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

CFVR2959 - Shell shear design according to SP63133 (compression & tension)

Design Shear plus axial force of a reinforced concrete shell according to SP63133

Element types used in the model: SHELL_181

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
Shell Vertex
Reinforced Concrete
Shear design
SP63133

Model Statistics
Number of elements 18
Number of nodes 28
Number of civil materials 2
Number of cross sections 0
Number of shell vertices 1

Log file: CFVR2959.DAT


   FINISH
  ~CFCLEAR,,1
   NomFile='CFVR2959'
  /TITLE, %NomFile%, Shell shear design with SP63133 (compression & tension)

! ------------------------------------------------------------------------------
! Model definition and solve
! ------------------------------------------------------------------------------

! Parameters
   A = 3
   B = 6
   T = 0.2         ! Thickness
   R = 0.03        ! Cover 
   AST = 3.0708e-3
   ASB = 6.296e-3
   EsIZ = 1.0
   S =0.1
   ASW = 0.00022831/S

   
! Code
~CODESEL,,SP63133,

! Units
 ~UNITS,,FORC,KN 
  
/PREP7
! Material
~CFMP,1,LIB,CONCRETE,SP63133,B25
~CFMP,2,LIB,REINF,SP63133,A400  
~CFMP,1,SP_C    ,GAMb    ,,18.5/14.5
~CFMP,1,SP_C    ,GAMbt   ,,1.55/1.05
  
! Element types
ET,1,181
KEYOPT,1,3,2
KEYOPT,1,8,2
   
! Shell Vertex 
~SHLRNF,1,T,1,2,R,AST,ASB,AST,ASB,0,0.000,45.00
~SHLSHR,1,2,,,,,,,ASW,ASW 

! Shell properties
~BMSHPRO,1,SHELL,1,1,1,1, 181,,, ,Shell 1   
  
! Geometry
RECTNG,0,A,0,B
  
! Mesh
AESIZE,1,EsIZ,  
AMESH,1
  
FINISH
/SOLU 
 
! BCs
NSEL,S,LOC,X,0
NSEL,A,LOC,X,A
D,ALL,UX  
NSEL,S,LOC,Y,0
D,ALL, , , , , ,UX,UY,UZ,ROTZ, ,   
NSEL,S,LOC,Y,B
D,ALL, , , , , ,UX,UZ,ROTZ, , ,
ALLSEL
  
! CASE 1: COMPRESSION 1
Q = 60
H = 100

! Vertical Pressure
!SFE,ALL,1,PRES, ,-Q

! Compression pressure to forces
!NSEL,S,LOC,Y,B
!*GET,NNODE,NODE,,COUNT
!F,ALL,FY,-H*A/NNODE

ALLSEL
SOLVE		! CASE 1

FDELE,ALL,ALL
SFEDELE,ALL,ALL,ALL

! CASE 2: COMPRESSION 2
Q = 60
H = 2000

! Vertical Pressure
!SFE,ALL,1,PRES, ,-Q

! Compression pressure to forces
!NSEL,S,LOC,Y,B
!*GET,NNODE,NODE,,COUNT
!F,ALL,FY,-H*A/NNODE

ALLSEL
SOLVE		! CASE 2

FDELE,ALL,ALL
SFEDELE,ALL,ALL,ALL

! CASE 3:  COMPRESSION
Q = 60
H = 2320

! Vertical Pressure
! SFE,ALL,1,PRES, ,-Q

! Compression pressure to forces
! NSEL,S,LOC,Y,B
! *GET,NNODE,NODE,,COUNT
! F,ALL,FY,-H*A/NNODE

ALLSEL
SOLVE		! CASE 3

FDELE,ALL,ALL
SFEDELE,ALL,ALL,ALL

! CASE 4:  TENSION
Q = 60
H = 200

! Vertical Pressure
! SFE,ALL,1,PRES, ,-Q

! Compression pressure to forces
! NSEL,S,LOC,Y,B
! *GET,NNODE,NODE,,COUNT
! F,ALL,FY,H*A/NNODE

ALLSEL
SOLVE		! CASE 4

FDELE,ALL,ALL
SFEDELE,ALL,ALL,ALL


FINISH
/POST1

! CASE 1
~CFSET,,1

ESEL,S, , ,     15  ! Lateral endside element
~FMREAD,CFVR2959,FM 
~DIMCON,SHELL,SHEAR

! CASE 2
~CFSET,,2

ESEL,S, , ,     16  ! Lateral endside element
~FMREAD,CFVR2959,FM 
~DIMCON,SHELL,SHEAR

! CASE 3
~CFSET,,3

ESEL,S, , ,     17  ! Lateral endside element
~FMREAD,CFVR2959,FM 
~DIMCON,SHELL,SHEAR

! CASE 4
~CFSET,,4

ESEL,S, , ,     18  ! Lateral endside element
~FMREAD,CFVR2959,FM 
~DIMCON,SHELL,SHEAR
 
 
 
!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp = 40
  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

! Labels
!--------------------------------------------------------------------------------------
 
! CASE 1 DESIGN
  LABEL(1 ) = 'D1_Q1_X'
  LABEL(2 ) = 'D1_QB_X'
  LABEL(3 ) = 'D1_QSW_X'
  LABEL(4 ) = 'D1_ASSH_X'
  LABEL(5 ) = 'D1_PHI_N_X'
  LABEL(6 ) = 'D1_Q1_Y'
  LABEL(7 ) = 'D1_QB_Y'
  LABEL(8 ) = 'D1_QSW_Y'
  LABEL(9 ) = 'D1_ASSH_Y'
  LABEL(10) = 'D1_PHI_N_Y'
    
! CASE 2 DESIGN
  LABEL(11) = 'D2_Q1_X'
  LABEL(12) = 'D2_QB_X'
  LABEL(13) = 'D2_QSW_X'
  LABEL(14) = 'D2_ASSH_X'
  LABEL(15) = 'D2_PHI_N_X'
  LABEL(16) = 'D2_Q1_Y'
  LABEL(17) = 'D2_QB_Y'
  LABEL(18) = 'D2_QSW_Y'
  LABEL(19) = 'D2_ASSH_Y'
  LABEL(20) = 'D2_PHI_N_Y'

! CASE 3 DESIGN
  LABEL(21) = 'D3_Q1_X'
  LABEL(22) = 'D3_QB_X'
  LABEL(23) = 'D3_QSW_X'
  LABEL(24) = 'D3_ASSH_X'
  LABEL(25) = 'D3_PHI_N_X'
  LABEL(26) = 'D3_Q1_Y'
  LABEL(27) = 'D3_QB_Y'
  LABEL(28) = 'D3_QSW_Y'
  LABEL(29) = 'D3_ASSH_Y'
  LABEL(30) = 'D3_PHI_N_Y'
  
! CASE 4 DESIGN
  LABEL(31) = 'D4_Q1_X'
  LABEL(32) = 'D4_QB_X'
  LABEL(33) = 'D4_QSW_X'
  LABEL(34) = 'D4_ASSH_X'
  LABEL(35) = 'D4_PHI_N_X'
  LABEL(36) = 'D4_Q1_Y'
  LABEL(37) = 'D4_QB_Y'
  LABEL(38) = 'D4_QSW_Y'
  LABEL(39) = 'D4_ASSH_Y'
  LABEL(40) = 'D4_PHI_N_Y'

! Correct values
  
! CASE 1 DESIGN
  VALUE(1 ,1)= 744.677
  VALUE(2 ,1)= 89.250
  VALUE(3 ,1)= 113.73
  VALUE(4 ,1)= 2.4e-3
  VALUE(5 ,1)= 1.007
  VALUE(6 ,1)= 764.64
  VALUE(7 ,1)= 89.250
  VALUE(8 ,1)= 292.28
  VALUE(9 ,1)= 6.18e-3
  VALUE(10,1)= 1.034
  
 
! CASE 2 DESIGN
  VALUE(11,1)= 841.551
  VALUE(12,1)= 89.250
  VALUE(13,1)= 90.36
  VALUE(14,1)= 1.9e-3
  VALUE(15,1)= 1.138
  VALUE(16,1)= 924.375
  VALUE(17,1)= 89.250
  VALUE(18,1)= 226.35
  VALUE(19,1)= 4.78e-3
  VALUE(20,1)= 1.25

! CASE 3 DESIGN
  VALUE(21,1)= 857.8
  VALUE(22,1)= 89.250
  VALUE(23,1)= 86.9
  VALUE(24,1)= 1.84e-3
  VALUE(25,1)= 1.16
  VALUE(26,1)= 739.5
  VALUE(27,1)= 89.250
  VALUE(28,1)= 305.25
  VALUE(29,1)= 6.45e-3
  VALUE(30,1)= 1.0
  
! CASE 4 DESIGN
  VALUE(31,1)= 669
  VALUE(32,1)= 89.250
  VALUE(33,1)= 136.61
  VALUE(34,1)= 2.88e-3
  VALUE(35,1)= 0.905
  VALUE(36,1)= 387.5
  VALUE(37,1)= 89.250
  VALUE(38,1)= 2**100
  VALUE(39,1)= 2**100
  VALUE(40,1)= 0.524
 
! Comparison
 ~CFSET,1

 ~CFGET,VALUE(1 ,2),ELEMENT,15,RESULT,Q1_X,I
 ~CFGET,VALUE(2 ,2),ELEMENT,15,RESULT,QB_X,I
 ~CFGET,VALUE(3 ,2),ELEMENT,15,RESULT,QSW_X,I
 ~CFGET,VALUE(4 ,2),ELEMENT,15,RESULT,ASSH_X,I
 ~CFGET,VALUE(5 ,2),ELEMENT,15,RESULT,PHI_N_X,I
 ~CFGET,VALUE(6 ,2),ELEMENT,15,RESULT,Q1_Y,I
 ~CFGET,VALUE(7 ,2),ELEMENT,15,RESULT,QB_Y,I
 ~CFGET,VALUE(8 ,2),ELEMENT,15,RESULT,QSW_Y,I
 ~CFGET,VALUE(9 ,2),ELEMENT,15,RESULT,ASSH_Y,I
 ~CFGET,VALUE(10,2),ELEMENT,15,RESULT,PHI_N_Y,I
 
 ~CFSET,2
 ~CFGET,VALUE(11,2),ELEMENT,16,RESULT,Q1_X,I
 ~CFGET,VALUE(12,2),ELEMENT,16,RESULT,QB_X,I
 ~CFGET,VALUE(13,2),ELEMENT,16,RESULT,QSW_X,I
 ~CFGET,VALUE(14,2),ELEMENT,16,RESULT,ASSH_X,I
 ~CFGET,VALUE(15,2),ELEMENT,16,RESULT,PHI_N_X,I
 ~CFGET,VALUE(16,2),ELEMENT,16,RESULT,Q1_Y,I
 ~CFGET,VALUE(17,2),ELEMENT,16,RESULT,QB_Y,I
 ~CFGET,VALUE(18,2),ELEMENT,16,RESULT,QSW_Y,I
 ~CFGET,VALUE(19,2),ELEMENT,16,RESULT,ASSH_Y,I
 ~CFGET,VALUE(20,2),ELEMENT,16,RESULT,PHI_N_Y,I
  
 ~CFSET,3
 ~CFGET,VALUE(21,2),ELEMENT,17,RESULT,Q1_X,I
 ~CFGET,VALUE(22,2),ELEMENT,17,RESULT,QB_X,I
 ~CFGET,VALUE(23,2),ELEMENT,17,RESULT,QSW_X,I
 ~CFGET,VALUE(24,2),ELEMENT,17,RESULT,ASSH_X,I
 ~CFGET,VALUE(25,2),ELEMENT,17,RESULT,PHI_N_X,I
 ~CFGET,VALUE(26,2),ELEMENT,17,RESULT,Q1_Y,I
 ~CFGET,VALUE(27,2),ELEMENT,17,RESULT,QB_Y,I
 ~CFGET,VALUE(28,2),ELEMENT,17,RESULT,QSW_Y,I
 ~CFGET,VALUE(29,2),ELEMENT,17,RESULT,ASSH_Y,I
 ~CFGET,VALUE(30,2),ELEMENT,17,RESULT,PHI_N_Y,I
 
 ~CFSET,4
 ~CFGET,VALUE(31,2),ELEMENT,18,RESULT,Q1_X,I
 ~CFGET,VALUE(32,2),ELEMENT,18,RESULT,QB_X,I
 ~CFGET,VALUE(33,2),ELEMENT,18,RESULT,QSW_X,I
 ~CFGET,VALUE(34,2),ELEMENT,18,RESULT,ASSH_X,I
 ~CFGET,VALUE(35,2),ELEMENT,18,RESULT,PHI_N_X,I
 ~CFGET,VALUE(36,2),ELEMENT,18,RESULT,Q1_Y,I
 ~CFGET,VALUE(37,2),ELEMENT,18,RESULT,QB_Y,I
 ~CFGET,VALUE(38,2),ELEMENT,18,RESULT,QSW_Y,I
 ~CFGET,VALUE(39,2),ELEMENT,18,RESULT,ASSH_Y,I
 ~CFGET,VALUE(40,2),ELEMENT,18,RESULT,PHI_N_Y,I

! Warning and error tolerances

! CASE 1 DESIGN
TOLER(1 ,1)= VALUE(1 ,1)*2E-2 $ TOLER(18,2)=VALUE(1 ,1)*2E-2 
TOLER(2 ,1)= VALUE(2 ,1)*1E-2 $ TOLER(19,2)=VALUE(2 ,1)*1E-2 
TOLER(3 ,1)= VALUE(3 ,1)*1E-2 $ TOLER(20,2)=VALUE(3 ,1)*1E-2
TOLER(4 ,1)= VALUE(4 ,1)*1E-2 $ TOLER(21,2)=VALUE(4 ,1)*1E-2
TOLER(5 ,1)= VALUE(5 ,1)*1E-2 $ TOLER(22,2)=VALUE(5 ,1)*1E-2
TOLER(6 ,1)= VALUE(6 ,1)*1E-2 $ TOLER(23,2)=VALUE(6 ,1)*1E-2
TOLER(7 ,1)= VALUE(7 ,1)*1E-2 $ TOLER(24,2)=VALUE(7 ,1)*1E-2
TOLER(8 ,1)= VALUE(8 ,1)*1E-2 $ TOLER(25,2)=VALUE(8 ,1)*1E-2 
TOLER(9 ,1)= VALUE(9 ,1)*1E-2 $ TOLER(26,2)=VALUE(9 ,1)*1E-2 
TOLER(10,1)= VALUE(10,1)*1E-2 $ TOLER(27,2)=VALUE(10,1)*1E-2

! CASE 2 DESIGN
TOLER(11,1)= VALUE(11,1)*2E-2 $ TOLER(11,2)=VALUE(11,1)*2E-2 
TOLER(12,1)= VALUE(12,1)*1E-2 $ TOLER(12,2)=VALUE(12,1)*1E-2 
TOLER(13,1)= VALUE(13,1)*2E-2 $ TOLER(13,2)=VALUE(13,1)*2E-2
TOLER(14,1)= VALUE(14,1)*1E-2 $ TOLER(14,2)=VALUE(14,1)*1E-2
TOLER(15,1)= VALUE(15,1)*1E-2 $ TOLER(15,2)=VALUE(15,1)*1E-2
TOLER(16,1)= VALUE(16,1)*1E-2 $ TOLER(16,2)=VALUE(16,1)*1E-2
TOLER(17,1)= VALUE(17,1)*1E-2 $ TOLER(17,2)=VALUE(17,1)*1E-2
TOLER(18,1)= VALUE(18,1)*1E-2 $ TOLER(18,2)=VALUE(18,1)*1E-2 
TOLER(19,1)= VALUE(19,1)*2E-2 $ TOLER(19,2)=VALUE(19,1)*2E-2 
TOLER(20,1)= VALUE(20,1)*1E-2 $ TOLER(20,2)=VALUE(20,1)*1E-2

! CASE 3 DESIGN                       
TOLER(21,1)= VALUE(21,1)*2E-2 $ TOLER(21,2)=VALUE(21,1)*2E-2 
TOLER(22,1)= VALUE(22,1)*1E-2 $ TOLER(22,2)=VALUE(22,1)*1E-2 
TOLER(23,1)= VALUE(23,1)*2E-2 $ TOLER(23,2)=VALUE(23,1)*2E-2
TOLER(24,1)= VALUE(24,1)*1E-2 $ TOLER(24,2)=VALUE(24,1)*1E-2
TOLER(25,1)= VALUE(25,1)*1E-2 $ TOLER(25,2)=VALUE(25,1)*1E-2
TOLER(26,1)= VALUE(26,1)*1E-2 $ TOLER(26,2)=VALUE(26,1)*1E-2
TOLER(27,1)= VALUE(27,1)*1E-2 $ TOLER(27,2)=VALUE(27,1)*1E-2
TOLER(28,1)= VALUE(28,1)*1E-2 $ TOLER(28,2)=VALUE(28,1)*1E-2 
TOLER(29,1)= VALUE(29,1)*2E-2 $ TOLER(29,2)=VALUE(29,1)*2E-2 
TOLER(30,1)= VALUE(30,1)*1E-2 $ TOLER(30,2)=VALUE(30,1)*1E-2

! CASE 4 DESIGN                       
TOLER(31,1)= VALUE(31,1)*2E-2 $ TOLER(31,2)=VALUE(31,1)*2E-2 
TOLER(32,1)= VALUE(32,1)*1E-2 $ TOLER(32,2)=VALUE(32,1)*1E-2 
TOLER(33,1)= VALUE(33,1)*2E-2 $ TOLER(33,2)=VALUE(33,1)*2E-2
TOLER(34,1)= VALUE(34,1)*1E-2 $ TOLER(34,2)=VALUE(34,1)*1E-2
TOLER(35,1)= VALUE(35,1)*1E-2 $ TOLER(35,2)=VALUE(35,1)*1E-2
TOLER(36,1)= VALUE(36,1)*1E-2 $ TOLER(36,2)=VALUE(36,1)*1E-2
TOLER(37,1)= VALUE(37,1)*1E-2 $ TOLER(37,2)=VALUE(37,1)*1E-2
TOLER(38,1)= VALUE(38,1)*1E-2 $ TOLER(38,2)=VALUE(38,1)*1E-2 
TOLER(39,1)= VALUE(39,1)*2E-2 $ TOLER(39,2)=VALUE(39,1)*2E-2 
TOLER(40,1)= VALUE(40,1)*1E-2 $ TOLER(40,2)=VALUE(40,1)*1E-2
 
!--------------------------------------------------------------------------------------
! Results comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC

Results

LabelTargetCivilFEMRatioTolerance
D1_Q1_X 744.68 744.61.000 0
D1_QB_X 89.25 89.251.000 0
D1_QSW_X 113.73 113.751.000 0
D1_ASSH_ 0.0024 0.00240460.998 0
D1_PHI_N 1.007 1.00691.000 0
D1_Q1_Y 764.64 7651.000 0
D1_QB_Y 89.25 89.251.000 0
D1_QSW_Y 292.28 292.11.001 0
D1_ASSH_ 0.00618 0.00617491.001 0
D1_PHI_N 1.034 1.03451.000 0
D2_Q1_X 841.55 841.51.000 16.83
D2_QB_X 89.25 89.251.000 0.8925
D2_QSW_X 90.36 90.3741.000 1.807
D2_ASSH_ 0.0019 0.00191050.9951.9e-005
D2_PHI_N 1.138 1.13791.000 0.01138
D2_Q1_Y 924.38 924.381.000 9.244
D2_QB_Y 89.25 89.251.000 0.8925
D2_QSW_Y 226.35 226.351.000 2.264
D2_ASSH_ 0.00478 0.0047850.9999.56e-005
D2_PHI_N 1.25 1.251.000 0.0125
D3_Q1_X 857.8 857.821.000 17.16
D3_QB_X 89.25 89.251.000 0.8925
D3_QSW_X 86.9 86.9570.999 1.738
D3_ASSH_ 0.00184 0.00183821.0011.84e-005
D3_PHI_N 1.16 1.161.000 0.0116
D3_Q1_Y 739.5 739.51.000 7.395
D3_QB_Y 89.25 89.251.000 0.8925
D3_QSW_Y 305.25 305.251.000 3.053
D3_ASSH_ 0.00645 0.00645291.0000.000129
D3_PHI_N 1 11.000 0.01
D4_Q1_X 669 669.071.000 13.38
D4_QB_X 89.25 89.251.000 0.8925
D4_QSW_X 136.61 136.671.000 2.732
D4_ASSH_ 0.00288 0.00288910.9972.88e-005
D4_PHI_N 0.905 0.904761.000 0.00905
D4_Q1_Y 387.5 387.361.000 3.875
D4_QB_Y 89.25 89.251.000 0.8925
D4_QSW_Y1.2677e+0301.2677e+0301.0001.268e+028
D4_ASSH_1.2677e+0301.2677e+0301.0002.535e+028
D4_PHI_N 0.524 0.523811.000 0.00524

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