Finish
/clear
! STARTEN PREPROCESSEOR
/prep7

/vup,-1,z       ! Vertikalachse
/view,1,1,1,1   ! Ansichtspunkte
/triad,lbot     ! Koordinatensystemansicht      ,orig=Ursprung
/plopts,info,1

! EINHEITEN
!mm,N,s,t

! VARIABLEN
LAENGE = 1000 !mm
Breite = 800 !mm
HALTERUNG = 50 !mm überstand der EPDM Halterung
DEPDM = 4 !mm Dicke der EPDM schicht
Hoehe1 = 12 !mm !Scheibendicke 1 
HOEHE2 = 0 !mm PVB-Folie
Hoehe3 = 0 !mm !Scheibendicke 2
PMAX = 0.058 ! N/mm²  1kpa=0.001N/mm²
TPLUS  = 0.022 ! Dauer des positiven Drucks in ms
AKONSTANTE = 0.65 !Zerfallskoeffizient
NETZ1 = HOEHE1 !Vernetzung Glasscheiben
NETZ3 = DEPDM !Vernetzung EPDM

! ELEMENTTYP
et,1,185 !Glasscheiben
et,2,185 !PVB-Folie
et,3,185 !EPDM

! Material
!Glas
mp,ex,1,70000   ! E-Modul
mp,nuxy,1,0.23  ! Querdehnzahl
mp,dens,1,2.5e-9  ! Dichte

!EPDM
mp,ex,3,1.3    ! E-Modul
mp,nuxy,3,0.49  ! Querdehnzahl
mp,dens,3,1.13e-9        ! Dichte


!KEYPOINTS,x,y,z
!Ausnutzung der Symetrie

!Scheibe 1
k,1,(0),0,0
k,91,(0),0,0
k,2,(0),LAENGE/2,0
k,92,(0),LAENGE/2,0
k,3,(0),LAENGE/2,BREITE/2
k,4,(0),0,Breite/2
k,94,(0),0,Breite/2

!SCHEIBE 2
k,105,(HOEHE1+HOEHE2+HOEHE3),0,0
k,106,(HOEHE1+HOEHE2+HOEHE3),LAENGE/2,0
k,107,(HOEHE1+HOEHE2+HOEHE3),LAENGE/2,BREITE/2
k,108,(HOEHE1+HOEHE2+HOEHE3),0,Breite/2
k,9105,(HOEHE1+HOEHE2+HOEHE3),0,0
k,9106,(HOEHE1+HOEHE2+HOEHE3),LAENGE/2,0
k,9108,(HOEHE1+HOEHE2+HOEHE3),0,Breite/2

!# Lange Einspannseite
!EPDM
k,34,(-DEPDM),-HALTERUNG,0
k,35,(-DEPDM),-HALTERUNG,-HALTERUNG
k,36,(-DEPDM),LAENGE/2,-HALTERUNG
k,37,(-DEPDM),LAENGE/2,0
k,80,(-DEPDM),0,-HALTERUNG
k,38,(0),-HALTERUNG,0
k,938,(0),-HALTERUNG,0
k,39,(0),-HALTERUNG,-HALTERUNG
k,939,(0),-HALTERUNG,-HALTERUNG
k,40,(0),LAENGE/2,-HALTERUNG
k,940,(0),LAENGE/2,-HALTERUNG
k,81,(0),0,-HALTERUNG
k,981,(0),0,-HALTERUNG

!Glas2
k,21,(HOEHE1+HOEHE2+HOEHE3),-HALTERUNG,0
k,22,(HOEHE1+HOEHE2+HOEHE3),-HALTERUNG,-HALTERUNG
k,23,(HOEHE1+HOEHE2+HOEHE3),LAENGE/2,-HALTERUNG
k,921,(HOEHE1+HOEHE2+HOEHE3),-HALTERUNG,0
k,922,(HOEHE1+HOEHE2+HOEHE3),-HALTERUNG,-HALTERUNG
k,923,(HOEHE1+HOEHE2+HOEHE3),LAENGE/2,-HALTERUNG
k,84,(HOEHE1+HOEHE2+HOEHE3),0,-HALTERUNG
k,984,(HOEHE1+HOEHE2+HOEHE3),0,-HALTERUNG

!EPDM
k,25,(HOEHE1+HOEHE2+HOEHE3+DEPDM),-HALTERUNG,0
k,26,(HOEHE1+HOEHE2+HOEHE3+DEPDM),-HALTERUNG,-HALTERUNG
k,27,(HOEHE1+HOEHE2+HOEHE3+DEPDM),LAENGE/2,-HALTERUNG
k,28,(HOEHE1+HOEHE2+HOEHE3+DEPDM),LAENGE/2,0
k,85,(HOEHE1+HOEHE2+HOEHE3+DEPDM),0,-HALTERUNG

!## Kurze Einspannseite
!EPDM
k,50,(-DEPDM),0,0
k,52,(-DEPDM),-HALTERUNG,Breite/2
k,53,(-DEPDM),0,Breite/2
k,56,(0),-HALTERUNG,Breite/2
k,956,(0),-HALTERUNG,Breite/2

!Glas2
k,68,(HOEHE1+HOEHE2+HOEHE3),-HALTERUNG,Breite/2
k,968,(HOEHE1+HOEHE2+HOEHE3),-HALTERUNG,Breite/2

!EPDM
k,70, (HOEHE1+HOEHE2+HOEHE3+DEPDM),0,0
k,72,(HOEHE1+HOEHE2+HOEHE3+DEPDM),-HALTERUNG,Breite/2
k,73,(HOEHE1+HOEHE2+HOEHE3+DEPDM),0,Breite/2


!VOLUMEN
!Fläche
v,1,2,3,4,105,106,107,108!V1 Glas Fläche 
!Lange Seite
v,80,36,37,50,981,940,92,91  !V4 EPDM
v,81,40,2,1, 84,23,106,105 !V5 Glas1
v,984,923,9106,9105,85,27,28,70  !V8 EPDM
!Kurze Seite
v,50,34,52,53,91,938,956,94  !V9 EPDM
v,1,38,56,4,105,21,68,108    !V10 Glas 
v,9105,921,968,9108,70,25,72,73 !V13 EPDM
!Quadrat im Eck
v,80,35,34,50,981,939,938,91 !v14 EPDM
v,81,39,38,1,84,22,21,105 !V15 Glas 
v,984,922,921,9105,85,26,25,70 !V18 EPDM

! Elementgröße und Vernetzen
esize,NETZ1
vsel,s,volu,,1
vsel,a,volu,,3
vsel,a,volu,,6
vsel,a,volu,,9
vatt,1,,1
vmesh,all
allsel,all

esize,NETZ3
vsel,s,volu,,2
vsel,a,volu,,4
vsel,a,volu,,5
vsel,a,volu,,7
vsel,a,volu,,8
vsel,a,volu,,10
vatt,3,,3
vmesh,all


!Lagerung für symetrischen Fall
nsel,s,loc,x,-DEPDM
nsel,a,loc,x,(HOEHE1+HOEHE2+HOEHE3+DEPDM)
d,all,all       

!Symmetrieflächen glas

DA,4 ,SYMM
DA,27 ,SYMM
DA,33 ,SYMM
DA,38,symm 
DA,3 ,SYMM
DA,9 ,SYMM
DA,15 ,SYMM
DA,20 ,SYMM 

!Festhaltung der Glasscheibe in X Richtung 
dk,39,ux

!###################### KONTAKTMANAGER LOG ###################
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
MP,MU,1,
MAT,1   
MP,EMIS,1,7.88860905221e-031
R,3 
REAL,3  
ET,4,170
ET,5,174
KEYOPT,5,9,0
KEYOPT,5,10,2   
R,3,
RMORE,  
RMORE,,0
RMORE,0 
ASEL,S,,,13 
CM,_TARGET,AREA 
TYPE,4  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
NSLE,A,CT2  
ESURF   
CMSEL,S,_ELEMCM 
ASEL,S,,,12 
CM,_CONTACT,AREA
TYPE,5  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,3  
ESEL,ALL
ESEL,S,TYPE,,4  
ESEL,A,TYPE,,5  
ESEL,R,REAL,,3  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
CMDEL,_TARGET   
CMDEL,_CONTACT  
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
MP,MU,1,0   
MAT,1   
R,4 
REAL,4  
ET,6,170
ET,7,174
KEYOPT,7,9,0
KEYOPT,7,10,2   
R,4,
RMORE,  
RMORE,,0
RMORE,0 
ASEL,S,,,45 
CM,_TARGET,AREA 
TYPE,6  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
NSLE,A,CT2  
ESURF   
CMSEL,S,_ELEMCM 
ASEL,S,,,44 
CM,_CONTACT,AREA
TYPE,7  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,6  
ESEL,A,TYPE,,7  
ESEL,R,REAL,,4  
ESEL,ALL
ESEL,S,TYPE,,6  
ESEL,A,TYPE,,7  
ESEL,R,REAL,,4  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
CMDEL,_TARGET   
CMDEL,_CONTACT  
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
MP,MU,1,0   
MAT,1   
R,5 
REAL,5  
ET,8,170
ET,9,174
KEYOPT,9,9,0
KEYOPT,9,10,2   
R,5,
RMORE,  
RMORE,,0
RMORE,0 
ASEL,S,,,48 
CM,_TARGET,AREA 
TYPE,8  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
NSLE,A,CT2  
ESURF   
CMSEL,S,_ELEMCM 
ASEL,S,,,49 
CM,_CONTACT,AREA
TYPE,9  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,8  
ESEL,A,TYPE,,9  
ESEL,R,REAL,,5  
ESEL,ALL
ESEL,S,TYPE,,8  
ESEL,A,TYPE,,9  
ESEL,R,REAL,,5  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
CMDEL,_TARGET   
CMDEL,_CONTACT  
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
MP,MU,1,0   
MAT,1   
R,6 
REAL,6  
ET,10,170   
ET,11,174   
KEYOPT,11,9,0   
KEYOPT,11,10,2  
R,6,
RMORE,  
RMORE,,0
RMORE,0 
ASEL,S,,,17 
CM,_TARGET,AREA 
TYPE,10 
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
NSLE,A,CT2  
ESURF   
CMSEL,S,_ELEMCM 
ASEL,S,,,18 
CM,_CONTACT,AREA
TYPE,11 
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,10 
ESEL,A,TYPE,,11 
ESEL,R,REAL,,6  
ESEL,ALL
ESEL,S,TYPE,,10 
ESEL,A,TYPE,,11 
ESEL,R,REAL,,6  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
CMDEL,_TARGET   
CMDEL,_CONTACT  
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
MP,MU,1,0   
MAT,1   
R,7 
REAL,7  
ET,12,170   
ET,13,174   
KEYOPT,13,9,0   
KEYOPT,13,10,2  
R,7,
RMORE,  
RMORE,,0
RMORE,0 
ASEL,S,,,30 
CM,_TARGET,AREA 
TYPE,12 
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
NSLE,A,CT2  
ESURF   
CMSEL,S,_ELEMCM 
ASEL,S,,,29 
CM,_CONTACT,AREA
TYPE,13 
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,12 
ESEL,A,TYPE,,13 
ESEL,R,REAL,,7  
ESEL,ALL
ESEL,S,TYPE,,12 
ESEL,A,TYPE,,13 
ESEL,R,REAL,,7  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
CMDEL,_TARGET   
CMDEL,_CONTACT  
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
MP,MU,1,0   
MAT,1   
R,8 
REAL,8  
ET,14,170   
ET,15,174   
KEYOPT,15,9,0   
KEYOPT,15,10,2  
R,8,
RMORE,  
RMORE,,0
RMORE,0 
ASEL,S,,,34 
CM,_TARGET,AREA 
TYPE,14 
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
NSLE,A,CT2  
ESURF   
CMSEL,S,_ELEMCM 
ASEL,S,,,35 
CM,_CONTACT,AREA
TYPE,15 
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,14 
ESEL,A,TYPE,,15 
ESEL,R,REAL,,8  
ESEL,ALL
ESEL,S,TYPE,,14 
ESEL,A,TYPE,,15 
ESEL,R,REAL,,8  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
CMDEL,_TARGET   
CMDEL,_CONTACT  

finish

!################################# SOLUTION ####################################
/solu
antype,trans
!betad,0.0001 !Dämpfung

/outp,ber,out

autots,off
lnsrch,off
nlgeom,on
sstiff,on
neqit,50
time,0
cnvtol,u,,1.0e-4,,1.0
outres,all,all
eqslv,sparse

!#Anfangsstellung
timint,on
nsubst,1,1,1	! mehrere Steps eingeben nsubst,10,10,10
time,0.000001
kbc,1
sfa,1,,pres, PMAX 
allsel
Lswrite,1

!# Explosionslast
timint,on
!Schleife
*do,ZEIT,0.001,TPLUS,0.001
TIME,ZEIT			
KBC,0
sfa,1, ,pres,PMAX*(1-ZEIT/TPLUS)*exp(-AKONSTANTE*(ZEIT/TPLUS))
nsubst,5,5,5
allsel
Lswrite,(ZEIT*1000+1)
*enddo

!# Ausschwingvorgang
sfadele,all,all,all
nsubst,10,10,10
time,0.05
allsel
Lswrite,(ZEIT*1000+2)

LSSOLVE,1,(ZEIT*1000+2),1

finish