ansys螺旋线画法的几个例子
1. EXAMPLE_1CSYS,1 *AFUN,DEG /prep7 ET,1,45 *do,i,0,36,1 x=10 y=10*i z=i/2 k,i+1,x,y,z *enddo *do,j,1,36,1 L,j,j+1 *enddo *do,i,0,35,1 LWPL,-1,i+1,0 CYL4,0,0,1 *enddo *do,i,1,36,1 VDRAG, i, , , , , , i, , , , , *enddo vglue,all TYPE,1 vmesh,all过程图：2. EXAMPLE_2/PREP7 wpoff,-10 CYL4, , ,1 ET,1,MESH200 ET,2,SOLID45 KEYOPT,1,1,6 KEYOPT,1,2,0 LESIZE,ALL,10,1 MSHAPE,0,2D MSHKEY,1 AMESH,ALL CSYS,0 WPAVE,0,0,0 CSYS,0 !* CSYS,1 KBETW,3,1,0,RATI,0.5, ！Creates a keypoint between two existing keypoints*do,i,1,72 k,5+i,10,180-i*10,i lstr,5+i,4+i *enddo *do,i,1,72 VDRAG,1+(i-1)*5, , , , , ,4+i eplot *enddo3. EXAMPLE_3finish /clear /prep7*AFUN,RADpi=2*ASIN(1) !获得pi常数 *afun,degR1=6.3 !轴外半径rs=1.65 !丝半径R2=R1+rs !螺旋外半径h0=64 !螺距H=h0 !螺旋总高度 K,1,R2,0,0theta=90-Atan(64/(2*pi*R2) !获得旋转角wprota,0,-theta,0KWPAVE, 1 ！Moves the working plane origin to the average location of keypoints.rad1=rs !圆半径 rad2=rad1*0.6 !四边形半径 *get,knum,kp,num,max CYL4, , ,rad1 CSYS,1 !选择柱坐标系circle_p=10*get,knum,kp,num,maxK,knum+1,R2,0,0nn=NINT(h*circle_p/h0)*DO,i,1,nn-1,1xx=R2yy=i*360/circle_pzz=i*h0/circle_pK,knum+i+1,xx,yy,zz*ENDDOK,knum+nn+1,xx,h/h0*360,hLSEL,U, , ,ALL*DO,i,1,nn,1L,knum+i,knum+i+1*ENDDOLCOMB,ALL,0*GET,lnum,LINE,0,NUM,MAX !获得最大线条号ET,2,SOLID45TYPE, 2LESIZE,lnum, , ,40, , , , ,1VDRAG,ALL,lnum4. EXAMPLE_4钢丝软轴的模型/prep7 *AFUN,DEG CSYS,1 *do,i,1,500 k,i,2,10*(i-1),0.01*(i-1) *enddo *do,i,1,499 l,i,i+1 *enddo FLST,2,499,4,ORDE,2 FITEM,2,1 FITEM,2,-499 LCOMB,P51X,0 CSYS,0 KWPAVE, 1 wprot,0,90 PCIRC,0.18,0 ,0,360, VDRAG, 1, , , , , 1 FINISH5. EXAMPLE_5近来算了一个螺旋形的钢桥,以下是一段用APDL建立螺旋线的代码,主要使用"BSPLIN"命令来绘制空间曲线,然后再进行网格划分。适当修改一下可以用做各种螺旋形模型的建立。/PREP7PI=3.1415926R=18.100 !螺旋线半径H=24.836 !螺旋线单圈高N=6 !曲线关键点数量ANG1=45.18 !曲线起始角度ANG2=135.71 !曲线终止角度HH=H*(ANG2-ANG1)/360 !曲线段高度DH=HH/(N-1) !每个关键点的高度增量DA=(ANG2-ANG1)/(N-1) !每个关键点的角度增量*DO,I,0,N-1XX=R*COS(ANG1+I*DA)*PI/180)YY=R*SIN(ANG1+I*DA)*PI/180)ZZ=I*DHK,I+1,XX,YY,ZZ*ENDDOBSPLIN,1,2,3,4,5,6*GET,LENG, LINE, 1, LENG !得到螺旋线的长度6. EXAMPLE_6Ansys生成螺旋线模型/prep7/title,helix heater/units,user,1e-3,273Rh=2.5 !螺“线”半径H=20 !总长dis=0.8 !螺距rs=0.5 !丝径*afun,degtheta=90-atan(dis/Rh)k,1,0,0,0k,2,0,0,Hk,3,Rh,0,0k,4,Rh,0,Hl,3,4 !可为曲线,但要修改参考柱面的编号arotat,1,1,2 !生成面上两点间最短线的参考面wprota,0,-theta,0 !旋转工作平面xc=Rhyc=0cyl4,xc,yc,rs !丝截面,用于拉伸成线wpcsys,0KPstr=500k,KPstr,Rh,0,0num=H/dis*dim,C,array,num,3,Number_i,Coordinate_icsys,1*do,i,1,num C(i,1)=Rh C(i,2)=ky(i+KPstr-1)+90 C(i,3)=kz(i+KPstr-1)+dis k,i+KPstr,C(i,1),C(i,2),C(i,3) !生成下一个关键点 NSarea=mod(i-1),4)+1 *get,Narea,area,0,num,max larea,i+KPstr-1,i+KPstr,NSarea !由关键点生成拉伸路径 *get,Nline,line,0,num,max vdrag,Narea,Nline !拉伸成体*enddocsys,0adele,1,4vadd,allnumcmp,all
