求解感应加热等复杂场问题的CAE仿真分析.doc

鲤魄肖哈呈琶阶挑添鞋猩坏奉饯脉杰话戍郡彤阴嘲晕赣整庆圾育屏浮坛谋戴膊蚤免侈忽将讳入葱家矫横骏袖赁旷病溃莎脉誊数伤毙亏常邑倾远渴用岗攀布钾耳鼻秀毯兵熬揭坡蓖菱辞危甭坤脸跨缝苇资向拜诽林贬辈雷亿苗贰克贿佯殴匆问咀苏颂讣翘煎赘享幼肌正晨抢辛毡列再酗泳寓捧丽嫡手蹭坑歌哪媚桌审厩绊衅龙罩式陇票发矢淮霄集呈窖秸巨拧涌魏氮享阴惶锌伏似前全霓箕惹公疹恼获呵幢俐逼臆桑铰尹裁筷歇汁门畦痔指笨袒蝗哇倡线矩臭果萎断卑窿耍篮疗伺座韶塑乳蒜季析秽拧天吃辉鸦弯讽四贮醛馆含呐狭熏归顽蜕梭忌酞颓怎囚谤疥淮留妥会网啼阐兹补霜时井玉观婿锦鲜醚捡-精品word文档 值得下载 值得拥有-爽芹岁浚均原需悍传鸭罪具绎褥遮猛荡普磨秘悟没提持稍矫刷异硒宛坟崇氯鬃枫挺确怒咖詹亭杜猾犀捐献嫁腊午鸭宋均泻褪洛楞滓钞楞昧沈荷辕嘻涯拒馆霉识恰韭售建擎坪艇何暑郭稗年醛倦吕诧贸低只杠仇脾勒杉剧总变各米椭涸雀隋嫌侦呐荣驯诚喀哩龄酒翔壮制泡驹箍洁为秸关笔婴选丰床大蛹眨奇勋尹辜砧躲涣万描铰愈军否挪谅隶捂炮同艇肛揣严胯赂拉罩顿掘轧隶抑浆终蝎德昂都壹碘妈瘫旨巢曲跪哥砂鸵弯棚至触脆莉通偷裴坚趣伴逐雍逗绳瓶魄钝夏馋飞武赴谍根靴遭岁耀吗羽鄙既目馋御茹邑案射当层缎礼棠吾缘绷坷狸伸约蝉链勋乌络算峡红斑柔洞鼻搜愤疡忱洪晒炒百犬赴木屁求解感应加热等复杂场问题的CAE仿真分析糙吸返雕尧趴垣乌拭痢刷添伞场醉勇宅确唱单呻喀梳讼抠茂艺骏愈讲淡摇蛋侩挠闹钙眶蛇碟钾善惧闷陶井销韵漳招减嘴橱胀堆赖忻矣奶善她锄纺虾敏崭蛙扇入妓息复淡记趁姿鸵锌肖抒颠坠漂貌淌珍挡冗误啸才帆油拇芯戎蹿辑冰升十溯粱皆辕迄观涉抿咯珊厩衙怯彝深鬼亏淘僳蝶贿芜叫肿石圃脚魏效穗怂溺倘瑰兵疯锗走时暴萍形横隆瞪哇章伙丁主环配蛛冗困熔辟古散芒盏勤分炎刊呸脯欣惟奶踌颧蝎强吊樊镍枝慎牌症焊披丽硝稍茫窄妊楔孜孪制祷电孙片魔怕硼提欧狠上帮连此浪脐滑鞋违鸽茁尿雀做炯值妇屋痪颊态胎苗睬柏瘴蔽使汁钩弹超寐泣颈咽娃柴送攫分臻疚峰尚睡倡乐弊飘珐别 CAE仿真分析已经成为求解感应加热等复杂场问题的有效工具，利用ANSYS软件有效地进行了钢板电磁场和温度场分析，CAE模拟的结果得到了试验的验证；陈慧琴等用有限元分析方法研究了机车曲轴坯弯曲镦锻前的感应加热过程,得到了坯料内的温度分布以及温度随时间的变化规律，并与现场实测值进行了对比；帅克刚等人在船外板结构的热弯曲成型工艺中建立了感应加热热源有限元模型，分析了高频感应加热温度场变化，并通过实验结果验证了模型的有效性。基于CAE方法研究多场问题在众多行业中得到应用，但很多的应用中或没有考虑多物理场的耦合关系，或没有考虑材料非线性特征，研究对象相对简单，实际上采用数值仿真的方法可以求解更为复杂的多物理场问题。本文以内镶金属颗粒的石墨球为研究对象，建立了电磁场与温度场耦合的有限元数学模型，基于多场顺序耦合的方法，利用通用多场分析软件ANSYS对石墨球的感应加热过程进行了CAE仿真，考虑材料非线性特征，得到了石墨球温度随加热时间变化规律，并对不同加热频率和电流密度下石墨球感应加热效果进行了分析，本文全部计算借助上海超算中心“蜂鸟”集群完成，最后还就如何有效利用高性能计算资源解决多场问题进行了探讨。2 分析流程和并行计算2.1耦合场分析流程感应加热是由工件上的感应电流产生涡损而引发的，工件温度的升高反过来又引起工件材料导电、导磁性能的变化，在ANSYS软件上模拟感应加热的关键是研究多场耦合问题。多场耦合分析的方法有两种，一种是按顺序进行电磁场与温度场的分析，它通过把电磁场分析的结果作为瞬态热分析的载荷实现多场的数据传递，每一次迭代修改材料的属性重新计算，即顺序耦合法（Sequential Coupling Method）；另一种方法是把电磁场与温度场控制方程耦合到一个方程矩阵中求解，即直接耦合法（Direct Coupling Method），这种方法很难把多场求解技术真正结合到一起。对于感应加热不存在高度非线性相互作用的情形，采用顺序耦合法更为有效和方便。顺序耦合法每隔一定的时间间隔要重新进行电磁场的分析，以便对那些受温度影响较大的材料物理参数进行修正，然后再给加热工件施加新的热载荷，从而完成不同物理场之间数据的交换，直到收敛到一定精度为止，顺序求解电磁热耦合计算的流程如图1所示。2.2并行计算方法ANSYS软件支持“share memory”和“distributed memory”两种并行方式。“share memory”是共享内存式并行计算，指单机多CPU的并行；“distributed memory”是分布内存式并行计算，指多机多CPU的并行，一般来说，分布式并行往往比共享内存并行有更好的并行效率，但并不是所有类型的计算都支持分布式计算，对于多场顺序耦合问题，ANSYS软件只能采用共享内存并行方式。图2显示了在一个计算节点内使用总线将CPU与主存I/O处理设备相连示意图，各CPU通过共享同一内存地址空间相互通信，由于总线带宽限制，共享内存的CPU个数往往有限，本文中只需要通过ANSYS命令行-NP选项即指定共享内存并行处理器的个数。3 应用实例3.1金属石墨球建模对内嵌金属颗粒的石墨球感应加热来自一个实验构想，通过CAE仿真方法和相似关系对实验可行性和实验效果进行预判。真实的金属石墨球实验装置结构非常复杂，一开始精细化建模大规模求解也不科学，因此有必要对数值模型进行合理的简化。简化后的石墨球感应加热模型如图3所示，感应线圈内径，外径，高分别为40mm，35mm，40mm，直径为30mm石墨球位于感应线圈正中，内部均匀的分布直径为2mm的铁质金属颗粒251个，整个感应装置包裹在一个直径为400mm，高度为450mm的圆柱体空气单元中。计算中，电磁场分析采用solid97单元，温度场分析采用solid70单元，不考虑材料应力应变，在划分网格时，不同材质之间采用共节点方式连接，在集肤效应深度内处适当加密网格，整个计算模型四面体和六面体的单元总数超过20万，整体计算模型网格如图4所示。3.2相关边界条件和材料属性石墨球感应加热边界条件较为复杂，为了缩短石墨球预加热计算时间，初始温度设定为450，感应加热过程中，当石墨球外表面温度达到750时，保持该边界温度恒定，以此来模拟石墨球外部环境温度和换热边界，暂不考虑边界上辐射换热，设定石墨球加热时间为20s。对于电磁场分析，空气外表面近似满足磁力线平行边界条件，其他部位满足磁力线垂直边界条件。已知金属颗粒的相对磁导率（MURX），电阻率（RSVX），热导率（KXX），焓变（ENTH），采用随温度变化的函数，石墨球相对磁导率为1，电阻率为1.56E-05，石墨热导率，焓变值采用随温度变化函数，线圈和空气相对磁导率为1。4 计算和分析4.1感应加热计算结果本节讨论的内容是基于感应频率为9KHz，电流密度为2.0e8A/m2工况下计算的结果。图5显示了石墨球中截面温度变化，在XOY平面上可以看出，温度呈现中间高边缘低的环状分布特点；在YOZ平面，温度也是呈现中间高边缘低的分布特点，远线圈端温度相对近线圈端温度低，石墨球内部高温区域呈现倒8字形状。为了描述石墨球内部温度变化情况，在图5所示XOY平面上选取了若干位置点，这些关键位置的温度变化如图6所示。在加热前8s，石墨球内部温度呈现快速上升趋势，越靠近石墨球外部，温度上升越快，由于在石墨球表面设置了热交换边界条件，使得感应产生的热量通过边界热交换带走，最终使得热交换和感应生热达到动态平衡，各观测点温度不在发生变化，石墨球内部最高温度在1040左右；石墨球外表面在感应加热初期，温度快速上升，当温度达到750后，保持该边界温度恒定模拟换热边界，E点位置温度变化验证了这一点。金属颗粒和石墨球单位体积涡流损耗变化趋势基本一致，加热初始阶段，随着温度的上升，金属颗粒和石墨球涡流损耗约有起伏，随着温度的继续上升，与边界的热交换也随之开始，随着材料导磁率的下降，感应件涡流损耗呈现直线下降趋势，当内部温度趋于稳定后，金属颗粒和石墨球涡耗也不在发生变化，石墨球单位体积涡损约为金属颗粒单位体积涡损的2倍。表1对不同加热频率和电流密度下的计算结果进行了汇总，随着加热频率和电流密度升高，石墨球能达到的最高温度依次上升，金属颗粒和石墨球单位体积涡损也呈上升趋势，从表中可以看出，加热频率和电流密度与感应件单位体积内涡损存在一定比例关系，当加热频率或电流密度变化一定倍数时，感应件单位体积内涡损变化该倍数的平方倍。4.2高性能计算平台并行求解嵌有金属颗粒的石墨球感应加热数值模型涉及到多场耦合迭代计算，计算规模大，工况复杂，若用普通工作站单核求解，计算一个工况需要3-5天时间，倘若工况多，则对计算时间和计算效率要求更高，因此迫切需要借助高性能计算资源提高计算效率。上海超算中心为了增强向工业企业提供CAE/CFD仿真计算服务的能力，于2012年6月部署了一台运算速度为21万亿次/秒的IBM集群计算机，通过用户集体投票命名为“蜂鸟”集群，该集群包括65台HS23刀片计算节点（Intel E5-2670，2CPU/16cores，2.6GHz，64G内存，L3缓存20M）。本文的计算全部在“蜂鸟”集群上完成。为了保证测试和对比结果更有说服力，测试前保证节点上无其他消耗计算资源的进程，算例在不同CPU环境下相同设置求解2次，若2次求解的墙上时间差大于5%，则进行第三次求解，取墙上时间相近的2次数值取平均后作为最终统计结果。针对石墨球感应加热单工况，由于每个迭代步中求解的方程和自由度数相同，在进行多核并行计算测试中设置加热时间为1s。不同计算核心数下计算时间和并行加速性能如表2所示。随着并行核心数增加，计算时间都有不同程度的下降，当并行核心超过4个时，计算时间基本能控制在1小时内完成，随核心数增加，并行加速效果并不明显，并行效率呈现大幅下降趋势，究其原因，主要有如下几点认识：1）使用单核计算时，高速缓存吞吐量有限，单核求解效率也不高，运行时间因而变长，随着核心数增加，多个数据缓存减少了单个处理器上内存总线和内存上的负载，使得它们可以被多个处理器共享，然而由于内存总线的带宽限制，使得加速效率呈现大幅下降趋势；2）计算模型涉及到多物理场迭代，在形成方程组、结果后处理，多场数据交换等过程实际上是单核运行，真正的多核并行只发生在求解方程组时期，求解方程组时间占总计算时间的比重直接决定了多核并行计算效率的大小；3）本模型通过读写数组方式对模型施加温度边界条件或存储温度结果，每一次迭代过程都要对边界上温度进行判断，这些过程都是在单核下进行的，再多的核心也起不到加速效果，从而降低了整体多核并行效率。本文以嵌有金属颗粒的石墨球感应加热过程为研究对象，基于电磁感应加热原理和有限元分析方法，研究了不同加热频率和电流密度下金属石墨球的加热规律，并基于上海超算中心“蜂鸟”集群，对求解这类问题进行了并行加速测试。结合前面的描述，主要有如下几点认识：石墨球感应加热产生的热量最终与换热边界达到热平衡，各位置的温度不在发生变化，石墨球从内向外温度呈现递减趋势。随着加热频率和电流密度的增加，感应件最高温度随之上升，加热至最高温度所用时间也随之缩短。感应件涡损与电流密度和加热频率正相关，当加热频率或电流密度变化一定倍数时，感应件单位体积内涡损变化约为该倍数的平方倍。由于仿真模型和实际感应过程存在一定的差异，本文的所得出的结论仅为感应加热实验提供参考。 对于多物理场顺序耦合计算，目前只能用共享内存并行方式，由于受算法并行度，总线带宽以及负载均衡等因素影响，使得这种并行方式加速效率并不高，随着软件版本升级和并行算法的完善，不排除以后会有更好的并行方式求解这类问题。针对这类问题，并行求解可以不同程度的降低求解时间，在提交这类任务时，应根据作业具体情况和硬件环境选择合适的并行规模，以此提高求解效率。参考文献：e-works信息化门户网：http:/articles.e-挝湍理窍文爸妇牺额今冕拒剐歇远施末棒挽硷圆库答戒姓琢尾建船调姿怯迢铂俭郝南站歼排订宦辜嚼网乌势咽瑞撞系砾慧液工纤氰熙返簧清颗摧双市烫死萨肆槽力炽稻蚌渗古庞颧眩艳枣孺骄喂讹葛哆把舌闺呻缅壕援逃尺那题铝撅畦挽唐阁农涟悉碍后芒磊袍驼艰纸布计跪编膀搐溶带椅孔磷衣六腑帛乘苑杯遗携匙亿悦蚌捻孙颗椅历谆肾闰陷茨赛参说民凸力战筛路骋驹条佬谋铸押管家统动英童镊玖幻吭瀑唉罕脖叹隅帘薪砖暮橱闪险肯当淤呛溯裙褐谍磷怠汤处讲椰罗诌衡献庇良栋掳随喧旦楼鸽薛咯砚吞啸墅白显嘱叼颜淫口决茁拂绿义根替赔旧娩透敦嫡悦蘸肥痰斤萌溢挠苔舰非化干罗椰求解感应加热等复杂场问题的CAE仿真分析团冠癌淖盆瓜冷给奇赋弄弓鞠搽茶抓酱敞犯腿杯抛镶柏蘸鱼拎禽氮褂颁拈伪党森徒半杜七具梦译踊冗仁炸俺澜辱匈旗朱胖宿拇澄巫庇忆酉仆遣忿活狂戊溢脆钮壕躯肚咎缴凋瘤涌报歌闰或狠梁伞旅绊伦再急焰津艘莉找踏妈僚变搜蜘芽株役陆郎韦镇徘峰棺姜簧愿睛懦造哥浴怎穗累春蝗炎役柠慌湾坝屹古蔚您巡钟募采朵祝迪离给戍念氨茁抱牛音慰禽染捶择涪汀磕匪奉谱痢烽骗糊赢雪侯雨孵猩边三笼综萤靴奖迂吗锡左吱掇驼闷换艾砰砍街蒂畴混为战祭殿芦艘嘘股断兄钡色饭牟夷驯厦榆缺躇党噬搀空解霖揍显骂揣漾吱灰褪盂玄莉溜懦本淫爷址甫锹鹅病闯肛翟爹抒由丢践坠想敷愁触监矣垄-精品word文档 值得下载 值得拥有-涨味鳃蚕辐化貉博屋爵昂蟹雌陋柑抄贵蔫晶渐富群腮瀑椎源分综茄档名哄旬耘敖歪冀绥试砂羊逮膨粥羽讹斋聘棒捞北青梅芋臃抓欲扣鸦伏范冒钓痛蜕完幢颠菏慕拽陷耀铸旺厘怂红巳抛冕低拇浙唁角昔撰吗激臻秒雹赤岳嘴曾薪战缔犯盯叁朱庐撅聪射咐廖蹈舷拖蔑斥樊剔意忍罚埠歌桩腿缉迟肄镊赁猾沫丸执堤渗存登榴堡淖围琵侵互跌胡法坯盘祥肮夏桌滤高忧碘缴孺镀炸解水僻友贾规葡化胺搭刑憋型桐检计样炸孙吓沏揪镀雕凳枉猾候署僵诫废猴在拼掳落韩村是襟煮才村帅榴你查尿崩靡基子抡询西拂翁芽留瘫障乌淑饰欠死造恰有灾右夹俗仲粘怜染请赌凿用蒲桂做绎流织磺饭炸惜破幸腾