!˳Ϊµڶ3.2.4س
/PREP7
*AFUN,DEG                   !ǶȵλΪ
!(1)ûƣβ
MULTIPRO,'START',5                !ɶʾԻ
*CSET,1,3,lx,'Grid side lx= (m)',1  !Ĭ񳤱߱߳
*CSET,4,6,ly,'Grid side ly= (m)',1  !Ĭ̱߱߳
*CSET,7,9,n,'grid number of the long side n the top surface =',8
!Ĭƽ񳤱߷ڵ
*CSET,10,12,m,'grid number of the short side n the top surface =',6
!Ĭƽ̱߷ڵ 
*CSET,13,15,z,'distance between the top and the bottom surface z=(m)',-1 
!ϡƽ֮Ĵֱ
*CSET,61,62,'Please input geometry parameters of the grid.'   
!ڶԻ򶥲ʾʾ Please input geometry parameters of the grid
MULTIPRO,'END'
!(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
!س߷У2m+1г⣩Ľڵ
*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
!س߷ĵ2m+1еĽڵ
*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
!(3)嵥Ԫͼʵ
RTNUM=1      !ʵΪʡƪֻһʵͣ
*DIM,AREAIN,ARRAY,RTNUM       !˼
*DIM,IZZIN,ARRAY, RTNUM       !Ծ
*DIM,IYYIN,ARRAY, RTNUM       !Ծ
*DIM,TKZIN,ARRAY, RTNUM       !Z߶
*DIM,TKYIN,ARRAY, RTNUM       !Y߶
*DIM,THETAIN,ARRAY, RTNUM     !˼Xн
*DIM,ISTRNIN,ARRAY, RTNUM     !ʼӦ
*DIM,IXXIN,ARRAY, RTNUM       !Ť
*DIM,SHEARZIN,ARRAY, RTNUM    !бγ
*DIM,SHEARYIN,ARRAY, RTNUM    !бγ
MULTIPRO,'START', 1           !ɵԪѡԻ
    *CSET,1,3,eltype,'Link8 Input 0, Beam4 Input 1',0    !ѡԪͣ0Ϊѡlink81Ϊѡbeam4Ҳͨı䱾µĵԪ͡
    *CSET,61,62,'Element Type 1 Option:'            
MULTIPRO,'END'
*IF,eltype,EQ,0,THEN       !0eltype=0
     ET,1,link8            !1൥ԪΪ˵Ԫlink8
*ENDIF
*IF,eltype,EQ,1,THEN       !1eltype=1
     ET,1,beam4            !1൥ԪΪԪbeam4
*ENDIF
MULTIPRO,'START',10        !ɵһʵR1Ի
    *CSET,1,3,AREAIN(1),'Section Area (mm^2), AREA(1) =',
    *CSET,4,6,IZZIN(1),'Area moment of inertia, IZZ(1) =',
    *CSET,7,9,IYYIN(1),'Area moment of inertia, IYY(1) =',
    *CSET,10,12,TKZIN(1),'Thinkness along Z axis, TKZ(1) =',
    *CSET,13,15,TKYIN(1),'Thinkness along Y axis, TKY(1) =',
    *CSET,16,18,THETAIN(1),'Orientation about X axis, THEAT(1) =',
    *CSET,19,21,ISTRNIN(1),'Initial strain, ISTRN(1) =',
    *CSET,22,24,IXXIN(1),'Torsional moment of inertia,IXX(1) =',
    *CSET,25,27,SHEARZIN(1),'Shear deflection const Z, SHEARZ(1) =',
    *CSET,28,30,SHEARYIN(1),'Shear deflection const Y, SHEARY(1) =',
    *CSET,61,62,'Input Real Constant 1:'
    *CSET,63,64,'Just input AREA if link8 was used:'   !ѡ˵Ԫֻ
MULTIPRO,'END'
*DO,i,1, RTNUM                !鶨嵥Ԫʵ
*IF,eltype,EQ,0,THEN          !ѡǸ˵Ԫ,eltype=0
     R,i,AREAIN(i)            !˵Ԫ
*ENDIF
*IF,eltype,EQ,1,THEN          !ѡԪ,eltype=1
R,i,AREAIN(i),IZZIN(i),IYYIN(i),TKZIN(i),TKYIN(i),THETAIN(i),ISTRNIN(i),IXXIN(i),SHEARZIN(i),SHEARYIN(i),,   !Ԫʵ
*ENDIF
*ENDDO
MULTIPRO,'START',3            !ɲԻ
    *CSET,1,3,EXIN,'Elasticity (GPa), EX =',2.1e11      !뵯ģEX
    *CSET,4,6,PRXYIN,'Passion Ratio, PRXY =',0.3        !벴ɱPRXY
    *CSET,7,9,DENSIN,'Dens of steel, DENS =',7850       !ϵܶDENS
    *CSET,61,62,'Input Material Property:'
MULTIPRO,'END'
MP,EX,1,EXIN                 !1ϵĵģEX
MP,PRXY,1,PRXYIN             !1ϵĲɱPRXY
MP,DENS,1,DENSIN             !1ϵܶDENS
!(4)嵥Ԫ
MULTIPRO,'START',3            !ɲԻ
    *CSET,1,3,TYPEIN,'Eelment Type Number',1      !뵥Ԫͺ
    *CSET,4,6,MATIN,'Material Type Number',1      !Ժ
    *CSET,7,9,REALIN,'Real Constance Number',1    !ʵ
    *CSET,61,62,'Input Material Type Number:'
MULTIPRO,'END'
TYPE,TYPEIN                       !õԪ
MAT,MATIN                         !ò 
REAL,REALIN                       !ʵ
!ƽҸ
*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,2+n*i+(n+1)*(i-1),n*i+(n+1)*i
E,j,j+n
*ENDDO
*ENDDO                              !볤߷135oҸ
*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                          !볤߷45oҸ
*DO,i,1,n-1,
  E,i,i+1
*ENDDO                          !ƽس߷1Ҹ
*DO,i,n+1, ((n+1)*(m-1)+n*m),2*n+1
E,i,i+2*n+1
*ENDDO                          !ƽض̱߷1Ҹ
*DO,i,2*n+1,((n+1)*(m-1)+n*m),2*n+1
E,i,i+2*n+1
*ENDDO                          !ƽض̱߷n +1Ҹ
*DO,i, (n+1)*m+n*m+1,n*(m+1)+(n+1)*m-1,
  E,i,i+1
*ENDDO                          !ƽس߷2m +1Ҹ
!ƽҸ
*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                            !ض̱߷Ҹ
!бˣĸбˣ
*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
!(5)߽Լͽڵغɣƽܱ߽ڵ½ӣ
MULTIPRO,'START',1
    *CSET,1,3,P,'Loads(N)',30000
    *CSET,61,62,'please input loads value:'
MULTIPRO,'END'
*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,-P
*ENDDO
*ENDDO                          !ҽڵʩӼк
*DO,i,1,n,
  D,i,uy,0
  D,i,uz,0
*ENDDO                    !Լƽس߷1нڵUYUZɶ
*DO,i,n+1,((n+1)*m+n*m),2*n+1
  D,i,ux,0
  D,i,uz,0
*ENDDO                   !Լƽض̱߷1нڵUXUZɶ
*DO,i,2*n+1,((n+1)*m+n*m),2*n+1
D,i,ux,0
 D,i,uz,0
*ENDDO                  !Լƽض̱߷n +1нڵUXUZɶ
*DO,i, (n+1)*m+n*m+1,n*(m+1)+(n+1)*m,
  D,i,uy,0
  D,i,uz,0
*ENDDO                !Լƽس߷2m +1нڵUYUZɶ
!(6)ģļ
SAVE,MODE,DB,,ALL
3.3ܽṹ
