!-----------------------------------------------
! RELY - Reliability analysis
!-----------------------------------------------
+PROG RELY urs:7
HEAD Reliability analysis


! Define random variables
!(TYPE = distibution type, PTYP = distribution defined in terms of moments or parameter,
! P1 = mean, P2 = standard deviation)

    ! Strength
    VAR 'FM'   TYPE LOGN     PTYP MOME P1 37.09[MPa] P2 9.27[MPa]  ! Bending strength               [MPa]
    VAR 'FC0'  TYPE LOGN     PTYP MOME P1 33.90[MPa] P2 6.78[MPa]  ! Compression strength           [MPa]
    VAR 'FT0'  TYPE LOGN     PTYP MOME P1 27.07[MPa] P2 8.12[MPa]  ! Tension strength               [MPa]
    VAR 'FVR'  TYPE LOGN     PTYP MOME P1 1.17[MPa]  P2 0.26[MPa]  ! Shear strength (rolling)       [MPa]

    ! Stiffness
    VAR 'E1'  TYPE LOGN      PTYP MOME P1 11600      P2  1508      ! Modulus of elasticity(para.)  [MPa]
    VAR 'E2'  TYPE LOGN      PTYP MOME P1 300        P2  39.0      ! Modulus of elasticity(perp.)  [MPa]
    VAR 'G1'  TYPE LOGN      PTYP MOME P1 650        P2  84.5      ! Shear moduls (para.)          [MPa]
    VAR 'G2'  TYPE LOGN      PTYP MOME P1 65         P2  16.3      ! Shear moduls (perp.)          [MPa]
    VAR 'CM'  TYPE LOGN      PTYP MOME P1 248        P2  24.8      ! Rot. spring constant          [kNm/rad/m]


    ! Model
    VAR 'X_R'  TYPE LOGN     PTYP MOME P1 1.00       P2 0.05       ! Model uncertainty              [-]

    ! Load
    VAR 'g_k1' TYPE NORM     PTYP MOME P1 1.37       P2 0.14       ! Dead load 1                    [kN/m2]
    VAR 'g_k2' TYPE NORM     PTYP MOME P1 1.26       P2 0.13       ! Dead load 2                    [kN/m2]
    VAR 'g_k3' TYPE NORM     PTYP MOME P1 1.28       P2 0.13       ! Dead load 3                    [kN/m2]
    VAR 'q_k'  TYPE GUMBMAX  PTYP MOME P1 0.988      P2 0.198      ! Live load                      [kN/m2]
    VAR 'w_k'  TYPE GUMBMAX  PTYP MOME P1 0.461      P2 0.185      ! Wind load                      [kN/m2]



! Define correlation between variables

    CORR VAR1 'FM'  VAR2 'FC0'  C 0.49  ! Original 0.8
    CORR VAR1 'FM'  VAR2 'FT0'  C 0.49  ! Original 0.8
    CORR VAR1 'FM'  VAR2 'FVR'  C 0.3
    CORR VAR1 'FM'  VAR2 'E1'   C 0.49  ! Original 0.8
    CORR VAR1 'FM'  VAR2 'E2'   C 0.49  ! Original 0.6
    CORR VAR1 'FM'  VAR2 'G1'   C 0.4
    CORR VAR1 'FM'  VAR2 'G2'   C 0.2
    CORR VAR1 'FT0' VAR2 'FC0'  C 0.49  ! Original 0.5
    CORR VAR1 'FT0' VAR2 'FVR'  C 0.4
    CORR VAR1 'FT0' VAR2 'E1'   C 0.49  ! Original 0.6
    CORR VAR1 'FT0' VAR2 'E2'   C 0.2
    CORR VAR1 'FT0' VAR2 'G1'   C 0.4
    CORR VAR1 'FT0' VAR2 'G2'   C 0.2
    CORR VAR1 'FC0' VAR2 'FVR'  C 0.2
    CORR VAR1 'FC0' VAR2 'E1'   C 0.49  ! Original 0.6
    CORR VAR1 'FC0' VAR2 'E2'   C 0.49
    CORR VAR1 'FC0' VAR2 'G1'   C 0.4
    CORR VAR1 'FC0' VAR2 'G2'   C 0.2
    CORR VAR1 'FVR' VAR2 'E1'   C 0.2
    CORR VAR1 'FVR' VAR2 'E2'   C 0.4
    CORR VAR1 'FVR' VAR2 'G1'   C 0.4
    CORR VAR1 'FVR' VAR2 'G2'   C 0.4
    CORR VAR1 'E1'  VAR2 'E2'   C 0.4
    CORR VAR1 'E1'  VAR2 'G1'   C 0.49  ! Original 0.6
    CORR VAR1 'E1'  VAR2 'G2'   C 0.4
    CORR VAR1 'E2'  VAR2 'G1'   C 0.49  ! Original 0.6
    CORR VAR1 'E2'  VAR2 'G2'   C 0.4
    CORR VAR1 'G1'  VAR2 'G2'   C 0.4

! Define a limit-state function

!---------- ! ULS - ULTIMATE LIMIT STATE !-----------------!

    ! Slab
   FUNC 'LS1'    ! Limit state function 1: Bending moment
   ! FUNC 'LS2'    ! Limit state function 2: Rolling shear

    ! Interior wall: wall 1
   ! FUNC 'LS3'    ! Limit state function 3: Combined compression and bending

    ! Exterior walls: wall 2
   ! FUNC 'LS4'    ! Limit state function 4: Combined compression and bending

    ! Exterior wall: wall 3
    $FUNC 'LS5'    ! Limit state function 5: Combined compression and bending
    $FUNC 'LS6'    ! Limit state function 6: Combined tension and bending

!---------- ! SLS - SERVICE LIMIT STATE !-----------------!

    ! Slab
    $FUNC 'LS7'    ! Limit state function 7: Vertical displacement failure

    ! Entire building
    $FUNC 'LS8'    ! Limit state function 8: Horisontal displacemnet


! Define the SOFiSTiK project that evaluates the limit-state function
! DAT = Name of the DAT-file of the project to call

    PROJ DAT 'TypoA_TEDDY'

! FORM – First order reliability method

    FORM

! CTRL – Control switches (SEED = ID of a seed value to use)

    CTRL SEED 1

END