!˳Ϊµڶ3.2.1س
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     !س߱߷ɵ1~mнڵ
*ENDDO
*ENDDO
 !±ڵλòڵ
*DO,i,2,m-1,
*DO,j,n*m+1+(i-2)*(n-2),(n-2)*(i-1)+n*m
  N,j,(j-n*m-(i-2)*(n-2))*lx,ly*(i-1),z   !س߱߷ɵ2~m-1нڵ
*ENDDO
*ENDDO
 !(3)嵥Ԫͼʵ
RTNUM=3     !ʵ
*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'
MULTIPRO,'START',1        !ѡǷʵR2
   *CSET,1,3,addreal,'Add real constant 2? YES(1)/NO(0)',0    !ĬΪR2
   *CSET,61,62,'whether add real constant 2'
MULTIPRO,'END'
*IF,addreal,eq,1,then      !ѡʵ
    MULTIPRO,'START',10    !ɵڶʵR2Ի
        *CSET,1,3,AREAIN(2),'Section Area (mm^2), AREA(2) =',1382e-6
        *CSET,4,6,IZZIN(2),'Area moment of inertia, IZZ(2) =',2093500e-12
        *CSET,7,9,IYYIN(2),'Area moment of inertia, IYY(2) =',2093500e-12
        *CSET,10,12,TKZIN(2),'Thinkness along Z axis, TKZ(2) =',114e-3
        *CSET,13,15,TKYIN(2),'Thinkness along Y axis, TKY(2) =',114e-3
        *CSET,16,18,THETAIN(2),'Orientation about X axis, THEAT(2) =',0
        *CSET,19,21,ISTRNIN(2),'Initial strain, ISTRN(2) =',0
        *CSET,22,24,IXXIN(2),'Torsional moment of inertia, IXX(2) =', 0
        *CSET,25,27,SHEARZIN(2),'Shear deflection const Z, SHEARZ(2) =',0
        *CSET,28,30,SHEARYIN(2),'Shear deflection const Y, SHEARY(2) =',0
        *CSET,61,62,'Input Real Constant 2:'
        *CSET,63,64,'Just input AREA if link8 was used:'
    MULTIPRO,'END'
    MULTIPRO,'START',1           !ѡǷʵR3
         *CSET,1,3,addreal,'Add real constant 3? YES(1)/NO(0)',0
         *CSET,61,62,'whether add real constant 3'
    MULTIPRO,'END'
      *IF,addreal,eq,1,then       !ѡʵR3
           MULTIPRO,'START',10    !ɵڶʵR3Ի
             *CSET,1,3,AREAIN(3),'Section Area (mm^2), AREA(3) =',1382e-6
             *CSET,4,6,IZZIN(3),'Area moment of inertia, IZZ(3) =',2093500e-12
             *CSET,7,9,IYYIN(3),'Area moment of inertia, IYY(3) =',2093500e-12
             *CSET,10,12,TKZIN(3),'Thinkness along Z axis, TKZ(3) =',114e-3
             *CSET,13,15,TKYIN(3),'Thinkness along Y axis, TKY(3) =',114e-3
             *CSET,16,18,THETAIN(3),'Orientation about X axis, THEAT(3) =',0
             *CSET,19,21,ISTRNIN(3),'Initial strain, ISTRN(3) =',0
             *CSET,22,24,IXXIN(3),'Torsional moment of inertia, IXX(3) =', 0
             *CSET,25,27,SHEARZIN(3),'Shear deflection const Z, SHEARZ(3) =',0
             *CSET,28,30,SHEARYIN(3),'Shear deflection const Y, SHEARY(3) =',0
             *CSET,61,62,'Input Real Constant 2:'
             *CSET,63,64,'Just input AREA if link8 was used:'
           MULTIPRO,'END'
            !Ҫʵɼ϶
       *ENDIF
*ENDIF                         !addreal3Ϊ3Ӧ
*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                        !ض̱߷1~nҸ
*ENDDO
*DO,i,1,m,
*DO,j,(i-1)*n+1,i*n-1,
 E,j,j+1                        !س߷1~mҸ
*ENDDO
*ENDDO 
 !±Ҹ
*DO,i,2,m-2,
*DO,j,n*m+1+(i-2)*(n-2),(n-2)*(i-1)+n*m
 E,j,j+n-2                       !ض̱߷2~n-1Ҹ
*ENDDO
*ENDDO
*DO,i,2,m-1,
*DO,j,n*m+1+(i-2)*(n-2),(n-2)*(i-1)+n*m-1,
 E,j,j+1                        !س߷2~m-1Ҹ
*ENDDO
*ENDDO 
 !
*DO,i,2,m-1,
*DO,j,(i-1)*n+2,i*n-1,
 E,j,j+n*m+1+(i-2)*(n-2)-(i-1)*n-2              !
*ENDDO
*ENDDO 
 !б
*DO,i,2,m-1,
*DO,j,(i-1)*n+1,n*i-(n-1)/2-1,
 E,j,j+n*m+1+(i-2)*(n-2)-(i-1)*n-1  !س߷2~m-1еб
*ENDDO
*ENDDO
*DO,i,2,m-1,
*DO,j,n*i-(n-1)/2+1,n*i,
 E,j,j+n*m+(i-2)*(n-2)+n-2-n*I      !س߷2~m-1еҰб
*ENDDO
*ENDDO
*DO,i,1,(m-1)/2,
*DO,j,n*(i-1)+2,n*i-1
 E,j,j+n*m-1+(i-1)*(n-2)-(i-1)*n   !ض̱߷2~m-1е°б
*ENDDO
*ENDDO
*DO,i,(m-1)/2+1,m-1,
*DO,j,n*i+2,n*(i+1)-1
 E,j,j+n*m+1+(i-2)*(n-2)-i*n-2     !ض̱߷2~m-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-2)*(m-2),
  F,i,fz,-p                         !ҽڵʩӼк
*ENDDO                     
*DO,i,1,n,
  D,i,all,0                   !Լس߷һнڵuxuyuzɶ
*ENDDO
*DO,i,n*(m-1)+1,n*m,
  D,i,all,0                  !Լس߷mнڵuxuyuzɶ
*ENDDO
*DO,i,n+1,(n*(m-2)+1),n
   D,i,all,0    !Լض̱߷һнڵ㣨ϵĽڵ⣩uxuyuzɶ
*ENDDO
*DO,i,n*2,n*(m-1),n
   D,i,all,0    !Լض̱߷nнڵ㣨ϵĽڵ⣩uxuyuzɶ*ENDDO
!(6)ģļ
SAVE,MODE,DB,,ALL
