《ABAQUS热传导与热应力分析课稿课件》由会员分享,可在线阅读,更多相关《ABAQUS热传导与热应力分析课稿课件(43页珍藏版)》请在金锄头文库上搜索。
1、ABAQUSABAQUS专题教程热传导和热应力分析1ABAQUS热传导与热应力分析课稿第一讲:固体热传导介绍概述介绍分析过程材料热性质ABAQUS/Standard中的热传导单元库边界条件和载荷稳态分析瞬态分析非线性分析2ABAQUS热传导与热应力分析课稿介绍-ABAQUS主要是用来进行应力分析的软件-但ABAQUS也有一个重要的特性:就是可以求解规模大的、复杂的和多组件模型的热传导问题。热传导求解能力是从求解热应力问题中发展出来的3ABAQUS热传导与热应力分析课稿ABAQUS中的热传导特性-稳态响应-瞬态响应,包括自适应时间步长-全套热传导边界条件-材料属性(和载荷)可以是温度相关-热“接
2、触”允许在“接触表面”有热流动-可以方便的将温度场导入热应力分析中-特性潜热项(由相变产生)强制对流应力-热传导耦合分析功能热传导壳单元(沿厚度方向温度梯度)空腔辐射(加热炉升温)功能介绍4ABAQUS热传导与热应力分析课稿介绍ABAQUS不能做什么ABAQUS不是专业热传导分析软件无流体分析无自由对流无浮力驱使流动对热冲击问题无自适应网格划分无逆传热分析5ABAQUS热传导与热应力分析课稿介绍力平衡与能量守恒之间的类比-在应力分析中,ABAQUS求解力平衡方程:Mu=PI-在热传导分析中,ABAQUS求解能率守恒方程并确定温度的分布。密度 比热 温度变化率外部热量内部热量6ABAQUS热传导
3、与热应力分析课稿介绍热传导分析中的基本物理量-温度Temperature单位-热能Heatenergy单位J-热率Heatratepower单位J/torW-热流量Heatflux=Powerperunitarea单位J/t/L2-热传导率k,衡量物质中热量流动的能力单位J/T/L/:热流量正比于热传导率和温度梯度:TaTbALQ7ABAQUS热传导与热应力分析课稿介绍-比热,衡量物质储存热的能力单位:J/M/时间增量温度增量比热-一维热传导公式热扩散率8ABAQUS热传导与热应力分析课稿介绍-类比StressHeat9ABAQUS热传导与热应力分析课稿分析过程在ABAQUS/Standard
4、中,热传导分析的执行是通过将几何体离散成扩散热传导单元,并且使用*HEATTRANSFER过程选项*HEATTRANSFER瞬态分析(默认)*HEATTRANSFER,STEADYSTATE稳态分析在ABAQUS/Explicit中,没有单纯的热传导分析选项,然而可以进行全耦合的热-应力分析。这个功能通过设定适当的边界条件,可以模拟纯热传导工程;除空腔辐射和利用用户子程序定义的不均匀热载荷之外,其他在ABAQUS/Standard中可以使用的热属性,都可以用在Explicit中。10ABAQUS热传导与热应力分析课稿材料热性质定义材料的热性质在inp中的*MATERIAL关键字定义*MATER
5、IAL,NAME=MATERIAL-1*CONDUCTIVITY1.0*DENSITY1.0*SPECIFICHEAT1.0热传导率:*CONDUCTIVITY,可以定义各向同性(默认)或各向异性(正交或完全)用TYPE参数:*CONDUCTIVITY,TYPE=ISO|ORTHO|ANISO-热传导率可以是温度的函数,这样就成了一个非线性问题。-热传导率也可以是任意数量预设的场变量的函数-预设场变量相关的材料性质不会涉及非线性,ABAQUS使用简单的插值方法确定材料性质。例如:*CONDUCTIVITY,DEPENDENCIES=163.0,20,16070.5,200,200*INITIA
6、LCONDITIOINS,TYPE=FIELD,VAR=1NALL,160*STEP*FIELD,VARIABLE=1,AMPLITUDE=TIMEVARNALL,180*ENDSTEP材料热性质定义温度 场变量设置包括的预设场变量数量比热:*SPECIFICHEAT,-比热可以定义为随温度与场变量变化-大多数材料的比热随温度平稳变化密度:*DENSITY,-密度可以定义为随温度与场变量变化11ABAQUS热传导与热应力分析课稿热传导单元定义连续单元:ABAQUS中连续扩散热传导单元库包括:一阶(线性)插值单元二阶(抛物线)单元用于一维,二维,轴对称和三维应用单元命名规则:DC3D20扩散di
7、ffusion连续体continuum节点数几何,3D单元-这些单元节点的基本变量(自由度)是温度标量q ABAQUS中用自由度11表示温度。节点温度输出变量为NT11.点单元热容单元HEATCAP模拟在一点的集中热容热容可以是温度或场变量的函数该单元可以在ABAQUS/Explicit中使用12ABAQUS热传导与热应力分析课稿壳单元一阶和二阶插值用于轴对称单元(DSAX1,DSAX2)和三维(DS3,DS4,DS6,DS8)应用的壳单元包含有单元库中。壳单元用于模拟承受热载荷的薄壁结构如:压力容器,管道系统和金属片元件等。热传导单元定义-壳单元表面下方的温度自由度为11(输出变量为NT11
8、)-在正表面的温度自由度为10+n,n为壳截面上使用截面点的数量-在单层(均质)壳中,截面点在厚度上均匀分布,默认为5个点-每层壳必须是奇数个截面点,这是由ABAQUS/standard在厚度上使用分段抛物线型插值方法决定的。nNT11NT12NT13-单元在每个壳节点的厚度方向的多个点上提供了温度自由度,这样温度不仅随着壳的参考平面变化,也随厚度方向变化。13ABAQUS热传导与热应力分析课稿热传导单元定义复合材料壳单元材料1材料1材料1材料1t1t2t3t4多层复合材料热壳可以被构建每一层可以是不同厚度,不同主方向的不同材料组成材料特性在*SHELLSECTION中定义:*SHELLSEC
9、TION,COMPOSITELAYER1的厚度,温度自由度数量(截面点数),材料名,材料方向参考的orientation名称LAYER2的厚度,温度自由度数量(截面点数),材料名,材料方向参考的orientation名称LAYER3的厚度,温度自由度数量(截面点数),材料名,材料方向参考的orientation名称多层复合材料热壳的默认截面点数量为3所有层的单层截面点数量必须相等14ABAQUS热传导与热应力分析课稿边界条件与载荷边界条件应力分析中,每个自由度都有一对共轭变量:位移-作用或反作用力默认情况下位移是未知的,力是已知的。热传导分析中,这对共轭变量是温度-热率(单位时间的能量流)默认
10、情况下温度是未知的,热率是已知的-已知的热率=0,相当于绝热边界条件;-没有外部的能量流进或流出节点。ABAQUS中的几种热边界条件和热载荷1.在某些节点上预设温度,*BOUNDARY,自由度112.在某些点上或者某些表面上或者体积内预设热率q*CFLUX,*DFLUX,*DSFLUX3.在某些点上或者某些表面上的边界层(薄膜)条件*CFILM,*FILM和*SFILM4.在某些点上或者某些表面上的辐射条件*CRADIATE,*RADIATE,和*SRADIATE5.自然边界条件(默认)15ABAQUS热传导与热应力分析课稿边界条件与载荷1.预设的温度*BOUNDARYTNODE,11,11,
11、500节点集第一个自由度温度最后个自由度温度值不变:变化的温度:*BOUNDARY,AMPLITUDE=amp-1TNODE,11,11,500温度幅值温度受幅值曲线amp-1控制11t0幅值曲线变化的温度5001t0T温度的共轭反作用是热率(热能进入一个已经预设温度值的节点的流通率)输出变量:RFLn16ABAQUS热传导与热应力分析课稿边界条件与载荷2.预设的热流量(热率)节点的集中热流量(与自由度11共轭)通过关键字*CFLUX施加*CFLUX,AMP=amp-1FNODE,11,30热率参考值输入可以参考一个AMPLITUDE曲线,使得输入的热率可以随时间变化。输出变量CFLn可以反映
12、节点*CFLUX的当前值。分布热流量(通过关键字*DFLUX或DSFLUX施加q*DFLUX可以施加在面或体上*DSFLUX只能施加在面上*DFLUX,AMP=amp-1ELHOL,S1,300*DSFLUX,AMP=amp-1SHOL,S,30017ABAQUS热传导与热应力分析课稿边界条件与载荷3.边界层(薄膜)条件-热传导中最常见的一种边界条件是一个自由表面被紧临的流体加热或降温-关键字*CFLIM,*FILM和*SFILM用于定义边界层条件。-边界层系数h是ABAQUS的一个输入参数,量纲:JL-2T-1q-1-边界层系数的重要性:热传导的结果严重依赖这个参数典型的,h是流体雷诺数和流
13、通温度的函数,但也与表面状况如粗糙度,脏污和方位强相关,因此很难去特征化。通常,需要用试验校准的方式来确定h的取值。Film,coefficienthq流体,温度q*FILMPROPERTY,NAME=H111.6E-6,4014.2E-6,6019.3E-6,80定义hh是温度q的函数18ABAQUS热传导与热应力分析课稿边界条件与载荷3.边界层(薄膜)条件*CFILM施加在节点上*CFILMNODESET,100.,450,2.3E-3*FILM二维情况下施加在单边上,三维情况下施加在单元面上*FILMELSET,F3.,450,2.3E-3面积温度h温度h*SFILM二维情况下施加在面上
14、*FILMSURSET,F.,450,2.3E-3温度h19ABAQUS热传导与热应力分析课稿边界条件与载荷4.向环境的辐射热传导中的另一种边界条件是黑体辐射q=-A(T4Te4)*CRADIATE施加在节点上*CRADIATENODESET,100.,450,0.1Emissivity(01)*RADIATE施加在单元上*RADIATEELSET,R1.,450,0.1*SRADIATE施加在面上*CRADIATESURSET,R.,450,0.1单元面编号定义辐射边界条件,需要定义Stefan-Boltzmann常数和绝对零度*PHYSICALCONSTANTS,ABSOLUTEZERO=
15、-273.16STEFANBOLTZMANN=5.6697E-820ABAQUS热传导与热应力分析课稿边界条件与载荷4.向环境的辐射辐射率emissivity是衡量一个表面有多接近理想黑体的指标一些常用材料的辐射率:Commercialaluminumsheet:0.09Heavilyoxidizedaluminumsheet:0.2Polishedgold:0.02Rustedironplate:0.6Polishedironplate:0.07Turned,heatedcastiron:0.44Type301stainlesssteel:0.58Redbrick:0.93Blackshin
16、ylacqueroniron:0.88Whitevamish:0.09Water:0.9521ABAQUS热传导与热应力分析课稿边界条件与载荷4.向环境的辐射是否需要考虑辐射边界条件FilmRadiationHeatfluxSurfacetemperature0200100Te=Roomtemp(23oC)h=10W/m2/oC辐射率=1温度越高,辐射现象越强22ABAQUS热传导与热应力分析课稿边界条件与载荷5.自然边界条件在任何温度下没有给定热流并没有外部热流的表面,默认条件是通过表面q=0,即没有通过表面的热流:理想绝热条件这是自然(无热载荷)边界条件,用于诸如施加对称边界条件的时候,如
17、外部热源温度800oC内部热源温度400oC23ABAQUS热传导与热应力分析课稿稳态分析实例二维热传导xy1.00.5ABCDE0.2Conductivity=52W/m/oCFilmcoefficient=750W/m2/oCBoundaryconditions:=100oCCalongABHeatflux=0alongDAConvectiontoambienttemperatureof0oCalongBCandCDObjective:FindqatETargetsolution:18.3oCatE24ABAQUS热传导与热应力分析课稿稳态分析实例定义热传导率定义薄膜换热系数换热条件边界条
18、件25ABAQUS热传导与热应力分析课稿稳态分析实例二维热传导26ABAQUS热传导与热应力分析课稿瞬态分析-有限元方法将问题在空间中离散化,对于瞬态传热问题,控制方程也必须通过时间积分进行求解-在ABAQUS中对瞬态固体传热进行时间积分的操作是利用后向差分算法:-后向差分算法是:相当的精确无条件稳定的-算法的稳定性非常重要,因为许多瞬态传热问题是在长的时间周期内进行分析的。(典型的是要到达到稳态条件)27ABAQUS热传导与热应力分析课稿瞬态分析-瞬态传热是扩散主导的过程在对一些对外界条件改变的响应中,开始时温度随时间的变化很快,然而到后期,可以看到温度的缓慢变化。-在ABAQUS传热分析中
19、,自动时间增量过程具有这种逻辑上内建的期望响应类型:指数衰减或增加。-这种结合精确设置DELTMX的方案,允许ABAQUS/Standard保持在所有分析阶段的整个过程中具有一致的精确性。28ABAQUS热传导与热应力分析课稿二维瞬态热传导例子xy1.00.5ABCDE0.2Conductivity=52W/m/oCSpecificheat=434J/kg/oCDensity=7832kg/m3Filmcoefficient=750W/m2/oCBoundaryconditions:=100oCCalongABHeatflux=0alongDAConvectiontoambienttemper
20、atureof0oCalongBCandCDObjective:FindqatETargetsolution:18.3oCatEatsteadystate.瞬态分析29ABAQUS热传导与热应力分析课稿瞬态分析在稳态算例基础上,增加密度和比热参数瞬态传热分析步设定DELTMAX30ABAQUS热传导与热应力分析课稿瞬态分析31ABAQUS热传导与热应力分析课稿瞬态分析-DELTMAX是一个时间积分精度参数在利用时间积分计算瞬态传热方程通过控制饿过程中,温度在每个时间最大允许的温度变化值,来控制求解的精度。配合使用自动时间增量方法,可以严格的控制时间增量步的大小,来满足DELTMAX的设定。-如
21、果计算过程中都能够满足DELTMAX,ABAQUS/Standard会尝试尽量增大时间增量步。自动时间增量步算法会尝试选择最优化的增量步时间,来兼顾计算精度和效率-瞬态传热分析可以通过设定当温度变化小于设定值时停止计算32ABAQUS热传导与热应力分析课稿瞬态分析-瞬态传热分析中的Initialconditions可以再瞬态传热分析之前,设定一个初始的温度分布如果没有给定初始值,abaqus的默认初始温度为033ABAQUS热传导与热应力分析课稿瞬态分析-最小可用时间增量步设置(仅ABAQUS/Standard适用)在对瞬态扩散过程的近似离散中,非常重要的一个问题是初始时间增量的选择。空间单元
22、的大小和时间增量步之间的关系是:如果一个时间增量步太小,将会产生很多无用信息,并且事实上还好经常出现一些虚假的震荡的结果。当使用二阶单元时,震荡会比较显著最小时间增量准则:Dt=时间增量r=密度c=比热k=热传导率Dl=在最大温度梯度区域靠近表面的单元尺度34ABAQUS热传导与热应力分析课稿非线性分析-一个典型的传热分析会包含以下的一些非线性:材料非线性:1.热传导率是温度的函数2.比热是温度的函数3.潜热效应,一种很强的非线性边界条件非线性:1.具有辐射边界条件,有时也是一种很强的非线性2.换热系数是温度的函数3.任何热流边界条件中,热流是温度的函数-ABAQUS/Standard使用牛-
23、拉迭代法,求解非线性问题的方程:35ABAQUS热传导与热应力分析课稿热“接触”-热量通过接触界面传导1.通过这些薄的界面进行传热是热分析的一个重要方面2.这些界面通常居于较低的导热率因此,在它们之间允许较大的温度差异3.然而,薄的界面具有可以忽略的“热质量”因此可以忽略界面内部的热能,假设它具有零比热36ABAQUS热传导与热应力分析课稿热“接触”-热界面的例子流体速度曲线流体边界层芯片核心载板表面1表面2在ABAQUS中,这种效应被模型化为:机遇面之间的相互关系,这些面(三维或二维)通常是物理上很相近,但两边具有不同的温度。这些面可以再不同物体上,也可以在同一物体上。37ABAQUS热传导
24、与热应力分析课稿热“接触”-热量可以通过以下方式穿过界面1.热传导,定义间隙热传导系数*GAPCONDUCTANCE2.热辐射,定义间隙热辐射率*GAPRADIATION-通常以上两种传热模式都存在,它们的相对重要关系取决于表面温度和界面间的媒介:1.界面间存在一定物质:热传导相对重要一些。因为热辐射不需要中间媒介,中间媒介的存在反而会吸收掉辐射热能。热“接触”如何定义1.定义表面2.定义定义接触对与机械接触相同。只是定义接触关系属性时需要用热的条件*SURFACEINTERACTION,NAME=THERM11.0*GAPCONDUCTANCE20,020,138ABAQUS热传导与热应力分
25、析课稿热应力分析39ABAQUS热传导与热应力分析课稿热应力分析过程-ABAQUS提供三种热应力分析程序:1.顺序耦合热应力分析,最常用的方法当应力是由热量场存在造成的,并且热求解过程与应力状态无关,也就是说应力依赖于热产生,而热并不依赖位移。需要跑两个分析:先分析热传导,再将温度结果导热应力分析热分析的结果,如温度(位置,时间的函数)被读入应力分析,作为一个预定义场。2.完全耦合热应力分析,最常用的方法应力依赖于温度场并且温度也依赖于应力场。只需要跑一个析。3.绝热分析这个分析的目的是模拟机械的变形产生局部的热量,但是时间很短,热传导不明显所有的温度增加是在材料局部发生的,并且只影响局部材料
26、点的机械属性 。40ABAQUS热传导与热应力分析课稿顺序耦合热应力分析-在热传导分析中,温度是未知量求解温度场-在应力分析中,位移是未知量节点的温度作为已知的外部载荷来产生热应变:-对于静力学分析计算应变:应力求解:这样,热场通过以下方式影响机械场:热膨胀(收缩)与温度相关的机械属性41ABAQUS热传导与热应力分析课稿顺序耦合热应力分析-热膨胀一个未约束体在温度增加的情况下,通常会产生体积膨胀,这种应变与温度之间的特征因子成为:热膨胀系数(CoefficientofThermalexpansion,CTE)CTE可以定义为与温度相关的,也可以是各向同性或各项异性使用*EXPANSION定义CTE*EXPANSION,TYPE=ISO,ZERO=201.0E-6,100,11.5E-6,200,12.0E-6,100,22.5E-6,200,2计算公式:q0:参考温度,qI,初始温度CTE会随温度的增加而缓慢增加42ABAQUS热传导与热应力分析课稿顺序耦合热应力分析-热膨胀算例43ABAQUS热传导与热应力分析课稿