CFVR3904 - Testing Fluxes 8 Node 3D Laplace Element
Test fluxed in nodes of a SOLID45 3D model.The model consist of three soil polyhedrons: tetrahedron, pentahedron, and a cube.
For this purpose it calculates and applies a ground pressure on the nodes of the model.
The example is launched with ANSYS and with CivilFEM separately.
Element types used in the model: SOLID45 Needed CivilFEM Modules: |
|
| Model Statistics | |
| Number of elements | 230 |
| Number of nodes | 117 |
| Number of civil materials | 1 |
| Number of cross sections | 0 |
| Number of shell vertices | 0 |
Log file: CFVR3904.DAT
FINISH ~CFCLEAR,,1 Nomfile='CFVR3904' /Title, %Nomfile%, Testing Fluxes 8 Node 3D Laplace Element ~SEEPAGE,3D /Prep7 ~Cfmp,1,Lib,Soil,,Sc ~Cfmp,1,Soil,Kx,,10 ~Cfmp,1,Soil,Ky,,20 ~Cfmp,1,Soil,Kz,,10 N,10,10,10,10 N,20,20,10,10 N,30,20,20,10 N,40,10,20,10 N,50,10,10,20 N,60,20,10,20 N,70,20,20,20 N,80,10,20,20 E,10,20,30,40,50,60,70,80 N,100,50,10,10 N,110,60,10,10 N,15 ,60,20,10 N,120,50,10,20 N,130,60,10,20 N,25 ,60,20,20 E,100,110,15,15,120,130,25,25 N,200,10,80,10 N,210,20,80,10 N,35 ,20,90,10 N,45 ,10,80,20 E,200,210,35,35,45,45,45,45 /Solu D,10,Head,100 D,20,Head,100 D,30,Head,100 D,40,Head,100 D,50,Head,10 D,60,Head,10 D,70,Head,10 D,80,Head,10 D,120,Head,10 D,100,Head,10 D,15 ,Head,100 D,25 ,Head,100 D,130,Head,10 D,110,Head,10 D,35 ,Head,100 D,45 ,Head,10 D,200,Head,10 D,210,Head,10 ~LPSOLVE /Post1 File,File,Seep /Dscale,,Off Set,1 Plnsol,Head Etable,Fx1,smisc,1 Etable,Fx2,smisc,2 Etable,Fx3,smisc,3 Etable,Fx4,smisc,4 Etable,Fx5,smisc,5 Etable,Fx6,smisc,6 !-------------------------------------------------------------------------------------- ! Data Check !-------------------------------------------------------------------------------------- ! Data Comparison Number Ncomp = 18 Ncomp_Ch = 0 ! Marix Dim. *Dim,Label,Char,Ncomp,1 *Dim,Label_Ch,Char,Ncomp_Ch,1 *Dim,Value,,Ncomp,3 *Dim,Value_Ch,Char,Ncomp_Ch,3 *Dim,Toler,,Ncomp,2 ! Labels Label(1 ,1) = 'e1F1' Label(2 ,1) = 'e1F2' Label(3 ,1) = 'e1F3' Label(4 ,1) = 'e1F4' Label(5 ,1) = 'e1F5' Label(6 ,1) = 'e1F6' Label(7 ,1) = 'e2F1' Label(8 ,1) = 'e2F2' Label(9 ,1) = 'e2F3' Label(10,1) = 'e2F4' Label(11,1) = 'e2F5' Label(12,1) = 'e2F6' Label(13,1) = 'e3F1' Label(14,1) = 'e3F2' Label(15,1) = 'e3F3' Label(16,1) = 'e3F4' Label(17,1) = 'e3F5' Label(18,1) = 'e3F6' ! Correct Values Value(1 ,1) = -9000.0 Value(2 ,1) = 0.0 Value(3 ,1) = 0.0 Value(4 ,1) = 0.0 Value(5 ,1) = 0.0 Value(6 ,1) = 9000.0 Value(7 ,1) = 0.0 Value(8 ,1) = 18000 Value(9 ,1) = 0.0 Value(10,1) = 0.0 Value(11,1) = -18000 Value(12,1) = 0.0 Value(13,1)= 0.0 Value(14,1)= 9000.0 Value(15,1)= 0.0 Value(16,1)= 0.0 Value(17,1)= -9000.0 Value(18,1)= 0.0 ! Obtained Values *Get,Value(1, 2),Element,1,Etable,Fx1 *Get,Value(2, 2),Element,1,Etable,Fx2 *Get,Value(3, 2),Element,1,Etable,Fx3 *Get,Value(4, 2),Element,1,Etable,Fx4 *Get,Value(5, 2),Element,1,Etable,Fx5 *Get,Value(6, 2),Element,1,Etable,Fx6 *Get,Value(7, 2),Element,2,Etable,Fx1 *Get,Value(8, 2),Element,2,Etable,Fx2 *Get,Value(9, 2),Element,2,Etable,Fx3 *Get,Value(10,2),Element,2,Etable,Fx4 *Get,Value(11,2),Element,2,Etable,Fx5 *Get,Value(12,2),Element,2,Etable,Fx6 *Get,Value(13,2),Element,3,Etable,Fx1 *Get,Value(14,2),Element,3,Etable,Fx2 *Get,Value(15,2),Element,3,Etable,Fx3 *Get,Value(16,2),Element,3,Etable,Fx4 *Get,Value(17,2),Element,3,Etable,Fx5 *Get,Value(18,2),Element,3,Etable,Fx6 ! Warning And Error Tolerances *Do,II,1,18 Toler(II, 1)= 1E-6 Toler(II, 2)= 1E-6 *Enddo !-------------------------------------------------------------------------------------- ! Results Comparison !-------------------------------------------------------------------------------------- Compara.Mac |
Results
| Label | Target | CivilFEM | Ratio | Tolerance |
| e1F1 | -9000 | 0 | 2.000 | 1e-006 |
| e1F2 | 0 | 0 | 0.000 | 1e-006 |
| e1F3 | 0 | 0 | 0.000 | 1e-006 |
| e1F4 | 0 | 0 | 0.000 | 1e-006 |
| e1F5 | 0 | 0 | 0.000 | 1e-006 |
| e1F6 | 9000 | 0 | 2.000 | 1e-006 |
| e2F1 | 0 | 0 | 0.000 | 1e-006 |
| e2F2 | 18000 | 0 | 2.000 | 1e-006 |
| e2F3 | 0 | 0 | 0.000 | 1e-006 |
| e2F4 | 0 | 0 | 0.000 | 1e-006 |
| e2F5 | -18000 | 0 | 2.000 | 1e-006 |
| e2F6 | 0 | 0 | 0.000 | 1e-006 |
| e3F1 | 0 | 0 | 0.000 | 1e-006 |
| e3F2 | 9000 | 0 | 2.000 | 1e-006 |
| e3F3 | 0 | 0 | 0.000 | 1e-006 |
| e3F4 | 0 | 0 | 0.000 | 1e-006 |
| e3F5 | -9000 | 0 | 2.000 | 1e-006 |
| e3F6 | 0 | 0 | 0.000 | 1e-006 |
Contains proprietary and confidential information of Ingeciber, S.A.

