Samcef 热烧蚀及热固耦合分析
Samcef Mecano 客户和应用案例
EADS Astrium 公司曾作为协助方之一与LMS SAMTECH 公司合作开发了Amaryllis 软件，而后不断帮助Amaryllis 软件在整个欧洲航天领域中的推广和使用，并且极其欣赏Amaryllis 软件分析结果的精确性给EADS Astrium 宇航研究带来的便利。如今Amaryllis 软件已在其热防护领域被大范围使用，而最初时EADS 集团利用Amaryllis 软件进行仿真实验相关性分析时就已经对软件的计算能力和精确性赞叹不已。 在AURORA 返回舱烧蚀热防护设计项目中，AURORA 返回舱结构呈倒置的伞状，返回舱在中心的球罐结构中被保护起来，再入大气层速度高达13km/s ，热流达6~9MW/m2
。
EADS ASTRIUM 采用SAMCEF Amaryllis 和BOSS Quattro 相结合的方式，进行了热防护材料的优化选择和厚度设计，在考虑重量、温度和烧蚀厚度等约束的情况下，实现了返回舱不同部位分布采用不同材料的优化设计方案
图中所示，样件在等离子束喷射的环境下进行试验，样件为轴对称结构，另外就是需要考虑等离子束的偏离入射。实验进行之后，用户测量了样件的几何尺寸变化、残渣、气体量及能量，并从中评估了此烧蚀材料的各项主要性能指标。
烧蚀样件实验示意图
而后开始在Amaryllis 环境下建模并完成热解烧蚀计算，下图中是建模时输入的烧蚀材料密度分布，以及计算结果烧蚀结束时的温度和材料密度分布。
样件
2D
轴对称初始模型
样件2D轴对称模型（网格已退化）计算结果温度分布（左）及密度分布（右）
而后，EADS Astrium使用Amaryllis软件负责欧洲AURORA返回舱项目中对烧蚀热防护结构的分析。AURORA返回舱结构呈倒置的伞状，返回舱在中心的球罐里被保护起来。
瞬态过程中随驻点压力变化的温度历程曲线及物质消耗率曲线与实验结果的对比
AURORA返回舱热防护结构示意图
AURORA
返回舱热流加载
AURORA 返回舱烧蚀结果温度分布
不同位置处温度历程曲线
出于保密原因，不便将此处各项实验结果展示。但是各项结果的精度均让EADS和其他用户非常满意。
网格退化过程不同位置表面结点位移曲线
ＭＢＤＡ案例－导弹弹头热防护设计
MBDA是欧洲著名的导弹及防务系统设计制造商，其对导弹的热防护研究水平处于世界前列。弹道导弹，尤其是弹头位置，热防护工作至关重要，MBDA导弹弹头热防护设计采用的分析工具也是SAMCEF Amaryllis。下面介绍MBDA集团利用SAMCEF Amaryllis进行导弹弹头热烧蚀建模和分析的案例（此处几何模型被简化）
MBDA
弹头热防护实验及建模示意此处使用的材料为碳化复合材料，密度大约1925kg/m3。对弹头位置，需要考虑对流换热、辐射换热，以及在很强压力下的化学和机械烧蚀作用。
有限元网格模型及材料化学烧蚀速率特性曲线
输入了复杂的材料特性参数以及其他边界条件后，Amaryllis 为MBDA 集团提供了所关心的完整的仿真结果。
瞬态模拟中不同位置的热交换系数和热焓随时间和位置变化的曲线
烧蚀结束时的温度场及不同位置处温度变化曲线
烧蚀结束时的几何形态及不同位置处位移变化曲线
欧盟第七科技框架计划（ＦＰ７）项目（ＡＥＲＯｃａｐｔｕｒｅ ｆｏｒ Ｆｕｔｕｒｅ ｓｐＡｃｅ ｔｒａｎＳｐｏｒＴａｔｉｏｎ）
AEROFAST 项目是欧盟第七科技框架计划项目，旨在设计一种航天器，能够用于人类探索太阳系，作为一种先进的运输工具，实现人或物在月球、火星或GEO 、LEO 卫星与地球间的往返运输。要设计这样的一种航天器，热防护系统就显得尤为重要。
在AEROFAST 项目中，热防护系统采用酚醛树脂复合材料，承受组分为CO2的火星大气层热流环境，采用SAMCEF Amaryllis 软件进行热烧蚀分析，分析模型采用3D 热烧蚀模型。
AURORA返回舱热防护设计
EADS ASTRIUM作为Samcef Supervisor软件的开发合作方，也是软件的忠实用户。Samcef Supervisor被广泛用于EADS ASTRIUM的多个火箭和再入航天器的热防护设计项目。在AURORA返回舱热防护设计项目中，AURORA返回舱结构呈倒置的伞状，返回舱在中心的球罐里被保护起来，再入大气层速度高达13km/s，热流达6-9MW/m2。EADS ASTRIUM 采用Samcef Supervisor和Boss Quattro相结合，进行热防护材料的优化选择和厚度设计，在考虑重量、温度和烧蚀厚度的约束情况下，最终实现了返回舱不同部位分布采用不同材料的优化设计方案。
太阳能帆板展开机构热固耦合分析
对于太阳能翻帆板展开机构的分析，如何保证机构的可靠性并防止展开过程中与周围部件产生碰撞尤为重要，因此需要非常精确的力学模型，充分考虑柔性、非线性及热等因素。ESA和AREVA等公司运用Samcef Supervisor成功模拟了太阳能帆板在温度因素影响下，非线性机构的展开过程，并且仿真和实验体现出良好的一致性。
火箭发动机热-固耦合分析
多种液体火箭发动机推进器的热-结构耦合分析。
如考虑热辐射和对流载荷以及碳环的非正常磨损
的液氧涡轮泵喷管时域非线性的瞬态分析。
另外还进行了考虑3D热传导、热对流和热辐
射在内的包括H2/O2供给、增压、管路排空和涡
轮泵冷却等多种分析。并且考虑燃烧过程中可能存在很大的振动和高频现象，结合
Samcef线性和非线性求解器也可以进行耦合分析。
ANSYS Example07热-结构耦合分析算例 (ANSYS)
07 热－结构耦合分析算例(ANSYS) 在土木工程结构中，温度应力在很多情况下对结构的影响很大。很多时候需要先对结构进行热传导分析，得到结构内部的温度应力分布，再进行结构分析，得到由于温度产生的结构内力。ANSYS提供了很方便的热分析－结构分析切换工具，本节将以一个圆环的热应力分析为例，介绍ANSYS提供的相关功能。 (1)首先进行热分析，进入ANSYS主菜单Preprocessor->Element Type->Add/Edit/Delete， 选择添加单元为Quad 4 node 55 号热分析单元 (2)进入ANSYS主菜单Preprocessor->Material Props->Material Models，添加热传导速率 参数Thermal->Conductivity->Isotropic，设定热传导速率为0.07。添加力学属性Structural->Linear->Elastic->Isotropic，设定弹性模量为30e9，泊松比为0.2。添加热膨胀系数Structural->Thermal Expansion->Secant Coefficient->Isotropic，设定热膨胀系数为1e-5。 (1)开始建立模型。还是按照ANSYS标准的点、线、面、体建立模型。首先建立关键点。 在ANSYS主菜单Preprocessor->Modeling->Create->Keypoints->In Active CS，输入以下关键点信息 (2)下面开始建立弧线。在ANSYS主菜单Preprocessor-> Modeling-> Create-> Lines-> Arcs-> By End KPs&Rad，首先点选关键点2和3，然后点选中心点1，最后输入半径为5，生成第一个圆弧。接着点选关键点4和5，然后点选中心点，输入半径8。生成第二个圆弧 (3)在ANSYS主菜单Preprocessor->Modeling->Create->Lines->Straight Line，连接关键 点2，4和3，5。组成圆环轮廓 (4)在ANSYS主菜单Preprocessor->Modeling->Create->Arbitrary->By Lines，点选圆环周 边轮廓线，生成圆环面。 (5)下面划分网格，由于本模型只有一种单元一种材料，所以不必复杂的设置属性。进入
热-电直接耦合实例分析
热-电直接耦合实例分析 上海交通大学（单位）江丙云（姓名）分析目的： 电导产品的温升 分析类型： 热-电耦合 分析对象： 电子连接器产品 知识要点： 热传、电导材料 热-电场直接耦合 稳态和瞬态求解 1.1 热-电直接耦合稳态分析实例 本节以图1-1所示的电子产品的端子通电升温为例，采用稳态求解，详细讲解多物理场的热-电直接耦合。 1.1.1 问题描述 图1-1（a）所示为电子产品用于导通电流1A的两个端子：端子1（Terminal-1）和端子2（Terminal-2）；端子1&2材料为铍铜C7025-TM02，其热传导系数为0.19W/mm- ℃，电传导IACS%为45%，即电传导系数为26100/(ohm-mm)， 工程实际中，端子1的B处与端子2的C处受力接触，造成接触阻抗R C（Contact Resistance），而端子自身阻抗称为体阻抗R b（Body Resistance）；简化有限元模型，接触用Beam连接代替，接触阻抗通过赋予Beam对应的材料实现，稳态求解即可。
整个电路导通如图1-1（b）所示，电流从A点导进，经Beam连接，至D点流出。 （a）端子几何模型 （b）电路导通示意 图 1-1导通端子有限元模型 1.1.2 导入几何部件 ●创建、保存模型 打开Abaqus/CAE，创建Model Database: With Standard/Explicit Model，应用File→Save as命令保存模型为9.3_T-E.cae。 ●重命名模型 如图9-9所示，右击树目录的Model-1，重命名（Rename）Model-1为Model-1_T-E。 ●导入几何部件 应用File→Import: Part命令，或右击树目录Parts→Import，导入光盘中本节几何文件：T-E.igs，其余，如图9-10默认选项即可。
第19章热-结构耦合分析
第19章热－结构耦合分析 热－结构耦合问题是结构分析中通常遇到的一类耦合分析问题。由于结构温度场的分布不均会引起结构的热应力，或者结构部件在高温环境中工作，材料受到温度的影响会发生性能的改变，这些都是进行结构分析时需要考虑的因素。为此需要先进行相应的热分析，然后在进行结构分析。热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度（热通量）等。本章主要介绍在ANSYS中进行稳态、瞬态热分析的基本过程，并讲解如何完整的进行热－结构耦合分析。 19.1 热－结构耦合分析简介 热－结构耦合分析是指求解温度场对结构中应力、应变和位移等物理量影响的分析类型。对于热－结构耦合分析，在ANSYS中通常采用顺序耦合分析方法，即先进行热分析求得结构中的温度场，然后再进行结构分析，且将前面得到的温度场作为体载荷加到结构中，求解结构的应力分布。为此，我们需要先了解热分析的基本知识，然后在学习耦合分析方法。 19.1.1 热分析基本知识 ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数。ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外，还可以分析相变、有内热源、接触热阻等问题。 热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热对流是指固体的表面与它周围接触的流体之间，由于温差的存在引起的热量的交换。热辐射指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。 如果系统的净热流率为０，即流入系统的热量加上系统自身产生的热量等于流出系统的热量：q流入+q生成-q流出=0，则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。 瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。
[热弹耦合稳定性问题探究,稳定性,结构,问题]关于结构的刚
关于结构的刚-热弹耦合稳定性问题探究 1 引言 高超声速飞行器因为其特殊的应用价值而成为当前国际研究的热点, 而随着航天飞行器速度的提高, 气动加热现象趋于严重. 为了得到更准确的分析结果, 航天飞行器设计中越来越多的开始计入热-弹耦合效应的影响[1,2]. 国内外学者已经进行了有益的研究, 可供我们参考和借鉴. Krishma 和Bainum[3]改进了模拟暴露于太阳辐射下的自由梁和板的偏差的计算方法, 考虑了表面反射率和入射角的影响. Thornton 和Kim[4]对哈勃太空望远镜的太阳帆板进行了简化, 考虑了支撑杆弹性变形与热辐射强度的耦合, 研究了结构的热诱发振动现象, 得到了系统稳定性的判据. Johnston 和Thornton[5]在温度变化确定的情况下, 用有限单元法研究了热效应对航天器姿态动力学的影响. Johnston在建模过程考虑了刚体运动和弹性变形的耦合, 但是基于线弹性假设, 忽视了几何非线性效应. 为了考虑温度变化对刚- 柔耦合动力学性态的影响, Oguamanam 等人[6]对受热载荷且存在平面运动的柔性体进行了研究, 采用16 节点的等参拉格朗日插值的单元, 在温度变化规律为已知的情况下分别建立了中心刚体-曲梁系统和中心刚体-薄板系统的动力学模型, 考虑了几何非线性效应, 但未考虑温度和变形的耦合, 也没有对温度变化规律未知的情况做进一步研究. Cannarozzi 和Ubertini[7], Kidawa-Kukla[8],Hosseini-Tehrani 和Eslami[9], Mahi 等人[10]都在这一研究领域做了有益的工作. 由于我国航天工程发展的需求, 我国学者在这一领域也做出重要贡献. 刘锦阳和洪嘉振[11]同时考虑了几何非线性和热效应, 用假设模态法对各柔性梁进行离散, 从虚功原理出发, 建立了带集中质量的柔性梁系统的动力学方程, 研究了几何非线性和热效应的综合影响. 王捷和刘锦阳[12]以哈勃天文望远镜(HST)为研究对象, 基于柔性多体系统动力学理论,考虑了柔性附件弹性变形引起的热辐射边界条件的变化, 建立了中心刚体和太阳能毯柔性附件多体系统的刚-柔-热耦合的动力学方程. 用假设模态法进行离散, 对HST 双侧太阳翼的振动特性进行分析, 研究系统各特征参数对于柔性附件热颤振的影响. 潘科琪和刘锦阳[13]研究了在热冲击下任意形状(仅一个方向有曲率)复合材料壳的非线性刚-柔耦合动力学响应. 根据Mindlin 理论, 建立了任意形状的复合材料壳的非线性应变-位移关系. 借助于数学理论以及几何关系, 描述了壳上任意点的变曲率. 用虚功原理建立了动力学变分方程, 并采用等参单元对壳的连续动力学方程进行离散, 建立了中心刚体-复合材料壳的刚-柔耦合动力学方程, 分析了在热冲击作用下复合材料壳的线性、非线性的动力学特性, 以及曲率、材料特性对动力学响应的影响. 蒋卓良[14]研究了太阳帆板在热载荷作用下的固有频率和动态响应.史晓鸣和杨炳渊[15], 王宏宏等人[16], 李忠学和严宗达[17]也在这一研究领域做了有益的工作. 将航天动力学处理为刚-弹耦合动力学问题取得成功. 由于高超声速飞行器的飞行马赫数高, 气动加热效应大, 在飞行过程中承受着严酷的气动力载荷和气动热载荷, 因此, 本文将之处理为刚-热弹耦合动力学问题. 本文应用刚-热弹耦合动力学变分原理来研究结构的刚-热弹耦合稳定性问题.
Maxwell与Fluent电磁热流耦合分析介绍
14.5耦合实例4——Maxwell和FLUENT电磁热流耦合 例， 14.5.1 析钢块在上述工况下的温度场分布情况、风的流线图及风的温度分布云图。 图14-164几何模型 14.5.2软件启动与保存 Step1：启动Workbench。如图14-165所示，在Windows XP下单击“开始”→“所有程序”→ANSYS14.0→Workbench 14.0命令，即可进入Workbench主界面。 图14-165 Workbench启动方法 Step2：保存工程文档。进入Workbench后，单击工具栏中的按钮，将文件保
存为“MagtoThemtoFluid”，单击Getting Started窗口右上角的（关闭）按钮将其关闭。 注意：本节算例需要用到ANSOFT Maxwell14.0软件，请读者进行安装； 由于ANSOFT Maxwell软件不支持保存路径中存在中文名，故在进行文档保存时，保存的路径不不能含有中文字符，否则会发生错误。 14.5.3导入几何数据文件 Step1：创建几何生成器。如图14-166所示，在Workbench左侧Toolbox（工具箱）的Analysis Systems中单击Maxwell 3D并按住左键不放将其拖到右侧的Project Schematic窗口中，此时即可创建一个如同EXCEL表格的项目A。 Step2：双击A2（Geometry）进入如图14-167所示的电磁分析环境，此时启动了Maxwell 3D软件。 图14-166项目A Step3：依次选择菜单Modeler→Import，在出现的Import File对话框中选择ThermaltoFluid.x_t几何文件，并单击打开按钮。 图14-167电磁分析环境 Step4：此时模型文件已经成功显示在Maxwell软件中，如图14-168所示，同时弹出Modal Analysis对话框，在对话框左侧的栏中显示的几何图形为Good表示数据读取无误，单击Close按钮。
换热器热应力耦合有限元讲解
第一章 课题相关知识介绍 2.1散热片知识 散热片是一种给电器中的易发热电子元件散热的装置，多由铝合金，黄铜或青铜做成板状，片状，多片状等，如电脑中CPU 中央处理器要使用相当大的散热片，电视机中电源管，行管，功放器中的功放管都要使用散热片。一般散热片在使用中要在电子元件与散热片接触面涂上一层导热硅脂，使元器件发出的热量更有效的传导到散热片上，在经散热片散发到周围空气中去。 2.1.1散热片的材质比较 就散热片材质来说，每种材料其导热性能是不同的，按导热性能从高到低排列，分别是银，铜，铝，钢。不过如果用银来作散热片会太昂贵，故最好的方案为采用铜质。虽然铝便宜得多，但显然导热性就不如铜好（大约只有铜的50%左右）。 目前常用的散热片材质是铜和铝合金，二者各有其优缺点。铜的导热性好，但价格较贵，加工难度较高，重量过大（很多纯铜散热器都超过了CPU 对重量的限制），热容量较小，而且容易氧化。而纯铝太软，不能直接使用，都是使用的铝合金才能提供足够的硬度，铝合金的优点是价格低廉，重量轻，但导热性比铜就要差很多。有些散热器就各取所长，在铝合金散热器底座上嵌入一片铜板。 对于普通用户而言，用铝材散热片已经足以达到散热需求了。 北方冬季取暖的暖气片也叫散热片。 散热片在散热器的构成中占有重要的角色，除风扇的主动散热以外，评定一个散热器的好坏，很大程度上取决于散热片本身的吸热能力和热传导能力 2.1.2散热片结构的设计 1. 肋片的散热量 肋基导入的热量向肋端传递，经肋片传给流体，因此肋片得热平衡方程为： 肋基导入的热量Φ=Φ流体带走的热量λ 所以肋片向流体的传热量恒等于肋基截面上导入的热量，根据傅立叶定律得 每片等截面直肋散热量的计算式为： )(1)(0mH th m h m h mH th m A H H λ λθλ++ =Φ （2—1）
热结构耦合
第21章热－结构耦合分析 热－结构耦合问题是结构分析中通常遇到的一类耦合分析问题。由于结构温度场的分布不均会引起结构的热应力，或者结构部件在高温环境中工作，材料受到温度的影响会发生性能的改变，这些都是进行结构分析时需要考虑的因素。为此需要先进行相应的热分析，然后在进行结构分析。热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度（热通量）等。本章主要介绍在ANSYS中进行稳态、瞬态热分析的基本过程，并讲解如何完整的进行热－结构耦合分析。 21.1 热－结构耦合分析简介 热－结构耦合分析是指求解温度场对结构中应力、应变和位移等物理量影响的分析类型。对于热－结构耦合分析，在ANSYS中通常采用顺序耦合分析方法，即先进行热分析求得结构的温度场，然后再进行结构分析。且将前面得到的温度场作为体载荷加到结构中，求解结构的应力分布。为此，首先需要了解热分析的基本知识，然后再学习耦合分析方法。 21.1.1 热分析基本知识 ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数。ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外，还可以分析相变、有内热源、接触热阻等问题。 热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热对流是指固体的表面和与它周围接触的流体之间，由于温差的存在引起的热量的交换。热辐射指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。 如果系统的净热流率为0，即流入系统的热量加上系统自身产生的热量等于流出系统的热量：q流入+q生成-q流出=0，则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。 瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。
加热弹性圆板的大振幅自由振动
收稿日期:1996-12-24 *机械部教育司科研基金资助项目 加热弹性圆板的大振幅自由振动* 李 世 荣 (甘肃工业大学基础课部,兰州 730050) 摘 要 基于v o n K árm án 理论和Hamilton 原理,导出了均匀加热弹性圆板用中面位移表示的大振幅自由振动动力学控制方程.并在调和振动模态假设下,采用Kan-to rovich 平均方法将所得混合初-边值问题转化为相应的非线性常微分方程两点边值问题,采用打靶法和解析延拓法,分别获得了不可移简支和夹紧加热圆板非线性振动的调和振动响应,绘出了不同加热温度下的幅-频特征曲线.得出:升温使圆板的固有频率降低,从而实现改变板的温度对其固有频率的控制. 关键词 加热圆板 大振幅振动 打靶法 解析延拓 固有频率分类号 O.343 受约束弹性构件在变温环境中的振动和热屈曲问题的研究对有效的工程设计具有十分重要的意义.关于薄板构件的热弯曲、热屈曲和热振动问题的研究状况可见Tauchert (1991)和 Thornto (1993)的综述文章 [1,2] ;文献[3]中采用Galerkin 法研究了面内热应力作用下的圆薄板在热屈曲构形附近的微幅振动;在文献[4]中采用混合方法分析了变温产生的内力对矩形板 和斜形板微幅振动的影响;Tani [5]采用有限差分格式研究了横向均布载荷作用下的加热环板在轴对称大变形弯曲平衡构形附近微幅振动的特征值问题;本文作者也曾采用摄动法与有限差分法有机结合的方法分别讨论过等厚度和变厚度环板在均匀变温场内的非线性振动问题 [6,7] .本文将以横向和径向位移作为基本未知量,讨论了周边受约束圆板在静态变温下的大 振幅自由振动. 1 控制方程 考虑一半径为a ,厚度为h 的周边不可移圆形薄板,设板从自然状态起的均匀升温为T ,分析该加热圆板的轴对称大振幅自由振动.忽略面内位移的惯性项,由Hamilto n [8] 原理可导出用中面位移表示的系统vo n K árm án 型的无量纲动力学控制方程: 4 w x 4+2x 3 w x 3-1x 2 2 w x 2+1x 3 w x - 2 w f 2+λx x x w x =12W 2 x x u x +12 w x 2 +_x u w x (1) 2 u x 2 +1x u x -u x 2+ 2 w x 2 w x +1-_2x w x 2 =0 (2) 第23卷第2期1997年6月甘 肃 工 业 大 学 学 报Jour nal of Gansu U niv er sity o f Technolog y V o l.23N o.2 Jun.1997
热-结构耦合(单元转换)
ANSYS热－结构耦合分析实例 在土木工程结构中，温度应力在很多情况下对结构的影响很大。很多时候需要先对结构进行热传导分析，得到结构内部的温度应力分布，再进行结构分析，得到由于温度产生的结构内力。ANSYS提供了很方便的热分析－结构分析切换工具，本节将以一个圆环的热应力分析为例，介绍ANSYS提供的相关功能。 (1) 首先进行热分析，进入ANSYS主菜单Preprocessor->Element Type->Add/Edit/Delete，选择添加单元为Quad 4 node 55 号热分析单元 (2) 进入ANSYS主菜单Preprocessor->Material Props->Material Models，添加热传导速率参数Thermal->Conductivity->Isotropic，设定热传导速率为0.07。添加力学属性Structural->Linear->Elastic->Isotropic，设定弹性模量为30e9，泊松比为0.2。添加热膨胀系数Structural->Thermal Expansion->Secant Coefficient->Isotropic，设定热膨胀系数为1e-5。 (1) 开始建立模型。还是按照ANSYS标准的点、线、面、体建立模型。首先建立关键点。在ANSYS主菜单 Preprocessor->Modeling->Create->Keypoints->In Active CS，输入以下关键点信息 (2) 下面开始建立弧线。在ANSYS主菜单Preprocessor-> Modeling-> Create-> Lines-> Arcs-> By End KPs&Rad，首先点选关键点2和3，然后点选中心点1，最后输入半径为5，生成第一个圆弧。接着点选关键点4和5，然后点选中心点，输入半径8。生成第二个圆弧 (3) 在ANSYS主菜单Preprocessor->Modeling->Create->Lines->Straight Line，连接关键点2，4和3，5。组成圆环轮廓 (4) 在ANSYS主菜单Preprocessor->Modeling->Create->Arbitrary->By Lines，点选圆环周边轮廓线，生成圆环面。 (5) 下面划分网格，由于本模型只有一种单元一种材料，所以不必复杂的设置属性。进入ANSYS主菜单Preprocessor->Meshing->Size Cntrls-> ManualSize->Global->Size，在Global Element Size窗口中设置单元尺寸为0.5。在ANSYS主菜单Preprocessor->Meshing->Mesh->Areas，点选圆环进行网格划分 (6) 下面首先进入热分析，进入ANSYS主菜单Solution->Analysis Type->New Analysis，设置分析类型为稳态分析Steady-state
热结构耦合分析的例子
这是两个同心圆，我画的不是很圆，请大家见谅。外圆外边温度70o 内圆内边温度200 求圆筒的温度分布，径向盈利，主环向应力 /batch,list /show /title,thermal stress in concentic cylinders-indirect method /prep7 et,1,plane77,,,1 mp,kxx,1,2.2 mp,kxx,2,10.8 rectng,0.1875,0.4,0.05 rectng,0.4,0.6,0,0.05 aglue,all numcmp,area asel,s,area,,1 aatt,1,1,1 asel,s,area,,2 aatt,2,1,1 asel,all esize,0.05 amseh,all esize,0.05 amesh,all nsel,s,loc,x,0.1875 d,all,temp,200 nsel,s,loc,x,0.6 d,all,temp,70 nsel,all finish /solu solve finish /post1 path,radial,2 !设置路径名和定义路径的点数 ppath,l,,,0.1875 ！通过坐标来定义路径 ppath,2,,0.6 pdef,temp,temp ！温度映射到路径上 T0
paget,path,points,radial ！用数组的形式保存路径 plpath,temp finish /prep7 et,1,82,,,1 mp,ex,1,30e6 mp,alpx,1,0.65e-5 mp,nuxy,1,0.3 mp,ex,2,10.6e6 mp,aplx,2,1.35e-5 mp,nuxy,2,0.33 nsel,s,loc,y,0.05 cp,1,uy,all nsel,s,loc,x,0.1875 cp,2,ux,all nsel,s,loc,y,0 d,all,uy,0 nsel,all finish /solu tref,70 ldread,temp,,,,,,rth solve finish /post1 paput,path,points,radial pmap,,mat ！设置路径映射来处理材料的不连续 pdef,sx,s,x ！映射径向应力 pdef,sz,s,z ！映射环向应力 plpath,sx,sz ！显示应力结果 plpagm,sx,,node ！在几何模型上显示径向应力 finish 这儿是一个在热结构耦合分析的例子，大家有兴趣可以看看，我想同时问一下，cp 这个命令是什么意思啊
耦合摆小振动问题研究
耦合摆小振动问题研究 唐志根 （华南师范大学物理与电信工程学院08物理学专业3班20082301053）摘要：利用拉格朗日方程，给出了耦合摆在平衡位置附近的动力学方程，求解出耦合摆在平衡位置附近作小振动时的本征频率。并对耦合摆的振动坐标作适当的线性组合，得出了简正坐标和简正振动频率。 关键词：拉格朗日方程；小振动；耦合摆；简正坐标； Key words：Lagrange equation；small oscillation；couple pendulum；normal coordinate Abstract：The dynamic equations of the couple pendulum near the balance location were given by using Lagrange equation．The eigenfrequency of small oscillation of the couple pendulum near the balanced location was also solved．The normal coordinates and normal oscillation frequency were gained by making adaptablely linear combination with oscillation coordinates of couple pendulum． 前言：力学系统大多是较为复杂的非线性系统，不太可能对其运动得出完整的一般解。但对某些小振动下的非线性系统能通过某些过程给出解析解。耦合摆由两个单摆通过一弹簧相连构成，关于耦合摆的小振动在大学力学教材中曾求解过，而且简正坐标的求法不少文献作了讨论，但求解过程较繁。本文利用分析力学的拉格朗日方程，求解了耦合摆在平衡位置附近作小振动时的情况，给出了处理耦合摆小振动的一种较简单适用的方法。虽然只是两个自由度的力学体系，但它包含了处理非线性力学系在平衡位置，附近作小振动的主要过程。 1、耦合摆 如上图所示，耦合摆由两个相同的单摆平行悬挂，两个摆锤之间用一轻质弹簧相连构成。两球摆锤的质量同为m：摆长为l；弹簧劲度系数为k；弹簧的质量忽略不计。
ANSYS电热耦合分析
一、Electric-Thermal Analysis ANSYS中电热耦合分析主要焦耳热效应（Joule heating）、塞贝克效应（Seebeck effect）、珀尔帖效应（Peltier effect）、珀尔帖效应(Thomson effect)。我们这里的分析主要是Joule heating分析，即通电产生热量，用于加热双层薄片。 1. ANSYS电-热耦合知识点 1.1、Element DOFs选项：UX, UY, UZ, and TEMP： 可用于Thermal-Electric Analysis 的单元类型如上表所示，其中LINK68, PLANE67, SOLID69, and SHELL157 是专用的thermal-electric elements，专用于Joule heating effects，SOLID5, SOLID98, PLANE223, SOLID226, and SOLID227 则需要选择DOFs选项为TEMP and VOLT。 For SOLID5 or SOLID98, set KEYOPT(1) to 1； For PLANE223, SOLID226, or SOLID227, set KEYOPT(1) to 110。 1.2、Material Properties设置： 对于Joule heating effects，需要设置材料参数： 电学参数：electric permittivity电阻率RSVX、RSVY、RSVZ 热学参数：thermal conductivity导热系数KXX, KYY, KZZ 若考虑瞬态热效应，需设置密度DENS、比热C或焓ENTH 1.3、Load载荷设置： 设置Applied Voltage or Current 设置对流、辐射、传热等边界条件
热力耦合分析单元简介
热力耦合分析单元简介！ SOLID5－三维耦合场实体 具有三维磁场、温度场、电场、压电场和结构场之间有限耦合的功能。本单元由8个节点定义，每个节点有6个自由度。在静态磁场分析中，可以使用标量势公式（对于简化的RSP，微分的DSP，通用的GSP）。在结构和压电分析中，具有大变形的应力钢化功能。与其相似的耦合场单元有PLANE13、SOLID62和SOLID98。 INFIN9－二维无限边界 用于模拟一个二维无界问题的开放边界。具有两个节点，每个节点上带有磁向量势或温度自由度。所依附的单元类型可以为PLANE13和PLANE53磁单元，或PLANE55和PLANE77和PLANE35热单元。使用磁自由度（ＡＺ）时，分析可以是线性的也可以是非线性的，静态的或动态的。使用热自由度时，只能进行线性稳态分析。 PLANE13－二维耦合场实体 具有二维磁场、温度场、电场和结构场之间有限耦合的功能。由4个节点定义，每个节点可达到4个自由度。具有非线性磁场功能，可用于模拟B－H曲线和永久磁铁去磁曲线。具有大变形和应力钢化功能。当用于纯结构分析时，具有大变形功能，相似的耦合场单元有SOLID5、SOLID98和SOLID62。 LINK31－辐射线单元 用于模拟空间两点间辐射热流率的单轴单元。每个节点有一个自由度。可用于二维（平面或轴对称）或三维的、稳态的或瞬态的热分析问题。 允许形状因子和面积分别乘以温度的经验公式是有效的。发射率可与温度相关。如果包含热辐射单元的模型还需要进行结构分析，辐射单元应当被一个等效的或（空）结构单元所代替。 LINK32－二维传导杆 用于两节点间热传导的单轴单元。该单元每个节点只有一个温度自由度。可用于二维（平面或轴对称）稳态或瞬态的热分析问题。 如果包含热传导杆单元的模型还需进行结构分析，该单元可被一个等效的结构单元所代替。 LINK33－三维传导杆 用于节点间热传导的单轴单元。该单元每个节点只有一个温度自由度。可用于稳态或瞬态的热分析问题。 如果包含热传导杆单元的模型还需进行结构分析，该单元可被一个等效的结构单元所代替。 LINK34－对流线单元 用于模拟节点间热对流的单轴单元。该单元每个节点只有一个温度自由度。热对流杆单元可用于二维（平面或轴对称）或三维、稳态或瞬态的热分析问题。 如果包含热对流单元的模型还需要进行结构分析，热对流单元可被一个等效（或空）的结构单元所代替。单元的对流换热系数可分为非线性，即对流换热系数是温度或时间的函数。
ANSYS热应力分析 精选实例x
ANSYS热应力分析实例 当一个结构加热或冷却时，会发生膨胀或收缩。如果结构各部分之间膨胀收缩程度不同，和结构的膨胀、收缩受到限制，就会产生热应力。 热应力分析的分类 ANSYS提供三种进行热应力分析的方法： 在结构应力分析中直接定义节点的温度。如果所以节点的温度已知，则可以通过命令直接定义节点温度。节点温度在应力分析中作为体载荷，而不是节点自由度间接法。首先进行热分析，然后将求得的节点温度作为体载荷施加在结构应力分析中。 直接法。使用具有温度和位移自由度的耦合单元，同时得到热分析和结构应 力分析的结果。 如果节点温度已知，适合第一种方法。但节点温度一般是不知道的。对于大多数问题，推荐使用第二种方法一间接法。因为这种方法叮以使用所有热分析的功能和结构分析的功能。如果热分析是瞬态的，只需要找出温度梯度最大的时间点，并将此时间点的节点温度作为荷载施加到结构应力分析中去。如果热和结构的耦合是双向的，即热分析影响结构应力分析，同时结构变形又会影响热分析（如大变形、接触等），则可以使用第三种直接法一使用耦合单元。此外只有第三种方法可以考虑其他分析领域（电磁、流体等）对热和结构的影响。 间接法进行热应力分析的步骤 首先进行热分析。可以使用热分析的所有功能，包括传导、对流、辐射和表 面效应单元等，进行稳态或瞬态热分析。但要注意划分单元时要充分考虑结构分
析的要求。例如，在有可能有应力集中的地方的网格要密一些。如果进行瞬态分析，在后处理中要找出热梯度最大的时间点或载荷步。 重新进入前处理，将热单元转换为相应的结构单元，表7-1是热单元与结构
单元的对应表。可以使用菜单进行转换:
第20章 热-应力耦合分析实例
第20章热－应力耦合分析实例 由于温度的分布不均在部件内部会产生热应力，在结构分析中常会遇到需要考虑温度场对应力分布影响的情况。特别在进行各类燃机部件，如航空发动机的涡轮盘、叶片等的强度计算分析时通常要考虑热问题。还有一些输送管道由于内外温度不同也会产生热应力。另外材料的性能和其温度是相关的，不同的温度下其性能通常不同，这也会造成部件应力的变化。本章将通过实例来讲解如何用ANSYS6.1来进行这类问题的分析。 20.1 问题描述 一无限长的截面形状和尺寸如图20.1所示的厚壁双层圆管，其内外层温度分别为Ti 和To，材料数据和边条如表20.1所示，利用ANSYS程序来求解圆管沿径向的温度分布情况，并求解圆管内沿径向和周向的应力情况。 图20.1 双层管道的截面图 从上面描述的问题可以看出，本实例属于轴对称问题，我们可以采用轴对称方法来进行分析。同时本问题为典型的热－应力耦合问题，可以采用间接法顺序耦合分析的一般步骤进行分析。因为管道为无限长，故建模时轴向尺寸可以是任意大于零的值，且将其一边
轴向约束，一边所有节点轴向自由度耦合。下面我们将首先建立有限元模型，进行稳态热分析，并观察分析其沿径向的温度分布情况。然后将模型中的热单元类型转换称对应的结构分析单元类型，从新定义材料的力学性能参数，并将热分析的结果以体载荷的形式施加到模型中，定义合理的边界条件，进行结构静力求解。最后，观察并分析整个结构沿径向和周向的应力分布情况。 20.2 建立模型 在ANSYS6.1中，首先我们通过完成如下工作来建立本算例的有限元模型，需要完成的工作有：指定分析标题，定义材料性能，定义单元类型，建立几何模型并划分有限元网格等。本节中定义的单元类型和材料属性都是针对热分析的。下面将详细讲解分析过程。 20.2.1指定分析标题并设置分析范畴 在这一步中我们将指定本实例的分析路径、数据库的名称、分析标题。另外为了得到适合热分析的菜单选项，需要将分析范畴指定为热分析。 1．选取菜单路径Utility Menu >File >Change Jobname，将弹出修改文件名(Change Jobname)对话框，如图20.2所示。在输入新文件名(Enter new jobname)文本框中输入文字“CH20”，为本分析实例的数据库文件名。单击按钮，完成文件名的修改。 图20.2 修改文件名(Change Jobname)对话框 2．选取菜单路径Utility Menu >File >Change Title，将弹出修改标题(Change Title)对话框，如图20.3所示。在输入新标题(Enter new title)文本框中输入文字“Thermal Stress in Concentric Cylinders-Indirect Method”，为本分析实例的标题名。单击按钮，完成对标题名的指定。 图20.3 修改标题(Change Title)对话框 3．选取菜单路径Main Menu >Preference，将弹出菜单过滤参数选择(Preference of GUI Filtering)对话框，如图20.4所示。单击对话框中的Themal(热)选择按钮，选中Thermal选项，以便ANSYS6.1的主菜单设置为与热分析相对应的菜单选项。单击按钮，完成
O形橡胶密封圈的热应力耦合分析
O形橡胶密封圈的热应力耦合分析 研究原油高温热采工具O 形橡胶密封圈在高温高压下的密封特性。借助于 大型有限元分析软件ANSYS，建立O 形橡胶密封圈及其边界的二维轴对称有限元模型，研究油压、装配间隙和摩擦因数对密封面最大接触应力、剪切应力和Von Mises 应力的影响，并采用热应力耦合分析方法，分析温度对O 形密封圈密封 性能的影响。结果表明：摩擦因数对应力影响不大，而油压和装配间隙对应力影响很大，过大的装配间隙会造成O 形橡胶密封圈最大接触应力下降和最大剪切 应力上升，造成密封失效；当温度升高时，密封圈最大剪切应力和接触应力相应减小，而最大Von Mises 应力明显减小，因此应使O 形密封圈在适当的温度下工作，以确保密封的可靠性。 随着国民经济的快速发展，陆地石油资源的不断减少，海上稠油开采将会成为热点之一。海上稠油开采通常采用高温热采作业模式，因此对开采工具的高温高压密封性能要求较高。当前原油开采一般都是在150 ℃以下进行，开采工具 的密封部件基本上都是采用简单高效、可靠性高的橡胶O 形圈，而橡胶O 形圈的耐用性、可靠性以及密封性能对温度比较敏感，同时对安装结构尺寸也有一定的要求。O 形密封圈简称O 形圈，是一种截面为圆形的橡胶密封圈，其成本低廉、安装和使用方便，被广泛应用于机械、汽车、动力及石油化工等领域。 对于常温下橡胶O 形圈的密封性能，目前研究报道较多，而对于高温(150 ℃以上) 下密封性能的研究，公开报道尚不多见。研究高温( 在橡胶材料可用最高温度范围内) 和高压( 稠油开采要求最高压强达到20 MPa) 下橡胶O 形圈的热应力耦合特性，对提升稠油高温热采工具密封性能有着至关重要的作用，可为密封部件的选择和结构优化提供理论依据和数据支撑。 利用ANSYS 有限元分析软件对O 形橡胶密封圈进行应力分析，并采用热应力耦合分析方法，分析温度对O 形密封圈密封性能的影响，得到温度与应力变 化的关系。 1、模型的建立 1. 1、几何模型的建立 O 形橡胶密封圈工作时依靠密封圈发生弹性变形，在密封接触面上产生接触应力，当接触应力大于密封介质的压力时，则不发生泄漏。如图1(a) 所示，当O 型密封圈装入密封槽后，其截面承受接触压缩应力而产生弹性变形，对接触面
热力耦合分析单元简介
共享：热力耦合分析单元简介！ 挑选了部分常用的，希望能方便大家的使用，其中自己翻译了一部分，不准确之处还望见谅，大家还可以继续补充哦！： SOLID5－三维耦合场实体 具有三维磁场、温度场、电场、压电场和结构场之间有限耦合的功能。本单元由8个节点定义，每个节点有6个自由度。在静态磁场分析中，可以使用标量势公式（对于简化的RSP，微分的DSP，通用的GSP）。在结构和压电分析中，具有大变形的应力钢化功能。与其相似的耦合场单元有PLANE13、SOLID62和SOLID98。 INFIN9－二维无限边界 用于模拟一个二维无界问题的开放边界。具有两个节点，每个节点上带有磁向量势或温度自由度。所依附的单元类型可以为PLANE13和PLANE53磁单元，或PLANE55和PLANE77和PLANE35热单元。使用磁自由度（ＡＺ）时，分析可以是线性的也可以是非线性的，静态的或动态的。使用热自由度时，只能进行线性稳态分析。 PLANE13－二维耦合场实体 具有二维磁场、温度场、电场和结构场之间有限耦合的功能。由4个节点定义，每个节点可达到 4个自由度。具有非线性磁场功能，可用于模拟B－H曲线和永久磁铁去磁曲线。具有大变形和应力钢化功能。当用于纯结构分析时，具有大变形功能，相似的耦合场单元有SOLID5、SOLID98和SOLID62。LINK31－辐射线单元 用于模拟空间两点间辐射热流率的单轴单元。每个节点有一个自由度。可用于二维（平面或轴对称）或三维的、稳态的或瞬态的热分析问题。 允许形状因子和面积分别乘以温度的经验公式是有效的。发射率可与温度相关。如果包含热辐射单元的模型还需要进行结构分析，辐射单元应当被一个等效的或（空）结构单元所代替。 LINK32－二维传导杆 用于两节点间热传导的单轴单元。该单元每个节点只有一个温度自由度。可用于二维（平面或轴对称）稳态或瞬态的热分析问题。 如果包含热传导杆单元的模型还需进行结构分析，该单元可被一个等效的结构单元所代替。 LINK33－三维传导杆 用于节点间热传导的单轴单元。该单元每个节点只有一个温度自由度。可用于稳态或瞬态的热分析问题。 如果包含热传导杆单元的模型还需进行结构分析，该单元可被一个等效的结构单元所代替。
热应力分析
ABAQUS可以求解以下类型的传热问题： 1.非耦合传热分析：温度场不受应力应变场或电场的影响。应用ABAQUS/Standard可以求 解导热问题、强制对流、边界辐射和空腔辐射问题，其分析类型可以是瞬态或稳态的，也可以是线性或非线性的。 2.顺序耦合热应力分析：应力应变场受温度场的影响，但温度场不受应力应变场的影响。 此类问题用ABAQUS/Standard求解的步骤为：先求解温度场，然后以其作为已知条件，进行热应力分析，得到应力应变场。分析传热问题和热应力分析可以使用不一样的网格，abaqus会自动进行差值处理（此类问题称为热应力分析）。 3.完全耦合热应力分析：温度场和应力应变场之间有着强烈的相互作用。 4.绝热分析：在此类分析中，力学变形会产生热，而且整个过程中时间极短，不发生热扩 散。 5.热电耦合分析：用来求解电流产生的温度场。 7.1热应力分析中的主要问题 设定线胀系数、模型的初始温度场，并可以修改分析步中的温度场。 7.2带孔平板的热应力分析 学习： 在LOAD功能模块中，使用预定义场（predefined field）来定义温度场。 在此模块中可以直接指定温度场或读入分析结果文件中的温度场，可以指定并精确读入某个分析步中某个增量步的温度场 7.3法兰盘感应淬火的残余应力模拟 学习： 使用热应力来模拟残余应力；在LOAD功能模块中，为模型的各个区域定义不同的温度场 表面感应淬火：常用的热处理工艺，使用感应器对工件表面进行局部加热，然后迅速冷却，在工件内部产生残余压应力。它可以提高工件的弯曲疲劳抗力和扭转疲劳抗力，工件表面的
马氏体具有良好的耐磨性。 Abaqus可以完整的模拟淬火的全过程，即通过分析工件和感应器之间以及工件与冷却液之间的热场过程来确定工件的温度场，从而得到相应的塑性应变场和冷却后的残余应变场。 比较简单的模拟方法：先设定整个模型的初始温度场，在分析过程中令淬硬层区域的温度升高至某个温度值，其余区域的温度保持不变。经过几次试算，找到合适的淬硬层温度值，使得法兰盘内圆角处的表面压应力与实验结果吻合。施加工作载荷，保持上述温度场不变，就可以模拟在残余应力作用下的应力场。 优点：通用性强，可以模拟不同工艺所产生的残余应力场 缺点：精确度不高 改进方法：参淬硬层的不同区域设定不同的温度值
热应力效应
激光加热热分析 本文使用ANSYS的参数化设计语言APDL建立热分析模型，模拟分析区域为圆形和方形的温度分布和应力/应变分布。采用直接耦合法同时计算温度场和应力场。单元类型选择SOLID226，设定自由度为结构和温度自由度。分析条件如下表一所示。 表一分析条件 一、圆柱形玻璃模型 采用ANSYS内部的建模工具建立几何模型，由于光斑内部热流较为集中，将几何模型分为内部加热区和外部扩散区两个部分，三维图形如图1。 图1 三维几何模型 采用高斯热源模拟激光对工件的加热情况，在光斑内热流密度随半径变化较大，需要细化加热区内的网格，提高求解精度。为控制网格数量以及网格质量，整体采用映射方式划分网格，高度方向上让网格大小线性变化，靠近热源处网格较小。划分后的有限元网格模型如图2所示。其中节点个数为44671个，单元个数为42240个。
图2 有限元网格模型 采用函数边界条件加载指数变化的热流密度边界条件，图3为ANSYS的函数定义窗口。定义好函数之后，图形显示其变化曲线，检查热流密度变化情况，图4为热流密度随中心距离变化曲线。 图3 函数定义窗口 图4 热流密度随中心距离变化曲线
在定义完函数之后，在圆柱形上表面施加热流密度载荷，打开热流密度载荷显示符号。图5为圆柱体上表面热流密度等值线云图。 图5 热流密度等值线分布云图 选取圆柱表面节点，施加固定温度载荷20℃，如图6所示 图6 施加温度载荷 完成上述设置之后，进入求解器，进行如下求解选项设置： ●设定分析类型为稳态； ●设定载荷步控制选项，打开自动时间步； ●选择PCG求解器； ●求解输出设置，保存每个载荷步最后子步的结果；