请教关于117后处理洛仑兹力计算的问题

fanyangyang

我的命令流如下，要求一段导体上的受自身磁场的受力。导体周围有铁磁物质
/PREP7
ET,1,117            
MP,MURX,1,1       !空气相对磁导率
MP,MURX,2,200   !铁磁材料的相对磁导率
MP,MURX,4,1       !导体的相对磁导率
MP,RSVX,4,3.57E-4      !导体的相对电导率
MP,MURX,5,1         !导体的相对磁导率
MP,RSVX,5,1.8E-7    !导体的相对电导率
MP,MURX,6,200   
MP,RSVX,6,1.8E-7  


BLOCK,0.005,0.025,-0.1,0.1,-0.01,-0.03             !建立铁磁体模型
BLOCK,-0.025,-0.005,-0.1,0.1,-0.01,-0.03
BLOCK,-0.025,0.025,-0.1,0.1,0,-0.01

CYLIND,0.0006,0,0,-0.01,0,360                          !建立导体模型
CYLIND,0.0006,0,-0.01,-0.0113,0,360
CYLIND,0.0006,0,-0.0113,-0.08,0,360

VGLUE,1,2,3                                 !进行粘接
VGLUE,4,5,6
BLOCK,-0.05,0.05,-0.2,0.2,0,-0.08                       !建立空气模型
   
VOVLAP,ALL         !布尔叠分
NUMCMP,VOLU       !压缩体号

VSEL,S,VOLU,,5     
VATT,5,1,1        !设置导体5属性
VSEL,S,VOLU,,3     
VATT,4,1,1        !设置导体3属性
VSEL,S,VOLU,,4     
VATT,6,1,1        !设置导体4上导体属性
VSEL,S,VOLU,,1     
VATT,2,1,1       !设置铁磁属性
VSEL,S,VOLU,,2     
VATT,2,1,1       !设置铁磁属性
VSEL,S,VOLU,,6     
VATT,2,1,1       !设置铁磁属性
VSEL,S,VOLU,,7   
VATT,1,1,1       !设置外层空气属性
VSEL,S,VOLU,,3
SMRt,1               !智能划分等级1
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
VSEL,S,VOLU,,5
SMRt,8                !智能划分等级8
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
VSEL,S,VOLU,,4
SMRt,8                !智能划分等级8
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
VSEL,S,VOLU,,2
SMRt,6                !智能划分等级6
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
VSEL,S,VOLU,,1
SMRt,6                !智能划分等级6
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL

VSEL,S,VOLU,,6
SMRt,6                !智能划分等级6
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL

VSEL,S,VOLU,,7
SMRt,9              !智能划分等级9
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
              
NSEL,EXT        
D,ALL,AZ,0   
ESEL,S,MAT,,4                 ! 选择材料4也就是导体上的单元
BFE,ALL,JS,,0,0,354000000       !施加电流密度
ESEL,S,MAT,,5                 ! 选择材料5也就是导体上的单元
BFE,ALL,JS,,0,0,354000000       !施加电流密度
ESEL,S,MAT,,6                 ! 选择材料6也就是导体上的单元
BFE,ALL,JS,,0,0,354000000       !施加电流密度
ALLSEL,ALL
/solu
SOLVE                           !计算
/POST1
ESEL,S,MAT,,4   
ETABLE,FMAGX,FMAG,X
SSUM


得到的结果是材料4的导体受力为-0.8E-3，因为模型是对称的，所以我感觉这时的受力应该为0啊。
而分别将导体侧面的两个导磁体用空气单元导体，得到的受力依次为0.7E-4，和 -0.170458E-2。感觉这个结果不对啊，请问错在哪里了啊。结果的符号代表力的方向吗

[ 本帖最后由 fanyangyang 于 2007-12-21 10:05 编辑 ]

fanyangyang

大家给点建议啊

lilulilu

我运行的结果是
SUM ALL THE ACTIVE ENTRIES IN THE ELEMENT TABLE
TABLE LABEL     TOTAL
FMAGX    -0.976470E-19

你给的apdl运行产生的模型和你帖的模型不是一个吧。。

fanyangyang

我给的模型是没有包空气的，完全按照我帖的命令运行的结果如下啊，感觉差距太大了啊, 版主的结果就是我要得到的结果啊

SUM ALL THE ACTIVE ENTRIES IN THE ELEMENT TABLE

TABLE LABEL     TOTAL
FMAGX    -0.800666E-03

我用的是8.0，难道是版本的事，不会吧。好心的朋友也帮忙看看，谢了

[ 本帖最后由 fanyangyang 于 2007-12-19 09:43 编辑 ]

zbhit

我用8.0运行也是FMAGX    -0.800666E-03

giggsun

听说117在用铁磁材料是应该设置什么选项的

fanyangyang

版主，是这样的吗，我做的是静态，直流的磁场分析啊

fanyangyang

我又试了一下最简单的一根导体被空气包围的例子，
!/BACH,LIST
KEYW,MAGEDG,1
/PREP7
ET,1,117            
MP,MURX,1,1       !空气相对磁导率
MP,MURX,2,1      !导体相对磁导率
MP,RSVX,2,1E-8
CYLIND,0.01,0,0,0.2,0,360,  
CYLIND,0.1,0,0,0.2,0,360,
VOVLAP,ALL         !布尔叠分
NUMCMP,VOLU       !压缩体号
VSEL,S,VOLU,,1     
VATT,2,1,1         !设置导体属性
VSEL,S,VOLU,,2     
VATT,1,1,1   
ALLSEL,ALL
SMRt,6                !智能划分等级6
MSHAPE,1,3D
MSHKEY,0
VMESH,ALL
              
NSEL,EXT        
D,ALL,AZ,0   
ESEL,S,MAT,,2                 ! 选择材料2也就是导体上的单元
BFE,ALL,JS,,0,0,100000       !施加电流密度
ALLSEL,ALL
/solu
SOLVE                           !计算

/POST1
ETABLE,FMAGX,FMAG,X
SSUM

算出的力也不小啊
SUM ALL THE ACTIVE ENTRIES IN THE ELEMENT TABLE

TABLE LABEL     TOTAL
FMAGX     0.318665E-06
真是搞不懂了

[ 本帖最后由 fanyangyang 于 2007-12-19 20:19 编辑 ]

fanyangyang

我换了9.0，和11.0版本，问题解决了，想不懂为什么差这么大啊

maomaoface

为什么一根单独的导体通电也会有洛伦兹力啊？我试了一下，如果把jt载荷改为cp节点后加电流，洛伦兹力为0，这是为什么啊？

llkg

我又试了一下最简单的一根导体被空气包围的例子

=========================
你提到的这个例子我用11版本为何无法显示结果。
提示：
The request FMAG data is not available
The etable command is ignored

bftdlsyj

这个例子与导体的一端volt 为零，另一端耦合，在加电流，算力有什么不同吗？
还有个问题是：
在求力的时候，不需要施加力标志吗？

[ 本帖最后由 bftdlsyj 于 2008-7-29 10:58 编辑 ]

llkg

用117单元分析，可以不施加力的标志。
If using elements solid117,plane121,solid122,solid123(static analyses only), use emft to summarize electromagnetic force and torque. If you do use FMAGSUM, you do not need to first set either the Maxwell or the virtual work force flags via fmagbc. Electromagnetic or electrostatic forces are reported as FMAG_X, FMAG_Y, and FMAG_Z.
如果用上述单元，用emft可以计算电磁力和扭矩。如果用fmagsum，你无需通过fmagbc施加麦克斯韦或者虚功力的标志。电磁或者静电力记录在fmag_x，fmag_y,fmag_z中。

[ 本帖最后由 llkg 于 2008-7-29 13:49 编辑 ]