Left

CivilFEM Online Help

Right

CFVR3903 - Testing fluxes 4 node 2D laplace element


Test fluxes in nodes of a PLANE42 2D model. The model consist of both soil triangle and square.
For this purpose it calculates and applies a ground pressure on both plane elements.

Element types used in the model: PLANE42

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 111
Number of civil materials 1
Number of cross sections 0
Number of shell vertices 0

Log file: CFVR3903.DAT


  FINISH
  ~CFCLEAR,,1
  NomFile='CFVR3903'
  /TITLE, %NomFile%, Testing fluxes 4 node 2D laplace element

  ~SEEPAGE,2D

/prep7

~cfmp,1,lib,soil,,sc
~cfmp,1,soil,kx,,10
~cfmp,1,soil,ky,,10

!      X   Y   Z
n,10, 10, 10,  0
n,20, 20, 10,  0
n,30, 20, 20,  0
n,40, 10, 20,  0
n,50, 50, 10,  0
n,60, 60, 10,  0
n,70, 55, 20,  0

e,10,20,30,40
e,50,60,70,70

/SOLU
d,10,head,100
d,20,head,100
d,30,head,10
d,40,head,10
d,70,head,100
d,60,head,10
d,50,head,10

~LPSOLVE
/post1
file,file,seep
set,1
/dscale,,off
plnsol,head

etable,fx1,smisc,1
etable,fx2,smisc,2
etable,fx3,smisc,3
etable,fx4,smisc,4

!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp    =  8
  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) = 'E2F1'
  LABEL(6,1) = 'E2F2'
  LABEL(7,1) = 'E2F3'
  LABEL(8,1) = 'E2F4'

! Correct values
  VALUE(1,1)= -900.0
  VALUE(2,1)= 0.0
  VALUE(3,1)= 900.0
  VALUE(4,1)= 0.0
  VALUE(5,1)= 900.0
  VALUE(6,1)= -450.0
  VALUE(7,1)= 0.0
  VALUE(8,1)= -450.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,2,etable,fx1
  *get,VALUE(6,2),element,2,etable,fx2
  *get,VALUE(7,2),element,2,etable,fx3
  *get,VALUE(8,2),element,2,etable,fx4

! Warning and error tolerances
  *DO,II,1,8
    TOLER(II, 1)= 1E-9
    TOLER(II, 2)= 1E-9
  *ENDDO

!--------------------------------------------------------------------------------------
! Results comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC

Results

LabelTargetCivilFEMRatioTolerance
E1F1 -900 02.000 1e-009
E1F2 0 00.000 1e-009
E1F3 900 02.000 1e-009
E1F4 0 00.000 1e-009
E2F1 900 02.000 1e-009
E2F2 -450 02.000 1e-009
E2F3 0 00.000 1e-009
E2F4 -450 02.000 1e-009


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