!˳Ϊµڶ3.2.2س
/PREP7
 !(1) û,β
MULTIPRO,'START',5                       !ɶʾԻ
*CSET,1,3,lx,'Grid side lx= (m)',1.2     !x귽45oҸ˳
*CSET,4,6,ly,'Grid side ly= (m)',1      !x귽135oҸ˳
*CSET,7,9,n,'grid number of the long side n the top surface =',6  
 !Ĭƽx귽ڵ
*CSET,10,12,m,'grid number of the short side n the top surface =',6  
 !Ĭƽy귽ڵ
*CSET,13,15,z,'distance between the top and the bottom surface z=(m)',-1  
 !ƽƽ֮Ĵֱ
MULTIPRO,'END'
 !(2) ڵλ,ڵ     
 !ƽڵλòڵ
*DO,i,1,m,
*DO,j,(2*n-1)*(i-1)+1,(2*n-1)*i
*IF,i,le,m-1,then
*IF,j,le,n*i+(i-1)*(n-1),then
N,j,(j-1-(2*n-1)*(i-1))*lx*2*sin(45),ly*2*sin(45)*(i-1),0 
!x귽ӵ1У2m-1г⣩ڵ
*ELSE
n,j,(2*(j-n*i-(i-1)*(n-1))-1)*lx*2*sin(45)/2,ly*2*sin(45)*(2*i-1)/2,0
!x귽ӵ2żнڵ
*ENDIF
*ELSE
*IF,j,le,n*i+(i-1)*(n-1),then
N,j,(j-1-(2*n-1)*(i-1))*lx*2*sin(45),ly*2*sin(45)*(i-1),0
!x귽2m-1нڵ
*ENDIF
*ENDIF
*ENDDO
*ENDDO
 !ƽڵλòڵ
*DO,i,1,m-1,
*DO,j,(n-1)*(m-1)+n*m+1+(2*n-3)*(i-1), (n-1)*(m-1)+n*m+1+(2*n-3)*i
*IF,i,le,(m-2),then
*IF,j,le,(n-1)*(m-1)+n*m+(n-1)*i+(i-1)*(n-2),then
N,j,(2*(j-(n-1)*(m-1)-n*m-(i-1)*(2*n-3))-1)*lx*2*sin(45)/2,ly*2*sin(45)*(2*i-1)/2,z
!x귽ӵ2пʼżУ2m-2г⣩ڵ
*ELSE
N,j,(j-(n-1)*(m-1)-n*m-(n-1)*i-(i-1)*(n-2))*lx*2*sin(45),ly*2*sin(45)*i,z
!x귽ӵ3пʼУ2m-1г⣩ڵ
*ENDIF
*ELSE
*IF,j,le,(n-1)*(m-1)+n*m+(n-1)*i+(i-1)*(n-2),then
N,j,(2*(j-(n-1)*(m-1)-n*m-(i-1)*(2*n-3))-1)*lx*2*sin(45)/2,ly*2*sin(45)*(2*i-1)/2,z
!x귽ӵ2m-2нڵ
*ENDIF
*ENDIF
*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) =',1382e-6
    *CSET,4,6,IZZIN(1),'Area moment of inertia, IZZ(1) =',2093500e-12
    *CSET,7,9,IYYIN(1),'Area moment of inertia, IYY(1) =',2093500e-12
    *CSET,10,12,TKZIN(1),'Thinkness along Z axis, TKZ(1) =',114e-3
    *CSET,13,15,TKYIN(1),'Thinkness along Y axis, TKY(1) =',114e-3
    *CSET,16,18,THETAIN(1),'Orientation about X axis, THEAT(1) =',0
    *CSET,19,21,ISTRNIN(1),'Initial strain, ISTRN(1) =',0
    *CSET,22,24,IXXIN(1),'Torsional moment of inertia, IXX(1) =', 0
    *CSET,25,27,SHEARZIN(1),'Shear deflection const Z, SHEARZ(1) =',0
    *CSET,28,30,SHEARYIN(1),'Shear deflection const Y, SHEARY(1) =',0
    *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-1
*DO,j,1+(2*n-1)*(i-1),n-1+(2*n-1)*(i-1)
E,j,j+n
*ENDDO
*ENDDO
*DO,i,1,m-1
*DO,j,2+(2*n-1)*i,n+(2*n-1)*i      
E,j,j-n
*ENDDO
*ENDDO                                     !x귽45oҸ
*DO,i,1,m-1
*DO,j,2+(2*n-1)*(i-1),n+(2*n-1)*(i-1)
E,j,j+n-1
*ENDDO
*ENDDO
*DO,i,1,m-1
*DO,j,1+(2*n-1)*i,n-1+(2*n-1)*i
E,j,j-n+1
*ENDDO
*ENDDO                                    !x귽135oҸ
*DO,i,1,n-1,
E,i,i+1
*ENDDO                                   !x귽1Ҹ
*DO,i,1,(n*(m-1)+(n-1)*(m-1)),2*n-1
E,i,i+2*n-1
*ENDDO                                   !y귽1Ҹ
*DO,i,n,(n*(m-1)+(n-1)*(m-1)),2*n-1
E,i,i+2*n-1
*ENDDO                                   !y귽2n-1Ҹ
*DO,i,n*(m-1)+(n-1)*(m-1)+1,n*m+(n-1)*(m-1)-1,
E,i,i+1
*ENDDO                                    !x귽2m-1Ҹ
 !ƽҸ
*DO,i,1,(m-2)
*DO,j,m*n+(m-1)*(n-1)+1+(2*n-3)*(i-1), m*n+(m-1)*(n-1)+n-2+(2*n-3)*(i-1)
E,j,j+n-1
*ENDDO
*ENDDO
*DO,i,1,(m-2)
*DO,j,m*n+(m-1)*(n-1)+2+(2*n-3)*i,m*n+(m-1)*(n-1)+n-1+(2*n-3)*i
E,j,j-n+1
*ENDDO
*ENDDO                                   !x귽45oҸ
*DO,i,1,(m-2)
*DO,j,m*n+(m-1)*(n-1)+2+(2*n-3)*(i-1), m*n+(m-1)*(n-1)+n-1+(2*n-3)*(i-1)
E,j,j+n-2
*ENDDO
*ENDDO
*DO,i,1,(m-2)
*DO,j,m*n+(m-1)*(n-1)+1+(2*n-3)*i, m*n+(m-1)*(n-1)+n-2+(2*n-3)*i
E,j,j-n+2
*ENDDO
*ENDDO                                     !x귽135oҸ
!
*DO,i,1,m-1,
*DO,j,((m*n)+(n-1)*(m-1)+1+(2*n-3)*(i-1)),((m*n)+(n-1)*(m-1)+(2*n-3)*i)
*IF,i,le,(m-2),then
*IF,j,le,((m*n)+(n-1)*(m-1)+n-1+(2*n-3)*(i-1)),then
E,j,j-(m*n+(n-1)*(m-1)-n-(2*i-2))
!x귽2пʼżУ2m-2г⣩
*ELSE
E,j,j-(m*n+(n-1)*(m-1)-n-(2*i-1))
!x귽3пʼУ2m-1г⣩
*ENDIF
*ELSE
*IF,j,le,((m*n)+(n-1)*(m-1)+n-1+(2*n-3)*(i-1)),then
E,j,j-(m*n+(n-1)*(m-1)-n-(2*i-2))
!x귽2m-2
*ENDIF
*ENDIF
*ENDDO
*ENDDO 
 !б
 !x귽135oб
*DO,i,1,m-1
*DO,j,n*(i-1)+(n-1)*(i-1)+1,n*i+(n-1)*(i-1)
*IF,j,gt,n*(i-1)+(n-1)*(i-1)+1+i-1,then
 E,j,j+m*n+(m-1)*(n-1)-(2*i-1)
*ENDIF
*ENDDO
*ENDDO
*DO,i,1,m-2
*DO,j,n*i+(n-1)*(i-1)+1,n*i+(n-1)*i
*IF,j,gt,n*i+(n-1)*(i-1)+1+i-1,then
 E,j,j+m*n+(m-1)*(n-1)-2*i
*ENDIF
*ENDDO
*ENDDO
*DO,i,1,m
*DO,j,n*(i-1)+(n-1)*(i-1)+1,n*i+(n-1)*(i-1)
*IF,j,lt,n*(i-1)+(n-1)*(i-1)+1+i-1,then
 E,j,j+(m*n+(m-1)*(n-1)+1-2*n-2*(i-2))
 *ENDIF
*ENDDO
*ENDDO
*DO,i,1,m-1
*DO,j,n*i+(n-1)*(i-1)+1,n*i+(n-1)*i
*IF,j,lt,n*i+(n-1)*(i-1)+1+i-1,then
 E,j,j+(m*n+(m-1)*(n-1)+1-2*n-(2*i-3))
*ENDIF
*ENDDO
*ENDDO
!x귽45oб
*DO,i,1,m-1
*DO,j,n*(i-1)+(n-1)*(i-1)+1,n*i+(n-1)*(i-1)
*IF,j,lt,n*i+(n-1)*(i-1)-(i-1),then
 E,j,j+m*n+(m-1)*(n-1)-(2*i-2)
*ENDIF
*ENDDO
*ENDDO
*DO,i,1,m-2
*DO,j,n*i+(n-1)*(i-1)+1,n*i+(n-1)*i
*IF,j,lt,n*i+(n-1)*i-(i-1),then
 E,j,j+m*n+(m-1)*(n-1)-(2*i-1)
*ENDIF
*ENDDO
*ENDDO
*DO,i,1,m
*DO,j,n*(i-1)+(n-1)*(i-1)+1,n*i+(n-1)*(i-1)
*IF,j,gt,n*i+(n-1)*(i-1)-(i-1),then
E,j,j+m*n+(m-1)*(n-1)-2*n-(2*i-4)
*ENDIF
*ENDDO
*ENDDO
*DO,i,1,m-1
*DO,j,n*i+(n-1)*(i-1)+1,n*i+(n-1)*i
*IF,j,gt,n*i+(n-1)*i-(i-1),then
 E,j,j+m*n+(m-1)*(n-1)-2*n-(2*i-3)
*ENDIF
*ENDDO
*ENDDO
 !(5)߽Լͽڵغ
MULTIPRO,'START',1
    *CSET,1,3,P,'Loads(N)',200
    *CSET,61,62,'please input loads value:'
MULTIPRO,'END'
*DO,i,(n-1)*(m-1)+m*n+1,(n-1)*(m-1)+m*n +(m-1)*(n-1)+(n-2)*(m-2),
  F,i,fz,-P                                  !ƽڵܼк
*ENDDO
*DO,i,1,n,
  D,i,all,0
*ENDDO               !ԼX귽һнڵUXUYUZɶ
*DO,i,n*(m-1)+(n-1)*(m-1)+1,n*m+(n-1)*(m-1),
  D,i,all,0             
*ENDDO               !ԼX귽2M-1нڵUXUYUZɶ
*DO,i,1,(n*(m-1)+(n-1)*(m-1)),2*n-1
*IF,i,gt,1,then
D,i,all,0                !ԼY귽һʣڵUXUYUZɶ
*ENDIF
*ENDDO
*DO,i,n,(n*(m-1)+(n-1)*(m-1)),2*n-1
*IF,i,gt,n,then
D,i,all,0             !ԼY귽2N-1ʣڵUXUYUZɶ
*ENDIF
*ENDDO
!(6)ģļ
SAVE,MODE,DB,,ALL
