97单元不能够用在空气+线圈+铁芯中么？

ponderadams

听说ANSYS公司的培训课件有这样一个说法


97单元是MVP，所以不能用在倒数第三行的情况下

在help里面并没有提到这一些，至少我没有看到，各位版主，不知道能否说明一下

xiaolanmao

这个我也看到了，好像用117和97 算出来的结果确实不一样阿，同问

小超

恩，我做出来的好像是不对

icemman

在铁芯-空气边界处,空气和铁芯的A值的切向分量相同,法向分量不同.可以通过复制边界的结点在加上CP来绕过这个问题,具体请看:
http://x.ansys.net/posting/50332

ponderadams

1) rotate nodal directions of nodes on the boundary so that one direction is
normal to the boundary
2) make a copy of the original mesh of the iron and delete the original mesh
3) couple the tangential components of the MVP at resulting coincident nodes
on the boundary, leave the normal component unspecified

怎么操作呀？

ponderadams

icemman  帮我解释一下好么，不胜感激

icemman

具体操作如下：
1）在空气-铁芯边界复制两场node,一层属于空气，另一层属于铁芯。
2）旋转node坐标系，使之与边界面的法向一致，用nora命令
3）用cpintf耦合AY,AZ

ponderadams

1) rotate nodal directions of nodes on the boundary so that one direction is normal to the boundary
2) make a copy of the original mesh of the iron and delete the original mesh
3) couple the tangential components of the MVP at resulting coincident nodes on the boundary, leave the normal component unspecified
从字面上翻译：
1.旋转空气-铁心交界边界的节点方向，使一个方向垂直于边界
2.复制铁心的原网格，并删除原网格。
3.对空气-铁心交界上的节点耦合节点法的正切分量，使不明确其垂直分量

楼上说的那种方法，我有些不明白
1.复制空气-铁心交界的节点，还是复制网格
2.另外这些操作都是在划分网格之后做吧？

各位高手，讨论一下吧

[ 本帖最后由 ponderadams 于 2007-9-18 09:12 编辑 ]

北极熊甲

XANSYS上你可能会搜索到一个叫 bill 的人写的一个宏。我的不知道放在哪里了，没找到。
内存的问题我也写道了FAQ里面。

icemman

Bill Bulat 用了3个macro来作这3步，我放在附件里了。其中旋转坐标系因为当时ANSYS的NROTATE功能还不完善，他用了比较麻烦的编程方法，速度较慢。我把我自己改进的方法也放在这里。

XXXX !选中边界的棱边，因为棱边AX,AY,AZ都要耦合
YYYY !选择边界面
nsll,s,1
cm,mubdry2n,node
nsla,s,all
cmse,u,mubdry2n
cm,mubdry1n,node
nora,all,1

XXXX和YYYY需要根据自己的问题写

ponderadams

icemman，你好，你的完整的程序是这样的么？
XXXX !选中边界的棱边，因为棱边AX,AY,AZ都要耦合
YYYY !选择边界面
！XXXX和YYYY需要根据自己的问题写
nsll,s,1
cm,mubdry2n,node
nsla,s,all
cmse,u,mubdry2n
cm,mubdry1n,node
nora,all,1
alls                ! UNSELECT "FLAT SURFACE" NODES, LEAVING CORNER
cmse,u,mubdry2n        ! AND EDGE NODES. COUPLE AX/Y/Z
cpintf,ax,1e-8
cpintf,ay,1e-8
cpintf,az,1e-8

alls                ! UNSELECT EDGE/CORNER NODES, LEAVING ONLY FLAT
cmse,u,mubdry1n        ! SURFACE NODES. COUPLE AY/Z ONLY
cpintf,ay,1e-8
cpintf,az,1e-8

我有几个疑问：
1.交界处的棱边上的节点mubdry2n和面上的节点mubdry1n都是空气和铁心的交界上的是属于空气的节点还是铁芯上的节点，我的意思是，这些交界上的点/单元的材料属性？

2.Bill Bulat 的3个宏，icemman好像没有调用，是吧？

3.在Bill Bulat 的copy_iron宏里定义的两个组件，mubdry0n和mubdry0e 在mbdryCP里面有说明
C***   1) COMPONENT mubdry0n (NODES ON BOUNDARY BETWEEN  DOMAINS OF DISSIMILAR mur EXCLUDING THOSE ON PERIODICITY SURFACES).交界面上相对磁导率不一样的那些节点
C***   2) COMPONENT mubdry0e ("IRON" ELEMENTS ON SURFACE)交界面上的单元？是么？

4.在Bill Bulat 的mbdryMVP宏里对mubdry0e和mubdry0n对应的交界面上重新定义mesh200单元并划分网格，把网格上的单元和节点分别组件为srfc_n和srfc_e，这里我想知道虽然有不同的组件名字，那mubdry0e和srfc_e代表相同的那些单元，节点也类似，所得明白一点，就是在esurf命令之后产生了新的单元和节点？是么？

5.在Bill Bulat 的mbdryMVP宏的后面一些程序是去似乎实现了nora的功能，是么？

6.另外最后一个宏mbdryCP，我猜大概的意思是耦合自由度了，但是我不是很理解这几句
alls  ! UNSELECT "FLAT SURFACE" NODES, LEAVING CORNER选中所有
cmse,u,mubdry2n ! AND EDGE NODES. COUPLE AX/Y/Z
cpintf,ax,1e-8
cpintf,ay,1e-8
cpintf,az,1e-8                                                          ！除了棱边上的节点耦合自由度
alls  ! UNSELECT EDGE/CORNER NODES, LEAVING ONLY FLAT选中所有
cmse,u,mubdry1n ! SURFACE NODES. COUPLE AY/Z ONLY
cpintf,ay,1e-8
cpintf,az,1e-8                                                          ！除了交界面上的节点耦合自由度                                               
a.为什么是除去了棱边上和交界面上的节点呢，不是正要对他们进行耦合自由度么
b.棱边是两个交界面的公共线，那么棱边上的节点同时也在交界面上，我想应该换一下顺序，相对交界面的节点耦合ay，az自由度，再对棱边上进行三个方向的自由度耦合。

6.Bill Bulat 和icemman的程序大体上是一样的思想，但是我觉得是有区别的，虽然都对mubdry0n进行nora操作，去别关键是在第4个问题上，是不是产生新的mesh200单元和节点，我不知道为什么Bill Bulat要这样做

7.这样做真得能使solid97单元能够弥补空气-铁芯求解的不准确的缺陷么？

[ 本帖最后由 ponderadams 于 2007-9-19 11:00 编辑 ]

icemman

你的代码没有复制节点，不会有效果。
我的程序调了bill的两个宏，mbdryMVP.mac用我列的代码代替了。mubdry2n，mubdry1n分别是边界面和棱边的点。边界面只耦合切向的自由度，棱边要耦合所有自由度。
其实只有SOLID97会有这个问题，如果可能还是用SILID98或SOLID117，不过都有各自的限制条件比如无法处理涡流，无法耦合电路。
至于准不准可以自己拿平面问题参照，平面问题无此误差。

lilulilu

So if you have shared nodes on the boundary between iron and air, single
valued MVP on that boundary forces both the normal and tangential components
of B to be continuous. Physically, only the normal component on such a
boundary is continuous. The tangential component is discontinuous.

“边界法向连续，切向不连续”
和
4楼“在铁芯-空气边界处,空气和铁芯的A值的切向分量相同,法向分量不同.”
一样？

98应该和96一个情况，都是标量势

icemman

B=curl(A), B法向连续等价于A切向连续

ponderadams

你说的没有复制节点，是因为缺少下面这一段么？
alls
*get,nmax,node,,num,max
esel,s,mat,,6
nsle
cm,scrap_e,elem
cm,scrap_n,node
ngen,2,nmax,all
egen,2,nmax,all
cmse,s,scrap_n
cmse,s,scrap_e
edel,all
ndel,all
alls
numc,node
大概了解你的思路了，另外，还是十分麻烦你，有空帮我回答一下我的7个问题（11楼）

[ 本帖最后由 ponderadams 于 2007-9-24 16:02 编辑 ]

ponderadams

顶一下，麻烦icemman 解释一下7个问题（11楼）


由于我铁芯比较复杂，经过这几天的努力，终于把棱边上的线、面选中了

对于自由度的问题，是针对哪些节点？看bill的程序，有点不清楚
对于复制节点/网格的问题，我想问一下，是复制整个铁芯上的节点/网格，还是仅仅那些交界面上的节点/网格

[ 本帖最后由 ponderadams 于 2007-9-27 09:49 编辑 ]

icemman

贴一下你的代码,或许可以帮上忙

ponderadams

XXXX !选中边界的棱边，因为棱边AX,AY,AZ都要耦合
YYYY !选择边界面
！XXXX和YYYY需要根据自己的问题写
nsll,s,1
cm,mubdry2n,node
nsla,s,all
norl,all,,1                                                 !这一条是我加上的，因为相对于下面nora的那一指令
cmse,u,mubdry2n
cm,mubdry1n,node
nora,all,1
alls
*get,nmax,node,,num,max
esel,s,mat,,6
nsle
cm,scrap_e,elem
cm,scrap_n,node
ngen,2,nmax,all
egen,2,nmax,all
cmse,s,scrap_n
cmse,s,scrap_e
edel,all
ndel,all
alls
numc,node

alls                ! UNSELECT "FLAT SURFACE" NODES, LEAVING CORNER
cmse,u,mubdry2n        ! AND EDGE NODES. COUPLE AX/Y/Z   ？？？？？这一条为什么不是cmse,s,mubdry2n
cpintf,ax,1e-8
cpintf,ay,1e-8
cpintf,az,1e-8

alls                ! UNSELECT EDGE/CORNER NODES, LEAVING ONLY FLAT
cmse,u,mubdry1n        ! SURFACE NODES. COUPLE AY/Z ONLY？？？？？这一条为什么不是cmse,s,mubdry1n
cpintf,ay,1e-8
cpintf,az,1e-8

可能我没有说明，bill的mbdryMVP.mac用你的一段程序代替了，我把代替的那段程序补上去了，剩下的步骤就是耦合自由度、复制和删除网格和节点了
但是我看了bill的程序，我不知道这两部是针对那些节点和网格

[ 本帖最后由 ponderadams 于 2007-9-27 15:59 编辑 ]

icemman

等有空贴个例子

[ 本帖最后由 icemman 于 2007-9-27 21:41 编辑 ]

ponderadams

请问各位版主：

norl和nora的用处