文档库 最新最全的文档下载
当前位置:文档库 › 基于ANSYS的发动机气门热_应力耦合分析_王景

基于ANSYS的发动机气门热_应力耦合分析_王景

基于ANSYS的发动机气门热_应力耦合分析_王景
基于ANSYS的发动机气门热_应力耦合分析_王景

ANSYS热应力分析经典例题

ANSYS热应力分析例题 实例1圆简内部热应力分折: 有一无限长圆筒,其核截面结构如图13—1所示,简内壁温度为200℃,外壁温度为20℃,圆筒材料参数如表13.1所示,求圆筒内的温度场、应力场分布。 该问题属于轴对称问题。由于圆筒无限长,忽略圆筒端部的热损失。沿圆筒纵截面取宽度为10M的如图13—2所示的矩形截面作为几何模型。在求解过程中采用间接求解法和直接求解法两种方法进行求解。间接法是先选择热分析单元,对圆筒进行热分析,然后将热分析单元转化为相应的结构单元,对圆筒进行结构分析;直接法是采用热应力藕合单元,对圆筒进行热力藕合分析。 /filname,exercise1-jianjie /title,thermal stresses in a long /prep7 $Et,1,plane55 Keyopt,1,3,1 $Mp,kxx,1,70 Rectng,0.1,0.15,0,0.01 $Lsel,s,,,1,3,2 Lesize, all,,,20 $Lsel,s,,,2,4,2 Lesize,all,,,5 $Amesh,1 $Finish /solu $Antype,static Lsel,s,,,4 $Nsll,s,1 $d,all,temp,200 lsel,s,,,2 $nsll,s,1 $d,all,temp,20 allsel $outpr,basic,all solve $finish /post1 $Set,last /plopts,info,on Plnsol,temp $Finish /prep7 $Etchg,tts Keyopt,1,3,1 $Keyopt,1,6,1 Mp,ex,1,220e9 $Mp,alpx,,1,3e-6 $Mp,prxy,1,0.28 Lsel,s,,,4 $Nsll,s,1 $Cp,8,ux,all Lsel,s,,,2 $Nsll,s,1 $Cp,9,ux,all Allsel $Finish /solu $Antype,static D,all,uy,0 $Ldread,temp,,,,,,rth Allsel $Solve $Finish /post1 /title,radial stress contours Plnsol,s,x /title,axial stress contours Plnsol,s,y /title,circular stress contours Plnsol,s,z /title,equvialent stress contours Plnsol,s,eqv $finish

ansys热分析

第三章稳态热分析 3.1稳态传热的定义 ANSYS/Multiphysics,ANSYS/Mechanical,ANSYS/FLOTRAN和 ANSYS/Professional这些产品支持稳态热分析。稳态传热用于分析稳定的热载荷对系统或部件的影响。通常在进行瞬态热分析以前,进行稳态热分析用于确定初始温度分布。也可以在所有瞬态效应消失后,将稳态热分析作为瞬态热分析的最后一步进行分析。 稳态热分析可以计算确定由于不随时间变化的热载荷引起的温度、热梯度、热流率、热流密度等参数。这些热载荷包括: 对流 辐射 热流率 热流密度(单位面积热流) 热生成率(单位体积热流) 固定温度的边界条件 稳态热分析可用于材料属性固定不变的线性问题和材料性质随温度变化的非线性问题。事实上,大多数材料的热性能都随温度变化,因此在通常情况下,热分析都是非线性的。当然,如果在分析中考虑辐射,则分析也是非线性的。 3.2热分析的单元 ANSYS和ANSYS/Professional中大约有40种单元有助于进行稳态分析。有关单元的详细描述请参考《ANSYS Element Reference》,该手册以单元编号来讲述单元,第一个单元是LINK1。单元名采用大写,所有的单元都可用于稳态和瞬态热分析。其中SOLID70单元还具有补偿在恒定速度场下由于传质导致的热流的功能。这些热分析单元如下: 表3-1二维实体单元 单元维数形状及特点自由度 PLANE35 二维六节点三角形单元温度(每个节点) PLANE55 二维四节点四边形单元温度(每个节点) PLANE75 二维四节点谐单元温度(每个节点) PLANE77 二维八节点四边形单元温度(每个节点) PLANE38 二维八节点谐单元温度(每个节点)

ANSYS热分析指南与经典案例

第一章简介 一、热分析的目的 热分析用于计算一个系统或部件的温度分布及其它热物理参数,如热量的获取或损失、热梯度、热流密度(热通量〕等。 热分析在许多工程应用中扮演重要角色,如内燃机、涡轮机、换热器、管路系统、电子元件等。 二、ANSYS的热分析 ?在ANSYS/Multiphysics、ANSYS/Mechanical、ANSYS/Thermal、ANSYS/FLOTRAN、ANSYS/ED五种产品中包含热分析功能,其中 ANSYS/FLOTRAN不含相变热分析。 ?ANSYS热分析基于能量守恒原理的热平衡方程,用有限元法计算各节点的温度,并导出其它热物理参数。 ?ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外,还可以分析相变、有内热源、接触热阻等问题。 三、ANSYS 热分析分类 ?稳态传热:系统的温度场不随时间变化 ?瞬态传热:系统的温度场随时间明显变化 四、耦合分析 ?热-结构耦合 ?热-流体耦合 ?热-电耦合 ?热-磁耦合 ?热-电-磁-结构耦合等

第二章 基础知识 一、符号与单位 W/m 2-℃ 3 二、传热学经典理论回顾 热分析遵循热力学第一定律,即能量守恒定律: ● 对于一个封闭的系统(没有质量的流入或流出〕 PE KE U W Q ?+?+?=- 式中: Q —— 热量; W —— 作功; ?U ——系统内能; ?KE ——系统动能; ?PE ——系统势能; ● 对于大多数工程传热问题:0==PE KE ??; ● 通常考虑没有做功:0=W , 则:U Q ?=; ● 对于稳态热分析:0=?=U Q ,即流入系统的热量等于流出的热量; ● 对于瞬态热分析:dt dU q = ,即流入或流出的热传递速率q 等于系统内能的变化。 三、热传递的方式 1、热传导 热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热传导遵循付里叶定律:dx dT k q -='',式中''q 为热流

ANSYS稳态热分析的基本过程和实例

ANSYS稳态热分析的基本过程 ANSYS热分析可分为三个步骤: ?前处理:建模、材料和网格 ?分析求解:施加载荷计算 ?后处理:查看结果 1、建模 ①、确定jobname、title、unit; ②、进入PREP7前处理,定义单元类型,设定单元选项; ③、定义单元实常数; ④、定义材料热性能参数,对于稳态传热,一般只需定义导热系数,它可 以是恒定的,也可以随温度变化; ⑤、创建几何模型并划分网格,请参阅《ANSYS Modeling and Meshing Guide》。 2、施加载荷计算 ①、定义分析类型 ●如果进行新的热分析: Command: ANTYPE, STATIC, NEW GUI: Main menu>Solution>-Analysis Type->New Analysis>Steady-state ●如果继续上一次分析,比如增加边界条件等: Command: ANTYPE, STATIC, REST GUI: Main menu>Solution>Analysis Type->Restart ②、施加载荷 可以直接在实体模型或单元模型上施加五种载荷(边界条件) : a、恒定的温度 通常作为自由度约束施加于温度已知的边界上。 Command Family: D GUI:Main Menu>Solution>-Loads-Apply>-Thermal-Temperature b、热流率 热流率作为节点集中载荷,主要用于线单元模型中(通常线单元模型不能施加对流或热流密度载荷),如果输入的值为正,代表热流流入节点,即单元获取热量。如果温度与热流率同时施加在一节点上则ANSYS读取温度值进行计算。 注意:如果在实体单元的某一节点上施加热流率,则此节点周围的单元要

一个经典的ansys热分析实例(流程序)

/PREP7 /TITLE,Steady-state thermal analysis of pipe junction /UNITS,BIN ! 英制单位;Use U. S. Customary system of units (inches) ! /SHOW, ! Specify graphics driver for interactive run ET,1,90 ! Define 20-node, 3-D thermal solid element MP,DENS,1,.285 ! Density = .285 lbf/in^3 MPTEMP,,70,200,300,400,500 ! Create temperature table MPDATA,KXX,1,,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12 ! 指定与温度相对应的数据材料属性;导热系数;Define conductivity values MPDATA,C,1,,.113,.117,.119,.122,.125 ! Define specific heat values(比热) MPDATA,HF,2,,426/144,405/144,352/144,275/144,221/144 ! Define film coefficient;除144是单位问题,上面的除12也是单元问题 ! Define parameters for model generation RI1=1.3 ! Inside radius of cylindrical tank RO1=1.5 ! Outside radius Z1=2 ! Length RI2=.4 ! Inside radius of pipe RO2=.5 ! Outside pipe radius Z2=2 ! Pipe length CYLIND,RI1,RO1,,Z1,,90 ! 90 degree cylindrical volume for tank WPROTA,0,-90 ! 旋转当前工作的平面;从Y到Z旋转-90度;;Rotate working plane to pipe axis CYLIND,RI2,RO2,,Z2,-90 ! 角度选择在了第四象限;90 degree cylindrical volume for pipe WPSTYL,DEFA ! 重新安排工作平面的设置;另外WPSTYL,STAT to list the status of the working plane;;Return working plane to default setting BOPT,NUMB,OFF ! 关掉布尔操作的数字警告信息;Turn off Boolean numbering warning VOVLAP,1,2 ! 交迭体;Overlap the two cylinders /PNUM,VOLU,1 ! 体编号打开;Turn volume numbers on /VIEW,,-3,-1,1

ansys热分析例题

问题描述:一个30公斤重、温度为70℃的铜块,以及一个20公斤重、温度为80℃的铁块,突然放入温度为20℃、盛满了300升水的、完全绝热的水箱中,如图所示。过了一个小时,求铜块与铁块的最高温度(假设忽略水的流动)。 材料热物理性能如下:热性能单位制 铜铁水 导热系数W/m℃ 383 37 密度Kg/m 8889 7833 996 比热J/kg℃ 390 448 4185 菜单操作过程: 一、设置分析标题 1、选择“Utility Menu>File>Change Jobname”,输入文件名Transient1。 2、选择“Utility Menu>File>Change Title”输入Thermal Transient Exercise 1。 二、定义单元类型 1、选择“Main Menu>Preprocessor”,进入前处理。 2、选择“Main Menu>Preprocesor>Element Type>Add/Edit/Delete”。选择热平面单元plane77。 三、定义材料属性 1、选择“Main Menu>Preprocessor>Material Props>Material Models”,在弹出的材料定义窗口中顺序双击Thermal选项。 2、点击Conductivity,Isotropic,在KXX框中输入383;点击Density,在DENS框中输入8898;点击Specific Heat,在C框中输入390。 3、在材料定义窗口中选择Material>New Model,定义第二种材料。 4、点击Conductivity,Isotropic,在KXX框中输入70;点击Density,在DENS框中输入7833;点击Specific Heat,在C框中输入448。 5、在材料定义窗口中选择Material>New Model,定义第三种材料。 6、点击Conductivity,Isotropic,在KXX框中输入.61;点击Density,在DENS框中输入996;点击Specific Heat,在C框中输入4185。 四、创建几何模型 1、选择“Main Menu>Preprocessor>-Modeling->Create>-Areas->Retangle>By Dimensions”,输入X1=0, Y1=0, X2=, Y2=, 点击Apply;输入X1=, Y1=, X2= ,Y2=, 点击Apply;输入X1= Y1=, X2= Y2=+, 选择OK。 2、选择“Main Menu>Preprocessor>-Modeling->Operate>Booleans>Overlap”,选择Pick All。 3、选择“Utility Menu>Plotctrls>Numbering>Areas, on”。 4、选择“Utility Menu>Plot>Areas”。 五、划分网格 1、选择“Main Menu>Preprocessor>-Attributes->Define->All Areas”,选择材料1。 2、选择“Main Menu>Preprocessor>Meshing->Size Cntrls->-Manualsize->-Global->Size”,输入单元大小。 3、选择“Main Menu>Preprocessor>Meshing->Mesh->-Areas->Mapped>3 or 4 sided”,选择铜块。 4、选择“Main Menu>Preprocessor>-Attributes->Define->All Areas”,选择材料2。 5、选择“Main Menu>Preprocessor>Meshing->Mesh->-Areas->Mapped>3 or 4 sided”,选

热结构耦合

第21章热-结构耦合分析 热-结构耦合问题是结构分析中通常遇到的一类耦合分析问题。由于结构温度场的分布不均会引起结构的热应力,或者结构部件在高温环境中工作,材料受到温度的影响会发生性能的改变,这些都是进行结构分析时需要考虑的因素。为此需要先进行相应的热分析,然后在进行结构分析。热分析用于计算一个系统或部件的温度分布及其它热物理参数,如热量的获取或损失、热梯度、热流密度(热通量)等。本章主要介绍在ANSYS中进行稳态、瞬态热分析的基本过程,并讲解如何完整的进行热-结构耦合分析。 21.1 热-结构耦合分析简介 热-结构耦合分析是指求解温度场对结构中应力、应变和位移等物理量影响的分析类型。对于热-结构耦合分析,在ANSYS中通常采用顺序耦合分析方法,即先进行热分析求得结构的温度场,然后再进行结构分析。且将前面得到的温度场作为体载荷加到结构中,求解结构的应力分布。为此,首先需要了解热分析的基本知识,然后再学习耦合分析方法。 21.1.1 热分析基本知识 ANSYS热分析基于能量守恒原理的热平衡方程,用有限元法计算各节点的温度,并导出其它热物理参数。ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。此外,还可以分析相变、有内热源、接触热阻等问题。 热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。热对流是指固体的表面和与它周围接触的流体之间,由于温差的存在引起的热量的交换。热辐射指物体发射电磁能,并被其它物体吸收转变为热的热量交换过程。 如果系统的净热流率为0,即流入系统的热量加上系统自身产生的热量等于流出系统的热量:q流入+q生成-q流出=0,则系统处于热稳态。在稳态热分析中任一节点的温度不随时间变化。 瞬态传热过程是指一个系统的加热或冷却过程。在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。

ANSYS热应力分析 精选实例x

ANSYS热应力分析实例 当一个结构加热或冷却时,会发生膨胀或收缩。如果结构各部分之间膨胀收缩程度不同,和结构的膨胀、收缩受到限制,就会产生热应力。 热应力分析的分类 ANSYS提供三种进行热应力分析的方法: 在结构应力分析中直接定义节点的温度。如果所以节点的温度已知,则可以通过命令直接定义节点温度。节点温度在应力分析中作为体载荷,而不是节点自由度间接法。首先进行热分析,然后将求得的节点温度作为体载荷施加在结构应力分析中。 直接法。使用具有温度和位移自由度的耦合单元,同时得到热分析和结构应 力分析的结果。 如果节点温度已知,适合第一种方法。但节点温度一般是不知道的。对于大多数问题,推荐使用第二种方法一间接法。因为这种方法叮以使用所有热分析的功能和结构分析的功能。如果热分析是瞬态的,只需要找出温度梯度最大的时间点,并将此时间点的节点温度作为荷载施加到结构应力分析中去。如果热和结构的耦合是双向的,即热分析影响结构应力分析,同时结构变形又会影响热分析(如大变形、接触等),则可以使用第三种直接法一使用耦合单元。此外只有第三种方法可以考虑其他分析领域(电磁、流体等)对热和结构的影响。 间接法进行热应力分析的步骤 首先进行热分析。可以使用热分析的所有功能,包括传导、对流、辐射和表 面效应单元等,进行稳态或瞬态热分析。但要注意划分单元时要充分考虑结构分

析的要求。例如,在有可能有应力集中的地方的网格要密一些。如果进行瞬态分析,在后处理中要找出热梯度最大的时间点或载荷步。 重新进入前处理,将热单元转换为相应的结构单元,表7-1是热单元与结构

单元的对应表。可以使用菜单进行转换:

Ansys模拟水结冰的热分析过程

Ansys模拟水结冰的热分析过程 一、问题描述: 对一茶杯水的结冰过程进行分析,水和茶杯的初始温度为0℃,环境温度为-10℃,杯子侧面和顶面的对流换热系数为12.5W/m^2·℃,杯子放在桌面上,假设桌面可以对杯子底面提供-10℃的温度载荷。计算3000s之后的温度分布。 模型如下: 茶杯底面外径54.41mm,内径50mm,高度85mm,顶面内径60mm,抽壳厚度为5mm(内部水的高度80mm)。 分析采用SI单位制,水的材料属性如下: 导热率:0.6 密度:1000 比热容:4200

焓值: 温度℃-10 -1 0 10 焓J/m^3 0 37.8e6 79.8e6 121.8e6 茶杯采用铁的材料属性: 导热率:70 密度:7833 比热容:448 二、问题分析: 本例采用70热单元进行分析,由于对称性,采用1/4模型进行建模分析。由于包含相变分析,因此水的焓值是必要的。假设温度0℃的水结成0℃的冰需要放出42000J/kg·℃的热量,通过定义焓值来实现。假设温度区间长度为1℃,因此温度低于-1℃,表示水已结成冰。本例通过apdl进行分析,方便输入及调试。 三、分析步骤: 1、定义工程名及标题 fini /cle !清除数据库 /filname.shuijiebing,1 !此处设置工作名 /title,lovz !此处设置标题 *afun deg !定义角度为度 2、进入前处理,定义单元及材料属性 /prep7 !进入前处理模块

et,1,70 !定义70单元 mp,kxx,1,0.6 !设置材料属性 mp,c,1,4200 mp,dens,1,1000 mptemp,1,-10,-1,0,10 mpdata,enth,1,1,0,37.8e6,79.8e6,121.8e6 !焓值定义mp,kxx,2,70 mp,dens,2,7833 mp,c,2,448 这里定义1号材料为水,2号材料为茶杯 3、定义参数 r1=50e-3 r2=60e-3 r3=54.41e-3 r4=65e-3 h1=80e-3 h2=85e-3 4、建模 wprot,,-90 !旋转工作平面 /pnum,volu,1 !打开体积显示 /view,1,1,1,1 !Iso视角 cone,r1,r2,0,h2,0,90 !建立水的1/4圆台模型

热力耦合分析单元简介

热力耦合分析单元简介! SOLID5-三维耦合场实体 具有三维磁场、温度场、电场、压电场和结构场之间有限耦合的功能。本单元由8个节点定义,每个节点有6个自由度。在静态磁场分析中,可以使用标量势公式(对于简化的RSP,微分的DSP,通用的GSP)。在结构和压电分析中,具有大变形的应力钢化功能。与其相似的耦合场单元有PLANE13、SOLID62和SOLID98。 INFIN9-二维无限边界 用于模拟一个二维无界问题的开放边界。具有两个节点,每个节点上带有磁向量势或温度自由度。所依附的单元类型可以为PLANE13和PLANE53磁单元,或PLANE55和PLANE77和PLANE35热单元。使用磁自由度(AZ)时,分析可以是线性的也可以是非线性的,静态的或动态的。使用热自由度时,只能进行线性稳态分析。 PLANE13-二维耦合场实体 具有二维磁场、温度场、电场和结构场之间有限耦合的功能。由4个节点定义,每个节点可达到4个自由度。具有非线性磁场功能,可用于模拟B-H曲线和永久磁铁去磁曲线。具有大变形和应力钢化功能。当用于纯结构分析时,具有大变形功能,相似的耦合场单元有SOLID5、SOLID98和SOLID62。 LINK31-辐射线单元 用于模拟空间两点间辐射热流率的单轴单元。每个节点有一个自由度。可用于二维(平面或轴对称)或三维的、稳态的或瞬态的热分析问题。 允许形状因子和面积分别乘以温度的经验公式是有效的。发射率可与温度相关。如果包含热辐射单元的模型还需要进行结构分析,辐射单元应当被一个等效的或(空)结构单元所代替。 LINK32-二维传导杆 用于两节点间热传导的单轴单元。该单元每个节点只有一个温度自由度。可用于二维(平面或轴对称)稳态或瞬态的热分析问题。 如果包含热传导杆单元的模型还需进行结构分析,该单元可被一个等效的结构单元所代替。 LINK33-三维传导杆 用于节点间热传导的单轴单元。该单元每个节点只有一个温度自由度。可用于稳态或瞬态的热分析问题。 如果包含热传导杆单元的模型还需进行结构分析,该单元可被一个等效的结构单元所代替。 LINK34-对流线单元 用于模拟节点间热对流的单轴单元。该单元每个节点只有一个温度自由度。热对流杆单元可用于二维(平面或轴对称)或三维、稳态或瞬态的热分析问题。 如果包含热对流单元的模型还需要进行结构分析,热对流单元可被一个等效(或空)的结构单元所代替。单元的对流换热系数可分为非线性,即对流换热系数是温度或时间的函数。

ANSYS传热分析实例汇总

实例1: 某一潜水艇可以简化为一圆筒,它由三层组成,最外面一层为不锈钢,中间为玻纤隔热层,最里面为铝层,筒内为空气,筒外为海水,求内外壁面温度及温度分布。 几何参数: 筒外径30 feet 总壁厚 2 inch 不锈钢层壁厚0、75 inch 玻纤层壁厚 1 inch 铝层壁厚0、25 inch 筒长200 feet 导热系数不锈钢8、27 BTU/hr、ft、o F 玻纤0、028 BTU/hr、ft、o F 铝117、4 BTU/hr、ft、o F 边界条件空气温度70 o F 海水温度44、5 o F 空气对流系数2、5 BTU/hr、ft2、o F 海水对流系数80 BTU/hr、ft2、o F 沿垂直于圆筒轴线作横截面,得到一圆环,取其中1度进行分析,如图示。 以下分别列出log文件与菜单文件。 /, Steady1 /title, Steady-state thermal analysis of submarine /units, BFT Ro=15 !外径(ft) Rss=15-(0、75/12) !不锈钢层内径ft) Rins=15-(1、75/12) !玻璃纤维层内径(ft)

Ral=15-(2/12) !铝层内径(ft) Tair=70 !潜水艇内空气温度 Tsea=44、5 !海水温度 Kss=8、27 !不锈钢的导热系数(BTU/hr、ft、oF) Kins=0、028 !玻璃纤维的导热系数(BTU/hr、ft、oF) Kal=117、4 !铝的导热系数(BTU/hr、ft、oF) Hair=2、5 !空气的对流系数(BTU/hr、ft2、oF) Hsea=80 !海水的对流系数(BTU/hr、ft2、oF) /prep7 et,1,plane55 !定义二维热单元 mp,kxx,1,Kss !设定不锈钢的导热系数 mp,kxx,2,Kins !设定玻璃纤维的导热系数 mp,kxx,3,Kal !设定铝的导热系数 pcirc,Ro,Rss,-0、5,0、5 !创建几何模型 pcirc,Rss,Rins,-0、5,0、5 pcirc,Rins,Ral,-0、5,0、5 aglue,all numcmp,area lesize,1,,,16 !设定划分网格密度 lesize,4,,,4 lesize,14,,,5 lesize,16,,,2 eshape,2 !设定为映射网格划分 mat,1 amesh,1 mat,2 amesh,2 mat,3 amesh,3 /SOLU SFL,11,CONV,HAIR,,TAIR !施加空气对流边界 SFL,1,CONV,HSEA,,TSEA !施加海水对流边界 SOLVE /POST1 PLNSOL !输出温度彩色云图

ANSYS_热分析(两个实例)有限元热分析上机指导书

第四讲 热分析上机指导书 CAD/CAM 实验室,USTC 实验要求: 1、通过对冷却栅管的热分析练习,熟悉用ANSYS 进行稳态热分析的基本过程,熟悉用直接耦合法、间接耦合法进行热应力分析的基本过程。 2、通过对铜块和铁块的水冷分析,熟悉用ANSYS 进行瞬态热分析的基本过程。 内容1:冷却栅管问题 问题描述:本实例确定一个冷却栅管(图a )的温度场分布及位移和应力分布。一个轴对称的冷却栅结构管内为热流体,管外流体为空气。冷却栅材料为不锈钢,特性如下: 导热系数: W/m ℃ 弹性模量:×109 MPa 热膨胀系数:×10-5 /℃ 泊松比: 边界条件: (1)管内:压力: MPa 流体温度:250 ℃ 对流系数 W/m 2℃ (2)管外:空气温度39℃ 对流系数: W/m 2℃ 假定冷却栅管无限长,根据冷却 栅结构的对称性特点可以构造出的有限元模型如图b 。其上下边界承受边界约束,管内部承受均布压力。 练习1-1:冷却栅管的稳态热分析 步骤: 定义工作文件名及工作标题 定义工作文件名:GUI: Utility Menu> File> Change Jobname ,在弹出的【Change Jobname 】对话框中输入文件名Pipe_Thermal ,单击OK 按钮。 定义工作标题:GUI: Utility Menu> File> Change Title ,在弹出的【Change Title 】对话框中2D Axisymmetrical Pipe Thermal Analysis ,单击OK 按钮。 关闭坐标符号的显示:GUI: Utility Menu> PlotCtrls> Window Control> Window Options ,在弹出的【Window Options 】对话框的Location of triad 下拉列表框中选择No Shown 选项,单击OK 按钮。 定义单元类型及材料属性 定义单元类型:GUI: Main Menu> Preprocessor> Element Type> Add/Edit>Delete 命令,弹出【Element Types 】对话框,单击Add 按钮,弹出【Library Type 】对话框,选择Thermal Solid Quad 8node 77选项,单击OK 按钮。 设置单元选项:单击【Element Type 】对话框的Options 按钮,弹出【Plane77 element type options 】对话框,在Element behavior 下拉列框中选择Axisymmetrical 选项,单击OK 按钮,单击Close 按钮。 设置材料属性:GUI: Main Menu> Preprocessor> Material Props> Material Models ,弹出【Define () ()

热应力分析

ABAQUS可以求解以下类型的传热问题: 1.非耦合传热分析:温度场不受应力应变场或电场的影响。应用ABAQUS/Standard可以求 解导热问题、强制对流、边界辐射和空腔辐射问题,其分析类型可以是瞬态或稳态的,也可以是线性或非线性的。 2.顺序耦合热应力分析:应力应变场受温度场的影响,但温度场不受应力应变场的影响。 此类问题用ABAQUS/Standard求解的步骤为:先求解温度场,然后以其作为已知条件,进行热应力分析,得到应力应变场。分析传热问题和热应力分析可以使用不一样的网格,abaqus会自动进行差值处理(此类问题称为热应力分析)。 3.完全耦合热应力分析:温度场和应力应变场之间有着强烈的相互作用。 4.绝热分析:在此类分析中,力学变形会产生热,而且整个过程中时间极短,不发生热扩 散。 5.热电耦合分析:用来求解电流产生的温度场。 7.1热应力分析中的主要问题 设定线胀系数、模型的初始温度场,并可以修改分析步中的温度场。 7.2带孔平板的热应力分析 学习: 在LOAD功能模块中,使用预定义场(predefined field)来定义温度场。 在此模块中可以直接指定温度场或读入分析结果文件中的温度场,可以指定并精确读入某个分析步中某个增量步的温度场 7.3法兰盘感应淬火的残余应力模拟 学习: 使用热应力来模拟残余应力;在LOAD功能模块中,为模型的各个区域定义不同的温度场 表面感应淬火:常用的热处理工艺,使用感应器对工件表面进行局部加热,然后迅速冷却,在工件内部产生残余压应力。它可以提高工件的弯曲疲劳抗力和扭转疲劳抗力,工件表面的

马氏体具有良好的耐磨性。 Abaqus可以完整的模拟淬火的全过程,即通过分析工件和感应器之间以及工件与冷却液之间的热场过程来确定工件的温度场,从而得到相应的塑性应变场和冷却后的残余应变场。 比较简单的模拟方法:先设定整个模型的初始温度场,在分析过程中令淬硬层区域的温度升高至某个温度值,其余区域的温度保持不变。经过几次试算,找到合适的淬硬层温度值,使得法兰盘内圆角处的表面压应力与实验结果吻合。施加工作载荷,保持上述温度场不变,就可以模拟在残余应力作用下的应力场。 优点:通用性强,可以模拟不同工艺所产生的残余应力场 缺点:精确度不高 改进方法:参淬硬层的不同区域设定不同的温度值

热力耦合分析单元简介

共享:热力耦合分析单元简介! 挑选了部分常用的,希望能方便大家的使用,其中自己翻译了一部分,不准确之处还望见谅,大家还可以继续补充哦!: SOLID5-三维耦合场实体    具有三维磁场、温度场、电场、压电场和结构场之间有限耦合的功能。本单元由8个节点定义,每个节点有6个自由度。在静态磁场分析中,可以使用标量势公式(对于简化的RSP,微分的DSP,通用的GSP)。在结构和压电分析中,具有大变形的应力钢化功能。与其相似的耦合场单元有PLANE13、SOLID62和SOLID98。 INFIN9-二维无限边界    用于模拟一个二维无界问题的开放边界。具有两个节点,每个节点上带有磁向量势或温度自由度。所依附的单元类型可以为PLANE13和PLANE53磁单元,或PLANE55和PLANE77和PLANE35热单元。使用磁自由度(AZ)时,分析可以是线性的也可以是非线性的,静态的或动态的。使用热自由度时,只能进行线性稳态分析。 PLANE13-二维耦合场实体    具有二维磁场、温度场、电场和结构场之间有限耦合的功能。由4个节点定义,每个节点可达到 4个自由度。具有非线性磁场功能,可用于模拟B-H曲线和永久磁铁去磁曲线。具有大变形和应力钢化功能。当用于纯结构分析时,具有大变形功能,相似的耦合场单元有SOLID5、SOLID98和SOLID62。LINK31-辐射线单元   用于模拟空间两点间辐射热流率的单轴单元。每个节点有一个自由度。可用于二维(平面或轴对称)或三维的、稳态的或瞬态的热分析问题。   允许形状因子和面积分别乘以温度的经验公式是有效的。发射率可与温度相关。如果包含热辐射单元的模型还需要进行结构分析,辐射单元应当被一个等效的或(空)结构单元所代替。 LINK32-二维传导杆   用于两节点间热传导的单轴单元。该单元每个节点只有一个温度自由度。可用于二维(平面或轴对称)稳态或瞬态的热分析问题。   如果包含热传导杆单元的模型还需进行结构分析,该单元可被一个等效的结构单元所代替。 LINK33-三维传导杆   用于节点间热传导的单轴单元。该单元每个节点只有一个温度自由度。可用于稳态或瞬态的热分析问题。   如果包含热传导杆单元的模型还需进行结构分析,该单元可被一个等效的结构单元所代替。

ANSYS热应力分析经典例题

ANSYS热应力分析例题 实例1——圆简内部热应力分折:有一无限长圆筒,其核截面结构如图13—1所示,简内壁温度为200℃,外壁温度为20℃,圆筒材料参数如表13.1所示,求圆筒内的温度场、应力场分布。 该问题属于轴对称问题。由于圆筒无限长,忽略圆筒端部的热损失。沿圆筒纵截面取宽度为10M的如图1 3—2所示的矩形截面作为几何模型。在求解过程中采用间接求解法和直接求解法两种方法进行求解。间接法是先选择热分析单元,对圆筒进行热分析,然后将热分析单元转化为相应的结构单元,对圆筒进行结构分析;直接法是采用热应力藕合单元,对圆筒进行热力藕合分析。 fini clear /filname,exercise1-jianjie /title,thermal stresses in a long /prep7 $Et,1,plane55 Keyopt,1,3,1 $Mp,kxx,1,70 Rectng,0.1,0.15,0,0.01 $Lsel,s,,,1,3,2 Lesize, all,,,20 $Lsel,s,,,2,4,2 Lesize,all,,,5 $Amesh,1 $Finish /solu $Antype,static Lsel,s,,,4 $Nsll,s,1 $d,all,temp,200 lsel,s,,,2 $nsll,s,1 $d,all,temp,20 allsel $outpr,basic,all solve $finish /post1 $Set,last /plopts,info,on Plnsol,temp $Finish /prep7 $Etchg,tts Keyopt,1,3,1 $Keyopt,1,6,1 Mp,ex,1,220e9 $Mp,alpx,,1,3e-6 $Mp,prxy,1,0.28 Lsel,s,,,4 $Nsll,s,1 $Cp,8,ux,all Lsel,s,,,2 $Nsll,s,1 $Cp,9,ux,all Allsel $Finish /solu $Antype,static D,all,uy,0 $Ldread,temp,,,,,,rth Allsel $Solve $Finish /post1 /title,radial stress contours Plnsol,s,x /title,axial stress contours Plnsol,s,y /title,circular stress contours Plnsol,s,z /title,equvialent stress contours

ANSYS传热分析实例汇总

实例1: 某一潜水艇可以简化为一圆筒,它由三层组成,最外面一层为不锈钢,中间为玻纤隔热层,最里面为铝层,筒内为空气,筒外为海水,求内外壁面温度及温度分布。 几何参数:筒外径30 feet 总壁厚 2 inch 不锈钢层壁厚0.75 i nch 玻纤层壁厚 1 inch 铝层壁厚0.25 i nch 筒长200 feet 导热系数不锈钢8.27 B TU/hr.ft.o F 玻纤0.028 BTU/hr.ft.o F 铝117.4 BTU/hr.ft.o F 边界条件空气温度70 o F 海水温度44.5 o F 空气对流系数 2.5 BTU/hr.ft2.o F 海水对流系数80 BTU/hr.ft2.o F 沿垂直于圆筒轴线作横截面,得到一圆环,取其中1度进行分析,如图示。 以下分别列出log文件和菜单文件。 /filename, Steady1 /title, Steady-state thermal analysis of submarine /units, BFT Ro=15 !外径(ft)

Rss=15-(0.75/12) !不锈钢层内径ft) Rins=15-(1.75/12) !玻璃纤维层内径(ft) Ral=15-(2/12) !铝层内径(ft) Tair=70 !潜水艇内空气温度 Tsea=44.5 !海水温度 Kss=8.27 !不锈钢的导热系数(BTU/hr.ft.oF) Kins=0.028 !玻璃纤维的导热系数(BTU/hr.ft.oF) Kal=117.4 !铝的导热系数(BTU/hr.ft.oF) Hair=2.5 !空气的对流系数(BTU/hr.ft2.oF) Hsea=80 !海水的对流系数(BTU/hr.ft2.oF) /prep7 et,1,plane55 !定义二维热单元 mp,kxx,1,Kss !设定不锈钢的导热系数 mp,kxx,2,Kins !设定玻璃纤维的导热系数 mp,kxx,3,Kal !设定铝的导热系数 pcirc,Ro,Rss,-0.5,0.5 !创建几何模型 pcirc,Rss,Rins,-0.5,0.5 pcirc,Rins,Ral,-0.5,0.5 aglue,all numcmp,area lesize,1,,,16 !设定划分网格密度 lesize,4,,,4 lesize,14,,,5 lesize,16,,,2 eshape,2 !设定为映射网格划分 mat,1 amesh,1 mat,2 amesh,2 mat,3 amesh,3 /SOLU SFL,11,CONV,HAIR,,TAIR !施加空气对流边界

第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的主菜单设置为与热分析相对应的菜单选项。单击按钮,完成

稳态热分析案例(ANSYS-15.0版)

本例题的主要部分为一个圆筒形罐,其上沿径向有一材料一样的接管(如图????所所示),罐内流动着450°F(232°C)的高温流体,接管内流动着100°F(38 °C)的低温流体,两个流体区域由薄壁管隔离。罐的对流换热系数为250Btu/hr-ft2-o F(1420watts/m2-°K),接管的对流换热系数随管壁温度而变,它的热物理性能如表???所示。要求计算罐与接管的温度分布。 表???? 6.5.1 预处理 Step 1: 确定分析标题 起动ANSYS后,开始一个分析,需要输入一个标题,按下面方法进行操作: 1.选择Utility Menu> File> Change Title,弹出相应对话框 2.输入 Steady-state thermal analysis of pipe junction。 3.点击OK。 Step 2: 设置分析单位系统 You need to specify units of measurement for the analysis. For this pipe junction example, measurements use the U. S. Customary system of units (based on inches).

To specify this, type the command /UNITS,BIN in the ANSYS Input window and press ENTER.在分析之前,需要为分析系统设定单位系统, Step 3: Define the Element Type The example analysis uses a thermal solid element. To define it, do the following: 1.Choose Main Menu> Preprocessor> Element Type> Add/Edit/Delete. The Element Types dialog box appears. 2.Click on Add. The Library of Element Types dialog box appears. 3.In the list on the left, scroll down and pick (highlight) "Thermal Solid." In the list on the right, pick "Brick20node 90." 4.Click on OK. 5.Click on Close to close the Element Types dialog box. Step 4: Define Material Properties To define material properties for the analysis, perform these steps: 1.Choose Main Menu> Preprocessor> Material Props> Material Models. The Define Material Model Behavior dialog box appears. 2.In the Material Models Available window, double-click on the following options: Thermal, Density. A dialog box appears. 3.Enter .285 for DENS (Density), and click on OK. Material Model Number 1 appears in the Material Models Defined window on the left. 4.In the Material Models Available window, double-click on the following options: Conductivity, Isotropic. A dialog box appears. 5.Click on the Add Temperature button four times. Four columns are added. 6.In the T1 through T5 fields, enter the following temperature values: 70, 200, 300, 400, and 500. Select the row of temperatures by dragging the cursor across the text fields. Then copy the temperatures by pressing Ctrl-c. 7.In the KXX (Thermal Conductivity) fields, enter the following values, in order, for each of the temperatures, then click on OK. Note that to keep the units consistent, each of the given values of KXX must be divided by 12. You can just input the fractions and have ANSYS perform the calculations. 8.35/12 8.90/12 9.35/12 9.80/12 10.23/12 8.In the Material Models Available window, double-click on Specific Heat. A dialog box appears. 9.Click on the Add Temperature button four times. Four columns are added.

相关文档