+prog aqua urs:1 kopf Material norm en 199X-200X !*!Label C24 Parallel zur Faser holz nr 1 ep 11000.00[N/mm2] g 690[N/mm2] e90 370[N/mm2] g90 50.0[N/mm2] gam 4.2 $$ fm 24.00[N/mm2] fm90 0.00[N/mm2] ft0 14.50[N/mm2] ft90 0.40[N/mm2] fc0 21.00[N/mm2] fc90 2.50[N/mm2] $$ fv 4.00[N/mm2] fvr 2.50[N/mm2] fvb 0.70[N/mm2] !*!Label C24 Senkrecht zur Faser holz nr 2 ep 370[N/mm2] g 50[N/mm2] e90 11000.00[N/mm2] g90 690.0[N/mm2] gam 4.2 $$ fm 0.00[N/mm2] fm90 24.00[N/mm2] ft0 0.40[N/mm2] ft90 14.50[N/mm2] fc0 2.50[N/mm2] fc90 21.00[N/mm2] $$ fv 4.00[N/mm2] fvr 2.50[N/mm2] fvb 0.70[N/mm2] mlay nr 11 t0 40[mm] nr0 1 $$ t1 40[mm] nr1 2 $$ t2 40[mm] nr2 1 $$ t3 40[mm] nr3 2 $$ t4 40[mm] nr4 1 ende $----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- $ Gruppe Nr 1 +PROG SOFIMSHA urs:5 KOPF Deckenplatte SYST 3D GDIR NEGZ GDIV 100000 $ GDIV: Gruppendivisor -> Die Gruppennummer errechnet sich, indem man die Elementnummer ganzzahlig durch denGruppendivisor teilt. !*!Label Eingabeparameter LET#PLATTEN_LAENGE 5.0 $ Länge der Platte in X-Richtung in [m] ; Haupttragrichtung (lang) ; nur gerade Nachkommastellen zulässig (6.0/6.2/6.4 ...) LET#PLATTEN_BREITE 4.0 $ Breite der Platte in X-Richtung in [m] ; Nebentragrichtung (kurz) ; nur gerade Nachkommastellen zulässig (5.0/5.2/5.4 ...) STO#Masse_MFS 0.06 $ Masse für das Masse-Feder-System in [t] LET#Federsteif_MFS 100$13324 $ Federsteifigkeit für das Masse-Feder-System in [kN/m] LET#Daempf_koef_MFS 3.11 $ Dämpfungskoeffizient für den TMD, der bei der Feder berücksichtigt wird in [kN*s/m] LET#Z_KOORD_FEDER_U -0.3[m] $ Länge der Feder LET#Stuetzeneinr_X 0.6 $ Abstand der Stütze vom Rand in X-Richtung in [m] ; nur Abstand mit max. 1 Dezimalstelle zul. (0.2m, 0.5m, 1.2m) ; außer #x1 wird geändert LET#Stuetzeneinr_Y 0.6 $ Abstand der Stütze vom Rand in Y-Richtung in [m] ; nur Abstand mit max. 1 Dezimalstelle zul. (0.2m, 0.5m, 1.2m) ; außer #y1 wird geändert LET#Daempfer_X 2.5 $ X-Koordinate des Dämpfers in [m] LET#Daempfer_Y 2.0 $ Y-Koordinate des Dämpfers in [m] !*!Label Variablen LET#x1 0.1 $ Abstand in x (Netz) in [m] ; am besten nicht verändern oder nur in 10er Potenzen LET#y1 0.1 $ Abstand in y (Netz) in [m] ; am besten nicht verändern oder nur in 10er Potenzen STO#SEG_X #PLATTEN_LAENGE/#x1+1 $ 6.4m ->64 Elemente+1 Knoten STO#SEG_Y #PLATTEN_BREITE/#y1+1 $ 5.4m ->54 Elemente+1 Knoten LET#KNR 1000 $ Knotennummer LET#X_Count 0 LET#Y_Count 0 LET#X_l 0 LET#Y_l 0 !*!Label Generieung Knoten $ Knotennummern: 1000 $ Erstellt erst in X-Richtung eine Reihe mit Knoten im Abstand der x1 Größe (hier im Abstand von 0,1m). Der erste Knoten ist im Punkt (0,0,0) und der letzte Knoten ist bei (Plattenlänge=(SEG_X-1)*0.1m,0,0). Anschließend wird eine weitere Reihe mit Knoten äquivalent um den y1 Abstand (hier 0,1m) erstellt. Der erste Knoten ist im Punkt (0,0.1m,0) und der letzte Knoten ist bei (Plattenlänge=(SEG_X-1)*0.1m,0.1m,0). Dieser Vorgang wird für die Anzahl der Elemente SEG_Y wiederholt. Der allerletzte Konten ist bei (Plattenlänge=(SEG_X-1)*0.1m,Plattenbreite=(SEG_Y-1)*0.1m,0). LOOP #SEG_Y LOOP #SEG_X KNOT NR #KNR X #X_l Y #Y_l Z 0 LET#KNR #KNR+1 LET#X_l #X_l+#x1 ENDLOOP LET#X_l 0 LET#Y_l #Y_l+#y1 ENDLOOP !*!Label Berechnung Mittelknoten STO#KNOT_NR_MITTE 1000+(#SEG_X*#SEG_Y+1)/2-1 $KNOT NR 9999 X #coor_X(#KNOT_NR_MITTE) Y #coor_Y(#KNOT_NR_MITTE) Z #coor_Z(#KNOT_NR_MITTE) PRT#KNOT_NR_MITTE $ Ausgabe der Variablen im Echoprint der Eingabe !*!Label Berechnung Knoten Dämpfer LET#Schritte_D_X #Daempfer_X/#x1 $ Berechnet, wie viele Schritte entlang des FE-Netzes in X-Richtung bis zum Dämpfer gemacht werden müssen LET#Schritte_D_Y #Daempfer_Y/#y1 $ Berechnet, wie viele Schritte entlang des FE-Netzes in Y-Richtung bis zum Dämpfer gemacht werden müssen LET#Daempferknoten 1000+#SEG_X*#Schritte_D_Y+#Schritte_D_X !*!Label Generierung Elemente $ Elementnummern: 100000 $ Erzeugt Flächenelemente durch jeweils 4 der gerade erstellten Knoten. Das erste Element umfasst Knoten 1000 und 1001 sowie 2 der um 0,1m in Y-Richtung verschobenen Punkte. Die Erstellung der Flächenelemente durchläuft zuerst in X-Richtung die Punkte und springt eine Zeile in Y-Richtung und läuft auch dort in X-Richtung durch. LET#X_Count 0 LET#Y_Count 0 LET#X 0 LET#Y 0 LET#KNR 1000 LET#Quad_Nr 100000 LOOP #SEG_Y-1 LOOP #SEG_X-1 quad nr #Quad_Nr k1 #KNR k2 #KNR+1 k3 #KNR+#SEG_X+1 k4 #KNR+#SEG_X mnr 11 nra 7 lage mitt t 20.0000[cm] drx 1.00000000 dry 0.00000000 drz 0.00000000 LET#KNR #KNR+1 LET#Quad_Nr #Quad_Nr+1 ENDLOOP LET#KNR #KNR+1 ENDLOOP !*!Label Generierung Feder $ Elementnummern Feder: 10000 $ Für Eigenwertberechnung diesen Block Feder auskommentieren LET#KNOT_NR_FEDER_O #Daempferknoten $ Knotennummer, an dem die Feder angebracht werden soll: Bsp. in der Mitte -> #KNOT_NR_MITTE oder #Daempferknoten (bei Eingabeparameter X- und Y-Wert festlegen) STO#KNOT_NR_FEDER_U #KNOT_NR_FEDER_O+1000000 $ Knotennummer, am unteren Ende der Feder: Knotennummer an dem Feder angebracht werden soll + 1.000.000 KNOT NR #KNOT_NR_FEDER_U X #coor_X(#KNOT_NR_FEDER_O) Y #coor_Y(#KNOT_NR_FEDER_O) Z #Z_KOORD_FEDER_U GRUP 2 fede NR 10000+#KNOT_NR_FEDER_O KA #KNOT_NR_FEDER_O KE #KNOT_NR_FEDER_U DX 0 DY 0 DZ 1 CP #Federsteif_MFS DP #Daempf_koef_MFS $ Dämpfung in [kN*s/m] als Parameter einführen! GRUP 3 MASS #KNOT_NR_FEDER_U MX #Masse_MFS MY #Masse_MFS MZ #Masse_MFS REF GLOB SELE PG !*!Label Auflagerbedingungen LET#Schritte_einr_X #Stuetzeneinr_X/#x1 $ Berechnet, wie viele Schritte entlang des FE-Netzes in X-Richtung bis zur eingerückten Stütze gemacht werden müssen LET#Schritte_einr_Y #Stuetzeneinr_Y/#y1 $ Berechnet, wie viele Schritte entlang des FE-Netzes in Y-Richtung bis zur eingerückten Stütze gemacht werden müssen LET#Stuetzknoten_1 1000+#SEG_X*#Schritte_einr_Y+#Schritte_einr_X LET#Stuetzknoten_2 1000+#SEG_X*(#Schritte_einr_Y+1)-1-#Schritte_einr_X LET#Stuetzknoten_3 1000+#SEG_X*(#SEG_Y-#Schritte_einr_Y-1)+#Schritte_einr_X LET#Stuetzknoten_4 1000+#SEG_X*(#SEG_Y-#Schritte_einr_Y)-1-#Schritte_einr_X knot nr #Stuetzknoten_1 fix pp $ Festhaltung gegen Verschiebung in x,y und z knot nr #Stuetzknoten_2 fix pp knot nr #Stuetzknoten_3 fix pp knot nr #Stuetzknoten_4 fix pp ende $----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- +prog sofiload urs:4 kopf Ständige Ausbaulast lf 2 $Lastfall 2 quad von grp bis 1 typ PG p 0.15 $ quad: Elementbelastung VON Anfangselement BIS Endelement; VON: Bei VON kann das Literal GRP eingegeben werden dann wird mit BIS eine Gruppennummer komplett ausgewählt. Die Elemente haben Nummer von 100000 bis 103574 und mit Gruppendivisor 100000 werden sie alle Gruppe 1 zugeordnet.; PG: Belastung in Eigengewichtsrichtung; P: Lastwert, hier 4 kN/m2 ende +prog sofiload urs:6 kopf Gehfunktion LF 70 TYP NONE BEZ 'FJMP: Gehlast' $ Es wird eine Einzellast von 2kN an Knoten in der Mitte der Platte aufgebracht, die eine Anregung durch Fußgänger simuliert KNOT #KNOT_NR_MITTE PG 2 FUNK S WALK ende $+prog sofiload urs:8 $kopf Gehfunktion $LF 92 TYP NONE BEZ 'Masse' $ Es wird eine Masse an Knoten in der Mitte der Platte aufgebracht am unteren Ende der Feder $ KNOT #KNOT_NR_FEDER_U typ PG p1 #Masse_MFS $ $ende $----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- +prog ase urs:3 kopf Aufstellung Lastfälle steu warn 240 syst prob line grup 1 voll grup 2 nein grup 3 nein lf alle ende +prog ase urs:10 kopf Aufstellung Lastfälle steu warn 240 syst prob nonl grup 1 voll grup 2 nein grup 3 nein lf 1 egz 1.0 bez 'g1: Eigengewicht Decke' ende +prog ase urs:2 kopf Aufstellung Lastfälle steu warn 240 syst prob nonl grup 1 voll grup 2 voll grup 3 voll lf 10 egz 1.0 bez 'g1: Eigengewicht Decke mit Extramasse' ende +prog ase urs:11 kopf Aufstellung Lastfälle steu warn 240 syst prob nonl grup 1 voll grup 2 nein grup 3 nein GRUP NR 1 MODD 0.025 lf 201 egz 1.35 bez 'uls : 1.35g + 1.5q' lc 2 fakt 1.35 ende +prog ase urs:7 kopf Lastfall EG, Ausbau, Gehlast steu warn 240 syst prob nonl GRUP NR 1 MODD 0.025 grup 2 nein grup 3 nein lf 202 egz 1.0 bez 'EG, Ausbaulast, Gehlast ohne masse ' lc 2 fakt 1.0 lc 70 fakt 1.0 ende +prog ase urs:12 kopf Lastfall EG, Ausbau, Gehlast steu warn 240 syst prob nonl lf 203 egz 1 bez 'EG, Ausbaulast, Gehlast mit masse ' GRUP NR 1 MODD 0.025 Grup 2 voll Grup 3 voll lc 2 fakt 1.0 lc 70 fakt 1.0 ende $----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- +PROG DYNA urs:19.1 KOPF Eigenfrequenzen ECHO ACCE EXTR mass 0 GRUP NR 1 modd 0 grup 2 nein grup 3 nein EIGE neig 50 etyp lanc lf 501; modd d 0.0 ende +PROG DYNA urs:14 KOPF Eigenfrequenzen 1 ECHO ACCE EXTR mass 0 grup nr 1 grup 2 nein grup 3 nein EIGE neig 50 etyp lanc lf 601; modd d 0.025 ende +PROG DYNA urs:15 KOPF Eigenfrequenzen 2 ECHO ACCE EXTR Grup 1 voll grup 2 voll grup 3 voll EIGE neig 50 etyp lanc lf 701; ende +PROG DYNA urs:16 KOPF Dynamische Berechnung Gehfunktion ECHO ACCE EXTR mass 0 LF 202 STEP N 50 DT 0.1 $ Berechnung mittels Newmark-Verfahren !*!Label Beschleunigung HIST TYP A-X VON 1000 BIS 1000+#SEG_X*#SEG_Y ERG DRUC LFSP 71 $ Knotennummer HIST TYP A-Y VON 1000 BIS 1000+#SEG_X*#SEG_Y ERG DRUC LFSP 72 $ Knotennummer HIST TYP A-Z VON 1000 BIS 1000+#SEG_X*#SEG_Y ERG DRUC LFSP 73 $ Knotennummer LF 202 !*!Label Verschiebung HIST TYP U-X VON 1000 BIS 1000+#SEG_X*#SEG_Y ERG DRUC LFSP 91 $ Knotennummer HIST TYP U-Y VON 1000 BIS 1000+#SEG_X*#SEG_Y ERG DRUC LFSP 92 $ Knotennummer HIST TYP U-Z VON 1000 BIS 1000+#SEG_X*#SEG_Y ERG DRUC LFSP 93 $ Knotennummer ENDE +PROG DYNA urs:17 KOPF Dynamische Berechnung Gehfunktion ECHO ACCE EXTR LF 203 STEP N 50 DT 0.1 $ Berechnung mittels Newmark-Verfahren !*!Label Beschleunigung HIST TYP A-X VON 1000 BIS 1000+#SEG_X*#SEG_Y ERG DRUC LFSP 74 $ Knotennummer HIST TYP A-Y VON 1000 BIS 1000+#SEG_X*#SEG_Y ERG DRUC LFSP 75 $ Knotennummer HIST TYP A-Z VON 1000 BIS 1000+#SEG_X*#SEG_Y ERG DRUC LFSP 76 $ Knotennummer ENDE $----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- +PROG DYNR urs:1080 kopf Darstellung Ergebnisse I ECHO VOLL WERT voll !*!Label Diagramm: Zeit-Beschleunigung HIST LFSP 71 TYP A-X NR0 #KNOT_NR_MITTE $ Beschleunigung in X-Richtung im Mittleren Punkt der Platte & dort wo Dämpfer hin sollen HIST LFSP 72 TYP A-Y NR0 #KNOT_NR_MITTE $ Beschleunigung in Y-Richtung im Mittleren Punkt der Platte & dort wo Dämpfer hin sollen HIST LFSP 73 TYP A-Z NR0 #KNOT_NR_MITTE $ Beschleunigung in Z-Richtung im Mittleren Punkt der Platte & dort wo Dämpfer hin sollen HIST LFSP 76 TYP A-Z NR0 #KNOT_NR_MITTE $ Beschleunigung in Z-Richtung im Mittleren Punkt der Platte & dort wo Dämpfer hin sollen !*!Label Diagramm: Zeit-Verschiebung HIST LFSP 91 TYP U-X NR0 #KNOT_NR_MITTE HIST LFSP 92 TYP U-Y NR0 #KNOT_NR_MITTE HIST LFSP 93 TYP U-Z NR0 #KNOT_NR_MITTE !*!Label Diagramm: Frequenz-Beschleunigung ERRE LFSP 73 NR #KNOT_NR_MITTE KOMP Z RESP AA 0 1000 0.01 SPEK LOGF $ Beschleunigung in Z-Richtung im Mittleren Punkt der Platte & dort wo Dämpfer hin sollen !*!Label Diagramm: Frequenz-Beschleunigung ERRE LFSP 76 NR #KNOT_NR_MITTE KOMP Z RESP AA 0 1000 0.01 SPEK LOGF $ Beschleunigung in Z-Richtung im Mittleren Punkt der Platte & dort wo Dämpfer hin sollen ENDE +PROG DYNR urs:9 !*!Label Diagramm: Frequenz-Verschiebung kopf Darstellung Ergebnisse II ECHO VOLL WERT voll ERRE LFSP 73 NR #KNOT_NR_MITTE KOMP Z RESP S 0 1000 0.01 $ Verschiebung in Z-Richtung im Mittleren Punkt der Platte & dort wo Dämpfer hin sollen SPEK LOGF ENDE