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: |
|
| 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
| Label | Target | CivilFEM | Ratio | Tolerance |
| DUC_OP_X | 1.2 | 1.2 | 1.000 | 0.01 |
| FSE | 1 | 1 | 1.000 | 0.05 |
Contains proprietary and confidential information of Ingeciber, S.A.

