Left

CivilFEM Online Help

Official Websites: CivilFEM | CivilFEM Licenses | Ingeciber
CivilFEM powered by Marc: Installation | Documentation | Python Guide | Theory Manual | Verification Manual | PDF | Class 3
CivilFEM for Ansys: Installation | Documentation | Verification Manual | Instalación | Documentación | PDF | PDF (ES) | Class 3
CivilFEM Apps for Ansys Workbench: Installation | Licensing | FAQ | Reinforced Concrete | Structural Steel | Theory Manual | PDF | Class 3
Right

CFVR0838 - Combination module test

Example to check the combinations utility

This example checks the functionality of the combinations module.
It checks the following:

  • Results handle on SHELL181 elements.
  • Date handle with unsorted numeration (elements, load steps...)

The model is a two elements horizontal slab, each element is 1.00m * 1.50m. All the nodes are restrained (no movement or rotation allowed).
The load cases are:
  • Load Step 1, Sub Step 1: Imposed vertical movement on the middle point of one side (node 50) of 1 meter.
  • Load Step 2, Sub Step 1: Imposed horizontal movement on the middle point of one side (node 50) of 1 meter (Y direction).
No Civil materials or shell properties have been used in this example, only ANSYS results will be checked.

The target for the combinations is:
  • Target 45; Maximum Axial force in the Y direction TY

The combination rule defined is:
  • Combination rule 10: ADDITION: H1+H2

Element types used in the model: SHELL181

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
Shell
Combinations

Model Statistics
Number of elements 2
Number of nodes 6
Number of civil materials 0
Number of cross sections 0
Number of shell vertices 0

Log file: CFVR0838.DAT


  FINISH
  ~CFCLEAR,,1 
 AnsLic='ansys'

  NomFile='CFVR0838'
  /TITLE, %NomFile%, Combination module test for SHELL181
! ---------------------------------------------------------------------------------
! Model definition and solve
! ---------------------------------------------------------------------------------
/PREP7
! Materials
  MP,EX,5,3E6
  MP,NUXY,5,0.3
! Real constants
  R,10,0.1,0.1,0.1,0.1
! Element types
  ET,2,SHELL181
! Nodes
  N,10,0.0,0.0
  N,20,1.0,0.0
  N,30,2.0,0.0
  N,40,0.0,1.5
  N,50,1.0,1.5
  N,60,2.0,1.5
! Elements
  MAT,5
  REAL,10
  TYPE,2
  EN,100, 10, 20, 50, 40
  EN,200, 20, 30, 60, 50
! Boundary conditions
  D,ALL,ALL,0.0
! Load cases
  /SOLU
  ! Load case 1:
    /TITLE, Load case 1
    OUTRES,ALL,ALL
    D,50,UZ,1
    SOLVE
    D,50,UZ,0.0
  ! Load case 2:
    /TITLE, Load case 2
    OUTRES,ALL,ALL
    D,50,UY,1
    SOLVE
    D,50,UZ,0.0
/POST1
! ---------------------------------------------------------------------------------
! Load case results
! ---------------------------------------------------------------------------------
  Dat=44 ! Number of data to check
  RES=
  *DIM,RES,,Dat,2
! Load case results data retrieval
  ! Load case 1
    SET,1,1
    ETABLE,TX ,SMISC,1
    ETABLE,TY ,SMISC,2
    ETABLE,TXY,SMISC,3
    ETABLE,MX ,SMISC,4
    ETABLE,MY ,SMISC,5
    ETABLE,MXY,SMISC,6
    ETABLE,NX ,SMISC,7
    ETABLE,NY ,SMISC,8
    *GET,RES( 1,1),ELEM,100,ETAB,TX
    *GET,RES( 2,1),ELEM,100,ETAB,TY
    *GET,RES( 3,1),ELEM,100,ETAB,TXY
    *GET,RES( 4,1),ELEM,100,ETAB,MX
    *GET,RES( 5,1),ELEM,100,ETAB,MY
    *GET,RES( 6,1),ELEM,100,ETAB,MXY
    SHELL,TOP
    *GET,RES( 7,1),NODE,20,S,X
    *GET,RES( 8,1),NODE,20,S,Y
    *GET,RES( 9,1),NODE,20,S,Z
    *GET,RES(10,1),NODE,20,S,XY
    *GET,RES(11,1),NODE,20,S,YZ
    *GET,RES(12,1),NODE,20,S,XZ
    *GET,RES(13,1),NODE,20,EPEL,X
    *GET,RES(14,1),NODE,20,EPEL,Y
    *GET,RES(15,1),NODE,20,EPEL,Z
    *GET,RES(16,1),NODE,20,EPEL,XY
    *GET,RES(17,1),NODE,20,EPEL,YZ
    *GET,RES(18,1),NODE,20,EPEL,XZ
    SHELL,MID
    *GET,RES(19,1),NODE,20,S,X
    *GET,RES(20,1),NODE,20,S,Y
    *GET,RES(21,1),NODE,20,S,Z
    *GET,RES(22,1),NODE,20,S,XY
    *GET,RES(23,1),NODE,20,S,YZ
    *GET,RES(24,1),NODE,20,S,XZ
    *GET,RES(25,1),NODE,20,EPEL,X
    *GET,RES(26,1),NODE,20,EPEL,Y
    *GET,RES(27,1),NODE,20,EPEL,Z
    *GET,RES(28,1),NODE,20,EPEL,XY
    *GET,RES(29,1),NODE,20,EPEL,YZ
    *GET,RES(30,1),NODE,20,EPEL,XZ
    SHELL,BOT
    *GET,RES(31,1),NODE,20,S,X
    *GET,RES(32,1),NODE,20,S,Y
    *GET,RES(33,1),NODE,20,S,Z
    *GET,RES(34,1),NODE,20,S,XY
    *GET,RES(35,1),NODE,20,S,YZ
    *GET,RES(36,1),NODE,20,S,XZ
    *GET,RES(37,1),NODE,20,EPEL,X
    *GET,RES(38,1),NODE,20,EPEL,Y
    *GET,RES(39,1),NODE,20,EPEL,Z
    *GET,RES(40,1),NODE,20,EPEL,XY
    *GET,RES(41,1),NODE,20,EPEL,YZ
    *GET,RES(42,1),NODE,20,EPEL,XZ
    *GET,RES(43,1),ELEM,100,ETAB,NX
    *GET,RES(44,1),ELEM,100,ETAB,NY
  ! Load case 2
    SET,2,1
    ETABLE,TX ,SMISC,1
    ETABLE,TY ,SMISC,2
    ETABLE,TXY,SMISC,3
    ETABLE,MX ,SMISC,4
    ETABLE,MY ,SMISC,5
    ETABLE,MXY,SMISC,6
    ETABLE,NX ,SMISC,7
    ETABLE,NY ,SMISC,8
    *GET,RES( 1,2),ELEM,100,ETAB,TX
    *GET,RES( 2,2),ELEM,100,ETAB,TY
    *GET,RES( 3,2),ELEM,100,ETAB,TXY
    *GET,RES( 4,2),ELEM,100,ETAB,MX
    *GET,RES( 5,2),ELEM,100,ETAB,MY
    *GET,RES( 6,2),ELEM,100,ETAB,MXY
    SHELL,TOP
    *GET,RES( 7,2),NODE,20,S,X
    *GET,RES( 8,2),NODE,20,S,Y
    *GET,RES( 9,2),NODE,20,S,Z
    *GET,RES(10,2),NODE,20,S,XY
    *GET,RES(11,2),NODE,20,S,YZ
    *GET,RES(12,2),NODE,20,S,XZ
    *GET,RES(13,2),NODE,20,EPEL,X
    *GET,RES(14,2),NODE,20,EPEL,Y
    *GET,RES(15,2),NODE,20,EPEL,Z
    *GET,RES(16,2),NODE,20,EPEL,XY
    *GET,RES(17,2),NODE,20,EPEL,YZ
    *GET,RES(18,2),NODE,20,EPEL,XZ
    SHELL,MID
    *GET,RES(19,2),NODE,20,S,X
    *GET,RES(20,2),NODE,20,S,Y
    *GET,RES(21,2),NODE,20,S,Z
    *GET,RES(22,2),NODE,20,S,XY
    *GET,RES(23,2),NODE,20,S,YZ
    *GET,RES(24,2),NODE,20,S,XZ
    *GET,RES(25,2),NODE,20,EPEL,X
    *GET,RES(26,2),NODE,20,EPEL,Y
    *GET,RES(27,2),NODE,20,EPEL,Z
    *GET,RES(28,2),NODE,20,EPEL,XY
    *GET,RES(29,2),NODE,20,EPEL,YZ
    *GET,RES(30,2),NODE,20,EPEL,XZ
    SHELL,BOT
    *GET,RES(31,2),NODE,20,S,X
    *GET,RES(32,2),NODE,20,S,Y
    *GET,RES(33,2),NODE,20,S,Z
    *GET,RES(34,2),NODE,20,S,XY
    *GET,RES(35,2),NODE,20,S,YZ
    *GET,RES(36,2),NODE,20,S,XZ
    *GET,RES(37,2),NODE,20,EPEL,X
    *GET,RES(38,2),NODE,20,EPEL,Y
    *GET,RES(39,2),NODE,20,EPEL,Z
    *GET,RES(40,2),NODE,20,EPEL,XY
    *GET,RES(41,2),NODE,20,EPEL,YZ
    *GET,RES(42,2),NODE,20,EPEL,XZ
    *GET,RES(43,2),ELEM,100,ETAB,NX
    *GET,RES(44,2),ELEM,100,ETAB,NY

! ---------------------------------------------------------------------------------
! Combination definition and performance
! ---------------------------------------------------------------------------------
! Combination module Initialization
  ~CMBCLR,
! Targets definition
  ~TRGDEF,45,SHELL,T,Y,MAX
! Combination definition
  ! Combination 10: ADDITION [H1+H2]
      ~CMBDEF, 10,ADD,2
      ~STSTDEF,10,1,LSTEP,1, ,2,1,1
      ~STSTCFT,10, ,1.00
! Combinations
  ~COMBINE

!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp = 44
  NComp_ch = 0

! Marix dim.
  *DIM,LABEL,CHAR,Ncomp,1
  *DIM,VALUE,,Ncomp,3
  *DIM,TOLER,,Ncomp,2

! Labels
! ---------------------------------------------------------------------------------
  *DO,I,1,NComp
    LABEL(I) ='TEST%I%'
  *ENDDO

! Correct values
! ---------------------------------------------------------------------------------
  *DO,I,1,NComp
    VALUE(I,1)=RES(I,1)+RES(I,2)
  *ENDDO

! Obtained values
! ---------------------------------------------------------------------------------
  ~CMBDAT,2
  ~CFSET,,10,45
  ETABLE,TX ,SMISC,1
  ETABLE,TY ,SMISC,2
  ETABLE,TXY,SMISC,3
  ETABLE,MX ,SMISC,4
  ETABLE,MY ,SMISC,5
  ETABLE,MXY,SMISC,6
  ETABLE,NX ,SMISC,7
  ETABLE,NY ,SMISC,8
  *GET,VALUE( 1,2),ELEM,100,ETAB,TX
  *GET,VALUE( 2,2),ELEM,100,ETAB,TY
  *GET,VALUE( 3,2),ELEM,100,ETAB,TXY
  *GET,VALUE( 4,2),ELEM,100,ETAB,MX
  *GET,VALUE( 5,2),ELEM,100,ETAB,MY
  *GET,VALUE( 6,2),ELEM,100,ETAB,MXY
  SHELL,TOP
  *GET,VALUE( 7,2),NODE,20,S,X
  *GET,VALUE( 8,2),NODE,20,S,Y
  *GET,VALUE( 9,2),NODE,20,S,Z
  *GET,VALUE(10,2),NODE,20,S,XY
  *GET,VALUE(11,2),NODE,20,S,YZ
  *GET,VALUE(12,2),NODE,20,S,XZ
  *GET,VALUE(13,2),NODE,20,EPEL,X
  *GET,VALUE(14,2),NODE,20,EPEL,Y
  *GET,VALUE(15,2),NODE,20,EPEL,Z
  *GET,VALUE(16,2),NODE,20,EPEL,XY
  *GET,VALUE(17,2),NODE,20,EPEL,YZ
  *GET,VALUE(18,2),NODE,20,EPEL,XZ
  SHELL,MID
  *GET,VALUE(19,2),NODE,20,S,X
  *GET,VALUE(20,2),NODE,20,S,Y
  *GET,VALUE(21,2),NODE,20,S,Z
  *GET,VALUE(22,2),NODE,20,S,XY
  *GET,VALUE(23,2),NODE,20,S,YZ
  *GET,VALUE(24,2),NODE,20,S,XZ
  *GET,VALUE(25,2),NODE,20,EPEL,X
  *GET,VALUE(26,2),NODE,20,EPEL,Y
  *GET,VALUE(27,2),NODE,20,EPEL,Z
  *GET,VALUE(28,2),NODE,20,EPEL,XY
  *GET,VALUE(29,2),NODE,20,EPEL,YZ
  *GET,VALUE(30,2),NODE,20,EPEL,XZ
  SHELL,BOT
  *GET,VALUE(31,2),NODE,20,S,X
  *GET,VALUE(32,2),NODE,20,S,Y
  *GET,VALUE(33,2),NODE,20,S,Z
  *GET,VALUE(34,2),NODE,20,S,XY
  *GET,VALUE(35,2),NODE,20,S,YZ
  *GET,VALUE(36,2),NODE,20,S,XZ
  *GET,VALUE(37,2),NODE,20,EPEL,X
  *GET,VALUE(38,2),NODE,20,EPEL,Y
  *GET,VALUE(39,2),NODE,20,EPEL,Z
  *GET,VALUE(40,2),NODE,20,EPEL,XY
  *GET,VALUE(41,2),NODE,20,EPEL,YZ
  *GET,VALUE(42,2),NODE,20,EPEL,XZ
  *GET,VALUE(43,2),ELEM,100,ETAB,NX
  *GET,VALUE(44,2),ELEM,100,ETAB,NY

! Warning and error tolerances
  TOLER( 1, 1)= 1E-03   $   TOLER( 1, 2)= 1E-02
  TOLER( 2, 1)= 1E-03   $   TOLER( 2, 2)= 1E-02
  TOLER( 3, 1)= 1E-03   $   TOLER( 3, 2)= 1E-02
  TOLER( 4, 1)= 1E-20   $   TOLER( 4, 2)= 1E-20
  TOLER( 5, 1)= 1E-20   $   TOLER( 5, 2)= 1E-20
  TOLER( 6, 1)= 1E-20   $   TOLER( 6, 2)= 1E-20
  TOLER( 7, 1)= 1E-02   $   TOLER( 7, 2)= 1E-01
  TOLER( 8, 1)= 1E-02   $   TOLER( 8, 2)= 1E-01
  TOLER( 9, 1)= 1E-15   $   TOLER( 9, 2)= 1E-14
  TOLER(10, 1)= 1E-15   $   TOLER(10, 2)= 1E-14
  TOLER(11, 1)= 1E-03   $   TOLER(11, 2)= 1E-02
  TOLER(12, 1)= 1E-13   $   TOLER(12, 2)= 1E-12
  TOLER(13, 1)= 1E-14   $   TOLER(13, 2)= 1E-13
  TOLER(14, 1)= 1E-09   $   TOLER(14, 2)= 1E-08
  TOLER(15, 1)= 1E-10   $   TOLER(15, 2)= 1E-09
  TOLER(16, 1)= 1E-15   $   TOLER(16, 2)= 1E-14
  TOLER(17, 1)= 1E-10   $   TOLER(17, 2)= 1E-09
  TOLER(18, 1)= 1E-15   $   TOLER(18, 2)= 1E-14
  TOLER(19, 1)= 1E-02   $   TOLER(19, 2)= 1E-01
  TOLER(20, 1)= 1E-02   $   TOLER(20, 2)= 1E-01
  TOLER(21, 1)= 1E-14   $   TOLER(21, 2)= 1E-13
  TOLER(22, 1)= 1E-12   $   TOLER(22, 2)= 1E-11
  TOLER(23, 1)= 1E-03   $   TOLER(23, 2)= 1E-02
  TOLER(24, 1)= 1E-13   $   TOLER(24, 2)= 1E-11
  TOLER(25, 1)= 1E-14   $   TOLER(25, 2)= 1E-13
  TOLER(26, 1)= 1E-09   $   TOLER(26, 2)= 1E-08
  TOLER(27, 1)= 1E-09   $   TOLER(27, 2)= 1E-08
  TOLER(28, 1)= 1E-14   $   TOLER(28, 2)= 1E-13
  TOLER(29, 1)= 1E-10   $   TOLER(29, 2)= 1E-09
  TOLER(30, 1)= 1E-14   $   TOLER(30, 2)= 1E-13
  TOLER(31, 1)= 1E-02   $   TOLER(31, 2)= 1E-01
  TOLER(32, 1)= 1E-01   $   TOLER(32, 2)= 1E-00
  TOLER(33, 1)= 1E-14   $   TOLER(33, 2)= 1E-13
  TOLER(34, 1)= 1E-14   $   TOLER(34, 2)= 1E-13
  TOLER(35, 1)= 1E-03   $   TOLER(35, 2)= 1E-02
  TOLER(36, 1)= 1E-13   $   TOLER(36, 2)= 1E-12
  TOLER(37, 1)= 1E-14   $   TOLER(37, 2)= 1E-13
  TOLER(38, 1)= 1E-09   $   TOLER(38, 2)= 1E-08
  TOLER(39, 1)= 1E-09   $   TOLER(39, 2)= 1E-08
  TOLER(40, 1)= 1E-14   $   TOLER(40, 2)= 1E-13
  TOLER(41, 1)= 1E-10   $   TOLER(41, 2)= 1E-09
  TOLER(42, 1)= 1E-14   $   TOLER(42, 2)= 1E-13
  TOLER(41, 1)= 1E-04   $   TOLER(41, 2)= 1E-03
  TOLER(42, 1)= 1E-06   $   TOLER(42, 2)= 1E-05

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

Results

LabelTargetCivilFEMRatioTolerance
TEST1 32967 329671.000 0.01
TEST21.0989e+0051.0989e+0051.000 0.01
TEST3 57692 576921.000 0.01
TEST4-2.8422e-014-2.8422e-0141.000 1e-020
TEST51.1369e-0131.1369e-0131.000 1e-020
TEST6-5.6843e-014-5.6843e-0141.000 1e-020
TEST73.2967e+0053.2967e+0051.000 0.1
TEST81.0989e+0061.0989e+0061.000 0.1
TEST9 0 00.000 1e-014
TEST10 0 00.000 1e-014
TEST11 0 00.000 0.01
TEST12 0 00.000 1e-012
TEST13 0 00.000 1e-013
TEST14 0.33333 0.333331.000 1e-008
TEST15 -0.14286 -0.142861.000 1e-009
TEST16 0 00.000 1e-014
TEST17 0 00.000 1e-009
TEST18 0 00.000 1e-014
TEST193.2967e+0053.2967e+0051.000 0.1
TEST201.0989e+0061.0989e+0061.000 0.1
TEST21 0 00.000 1e-013
TEST22 0 00.000 1e-011
TEST23 0 00.000 0.01
TEST24 0 00.000 1e-011
TEST25 0 00.000 1e-013
TEST26 0.33333 0.333331.000 1e-008
TEST27 -0.14286 -0.142861.000 1e-008
TEST28 0 00.000 1e-013
TEST29 0 00.000 1e-009
TEST30 0 00.000 1e-013
TEST313.2967e+0053.2967e+0051.000 0.1
TEST321.0989e+0061.0989e+0061.000 1
TEST33 0 00.000 1e-013
TEST34 0 00.000 1e-013
TEST35 0 00.000 0.01
TEST36 0 00.000 1e-012
TEST37 0 00.000 1e-013
TEST38 0.33333 0.333331.000 1e-008
TEST39 -0.14286 -0.142861.000 1e-008
TEST40 0 00.000 1e-013
TEST41 0 00.000 0.001
TEST42 0 00.000 1e-005
TEST43 48032 480321.000 0
TEST44 32021 320211.000 0

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