***************************************************************************** ***************************************************************************** *preprint,model=no ***************************************************************************** ***************************************************************************** ** **************************************************************************** *********************************** ** PARAMETER *********************************** *parameter t=10.25 ts=4.54 tpvb=1.52 h=251 l=763.04 *********************************** Eg=70000 vg=0.23 Ep=3.892 vp=0.39 Es=2100000 vs=0.3 *********************************** t1=ts+tpvb k=10 ku=-k t2=ts+ts+t+tpvb+tpvb ku2=-0.5*k ***************************************************************************** ***************************************************************************** ** ** KNOTEN *********************************** ** Scheibe *********************************** *node 1,0,0,0 46001,0,,0 *ngen,nset=Al 1,46001,1000 *ncopy,change number=12,old set=Al,new set=Ar,shift ,0,0 0,0,0,0,0,0,0 *nfill,nset=A Al,Ar,12,1 *ncopy,change number=100,old set=A,new set=B,shift 0,0, 0,0,0,0,0,0,0 *nfill,nset=glas1 A,B,2,50 *ncopy,change number=100,old set=B,new set=P1,shift 0,0, 0,0,0,0,0,0,0 *nfill,nset=pvb1 B,P1,2,50 *ncopy,change number=100,old set=P1,new set=S1,shift 0,0, 0,0,0,0,0,0,0 *nfill,nset=glas2 P1,S1,2,50 *ncopy,change number=100,old set=S1,new set=P2,shift 0,0, 0,0,0,0,0,0,0 *nfill,nset=pvb2 S1,P2,2,50 *ncopy,change number=100,old set=P2,new set=S2,shift 0,0, 0,0,0,0,0,0,0 *nfill,nset=pvb2 P2,S2,2,50 *** Lasteinleitung oben ********************************* *nset,nset=le1 46001 46051 46101 46151 46201 46251 46301 46351 46401 46451 46501 *ncopy,change number=12,old set=le1,new set=le2,shift ,0,0 0,0,0,0,0,0,0 *nfill,nset=le3 le1,le2,12,1 *ncopy,change number=2000,old set=le3,new set=le4,shift 0,,0 0,0,0,0,0,0,0 *nfill,nset=le5 le3,le4,2,1000 *********************************** ***Lasteinleitung unten ********************* *nset,nset=le6 1 51 101 151 201 251 301 351 401 451 501 *ncopy,change number=12,old set=le6,new set=le7,shift ,0,0 0,0,0,0,0,0,0 *nfill,nset=le8 le6,le7,12,1 *ncopy,change number=200000,old set=le8,new set=le9,shift 0,,0 0,0,0,0,0,0,0 *nfill,nset=le10 le8,le9,2,100000 *********************************** ** Sets *********************************** *imperfection,file=Eigenform_vsg1_Lasteinleitung,step=1 1,0.36 *************************************** ** ** beidseitig gekenkiges Lager *********************************** *nset,nset=Ga_X 200251 48251 *nset,nset=Ga_Y 200251 200252 200253 200254 200255 200256 200257 200258 200259 200260 200261 200262 200263 *nset,nset=Ga_Z 200251 200257 200263 48251 48263 48257 *********************************** ** Output *********************************** *nset,nset=Last 48251 48252 48253 48254 48255 48256 48257 48258 48259 48260 48261 48262 48263 *nset,nset=Last_V 23501 23502 23503 23504 23505 23506 23507 23508 23509 23510 23511 23512 23513 *nset,nset=mitte 23007 23257 23507 ***************************************************************************** ***************************************************************************** ** ** ELEMENTE ************************************ ** Scheibe ************************************ *element,type=C3D8I 1,1,51,52,2,1001,1051,1052,1002 *elgen,elset=glas1 1,12,1,1,2,50,50,46,1000,1000 *element,type=C3D8IH 101,101,151,152,102,1101,1151,1152,1102 *elgen,elset=pvb1 101,12,1,1,2,50,50,46,1000,1000 *element,type=C3D8I 201,201,251,252,202,1201,1251,1252,1202 *elgen,elset=glas2 201,12,1,1,2,50,50,46,1000,1000 *element,type=C3D8IH 301,301,351,352,302,1301,1351,1352,1302 *elgen,elset=pvb2 301,12,1,1,2,50,50,46,1000,1000 *element,type=C3D8I 401,401,451,452,402,1401,1451,1452,1402 *elgen,elset=glas3 401,12,1,1,2,50,50,46,1000,1000 *** Lasteinleitung oben *************** *element,type=C3D6,elset=last1 47001,46052,47052,47002,46051,47051,47001 47201,46252,47252,47202,46251,47251,47201 47401,46452,47452,47402,46451,47451,47401 *elgen,elset=last2 47001,12,1,1 47201,12,1,1 47401,12,1,1 *element,type=C3D6,elset=last3 47051,47102,47052,46052,47101,47051,46051 47251,47302,47252,46252,47301,47251,46251 47451,47502,47452,46452,47501,47451,46451 *elgen,elset=last4 47051,12,1,1 47251,12,1,1 47451,12,1,1 *element,type=C3D8I 48001,47001,47051,47052,47002,48001,48051,48052,48002 *elgen,elset=last5 48001,12,1,1,10,50,50 *elset,elset=lastoben last1 last2 last3 last4 last5 ***************************** *** Lasteinleitung unten ***************************** *************** *element,type=C3D6,elset=last6 100001,100002,100052,52,100001,100051,51 100201,100202,100252,252,100201,100251,251 100401,100402,100452,452,100401,100451,451 *elgen,elset=last7 100001,12,1,1 100201,12,1,1 100401,12,1,1 *element,type=C3D6,elset=last8 100051,52,100052,100102,51,100051,100101 100251,252,100252,100302,251,100251,100301 100451,452,100452,100502,451,100451,100501 *elgen,elset=last9 100051,12,1,1 100251,12,1,1 100451,12,1,1 *element,type=C3D8I 200001,200001,200051,200052,200002,100001,100051,100052,100002, *elgen,elset=last10 200001,12,1,1,10,50,50 *elset,elset=lastunten last6 last7 last8 last9 last10 ************************************ ** sets ************************************ *elset,elset=Spannung 23456 23007 *elset,elset=glas glas1 glas2 glas3 *elset,elset=pvb pvb1 pvb2 *elset,elset=metall lastoben lastunten ***************************************************************************** ***************************************************************************** ** ** MATERIAL ************************************ ** Glas ************************************ *solid section,elset=glas,material=glas *material,name=glas *elastic , *solid section,elset=pvb,material=pvb *material,name=pvb *elastic , *solid section,elset=metall,material=metall *material,name=metall *elastic , ***************************************************************************** ***************************************************************************** ** ** LAGERUNG ************************************ *boundary Ga_X,1 Ga_Y,2 Ga_Z,3 ************************************ ** BELASTUNG Vertikal *************************************************************************** *step,amplitude=RAMP,nlgeom *static 0.01,20, 1e-30,10 *CLOAD Last,2,-3107.69 *NODE FILE,FREQ=0 *EL FILE,FREQ=0 *restart,write, FREQ=1 *monitor,dof=3,node=45213 *output,field,frequency=1,variable=all *output, history,frequency=1 *el print, freq=1,elset=Spannung,summary=no,totals=no,position=nodes s22 *node print,freq=1,nset=mitte,summary=no,totals=no u3 *node print,freq=1,nset=Last,summary=no,totals=no CF2 *output,history,frequency=1 *node output,nset=Last CF2 u2 *node output,nset=mitte u3 *end step ************************************************************************ ** Belastung Horizontal ************************************************************************ *step, amplitude=RAMP, nlgeom *static 0.01,20, 1e-30,10 *CLOAD Last_V, 3,3000 Ga_X,1 Ga_Y,2 Ga_Z,3 *NODE FILE,FREQ=0 *EL FILE,FREQ=0 *restart,write, FREQ=1 *monitor,dof=3,node=45213 *output,field,frequency=1,variable=all,OP=Add *output,history,frequency=1,OP=Add *el print, freq=1,elset=Spannung,summary=no,totals=no,position=nodes s22 *node print,freq=1,nset=mitte,summary=no,totals=no u3 *node print,freq=1,nset=Last_V,summary=no,totals=no CF3 *output,history,frequency=1,OP=Add *node output,nset=Last_V CF3 *node output,nset=mitte u3 *end step