1、finish/clear/Prep7! basic geometrical parametersmultipro,start,5*cset,1,3,lx,Grid side lx= (m),3.3*cset,4,6,ly,Grid side ly= (m),3.3*cset,7,9,n,Grid number of the long side on the top surface =,10*cset,10,12,m,Grid number of the short side on the top surface =,10*cset,13,15,z,distance between the to
2、p and the bottom surface z=(m),-2.3*cset,61,62,Please input geometry parameters of the grid.multipro,end! nodes definition*do,i,1,m*do,j,n*(i-1)+1,n*i n,j,(j-1-n*(i-1)*lx,(i-1)*ly,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)*lx/2,(2*i-1)*ly/2,z*enddo*enddo
3、!element type and real constantsRTNUM=1*DIM,areain,array,rtnum*DIM,izzin,array,rtnum*DIM,iyyin,array,rtnum*DIM,tkzin,array,rtnum*DIM,tkyin,array,rtnum*DIM,thetain,array,rtnum*DIM,istrnin,array,rtnum*DIM,ixxin,array,rtnum*DIM,shearzin,array,rtnum*DIM,shearyin,array,rtnummultipro,start,1*cset,1,3,elty
4、pe,Link8 Input 0, Beam4 Input 1,0*cset,61,62,Element Type 1 Option:multipro,end*if,eltype,eq,0,then et,1,link8*endif*if,eltype,eq,1,then et,1,beam4*endifmultipro,start,10*cset,1,3,areain(1),Section Area (mm2), AREA(1)=,*cset,4,6,izzin(1),Area moment of inertia(mm4), IZZ(1)=,*cset,7,9,iyyin(1),Area m
5、oment of inertia(mm4), IYY(1)=,*cset,10,12,tkzin(1),Thickness along Z axis, TKZ(1)=,*cset,13,15,tkyin(1),Thickness 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(mm4),
6、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,Please 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 r,i,areain(i)*endif*if,eltype,eq,1,th
7、en r,i,areain(i),izzin(i),iyyin(i),tkzin(i),tkyin(i),thetain(i), rmore,istrnin(i),ixxin(i),shearzin(i),shearyin(i),*endif*enddomultipro,start,3*cset,1,3,exin,Elasticity (GPa), EX=,2.1E11*cset,4,6,prxyin,Poisson Ratio, PRXY=,0.3*cset,7,9,densin,Dens of steel, DENS(1)=,7850*cset,61,62,Please input mat
8、erial property:multipro,endmp,ex,1,exinmp,prxy,1,prxyinmp,dens,1,densin! multipro,start,3*cset,1,3,typein,Element type number,1*cset,4,6,matin,Material type number,1*cset,7,9,realin,Real constant type number,1*cset,61,62,Please input element type number:multipro,endtype,typeinmat,matinreal,realin!*d
9、o,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,(i-1)*n+n-1 e,j,j+n*m-(i-1)*enddo*enddo*do,i,1,m-1*do,j
10、,(i-1)*n+2,(i-1)*n+n e,j,j+n*m-(i-1)-1*enddo*enddo*do,i,1,m-1*do,j,i*n+1,i*n+n-1 e,j,j+n*m-(i-1)-n*enddo*enddo*do,i,1,m-1*do,j,i*n+2,i*n+n e,j,j+n*m-(i-1)-n-1*enddo*enddo!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*enddo*do,i,n*(m-1)+1,n*m d,i,all,0*enddo*do,i,n+1,(n*(m-2)+1),n d,i,all,0*enddo*do,i,n*2,n*(m-1),n d,i,all,0*enddo!save,model,db,all
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1