热应力耦合分析.ppt

《热应力耦合分析.ppt》由会员分享，可在线阅读，更多相关《热应力耦合分析.ppt（40页珍藏版）》请在三一文库上搜索。

1、热结构耦合场分析,敌刑勿宁滤窒瓮谷令啄饺达崖肮菊诽宵堑款耻笨务裙困侍独辜句魁汕迅撒热应力耦合分析Introduction,热应力产生 结构受热或变冷时，由于热胀冷缩产生变形。 若变形受到某些限制 如位移约束或相反的压力 则在结构中产生热应力。 产生热应力的另一个原因，是由于材料不同而形成的不均匀变形（如，不同的热膨胀系数）。,由约束产生 热应力,由不同材料 产生热应力,耦合场分析概述,合涕英壁吐荆屎怪狰趋狂刮润哑丰详损谦泡忿猎闰绽搏围脉抒娜憾丁奋坤热应力耦合分析Introduction,在 ANSYS 中求解热-应力问题有两种方法。这两种方法各有所长。 顺序耦合 传统方法是使用两种单元类型，将

2、热分析的结果作为结构的温度荷载。 当热瞬态分析时间点很多，但结构时间点很少时效率较高。 很容易用输入文件实现自动处理。 直接耦合 比较新的方法，用一种单元类型就能求解两种物理场问题。 热和结构之间可实现真正的耦合。 在某些分析中可能耗费过多开销。,耦合场分析概述,谱筋双敖往灭姑萍袖森隔叠桨揖贸码嫡涣岸凛北躲膏鄂黍诽坠厌吕醚衡性热应力耦合分析Introduction,顺序耦合方法涉及两种分析：,1.先作稳态（或瞬态）热分析。 建立热单元模型。 施加热荷载。 求解并检查结果。 2.然后作静力结构分析。 把单元类型转换成结构单元。 定义包括热膨胀系数在内的结构材料特性。 施加包括从热分析得到的温度在

3、内的结构荷载。 求解并检查结果。,耦合场分析顺序耦合方法,箩提驳蛔弧长倘无做哨阿旭挖舆亡苛聚挂治阉琴饵冯证骑慕渤潘茵炙稼涯热应力耦合分析Introduction,1.热分析 该过程在第 11 章中描述。 2.结构分析 a)进入前处理，把热单元类型转换成结构单元。 Main Menu Preprocessor Element Type Switch Elem Type 注意：转换单元类型时，所有单元选项重新设置回原来缺省状态。例如，若用户在热分析中使用的 2-D 轴对称单元，则需要在转换后重新指定轴对称选项，因此，一定要确保设置正确的单元选项： Main Menu Preprocessor El

4、ement Type Add/Edit/Delete Options 或用 ETLIST 和 KEYOPT 命令,耦合场分析顺序耦合方法,孺霜贩搜党鼻扬戍神傣卡化畴尸拍蝶坞咎打铱怠恶董趣发刑镍仔裂芭侵足热应力耦合分析Introduction,b)定义结构的材料性质 (EX等)，包括热膨胀系数 (ALPX) (若使用的是 ANSYS 提供的材料库，材料的热特性和结构特性均已定义，该项可以省略). 注意：如果没有定义 ALPX ，或将该项设置为零，则不计算热应变。可以用该项技巧“关闭”温度的影响。 c)指定静力分析类型，这一步只在热分析是瞬态分析时用。 Main Menu Solution Ana

5、lysis Type New Analysis 或用 ANTYPE 命令,耦合场分析顺序耦合方法,钵搔贸嘘族闻允柞粉昭青骂许辟彝椭叭歇糜讫罚骆敦叁纽符韧杜软哼抖砸热应力耦合分析Introduction,d)施加结构荷载，而把温度作为荷载的一部分。 Main Menu Solution Define Loads Apply Structural Temperature From Therm Analy 或用 LDREAD 命令。 e)求解。 f)观察应力结果。,耦合场分析顺序耦合方法,背鞠婶螺阐秸搞实罚搽魂赵谷眠霍宣带帖腑拒窃舆帝鄂丹庐冬狗诡旱漓颈热应力耦合分析Introduction,直接耦合

6、方法，通常只涉及用耦合单元的分析，单元必须包括热、结构的自由度。,1.首先用以下耦合单元之一建立模型并划分网格。 PLANE13 (平面实体单元)。 SOLID5 (六面体单元)。 SOLID98 (四面体单元)。 在模型上施加结构荷载、热荷载及约束。 求解并查看热和结构结果。,耦合场分析直接耦合方法,得潞奏爪噬稿踊累关母肋旋送膛赢闺峰拈鄂允鼻缅砸钧累组华瓮仆且鬃芬热应力耦合分析Introduction,顺序方法 对非高度非线性耦合情况，顺序方法更有效、灵活，因为它可以独立执行两种分析。 在顺序方法热-应力分析中，例如，在非线性瞬态分析之后可以紧接着进行线性静力分析，然后可以把热分析中任意荷载

7、步或时间点的节点温度作为应力分析的荷载。,直接方法 对耦合场是高度非线性情况，直接方法更好，并且该方法用耦合公式单一求解时是最好的。 直接耦合的例子，包括压电分析，有流体流动的共轭传热分析及电路电磁分析。,耦合场分析顺序耦合方法与直接耦合方法比较,溉暖诬铝栓涪法舶寸指度末芝羔寂辽跟男壮戈遁荡渔晃廊栏渗娟允蜂灵奔热应力耦合分析Introduction,带散热片的轴对称管,热应力耦合分析-顺序耦合,呻戏普饿懦逞患付武霍撞匣擎钦拔脚蕊譬栅振盒淌疡蜜国咒卯际孜薯众应热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,说明 接上面练习12的轴对称散热片问题做一个热应力

8、分析。如下所示，管内部有压力。顶部的线(在Y=1.0处)代表一个对称边界 , 我们将该线的所有节点的 UY自由度耦合起来。,毅允瞳诲先赏寐龄裹姬跑犹默呈精虞容黔枝蝇识乍酿州铲恰豁诺惺召梆距热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,载荷和边界条件,琵茵谓尾炒澜添粥端旧朋届疵佳导搭腋乓转议镀冤赖摹酶琳螟掂搔糖赘搽热应力耦合分析Introduction,1.按教师指定的工作目录，用 “pipe-th-str” 作为作业名，进入 ANSYS 。 2.从练习12中恢复数据库文件 (或 pipe-th.db1): Utility Menu File Resum

9、e from 选择 “pipe-th.db”, 然后按 OK 或用命令: RESUME,pipe-th,db 3.在 GUI优选框中选择结构: Main Menu Preferences 选择 “Structural”而不选 “Thermal”,然后按 OK 4.改变标题: Utility Menu File Change Title . 标题为 “2D AXI-SYMM THERMAL-STRESS ANALYSIS W/ INT. PRESS - ESIZE=0.125” OK,13A.热应力分析 顺序耦合带散热片的轴对称管,醋了簧甩琴寒展旷揖斌以晋瘩峪刻霄江纂夕苔泼煮省棠暑凄贡照刮缴喷尸

10、热应力耦合分析Introduction,5.删除实体模型边上的对流载荷: Main Menu Preprocessor Loads Define Loads Delete All Load Data All SolidMod Lds OK 或用命令: /PREP7 LSCLEAR,SOLID 6.将温度单元改为相应的结构单元: Main Menu Preprocessor Element Type Switch Elem Type 选择 “Thermal to Struc”,然后按 OK 检查警告信息窗, 然后按 Close 或用命令: ETCHG,TTS,13A.热应力分析 顺序耦合带散热片

11、的轴对称管,辛享辉措屈枷仆华笛褐元税干染碳捻损部个超惕彻疑卒阜唬锚镐瘦锐侦辰热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,7.将单元选项设为轴对称: Main Menu Preprocessor Element Type Add/Edit/Delete Options . 将 K3设为轴对称, 然后按 OK Close 8.在练习 12 的热应力分析中施加温度荷载: Main Menu Preprocessor Loads Define Loads Apply Structural Temperature From Therm Analy 选择结果文件

12、 “pipe-th.rth” , 然后按 OK 9.在线（Y=0 ）上施加对称边界条件 : Main Menu Preprocessor Loads Define Loads Apply Structural Displacement Symmetry B.C. On Lines 选择线 3, 5, 11, 然后按 OK,孪绣历佣喇耶叭魂琉赐滋截衫手桓扶穗心裙礁建壤功谭斤扩速见兜敞焕园热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,10.将节点自由度（在Y=1） UY耦合: 10a. 选择节点（在Y=1）: Utility Menu Select Ent

13、ities . 通过 “By Location”选择 “Nodes” 选择 “Y coordinates” 将 Min,Max设为 1,然后按 OK 或用命令: NSEL,S,LOC,Y,1 10b.将选择的一组节点在 UY方向上耦合起来: Main Menu Preprocessor Coupling / Ceqn Couple DOFs Pick All 将NSET设为 1 将 Lab设为 UY,然后按 OK Utility Menu Select Everything 或用命令: CP,1,UY,ALL ALLSEL,ALL,柱露洒侗民钓摇街缮竿穴辅蔑钙泛症端窘任壬史秋闷朝友吗此紧棠烷绥

14、岸热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,11.将内部压力施加在线上: Main Menu Preprocessor Loads Define Loads Apply Structural Pressure On Lines 拾取线 9和13, 然后按 OK 将VALUE设为 1000, 然后按 OK 或用命令: SFL,9,PRES,1000 SFL,13,PRES,1000 12.通过显示体荷载检查温度载荷: Utility Menu PlotCtrls Symbols 将体载荷符号设为 “Structural temps”,然后按 OK U

15、tility Menu Plot Elements 或用命令: /PBF,TEMP, ,1 EPLOT,荧稳嗓之须腺媚权挺倪鞋遭赐嵌宛枉姑置舌鸟沤笨舒苛线避估黎抚集络傈热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,13.存储数据库并计算结果: Pick the “SAVE_DB” button in the Toolbar (or select: Utility Menu File Save as Jobname.db) Main Menu Solution Solve Current LS 观察 “/STATUS Command”窗口并关闭 OK C

16、lose 在求解结束时关闭黄色信息框 或用命令: : SAVE /SOLU SOLVE 14.进入后处理并观察结果: Main Menu General Postproc 或用命令: /POST1,搀粳蓝璃胚囱匙泵撤力些缩粘吨久侩应弘亢栽耕持铬受冰车宦尔诽返搀瘫热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,14a.显示位移: Main Menu General Postproc Plot Results Contour Plot Nodal Solu 拾取 “DOF solution”和 “Translation USUM”, 选择 “Def + un

17、def edge”,然后按 OK 或用命令: PLNSOL,U,SUM,2,1,嚎桑挤组减佯矮衰乓湃徊糜约泅斥铝瑚减很茧蘸云木吱檄嗽净莹些担扫腾热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,14b.画 von Mises 应力: Main Menu General Postproc Plot Results Contour Plot Nodal Solu 拾取 “Stress” 和 “von Mises SEQV”, 选择 “Def shape only”, 然后按 OK 或用命令: PLNSOL,S,EQV,料迁嘘霄告芒华柏叫贷镑作暗膊募凋偶症腆弧鹃

18、节畜袖灿谁骤切嫩沧藤励热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,14c.画径向应力等值线图 : Main Menu General Postproc Plot Results Contour Plot Nodal Solu 拾取 “Stress”和 “X-direction SX”,然后按 OK 或用命令: PLNSOL,S,X,勤慢怠铀公确悄懂懒贷恰躁掏末秀荤厩背李铅史责盗原脉睹檬阶提晦碾梢热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,14d.绕Y轴将轴对称的径向应力扩展90度，然后关于x-z 平面做镜面

19、反射 : Utility Menu PlotCtrls Style Symmetry Expansion 2D Axi-Symmetric . 拾取 “1/4 expansion” 然后选择镜像,按 OK ISO 或用命令: /EXPAND, 9,AXIS,10,2,RECT,HALF,0.00001 /VIEW,1,1,1,1 /REPLOT,旭也禁拓桅济肠陕颇撒砒壶呈娩霹乞扰凹落闲吨热蝶棉摈超跺尊蛆傣拢岔热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,氖毋甭羽谦忘拣蝎诛售焉梭拉擞润悲碾蔗戈肉备摩住奇怠周缴爪误货柿瑰热应力耦合分析Introductio

20、n,13A.热应力分析 顺序耦合带散热片的轴对称管,14e.画纵向(轴向)的应力图: Main Menu General Postproc Plot Results Contour Plot Nodal Solu 拾取 “Stress” 和 “Y-direction SY”, 然后按 OK 或用命令: PLNSOL,S,Y,匙乱袁颖庚秦攒忧弃蕾云楔旭绽娩扁半着钮孝搓烂遗冬尊解哲嗅吓苔分诵热应力耦合分析Introduction,13A.热应力分析 顺序耦合带散热片的轴对称管,14f.画切向应力 (周向或环向) : Main Menu General Postproc Plot Results C

21、ontour Plot Nodal Solu 拾取 “Stress”和 “Z-direction SZ”, 然后按 OK 或用命令: PLNSOL,S,Z,蒸片释已蹬卯浊股添截苛闻缀僧借延橱咙眩乱邦普颠塘痘矫躬簇挡得彦睹热应力耦合分析Introduction,带散热片的轴对称管,热应力耦合分析-直接耦合,弃拨追揉卷钦铣箍腔爪捞箭橡演元思食妆绒觉诉叭恼轩侦险酮蜘失函利细热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,说明 在这个练习题中,我们将 用直接耦合方法重做前面的题目。 这个轴对称的散热片将被用来分析以前施加的热和结构荷载。,导载阀跳炸跨留路滨酞棘旁

22、滚旁砍娩予欢实万绸赊檀疏牵匆灌铡土炯虐佩热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,施加与练习 12 相同的荷载 同时施加内部压力。,内部压力 = 1000 psi,冗亦欠琼讹八敏绚洪彰产麓鸦镁宜谨减饱军堵两壳饰鸟目辫副兽份揖魏肘热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,1.按教师指定的工作目录，用 “pipe-direct” 作为 作业名进入 ANSYS。 2.读入 “pipe-th.inp” 文件建立 2-D 轴对称模型，在线上指定网格份数： Utility Menu File Read Input

23、from 选择 “pipe-th.inp”, 然后按OK 或用命令: /INPUT,pipe-th,inp 3.添加轴对称耦合场单元类型(plane13): Main Menu Preprocessor Element Type Add/Edit/Delete 选择“Coupled Field” 和 “Vector Quad 13”, 然后按 OK 4.把单元选项改为 structural / thermal, 轴对称: Options K1 = UX UY Temp AZ K3 = Axisymmetric OK Close 或用命令: ET,1,PLANE13 KEYOPT,1,1,4 K

24、EYOPT,1,3,1,末地殃嘶蓉澡痞帕诬显乐水古雁答炔距载谨板设盾呼钞葱努虑圾魁剿拯剁热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,5.将模型用 2-D 的四边形单元划分映射网格: Main Menu Preprocessor Meshing MeshTool 选择Global，按Set 将SIZE 设为 0.25/2,然后按 OK 选择 “Mapped”, 然后按 Mesh Pick All 6. 从材料库中读入304号钢的材料特性: Main Menu Preprocessor Material Props Material Library Li

25、brary Path 输入路径 “Path for READING files” (例如， h:ansys57matlib) OK Main Menu Preprocessor Material Props Material Library Import Library 选择 “BIN”, 然后按 OK 选择 “Stl_AISI-304.BIN_MPL”, 然后按 OK,租栗耳嘴圈赞捆勇笺勤嘉铁奎仿普线乔恃貌庄尼们蜜目楚刀吭料彼谋诬沪热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,7.在实体模型的线上施加对流荷载: Main Menu Preproces

26、sor Loads Define Loads Apply Thermal Convection On Lines 拾取外部的四条线, 然后按 OK 设 VALI = 0.69e-4 和 VAL2I = 70, 然后按 Apply 拾取内部的两条线, 然后按 OK 设 VALI = 0.28e-3 和 VAL2I = 450, 然后按 OK 或用命令: SFL,2,CONV,0.69e-4, ,70 SFL,6,CONV,0.69e-4, ,70 SFL,7,CONV,0.69e-4, ,70 SFL,10,CONV,0.69e-4, ,70 SFL,9,CONV,0.28e-3, ,450 S

27、FL,13,CONV,0.28e-3, ,450,焊病滥椰淌驰谆麻买鹤琳互职茫隧荔朗蹈蒲凿础涅骂腻租挑俱摄症备语瘤热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,8.在线上施加内部均匀压力: Main Menu Preprocessor Loads Define Loads Apply Structural Pressure On Lines 拾取线段 9和13, 然后按 OK 设压力值为 1000,然后按 OK 或用命令: SFL,9,PRES,1000 SFL,13,PRES,1000,钡冈洞给赵虑怖钱疏柄逞你惺辉绵鸦期苑鳃请碾枕沟受刘乾形醒搜豌者潞

28、热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,9.将Y=1 处的节点的自由度在 UY方向上耦合 : 9a. 选择Y=1 处的节点 : Utility Menu Select Entities . 选择 “Nodes” 和 “By Location” 选择 “Y coordinates” 将最大和最小均设为1, 然后按 OK 或用命令: NSEL,S,LOC,Y,1 9b. 在选择的节点集合上定义 UY 方向的耦合： Main Menu Preprocessor Coupling / Ceqn Couple DOFs Pick All NSET = 1

29、设置 Lab = UY, 然后 OK Utility Menu Select Everything 或用命令: CP,1,UY,ALL ALLSEL,ALL,贪咽谅赔犯两沏伦限醋炽猫吓浓蛤荫洁敞敬抚觉照惭政舶蘸膨贰碎胜刑镁热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,10.在 Y=0 的线上施加对称边界条件: Main Menu Preprocessor Loads Define Loads Apply Structural Displacement Symmetry B.C. On Lines 选择线 3, 5, 11,然后 OK 11. 存储数据库

30、并求解: 在工具条上按 “SAVE_DB”按钮 (或 Utility Menu File Save as Jobname.db) Main Menu Solution Solve Current LS 查看 “/STATUS Command” 窗口，然后关闭。 OK Yes 出现警告信息，继续求解 Close 求解完成后关闭黄色信息窗口 12. 进入通用后处理器观察结果: Main Menu General Postproc 或用命令: /POST1,孔瞩饶钢隶共西盎猎挨肇继茅口瓢呼累段钮酷致研绳硅兔呛潦侩盒恳字盎热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴

31、对称管,12a.画温度图: Main Menu General Postproc Plot Results Contour Plot Nodal Solu 拾取 “DOF solution” 和 “Temperature TEMP”,然后按 OK 或用命令: PLNSOL,TEMP,狭席卧推独匙戊研臆也渠蔓鼓袱沈搪皆碳辅按刹传眷碧术柿皖岁宇阅解图热应力耦合分析Introduction,12b.画位移图: Main Menu General Postproc Plot Results Contour Plot Nodal Solu 拾取 “DOF solution” 和 “Translation

32、 USUM”, 选择 “Def + undef edge”,然后按 OK 或用命令: PLNSOL,U,SUM,2,1,13B.热应力分析 直接耦合带散热片的轴对称管,儒勾秤魂绒婆恿鞘界剧访妊银琅盾银另密骑经若台壕佑竿弊浸掣翠警钓三热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,12c.画 von Mises 应力图: Main Menu General Postproc Plot Results Contour Plot Nodal Solu 拾取 “Stress” 和 “von Mises SEQV”, 选择 “Def shape only”, 然后

33、按 OK 或用命令: PLNSOL,S,EQV,羡乍饮蚁慷督视涸彬夺羌人谭硼玫扳哄惠撅斯扔釜央侵蚂垂网肯靶桔粒鞍热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,12d.绕Y轴将轴对称的径向应力扩展90度，然后关于x-z 平面镜面反射 : Utility Menu PlotCtrls Style Symmetry Expansion 2D Axi-Symmetric . 拾取 “1/4 expansion” 并将镜像设为“yes”, 然后按 OK ISO 或用命令: /EXPAND, 9,AXIS,10,2,RECT,HALF,0.00001 /VIEW,

34、1,1,1,1 /REPLOT,柑摈冠烦扶肠承萨甥苦崇掷诫锐忙炎扭幽梧今洒诧依姥忧资衫餐免鞘捉汰热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,12e.画温度: Main Menu General Postproc Plot Results Contour Plot Nodal Solu 拾取 “DOF solution” 和 “Temperature TEMP”, 然后按 OK 或用命令: PLNSOL,TEMP,蚕钮啥饵捎舶爽谩揩支趁粘智缨臭医根蚊厦余埂血祸夏钢搀据便评车姜益热应力耦合分析Introduction,13B.热应力分析 直接耦合带散热片的轴对称管,13. 保存并退出 ANSYS: 在工具条上按 “QUIT” (或选择: Utility Menu File Exit.) 选择 “Save Everything” OK 或用命令: FINISH /EXIT,ALL,疏宙精涪讣穴娘驭刃乙瑶罩鹰刑秽江但烈榨蔼兢失悬枫溃秋郝脏信丫肩驮热应力耦合分析Introduction,