水轮发电机瞬态分析求感应电动势

cqlyj2008

回复 20# wah_zyy
wah_zyy ：
你在第一步用timint，off的时候，耦合了电感单元两端的节点电压没？我这里用timint，off的时候，出现了需要耦合电感单元两端需要耦合电压的提示，说是因为timint，off相当于静态分析，电感中的电流为0，所以就要耦合。但是如果需要耦合的话，我发现很麻烦，而且在第一步进行耦合后，第二步就要用antype，trans，rest，这样就还得删掉先前的耦合，然后在计算，工程量很大啊，电脑慢。你有没有遇到过这样的问题，交流一下。版主，也麻烦您，能不能讲讲，这种情况是不是像上面我说的那样做？

cqlyj2008

回复 20# wah_zyy
wah_zyy ：
如果我在第一步就用瞬态求解的话，在旋转完转子后，出现了The value of VOLT at node 82540 is 1.914572563E+43.  It is greater than
the current limit of 1.E+25.  This generally indicates rigid body      
motion as a result of an unconstrained model.  Verify that your model   
is properly constrained.  
这样的错误，我看了下，这些点都位于转子线圈面，转子线圈时绞线圈，也耦合了emf和volt，不知道为什么会出现这样的情况？

wah_zyy

回复 22# cqlyj2008


    你好，我没出现你说的这些问题啊，你说的第二个问题我想问一下，励磁线圈我是耦合的emf 和curr自由度，没耦合volt，我定子侧也是这样耦合的，还有你的励磁加的是电流源还是电压源。我以前是加的电压，在计算中励磁侧电流密度经常变化，感觉不对，就改成加电流源了。但是现在还是无缘无故的跳出，没提示，好愁人

cqlyj2008

回复 23# wah_zyy
恩， 刚刚说错了。励磁线圈上面是用电流源的，也是耦合emf和curr。我想问的是，在第一步求解时用了timint，off没，还有就是在绕组连线的时候，用了电感单元没有？

wah_zyy

回复 24# cqlyj2008


     第一步求解时用timint,off了；
   我在绕组连接时用到电感单元了，是在绕组端部用的，模拟端部漏抗。到没有出现你说的问题，希望斑竹能关注一下，给予帮助。

cqlyj2008

回复 25# wah_zyy 这样的话，你耦合了电感两端的节点volt了吗？在静态分析的时候听说电感相当于短路，所以要耦合。现在我的模型上就出现了这个问题。下面是最后几行代码，您帮我看看，同时期待版主指导一下》
finish
/solu
*if,i,eq,0,then
  antype,trans,new     !求解静态
  timint,off
*else
  antype,trans,rest    !求解瞬态
  timint,on
*endif
!trnopt,full      !全波求解
eqslv,front       !选择稀疏矩阵求解器
!!!!!!
!!定义边界条件
!!!!!!
csys,1
nsel,s,loc,x,r1-delta,r1+delta
d,all,az,0
allsel
csys,0
/pnum,type,1
/number,1
NROPT,AUTO
ALLSEL
time,(i+1)*0.02/steps
NSUBST,1
KBC,0   
NEQIT,40
CNVTOL,A,,1E-2  
outres,all,all
solve

wah_zyy

回复 26# cqlyj2008


   不好意思，忘了回答你了，我的程序中没有对电感单元进行volt自由度耦合，也没有出现你说的问题，可能是我还没有遇到吧。
  看了你的程序，你应该是进行的非线性分析，我是的线性分析，所以你遇到的问题我还没遇到，不过我以后估计也会做非线性的，希望你的问题能早日解决，到时候还得请教你，呵呵
  我是用的阶跃加载方式，KBC,1，你用的是斜坡，这个会不会有影响呢？你试一下吧，我学的也不深入，刚学一段时间，不能给你太大帮助，不好意思啊。

沙隆巴斯

楼上的同志们讨论很热烈，最近实在很忙，现在让话题延续下去。
   之所以在场路耦合分析的时候需要求一个静态，就是timint,off，为得是让所有的数值有一个初始值，姑且称为a0，这样让计算有一个起点，当开始计算以后会逐渐计算到稳定的值。
   电感确实需要耦合，尤其是在初始的时候，个人不建议使用电感单元，这会带来计算上的复杂，当然，这也和使用者的水平有关。如果要求不是那么精确的话，端部电感就可以忽略，倒是电阻可以使用。
   至于上面数值计算过大的问题，肯定有地方不对
   上面给出代码的同学，一般边界条件等等约束都是在求解之前给出就不做修改的
   关于ce的删除，确实很多需要删除掉，这个没问题，只是记得重新耦合就是，这样才能旋转呐。

cqlyj2008

回复 28# 沙隆巴斯
版主你好：
    你的意思是边界条件应该放在循环前面，是吗?以前确实不知道这个问题。谢谢您。
    关于计算数值过大的问题解决了，是因为我的模型出了问题。现在又有新的问题来了。在第一步进行静态求解的时候，我把电感耦合了volt。收敛也挺快的，但是到第二步瞬态求解的时候，出现了不收敛的情况。我看了下命令行，主要是因为我设置的了收敛值。如果不对CURR和EMF设置，那就没有问题，但是我觉得这样结果不对，如果都设置的话，那么EMF老是达不到要求，求解就不收敛。不只何解。我的步长每次1个机械角。敬请版主指导：
CNVTOL,A,,1E-3
CNVTOL,CURR,1E-3
CNVTOL,EMF,1E-3

wah_zyy

版主你好，我的问题还是没有解决，我想问一下，什么原因会使得程序跳出而且没有错误提示呢？我的程序中间气隙还是分了200段，如果按每步走2段的话，可以计算到转子转动一个周期，即100步，但是下一周期就会出错；当按每步走1段时，转子转4步程序过一会跳出，而且都没有错误提示。希望版主能多多指教，万分感谢！

cqlyj2008

回复 30# wah_zyy 换了电脑还是这样吗？这样就比较难找原因了。你有用CNVTOL,A,,1E-3
CNVTOL,CURR,1E-3
CNVTOL,EMF,1E-3
这几条命令吗？我用了，结果到第二步emf就会不收敛，不知道为什么/

wah_zyy

回复 31# cqlyj2008


    我还没换电脑试，想找找自身原因，可是总找不到，等我试试看。还有我是做的线性分析，还用定义收敛容差吗？我看书上针对非线性有收敛容差的设置，你这样做应该没错吧

沙隆巴斯

容差这东西实在不好说，如果不设定软件有本身的默认值，一般设定小是为了求得更精确的解，但是有的时候会由于太小而找不到解，这就需要经验了。电感，我建议别一上来就加上,可以先算不带电感的，算的差不多了再把电感考虑进去。一个端部效应，个人觉得不会有决定性影响，但是对于计算过程来说增加的麻烦就大了。

吕鑫明

回复 32# wah_zyy
如果你写过数值计算程序就会知道，容差还是必须要设的，如果不设，默认值太小，程序反复迭代怎么也收敛不到允许的误差范围内，系统就会瘫掉，就自己跳出了。一开始把容差设大一点，收敛之后再减小一点。

wah_zyy

回复 34# 吕鑫明


    嗯，谢谢指导，我试试看。我算的场路耦合程序转子能旋转一个周期，但是第二个周期就会跳出，而且算的数值很大，比正常值大2个数量级。

wah_zyy

wah_zyy

上面的图形解释：为得到定子槽中一个导体感应电动势波形，只简单建立了一个槽中导体的场路耦合，即把该导体与一个无穷大电阻相连，取无穷大电阻上的voltage（ETABLE,voltage,SMISC,1），求得如上波形，横坐标为步数，纵坐标为电压值。转子一个周期共转100步，第一个周期能顺利完成，在第二个周期的第20几步（上图只记录到120步），程序跳出，且无提示。
问：
    2为什么感应电势波形不是上下对称的正弦波？
    3问什么求得的结果这么大？与正常值不符，大2个数量级？
希望大家能给予帮助，谢谢！

沙隆巴斯

1都酸了一个周期还跳出来有点不可理解
2感应电势里面可能会有谐波分量
3跳出来的问题就是由于值超大。

wah_zyy

回复 38# 沙隆巴斯


    版主好，谢谢您的回答。我也不知道该怎么办了，很是困惑，我在改改吧。版主，您有没有时间帮我看看程序，是不是我的程序有问题，有时间的话我发到您的QQ邮箱里吧。没有时间也没事，总之很感谢您的指导，谢谢！

沙隆巴斯

程序一看就晕，还是随便聊聊吧。