+SYS DEL $(projekt).cdb $$$$$$$$$$$$$$$$$$$$$$$$$ !#!Kapitel Eingangsgrößen $$$$$$$$$$$$$$$$$$$$$$$$$ !*!Label Variablendefinition $ Geometrie #define l=1.0 $ Laenge Pruefkoerper #define lm=0.5 $ Halbe Laenge Pruefkoerper #define b=0.6 $ Breite Pruefkoerper #define bm=0.3 $ Halbe Breite Pruefkoerper #define tu=0.3 $ Dicke Pruefkoerper unten #define to=0.3 $ Dicke Pruefkoerper oben #define df=0.01 $ Dicke Folie #define db=0.02 $ Durchmesser Duebel #define ds1=0.0025 $ Schlupf zwischen Duebel und Beton-RB #define ds2=0.001 $ Horizontaler Schlupf #define dl1=0.15 $ Duebellaenge oben #define dl2=0.15 $ Duebellaenge unten #define A=($(l)*$(to)) $ Flaeche Querschnitt #define ys=($(lm)) $ Flaechenschwerpunkt in Y #define zs=($(to)/2) $ Flaechenschwerpunkt in Z $ Netzgenerierung #define gdiv=10000 $ Gruppendivisor !*!Label Materialien $ Material $ Festigkeitseigenschaften Randbalken #define EC1=31000 #define E28=31000 #define aT=10.0*10**(-6) $ Temperaturausdehnungskoeffizient #define fct1=3.00 $ 28-Tage Zugfestigkeit #define fctm=3.00 $ mittlere Zugfestigkeit #define mue=0.2 $ Exponent nach MC90; Querdehnzahl $ Festigkeitseigenschaften Bruecke #define EC1_alt=31000 #define E28_alt=31000 $$$$$$$$$$$$$$$$$$$$$$$$$$$$ !#!Kapitel Systemgenerierung $$$$$$$$$$$$$$$$$$$$$$$$$$$$ +PROG AQUA urs:1 KOPF Materialeigenschaften !*!Label Materialbezeichnungen NORM EN 1992-2004 UNIT 0 $ Norm und Einheiten MATE 1 E $(E28) ALFA ($(at)) MUE $(mue) BEZ 'Randbalken' MATE 2 E $(E28_alt) ALFA ($(at)) MUE $(mue) gam 0 BEZ 'Bruecke' MATE 3 E 1 BEZ 'Abdichtung' MATE 4 TYP PS BEZ 'Daemmung' STAH 5 ART S GUET 450 BEZ 'Duebel' STAH 6 ART B GUET 550B BEZ 'Bewehrung' BETO 7 ART C FCN 50 BEZ 'Beton-Versuch' ENDE +PROG SOFIMSHC urs:2 KOPF Systemgenerierung SYST RAUM GDIR POSZ GDIV $(gdiv) STEU MESH 2 STEU HMIN 0.5 STEU FEIN 0.25 STEU TOLN 0.00001 let#r1 $(db)/2 $ Radius Duebel let#r2 $(db)/2+$(ds1) $ Radius Beton let#dxy1 #r1/(2**(1/2)) let#dxy2 #r2/(2**(1/2)) sto#r1 sto#r2 sto#dxy1 sto#dxy2 PRT#r1 PRT#r2 PRT#dxy1 PRT#dxy2 !*!Label Strukturpunkte $ Strukturpunkte für Flaechenelemente als Liste SPT NR X Y Z 1 0 0 $(to) 2 $(bm) 0 $(to) 3 $(bm) ($(lm)-#r2) $(to) 4 ($(bm)-#dxy2) ($(lm)-#dxy2) $(to) 5 ($(bm)-#r2) $(lm) $(to) 6 ($(bm)-#dxy2) ($(lm)+#dxy2) $(to) 7 $(bm) ($(lm)+#r2) $(to) 8 $(bm) $(l) $(to) 9 0 $(l) $(to) 11 $(bm) ($(lm)-#r1) $(to) 12 ($(bm)-#dxy1) ($(lm)-#dxy1) $(to) 13 ($(bm)-#r1) $(lm) $(to) 14 ($(bm)-#dxy1) ($(lm)+#dxy1) $(to) 15 $(bm) ($(lm)+#r1) $(to) 21 $(bm) ($(lm)-#r2) 0 22 ($(bm)-#dxy2) ($(lm)-#dxy2) 0 23 ($(bm)-#r2) ($(lm)) 0 24 ($(bm)-#dxy2) ($(lm)+#dxy2) 0 25 ($(bm)) ($(lm)+#r2) 0 !*!Label Strukturlinien $ Strukturlinien für die Flaechenelemente SLN NR NPA NPE 12 1 2 23 2 3 34 3 4 45 4 5 56 5 6 67 6 7 78 7 8 89 8 9 91 9 1 1112 11 12 1213 12 13 1314 13 14 1415 14 15 1511 15 11 2122 21 22 2223 22 23 2324 23 24 2425 24 25 2521 25 21 !*!Label Strukturflaechen $ Strukturflaechen für Flaechenelemente SAR 1 MNR 7 NRA 0 GRP 10 T 0 SARB NL 12 SARB NL 23 SARB NL 34 SARB NL 45 SARB NL 56 SARB NL 67 SARB NL 78 SARB NL 89 SARB NL 91 SAR 2 MNR 5 NRA 0 GRP 11 T 0 SARB NL 1112 SARB NL 1213 SARB NL 1314 SARB NL 1415 SARB NL 1511 SAR 3 MNR 7 NRA 0 GRP 12 T 0 SARB NL 2122 SARB NL 2223 SARB NL 2324 SARB NL 2425 SARB NL 2521 ENDE +PROG SOFIMSHA urs:3 KOPF Extrusion ECHO voll voll SYST REST GDIR POSZ GDIV $(gdiv) STEU gtol 0.001 $ Grundflaeche durch Spiegelung erzeugen GRUP 10 SPIE TYP QUAD VON 10 DELT GRP A 1 B 0 C 0 DIST $(bm) GRUP 11 SPIE TYP QUAD VON 11 DELT GRP A 1 B 0 C 0 DIST $(bm) GRUP 12 SPIE TYP QUAD VON 12 DELT GRP A 1 B 0 C 0 DIST $(bm) !*!Label Generierung der Volumenelemente $ Erzeugung von Volumenelementen - OBEN BRIC PROP MNR 7 GRUP 1 EXTR TYP QUAD VON GRP 10 PATH DIR 0 0 (-$(to)) DIV 6 BRIC PROP MNR 7 GRUP 1 EXTR TYP QUAD VON GRP 12 PATH DIR 0 0 ($(to)-$(dl1)-$(ds2)) DIV 3 $ Erzeugung der Volumenelemente - UNTEN GRUP 20 TRAN TYP QUAD VON 10 DELT GRP DZ ($(df)) BRIC PROP MNR 7 GRUP 2 EXTR TYP QUAD VON GRP 20 PATH DIR 0 0 ($(tu)) DIV 6 GRUP 22 TRAN TYP QUAD VON 12 DELT GRP DZ ($(to)+$(df)+$(tu)) BRIC PROP MNR 7 GRUP 2 EXTR TYP QUAD VON GRP 22 PATH DIR 0 0 -($(tu)-$(dl2)-$(ds2)) DIV 3 $ Erzeugung der Volumenelemente - DUEBEL BRIC PROP MNR 5 GRUP 3 EXTR TYP QUAD VON GRP 11 PATH DIR 0 0 (-$(dl1)) DIV 3 EXTR TYP QUAD VON GRP 11 PATH DIR 0 0 ($(df)) DIV 2 GRUP 21 TRAN TYP QUAD VON 11 DELT GRP DZ ($(df)) BRIC PROP MNR 5 GRUP 3 EXTR TYP QUAD VON GRP 21 PATH DIR 0 0 ($(dl2)) DIV 3 !*!Label Generierung der Federelemente $ Federkopplung zwischen Duebel und Beton - OBEN QUAD PROP MNR 0 NRA 0 T 0 GRUP 31 EXTR TYP SLN VON SLN 34 PATH DIR 0 0 -$(dl1) DIV 3 EXTR TYP SLN VON SLN 45 PATH DIR 0 0 -$(dl1) DIV 3 EXTR TYP SLN VON SLN 56 PATH DIR 0 0 -$(dl1) DIV 3 EXTR TYP SLN VON SLN 67 PATH DIR 0 0 -$(dl1) DIV 3 GRUP 32 EXTR TYP SLN VON SLN 1112 PATH DIR 0 0 -$(dl1) DIV 3 EXTR TYP SLN VON SLN 1213 PATH DIR 0 0 -$(dl1) DIV 3 EXTR TYP SLN VON SLN 1314 PATH DIR 0 0 -$(dl1) DIV 3 EXTR TYP SLN VON SLN 1415 PATH DIR 0 0 -$(dl1) DIV 3 GRUP 31 SPIE TYP QUAD VON 31 DELT GRP A 1 B 0 C 0 DIST $(bm) GRUP 32 SPIE TYP QUAD VON 32 DELT GRP A 1 B 0 C 0 DIST $(bm) GRUP 30 FEDE NR QGRP KA 31 KE 32 CP 10E6 CQ 10E6 RISS 0 $ Federkopplung zwischen Duebel und Beton - UNTEN QUAD PROP MNR 0 NRA 0 T 0 GRUP 41 EXTR TYP SLN VON SLN 34 PATH DIR 0 0 $(dl2) DIV 3 EXTR TYP SLN VON SLN 45 PATH DIR 0 0 $(dl2) DIV 3 EXTR TYP SLN VON SLN 56 PATH DIR 0 0 $(dl2) DIV 3 EXTR TYP SLN VON SLN 67 PATH DIR 0 0 $(dl2) DIV 3 GRUP 42 EXTR TYP SLN VON SLN 1112 PATH DIR 0 0 $(dl2) DIV 3 EXTR TYP SLN VON SLN 1213 PATH DIR 0 0 $(dl2) DIV 3 EXTR TYP SLN VON SLN 1314 PATH DIR 0 0 $(dl2) DIV 3 EXTR TYP SLN VON SLN 1415 PATH DIR 0 0 $(dl2) DIV 3 GRUP 43 TRAN TYP QUAD VON 41 DELT GRP DZ ($(df)) GRUP 44 TRAN TYP QUAD VON 42 DELT GRP DZ ($(df)) GRUP 43 SPIE TYP QUAD VON 43 DELT GRP A 1 B 0 C 0 DIST $(bm) GRUP 44 SPIE TYP QUAD VON 44 DELT GRP A 1 B 0 C 0 DIST $(bm) GRUP 40 FEDE NR QGRP KA 43 KE 44 CP 10E6 CQ 10E6 $ Gruppierung der Federelemente am Rand - OBEN Let#xmin ($(bm)-#r2-0.001) Let#xmax ($(bm)+#r2+0.001) Let#ymin ($(lm)-#r2-0.001) Let#ymax ($(lm)+#r2+0.001) Let#zmin ($(to)-$(dl1)-0.001) Let#zmax ($(to)-$(dl1)+0.001) DEL TYP FEDE VON 30 DELT GRP XMIN #xmin XMAX #xmax YMIN #ymin YMAX #ymax ZMIN #zmin ZMAX #zmax GRUP 35 Let#xmin ($(bm)-#r2-0.001) Let#xmax ($(bm)+#r2+0.001) Let#ymin ($(lm)-#r2-0.001) Let#ymax ($(lm)+#r2+0.001) Let#zmin ($(to)-0.001) Let#zmax ($(to)+0.001) TRAN TYP FEDE VON 30 DELT GRP XMIN #xmin XMAX #xmax YMIN #ymin YMAX #ymax ZMIN #zmin ZMAX #zmax DZ (-$(dl1)) !*!Label Generierung der Auflager $ Generierung der Auflager $ in z-Richtung bei z = tu + df + to Let#xmin (-0.01) Let#xmax ($(b)+0.01) Let#ymin (-0.01) Let#ymax ($(l)+0.01) Let#zmin ($(to)+$(df)+$(tu)-0.01) Let#zmax ($(to)+$(df)+$(tu)+0.01) MOD TYP KNOT XMIN #xmin YMIN #ymin ZMIN #zmin XMAX #xmax YMAX #ymax ZMAX #zmax FIX PZ $ in z-Richtung bei z = tu + df Let#xmin (-0.001) Let#xmax ($(bm)-#r2+0.001) Let#ymin (-0.001) Let#ymax ($(l)+0.001) Let#zmin ($(to)-0.001) Let#zmax ($(to)+0.001) MOD TYP KNOT XMIN #xmin YMIN #ymin ZMIN #zmin XMAX #xmax YMAX #ymax ZMAX #zmax FIX PZ Let#xmin ($(bm)+#r2-0.001) Let#xmax ($(b)+0.001) Let#ymin (-0.01) Let#ymax ($(l)+0.001) Let#zmin ($(to)-0.001) Let#zmax ($(to)+0.001) MOD TYP KNOT XMIN #xmin YMIN #ymin ZMIN #zmin XMAX #xmax YMAX #ymax ZMAX #zmax FIX PZ Let#xmin ($(bm)-#r2-0.001) Let#xmax ($(bm)+#r2+0.001) Let#ymin (-0.001) Let#ymax ($(lm)-#r2+0.001) Let#zmin ($(to)-0.001) Let#zmax ($(to)+0.001) MOD TYP KNOT XMIN #xmin YMIN #ymin ZMIN #zmin XMAX #xmax YMAX #ymax ZMAX #zmax FIX PZ Let#xmin ($(bm)-#r2-0.001) Let#xmax ($(bm)+#r2+0.001) Let#ymin ($(lm)+#r2+0.001) Let#ymax ($(l)+0.001) Let#zmin ($(to)-0.001) Let#zmax ($(to)+0.001) MOD TYP KNOT XMIN #xmin YMIN #ymin ZMIN #zmin XMAX #xmax YMAX #ymax ZMAX #zmax FIX PZ $ in x-Richtung bei x = 0 Let#xmin (-0.01) Let#xmax (+0.01) Let#ymin (-0.01) Let#ymax ($(l)+0.01) Let#zmin ($(to)+$(df)-0.001) Let#zmax ($(to)+$(df)+$(tu)+0.001) MOD TYP KNOT XMIN #xmin YMIN #ymin ZMIN #zmin XMAX #xmax YMAX #ymax ZMAX #zmax FIX PX $ in x-Richtung bei x = b Let#xmin ($(b)-0.01) Let#xmax ($(b)+0.01) Let#ymin (-0.01) Let#ymax ($(l)+0.01) Let#zmin ($(to)+$(df)-0.001) Let#zmax ($(to)+$(df)+$(tu)+0.001) MOD TYP KNOT XMIN #xmin YMIN #ymin ZMIN #zmin XMAX #xmax YMAX #ymax ZMAX #zmax FIX PX $ in y-Richtung bei y = 0 Let#xmin (-0.01) Let#xmax ($(b)+0.01) Let#ymin (-0.01) Let#ymax (+0.01) Let#zmin ($(to)+$(df)-0.001) Let#zmax ($(to)+$(df)+$(tu)+0.001) MOD TYP KNOT XMIN #xmin YMIN #ymin ZMIN #zmin XMAX #xmax YMAX #ymax ZMAX #zmax FIX PY Let#xmin (-0.01) Let#xmax ($(b)+0.01) Let#ymin (-0.01) Let#ymax (+0.01) Let#zmin (-0.001) Let#zmax ($(to)+0.001) MOD TYP KNOT XMIN #xmin YMIN #ymin ZMIN #zmin XMAX #xmax YMAX #ymax ZMAX #zmax FIX PY $ in y-Richtung bei y = l Let#xmin (-0.01) Let#xmax ($(b)+0.01) Let#ymin ($(l)-0.01) Let#ymax ($(l)+0.01) Let#zmin ($(to)+$(df)-0.001) Let#zmax ($(to)+$(df)+$(tu)+0.001) MOD TYP KNOT XMIN #xmin YMIN #ymin ZMIN #zmin XMAX #xmax YMAX #ymax ZMAX #zmax FIX PY Let#xmin (-0.01) Let#xmax ($(b)+0.01) Let#ymin ($(l)-0.01) Let#ymax ($(l)+0.01) Let#zmin (-0.001) Let#zmax ($(to)+0.001) MOD TYP KNOT XMIN #xmin YMIN #ymin ZMIN #zmin XMAX #xmax YMAX #ymax ZMAX #zmax FIX PY $ Entfernen der Flaechenelemente & Kopplungen DEL QUAD VON 1 BIS 99 DELT GRP DEL KOPP VON 1 BIS 99 DELT GRP ENDE +PROG SOFILOAD urs:4 KOPF Test-Verformung ECHO Last Voll; ECHO Act Voll $ Verschiebung - OBEN sto#mm 10 $ Anzahl Lastschritte let#dL_max 0.5[mm] $ maximale Verschiebung in mm let#lf 1 $ Nummer Lastfall LOOP #mm LF #lf TYP none BEZ 'Verschiebung in X' AREA TYP WZ P1 ((#dL_max/#mm)*#lf) X1 $(b) Y1 0 Z1 0 X2 $(b) Y2 0 Z2 $(to) X3 $(b) Y3 $(l) Z3 $(to) X4 $(b) Y4 $(l) Z4 0 let#lf #lf+1 ENDLOOP ENDE +PROG ASE urs:5 KOPF Spannungen-Elastisch ECHO voll nein; ECHO grup ja; ECHO mat ja; ECHO fede ja STEU SOLV 4 let#start 0 $ Abfrage für Beginn der Primärlastfälle let#lf 1 $ Nummer erster Lastfall für Verschiebungseinwirkung $-------------------------------------- $ Gruppensteuerung fuer Verschiebung in Laengsrichtung $ erster Lastfall ohne Eigengewicht SYST PROB line ITER 1000 $ Rechnung linear $ Gruppensteuerung grup wert line lf 1000 egz 0 ende $ alle Folgelastfaelle fuer die Horizontalverschiebung loop #mm SYST PLF 1000+#lf-1 PROB line ITER 1000 $ Lasten auf zuvor berechneten Lastfall aufbringen $ und ab "prob" lineare Berechnung $ Gruppensteuerung grup wert line lf 1000+#lf BEZ 'u-max' lc #lf ENDE let#lf #lf+1 endloop ENDE