本帖最后由 i我飞翔 于 2016-12-14 00:11 编辑
大家好,本人水货加小鲜肉,初学ansys 这个高大上的软件,想建立一个关于管道用线圈磁化的三维局部1/4模型,但是不知道怎么加载边界条件,加载电流部分已经解决,但是边界条件设置缺遇到问题,用的是静态分析和97号单元,请大神和朋友们帮助,先谢谢大家了,以下是个人仿真的程序,希望大家能够指点一二,以求共同进步,
/clear
/filname,xxx
r1=0.5 !线圈内径
th=1.85 !线圈厚度
hc=5.1 !线圈高度
rd=9 !导体半径
gao=12 !导体高
d=0.5 !线圈距导体距离
turn=400 !线圈匝数
fill=0.95473 !线圈填充因子
/prep7
et,1,solid97,1 !导体
et,2,solid97,2 !线圈
et,3,solid97 !空气
pi=3.1415926
mp,murx,1,1
mp,rsvx,1,2.7e-8
mp,murx,2,1
mp,rsvx,2,1.7241e-8
mp,murx,3,1
area=th*hc*1e-6 !线圈截面积
a_coil=pi*(r1+th)*(r1+th)-pi*r1*r1
v_coil=(a_coil*hc*1e-9)/4 !线圈体积
R,1,area,turn,v_coil,0,1,0,
RMORE,4,fill
cyl4,0,0,0,0,rd,90,gao
wpoffs,0,0,gao+d
cyl4,0,0,r1,0,r1+th,90,hc !线圈
wpoffs,0,0,-12.5
SPHERE,25,0,0,90 !空气场
SPHERE,50,0,0,90
/PNUM,VOLU,1 !显示实体编号
/REP,FAST
NUMCMP,VOLU
/REP,FAST
CSWPLA,11,1,1,1
CSWPLA,12,2,1,1
csys,0
alls
vovlap,all
numcmp,all
vsel,s,,,1 !划分被测体
vatt,1,,1,0 !vatt-定义所选单元属性 mat,real,type,esys
esize,0.75 !mat(material number)材料号,real(real constant set number)实常数,
mshkey,1 !type(element type number)单元类型编号,esys-单元坐标类型
mshape,0,3d
vmesh,all
vsel,s,,,2
vatt,2,1,2,11
esize,0.75
mshkey,1
mshape,0,3d
vmesh,all
vsel,s,,,4 !划分近场空气
vatt,3,,3,12
esize,1.5
mshkey,0
mshape,1,3d
vmesh,all
vsel,s,,,3 !划分远场空气
vatt,3,,3,12
esize,4.5
mshkey,0
mshape,1,3d
vmesh,all
alls
VLSCAL,all, , ,0.001,0.001,0.001, ,0,1 !模型缩小,米制
!******************************************扫掠划分远场空气测试
!vsel,s,,,3 !划分远场空气
!vatt,3,,3,12
!vatt,4,,4,12
!esize,4.5
!esize,40 !infin111作为远场单元,只能划分一层单元,不能为多层单元
!vsel,s,,,3
!vsweep,3
!***********************************************************************开始加载求解
/solu
ALLSEL,ALL
CSYS,0
!!!!!!!!!!=========================AZ,AY,AZ=0 1.7,1.5
!Asel,s,ext !!!!!!!
!da,all,ax,0 !!!!!!!
!da,all,ay,0 !!!!!!!
!da,all,az,0 !!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!===================== 0.029,0.15
!ASEL,S,AREA,,13
!da,all,ax,0
!da,all,ay,0
!da,all,az,0
!ASEL,S,LOC,X,0
!DA,ALL,AZ,0
!ASEL,S,LOC,Y,
!DA,ALL,AZ,0
!!!!!!!!!===================== 1.09,1.44
!ASEL,S,AREA,,13
!da,all,ax,0
!da,all,ay,0
!da,all,az,0
!ASEL,S,LOC,X,0
!DA,ALL,AX,0
!ASEL,S,LOC,Y,
!DA,ALL,AY,0
!!!!!!!!!========================= 0.016,0.1
!ASEL,S,AREA,,13
!da,all,ax,0
!da,all,ay,0
!da,all,az,0
!lsel,s,loc,,x,0
!lsel,r,loc,,y,0
!nsll,s,
!d,all,az,0
!!!!!!!!!=====================
!ASEL,S,LOC,X,0
!da,all,aX,0
!ASEL,S,LOC,Y,0
!da,all,aY,0
esel,s,mat,,1
nsle,s
cp,1,volt,all
!vsel,s,,,1
esel,s,mat,,2
!eslv,s !选中依附于体积的单元
nsle,s !选中附着在单元上的节点
cp,2,curr,all !耦合线圈电流自由度
BFE,all,VLTG,1,24*1.414 !加电压,峰值
alls
antype,harm
harfrq,10000
solve
/post1
esel,s,mat,,2
SET,1,LAST,1,1, , ,
PRNSOL,CURR !电流虚部值,保存
SET,1,LAST,1,0, , ,
PRNSOL,CURR !电流实部值,保存
请大家指点,谢谢
|