一、正弦曲线
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 save
2、直接用解析法算出切点
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