/PREP7
*AFUN,DEG

!SETTINGS
!--------

*SET,AUFLAGERHOEHE,1.5
*SET,AUFLAGERLINIE,2
*SET,RANDLINIE,29
*SET,ZENTRALLINIE,20
*SET,DIAGONALLINIE,39
*SET,X,0				!VERTIKALER ABSTAND ZWISCHEN FUß UND URSPRUNG DES SYSTEMS 
*SET,D,0.05				!NORMALES BETON

*SET,HSCHALE,4.671779316 		!Zu abmessen A
*SET,AREA,324.952811			!Zu abmessen R
*SET,PROJAREA,324.952811		!Zu abmessen R
*SET,RMAXSCHALE,14.32881511		!Zu abmessen A
*SET,ALPHA,24.1409			!Zu abmessen A
*SET,DXMAX,10.13200233			!Zu abmessen A
*SET,FUSSNEIGUNG,24.1409   		!Zu abmessen R+A

*SET,DICHTE,2550			!Kg/m³
*SET,EX,30e+09				!N/m²
*SET,PRXY,0.15				!()

*SET,BELASTUNG,2			!1=EIGG,2=SCHNEELAST,3=WINDLAST,4=ADAPTION

*SET,SCHNEELASTFALL,1			!1=SYM,2=ANTISYM_HAELFTE,3=ANTISYM_45_UMRICHTET,4=ANTISYM_VIERTEL
*SET,SCHNEEDRUCK_K,1858.4		!N/m² VAIHINGEN - 500M
*SET,SCHNEEDRUCK_I_GL,0.8*(1858.4)	!N/m² VAIHINGEN - 500M
*SET,SCHNEEDRUCKO,1858.4		!N/m²
*SET,SCHNEEDRUCKU,1383.5		!N/m²
*SET,WINDLASTFALL,2			!1=X_RICHT,2=X_RICHT_45_GEDREHT
*SET,WINDDRUCK,840			!N/m²
*SET,STOERUNG,0				!1=STOERUNG,0=KEINE STOERUNG
*SET,XCYL,4
*SET,YCYL,4
*SET,RCYL,1.5
				
*SET,AVX1_1,0.001				!>0,NACHAUSSN,<0,NACHINN
*SET,AVX1_2,0.001				!>0,NACHAUSSN,<0,NACHINN
*SET,AVZ1_1,0.001				!>0,OBEN
*SET,AVZ1_2,0.001				!>0,OBEN

!GEOMETRIE
!---------

LOCAL,11,1,0,0,X+AUFLAGERHOEHE
LOCAL,12,0,0,0,0,45
CSYS,0

K,101,50,50,AUFLAGERHOEHE+X
K,102,50,-50,AUFLAGERHOEHE+X
K,103,-50,50,AUFLAGERHOEHE+X
K,104,-50,-50,AUFLAGERHOEHE+X
A,101,102,104,103
ASBA,1,2
APLOT
ADELE,3,6,1,1

K,105,,50,50
K,106,,50,-50
K,107,,-50,50
K,108,,-50,-50
A,105,106,108,107
ASBA,ALL,1

K,105,50,,50
K,106,50,,-50
K,107,-50,,50
K,108,-50,,-50
A,105,106,108,107
ASBA,ALL,1

K,105,50,50,50
K,106,50,50,-50
K,107,-50,-50,50
K,108,-50,-50,-50
A,105,106,108,107
ASBA,ALL,1

K,105,50,-50,50
K,106,50,-50,-50
K,107,-50,50,50
K,108,-50,50,-50
A,105,106,108,107
ASBA,ALL,1

!STOERUNG
!--------

*IF,STOERUNG,EQ,1,THEN
	CYL4,XCYL,YCYL,RCYL,,,,99
	ASBV,ALL,ALL
*ENDIF

!ELEMENTS
!--------

ET,1,SHELL93
KEYOPT,1,8,2
ET,2,SURF154          

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,DICHTE
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,EX
MPDATA,PRXY,1,,PRXY

R,1,D
R,2

MAT,1
REAL,1
TYPE,1
ESYS,0

LSEL,S,LOC,Z,AUFLAGERHOEHE+X
LESIZE,ALL,,,AUFLAGERLINIE
LSEL,S,LOC,Z,AUFLAGERHOEHE+X+0.1,99
LSEL,U,LOC,X,0,0.001
LSEL,U,LOC,Y,0,0.001
CSYS,12
LSEL,U,LOC,X,0,0.001
LSEL,U,LOC,Y,0,0.001
LESIZE,ALL,,,RANDLINIE
LSEL,S,LOC,X,0,0.001
LSEL,A,LOC,Y,0,0.001
LESIZE,ALL,,,DIAGONALLINIE
CSYS,0
LSEL,S,LOC,X,0,0.001
LSEL,A,LOC,Y,0,0.001
LESIZE,ALL,,,ZENTRALLINIE
ALLS
LPLOT

MSHAPE,1,2D
MSHKEY,0
AMESH,ALL

MAT,2
REAL,2
TYPE,2
NSEL,ALL
ESURF
ALLS

NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
CSYS,11
NROT,ALL
NMODIF,ALL,,,,,,FUSSNEIGUNG
/PSYMB,NDIR,1
NPLOT


!RANDBEDINGUNGEN
!---------------

/SOLU
CSYS,0
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
D,ALL,UX
D,ALL,UZ
ALLS
CSYS,0

NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
CSYS,2
NSEL,U,LOC,Y,45.02,67.5
NSEL,U,LOC,Y,27.5,44.98
NSEL,U,LOC,Y,100,134.98
NSEL,U,LOC,Y,135.02,140
NSEL,U,LOC,Y,200,224.98
NSEL,U,LOC,Y,225.02,250
NSEL,U,LOC,Y,285,314.98
NSEL,U,LOC,Y,315.02,360
D,ALL,UY
ALLS
CSYS,0

!BELASTUNG
!---------

ACEL,,,9.81						!1_EIGENGEWICHT	
LSWRITE,1	
LSCLEAR,INER								

/PSF,PRES,2						!2_SCHNEE VOLLSTÄNDIG
ESEL,S,TYPE,,2
SFE,ALL,5,PRES,,SCHNEEDRUCK_I_GL,0,0,-1
ALLS
LSWRITE,2
SFEDELE,ALL,ALL,ALL
	
									
/PSF,PRES,2						!3_SCHNEE_HAELFTE	
ESEL,S,TYPE,,2
ESEL,U,,,4773,5372,1
ESEL,U,,,7139,7742,1
ESEL,U,,,8345,8940,1
ESEL,U,,,5969,6558,1
SFE,ALL,5,PRES,,SCHNEEDRUCK_K,0,0,-1
ALLS
LSWRITE,3
SFEDELE,ALL,ALL,ALL


	!SCHNEELAST_GES = PROJAREA/2*SCHNEEDRUCK_I_GL		!2_SCHNEE_HAELFTE_45GRAD_GEDREHT
	!CSYS,12		
	!ASEL,S,LOC,X,0,99						
	!NSLA,,1	
	!*GET,NODENUMBER,NODE,,COUNT				
	!NSEL,U,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001			
	!F,ALL,FZ,SCHNEELAST_GES/NODENUMBER
	

	!SCHNEELAST_GES = PROJAREA/4*SCHNEEDRUCK_I_GL		!2_SCHNEE_VIERTEL			
	!ASEL,S,LOC,X,0,99
	!ASEL,R,LOC,Y,0,99						
	!NSLA,,1	
	!*GET,NODENUMBER,NODE,,COUNT				
	!NSEL,U,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001			
	!F,ALL,FZ,SCHNEELAST_GES/NODENUMBER
	
	!SFGRAD,PRES,0,X,0,(WINDDRUCK*COS(ALPHA)/DXMAX)		!3_WIND
	!CSYS,0	
	!SF,ALL,PRES,0
	!/PSF,PRES,2

		
/PSF,PRES,2						!4_WIND_45GRAD_GEDREHT
SFGRAD,PRES,12,X,0,(WINDDRUCK*COS(ALPHA)/RMAXSCHALE)		
ESEL,S,TYPE,,2
!ESEL,U,,,4773,5372,1	Retirer les noeuds de la moitie sur 45 degre tournee
!ESEL,U,,,7139,7742,1
!ESEL,U,,,8345,8940,1
!ESEL,U,,,6559,7138,1
SFE,ALL,1,PRES,,,,,
ALLS
LSWRITE,4	
SFEDELE,ALL,ALL,ALL
DDELE,ALL,ALL


!ADAPTION LOAD CASE
!------------------

!AVX1_1
!------


CSYS,0
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
CSYS,12
NSEL,U,LOC,X,2,99
D,ALL,UX
ALLS
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
D,ALL,UZ
CSYS,0
ALLS
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
CSYS,2
NSEL,U,LOC,Y,45.02,67.5
NSEL,U,LOC,Y,27.5,44.98
NSEL,U,LOC,Y,100,134.98
NSEL,U,LOC,Y,135.02,140
NSEL,U,LOC,Y,200,224.98
NSEL,U,LOC,Y,225.02,250
NSEL,U,LOC,Y,285,314.98
NSEL,U,LOC,Y,315.02,360
D,ALL,UY
CSYS,0
ALLS
NSEL,S,LOC,Z,X+AUFLAGERHOEHE,X+AUFLAGERHOEHE+0.001
NSEL,R,LOC,X,0,99
NSEL,R,LOC,Y,0,99
D,ALL,UX,AVX1_1
ALLS
LSWRITE,5
DDELE,ALL,ALL


!AVX1_2
!------

CSYS,0
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
CSYS,12
NSEL,U,LOC,Y,-2,-99
D,ALL,UX
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
D,ALL,UZ
CSYS,0
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
CSYS,2
NSEL,U,LOC,Y,45.02,67.5
NSEL,U,LOC,Y,27.5,44.98
NSEL,U,LOC,Y,100,134.98
NSEL,U,LOC,Y,135.02,140
NSEL,U,LOC,Y,200,224.98
NSEL,U,LOC,Y,225.02,250
NSEL,U,LOC,Y,285,314.98
NSEL,U,LOC,Y,315.02,360
D,ALL,UY
CSYS,0
NSEL,S,LOC,Z,X+AUFLAGERHOEHE,X+AUFLAGERHOEHE+0.001
NSEL,R,LOC,X,0,99
NSEL,R,LOC,Y,0,-99
D,ALL,UX,AVX1_2
ALLS
LSWRITE,6
DDELE,ALL,ALL


!AVZ1_1
!------

CSYS,0
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
CSYS,12
NSEL,U,LOC,X,2,99
D,ALL,UZ
ALLS
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
D,ALL,UX
CSYS,0
ALLS
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
CSYS,2
NSEL,U,LOC,Y,45.02,67.5
NSEL,U,LOC,Y,27.5,44.98
NSEL,U,LOC,Y,100,134.98
NSEL,U,LOC,Y,135.02,140
NSEL,U,LOC,Y,200,224.98
NSEL,U,LOC,Y,225.02,250
NSEL,U,LOC,Y,285,314.98
NSEL,U,LOC,Y,315.02,360
D,ALL,UY
CSYS,0
ALLS
NSEL,S,LOC,Z,X+AUFLAGERHOEHE,X+AUFLAGERHOEHE+0.001
NSEL,R,LOC,X,0,99
NSEL,R,LOC,Y,0,99
D,ALL,UZ,AVZ1_1
ALLS
LSWRITE,7
DDELE,ALL,ALL


!AVZ1_2
!------

CSYS,0
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
CSYS,12
NSEL,U,LOC,X,-2,-99
D,ALL,UZ
ALLS
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
D,ALL,UX
CSYS,0
ALLS
NSEL,S,LOC,Z,AUFLAGERHOEHE+X,AUFLAGERHOEHE+X+0.001
CSYS,2
NSEL,U,LOC,Y,45.02,67.5
NSEL,U,LOC,Y,27.5,44.98
NSEL,U,LOC,Y,100,134.98
NSEL,U,LOC,Y,135.02,140
NSEL,U,LOC,Y,200,224.98
NSEL,U,LOC,Y,225.02,250
NSEL,U,LOC,Y,285,314.98
NSEL,U,LOC,Y,315.02,360
D,ALL,UY
CSYS,0
ALLS
NSEL,S,LOC,Z,X+AUFLAGERHOEHE,X+AUFLAGERHOEHE+0.001
NSEL,R,LOC,X,0,99
NSEL,R,LOC,Y,0,-99
D,ALL,UZ,AVZ1_2
ALLS
LSWRITE,8
DDELE,ALL,ALL




CSYS,0
ALLS
NPLOT
SAVE
FINISH


!BERECHNUNG
!----------

/SOLU
LSSOLVE,1,8,1 

!POSTPROCESSING
!--------------

/POST1
/GRA,FULL
/EDGE,1
/PLOPTS,TITLE,0
/PLOPTS,DATE,0
SHELL,MID


!LOAD CASES
!----------

SET,1			!LOAD CASE 10 - EIGG					
LCDEF,10,1			
LCWRITE,10
SAVE
SET,2			!LOAD CASE 21 - SCHNEE VOLL						
LCDEF,21,2				
LCWRITE,21
SAVE
SET,3			!LOAD CASE 22 - SCHNEE HAELFTE					
LCDEF,22,3				
LCWRITE,22
SAVE
SET,4			!LOAD CASE 23 - WIND 45GRAD GEDREHT
LCDEF,23,4				
LCWRITE,23
SAVE
SET,5			!LOAD CASE 24 - AVX1_1
LCDEF,24,5				
LCWRITE,24
SAVE
SET,6			!LOAD CASE 25 - AVX1_2
LCDEF,25,6				
LCWRITE,25
SAVE
SET,7			!LOAD CASE 26 - AVZ1_1
LCDEF,26,7				
LCWRITE,26
SAVE
SET,8			!LOAD CASE 27 - AVZ1_2
LCDEF,27,8				
LCWRITE,27
SAVE


LCASE,21		!KOMBI 31 ( 10 + 21 )
LCOPER,ADD,10,,,
LCWRITE,31
LCASE,22		!KOMBI 32 ( 10 + 22 )
LCOPER,ADD,10,,,
LCWRITE,32
LCASE,23		!KOMBI 33 ( 10 + 23 )
LCOPER,ADD,10,,,
LCWRITE,33
LCASE,24		!KOMBI 59 ( 24 + 25 ) - AVX1_1 + AVX1_2
LCOPER,ADD,25,,,
LCWRITE,59


!OUTPUT
!------

!/SHOW,PNG
!/RGB,INDEX,100,100,100,0
!/RGB,INDEX,0,0,0,15

!ETABLE,MX,SMISC,4   
!ETABLE,MY,SMISC,5   
!ETABLE,MXY,SMISC,6
!SADD,MX+MY,MX,MY,1,1,0			

!BIEGESPANNUNGEN
!---------------

!LCASE,10
!GRA,FULL
!*DIM,SZ_T,,9793
!*DIM,SZ_M,,9793
!*DIM,SZ_B,,9793
!SHELL,MID
!*VGET,SZ_M(1),NODE,1,S,Z
!SHELL,TOP
!*VGET,SZ_T(1),NODE,1,S,Z
!*VOPER,SZ_B(1),SZ_T(1),SUB,SZ_M(1)
!*VPUT,SZ_B(1),NODE,1,S,Z


!TOPOLOGY OPTIMIZATION
!---------------------

!/OPT
!TOCOMP,COMP,SING
!TOVAR,COMP,OBJ
!TOVAR,VOLUME,CON,,50
!TOTYPE,OC
!TODEF,0.0001
!TOLOOP,10,1
!TOPLOT
!PLNS,TOPO