资源描述
!EX8.26 钢筋混凝土简支梁数值分析!分离式模型,关闭压碎,keyopt(l)=0,keyopt(7)=l !力加载,位移收敛准则,误差1.5%,1/4 模型分析 !finish /clear /config,nres,2000/prep7!l. 定义单元与材料性质et,l,solid65,l!Kl=0,k7=let,2,link8 mp,ex,l,l3585 mp,prxy,l,0.2fc=l4.3 ft=l.43 tb,concr,l tbdata,0.5,0.95,ft,-l tb,miso,l,ll tbpt,0.0002,fc*0.l9 tbpt,0.0004,fc*0.36 tbpt,0.0006,fc*0.5l tbpt,0.0008,fc*0.64 tbpt,0.00l,fc*0.75 tbpt,0.00l2,fc*0.84 tbpt,0.00l4,fc*0.9l tbpt,0.00l6,fc*0.96 tbpt,0.00l8,fc*0.99 tbpt,0.002,fc tbpt,0.0033,fc*0.85 mp,ex,2,2.0e5 mp,prxy,2,0.3tb,biso,2 tbdata,300,0 pi=acos(-l) r,l,0.25*pi*22*22 r,2,0.25*pi*22*22/2 r,3,0.25*pi*l0*l0 r,4,0.25*pi*l0*l0/2 !2.创建几何模型 blc4,l50/2,300,2000/2*do,i,l,9 wpoff,l00 vsbw,all*enddo wpcsys,-1 wpoff,50 vsbw,all wpcsys,-1 wprota,-90 wpoff,30 vsbw,all wpoff,240 vsbw,all wpcsys,-1 wpoff,30 wprota,90 vsbw,all wpcsys,-1 !3.划分钢筋网格 elemsiz=50 lsel,s,loc,x,30 lsel,r,loc,y,30 cm,zj,line latt,2,1,2 lesize,all,elemsiz lsel,s,loc,x,75 lsel,r,loc,y,30 cm,zjb,line latt,2,2,2 lesize,all,elemsiz lsel,s,loc,x,30 lsel,r,loc,y,270 cm,jlj,line latt,2,3,2 lesize,all,elemsiz lsel,s,tan1,z lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 lsel,u,loc,z,50 cm,gj,line latt,2,3,2 lesize,all,elemsiz lsel,s,loc,z,0 lsel,r,loc,y,30,270 lsel,r,loc,x,30,70 cm,gjb,line latt,2,4,2lesize,all,elemsiz lsel,all cmsel,s,zj cmsel,a,zjb cmsel,a,jlj cmsel,a,gj cmsel,a,gjb cm,gj,line lmesh,all lsel,all !4.划分混凝土网格 vatt,1,1 mshkey,1 esize,elemsiz vmesh,all allsel,all !5.施加荷载和约束 lsel,s,loc,y,0 lsel,r,loc,z,900 dl,all,uy asel,s,loc,z,0 da,all,symm asel,s,loc,x,75 da,all,symm !p0=180000 q0=p0/150/100 asel,s,loc,z,0,50 asel,r,loc,y,300 sfa,all,1,pres,q0 allsel,all !6.求解控制设置 /solu antype,0 nsubst,80 outres,all,all autos,on neqit,50 cnvtol,u,0.015 solve!7.进入 POST1 查看结果 /post1set,lastpldisp,1esel,s,type,2 etable,saxl,ls,1 plls,saxl,saxlesel,s,type,1 /device,vector,on plcrackeall !8.进入时程后处理 /post26 nsol,2,205,u,y prod,3,1,p0/1000 prod,4,2,-1 xvar,4 plvar,3!EX2.6D 60kg/m 钢轨实体建模FINISH/CLEAR/PREP7!1.定义参数H1=30.5H2=48.5H3=48.5B1=45.75B2=29.25B3=16.5B4=36.5B5=10XD1=1/3XD2=1/9XD3=1/3XD4=1/20R0=400R1=2R2=4R3=40R4=20R5=25R6=8R7=13R8=80R9=300R10=2FAI=43H=H1+2*H2+H3 !2.创建轨底主要组成部分的线K,1K,2,B1+B2K,3,B1+B2,H1-B1*XD1-B2*XD2K,4,B1,H1-B1*XD1K,5,0,H1*DO,I,1,4L,I,I+1*ENDDO!3.计算确定轨头下三个圆弧的位置CTA=ATAN(XD3)BTA=(90-FAI)*ACOS(-1)/180-CTAA=1/TAN(CTA)C1=B3/2+R0-R6*SIN(CTA)-(R5-R6)*COS(BTA) C2=R5*SIN(BTA)+R6*COS(CTA)-R6*SIN(BTA)-H2C=(C1-C2/TAN(CTA)/(R0-R5)C3=A*C-SQRT(A*A-C*C+1)C3=C3/(A*A+1)REFA=ASIN(C3)C1=R6*SIN(CTA)+(R5-R6)*COS(BTA)+(R0-R5)*COS(REFA) !4.创建轨腰和轨头结构线ROUX=B3/2+R0-C1ROUY=ROUX*TAN(CTA)ANGI=ASIN(H1+H2)/R0)K,6,B3/2+R0-R0*COS(ANGI) K,7,B3/2+R0-R0*COS(REFA),H1+H2+R0*SIN(REFA)LARC,6,7,2,R0DX1=R6*(COS(BTA)-SIN(CTA)DY1=R6*(COS(CTA)-SIN(BTA)K,8,ROUX-DX1,H1+2*H2+ROUY-DY1K,9,ROUX,H1+2*H2+ROUYLARC,7,8,2,R5LARC,8,9,2,R6K,10,B4,H1+2*H2+B4*XD3K,11,B4-(H3-B4*XD3)*XD4,HL,9,10L,10,11ANGI=ASIN(B5/R9)K,12,0,HK,13,B5,H-(R9-R9*COS(ANGI)LARC,12,13,1,R9DX1=B4-B5+R8*SIN(ANGI)ANGI1=ASIN(DX1/R8)Y1=R8*COS(ANGI)-R8*COS(ANGI1)K,14,B4,KY(13)-Y1LARC,13,14,1,R8!5.倒角或弧线连接LFILLT,1,2,R1LFILLT,2,3,R2LFILLT,3,4,R3LPTN,4,5LFILLT,17,18,R4LFILLT,8,9,R10LPTN,9,11LFILLT,20,22,R7LDELE,15,16,1,1LDELE,19,21,2,1LSYMM,X,ALLNUMMRG,ALLNUMCMP,ALLAL,ALLASUMVOFFST,1,300!EX2.9B 斜向圆台建模 finish/clear/prep7r0=3 r1=2 h=4 cyl4,r0 wpoff,r1,h cyl4,r1 *do,i,1,4 askin,i,i+4*enddova,all!斜向延伸创建圆台 vgen,1,1,3*r0,1 wpcsyscyl4,r0vext,7,r1,h,r1/r0,r1/r0!EX2.9B 已知函数方程时的建模!玫瑰线 1=!极坐标方程为P =asin(4/3)finish/clear/prep7a=20csys,1*do,i,0,360*3fei=i*acos(-1)/180r0=a*sin(fei*4/3)k,r0,i*enddo*do,i,1,360*3l, i,i+1*enddo!玫瑰线 2=!极坐标方程为P =asin(4)finish/clear/prep7a=20csys,1*do,i,0,360fei=i*acos(-1)/180r0=a*sin(fei*4)k,r0,i*enddo*do,i,1,360l, i,i+1*enddo!长辐圆内旋轮线=!参数方程 x=(a-b)cost+ 入 cos(a-b)/bt!参数方程 y=(a-b)sint-入 sin(a-b)/btfinish/clear/prep7a=50b=5lmda=15*afun,deg*do,i,0,360x1=(a-b)*cos(i)+lmda*cos(a-b)/b*i) y1=(a-b)*sin(i)-lmda*sin(a-b)/b*i)k,x1,y1*enddo*do,i,1,360l, i,i+1*enddo!椭圆抛物面=!方程为 z=x*x/a/a+y*y/b/b! 使用蒙皮方法finish/clear/prep7a=20b=10n=20*do,i,1,niz=i/10lsel,nonewpoff,izcyl4,1adele,alllsscale,all,a*sqrt(iz),b*sqrt(iz),1*enddoallsel*do,i,1,4*(n-1),4askin,i,i+4askin,i+1,i+5askin,i+2,i+6askin,i+3,i+7*enddo!EX2.9C 旋转图素建模 !类花瓶建模finish/clear /prep7 K,1 K,2,100 K,3,300,400 K,4,150,550 k,5,130,880 K,6,300,1150 k,7,400,1100 k,8,500,1150 bsplin,all k,1001,1000 arotat,all,1,1001,360,5 !类轮建模finish/clear /prep7 !先创建类轮结构的一个断面,然后旋转之 blc4,4,6cyl4,2,6,1.5 asba,1,2 k,100,-3 k,101,-3,10 vrotat,all,100,101 !EX3.1 自定义箱形截面finish/clear/prep7k,1k,2,2k,3,2,2.2k,4,3,2.3k,5,3,2.5k,6,0,2.5a,1,2,3,4,5,6blc4,0.2,1.7,2asba,1,2wpoff,1.7wprota,90asbw,all wpoff,0.3 asbw,all wpoff,0.2 wprota,90 asbw,all wpoff,-2 asbw,all arsym,x,all wpcsys aglue,all et,1,plane82 lesize,all,1 amesh,all secwrite,myboxfinish /clear /prep7 et,1,beam189 mp,ex,1,3.0e10 mp,prxy,1,0.167 sectype,1,beam,mesh secread,mybox,mesh k,1k,2,10 k,100,5,5l,1,2 lesize,all,20 latt,1,1,100,1 lmesh,all /eshape,1 Eplot!EX3.2 自定义多种材料截面 finish /clear/prep7Ro=1.5Ri=1.0 csys,1 cyl4,ricyl4,ro aptn,all wprota,90 asbwa,all wprota,90 asbw,all wpcsys et,1,plane82 mymat1=4 mymat2=7 mp,ex,mymat1,1.0 mp,ex,mymat2,2.0 asel,s,loc,x,0,ri aatt,mymat1,1 asel,s,loc,x,ri,ro aatt,mymat2,1 allsel esize,0.25 mshape,0,2d mshkey,1 amesh,all secwrite,mycsolid,sectfinish /clear /prep7 et,1,beam189 mym1=4 mym2=7 mp,ex,mym1,3.0e10 mp,prxy,mym1,0.167 mp,ex,mym2,2.1e11 mp,prxy,mym2,0.3 sectype,1,beam,mesh secread,mycsolid,sect,mesh k,1k,2,10l,1,2 lesize,all,20 latt,1,1 lmesh,all /eshape,1 /pnum,mat,1 eplot!EX3.3A 单个方位关键点示例 FINISH/CLEAR /PREP7ET,1,BEAM189 MP,EX,1,2.1E5 MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8 K,1K,2,1000L,1,2K,100,500,500 LATT,1,1,100,1 LGEN,4,1,500LESIZE,ALL,10LMESH,ALL /ESHAPE,1 EPLOT!EX3.3B同时使用KB和KE方位关键点示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,L0L, 1,2K,100,DLK,200,DXC,-DLK,300,2*DXC,DLK,301,2*DXC+DLK,400K,500,8*DXCLGEN,5,1,DXCLSEL,S,1LATT,1,1,100,1LSEL,S,2LATT,1,1,200,1LSEL,S,3LATT,1,1,300,301,1LSEL,S,4LATT,1,1,400,1LSEL,S,5LATT,1,1,500,1LSEL,ALLLESIZE,ALL,50LMESH,ALL/ESHAPE,1EPLOT!EX3.3C 双方位关键点翻转示例FINISH/CLEAR/PREP7ET,1,BEAM189MP,EX,1,2.1E5MP,PRXY,1,0.3SECTYPE,1,BEAM,ISECDATA,100,40,160,10,10,8L0=1000DL=500DXC=400K,1K,2,L0K,3,DXCK,4,DXC,L0/2K,5,DXC,L0L,1,2L,3,4L,4,5K,100,DLk,101,-DLK,200,DXC,DLk,201,2*DXCK,202,DXC,-DLLSEL,S,1LATT,1,1,100,101,1 lesize,all,100LSEL,S,2LATT,1,1,200,201,1 lesize,all,50LSEL,S,3LATT,1,1,201,202,1 lesize,all,50LSEL,ALLLMESH,ALL/ESHAPE,1EPLOT!EX3.4A 两种单元形状和两种网格划分比较finish/clear/prep7et,1,plane82k,1k,2,8k,3,7,6k,4,1,6a,1,2,3,4esize,1mshape,0mshkey,0!mshape,0!mshkey,1!mshape,1!mshkey,0!mshape,1!mshkey,1amesh,all!EX3.4B 中间节点位置控制网格划分比较finish/clear/prep7et,1,plane82 cyl4,4,8,60 lesize,all,2 mshape,0 mshkey,1 mshmid,0 !mshmid,1 !mshmid,2 amesh,all!EX3.5 线上单元尺寸设置示例 !下边密上边稀 finish /clear/prep7 et,1,plane82 blc4,10,10 lsel,s,tan1,y lesize,all,10 lsel,s,loc,x,0 lesize,all,9,1/8 lsel,s,loc,x,10 lesize,all,9,8 lsel,all mshape,0 mshkey,1 amesh,all !中间密外边稀 finish /clear /prep7 et,1,plane82 blc4,10,10 lsel,s,tan1,y lesize,all,10,-1/5 lsel,s,tan1,x lesize,all,9,-1/8 lsel,all mshape,0 mshkey,1 amesh,all!EX3.6 合并线和连接线以进行映射网格划分 finish/clear /prep7 et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,4l,1,2l,2,3l,3,4 larc,4,5,3,10l,5,6l,6,1al,all esize,3 mshape,0 mshkey,1 lccat,1,2 lcomb,4,5 amesh,all!EX3.7 合并线或连接线的网格划分设置finish/clear/prep7et,1,plane82k,1,5k,2,10k,3,11,6k,4,6,15k,5,-1,8k,6,4l,1,2l,2,3l,3,4larc,4,5,3,10l,5,6l,6,1 al,all esize,10 lesize,6,8 lesize,1,4 lesize,2,3 lesize,4,5 lesize,5,2 mshape,0 mshkey,1 lcomb,1,2 !lccat,1,2!如采用lcatt,l,2则采用7个划分数 lcomb,4,5 amesh,all!EX3.7B 简化映射网格划分 finish/clear /prep7 et,l,plane82 k,l,5 k,2,l0 k,3,ll,6 k,4,6,l5k,5,-l,8 k,6,4l, l,21.3.4 larc,4,5,3,l0l,5,6l,6,lal,all esize,6 mshape,0 amap,l,2,5,3,4!EX3.8A 过渡四边形映射网格finish/clear/prep7et,1,plane42k,1 k,2,10,-1 k,3,8,6 k,4,1,3 a,1,2,3,4 lesize,1,8 lesize,3,3 lesize,4,7 lesize,2,2 mshape,0,2d mshkey,1 amesh,all !EX3.8B 过渡四边形映射网格 finish /clear /prep7 et,1,plane42k,1 k,2,10,-1 k,3,8,6 k,4,1,3 a,1,2,3,4 lesize,1,11 lesize,3,3 lesize,4,2 lesize,2,2 mshape,0,2d mshkey,1 amesh,all!EX3.9A 过渡六面体映射网格划分 finish/clear/prep7et,1,95 blc4,8,8,8 lesize,all,3 lesize,5,15 allsel mshape,0,3dmshkey,1 vmesh,all !EX3.9B 过渡六面体映射网格划分 finish /clear/prep7 et,1,95 blc4,8,8,8 lesize,all,4 lesize,7,12 mshape,0,3d mshkey,1 vmesh,all !EX3.9C 过渡六面体映射网格划分 finish /clear/prep7 et,1,95 blc4,8,8,8 lsel,s,1,8,7 lesize,all,6 lsel,s,2,7,5 lesize,all,3 lsel,s,4,5 lesize,all,7 lsel,s,3,6,3 lesize,all,2 lsel,s,9,12,1 lesize,all,5 allsel mshape,0,3d mshkey,1 vmesh,all!EX3.9D 过渡六面体映射网格划分 finish /clear/prep7 et,1,95 blc4,8,8,8 lesize,2,2 lesize,4,2 lesize,5,2 lesize,1,3 lesize,3,3lesize,8,3 lesize,7,6 lesize,6,7 lesize,11,7 lesize,9,3 lesize,10,3 lesize,12,3 allsel mshape,0,3d mshkey,1 vmesh,all!EX3.10 带两孔的长方体的扫掠网格划分 finish /clear/prep7 a=10 r=2 et,1,mesh200,6 et,2,solid45 blc4,2*a,a,a cyl4,a/2,a/2,r,a wprota,90 cyl4,1.5*a,a/2,r,-a vsbv,1,2 vsbv,4,3 wprota,90 wpoff,a/2 vsbw,all wpoff,a/2 vsbw,all wpoff,a/2 vsbw,all wpcsys wpoff,a/2,a/2 vsbw,all wprota,90 vsbw,all wpcsys esize,1 amap,105,15,16,26,63 amap,107,16,13,60,26amap,108,13,14,28,60 amap,103,15,14,28,63 lesize,94,4 lesize,79,5 vsel,s,loc,x,0,a vsweep,all asel,s,loc,y,a asel,r,loc,x,a,2*a lccat,2,45 lccat,57,71 lccat,68,78 lccat,65,72 mshape,0,2d mshkey,1 amesh,all vsel,s,loc,x,a,2*a vsweep,all allsel /view,1,1,2,3 Eplot!EX3.11 拉伸类命令生成体单元网格 finish/clear /prep7 et,1,82 et,2,95 blc4,4,4 blc4,6,4,4 blc4,12,4,4 blc4,18,4,4 esize,1 amesh,all esize,8 vrotat,1,1,4,90 vext,2,10,0.5,0.5 voffst,3,10 vdrag,4,35!EX3.12 圆的网格划分finish/clear/prep7 et,1,plane82 r0=10 cyl4,r0 cyl4,3*r0,r0 wprota,90 asbw,all wprota,90 asbw,all wpoff,3*r0 asbw,all wpcsys,-1 asel,s,loc,x,-r0,r0 lsla,s lesize,all,8 mshape,0,2d mshkey,1 amesh,all allsel asel,s,loc,x,2*r0,4*r0 lesize,all,8 lsel,r,length,r0 lesize,all,8,0.1,1 amesh,all allsel!EX3.13 圆环的网格划分 finish/clear /prep7 et,1,plane82 r0=10 cyl4,r0/3,r0,90 cyl4,2*r0,r0/10,r0,90 asel,s,loc,x,-r0,r0 lsla,s lesize,all,8 lsel,r,length,r0*2/3 lesize,all,3,1 mshape,0,2d mshkey,1 amesh,all asel,s,loc,x,2*r0,4*r0lsla,s lesize,5,12 lesize,7,6 lsel,s,length,r0*9/10 lesize,all,7 amesh,all!EX3.14 圆柱面和圆柱体的网格划分 !圆柱面 finish /clear/prep7 r0=10 h0=50 et,1,shell63 cyl4,r0 adele,1 cm,l1cm,line k,50 k,51,h0 l,50,51 adrag,l1cm,5 lsel,s,loc,z,0 lesize,all,6 lsel,s,length,h0 lesize,all,8 mshape,0,2d mshkey,1 amesh,all !圆柱体 finish /clear /prep7 r0=10 h0=50 et,1,solid95 cyl4,r0,h0 wprota,90 vsbw,all wprota,90 vsbw,all mshape,0,3d mshkey,1 lsel,s,loc,z,0lesize,all,6 lsel,s,length,h0 lesize,all,8 vmesh,all!EX3.15 圆锥的网格划分 !扫掠网格划分 finish /clear/prep7 cone,10,15,90 et,1,200,7 et,2,95 lsel,s,5,6 lesize,all,12,0.5 lsel,all lesize,3,6 mshape,0,2d mshkey,1 amesh,3 vsweep,1,3,4 vsymm,x,all vsymm,y,all vglue,all !六面体映射网格划分 finish /clear/prep7 cone,10,15,90 et,2,95 lesize,all,8 mshape,0,3d mshkey,1 vmesh,all vsymm,x,all vsymm,y,all vglue,all!EX3.16 正多边形面的通用网格划分 finish /clear/prep7 ns=11et,1,plane82 rpr4,ns,10 kp0=100+nsk,kp0 *if,mod(ns,2),eq,0,then *do,i,1,ns/2 l,kp0,2*i-1*enddo*else*do,i,1,nsl,kp0,i*enddo*endif lsel,s,ns+1,2*ns cm,l1cm,line allsel asbl,1,l1cm *if,mod(ns,2),eq,0,then cmsel,s,l1cm lesize,all,8 lsel,inve lesize,all,4*elselesize,all,8*endifallselmshape,0mshkey,1amesh,all!EX3.17 球体及球面网格划分 finish /clear/prep7 r0=10 sphere,r0,90 vsbw,all vdele,2,1 esize,2 et,1,solid95 mshape,0,3d mshkey,1vmesh,allvsymm,x,all vsymm,y,all vsymm,z,all vglue,all finish /clear /prep7 r0=10 sphere,r0,90 vsbw,all vdele,2,1 vdele,all asel,s,loc,x,0 asel,a,loc,y,0 asel,a,loc,z,0 adele,all,1 allsel esize,2 et,1,shell63 mshape,0,2d mshkey,1 amesh,all arsym,x,all arsym,y,all arsym,z,all nummrg,all!EX3.18 孔板网格划分 finish/clear /prep7 a0=300 b0=800 r0=15 blc4,a0,b0 cyl4,a0/4,b0/8,r0 agen,2,2,a0/2 agen,2,2,3,1,b0/8 agen,2,2,5,1,b0*5/8 asel,s,2,9,1 cm,a2cm,area allselasba,1,a2cm wprota,-90 *do,i,1,5 wpoff,b0/16 asbw,all *enddo wpoff,b0*5/16 *do,i,1,5 wpoff,b0/16 asbw,all *enddo wprota,90 *do,i,1,3 wpoff,a0/4 asbw,all *enddo wpcsys,-1 numcmp,alllsel,s,radius,r0 lesize,all,8 lsel,inve lesize,all,4 lsel,all et,1,82 mshape,0,2d mshkey,1 asel,u,loc,y,b0/16,b0*5/16 asel,u,loc,y,b0*11/16,b0*15/16 lsla,s lsel,r,tan1,x lesize,all,50,1 amesh,all allselamap,21,30,31,54,62 amap,22,31,32,54,64 amap,32,29,30,62,76 amap,33,29,32,64,76amap,42,34,35,76,82 amap,43,35,36,76,84 amap,41,33,34,53,82 amap,44,33,36,53,84amap,19,22,23,49,64amap,20,23,24,49,60amap,9,21,22,64,74amap,31,21,24,60,74amap,38,26,27,74,84amap,39,27,28,74,80amap,37,25,26,50,84amap,40,25,28,50,80amap,17,14,15,43,59amap,18,15,16,43,63amap,7,13,14,59,73amap,8,13,16,63,73amap,34,18,19,73,79amap,30,17,18,79,44amap,35,19,20,73,83amap,36,17,20,44,83amap,15,6,7,39,63amap,16,7,8,39,58amap,4,6,5,63,71amap,5,5,8,58,71amap,27,10,11,71,83amap,28,11,12,71,78amap,26,9,10,83,40amap,29,9,12,40,78!EX3.19 角支架的网格划分 finish/clear/prep7blc4,150,50 blc4,100,50,-50 cyl4,25,25 cyl4,125,-50,25 aadd,all numcmp,all cyl4,25,10cyl4,125,-50,10 asel,s,2,3 cm,a1cm,area asel,all asba,1,a1cm lcomb,1,6 lfillt,1,2,20 asbl,4,6 adele,1,1 lfillt,3,4,20 al,18,19,20 aadd,all numcmp,all !wprota,90 wpoff,-25 asbw,all wpoff,75 asbw,all wprota,90 asbw,all wpoff,125 asbw,all wpcsys,-1 wpoff,25 wprota,90 asbw,all kwpave,18 asbw,all kwpave,3 wprota,90 asbw,all kwpave,21 wprota,-45 asbw,8 wpcsys,-1 !et,1,plane82 mshape,0,2d mshkey,1 esize,6 lesize,33,6 lesize,37,6 lesize,42,6 amesh,3,5,2 amap,6,9,10,4,23 amap,7,9,12,1,23 amesh,11,13,2 amesh,1,9,8 amesh,2,4,2 amap,16,3,18,26,28 amap,12,14,15,5,28 amap,15,13,14,6,28 lcomb,19,22 lcomb,27,46 amesh,10,14,4 != !EX3.20 大板小孔的网格划分 A finish/clear /prep7 a0=100 blc4,a0,a0 cyl4,a0/100 asba,1,2 csys,1 k,50,a0/20 k,51,a0/20,90l,50,51 asbl,all,1 wprota,90,45 asbw,all et,1,plane82 mshape,0,2d mshkey,1 lsel,s,loc,x,0,a0/20 lesize,all,8 lsel,all lesize,7,12,0.1 lesize,8,12,0.1 lesize,9,12,10 amesh,all csys,0 arsym,x,all arsym,y,allnummrg,all != !EX3.20 大板小孔的网格划分 B finish/clear /prep7 a0=100 blc4,a0,a0 cyl4,a0/100 asba,1,2 csys,1 k,50,a0/20 k,51,a0/20,90l,50,51 asbl,all,1 et,1,plane82 mshape,0,2d mshkey,1 lesize,5,8 lesize,1,10 lesize,4,8,5 lesize,6,8,5 amesh,1 lesize,7,20,0.1 lesize,8,20,0.1 amap,2,50,51,2,4 csys,0 arsym,x,all arsym,y,all nummrg,all!EX3.21 Z形面网格划分及局部细分 finish /clear/prep7 blc4,15,10 blc4,10,6,14,12 aadd,all wprota,-90 wpoff,6 asbw,all wpoff,4asbw,allwprota,90 wpoff,10 asbw,all wpoff,5 asbw,all wpcsys,-1 et,1,82 esize,2 mshape,0,2d mshkey,1 amesh,all krefine,9,10,1,1,off !不光滑处理(节点位置不变) aclear,all amesh,allkrefine,9,10,1,1,smooth !光滑处理,节点位置可改变,但既有单元不删除 aclear,allamesh,allkrefine,9,10,1,1 !光滑处理,可删除既有单元重细分,节点位置改变!EX3.22 分布小孔环形面的网格划分 finish/clear /prep7 r1=15 r2=25 r3=3 n=8 cyl4,r1,r2,180/n cyl4,0.5*(r1+r2),r3 asba,1,2kl,2 larc,6,8,4,0.5*(r1+r2) asbl,all,4 lesize,4,4 lesize,10,4 lesize,6,8 lesize,3,6 lesize,8,4 lesize,5,8 lesize,9,4lesize,1,6 lesize,7,4 et,1,plane82 mshape,0,2d amap,2,6,7,4,8 amap,1,5,6,8,1 !lrefine,5,6,1,1,smooth arsym,y,all csys,1 agen,n,all,360/n nummrg,all!EX3.23 同心圆轴体的网格划分 finish /clear/prep7 cylind,50,0,100,150 cylind,40,0,200,250 cylind,30,0,100,250 cylind,20,0,50,250 cylind,10,0,0,250 vptn,all wprota,90 vsbw,all wprota,90 vsbw,all et,1,95 csys,1 lsel,s,length,50 lesize,all,5 lsel,inve lsel,u,radius,10,50 lesize,all,4 n1=6 n2=6 n3=8 n4=10 n5=12 lsel,r,loc,x,0,10 lesize,all,n1,1 lsel,s,radius,10 lesize,all,n1 lsel,s,radius,20lesize,all,n2 lsel,s,radius,30 lesize,all,n3 lsel,s,radius,40 lesize,all,n4 lsel,s,radius,50 lesize,all,n5allsel mshape,0,3d mshkey,1 vmesh,all!EX3.24A 圆柱与长方体组合模型的网格划分 !采用六面体映射网格划分 finish /clear/prep7 a0=30 h1=15 h2=25 r0=7 blc4,a0/2,a0/2,h1 cyl4,r0,90,h1+h2 vptn,all accat,4,6 et,1,95 esize,3 mshape,0,3d mshkey,1 vmesh,all asel,s,acca adele,all lsel,s,lcca ldele,all allsel vsymm,x,all vsymm,y,all nummrg,all !EX3.24B 圆柱与长方体组合模型的网格划分 !采用过渡六面体映射网格划分 finish /clear/prep7 a0=30 h1=15 h2=25 r0=7 blc4,a0/2,a0/2,h1 cyl4,r0,90,h1+h2 vptn,all lsel,s,radius,r0 lsel,a,length,r0 lesize,all,6 lsel,s,loc,z,0 lsel,a,loc,z,h1 lesize,all,7 asel,s,loc,x,a0/2 asel,a,loc,y,a0/2 accat,all allsel et,1,95 esize,3 mshape,0,3d mshkey,1 vmesh,all asel,s,acca adele,all lsel,s,lcca ldele,all allsel vsymm,x,all vsymm,y,all nummrg,all!EX3.25-任意圈数圆柱形螺旋弹簧的扫略网格划分FINISH /CLEAR /PREP7 !1.定义弹簧参数(同上) D=4C=8 DZ=C*D T=DZ/6 N=4.7 *IF,T,LT,D,THENT=D *ENDIF TKPD=90 !2.创建全部螺旋线 CSYS,1 TDEG=N*360 TDEG1=MOD(TDEG,TKPD) N0=(TDEG-TDEG1)/TKPD+1 *AFUN,DEG *DO,I,1,N0 CTA=(I-1)*TKPD Z=T/360*CTA K,I,DZ/2,CTA,Z *ENDDO *IF,TDEG1,LT,1.0E-2,THEN *ELSEN0=N0+1 CTA=CTA+TDEG1 Z=T/360*CTA K,N0,DZ/2,CTA,Z *ENDIF *DO,I,1,N0-1 L,I,I+1 *ENDDO CM,L1,LINE !3.在螺旋线端部创建簧丝截面 CSYS,0WPOFF,DZ/2 WPROTA,90 CYL4,D/2wprota,90 asbw,all wprota,90 asbw,all cm,a1cm,area VDRAG,a1cm,L1 et,1,mesh200,7 et,2,solid95 cmsel,s,a1cm lsla,s lesize,all,6 amesh,all allselesize,d/2 vsweep,all!EX3.26 具坑缺陷圆柱的网格划分 finish /clear/prep7 l0=50 r0=10 r1=5 r2=1 h0=5 et,1,solid45 cyl4,r0,l0 wpave,0,r0,l0/2 wprota,90 con4,r1,r2,h0 vsbv,1,2 wprota,90 vsbw,all wpoff,r1+3 vsbw,all wpoff,-2*(r1+3) vsbw,all wprota,90 vsbw,all wpoff,r1+1 vsbw,all wpoff,-2*(r1+1) vsbw,all wpcsys,-1 wprota,90 vsbw,all wpoff,-r1 vsbw,all wpoff,2*r1 vsbw,all wpcsys,-1 numcmp,all ksel,s,loc,z,l0/2 ksel,r,loc,x,0 ksel,u,loc,y,0*get,kp1,kp,0,num,min kp2=kpnext(kp1) ksel,s,loc,z,l0/2+r2*get,kp3,kp,0,num,min allsel k,1000,-r0,l0/2+5 l,kp3,1000 asel,none*get,l1,line,0,num,max arotat,l1,kp1,kp2 cm,a1cm,areaallselvsba,all,a1cm!vsel,s,loc,z,17,33!vsel,r,loc,x,-6,6!aslv,s!aplotaccat,128,61accat,160,138accat,224,202accat,183,35accat,31,134accat,151,164accat,215,231accat,82,195accat,130,5accat,162,145accat,229,209accat,192,89accat,12,123accat,117,155accat,109,219accat,14,178mshape,0,3dmshkey,1esize,1vmesh,all!EX3.27 长方体开圆柱槽的网格划分 finish /clear/prep7 blc5,20,10,30 wpoff,5,5 cyl4,8,10 wpoff,10 cyl4,4,12 vsel,s,2,3 cm,v1cm,volu allsel wpcsys,-1 vsbv,1,v1cm voffst,15,5 voffst,16,3 voffst,17,15 vptn,all wprota,90 vsbw,allwpcsys,-1 wpoff,5 vsbw,all wpoff,10 vsbw,all wpoff,12 vsbw,allet,1,mesh200,7 esize,2 mshape,0,2d mshkey,1 asel,s,loc,z,30 lccat,6,74 lccat,8,68 amesh,all asel,s,loc,z,0 lccat,3,73 lccat,1,67 amesh,all allsel et,2,solid95 mshape,0,3d mshkey,1 vsweep,all !ex4.
展开阅读全文