ansys如何实现时步有限元仿真电机旋转暂态

acl10

这些天一直在学习用ansys的APDL语言实现异步电机的时步有限元，也就是用ansys仿真电机暂态运动，遇到了一些问题想请教一下各位高手。我在网上找到了几篇关于ansys仿真电机的资料，分析《异步电动机综合物理场分析与Ansys实例计算》文章发现其只有建模和剖网格的APDL指令流，而关于求解过程代码并未贴出。后来我又在网上检索了很长时间，找到的资料都没有具体讲如何实现时步有限元的求解。
在ANSYS中实现时步有限元主要有以下几个步骤（参照《异步电机时步有限元模型在ansys中的实施》）：
1、首先计算场分析中无法自动计及的集中参数,即计算定子绕组的直流电阻、绕组端部的漏电感，以及转子端环电阻和漏感。
2、把计算出的集中参数加入到场路耦合模型中的电路方程式中。
3、完成一次某个时步tn下的非线性二维瞬态场的计算。
4、对每个时步的计算结果进行必要的后处理,电磁转矩的求解可以用麦克斯韦应力张量或虚位移的方法。
5、利用机械运动方程,可求得转子角加速度、角速度以及角位移,为转子转动到下一个位置做好准备。
6、利用上一步的计算结果将转子单元网格整体旋转到一个新的位置,并重新建立耦合与约
束方程。转子的旋转用耦合与约束方程的方法,其基本思路是:定子及气隙的一半组成一个整体,转子及气隙的另一半组成另一个整体,两个整体分别独立进行剖分,目的是使定、转子之间可以相互移动而不影响对方的单元网格。但是,这将在共同半径的圆周上形成不一致的网格,为了正确解算,还必须利用耦合与约束方程以及有限元的形函数,把气隙共同半径圆周上节点的自由度进行关联。
7、时间变量tn增加Δt,返回第三步,继续求解下一时步tn+1,直至时间变量到达设定的终值tend
个人认为，实施的关键为第六点，如何在ansys中实现插值独立坐标系法？其中，将气隙剖分成定子侧和转子侧容易实现，但是旋转之后的节点如何耦合不太明白，而且如何施加约束方程？希望各位高手给予帮助！！！

pengchangjun

论文中的内容是忽悠人的吧，

pengchangjun

论文前面写一堆自己都没整明白的公式，后面计算也没应用这些公式，再用商业仿真软件弄弄，自欺欺人，掩耳盗铃