!*****************************************************************************************************!
! TITLE: CFVR5226 - Time-dependent losses.Coupled calculation method (EHE)
! SUBTITLE:
!
! DESCRIPTION: The model is a prestressed beam with a length of 1600 cm.
! DESCRIPTION:
! DESCRIPTION: The beam is made of HA-40 concrete and the tendons are made of Y1960S3 steel.
! DESCRIPTION:
! DESCRIPTION: The section is rectangular 100 x 80, with reinforcement of s500 steel.
! DESCRIPTION:
! DESCRIPTION: The aim of the example is to verify the calculation of prestressing losses
! DESCRIPTION: for active tendons.
!
! ELEMENT TYPE: BEAM44
! MODULES: PRSC NLBR
! UNITS: User
! KEYWORD1: Beam
! KEYWORD2: Tendon
! KEYWORD3: Cross Section
! KEYWORD4: Prestressed Concrete
!
!
!*************************************************************************************!
  FINISH
  ~CFCLEAR,,1
  ~CFACTIV,NLBR,Y
  ~CFACTIV,PRSC,Y
  NomFile='CFVR5226'
  /TITLE, %NomFile%, Time-dependent losses. Coupled calculation method (EHE)

! -------------------------------------------------------------------------------------
! Initial data
! -------------------------------------------------------------------------------------
! Parameters
  DEPTH = 100
  WIDTH = 80
  L = 1600
! Setup
  ~CODESEL,EC3-92,EC2-91,EHE,,EC8-94
  ~UNITS,,LENG,CM
  ~UNITS,,TIME,S
  ~UNITS,,FORC,KP
  /PREP7
! Materials
  ~CFMP,1,LIB,CONCRETE,EHE,HA-40
  ~CFMP,2,LIB,REINF,EHE,B 500 S
  ~CFMP,3,LIB,PREST,EHE,Y1960S3
! Model definition
  ~CSECDMS,1,REC,1,DEPTH,WIDTH
  ~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
  ~PCEPDEF,2,11,0,0
! Plan points
  ~PCPPDEF,1,1,0,0
  ~PCPPDEF,2,11,0,0
! Tendons
  ~PCTNDEF,1,3,100,15
  ~PCTNMDF,1,EADD,1.
  ~PCTNMDF,1,EADD,2.
  ~PCTNMDF,1,PADD,1.
  ~PCTNMDF,1,PADD,2.
  ~PCTNMDF,1,METHOD,1
! Stress losses calculation
  ~PCLOSS,1E+06,1,EHE,0.8

!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp    = 9
  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)  = 'C01_SHR'
  LABEL( 2)  = 'C01_CRE'
  LABEL( 3)  = 'C01_RLX'
  LABEL( 4)  = 'C06_SHR'
  LABEL( 5)  = 'C06_CRE'
  LABEL( 6)  = 'C06_RLX'
  LABEL( 7)  = 'C11_SHR'
  LABEL( 8)  = 'C11_CRE'
  LABEL( 9)  = 'C11_RLX'

! Correct data, obtained data
!--------------------------------------------------------------------------------------
  VALUE( 1,1) =  7.750E2/1.1788      $  ~CFGET,value( 1,2),PRSCONC,,TENDON,LOSS,,1, 1,4
  VALUE( 2,1) =  5.611E2/1.1788      $  ~CFGET,value( 2,2),PRSCONC,,TENDON,LOSS,,1, 1,5
  VALUE( 3,1) =  4.910E2/1.1788*0.8  $  ~CFGET,value( 3,2),PRSCONC,,TENDON,LOSS,,1, 1,6
  VALUE( 4,1) =  7.750E2/1.1788      $  ~CFGET,value( 4,2),PRSCONC,,TENDON,LOSS,,1, 6,4
  VALUE( 5,1) =  5.611E2/1.1788      $  ~CFGET,value( 5,2),PRSCONC,,TENDON,LOSS,,1, 6,5
  VALUE( 6,1) =  7.870E2/1.1788*0.8  $  ~CFGET,value( 6,2),PRSCONC,,TENDON,LOSS,,1, 6,6
  VALUE( 7,1) =  7.750E2/1.1788      $  ~CFGET,value( 7,2),PRSCONC,,TENDON,LOSS,,1,11,4
  VALUE( 8,1) =  5.611E2/1.1788      $  ~CFGET,value( 8,2),PRSCONC,,TENDON,LOSS,,1,11,5
  VALUE( 9,1) =  7.220E2/1.1788*0.8  $  ~CFGET,value( 9,2),PRSCONC,,TENDON,LOSS,,1,11,6

! Warning and error tolerances
  *DO,II,1,NComp
    TOLER(II,1)  = 1E1
    TOLER(II,2)  = 1E1
  *ENDDO

  ~PLTEND ! Check Tendon Editor Window

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