CFVR0836 - Combination module test for CivilFEM results
Example to check the combinations utility
This example checks the functionality of the combinations module.
It checks the following:
- Results handle on BEAM3 elements.
- Date handle with unsorted numeration (elements, load steps...)
- Handle of all combination rules.
- Handle of nested combinations.
The load cases are:
- Load Step 2, Sub Step 1, Time 10: Surface load (1 N/m) on Span 1.
- Load Step 3, Sub Step 3, Time 30: Surface load (2 N/m) on Span 2.
- Load Step 6, Sub Step 2, Time 50: Surface load (3 N/m) on Span 3.
- Target 12: Absolute value MZ.
- Target 23: Minimum value of SX stress at point 26 of the cross section
- Target 25: Maximum value of X strain at point 57 of the cross section
- Combination rule 18: Addition [3.0*H1 + 3.0*H2]
- Combination rule 28: Addition [Combination 18 + 5.0*H3]
Element types used in the model: BEAM3 Needed CivilFEM Modules: |
|
| Model Statistics | |
| Number of elements | 6 |
| Number of nodes | 12 |
| Number of civil materials | 13 |
| Number of cross sections | 1 |
| Number of shell vertices | 0 |
Log file: CFVR0836.DAT
FINISH ~CFCLEAR,,1 AnsLic='ansys' NomFile='CFVR0836' /TITLE, %NomFile%, Combination module test for .RCV ! --------------------------------------------------------------------------------- ! Model definition and solve ! --------------------------------------------------------------------------------- ! Load scheme: ! LS2 SB1 TM10 Q=1 LS3 SB3 TM30 Q=2 LS6 SB2 TM50 Q=3 ! |||||||||||||||||| ||||||||||||||||| ||||||||||||||||||| ! |||||||||||||||||| ||||||||||||||||| ||||||||||||||||||| ! vvvvvvvvvvvvvvvvvv vvvvvvvvvvvvvvvvv vvvvvvvvvvvvvvvvvvv ! 2 4 6 8 10 12 14 : Nodes LS8 ! o=========o=========o=========o=========o=========o=========o ------> N=100 SB1 ! ^ 5 7 ^ 9 11 ^ 13 15 ^ : Elements ! 3 5 7 9 11 13 : Material ! 2 4 6 8 10 12 : Section ! 4 4 6 6 8 8 : Types ! 3 3 3 3 3 3 : Ename ! |------ 10 m -------|------ 16 m -------|------ 10 m -------| ! --------------------------------------------------------------------------------- ~UNITS,,LENG,M ~UNITS,,TIME,S ~UNITS,,FORC,KN /PREP7 ! Materials *DO,IMAT,2,14 ~CFMP,IMAT,LIB,STEEL,EA,A42 *ENDDO ! Element types ET,1,BEAM3 *DO,ITYP,2,10,4 ET,ITYP,BEAM3 *ENDDO *DO,ITYP,4,12,4 ET,ITYP,BEAM3 *ENDDO ! Sections ~SSECLIB,1,2,1,6 ! IPE 180 ! Beam & Shell Property *DO,ISEC,1,13 ~BMSHPRO,ISEC,BEAM,1,1,,,3,1,0, *ENDDO ! Nodes N, 1,-1 N, 2, 0 N, 4, 5 N, 6,10 N,10,26 FILL,6,10 N,14,36 FILL,10,14 ! Elements MAT, 3 $ REAL, 2 $ TYPE, 4 $ EN, 5, 2, 4 MAT, 5 $ REAL, 4 $ TYPE, 4 $ EN, 7, 4, 6 MAT, 7 $ REAL, 6 $ TYPE, 6 $ EN, 9, 6, 8 MAT, 9 $ REAL, 8 $ TYPE, 6 $ EN,11, 8,10 MAT,11 $ REAL,10 $ TYPE, 8 $ EN,13,10,12 MAT,13 $ REAL,12 $ TYPE, 8 $ EN,15,12,14 ! Boundary conditions D,2,UX D,2,UY D,6,UY D,10,UY D,14,UY ! Plot /ESHAPE,1 EPLOT ! Initial hypothesis solve /SOLU SOLCONTROL,OFF AUTOTS,OFF ! Hypothesis 0: LS1 Ignored TIME, 2 OUTRES,ALL,NONE SOLVE ! Hypothesis 1: LS2 SB1 TM10 Q=1 (1st Span) /TITLE, H1: LS2 SB1 TM10 Q=1 (1st Span) TIME,10 OUTRES,ALL,ALL SFBEAM, 5,1,PRES,1,1 SFBEAM, 7,1,PRES,1,1 SOLVE SFEDELE,ALL,ALL,PRES ! Hypothesis 2: LS3 SB3 TM30 Q=2 (2nd Span) /TITLE, H2: LS3 SB3 TM30 Q=2 (2nd Span) TIME,30 NSUBST,3 OUTRES,ALL,LAST SFBEAM, 9,1,PRES,2,2 SFBEAM,11,1,PRES,2,2 SOLVE SFEDELE,ALL,ALL,PRES ! Hypothesis 0: LS4 Ignored TIME,31 OUTRES,ALL,NONE SOLVE ! Hypothesis 0: LS5 Ignored TIME,32 OUTRES,ALL,NONE SOLVE ! Hypothesis 3: LS6 SB2 TM50 Q=3 (3rd Span) /TITLE, H3: LS6 SB2 TM50 Q=3 (3rd Span) TIME,50 NSUBST,2 OUTRES,ALL,ALL SFBEAM,13,1,PRES,3,3 SFBEAM,15,1,PRES,3,3 SOLVE SFEDELE,ALL,ALL,PRES ! Hypothesis 0: LS6 Ignored TIME,52 OUTRES,ALL,NONE SOLVE ! Hypothesis 4: LS8 SB1 TM60 N=100 /TITLE, H4: LS8 SB1 TM60 N=100 TIME,60 NSUBST,1 OUTRES,ALL,ALL F,14,FX,100 SOLVE FDELE,ALL,ALL /POST1 ! --------------------------------------------------------------------------------- ! Definition and solving of combinations ! --------------------------------------------------------------------------------- ! Targets: ! 12: MZ ABS in BEAMS ! 23: SBYT MIN in BEAMS ! 25: EPELDIR MAX in BEAMS ! Rules: ! 18 ADDITION: 3H1+2H2 ! --------------------------------------------------------------------------------- ! Reset combinations module ~CMBCLR, ! Targets definition ~TRGDEF,12,BEAM ,M ,Z ,ABS ~TRGDEF,23,BEAM ,S ,BYT ,MIN ~TRGDEF,25,BEAM ,EPEL ,BYT ,MAX ! Combinations definition ! Combination 18: ADDITION [3H1+2H2] ! ~CMBDEF , ICMB, Type, NSTST, KEYTITL, ITEM1 ~CMBDEF , 18, ADD, 2, ! ~STSTDEF, ICMB, ISTST, Type , ITEM1, ITEM2, ITIME, INCSTST, INCITEM1, INCITEM2 ~STSTDEF, 18, 1, LSTEP , 2 ! H1: LS2 SB1 TM10 Q=1 (1st Span) ~STSTDEF, 18, 2, LSTEP , 3 ! H2: LS3 SB3 TM30 Q=2 (2nd Span) ! ~STSTCFT, ICMB, ISTST, CFT1, CFT2 ~STSTCFT, 18, 1, 3.00, ~STSTCFT, 18, 2, 2.00, ! Combination 28: ADDITION CMB18+5H3 = [3H1+2H2]+5H3 ! ~CMBDEF , ICMB, Type, NSTST, KEYTITL, ITEM1 ~CMBDEF , 28, ADD, 2, ! ~STSTDEF, ICMB, ISTST, Type , ITEM1, ITEM2, ITIME, INCSTST, INCITEM1, INCITEM2 ~STSTDEF, 28, 1, CMB , 18 ! Cmb18 ~STSTDEF, 28, 2, LSTEP , 6 ! H3: LS6 SB2 TM50 Q=2 (3rd Span) ! ~STSTCFT, ICMB, ISTST, CFT1, CFT2 ~STSTCFT, 28, 1, 1.00, ~STSTCFT, 28, 2, 5.00, ! Combine ~COMBINE !-------------------------------------------------------------------------------------- ! DATA CHECK !-------------------------------------------------------------------------------------- ! Data comparison number NComp = 6 NComp_ch = 0 ! Matrix 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 ! Correct values !-------------------------------------------------------------------------------------- VALUE( 1,1) = 62.250 VALUE( 2,1) = -425398.63 VALUE( 3,1) = -0.202570778E-2 VALUE( 4,1) = 34.6764706 VALUE( 5,1) = -236969.044 VALUE( 6,1) = -0.112842403E-2 ! Obtained values !-------------------------------------------------------------------------------------- ~CMBDAT,2 ! Test 1: ~CFSET,,18,12 ETABLE,MZ_I,SMISC,6 ETABLE,MZ_J,SMISC,12 ! Comparison *GET,VALUE(1,2),ELEM,9,ETAB,MZ_J ! Test 2: ~CFSET,,18,23 ETABLE,SX_I,LS ,2 ETABLE,SX_J,LS ,5 ! Comparison *GET,VALUE(2,2),ELEM,9,ETAB,SX_J ! Test 3: ~CFSET,,18,25 ! Comparison ETABLE,EP_I,LEPEL,2 ETABLE,EP_J,LEPEL,5 *GET,VALUE(3,2),ELEM,9,ETAB,EP_J ! Test 4: ~CFSET,,28,12 ! Comparison ETABLE,MZ_I,SMISC,6 ETABLE,MZ_J,SMISC,12 *GET,VALUE(4,2),ELEM,9,ETAB,MZ_J ! Test 5: ~CFSET,,28,23 ! Comparison ETABLE,SX_I,LS ,2 ETABLE,SX_J,LS ,5 *GET,VALUE(5,2),ELEM,9,ETAB,SX_J ! Test 6: ~CFSET,,28,25 ! Comparison ETABLE,EP_I,LEPEL,2 ETABLE,EP_J,LEPEL,5 *GET,VALUE(6,2),ELEM,9,ETAB,EP_J ! Labels !-------------------------------------------------------------------------------------- LABEL( 1) ='TEST1' LABEL( 2) ='TEST2' LABEL( 3) ='TEST3' LABEL( 4) ='TEST4' LABEL( 5) ='TEST5' LABEL( 6) ='TEST6' ! Warning and error tolerances TOLER( 1, 1)= 1E-5 $ TOLER( 1, 2)= 1E-4 TOLER( 2, 1)= 1E-1 $ TOLER( 2, 2)= 1E-1 TOLER( 3, 1)= 1E-5 $ TOLER( 3, 2)= 1E-4 TOLER( 4, 1)= 1E-5 $ TOLER( 4, 2)= 1E-4 TOLER( 5, 1)= 1E-1 $ TOLER( 5, 2)= 1E-1 TOLER( 6, 1)= 1E-5 $ TOLER( 6, 2)= 1E-4 !-------------------------------------------------------------------------------------- ! Results comparison !-------------------------------------------------------------------------------------- COMPARA.MAC |
Results
| Label | Target | CivilFEM | Ratio | Tolerance |
| TEST1 | 62.25 | 62.25 | 1.000 | 0.0001 |
| TEST2 | -4.254e+005 | -4.254e+005 | 1.000 | 0.1 |
| TEST3 | -0.0020257 | -0.0020257 | 1.000 | 0.0001 |
| TEST4 | 34.676 | 34.676 | 1.000 | 0.0001 |
| TEST5 | -2.3697e+005 | -2.3697e+005 | 1.000 | 0.1 |
| TEST6 | -0.0011284 | -0.0011284 | 1.000 | 0.0001 |
Contains proprietary and confidential information of Ingeciber, S.A.