!*************************************************************************************!
! TITLE: CFVR5262 - 2D bending EHE with several tendons
! SUBTITLE: Rectangular section and HA-25
!
! DESCRIPTION: The model is a prestressed fixed-pinned beam 1600 cm length.
! DESCRIPTION:
! DESCRIPTION: The beam is made of HA-25 concrete with tendons that are made of Y1860S3 steel.
! DESCRIPTION:
! DESCRIPTION: The section is box one, with reinforcement of B500S steel.
! DESCRIPTION: It has the following dimensions:
! DESCRIPTION:
! DESCRIPTION: - h = 80 cm
! DESCRIPTION:
! DESCRIPTION: - b = 50 cm
! DESCRIPTION:
! DESCRIPTION: Loads:
! DESCRIPTION:
! DESCRIPTION: - Prestressing loads
! DESCRIPTION:
! DESCRIPTION: - Compression load Fy= -60000 kp
! DESCRIPTION:
! DESCRIPTION: - Surface load of 5 kp
! DESCRIPTION:
! DESCRIPTION: The aim of the example is to verify bending checking for prestressed concrete
! DESCRIPTION: according to EHE.
!
! ELEMENT TYPE: BEAM44
! MODULES: PRSC NLBR
! UNITS: User
! KEYWORD1: Beam
! KEYWORD2: Cross Section
! KEYWORD3: Prestressed Concrete
! KEYWORD4: Code Checking/Design
!
!*************************************************************************************!
  FINISH
  ~CFCLEAR,,1
  ~CFACTIV,NLBR,Y
  ~CFACTIV,PRSC,Y
  NomFile='CFVR5262'
  /TITLE, %NomFile%, 2D bending EHE. Rectangular section

! -------------------------------------------------------------------------------------
! Initial data
! -------------------------------------------------------------------------------------
! Parameters
  DEPTH = 80
  WIDTH = 50
  L = 800
! Setup
  ~CODESEL,EC3-92,EHE-98,EHE-98,,EC8-94
  ~UNITS,,LENG,CM
  ~UNITS,,TIME,S
  ~UNITS,,FORC,KP
  /PREP7
! Materials
  ~CFMP,1,LIB,CONCRETE,EHE,HA-25
  ~CFMP,2,LIB,REINF,EHE,B 500 S
  ~CFMP,3,LIB,PREST,EHE,Y1860S3
! Model definition
  ~CSECDMS,1,REC,1,DEPTH,WIDTH
  ~RNFDEF,1,1,2,4,0,4.0, , ,25.0,5, , ,4
  ~RNFDEF,1,2,2,2,0,4.0, , ,25.0,5, , ,4
  ~BMSHPRO,1,BEAM,1,1,,,4,1,0,,Beam 1
  ET,1,BEAM4
! Nodes
  N, 1,0,0
  N,11,L,0
  FILL
! Elements
  TYPE,1
  REAL,1
  MAT,1
  E,1,2
  EGEN,10,1,-1
! Capture support beam.
  ~SBBMDEF,1
! Elevation points
  ~PCEPDEF,1,1,0,-0.1
  ~PCEPDEF,2,11,0,0.1
  ~PCEPDEF,3,6,-25,0
! Plan points
  ~PCPPDEF,1, 1,-18,0
  ~PCPPDEF,2,11,-18,0
  ~PCPPDEF,3, 1, -6,0
  ~PCPPDEF,4,11, -6,0
  ~PCPPDEF,5, 1,  6,0
  ~PCPPDEF,6,11,  6,0
  ~PCPPDEF,7, 1, 18,0
  ~PCPPDEF,8,11, 18,0
! Tendon parameters
  DiV=5
  AsT=9.8
  AsT=11.6
! Tendon 1
  ~PCTNDEF,1,3,AsT,DiV
  ~PCTNMDF,1,EADD,1.
  ~PCTNMDF,1,EADD,3.
  ~PCTNMDF,1,EADD,2.
  ~PCTNMDF,1,PADD,1.
  ~PCTNMDF,1,PADD,2.
! Tendon 2
  ~PCTNDEF,2,3,AsT,DiV
  ~PCTNMDF,2,EADD,1.
  ~PCTNMDF,2,EADD,3.
  ~PCTNMDF,2,EADD,2.
  ~PCTNMDF,2,PADD,3.
  ~PCTNMDF,2,PADD,4.
! Tendon 3
  ~PCTNDEF,3,3,AsT,DiV
  ~PCTNMDF,3,EADD,1.
  ~PCTNMDF,3,EADD,3.
  ~PCTNMDF,3,EADD,2.
  ~PCTNMDF,3,PADD,5.
  ~PCTNMDF,3,PADD,6.
! Tendon 4
  ~PCTNDEF,4,3,AsT,DiV
  ~PCTNMDF,4,EADD,1.
  ~PCTNMDF,4,EADD,3.
  ~PCTNMDF,4,EADD,2.
  ~PCTNMDF,4,PADD,7.
  ~PCTNMDF,4,PADD,8.
! Stress losses calculation
  ~PCTNMDF,1,METHOD,1
  ~PCTNMDF,2,METHOD,1
  ~PCTNMDF,3,METHOD,1
  ~PCTNMDF,4,METHOD,1
  ~PCLOSS,1E6,0,,0
! Boundary conditions
  D, 1,UX,0
  D, 1,UY,0
  D, 1,UZ,0
  D, 1,ROTX,0
  D, 6,ROTY,0	! se restringe el nodo 6 para que resuelva la version 16.1
  D,11,UY,0

/SOLU
! Solution
  SFBEAM,ALL,2,PRES,5
  F,6,FY,-60000
  ~PCPL,0,
  SOLVE

/POST1
! Postprocessor
  ~CFSET,,1
  ~CHKPRS,2DB, ,0

!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp    = 3
  NComp_ch = 0

! Marix 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)  = 'CRT_TOT1'
  LABEL( 2)  = 'CRT_TOT2'
  LABEL( 3)  = 'CRT_TOT3'

! Correct data
  VALUE(1,1)=0.188
  VALUE(2,1)=0.404
  VALUE(3,1)=0.586

! Obtained data
  ~PLLSPRS,CRT_TOT,1,
   *GET,VALUE(1,2),ELEM,2,ETAB,CFETAB_I
  ~PLLSPRS,CRT_TOT,1,
   *GET,VALUE(2,2),ELEM,4,ETAB,CFETAB_I
  ~PLLSPRS,CRT_TOT,1,
   *GET,VALUE(3,2),ELEM,6,ETAB,CFETAB_I

! 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

!--------------------------------------------------------------------------------------
! Result comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC