I am trying to calculate this script of a bridge, with a train load. The script may not be complete, but anyways, I receive a calculation issue for the eigenvalue part. What might be the reason to this?
$SOFiSTiK Export Version 0.1.10
$Automatically generated File of ASCII Export
$ Exported by AQUA Version 19.00-70
PROG AQUA urs:1
HEAD Material and Cross Section
UNIT 5
NORM ‘NS’ ‘en199x-200x-bridge’ COUN 47 CAT ‘D’ SNOW ‘1’ UNIT 5
$Create Materials
CONC 1 C ‘40’ TYPR B TITL “=C 40/50 (EN 1992)”
STEE 2 B ‘500B’ TMAX 32 TITL “=B 500 B (EN 1992)” $Not needed
$MATE NO 101 E 2870 GAM 23.0 $Modulus of Elasticity E=2.87e6 kN/m2 = 2870 MPa; Gamma is unit weight in kN/m2 CTRL
CTRL
CTRL RFCS 1
CTRL FACE -1
CTRL REFD 0
CTRL STYP FEM
CTRL SCUT 17
CTRL PLAS 1
$Create the Cross Section
SECT 1 MNO 1 MRF 2 FSYM YES BTYP BEAM
SV IT 100[o/o]
LAY 1 TYPE MIN MRF 2
LAY 2 TYPE MIN MRF 2
LAY 3 TYPE MIN MRF 2
LAY 4 TYPE MIN MRF 2
POLY TYPE O MNO 1
VERT ‘0100’ Y -3300 Z 0 EXP 1
VERT ‘0101’ Y 3300 Z 0 EXP 1
VERT ‘0102’ Y 3300 Z 350 EXP 1
VERT ‘0103’ Y 1800 Z 350 EXP 1
VERT ‘0104’ Y 1800 Z 2850 EXP 1
VERT ‘0105’ Y 900 Z 2850 EXP 1
VERT ‘0106’ Y 900 Z 350 EXP 1
VERT ‘0107’ Y -900 Z 350 EXP 1
VERT ‘0108’ Y -900 Z 2850 EXP 1
VERT ‘0109’ Y -1800 Z 2850 EXP 1
VERT ‘0110’ Y -1800 Z 350 EXP 1
VERT ‘0111’ Y -3300 Z 350 EXP 1
VERT ‘0100’ Y -3300 Z 0 EXP 1
CUT ‘ZS’ ZB ‘S’ NS 0 MS 0 MNO 1 MRF 2 LAY 1 TYPE WEB VYFK 0 INCL 90
$SVAL NO 1 MNO 101 A 1.0 IY 2.90 TITL “CrossSection”
END
$ Exported by SOFIMSHC Version 23.00-70
PROG SOFIMSHC urs:2
HEAD System
SYST 3d GDIV 10000 GDIR NEGZ
CTRL MESH 1
ctrl node 1
$Create the Structural Points
SPT 1 X 0.0 0.0 0.0 FIX PP
SPT 2 X 15.0 0.0 0.0 FIX XP
SPT 3 X 40.0 0.0 0.0 FIX XP
SPT 4 X 75.0 0.0 0.0 FIX XP
SPT 5 X 95.0 0.0 0.0 FIX XP
SPT 6 X 155.0 0.0 0.0 FIX XP
$Create the Beam 1
SLN 1 NPA 1 2 GRP 1 STYP B SDIV 0.5 SNO 1 DRX 0 0 -1
SLN 2 NPA 2 3 GRP 1 STYP B SDIV 0.5 SNO 1 DRX 0 0 -1
SLN 3 NPA 3 4 GRP 1 STYP B SDIV 0.5 SNO 1 DRX 0 0 -1
SLN 4 NPA 4 5 GRP 1 STYP B SDIV 0.5 SNO 1 DRX 0 0 -1
SLN 5 NPA 5 6 GRP 1 STYP B SDIV 0.5 SNO 1 DRX 0 0 -1
$Create a dummy structure line with number 101 to be selected later in Module Dyna
$to indicate the travel direction
SLN 101 NPA 1 2
END
$ Exported by SOFILOAD Version 17.20-70
$PROG SOFILOAD urs:3
$HEAD EXPORT FROM DATABASE
UNIT TYPE 5
ACT ‘Q’ GAMU 1.500000 0 PSI0 0.700000 0.600000 0.500000 0.800000 PART ‘Q’ SUP COND TITL “øvrige variable laster”
ACT ‘G_1’ GAMU 1.350000 1 PSI0 1 1 1 1 PART ‘G’ SUP PERM TITL “dead load g1”
END
$ Exported by SOFILOAD Version 17.20-70
+PROG SOFILOAD urs:4
HEAD Vehicle
$UNIT TYPE 5
$LC 1 ‘G_1’ 1 FACD 1 TITL “Loadcase 1”
$Create the trains
LC 101 TYPE none titl ‘Loadtrain HSML A1’
TRAI HSLM P1 A3
end
$Start with the Dynamic Analysis
+PROG DYNA urs:5
HEAD Eigen Values of the Bridge
grp 1
EIGE 25 LANC LC 2001
end
+PROG template urs:6
head Definition of Modeling Parameters
LET#V 100.0 $Speed in km/h
sto#V #V/3.6 $Speed in m/s
end
+prog dyna urs:7
head Train Analysis
mass 0 $Refreshes the masses in the system
ctrl mcon 3 $Applies the rotational masses
grp - $Includes all the groups in the analysis
$DAMPING MODEL*
$********************
$Activate for Mass-Prop Damping
grp 1 rada 0.694802 radb 0.000122
$********************
$********************
CTRL RLC 10001
step N 2500 DT 0.001 $NSteps dt MassPropDamp StiffPropDamp
LC 101 $Number of the load case that defines the train loading
CONT NR 101 V #V 0.0 $Automatic timevalues in nodes from beams; CONT Creates the contact and the moving load
$$$Recorder to Dump the Values
HIST A-Z from 9 to 45 inc 6 lcst 3002 dump ‘bridgeAcc__QuartSpan_10Perc.txt’
HIST u-z from 9 to 45 inc 6 lcst 3003 dump ‘bridgeDisp_QuartSpan_10Perc.txt’
end
end