Left

CivilFEM Online Help

Right

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:
Geotechnical Module
Bridges and Civil Non Linearities Module
Advanced Prestressed Concrete Module

The example can be launched on an educational license.
KEYWORDS
Material Models
Materials

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

LabelTargetCivilFEMRatioTolerance
e1F1 -9000 02.000 1e-006
e1F2 0 00.000 1e-006
e1F3 0 00.000 1e-006
e1F4 0 00.000 1e-006
e1F5 0 00.000 1e-006
e1F6 9000 02.000 1e-006
e2F1 0 00.000 1e-006
e2F2 18000 02.000 1e-006
e2F3 0 00.000 1e-006
e2F4 0 00.000 1e-006
e2F5 -18000 02.000 1e-006
e2F6 0 00.000 1e-006
e3F1 0 00.000 1e-006
e3F2 9000 02.000 1e-006
e3F3 0 00.000 1e-006
e3F4 0 00.000 1e-006
e3F5 -9000 02.000 1e-006
e3F6 0 00.000 1e-006


Contains proprietary and confidential information of Ingeciber, S.A.