!˳Ϊµڶ3.2.3س
FINISH
/CLEAR
/PREP7
!(1)ûƣβ
MULTIPRO,'START',5                  !ɶʾԻ
*CSET,1,3,lx,'Grid side lx= (m)',1         !Ĭ񳤱߷Ҹ˳
*CSET,4,6,ly,'Grid side ly= (m)',1.2        !Ĭ̱߷Ҹ˳
*CSET,7,9,n,'grid number of the long side n the top surface =',7  
!Ĭƽ񳤱߷ڵ
 *CSET,10,12,m,'grid number of the short side n the top surface =',7  
 !Ĭƽ̱߷ڵ 
*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,
*DO,j,n*(i-1)+1,n*i
  N,j,(j-1-n*(i-1))*lx,ly*(i-1),0
*ENDDO
*ENDDO
!ƽڵλòڵ
*DO,i,1,m-1,
*DO,j,n*m+1+(i-1)*(n-1),(n-1)*i+n*m
  N,j,(2*(j-n*m-(i-1)*(n-1))-1)/2*lx,(2*i-1)/2*ly,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',1   !ѡԪͣ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,j,1,(m-1)*n,
  E,j,j+n
*ENDDO                           !ض̱߷Ҹ
*DO,i,1,m,
*DO,j,(i-1)*n+1,i*n-1,
 E,j,j+1
*ENDDO
*ENDDO                            !س߷Ҹ
!ƽҸ
*DO,i,1,m-2,
*DO,j,n*m+1+(i-1)*(n-1),(n-1)*i+n*m
 E,j,j+n-1                            
*ENDDO
*ENDDO                          !ض̱߷Ҹ
*DO,i,1,m-1,
*DO,j,n*m+1+(i-1)*(n-1),(n-1)*i+n*m-1,
 E,j,j+1                               
*ENDDO
*ENDDO                          !س߷Ҹ
!бˣĸбˣ
*DO,i,1,m-1,
*DO,j,(i-1)*n+1,n*(i-1)+n-1
 E,j,j+n*m-(i-1)
*ENDDO
*ENDDO
*DO,i,1,m-1,
*DO,j,(i-1)*n+2,n*(i-1)+n
 E,j,j+n*m-1-(i-1)
*ENDDO
*ENDDO
*DO,i,1,m-1,
*DO,j,i*n+1,n*i+n-1
 E,j,j+n*m-(i-1)-n
*ENDDO
*ENDDO
*DO,i,1,m-1,
*DO,j,i*n+2,n*i+n
 E,j,j+n*m-(i-1)-n-1
*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*m+1,n*m+(n-1)*(m-1),
  F,i,fz,-P         
*ENDDO                          !ƽнڵʩӼк
*DO,i,1,n,
  D,i,all,0            !Լƽس߷1нڵuxuyuzɶ
*ENDDO
*DO,i,n*(m-1)+1,n*m,
  D,i,all,0
*ENDDO             !Լƽس߷mнڵuxuyuzɶ
*DO,i,n+1,(n*(m-2)+1),n
   D,i,all,0
*ENDDO            !Լƽض̱߷1ʣڵuxuyuzɶ
*DO,i,n*2,n*m,n
   D,i,all,0
*ENDDO            !Լƽض̱߷nʣڵuxuyuzɶ

!(6)ģļ
SAVE,MODE,DB,,ALL
