在仿真中我建立了两平行半径r=0.25m的单位长导体,两导体相距2m,通流计算载流导体的受力和变形。
仿真的时候我先试用solid117单元算出电磁力,再耦合到结构单元计算出变形;
但在瞬态电磁仿真时,试用solid17单元可以算出电磁力,但是这个力该如何耦合到结构
单元呢?我查看了ldread的相关解释,其如下:
When transferring static magnetic forces from a SOLID117 model, ensure that all nodes and elements of the bodies intended to receive the forces are selected prior to transfer. See the discussion on force computation in the Low-Frequency Electromagnetic Analysis Guide.
For a full harmonic magnetic analysis, FORC represents the time-averaged force (use in conjunction with KIMG = 2). Values are in the nodal coordinate system for the force loads (F).
对静态和谐波分析均有说明,但是都瞬态分析,如果做到电磁结构耦合却没有说明,这一部分的使用该如何进行呢?
我尝试着使用solid62单元来进行瞬态仿真;我仿真vm172做了一个例子做了一个静态仿真分析,但是发现这个的计算结果不太对,具体流文件如下:
/prep7
et,1,solid62
et,2,solid97
mp,rsvx,1,2.6e-8
mp,murx,1,1
mp,ex,1,1e10
mp,prxy,1,0.3
mp,murx,2,1
wpoffs,-1
csys,4
cylind,0.25,,0,1
wpcsys,-1,0
wpoffs,1
csys,4
cylind,0.25,,0,1
wpcsys,-1,0
cylidn,2,,0,1
allsel,all
vovlap,all
numcmp,all
esize,0.05
mshkey,0
mshape,0
vatt,1,1,1
vsweep,1
vsweep,2
esize,0.5
mshkey,0
mshape,0
vatt,2,2,2
vsweep,3
allsel,all
pi=2*asin(1)
vSEL,S,volu,,1
eslv,s
BFE,ALL,JS,,,,16000/pi !保证施加电流为1000A,方便后期验证
vSEL,S,volu,,2
eslv,s
BFE,ALL,JS,,,,16000/pi
allsel,all
asel,s,ext
da,all,az,0 ! 可能不是这样
allsel,all
asel,s,area,,2
asel,a,area,,6
nsla,s,1
d,all,ux,0
d,all,uy,0
d,all,uz,0
allsel,all
/solu
allsel,all
antype,static
eqslv,sprase
cnvtol,f,1e-3
solve
/post1
vsel,s,volu,,1
eslv,s
etable,jsz,smisc,3 ! 获得Z方向的电流密度
ssum
/post1
vsel,s,volu,,1
eslv,s
etable,,nmisc,4 ! 获得X方向电磁力
etable,,nmisc,5 ! 获得y方向电磁力
etable,,nmisc,6 ! 获得z方向电磁力
ssum
通过后处理的JZ与我输入的JZ不一样,而且也无法获得各方向电磁力,计算有计算结果,电流密度分布正确,但是磁感应强度值不太对劲(根据安培环路定律推得),而且这个变形量也不太正确(与同条件下使用solid117计算结果比较),这么多问题,他出现的原因是什么啊?
我也一直在找,但是没找到原因,还恳请各位大侠给点提示。 |