fangdian 发表于 2011-8-29 15:36:54

还是关于电弧模拟的

大家好!
    还是发一个电弧模拟的问题,我用ansys做电弧模拟,想先模拟电场,再将电场计算的电流密度导入磁场进行计算,得出焦耳热,再将此焦耳热以生热率的方式加载到fluid142单元的flotran计算中,得出电弧的温度场,有关建模的帖子见我以前发的:
http://forum.simwe.com/thread-994870-1-1.html
我电场计算的结果感觉是对的,如下:

局部放大图:


磁场计算结果似乎也没有太大的问题:

磁场计算结果俯视图


但是流场计算后,温度场差距很大,最低温度和最高温度都太差,如下:


流场计算的速度场为:

流场速度场中的一个面局部放大为:


感觉差距太大,我看过人家文献中计算的温度场为:

我的问题可能在哪些方面?
感觉无从查起。

fangdian 发表于 2011-8-29 15:44:05

在ansys中,我感觉我没有好办法处理电弧的辐射问题,我查到了我所用其它的辐射损失(单位W/M3,随温度变化)


所以辐射我采用表格加载生热率的方式(将辐射损失值变为负的生热率,以bf方式加载到点上),如下:
*dim,radiat,table,73,1,1,temp,,   !等离子体辐射项

radiat(1,0,1)=5000,8098.290598,8397.435897,8696.581197,8995.726496,9294.871795,9604.700855,9903.846154,10202.99145,10502.13675,
radiat(11,0,1)=10801.28205,11100.42735,11399.57265,11698.71795,11997.86325,12297.00855,12596.15385,12895.29915,13205.12821,13504.2735,
radiat(21,0,1)=13803.4188,14102.5641,14401.7094,14700.8547,15000,15299.1453,15598.2906,15897.4359,16196.5812,16495.7265,
radiat(31,0,1)=16794.87179,17104.70085,17403.84615,17702.99145,18002.13675,18301.28205,18600.42735,18899.57265,19198.71795,19497.86325,
radiat(41,0,1)=19797.00855,20096.15385,20395.29915,20705.12821,21004.2735,21303.4188,21602.5641,21901.7094,22200.8547,22500,
radiat(51,0,1)=22799.1453,23098.2906,23397.4359,23696.5812,23995.7265,24294.87179,24604.70085,24903.84615,25202.99145,25502.13675,
radiat(61,0,1)=25801.28205,26100.42735,26399.57265,26698.71795,26997.86325,27297.00855,27596.15385,27895.29915,28205.12821,28504.2735,
radiat(71,0,1)=28803.4188,29102.5641,29401.7094

radiat(1,1,1)=0,-1335566.76,-2268284.633,-4756815.357,-12623258.01,-26864575.27,-54170060.-22,-110848146.2,-191051371.5,-311991459,
radiat(11,1,1)=-517041052.1,-784456182.4,-1196030567,-1677672305,-2296266893,-2877391359,-3587943375,-4745165140,-5417006022,-6526742603,
radiat(21,1,1)=-8299708372,-9290757363,-10148217728,-10976618837,-12227197960,-12287311518,-12968389593,-13553622272,-13487313507,-13487313507,
radiat(31,1,1)=-13822133392,-13687219403,-13487313507,-13553622272,-13620257037,-13620257037,-12968389593,-12841808759,-12904943977,-12904943977,
radiat(41,1,1)=-12968389593,-13225306655,-13553622272,-13421329145,-13890088260,-14804509950,-15472602389,-16170844255,-16817912863,-18190762563,
radiat(51,1,1)=-19293455187,-20462991134,-22461445394,-24776286181,-27062933282,-29851996081,-32607089391,-36144354123,-39869333689,-43123875632,
radiat(61,1,1)=-48749023556,-52987649376,-60787259812,-66397423113,-72170540429,-78831286752,-86106765088,-96862457770,-1.04257E+11,-1.12767E+11,
radiat(71,1,1)=-1.24389E+11,-1.37208E+11,-1.51349E+11

esel,s,mat,,1   
nsle,s
bf,all,hgen,%radiat%

不知道这样加载能否合适?
我采用的是bf加载到点上,也试过bfe加载到单元上,只是也是有问题(由于我先将磁场计算的生热率加载到单元上,如果再用bfe将辐射率加载到单元上会将原来的值替换掉,所以最终选择了bf)
唉,搞了好久了,就是没有出来,谢谢啊

wgq521 发表于 2011-8-29 19:44:06

你的阴极没有尖端,电弧在阳极与阴极之间,你这个是不是没有起弧?

wgq521 发表于 2011-8-29 19:46:37

而且你的思路,弧是由电极间的压差产生,产生电弧以后,温度升高,阳极与阴极之间导通,你单独计算电流的话,那温度怎么办呢,这是一个电磁场与流场耦合计算的问题,不能分别计算。

fangdian 发表于 2011-9-5 11:21:10

我初次计算的时候给氩气一个初始高温6000℃,就能引弧了,至于温度场,我打算通过电磁场计算出生热率导入流场计算温度,再把温度返回到电磁场,这样的话就能考虑温度了。

huangwei30 发表于 2011-9-8 14:18:34

楼主我也是做焊接电弧模拟的,只是现在懂的东西太少,希望能与楼主交流

taojunhust 发表于 2011-9-9 15:21:32

有些复杂,先学习了先。

fangdian 发表于 2011-9-8 20:44:09

huangwei30 发表于 2011-9-8 14:18 static/image/common/back.gif
楼主我也是做焊接电弧模拟的,只是现在懂的东西太少,希望能与楼主交流 ...

可以啊,我qq59158681

andylove55 发表于 2013-11-3 18:07:08

楼主我是个新手,能把你提到的文献名跟我说下吗,谢谢楼主

小七工作室 发表于 2013-11-4 10:22:52

我们计算过不少电弧 包括ansys高斯热源fluent激光焊接vof焊接等

simonsnow 发表于 2014-7-1 01:10:15

有些复杂,先学习了先。

fishgo 发表于 2014-7-4 09:46:03

到现在对电类还不怎么了解,好难的啊...楼主人才.

chenqian92 发表于 2014-7-8 15:16:40

支持一个,楼主我最近也在做电弧模拟,新手一个,以后希望多多指教

chenqian92 发表于 2014-7-8 15:23:01

楼主,我现在用热电耦合来仿电弧产生的整个过程,我想在这个热电场中加入外电路,请教一下如何实现?
或者说因为我在边界处加载的电流跟我计算过程中的场变量有关,我如何将这个场变量即时反馈过来?

神愤业炎 发表于 2015-6-10 23:01:38

好复杂的样子哦
页: [1]
查看完整版本: 还是关于电弧模拟的