CFVR5241 - Primary and secondary forces and moments

T cross section.

The model is a prestressed beam 15 m length.
The beam is made of C16/20 concrete and the tendons are made of Gr250 steel.
The section is a T one with reinforcement of S400 steel. It has the following dimensions:

h = 1 m
tw = 0.5 m
b = 1 m
tf = 0.3 m

Loads:

Self weight
Prestressing loads

The aim of the example is to verify the calculation of primary and secondary forces and moments.
This example is based on "Prestressed concrete, a fundamental approach", 4^th edition, page 16.

Element types used in the model: BEAM44

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

Beam

Tendon

Cross Section

Prestressed Concrete

Prestressing Losses

Model Statistics
Number of elements	8
Number of nodes	9
Number of civil materials	3
Number of cross sections	1
Number of shell vertices	0

Log file: CFVR5242.DAT

   FINISH
   ~CFCLEAR,,1
   ~CFACTIV,PRSC,Y
   NomFile='CFVR5242'
 /TITLE,%NomFile%, I&H CHECKING, ISOST BEAM, T SECTION

! -------------------------------------------------------------------------------------
! Initial data
! -------------------------------------------------------------------------------------
! CivilFEM SETUP
  ~UNITS,SI
  ~CODESEL,,EC2-91,EC2-91
  ~CFACTIV,PRSC,Y

! Preprocessor
  /PREP7

! Material Definition
  ~CFMP,1,LIB,CONCRETE,EC2,C16/20
  ~CFMP,2,LIB,REINF,CEB,S400
  ~CFMP,3,LIB,PREST,ASTMA416,Gr250
! Parameters for not considering prestressing losses.
  ~CFMP,3,Prest,MU,,0
  ~CFMP,3,Prest,K ,,0
  ~CFMP,3,Prest,A ,,0

! Element Type Definition
  ET,1,Beam44

! Cross-section definition
  ~CSECDMS,1,T,1,1 ,0.5,1 ,0.3,0 ,0 ,0 ,0 ,2 ,2
! Beam & Shell property
  ~BMSHPRO,1,BEAM,1,1,,,44,1,1,,Beam 1

! Model Construction
  k,1
  k,2,15
  l,1,2
  dk,1,ux,0,,,uy,uz,rotx
  dk,2,uy,0,,,uz
  h_elems = 4
  medcut = h_elems+1
  endcut   = h_elems*2+1
  lesize,all,,,h_elems*2
  lmesh,all

! Support beam construction
  area = 0.01
  c_diameter = 0.15
  ~SBBMDEF

  ~PCEPDEF,1,1,0,0
  ~PCEPMDF,1,MOVE,CG
  ~PCEPDEF,2,medcut,-0.900,0
  ~PCEPDEF,3,endcut,0,0
  ~PCEPMDF,3,MOVE,CG

  ~PCPPDEF,1,1,0,0
  ~PCPPDEF,3,endcut,0,0

   Papp = 934057.6
  ~PCTNDEF,1,3,area,c_diameter,Papp ,Papp ,0

  ~PCTNMDF,1,EADD,1.
  ~PCTNMDF,1,EADD,2.
  ~PCTNMDF,1,EADD,3.
  ~PCTNMDF,1,PADD,1.
  ~PCTNMDF,1,PADD,3.
  ~PCTNMDF,1,METHOD,1

  ~PCTNMDF,1,EPAROUT,1,0
  ~PCTNMDF,1,EPARIN,2,0

  ~PCLOSS,0,0,0,0.8,

/SOLU
  ~PCPL               ! Prestressing forces
  SOLVE               ! Calculation
  ACEL,,9.81
  SOLVE

/POST1

!--------------------------------------------------------------------------------------
! 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, Correct data, obtained data
!--------------------------------------------------------------------------------------
  ~CFSET,,1
  LABEL( 1) = 'KPRES1'
  VALUE( 1,1) =  1
  ~CFGET,VALUE( 1,2),PRSCONC,,KPRES
  ~CFGET,aux1,ELEMENT,h_elems,FORCE, MZ,J,
  ~CFGET,aux2,ELEMENT,h_elems,FORCE,IMZ,J,
  LABEL( 2) = 'MZ/IMZ'
  VALUE( 2,1) =  1
  VALUE( 2,2)=aux1/aux2
  ~CFGET,aux1,ELEMENT,h_elems,FORCE, FX,J,
  ~CFGET,aux2,ELEMENT,h_elems,FORCE,IFX,J,
  LABEL( 3) = 'FX/IFX'
  VALUE( 3,1) =  1
  VALUE( 3,2)=aux1/aux2
  ~CFGET,aux1,ELEMENT,h_elems,FORCE,HFY,J,
  ~CFGET,aux2,ELEMENT,h_elems,FORCE,IFY,J,
  LABEL( 4) = 'HFY'     $ VALUE( 4,1) =  0.0             $  VALUE( 4,2)=ABS(aux1)
  LABEL( 5) = 'RFX1'    $ VALUE( 5,1) =  0.0             $  *GET,VALUE( 5,2),NODE,1,RF,FX
  LABEL( 6) = 'RFY1'    $ VALUE( 6,1) =  0.0             $  *GET,VALUE( 6,2),NODE,1,RF,FY
  LABEL( 7) = 'RFZ1'    $ VALUE( 7,1) =  0.0             $  *GET,VALUE( 7,2),NODE,1,RF,FZ
  LABEL( 8) = 'RFY2'    $ VALUE( 8,1) =  0.0             $  *GET,VALUE( 8,2),NODE,2,RF,FY
  LABEL( 9) = 'RFZ2'    $ VALUE( 9,1) =  0.0             $  *GET,VALUE( 9,2),NODE,2,RF,FZ
 ~CFSET,,2
  LABEL(10) = 'KPRES2' $  VALUE(10,1) =  1
  ~CFGET,VALUE(10,2),PRSCONC,,KPRES
  ~CFGET,aux1,ELEMENT,h_elems,FORCE,HMZ,J,
  ~CFGET,aux2,ELEMENT,h_elems,FORCE,IMZ,J,
  LABEL(11) = 'HMZ/IMZ'$  VALUE(11,1) = -1
  VALUE(11,2)=aux1/aux2

! Warning and error tolerances
  TOLER(  1,1) = 1E-02   $  TOLER(  1,2) = 1E-02 
  TOLER(  2,1) = 1E-02   $  TOLER(  2,2) = 1E-02 
  TOLER(  3,1) = 1E-02   $  TOLER(  3,2) = 1E-02   
  TOLER(  4,1) = 1E-02   $  TOLER(  4,2) = 1E-02 
  TOLER(  5,1) = 1E-02   $  TOLER(  5,2) = 1E-02 
  TOLER(  6,1) = 1E-02   $  TOLER(  6,2) = 1E-02 
  TOLER(  7,1) = 1E-00   $  TOLER(  7,2) = 1E-00 
  TOLER(  8,1) = 1E-02   $  TOLER(  8,2) = 1E-02 
  TOLER(  9,1) = 1E-02   $  TOLER(  9,2) = 1E-02 
  TOLER( 10,1) = 1E-02   $  TOLER( 10,2) = 1E-02 
  TOLER( 11,1) = 1E-02   $  TOLER( 11,2) = 1E-02 

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

Results

Label	Target	CivilFEM	Ratio	Tolerance
KPRES1	1	1	1.000	0.01
MZ/IMZ	1	1	1.000	0.01
FX/IFX	1	1	1.000	0.01
HFY	0	1.9299e-009	0.000	0.01
RFX1	0	-8.1491e-010	0.000	0.01
RFY1	0	3.638e-010	0.000	0.01
RFZ1	0	0	0.000	1
RFY2	0	-9.8953e-010	0.000	0.01
RFZ2	0	0	0.000	0.01
KPRES2	1	1	1.000	0.01
HMZ/IMZ	-1	-0.9984	1.002	0.01

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

CivilFEM Online Help

CFVR5241 - Primary and secondary forces and moments

T cross section.

Results