finish
/clear



!===========================Parameter============================
ri=10.2
alpha1=0.93419954
alpha2=5.59933934


!Funktionen
b31=0.01969627                            !1. Funktion
b21=0.1930295
b11=0.30908763
b01=9.33646544

b42=0.39372145                            !2. Funktion
b22=-0.68722407
b02=9.5

b33=-0.01969627                           !3. Funktion
b23=0.1930295
b13=-0.30908763
b03=9.33646544




!=========================Preprocessing==========================
/prep7                                    !Preprozessor starten

et,1,plane42                              !Elementtyp festlegen
type,1

local,11,1,0,0,0,0,0,0,0                  !Koordinatensystem festlegen
csys,11



!=======================Modellgenerierung========================
!--------------------------Keypoints-----------------------------
*get,kp,kp,,num,max

winkel=-9.5
*do,ii,1,39
   *if, winkel,ge,-alpha2,and,winkel,le,-alpha1,then
      rad=(b31*winkel**3)+(b21*winkel**2)+(b11*winkel)+b01
      k,kp+ii,rad,winkel
   *elseif,winkel,gt,-alpha1,and,winkel,lt,alpha1,then
      rad=(b42*winkel**4)+(b22*winkel**2)+b02
      k,kp+ii,rad,winkel
   *elseif,winkel,ge,alpha1,and,winkel,le,alpha2,then
      rad=(b33*winkel**3)+(b23*winkel**2)+(b13*winkel)+b03
      k,kp+ii,rad,winkel
   *elseif,winkel,lt,-alpha2,or,winkel,gt,alpha2,then
      k,kp+ii,ri,winkel
   *endif
   winkel=winkel+0.5
*enddo



!--------------------------Linien--------------------------------

larc,1,3,2,-ri
*do,ii,1,35
   ltan,ii,ii+3
*enddo

*get,lmax,line,,num,max
*do,ii,2,lmax
   lcomb,1,ii
*enddo