CFVR9174 - Shell shear reinforcement check according to ACI 318 19

Shell shear reinforcement check according to ACI 318 19.

The model is made up of two square shell elements of 1.0m x 1.0m. Boundary conditions: + Node located at origin has its vertical movement restrained(perpendicular to the shells). + Node located at X=0 and Y=100 cm are completely restrained. + Node located at X=100 and Y=100 cm has its vertical movement restrained(perpendicular to the shells). Applied loads: + Pressure of 3 kg/cm² on one of the elements. + Punctual tension force of 50 kp (X=200,Y=0). + Punctual compression force of 100 kp (X=200,Y=100).
There is no bending reinforcement in the slab. The reinforcement checking code is the ACI 318 19.

Element types used in the model: SHELL43

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

ACI 318 19

Model Statistics
Number of elements	200
Number of nodes	231
Number of civil materials	2
Number of cross sections	0
Number of shell vertices	1

Log file: CFVR9174.DAT

   FINISH
  ~cfclear,,1
   NomFile='CFVR9174'
  /TITLE, %NomFile%, Shell shear reinforcement check according to ACI 318 19
! ---------------------------------------------------------------------------------
! Model definition and solve
! ---------------------------------------------------------------------------------

! Setup
 ~UNITS,,MONE,EURO
 ~UNITS,,LENG,CM
 ~UNITS,,TIME,S
 ~UNITS,,FORC,KP
 ~CODESEL,EC3-05,ACI_19 
 ~CFCONFG,PLRES   ,SHELL   ,EPSW 
 ~CFCONFG,RESULT,CFSHF,1
   /PREP7
   
! Materials
  ~CFMP,5,LIB,CONCRETE,EC2,C40/50
  ~CFMP,2,LIB,REINF,EC2,S400
  ~CFMP,5,ACI_14_C,Lambda  ,,0.8  
  
! Element types
  ET,2,SHELL43
  
! Shell Vertex
  t=25     ! Thickness
  r=4      ! Cover
  th = 0   ! Reinforcement angle
 ~SHLRNF,1,t,5,2,r,0.0,0.0,0.0,0.0,0,th,45
 ~SHLSHR,1,2,,12,,,14,14
 
! Beam & Shell property
 ~BMSHPRO,10,SHELL,1,1,1,1,43,,,,Shell property
  
! Elements
  MAT,5
  REAL,10
  TYPE,2
  RECTNG,0,200,0,100, 
  AESIZE,ALL,10,  
  AMESH,All
  
! Boundary conditions
  DL,4, ,ALL,  

! Solve
  /SOLU
  SFE,All,2,PRES,,3.0
  FK,3,FX,-100 
  FK,2,FY,50 
  SOLVE

  /POST1
! ---------------------------------------------------------------------------------
! Check with CivilFEM
! ---------------------------------------------------------------------------------
 ~CFSET,,1
 ~CHKCON,SHELL,SHEAR, , , , , , , , , ,0   ! Phi=.75 ; fy unlimited 
 ~CFSET,1
 
!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp = 11
  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
!--------------------------------------------------------------------------------------
  LABEL( 1,1) = 'VC_X'
  LABEL( 2,1) = 'VS_X'
  LABEL( 3,1) = 'VN_X'
  LABEL( 4,1) = 'VFI_X'
  LABEL( 5,1) = 'CRT_X'
  LABEL( 6,1) = 'VC_Y'
  LABEL( 7,1) = 'VS_Y'
  LABEL( 8,1) = 'VN_Y'
  LABEL( 9,1) = 'VFI_Y'
  LABEL( 10,1) = 'CRT_Y'
  LABEL( 11,1) = 'CRT_TOT'

! Correct values
  VALUE( 1,1)=178.301
  VALUE( 2,1)=247.136
  VALUE( 3,1)=424.8682
  VALUE( 4,1)=318.651134
  VALUE( 5,1)=1.36452
  VALUE( 6,1)=179.817
  VALUE( 7,1)=247.136
  VALUE( 8,1)=426.38426
  VALUE( 9,1)=319.7882
  VALUE(10,1)=0.0
  VALUE(11,1)=1.36452

! Comparison
 ESEL,S,ELEM,,1
 ~PLSHCON,VC_X 
 *GET,VALUE(1,2),NODE,1,EPSW
~PLSHCON,VS_X 
 *GET,VALUE(2,2),NODE,1,EPSW
~PLSHCON,VN_X  
 *GET,VALUE(3,2),NODE,1,EPSW
~PLSHCON,VFI_X 
 *GET,VALUE(4,2),NODE,1,EPSW
~PLSHCON,CRT_X  
 *GET,VALUE(5,2),NODE,1,EPSW
~PLSHCON,VC_Y 
 *GET,VALUE(6,2),NODE,1,EPSW
~PLSHCON,VS_Y 
 *GET,VALUE(7,2),NODE,1,EPSW
~PLSHCON,VN_Y   
 *GET,VALUE(8,2),NODE,1,EPSW
~PLSHCON,VFI_Y 
 *GET,VALUE(9,2),NODE,1,EPSW
~PLSHCON,CRT_Y  
 *GET,VALUE(10,2),NODE,1,EPSW
 *IF,ABS(VALUE(10,2)),LT,1.E-12,THEN
	VALUE(10,2)=0.0
 ~PLSHCON,CRT_TOT  
 *GET,VALUE(11,2),NODE,1,EPSW
 
! Warning and error tolerances
  TOLER(1,1)=  0.01*VALUE(1,1) $ TOLER(1,2)=  0.01*VALUE(1,1)
  TOLER(2,1)= 0.005*VALUE(2,1) $ TOLER(2,2)= 0.005*VALUE(2,1) 
  TOLER(3,1)= 0.005*VALUE(3,1) $ TOLER(3,2)= 0.005*VALUE(3,1) 
  TOLER(4,1)= 0.005*VALUE(4,1) $ TOLER(4,2)= 0.005*VALUE(4,1) 
  TOLER(5,1)= 0.005*VALUE(5,1) $ TOLER(5,2)= 0.005*VALUE(5,1) 
  TOLER(6,1)= 0.005*VALUE(6,1) $ TOLER(6,2)= 0.005*VALUE(6,1) 
  TOLER(7,1)= 0.005*VALUE(7,1) $ TOLER(7,2)= 0.005*VALUE(7,1)
  TOLER(8,1)= 0.005*VALUE(8,1) $ TOLER(8,2)= 0.005*VALUE(8,1) 
  TOLER(9,1)= 0.005*VALUE(9,1) $ TOLER(9,2)= 0.005*VALUE(9,1) 
  TOLER(10,1)= 0.005*VALUE(10,1) $ TOLER(10,2)= 0.005*VALUE(10,1) 
  TOLER(11,1)= 0.005*VALUE(11,1) $ TOLER(11,2)= 0.005*VALUE(11,1) 


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

Results

Label	Target	CivilFEM	Ratio	Tolerance
VC_X	178.3	179.34	0.994	1.783
VS_X	247.14	247.13	1.000	1.236
VN_X	424.87	426.47	0.996	2.124
VFI_X	318.65	319.85	0.996	1.593
CRT_X	1.3645	1.3594	1.004	0.006823
VC_Y	179.82	180.37	0.997	0.8991
VS_Y	247.14	247.13	1.000	1.236
VN_Y	426.38	427.5	0.997	2.132
VFI_Y	319.79	320.63	0.997	1.599
CRT_Y	0	0	0.000	0
CRT_TOT	1.3645	1.3594	1.004	0.006823

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

CivilFEM Online Help

CFVR9174 - Shell shear reinforcement check according to ACI 318 19

Shell shear reinforcement check according to ACI 318 19.

Results