finish
/clear, start

!**********Tabellenwerte*********************

	Ex_cfk=32.71e9 !in Pa
	Ey_cfk=37.71e9
	Ez_cfk=4e9
	pr_cfk=0.1
	GXY_cfk=3.26e9
	dens_cfk=1487.8
	Winkel_cfk=0

!***********Semperit**********************************************

	Ex_semperit=0.05e9 !Startwert 0.05e9
	pr_semperit=0.5
	dens_semperit=873.84	!aus Experiment, für jede Konfig anders weil große Schwankung, Startwert 1221.62
	


div_=5

!*DO,Winkel_,0,45,15


parsav,,'parameters_systematik',par !saves parameters to file
finish
/clear, start !clears all data including parameters
parres,,'parameters_systematik',par !reads parameters from file	
	
/prep7


!*******Material 1 - CFK********

mp, ex,1,Ex_cfk !in SI-Units, kg,m,Pa, erster wert 1.392e11
mp, ey,1,Ey_cfk
mp, ez,1,Ez_cfk	
mp, prxy,1,pr_cfk !alter Wert 0.26
mp, pryz,1,pr_cfk
mp, prxz,1,pr_cfk
mp, gxy,1,GXY_cfk !alter Wert 5.6e9
mp, gyz,1,GXY_cfk
mp, gxz,1,GXY_cfk
mp, dens,1,dens_cfk

!*******Material 2 - Gummi*******

mp,ex,2,Ex_semperit
mp,dens,2,dens_semperit
mp,prxy,pr_semperit

!********Geometrie*******
L_1=0.180 
L_2=0.02
L_3=0.00184
Dicke_1=0.00043
Dicke_2=0.000179
block,0,L_1,0,L_2,0,L_3 !block,x1,x2,y1,y2,z1,z2


et,1,190
keyopt,1,8,1
keyopt,1,2,1

sectype,1,shell,,LAM_1
secdata, Dicke_1,1,Winkel_,5 !secdata, dicke, mat, winkel, integrationpoints, layername
secdata, Dicke_1,1,Winkel_,5

sectype,2,shell,,Elast_1
secdata, Dicke_2,2,0,5

sectype,3,shell,,LAM_2
secdata, Dicke_1,1,Winkel_,5
secdata, Dicke_1,1,Winkel_,5


!*******Vernetzung********
elemsi=L_3 

nummrg,kp !merges items, kp->keypoints
esize,elemsi !number of line divisions
lesize,9,,,div_
lesize,10,,,div_
lesize,11,,,div_
lesize,12,,,div_

veorient,1,thin !elementorientierung
/eshape,0,0

vmesh,1,1,0.1

!*******Randeinspannung bei x=0*********
csys,0
nsel,s,loc,x,,0
d,all,all

!********Lösung**********

/solu
alls
Antype,MODAL
modenanz=1
modopt, lanb, modenanz,0.1,,,off !mode-extraction
mxpa,modenanz,,,yes
solve

/Post1
SET,First
set,list

!*Enddo