Left

CivilFEM Online Help

Official Websites: CivilFEM | CivilFEM Licenses | Ingeciber | Current downloads
CivilFEM for Ansys: Install | Doc | Verif | Instala | Doc(ES) | PDF | PDF(ES) | Support | Soporte | Class3
CivilFEM Apps for Ansys Workbench: Install | Lic | FAQ | Concrete | Steel | Theory | PDF | Support | Soporte | Class3
CivilFEM powered by Marc: Install | Capabilities | Cap(ES) | Doc | Python | Theory | Verif | PDF | Support | Soporte | Class3
Right

CFVR0911 - Shell reinforcement design according to EHE

Shell reinforcement design according to the Spanish code EHE using the CEB-FIP method.

The model is made up of a single squared shell element of 1.0m x 1.0m. The nodes of one side are completely restrained and one of the nodes of the other side has its vertical movement restrained (perpendicular to the shell).
The applied loads are horizontal tension loads (X direction) at the free side and a pressure of 30 kN/m2 on one of the elements.

The reinforcement design method is the CEB-FIP method.

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
EHE

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

Log file: CFVR0911.DAT


   FINISH
  ~CFCLEAR,,1
   NomFile='CFVR0911'
  /TITLE, %NomFile%, Shell reinforcement design according to EHE
! ---------------------------------------------------------------------------------
! Model definition and analysis
! ---------------------------------------------------------------------------------

! Setup
 ~UNITS,SI
 ~CODESEL,,EHE-08
 ~CFCONFG,PLRES   ,SHELL   ,EPSW 
 ~CFCONFG,RESULT,CFSHF,1

/PREP7
! Materials
 ~CFMP,5,LIB,CONCRETE,EHE,HA-35
 ~CFMP,2,LIB,REINF,EHE,B 500 S
! Element types
  ET,2,SHELL43
! Shell vertex
  t=0.12     ! Thickness
  r=0.03     ! Cover
  th = 0     ! Reinforcement angle
 ~SHLRNF,1,t,5,2,r,0.0,0.0,0.0,0.0,0,th,45
! Beam & Shell property
 ~BMSHPRO,10,SHELL,1,1,1,1,43,,,,Shell property

! Nodes
  N,10,0.0,0.0
  N,20,1.0,0.0
  N,40,0.0,1.0
  N,50,1.0,1.0
! Elements
  MAT,5 $ REAL,10 $ TYPE,2
  EN,100, 10, 20, 50, 40
! Boundary conditions
  D,10,ALL,0.0
  D,40,ALL,0.0
  D,50,UZ,0.0
! Solve
  /SOLU
  F,20,FX,30000
  F,50,FX,30000
  SFE,100,2,PRES, ,30000
  SOLVE

/POST1
! ---------------------------------------------------------------------------------
! Forces and moments calculation on nodes of element 100
! ---------------------------------------------------------------------------------
  *DIM,S  ,,6,3,4 ! Stresses(direction,layer,node)
  *DIM,ESF,,8,4   ! Forces & Moments(component,node)   (Tx,Ty,Txy,Mx,My,Mxy,Nx,Ny)
  ESEL,S,ELEM,,100
! Stresses
  SHELL,TOP
  *VGET,S(1,1,1),NODE,10,S,X,,,4
  *VGET,S(1,1,2),NODE,20,S,X,,,4
  *VGET,S(1,1,3),NODE,50,S,X,,,4
  *VGET,S(1,1,4),NODE,40,S,X,,,4
  SHELL,BOT
  *VGET,S(1,3,1),NODE,10,S,X,,,4
  *VGET,S(1,3,2),NODE,20,S,X,,,4
  *VGET,S(1,3,3),NODE,50,S,X,,,4
  *VGET,S(1,3,4),NODE,40,S,X,,,4
! Forces & Moments
  *DO,I,1,4
    ESF(1,I)=t*(S(1,1,I)+S(1,3,I))/2
    ESF(2,I)=t*(S(2,1,I)+S(2,3,I))/2
    ESF(3,I)=t*(S(4,1,I)+S(4,3,I))/2
    ESF(4,I)=-t**2*(S(1,1,I)-S(1,3,I))/12
    ESF(5,I)=-t**2*(S(2,1,I)-S(2,3,I))/12
    ESF(6,I)=-t**2*(S(4,1,I)-S(4,3,I))/12
    ESF(7,I)=t*(S(6,1,I)+S(6,3,I))/2
    ESF(8,I)=t*(S(5,1,I)+S(5,3,I))/2
  *ENDDO

! Designed reinforcement
  *DIM,Axi,,4
  *DIM,Ayi,,4
  *DIM,Axs,,4
  *DIM,Ays,,4

I=1
   Axi(I)= 0.00000000D+00
   Ayi(I)= 0.57603764D-04
   Axs(I)= 0.23519910D-03
   Ays(I)= 0.67195337D-04

I=2
   Axi(I)= 0.12076100D-04
   Ayi(I)= 0.64760552D-04
   Axs(I)= 0.28619516D-03
   Ays(I)= 0.81710711D-04

I=3
   Axi(I)= 0.15731954D-03
   Ayi(I)= 0.10012152D-03
   Axs(I)= 0.15673242D-03
   Ays(I)= 0.62130437D-04

I=4
   Axi(I)= 0.98433127D-04
   Ayi(I)= 0.46154418D-04
   Axs(I)= 0.10573635D-03
   Ays(I)= 0.47615063D-04

  ESEL,ALL

! ---------------------------------------------------------------------------------
! Design with CivilFEM
! ---------------------------------------------------------------------------------
 ~CFSET,,1,1
 ~DIMCON,SHELL,CEB

  ESEL,S,ELEM,,100
! Forces & Moments with CivilFEM
  *DIM,ESF_CF,,8,4   ! Forces & Moments(component,node)   (Tx,Ty,Txy,Mx,My,Mxy,Nx,Ny)
~PLSHFOR,TX
  *GET,ESF_CF(1,1),NODE,10,EPSW
  *GET,ESF_CF(1,2),NODE,20,EPSW
  *GET,ESF_CF(1,3),NODE,50,EPSW
  *GET,ESF_CF(1,4),NODE,40,EPSW
~PLSHFOR,TY
  *GET,ESF_CF(2,1),NODE,10,EPSW
  *GET,ESF_CF(2,2),NODE,20,EPSW
  *GET,ESF_CF(2,3),NODE,50,EPSW
  *GET,ESF_CF(2,4),NODE,40,EPSW
~PLSHFOR,TXY
  *GET,ESF_CF(3,1),NODE,10,EPSW
  *GET,ESF_CF(3,2),NODE,20,EPSW
  *GET,ESF_CF(3,3),NODE,50,EPSW
  *GET,ESF_CF(3,4),NODE,40,EPSW
~PLSHFOR,MX
  *GET,ESF_CF(4,1),NODE,10,EPSW
  *GET,ESF_CF(4,2),NODE,20,EPSW
  *GET,ESF_CF(4,3),NODE,50,EPSW
  *GET,ESF_CF(4,4),NODE,40,EPSW
~PLSHFOR,MY
  *GET,ESF_CF(5,1),NODE,10,EPSW
  *GET,ESF_CF(5,2),NODE,20,EPSW
  *GET,ESF_CF(5,3),NODE,50,EPSW
  *GET,ESF_CF(5,4),NODE,40,EPSW
~PLSHFOR,MXY
  *GET,ESF_CF(6,1),NODE,10,EPSW
  *GET,ESF_CF(6,2),NODE,20,EPSW
  *GET,ESF_CF(6,3),NODE,50,EPSW
  *GET,ESF_CF(6,4),NODE,40,EPSW
~PLSHFOR,NX
  *GET,ESF_CF(7,1),NODE,10,EPSW
  *GET,ESF_CF(7,2),NODE,20,EPSW
  *GET,ESF_CF(7,3),NODE,50,EPSW
  *GET,ESF_CF(7,4),NODE,40,EPSW
~PLSHFOR,NY
  *GET,ESF_CF(8,1),NODE,10,EPSW
  *GET,ESF_CF(8,2),NODE,20,EPSW
  *GET,ESF_CF(8,3),NODE,50,EPSW
  *GET,ESF_CF(8,4),NODE,40,EPSW

! Designed reinforcement obtained with CivilFEM
  *DIM,Axi_CF,,4
  *DIM,Ayi_CF,,4
  *DIM,Axs_CF,,4
  *DIM,Ays_CF,,4
~PLSHCON,ASBX
  *GET,Axi_CF(1),NODE,10,EPSW
  *GET,Axi_CF(2),NODE,20,EPSW
  *GET,Axi_CF(3),NODE,50,EPSW
  *GET,Axi_CF(4),NODE,40,EPSW
~PLSHCON,ASBY
  *GET,Ayi_CF(1),NODE,10,EPSW
  *GET,Ayi_CF(2),NODE,20,EPSW
  *GET,Ayi_CF(3),NODE,50,EPSW
  *GET,Ayi_CF(4),NODE,40,EPSW
~PLSHCON,ASTX
  *GET,Axs_CF(1),NODE,10,EPSW
  *GET,Axs_CF(2),NODE,20,EPSW
  *GET,Axs_CF(3),NODE,50,EPSW
  *GET,Axs_CF(4),NODE,40,EPSW
~PLSHCON,ASTY
  *GET,Ays_CF(1),NODE,10,EPSW
  *GET,Ays_CF(2),NODE,20,EPSW
  *GET,Ays_CF(3),NODE,50,EPSW
  *GET,Ays_CF(4),NODE,40,EPSW

!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp = 48
  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) ='TX_I'
  LABEL( 2) ='TX_J'
  LABEL( 3) ='TX_K'
  LABEL( 4) ='TX_L'
  LABEL( 5) ='TY_I'
  LABEL( 6) ='TY_J'
  LABEL( 7) ='TY_K'
  LABEL( 8) ='TY_L'
  LABEL( 9) ='TXY_I'
  LABEL(10) ='TXY_J'
  LABEL(11) ='TXY_K'
  LABEL(12) ='TXY_L'
  LABEL(13) ='MX_I'
  LABEL(14) ='MX_J'
  LABEL(15) ='MX_K'
  LABEL(16) ='MX_L'
  LABEL(17) ='MY_I'
  LABEL(18) ='MY_J'
  LABEL(19) ='MY_K'
  LABEL(20) ='MY_L'
  LABEL(21) ='MXY_I'
  LABEL(22) ='MXY_J'
  LABEL(23) ='MXY_K'
  LABEL(24) ='MXY_L'
  LABEL(25) ='NX_I'
  LABEL(26) ='NX_J'
  LABEL(27) ='NX_K'
  LABEL(28) ='NX_L'
  LABEL(29) ='NY_I'
  LABEL(30) ='NY_J'
  LABEL(31) ='NY_K'
  LABEL(32) ='NY_L'
  LABEL(33) ='ASxi_I'
  LABEL(34) ='ASxi_J'
  LABEL(35) ='ASxi_K'
  LABEL(36) ='ASxi_L'
  LABEL(37) ='ASyi_I'
  LABEL(38) ='ASyi_J'
  LABEL(39) ='ASyi_K'
  LABEL(40) ='ASyi_L'
  LABEL(41) ='ASxs_I'
  LABEL(42) ='ASxs_J'
  LABEL(43) ='ASxs_K'
  LABEL(44) ='ASxs_L'
  LABEL(45) ='ASys_I'
  LABEL(46) ='ASys_J'
  LABEL(47) ='ASys_K'
  LABEL(48) ='ASys_L'

! Correct values
! ---------------------------------------------------------------------------------
  ! Forces & Moments
    *DO,I,1,8
      *DO,J,1,4
        VALUE(4*(I-1)+J,1)= ESF(I,J)
      *ENDDO
    *ENDDO
  ! Reinforcement
    *DO,I,1,4
      VALUE(32+I,1)=Axi(I)
    *ENDDO
    *DO,I,1,4
      VALUE(36+I,1)=Ayi(I)
    *ENDDO
    *DO,I,1,4
      VALUE(40+I,1)=Axs(I)
    *ENDDO
    *DO,I,1,4
      VALUE(44+I,1)=Ays(I)
    *ENDDO

! Obtained values
! ---------------------------------------------------------------------------------
  ! Forces & Moments
    *DO,I,1,8
      *DO,J,1,4
        VALUE(4*(I-1)+J,2)= ESF_CF(I,J)
      *ENDDO
    *ENDDO
  ! Reinforcement
    *DO,I,1,4
      VALUE(32+I,2)=Axi_CF(I)
    *ENDDO
    *DO,I,1,4
      VALUE(36+I,2)=Ayi_CF(I)
    *ENDDO
    *DO,I,1,4
      VALUE(40+I,2)=Axs_CF(I)
    *ENDDO
    *DO,I,1,4
      VALUE(44+I,2)=Ays_CF(I)
    *ENDDO

! Warning and error tolerances
  TOLER( 1, 1)= 1E-01   $   TOLER( 1, 2)= 1E-00
  TOLER( 2, 1)= 1E-01   $   TOLER( 2, 2)= 1E-00
  TOLER( 3, 1)= 1E-02   $   TOLER( 3, 2)= 1E-01
  TOLER( 4, 1)= 1E-02   $   TOLER( 4, 2)= 1E-01
  TOLER( 5, 1)= 1E-02   $   TOLER( 5, 2)= 1E-01
  TOLER( 6, 1)= 1E-02   $   TOLER( 6, 2)= 1E-01
  TOLER( 7, 1)= 1E-02   $   TOLER( 7, 2)= 1E-01
  TOLER( 8, 1)= 1E-02   $   TOLER( 8, 2)= 1E-01
  TOLER( 9, 1)= 1E-11   $   TOLER( 9, 2)= 1E-11
  TOLER(10, 1)= 1E-11   $   TOLER(10, 2)= 1E-11
  TOLER(11, 1)= 1E-11   $   TOLER(11, 2)= 1E-11
  TOLER(12, 1)= 1E-11   $   TOLER(12, 2)= 1E-11
  TOLER(13, 1)= 1E-02   $   TOLER(13, 2)= 1E-01
  TOLER(14, 1)= 1E-03   $   TOLER(14, 2)= 1E-02
  TOLER(15, 1)= 1E-02   $   TOLER(15, 2)= 1E-01
  TOLER(16, 1)= 1E-03   $   TOLER(16, 2)= 1E-02
  TOLER(17, 1)= 1E-02   $   TOLER(17, 2)= 1E-01
  TOLER(18, 1)= 1E-04   $   TOLER(18, 2)= 1E-03
  TOLER(19, 1)= 1E-04   $   TOLER(19, 2)= 1E-03
  TOLER(20, 1)= 1E-03   $   TOLER(20, 2)= 1E-02
  TOLER(21, 1)= 1E-03   $   TOLER(21, 2)= 1E-02
  TOLER(22, 1)= 1E-03   $   TOLER(22, 2)= 1E-02
  TOLER(23, 1)= 1E-03   $   TOLER(23, 2)= 1E-02
  TOLER(24, 1)= 1E-03   $   TOLER(24, 2)= 1E-02
  TOLER(25, 1)= 1E-03   $   TOLER(25, 2)= 1E-02
  TOLER(26, 1)= 1E-03   $   TOLER(26, 2)= 1E-02
  TOLER(27, 1)= 1E-03   $   TOLER(27, 2)= 1E-02
  TOLER(28, 1)= 1E-03   $   TOLER(28, 2)= 1E-02
  TOLER(29, 1)= 1E-11   $   TOLER(29, 2)= 1E-11
  TOLER(30, 1)= 1E-03   $   TOLER(30, 2)= 1E-02
  TOLER(31, 1)= 1E-03   $   TOLER(31, 2)= 1E-02
  TOLER(32, 1)= 1E-13   $   TOLER(32, 2)= 1E-13
  TOLER(33, 1)= 1E-06   $   TOLER(33, 2)= 1E-05
  TOLER(34, 1)= 1E-06   $   TOLER(34, 2)= 1E-05
  TOLER(35, 1)= 1E-06   $   TOLER(35, 2)= 1E-05
  TOLER(36, 1)= 1E-06   $   TOLER(36, 2)= 1E-05
  TOLER(37, 1)= 1E-06   $   TOLER(37, 2)= 1E-05
  TOLER(38, 1)= 1E-06   $   TOLER(38, 2)= 1E-05
  TOLER(39, 1)= 1E-06   $   TOLER(39, 2)= 1E-05
  TOLER(40, 1)= 1E-06   $   TOLER(40, 2)= 1E-05
  TOLER(41, 1)= 1E-06   $   TOLER(41, 2)= 1E-05
  TOLER(42, 1)= 1E-06   $   TOLER(42, 2)= 1E-05
  TOLER(43, 1)= 1E-06   $   TOLER(43, 2)= 1E-05
  TOLER(44, 1)= 1E-06   $   TOLER(44, 2)= 1E-05
  TOLER(45, 1)= 1E-06   $   TOLER(45, 2)= 1E-05
  TOLER(46, 1)= 1E-06   $   TOLER(46, 2)= 1E-05
  TOLER(47, 1)= 1E-06   $   TOLER(47, 2)= 1E-05
  TOLER(48, 1)= 1E-06   $   TOLER(48, 2)= 1E-05

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

Results

LabelTargetCivilFEMRatioTolerance
TX_I 60000 600001.000 1
TX_J 60000 600001.000 1
TX_K 60000 600001.000 0.1
TX_L 60000 600001.000 0.1
TY_I 12000 120001.000 0.1
TY_J -6000 -60001.000 0.1
TY_K -6000 -60001.000 0.1
TY_L 12000 120001.000 0.1
TXY_I6.9849e-0126.9849e-0121.000 1e-011
TXY_J-1.397e-011-1.397e-0111.000 1e-011
TXY_K2.794e-0112.794e-0111.000 1e-011
TXY_L-1.7462e-011-1.7462e-0111.000 1e-011
MX_I -4904.5 -4904.51.000 0.1
MX_J -4767.3 -4767.31.000 0.01
MX_K 10.211 10.2111.000 0.1
MX_L -127.01 -127.011.000 0.01
MY_I -980.9 -980.91.000 0.1
MY_J -294.79 -294.791.000 0.001
MY_K 660.71 660.711.000 0.001
MY_L -25.403 -25.4031.000 0.01
MXY_I 876.33 876.331.000 0.01
MXY_J 2787.3 2787.31.000 0.01
MXY_K 3061.8 3061.81.000 0.01
MXY_L 1150.8 1150.81.000 0.01
NX_I -23397 -233971.000 0.01
NX_J -23397 -233971.000 0.01
NX_K 13608 136081.000 0.01
NX_L 13608 136081.000 0.01
NY_I6.3301e-0126.3301e-0121.000 1e-011
NY_J 5907.1 5907.11.000 0.01
NY_K 5907.1 5907.11.000 0.01
NY_L6.3301e-0126.3301e-0121.000 1e-013
ASxi_I 0 00.000 1e-005
ASxi_J1.2076e-0051.2076e-0051.000 1e-005
ASxi_K0.000157320.000157321.000 1e-005
ASxi_L9.8433e-0059.8433e-0051.000 1e-005
ASyi_I5.7604e-0055.7604e-0051.000 1e-005
ASyi_J6.4761e-0056.4761e-0051.000 1e-005
ASyi_K0.000100120.000100121.000 1e-005
ASyi_L4.6154e-0054.6154e-0051.000 1e-005
ASxs_I 0.0002352 0.00023521.000 1e-005
ASxs_J 0.0002862 0.00028621.000 1e-005
ASxs_K0.000156730.000156731.000 1e-005
ASxs_L0.000105740.000105741.000 1e-005
ASys_I6.7195e-0056.7195e-0051.000 1e-005
ASys_J8.1711e-0058.1711e-0051.000 1e-005
ASys_K6.213e-0056.213e-0051.000 1e-005
ASys_L4.7615e-0054.7615e-0051.000 1e-005

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