一、正弦曲线 finish /clear /FILNAME,sin,1 /prep7 k,100 CIRCLE,100,1,90 !Circle,pcent,rad,paxis,pzero,arc,nseg 产生圆弧线。 !该圆弧线为圆的一部分，依参数状况而定，与目前所在的坐标系统无关，点的号码和圆弧的线段号码会自动产生。 !Pcent为圆弧中心坐标点的号码； !RAD :圆的半径，若此值不输，则半径的定义为PCENT到PZERO的距离 !paxis 定义圆心轴正方向上任意点的号码； !Pzero定义圆弧线起点轴上的任意点的号码，此点不一定在圆上； !ARC :弧长（以角度表示），若输入为正值，则由开始轴产生一段弧长，若没输和，产生一个整圆。 !NSEG :圆弧欲划分的段数，此处段数为线条的数目，非有限元网格化时的数目。默认为4。 csys,1 KFILL,1,2,4,3,1 !KFILL,NP1,NP2,NFILL,NSTRT,NINC,SPACE 点的填充命令 !是自动将两点NP1,NP2间，在现有的坐标系下填充许多点，两点间填充点的个数（NFILL）及分布状态视其参数（NSTRT,NINC,SPACE）而定， !系统设定为均分填充。如语句 FILL,1,5，则平均填充3个点在1 和5 之间。 csys,0 k,7,1+3.1415926/2 /pnum,kp,on /pnum,line,on KFILL,1,7,4,8,1 KGEN,2,7,11,1,1 !Kgen,itime,Np1,Np2,Ninc,Dx,Dy,Dz,kinc,noelem,imove !注意：MAT,REAL,TYPE 将一起拷贝，不是当前的MAT,REAL,TYPE !Itime：拷贝份数 !Np1,Np2,Ninc：所选关键点 !Dx,Dy,Dz：偏移坐标 !Kinc：每份之间关键点号增量 !noelem: “0” 如果附有节点及单元，则一起拷贝。 ! “1”不拷贝节点和单元 !imove： “0” 生成拷贝 ! “1”移动原关键点至新位置，并保持号码，此时（itime,kinc,noelem）被忽略 LSTR,8,13 LSTR,9,14 LSTR,10,15 LSTR,11,16 LANG,2,3,90 !Lang,nl1,p3,ang,phit,locat 产生一新的线段，此新的线段与已存在的线段nl1的夹角为ang,phit为新产生点的号码 LANG,3,4,90 LANG,4,5,90 LANG,5,6,90 BSPLIN,1,17,18,19,20,12 !Bsplin,p1,p2,p3,p4,p5,p6,xv1,yv1,zv1,xv6,yv6,zv6 通过6点拟合一条曲线，并定义两端点的斜率。 !SPLINE,1,17,18,19,20,12 ！spline,p1,p2,p3,p4,p5,p6,xv1,yv1,zv1,xv6,yv6,zv6 通过6点拟合曲线，在每点之间形成一新线段，并可以定义两端点的斜率。 lsel,all lseL,u,14 LDELE,all,1 allsel,all NUMCMP,line KWPAVE,12 csys,4 !csys,kcn 声明坐标系统，系统默认为卡式坐标（csys,0）。 !kcn = 0 笛卡尔坐标系 (x,y,z) ! 1 柱坐标系 (r,z) ! 2 球坐标系 (r,) ! 4 工作平面坐标系 (wx,wy,wz) ! 5 柱坐标系（以Y轴为轴心）(r,y,) ! n 已定义的局部坐标系 lsymm,x,1 NUMMRG, KP,LOW !Nummrg,label,toler, Gtoler,action,switch 合并相同位置的item !label: 要合并的项目 ! node: 节点， Elem,单元，kp: 关键点 Line：线，area：面，volu：体 ! mat: 材料，type: 单元类型，Real: 实常数 ! cp：耦合项，CE：约束项，CE: 约束方程，All：所有项 !toler: 公差 !Gtoler：实体公差 !Action: sele 仅选择不合并 ! 空 合并 !switch: 较低号还是较高号被保留（low, high） !注意：可以先选择一部分项目，再执行合并。如果多次发生合并命令，一定要先合并节点，再合并关键点。 !合并节点后，实体荷载不能转化到单元，此时可合并关键点解决问题。 LCOMB, ALL,0 !Lcomb,nl1,nl2,keep !将两条相邻的线合并为一条线，keep=0时原线段删除，keep=1时保留。 KWPAVE,1 LGEN,all,-1,1习题（自选）、过一点作圆的切线 1、以圆心与点的连线为直径作圆与初始圆产生交点（即切点）fini /cle /filname,tangent1,1 /title,tangent1 /prep7 r1=10 xb=30 yb=20 k,1 k,2,xb,yb circle,1,r1 cyl5,0,0,xb,yb ADELE,1,0 !Adele,na1,na2,ninc,kswp !kswp=0时只删除掉面积本身，1时低单元一并删除。!KFILL,1,2,1,100!CIRCLE,100,1 /pnum,line,on /pnum,kp,on LOVLAP,4,7 lstr,2,11 NUMCMP,line LDELE,5,7,1,1 ldele,10,11,1,1 save2、直接用解析法算出切点 fini /cle /filname,tangent2,1 /title,tangent2 /prep7 r1=10 xb=30 yb=20 r2=sqrt(xb*xb+yb*yb) theta=acos(r1/r2) theta1=asin(yb/r2) xb1=r1*cos(theta-theta1) yb1=-r1*sin(theta-theta1) k,1 k,2,xb,yb k,3,xb1,yb1 CIRCLE,1,r1 lstr,2,3 save二、相交圆柱体 fini /cle /filname,cylinder,1 /prep7 /pnum,volu,on CYLIND,15,30,0,80 wpoffs,0,50,50 WPROTA,0,-60 csys,4 CYLIND,6,12,0,-50 vovlap,all VDELE,3,5,2,1 BLOCK,-2,2,-9,-12,0,-8 WPSTYL,1 !设置工作平面坐标系为极坐标系（即圆柱坐标系） csys,4 VGEN,3,1,120 vsbv,6,1 vsbv,5,2 vsbv,1,3 wprota,0,0,90 vsbw,all !用工作平面切割体 VDELE,6,9,3 vdele,1 vadd,all fini习题2.1（a）、轴习题2.1（b）、底座 /prep7 CYLIND,15,34,0,8 CYL4,27,0,2.5,8 WPSTYL,1 csys,4 VGEN,6,2,60 wpstyl wpoffs,8 csys,4 CYLIND,15,22,0,50 BLOCK,-4,4,10,16.65,-8,50 WPoffs,0,0,40 wprota,0,0,90 CYLIND,3,0,22 vs