用以下程序仿真
/batch,list
/filename,cichang
/prep7
shpp,off
/title,cichang
/com,
/com,
/com,
x=0.6
y=0.0005
z=0.3
d1=1.8
d2=0.15
d=0.002 !簧丝直径
a=0.03
c=1 !旋绕比
n=5 !圈数
dz=c*d1 !螺旋线的直径
t=dz/50 !螺距
et,1,117,1
et,2,117
et,3,117
emunit,mks
mp,murx,1,1
mp,murx,2,1
mp,murx,3,1
mp,rsvx,3,1.22e-6
CSWPLA,11,1,1/12,1 !根据当前定义的工作平面定义局部坐标系
!csys,1 !设置当前坐标系为柱坐标系
k,1,dz/2,0,-t/2 !创建两个关键点
k,2,dz/2,180
L,1,2 !创建半圈螺旋线
CSYS,0 !设置直角坐标系
LSYMM,z,1 !利用对称性生成另外半圈螺旋线
LSYMM,Y,2,,,,,1 !合并关键点,并将此两条线定义为原件l1
nummrg,all !合并关键点,
cm,l1,line
kwpave,1 !移动工作平面并旋转
wprota,,90
recteg,0,d,0,a !创建直径为d的圆面(簧丝截面)
vdrag,1,,,,,,l1 !拖拉面创建体
vgen,n,all,,,,,t !复制体n次
nummrg,kp !合并关键点,并将工作平面归位
wpcsys
vglue,all
block,-x,x,-y,y,-z/2,z
block,-(x+0.001),x+0.001,-(y+0.001),y+0.001,-(z/2+0.001),z+0.001
block,-(x+0.005),x+0.005,-(y+0.005),y+0.005,-(z/2+0.005),z+0.005
block,-(x+0.01),x+0.01,-(y+0.01),y+0.01,-(z/2+0.01),z+0.01
block,-(d1/2+0.05),d1/2+0.05,-(d2/2+0.05),d2/2+0.05,-(z/2+0.01),z+0.01
block,-(d1/2+0.2),d1/2+0.2,-(d2/2+0.2),d2/2+0.2,-(z/2+0.01),z+0.01
ALLSEL,ALL
VOVLAP,ALL
NUMCMP,VOLU
/TRLCY,VOLU,1,12,16,1
ALLSEL,ALL
vsel,s,volu,,2
vatt,3,1,1
vsel,s,volu,,1
vsel,a,volu,,3,11,1
vatt,2,2,2
vsel,s,volu,,12,16,1
vatt,1,3,3
alls
vsel,s,volu,,2
esize,0.01
vmesh,all
vsel,s,volu,,1
vsel,a,volu,,3,11,1
esize,0.02
vmesh,all
esha,1
vsel,s,volu,,13
esize,0.02
vmesh,all
esha,1
vsel,s,volu,,14
esize,0.03
vmesh,all
esha,1
vsel,s,volu,,15
esize,0.05
vmesh,all
esha,1
vsel,s,volu,,16
esize,0.07
vmesh,all
vsel,s,volu,,12
esize,0.09
vmesh,all
alls
csys,0
nsel,s,ext
d,all,az,0
nsel,all
allsel
esel,s,type,,2
emode,all,esys,11
esel,s,real,,2
bfe,all,js,,,1.67e7
alls
/solu
antype,harm
harfrq,30000
solve
save
finish
问题1:总是出现一个错误提示“Element type 3 is UNKNOWN, which cannot be used with the VMES command. Meshing of volume 13 aborted.”,但是划分网格之前,通过List中的命令可以看到,单元类型1,2,3和材料类型分别赋给了相应的体;
问题2:如果把“et,3,117”去掉,就可以正常划分网格,但是求解过程中又出现了这样的错误:“Nonlinear BH data or linear MUR data must be input for magnetic analysis.”可是通过List中相应的命令可以看到,单元类型1,2和材料类型1,2,3分别赋给了相应的element;
各位大神有知道原因的吗?请大家不吝赐教,谢谢!
|