$â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â–  +PROG TEMPLATE urs:15.3 HEAD Quantities sto#spring_cst 4 $ [kN/m³] spring constant sto#damping .005 $ [kNs/m³] damping sto#sdof_mass 1 $ [kg] sto#sdof_width 0.2 $ [m] sto#sdof_height 0.2 $ [m] sto#sdof_volume #sdof_width^2*#sdof_height $ [m³] sto#sdof_density #sdof_mass/#sdof_volume $ [kg/m³] sto#sdof_unitweight #sdof_density*9.81/10^3 $ [kN/m³] sto#sdof_mat 1 $ material ID of the SDOF mass sto#sdof_cs 1 $ cross section ID of the SDOF mass sto#spring_height 0.5 $ [m] sto#force 1 $ [kN] sto#lc_dl_sw 100 sto#lc_dl 200 sto#dt 0.005 $ [s] time step sto#nr_of_ts 1000 $ number of time steps END $â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â–  +PROG AQUA urs:5.3 HEAD Creation of materials and cross sections UNIT 5 MATE NO #sdof_mat GAM #sdof_unitweight TITL "SDOF" $ Material definition of the SDOF UNIT 5 SREC NO #sdof_cs H #sdof_width*1000 B #sdof_width*1000 MNO #sdof_mat MRF 0 REF C IT 0[o/o] AY 0[o/o] AZ 0[o/o] BCYZ '0' SPT 0 TITL "Cross section of SDOF" END $â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â–  +PROG SOFIMSHA urs:5.2 HEAD Creation of the bridge model SYST type 3D GDIR NEGZ GDIV 5000 NODE NO 1 X 0 Y 0 Z 0 FIX F NODE NO 2 X 0 Y 0 Z #spring_height FIX ZP NODE NO 3 X 0 Y 0 Z #spring_height+#sdof_height FIX ZP TRUS NO 1 NA 2 NE 3 NCS #sdof_cs SPRI NO 1 NA 1 NE 2 DZ 1 CP #spring_cst DP #damping $SPRI NO 2 NA 2 DX 1 CP 1 $ horizontally stabilizing spring for bottom node in x-direction $SPRI NO 3 NA 2 DY 1 CP 1 $ horizontally stabilizing spring for bottom node in y-direction END +PROG SOFILOAD urs:6.2 HEAD Loads on the bridge $__________________________________________________________________________________________ !*! Actions ACT G gamu 1 gamf 1 $ dead load ACT Q $ variable load $__________________________________________________________________________________________ !*! Loadcases LC 1 TYPE G TITL 'Impulse load' NODE 3 TYPE PG P1 #force FUNC T 0,1,1.1,1.2,20 F 0,0,1,0,0 END $â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â–  +PROG ASE urs:9 $ Linear Analysis HEAD Calculation of forces and moments PAGE UNII 0 CTRL OPT WARP VAL 0 LC ALL END $â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â–  +PROG ASE URS:5 HEAD Loadcases in ASE ECHO DISP,REAC,FORC,NOST,BEDD NO $ Create load case which contains the dirac impulse load LC 1 $ this load case is already defined in Sofiload, but has to be defined again in ASE in order to use it for the impulse load $ Create load case which contains the selfweight and the two dead loads: wearing layer and foam LC #lc_dl_sw FACD 1 TITL 'Primary LC: dead weight + SW' $ FACD 1 sets factor of structural dead weight (= selfweight) to 1, because this is 0 by default $ Create load case which contains the two dead loads: wearing layer and foam LC #lc_dl FACD 0 TITL 'Only additional masses' $ this time FACD is set to 0, so selfweight is not taken into account END $â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â–  +PROG ASE URS:6 HEAD Eigenfrequencies [without damping] $MASS LC #lc_dl PRZ 100 EIGE NEIG 2 LC #lc_dl_sw $ LC 5001 $ NEIG determines the number of eigenfrequencies that will be calculated END END $â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â– â–  -PROG DYNA urs:25.1 HEAD Dynamic time step analysis of response to impulse load SYST TYPE REST PHYS NONL CTRL RLC #lc_dyna_first $ RLC: "Result load case for Histories" [source: manual] CTRL STYP MAX $ STYP: "Superposition of results of spectra and stationary response" [source: manual] $ STYP MAX: "The functions are evaluated separately followed by a max/min selection" [source: manual] GRP NO - VAL FULL RADA #alpha RADB #beta $ calculation with Rayleigh damping, determined from predefined damping factor and two natural frequencies CTRL MCON 2 $ MCON: "Formulation of the mass matrix" $ MCON 2: "Consistent mass matrix for system, consistent translatoric element matrices" STEP DT #dt N #nr_of_ts BET 0.25 DEL 0.5 THE 1 $ BET, DEL and THE are parameters related to the integration method, they should be set to 0.25, 0.5 and 1 respectively in order to avoid numerical damping LC #lc_imp MASS LC #lc_dl $CONT NR 101 V 140/3.6 $ Movement of load HIST TYPE AZ LCST 802 HIST TYPE AZ FROM #node_middle LCST 800 $ save time history of the vertical acceleration (AZ) in the middle node to the database HIST TYPE AZ FROM #node_neighbor LCST 801 $ save time history of the vertical acceleration (AZ) in the node neighboring to the middle node to the database END