<P>一个combin39单元应用实例</P>
<P>橡胶垫选用combin39; <BR>若考虑叠层橡胶垫隔震,则可考虑combin14,该算例全部体现; <BR><BR>/prep7 <BR>local,11,1 <BR>cscir,11,1 <BR>f=8.892 <BR>h=2 <BR>w=39 <BR>r=(f*f+(w*w/4))/(2*f) <BR>b=(atan(2*(r-f)/w))*180/3.14159 <BR>n,1,r,b <BR>n,14,r,180-b <BR>fill,1,14 <BR>/view,1,1,1,1 <BR>/ang,1 <BR>/rep,fast <BR>local,11, <BR>ngen,22,14,1,14,1,,,3 <BR>local,11,1 <BR>n,309,r-h,(90+12*b)/13,1.5 <BR>ngen,13,1,309,309,,,(180-2*b)/13 <BR>local,11 <BR>ngen,21,13,309,321,1,,,3 <BR>et,1,link8 <BR>mp,ex,1,2.06e11 <BR>mp,dens,1,7890 <BR>mp,nuxy,1,0.3 <BR>MP,DAMP,1,0.02 <BR>r,1,703e-6 <BR>type,1 <BR>mat,1 <BR>real,1 <BR>*do,i,1,14,1 <BR>e,i,i+14 <BR>*enddo <BR>egen,21,14,1,14,1 <BR>*do,i,1,13,1 <BR>e,i,i+1 <BR>*enddo <BR>egen,22,14,295,307,1 <BR>r,2,543e-6 <BR>type,1 <BR>mat,1 <BR>real,2 <BR>*do,i,309,321,1 <BR>e,i,i+13 <BR>*enddo <BR>egen,20,13,581,593,1 <BR>*do,i,309,320,1 <BR>e,i,i+1 <BR>*enddo <BR>egen,21,13,841,852,1 <BR>r,3,433e-6 <BR>type,1 <BR>mat,1 <BR>real,3 <BR>nx=13 <BR>ny=21 <BR>nnup=308 <BR>*do,i,1,ny <BR>*do,j,1,nx <BR>k1=(i-1)*nx+j+nnup <BR>k2=k1-nnup+i-1 <BR>e,k1,k2 <BR>e,k1,k2+1 <BR>e,k1,k2+nx+1 <BR>e,k1,k2+nx+2 <BR>*enddo <BR>*enddo <BR>allsel,all <BR>acel,,0 <BR>et,4,mass21 <BR>r,6,850,850,850 <BR>type,4 <BR>real,6 <BR>*do,i,1,308,1 <BR>e,i <BR>*enddo <BR>local,11,1 <BR>nsel,s,loc,x,r <BR>local,11 <BR>!f,all,fy,-1000 !动力分析,恒载不于考虑 <BR>local,11,1 <BR>nsel,s,loc,z,0 <BR>nsel,a,loc,z,63 <BR>nsel,a,loc,y,b <BR>nsel,a,loc,y,180-b <BR>local,11 <BR>!f,all,fy,-500 <BR>/REP <BR>local,11,1 <BR>nsel,s,loc,x,r-h <BR>local,11 <BR>!f,all,fy,-500 <BR>local,11,1 <BR>nsel,s,loc,z,1.5 <BR>nsel,a,loc,z,61.5 <BR>nsel,a,loc,y,(90+12*b)/13 <BR>nsel,a,loc,y,(2250-12*b)/13 <BR>local,11, <BR>!f,all,fy,-250 <BR>allsel,all <BR>GPLOT !进行总绘图 <BR>/PBC,ALL, ,1 !显示边界条件及值 <BR>/rep <BR>/solu <BR>nsel,s,node,,2,13,1 <BR>nsel,a,node,,296,307,1 <BR>d,all,uy,0 <BR>nsel,all <BR>nsel,s,node,,1,295,14 <BR>local,11,1 <BR>nsel,a,loc,y,180-b <BR>d,all,uy <BR>allsel <BR>fini <BR>/prep7 <BR>local,11,0 <BR>ngen,2,581,1,295,14, <BR>ngen,2,876,14,308,14, <BR>et,2,combin14,,,2 !叠层橡胶!(ke bu laolui) <BR>r,4,0.41e6 <BR>type,2 <BR>mat,1 <BR>real,4 <BR>*do,i,1,295,14 <BR>e,i,i+581 <BR>*enddo <BR>*do,i,14,308,14 <BR>e,i,i+876 <BR>*enddo !(ke bu laok) <BR>et,3,combin39,1,,, <BR>r,5,0,0,0.00002,0.009e4,0.0008,0.0225e4 !夹芯橡胶 <BR>type,3 <BR>mat,1 <BR>real,5 <BR>*do,i,1,295,14 <BR>e,i,i+581 <BR>*enddo <BR>*do,i,14,308,14 <BR>e,i,i+876 <BR>*enddo <BR>/solu <BR>d,1,uy,,,14,13,uz, <BR>d,295,uy,,,308,13,uz, <BR>nsel,s,node,,582,876,14 <BR>nsel,a,node,,890,1184,14 <BR>d,all,all,0 <BR>allsel,all <BR>fini <BR></P> |