comsol电磁分析加deform ALE动网格求助
本帖最后由 flyingmind 于 2009-3-13 19:00 编辑做个个简单的模型,就是一个载流导体,让它在空气中运动,然后观察空间磁场分布的变化,加了ALE的动网格,但发现物体运动了,但磁场并不对,场图是相对于直接将初始时刻磁场的保持不变,然后将导线的区域移动,留下空白,且空气磁场分布并未根据导线位置变化
有人用过ALE与电磁瞬态运动(最好直线运动)的麻烦指导下,谢谢 本帖最后由 flyingmind 于 2009-3-13 19:09 编辑
呵呵,后来仔细看了自带的旋转电机的例子,后来正确设置了identity pairs之后,网格可以动了,场分布也比较正常了,但似乎加了ALE之后,瞬态分析收敛比较困难? 空欢喜一场,上面的对旋转运动是适用的,对直线运动还有些问题,搞不定,有达人能指点下么? 现在发现其实如果做非旋转运动的动网格,需要将运动体用空气包围,然后空气包设置为自由运动,运动体设置自己想要的运动,然后在ALE边界条件里将运动的边界运动设置与运动domain一样,空气包边界勾选dx,dy的选项,且使其值均为0,同时还需要在绝对容差里分别设置Az,lmx的容差,其中x为数值,在solver manger中看求解量里可以看到,这样就没有时间步长逐渐变得很小的问题,可以正常计算
但是现在有个问题,就是经常运动一段后,会碰到inverted mesh的问题,用了参考文档里的解决方法,还是解决不了,苦闷中 关于inverted mesh,看了帮助,在较大变形时,网格发生畸变,目前我还没有解决研究了两天了,前几天一直在正直线运动的动网格,现在网格可以动了,我让网格以一定速度运行,一定速度下,如果位移很小就没有问题,位移一大就有问题了,运行一会总是出现inverted mesh,然后磁力线就在那附近严重畸变,看了帮助里说打开allow remesh,同时设置stop condition,停止后,然后创建最后正常时刻的deform mesh在当前的geom里,然后,更改起始时间为最后正常的时刻,再用restart求解,可以正常运行一会,再出现,再重复这个步骤,这样的话,在restart之前,要把之前的仿真结果自己导出,只是这样似乎太麻烦了些
看了有些人说网格不能剖太细,也不能太疏,太难弄了 后来仔细研究了问题所在,在5楼打开allow remesh on的基础上,将ALE的smooting method由laplace改为winslow,这样至少不会出现inverted mesh,它允许网格有较大的位移,但在位移非常大时,在移动方向的反方向,网格会被拉的很长,还有在运动体边沿的空气包网格也被拉长,这导致场分布精度不够了(网格太大),所以最终解决方法可能还是要winslow+运行一定位移后,停下来remesh,然后再restart接着之前时刻计算,当然每次restart之前的结果要自己保存,否则会丢失
发现smoothing method为laplace时,由于与运动方向一致的靠近运动体的网格几乎不动,所以导致运动体网格与相邻空气的网格交叉(这就产生的inverted mesh),而且我发现网格剖的越稀疏,产生inverted mesh的时刻就晚一起,刚开始不知道什么原因,原来网格大了,运动体要运动较长时间才会和相邻空气包形成inverted mesh 做的一个方形导线电磁场分析+直线运动的例子
(仅仅为了测试直线运动的deform ALE)是参考winslow smoothing method方法做的,由于位移较大,在位移大时,网格畸变较为严重,如下图所示:
所以应该用+stop condition限制最小mineqal1的值,在每次条件达到停下来时,重新划分网格(需要建立deformed mesh,然后重新划分网格),然后将起始时间设置为上次停止的时刻,用restart求解,这样效果会更好,每次在停止后,将结果导出为fem结构,名字取得不一样,最后可以用postmovie命令处理多个fem structure 具体可以参考附件的comsolrelease32a(有关于ALE ..),当然对高版本也使用,我用的是3.5 楼主厉害啊,我在研究永磁体直线运动的情况,整好有参考。谢谢! 我现在在做一点数值的介电常数,在某个特定区域运动。也是弄不明白移动网格怎么设置。下载看看,或许对自己有帮助。谢谢楼主了。 楼主你太自强不息了,好敬佩你,先下载学习下,回评感受 楼主很强大 不错,学习了 多谢楼主,对我的启发很大。类似的运动在ansoft中的实现非常简单,但在comsol中就有点复杂。ansoft是每步强制remesh的,comsol中麻烦就在它太灵活。呵呵,我刚学comsol,不知道能不能有机会和楼主多交流一下 本帖最后由 f418063285 于 2012-1-31 22:42 编辑
回复 7# flyingmind
您好,看了您关于动网格的例子,为什么我总是做不出你的结果呢?在create geometry from mesh 时总是报:Varying geometry shape order not supported。另外一设置绝对容差Az lmx根本计算不了。不设置可以。请您帮忙看一下,十分感谢。 看过学习就要顶 回复 7# flyingmind
flyingmind,你好,能不能看看您摸你的 这个文件呢,我现在一直在高动网格 ,几个星期了 搞不定,求助啊~~~~~ 回复 16# f418063285
哥们儿,你好啊,你的问题是否解决了呢? 回复 5# flyingmind
想问一下,网格边界运动和子域运动要设置相同,为什么我设的位移为0.2*t运行时弹出 矩阵错误啊,你是怎么设的啊?求教~~~!!!
页:
[1]
2