!˳ΪµĽеس
/Prep7
*AFUN,DEG    !ýǶȵλΪ
lx=3  !Ĭ񳤱߱߳
ly=3  !Ĭ̱߱߳
n=8  !Ĭϱ񳤱߷ڵ
m=6 !Ĭϱ̱߷ڵ 
z=-2.5 !±֮Ĵֱ
Et,1,link8         !嵥Ԫ
D1=0.03             !Ҹֱ
AREA1=3.14*d1*d1/4   !Ҹ
R,1,area1          !˵Ԫ
d2=0.04            !Ҹֱ
AREA2=3.14*d2*d2/4  !Ҹ
R,2,area2         !˵Ԫ
D3=0.02           !бֱ
AREA3=3.14*d3*d3/4  !б
R,3,area3        !˵Ԫ
MP,EX,1,2.1e11                !1ϵĵģEX
MP,PRXY,1,0.3          !1ϵĲɱPRXY
MP,DENS,1,7850          !1ϵܶDENS

!(2)ڵλ,ڵ
!ϱڵλòڵ
*DO,i,1,(m+1),
*DO,j,(2*n+1)*(i-1)+1,(2*n+1)*i
*IF,i,le,m,then
*IF,j,le,n*i+(i-1)*(n+1),then
N,j,(2*(j-(i-1)*(2*n+1))-1)*lx*sin(45),ly*sin(45)*2*(i-1),0
*ELSE
N,j,(j-(n+1)*(i-1)-n*i-1)*lx*2*sin(45),ly*(2*i-1)*sin(45),0
*ENDIF
*ELSE
*IF,j,le,n*i+(i-1)*(n+1),then
N,j,(2*(j-(i-1)*(2*n+1))-1)*lx*sin(45),ly*2*sin(45)*(i-1),0
*ENDIF
*ENDIF
*ENDDO
*ENDDO
!±ڵλòڵ
*DO,i,1,m
*DO,j,n*(m+1)+(n+1)*m+n*(i-1)+1, n*(m+1)+(n+1)*m+n*i
N,j,(2*(j-(n*(m+1)+(n+1)*m+n*(i-1)))-1)*lx*sin(45),(2*i-1)*ly* sin(45),z
*ENDDO
*ENDDO
!(4)嵥Ԫ
TYPE,1                   !õԪ
MAT,1                       !ò 
REAL,1                    !ʵ
!ϱҸ
*DO,i,1,m
*DO,j,1+n*(i-1)+(n+1)*(i-1),n*i+(n+1)*(i-1)
E,j,j+n
*ENDDO
*ENDDO
*DO,i,1,m
*DO,j,1+n*(i-1)+(n+1)*(i-1),n*i+(n+1)*(i-1)
E,j,j+n+1
*ENDDO
*ENDDO
*DO,i,1,m
*DO,j,1+n*i+(n+1)*(i-1),n*(i+1)+(n+1)*(i-1)
E,j,j+n+1
*ENDDO
*ENDDO
*DO,i,1,m
*DO,j,2+n*i+(n+1)*(i-1),n*i+(n+1)*i
E,j,j+n
*ENDDO
*ENDDO
*DO,i,1,n-1,
  E,i,i+1
*ENDDO
*DO,i,n+1, ((n+1)*(m-1)+n*m),2*n+1
E,i,i+2*n+1
*ENDDO
*DO,i,2*n+1,((n+1)*(m-1)+n*m),2*n+1
E,i,i+2*n+1
*ENDDO
*DO,i, (n+1)*m+n*m+1,n*(m+1)+(n+1)*m-1,
  E,i,i+1
*ENDDO
TYPE,1                   !õԪ
MAT,1                       !ò 
REAL,2                   !ʵ
!±Ҹ
*DO,i,1,m
*DO,j,n*(m+1)+(n+1)*m+n*(i-1)+1, n*(m+1)+(n+1)*m+n*i-1
E,j,j+1
*ENDDO
*ENDDO
*DO,i,1,m-1
*DO,j,n*(m+1)+(n+1)*m+n*(i-1)+1, n*(m+1)+(n+1)*m+n*i
E,j,j+n
*ENDDO
*ENDDO
TYPE,1                   !õԪ
MAT,1                       !ò 
REAL,3                   !ʵ
!б
*DO,i,1,m
*DO,j,n*(m+1)+(n+1)*m+n*(i-1)+1,n*(m+1)+(n+1)*m+n*i
E,j,j-(n*(m+1)+(n+1)*m-(n+1)*(i-1))
*ENDDO
*ENDDO
*DO,i,1,m
*DO,j,n*(m+1)+(n+1)*m+n*(i-1)+1,n*(m+1)+(n+1)*m+n*i
E,j,j-(n*(m+1)+(n+1)*m-n-(n+1)*(i-1))
*ENDDO
*ENDDO
*DO,i,1,m
*DO,j,n*(m+1)+(n+1)*m+n*(i-1)+1,n*(m+1)+(n+1)*m+n*i
E,j,j-(n*(m+1)+(n+1)*m-n-1-(n+1)*(i-1))
*ENDDO
*ENDDO
*DO,i,1,m
*DO,j,n*(m+1)+(n+1)*m+n*(i-1)+1,n*(m+1)+(n+1)*m+n*i
E,j,j-(n*(m+1)+(n+1)*m-(n+1)*i-n)
*ENDDO
*ENDDO
*DO,i,1,m
*DO,j,n*(m+1)+(n+1)*m+n*(i-1)+1, n*(m+1)+(n+1)*m+n*i
F,j,fz,-30000
*ENDDO
*ENDDO
*DO,i,1,n,
  d,i,uy,0
  d,i,uz,0
*ENDDO
*DO,i,n+1,((n+1)*m+n*m),2*n+1
  D,i,ux,0
  D,i,uz,0
*ENDDO
*DO,i,2*n+1,((n+1)*m+n*m),2*n+1
D,i,ux,0
D,i,uz,0
*ENDDO
*DO,i, (n+1)*m+n*m+1,n*(m+1)+(n+1)*m,
  D,i,uy,0
  D,i,uz,0
*ENDDO
/SOLU
ACEL,,,9.8           !ʩٶȣԿǸ˼
ALLSELALL
SOLVE
FINISH              !˳ģ
/POST1              !ͨú
ETABLE,EVOLUME,VOLU        !嵥Ԫб
SSUM                       !Ԫб
*GET,VTOT,SSUM,,ITEM,EVOLUME   !volume洢ܽṹĲ
RO=7850
WT=VTOT*RO
ETABLE,SMAX_E,LS, 1 
ESORT,ETAB,SMAX_E,0,1, ,
*GET,SMAXE,SORT,,MAX
*GET,DZMAX,NODE,55,u,z    
*SET,DMAX,-1*DZMAX        
LGWRITE,'OPT1','lgw','e:\aa',comment   !Żļ
FINISH                       !˳ͨúģ
/OPT                         !Ż
OPANL,'opt1','lgw',''!ָŻļ
OPVAR,D1,DV,0.02,0.07, 
OPVAR,D2,DV,0.03,0.06,
OPVAR,D3,DV,0.01,0.05,
OPVAR,DMAX,SV,0.01,0.07,     
OPVAR,SMAXE,SV,180000000,210000000,
OPSAVE,'OPT1','opt',' '       !洢Żݿ
OPVAR,WT,OBJ,,,     !Ŀwt
OPTYPE,FIRS             !ָŻ
OPFRST,100,100,0.2       !Ż趨
OPEXE          !Ż
*STATUS         !鿴Ż
