Left

CivilFEM Online Help

Right

CFVR0972 - Seismic margin in reinforced shell according to Eurocode 2 (2008)

Shear seismic margin check.

The model is made up of two square shell elements of 1.0m x 1.0m. The nodes of one side are completely restrained and one of the nodes between both element has its vertical movement restrained (perpendicular to the shells).
The initial shear reinforcement are 12 mm bars each 14 cm in both directions.
The applied nonseismic pressure is 1 kg/cm2 on one of the elements. The applied seismic pressure is 3 kg/cm2 on same element.
The aim of this example is to obtain the scale factor for elastic (FSE) and permissible inelastic response (FSI) according to EC2-08 code.
Inelastic energy absorption factors for X and Y direction are 1.2 and 1.3.

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
Eurocode 2

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

Log file: CFVR0972.DAT


   FINISH
  ~CFCLEAR,,1
   NomFile='CFVR0972'
  /TITLE, %NomFile%, Seismic margin in reinforced shell according to Eurocode 2 (2008)
! ---------------------------------------------------------------------------------
! Model definition and solve
! ---------------------------------------------------------------------------------

! Setup
 ~UNITS,,MONE,EURO
 ~UNITS,,LENG,CM
 ~UNITS,,TIME,S
 ~UNITS,,FORC,KP
 ~CODESEL,EC3-92,EC2-08,EC2-08,,EC8-94
 ~CFCONFG,RESULT,CFSHF,0

  /PREP7
! PArameters
 P_NS = 1
 P_S = 3
! Materials
  ~CFMP,5,LIB,CONCRETE,EC2,C40/50
  ~CFMP,2,LIB,REINF,EC2,S400
! Element types
  ET,2,SHELL43

! Shell Vertex
  t=25     ! Thickness
  r=4      ! Cover
 ~SHLRNF,1,t,5,2,r,0.0,0.0,0.0,0.0,0,0,45.0
 ~SHLSHR,1,2,,12,,,14,14
 ~SHLMDF,1,CP,DUC_OP_X,,1.200E+00
 ~SHLMDF,1,CP,DUC_OP_Y,,1.300E+00
 
! Beam & Shell property
 ~BMSHPRO,10,SHELL,1,1,1,1,43,,,,Shell property

! Nodes
  N,10,0,0
  N,20,100,0
  N,30,200,0
  N,40,0,100
  N,50,100,100
  N,60,200,100

! Elements
  MAT,5
  REAL,10
  TYPE,2
  EN,100, 10, 20, 50, 40
  EN,200, 20, 30, 60, 50
  
! Boundary conditions
  D,10,ALL,0.0
  D,40,ALL,0.0
  D,50,UZ,0.0
 
  ~CFSAVE,%NomFile%
  FINISH

! Nonseismic loads
 /FILNAME,CNS
 ~CFRESUM,%NomFile%
 FINISH
 /SOLU
 
 SFE,100,2,PRES,,P_NS
 SOLVE
  
! Seismic loads
 FINISH
 /FILNAME,SPECTRAL
 ~CFRESUM,%NomFile%
 FINISH
 /SOLU
 
 SFE,100,2,PRES,,P_S 
 SOLVE
  
 FINISH
  
 /FILNAME,CNS 	! Point to the ns loads
  

 FINISH
 /POST1
! ---------------------------------------------------------------------------------
! Checking with CivilFEM
! ---------------------------------------------------------------------------------

 ~CFSET,,1
 ~CFFILE2,4,SPECTRAL.RCV !  Specify different RCV file from which to read the load steps for seismic margin check
 ~HCLPFCN,SHELL,SHEAR,,,,,,,,,,,,1,1
 ~CFGET,FSE, ELEMENT, 100 ,RESULT ,  FSE, I, , , 1 
 ~CFSAVE,%NomFile%
 
 FINISH
  
  /FILNAME,SPECTRALSCALED
  ~CFRESUM,%NomFile%
  FINISH
  /SOLU
! Apply Load
  SFE,100,2,PRES,,P_S*FSE
  
  ! Solve
  SOLVE
   FINISH
   
/FILNAME,CNS 
/POST1
! Postprocessor
  ~CFSET,,1,1       ! Read first load step
  ~CFFILE2,4,SPECTRALSCALED.RCV ! Indicamos el sismico
  ~HCLPFCN,SHELL,SHEAR,,,,,,,,,,,,1,1
  ~PLHCLPF,SHELL,FSE
 
 !--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp = 2
  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) = 'DUC_OP_X'
  LABEL(2,1) = 'FSE'

! Correct values
  VALUE(1,1)=min(1.2,1.3)
  VALUE(2,1)=1.

! Comparison
! CivilFEM data
 ~CFGET,VALUE(1,2),SHLVERT,1,CP,DUC_OP_X
 ~CFGET,VALUE(2,2), ELEMENT, 100 ,RESULT ,  FSE, I, , , 1 
 
! Warning and error tolerances
  TOLER( 1, 1)= 1E-02 $ TOLER( 1, 2)= 1E-02
  TOLER( 2, 1)= 5E-02 $ TOLER( 2, 2)= 5E-02
  
!--------------------------------------------------------------------------------------
! Results comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC

Results

LabelTargetCivilFEMRatioTolerance
DUC_OP_X 1.2 1.21.000 0.01
FSE 1 11.000 0.05


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