+prog sofiload urs:26.2
head
!*! Cas de charges et raccourcissement des câbles
HEAD Definition of actions and loads
UNIT 5 $ units: sections in mm, geometry+loads in m
ECHO ACT Full
ACT G_1
ACT G_2
ACT B   $ to be able to give a separate safety factor to this cable prestress
ACT P
ACT C

$ Loads used in CSM should get TYPE none because construction stages get an own type in CSM-CS
LC   200 TYPE NONE  TITL 'g1'            ;  beam grp 0 TYPE pg  pa 120
$LC  202 TYPE NONE  TITL 'q'             ;  BEAM GRP 0 TYPE PG  PA 30


let#i 1
loop 116
LC #i TYPE NONE TITL 'cable WX' ; CABL GRP #i TYPE WX PA -10 $ cable shortening 10 mm - scalable loadcases
let#i #i+1
endloop

end

+prog ase urs:26.3
!*! Pré-tension dans les câbles
head  first estimate
$ctrl iter 3 v2 1
syst prob th3 iter -
grp -         yes
sto#prestress 10500
grp (1  29 1)  full  prex #prestress/2
grp (59 87 1)  full  prex #prestress/2
grp (30 40 1)  full  prex #prestress
grp (41 58 1)  full  prex #prestress
grp (88 98 1)  full  prex #prestress
grp (99 116 1) full  prex #prestress


lc 900 facd 1 titl 'first estimate'
   lcc 200
end

+prog ase urs:26.4
!*! Effet de chaînette

head With this frozen stiffness - calculate all scalable loadcases
sfix lc 900     $ frozen stiffness
syst prob line  $ but with frozen stiffness of sfix
LC (1  116 1)
end

+PROG CSM urs:26.5
HEAD
EQLC  900   G;                $ target state g_1 and additional dead load g_2

eqlc  ( 30 58 1);
equu  no    uz
    145     0
    1116    0
    169     0
    1113    0
    1110    0
    1109    0
    1106    0
    1105    0
    1102    0
    1101    0
    185     0
    214     0
    155     0
    216     0
    209     0
    210     0
    223     0
    171     0
    157     0
    211     0
    222     0
    164     0
    181     0
    159     0
    196     0
    178     0
    173     0
    208     0
    200     0


eqlc  ( 88 116 1);
equu        no  uz
            182 0
            175 0
            186 0
            201 0
            149 0
            144 0
            163 0
            194 0
            204 0
            220 0
            203 0
            138 0
            221 0
            160 0
            191 0
            176 0
            202 0
            224 0
            183 0
            150 0
            226 0
            193 0
            205 0
            206 0
            146 0
            151 0
            188 0
            177 0
            190 0

$ superstructure uz positions:
eqlc (1 29 1)
equu no   ux
    1161    0
    1157    0
    1156    0
    1153    0
    1152    0
    1149    0
    1148    0
    1145    0
    1144    0
    1141    0
    1140    0
    1137    0
    1136    0
    1133    0
    1132    0
    1129    0
    1128    0
    1125    0
    1124    0
    1121    0
    1119    0
    1118    0
    1084    0
    1086    0
    1080    0
    1079    0
    1076    0
    1070    0
    1068    0

$ pylon horizontal position:

eqlc (59 87 1)
equu no   ux
    1064    0
    1062    0
    1061    0
    1060    0
    1059    0
    1058    0
    1057    0
    1056    0
    1055    0
    1054    0
    1053    0
    1052    0
    1051    0
    1050    0
    1049    0
    1048    0
    1047    0
    1046    0
    1045    0
    1044    0
    1043    0
    1042    0
    1020    0
    1021    0
    1018    0
    1017    0
    1012    0
    1011    0
    1010    0
END

+prog ase urs:11
head Formfinding first iteration -> LC 910
ctrl warn 256 $ Warning 265: prestress loads act aditionally is correct here
syst prob th3 tol8  -5.0 plc 900
LC 910 FACD 1 titl 'Formfinding first iteration'
  lcc 200
  lcc 'iter' fact 'iter' $ all loadcase and factors from last CSM-Equation run - here loadcases 101:104
end

+PROG CSM URS:3
HEAD 2. CSM-Equation iteration
EQLC cont   $ continue optimizing targets -> set new factors
CTRL ADDF 0 $ do not add factors but use last ASE as PLC in next ASE = best method for cable stressing
            $ only if you use real loads, you must use CTRL ADDF 1 = add factors to last values - and use PLC 91 in next ASE

END

+prog ase urs:13
head Formfinding 2.  final iteration
ctrl warn 265 $ Warning 265: prestress loads act aditionally is correct here
syst prob th3 tol -5.0 plc 910
LC 920 FACD 1 titl 'Formfinding 2. iteration'            $ new loadcase number !
  lcc 200
  lcc 'iter' fact 'iter' $ all loadcase and factors from last CSM-Equation run - here loadcases 101:104
end


+PROG WING urs:26.6
HEAD HEAD Bending moments after first CSM construction stage run
SIZE TYPE WIN SC 1100 SPLI 2x1 ; ctrl warn 873,95
VIEW TYPE DIRE X 0 Y 1 Z 0 AXIS NEGZ
$BOX  XMIN -139.9311 YMIN 12.00012 ZMIN -3.863609 XMAX 133.5521 YMAX -12.00012 ZMAX 59.16440 TYPE BWIN OPER AND BGRP 0
$MOVE X DEL DTYP TEXT ; MOVE X 1 14 UNIT PR SCHH .4 C 1001 DTYP TEXT
   TEXT after first CSM construction stage run - estimated prestress
LC 900 ; grp lc yes ; defo yes 25 ; BEAM N UNIT 10000 SCHH 0.35 STYP cabl ; defo no ; BEAM MY UNIT 20000 SCHH 0.30 ; node uz
LC 910 ; grp lc yes ; defo yes 25 ; BEAM N UNIT 10000 SCHH 0.35 STYP cabl ; defo no ; BEAM MY UNIT 20000 SCHH 0.30 ; node uz
LC 920 ; grp lc yes ; defo yes 25 ; BEAM N UNIT 10000 SCHH 0.35 STYP cabl ; defo no ; BEAM MY UNIT 20000 SCHH 0.30; node uz
END

$ Results are good now - no further iteration necessary

!+!Kapitel Final System Check
+prog ase urs:26.7
head Check LC 930 if everything is in equilibrium in final shape sw+g1

syst post 1  prob th3 tol8 -5 PLC 920 fmax 2
LC 930 FACD 1 titl 'Check final shape -> v=0'
  lcc  200 $ sdl
end

$+prog csm urs:26.8
$head
$cs 1 type G_1
$    grp 1
$end

$-apply "$(NAME)_csm.dat"

$-prog csm urs:26.10
$head
$camb cs 1
$end

$-prog csm urs:26.9
$head
$page line 9999
$echo crep 250
$camb tabn 182,200
$camb tabf 2
$end
$page line 68
$end
$$ Please notice that deformation in the superstructure is very very small !
$-sys del $(project).$d?