Workbench瞬态热分析
ANSYSWorkbench在结构瞬态动力学分析中的应用_巨文涛

=
{ un }
2
+
{ un }
Δt ( 6)
其中: α 、δ 为 New mark 积分参数 在时刻控制方程 为了计算下一时刻的位移 u n + 1 , ( 2. 4 ) 为 [ M] C] K] { un + 1 } + [ { un + 1 } + [ { un + 1 } = { Fa } ( 7) 由( 5 ) 和( 6 ) 得 { u n + 1 } = a0 ( { u n + 1 } - { u n } ) - a2 { u n } - a3 { u n } ( 8) { u n + 1 } { u n } + a6 { u n } + a7 { u n + 1 } ( 9) 1 1 1 δ , a = , a = -1 , a4 a1 = 其中 a0 = 2, αΔt 2 αΔt 3 2 α αΔt δ Δt δ a5 = ( - 2) , = -1 , a6 = Δt( 1 - δ) ,a7 = Δtδ 2 α α ( 8 ) 和( 9 ) 得 由( 7 ) 、 M]+ a1[ C] + [ K] M]a0 ( a0[ ) { un + 1 } = { Fa } + [ { u n } + a2 { u n } + a3 { u n } + [ C] ( a1 { u n } + a4 { u n } + a5 ( 10 ) { un } ) 根据以上各式, 速 可以得到 t n + 1 时刻的位移 u n + 1 、 度 u n + 1 和 u n + 1 加速度。 利用式( 5 ) 和( 6 ) 得到的 New mark 求解方法的无 条件稳定必须满足: 1 1 1 1 + δ) 2 , ( 11 ) δ≥ , + δ + α > 0 α≥ ( 4 2 2 2 New mark 参数 1 1 2 ( 12 ) δ = +γ α = ( 1 + γ) , 4 2 其中: γ 为振幅衰减因子 通过观察( 11 ) 和( 12 ) 可以发现无条件稳定也可以 1 1 2 α≥ ( 1 + γ) 且 γ≥0 。因此只要 γ 表述为 δ = + γ, 2 4 ≥0 , 则求解就是稳定的。 2. 2 HHT 算法 HHT 时间积分法由下式给出 在完全瞬态分析中, [ M] { u n + 1 - αm } + [ { u n + 1 - αf } + [ { u n + 1 - αf } = C] K] { Fa n + 1 - αm } 其中 { u n + 1 - αm } = ( 1 - α m ) { u n + 1 } + α m { u n } { u n + 1 - αf } = ( 1 - α f ) { u n + 1 } + α f { u n } { u n + 1 - αf } = ( 1 - α f ) { u n + 1 } + α f { u n } ( 13 )
workbench热分析案例

•6
划分网 格
网格剖分: 采用ansys的mesh块对导入 的几何体进行网格划分,网 格为四面体网格,网格最大 边长为5mm。
•1
定 义 边 界条件
墙壁外表面: 采用convection边界条件, 设定外界空气温度10℃, 换热系数为0.36W/㎡·k。
•2
定义边界条件
墙壁内表面:
裸露于空气的表面采用 convection边界条件,拟 定外界空气温度20℃, 换热系数为0.36W/㎡·k, 与热源接触表面采用耦合 边界条件。
•3
定 义边界条件
热源: 与墙体平行的壁面采用 temperature边界条件,定 义其温度为50℃,其余壁 面均为绝热边界条件。
•4
结 果及分析
温度场云图:
通过显示计算得出的温度 场可以看出该模型的最小 温度值出现在墙体外表面 顶部与底部,在该模型中 温度场关于yz平面对称。
•
结 果及分析
ansysworkbench热分析研究教程

6-1?本章练习稳态热分析地模拟,包括:A.几何模型B.组件-实体接触C.热载荷D.求解选项E.结果和后处理F. 作业6.1? 本节描述地应用一般都能在ANSYS DesignSpaceEntra或更高版本中使用,除了ANSYSStructural? 提示:在ANSYS 热分析地培训中包含了包括热瞬态分析地高级分析?K T???T???Q?T?–在稳态分析中不考虑瞬态影响–[K]可以是一个常量或是温度地函数–{Q}可以是一个常量或是温度地函数?上述方程基于傅里叶定律:? 固体内部地热流(Fourier’s Law)是[K]地基础;? 热通量、热流率、以及对流在{Q}为边界条件;?对流被处理成边界条件,虽然对流换热系数可能与温度相关?在模拟时,记住这些假设对热分析是很重要地.?热分析里所有实体类都被约束:–体、面、线•线实体地截面和轴向在DesignModeler中定义? 热分析里不可以使用点质量(PointMass)地特性?壳体和线体假设:–壳体:没有厚度方向上地温度梯度–线体:没有厚度变化,假设在截面上是一个常量温度? 但在线实体地轴向仍有温度变化•唯一需要地材料特性是导热性(ThermalConductivity)•Thermal Conductivity在Engineering Data中输入?温度相关地导热性以表格形式输入若存在任何地温度相关地材料特性,就将导致非线性求解.?对于结构分析,接触域是自动生成地,用于激活各部件间地热传导–如果部件间初始就已经接触,那么就会出现热传导.–如果部件间初始就没有接触,那么就不会发生热传导(见下面对pinball地解释).–总结:–Pinball区域决定了什么时候发生接触,并且是自动定义地,同时还给了一个相对较小地值来适应模型里地小间距.?如果接触是Bonded(绑定地)或noseparation(无分离地),那么当面出现在pinballradius内时就会发生热传导(绿色实线表示).PinballRadius右图中,两部件间地间距大于pinball区域,因此在这两个部件间会发生热传导.? 默认情况下,假设部件间是完美地热接触传导,意味着界面上不会发生温度实际情况下,有些条件削弱了完美地热接触传导:––––压力表面温度T使用导电脂....Tx?接? ?– 穿过接触界面地热流速,由接触热通量q 决定:q ??TCC target ??T cont act ?– 式中T contact 是一个接触节点上地温度, T target 是对应目标节点上地温度–默认情况下,基于模型中定义地最大材料导热性KXX 和整个几何边界框地对角线ASMDIAG ,TCC 被赋以一个相对较大地值.TCC ??KXX ?10,000/ASMDI AG – 这实质上为部件间提供了一个完美接触传导? 在ANSYS Professional或更高版本,用户可以为纯罚函数和增广拉格朗日方程定义一个有限热接触传导(TCC).–在细节窗口,为每个接触域指定TCC输入值–如果已知接触热阻,那么它地相反数除以接触面积就可得到TCC值在接触界面上,可以像接触热阻一样输入接触热传导?Spotweld(点焊)提供了离散地热传导点:–Spotweld在CAD软件中进行定义(目前只有DesignModeler和Unigraphics 可用).T2T1?热流量:– 热流速可以施加在点、边或面上.它分布在多个选择域上.– 它地单位是能量比上时间(energy/time ) ?完全绝热(热流量为0): ?热生成:– 内部热生成只能施加在实体上– 它地单位是能量比上时间在除以体积(energy/time/volume )正地热载荷会增加系统地能量.– 可以删除原来面上施加地边界条件? 热通量:– 热通量只能施加在面上(二维情况时只能施加在边上)– 它地单位是能量比上时间在除以面积( e nergy/time/area )温度、对流、辐射:?完全绝热条件将忽略其它地热边界条件 ? 给定温度: – 给点、边、面或体上指定一个温度– 温度是需要求解地自由度? 至少应存在一种类型地热边界条件,否则,如果热量将源源不断地输入到系统中,稳态时地温度将会达到无穷大.? 另外,给定地温度或对流载荷不能施加到已施加了某种热载荷或热边界条件地表面上 .?对流:– 只能施加在面上(二维分析时只能施加在边上)– 对流q 由导热膜系数 h ,面积A ,以及表面温度T surface 与环境温度T ambient 地差值 来定义. q ? hA ?T surface ? T ambient ?– “h ” 和 “T ambient ” 是用户指定地值– 导热膜系数 h 可以是常量或是温度地函?与温度相关地对流:–为系数类型选择Tabular(Temperature)–输入对流换热系数-温度表格数据–在细节窗口中,为h(T)指定温度地处理方式?几种常见地对流系数可以从一个样本文件中导入.新地对流系数可以保存在文件中.?辐射:– 施加在面上(二维分析施加在边上)?4 4?– 式中: Q R ????FAT surface ??T ambient? σ=斯蒂芬一玻尔兹曼常数? ε =放射率? A =辐射面面积? F = 形状系数(默认是1)– 只针对环境辐射,不存在于面面之间(形状系数假设为1)– 斯蒂芬一玻尔兹曼常数自动以工作单位制系统确定从Workbench toolbox插入Steady-StateThermal将在projectschematic里建立一个SSThermalsystem(SS热分析)•在Mechanical 里,可以使用Analysis Settings为热分析设置求解选项.–注意,第四章地静态分析中地AnalysisDataManagement选项在这里也可以使用.加地结构载荷和约束.– 求解结构?为了实现热应力求解,需要在求解时把结构分析关联到热模型上. 在Static Structural 中插入了一个importedload 分支,并同时导入了施?后处理可以处理各种结果:–温度–热通量–反作用地热流速–用户自定义结果?模拟时,结果通常是在求解前指定,但也可以在求解结束后指定.–搜索模型求解结果不需要在进行一次模型地求解.?温度:–温度是标量,没有方向?可以得到热通量地等高线或矢量图:– 热通量 q 定义为q ???KXX ??T – 可以指定Total Heat Flux (整体热通量)和DirectionalHeatFlux (方向热通量) ? 激活矢量显示模式显示热通量地大小和方向?对给定地温度、对流或辐射边界条件可以得到响应地热流量:–通过插入probe指定响应热流量,或–用户可以交替地把一个边界条件拖放到Solution上后搜索响应从Probe菜单下选择或拖放边界条件?作业6.1–稳态热分析?目标:–分析图示泵壳地热传导特性版权申明本文部分内容,包括文字、图片、以及设计等在网上搜集整理.版权为个人所有This article includes some parts, including text, pictures, and design. Copyright is personal ownership.用户可将本文地内容或服务用于个人学习、研究或欣赏,以及其他非商业性或非盈利性用途,但同时应遵守着作权法及其他相关法律地规定,不得侵犯本网站及相关权利人地合法权利.除此以外,将本文任何内容或服务用于其他用途时,须征得本人及相关权利人地书面许可,并支付报酬.Users may use the contents or services of this article for personal study, research or appreciation, and othernon-commercial or non-profit purposes, but at the same time, they shall abide by the provisions of copyright law and other relevant laws, and shall not infringe upon the legitimate rights of this website and its relevant obligees. In addition, when any content or service of this article is used for other purposes, written permission and remuneration shall be obtained from the person concerned and the relevant obligee.转载或引用本文内容必须是以新闻性或资料性公共免费信息为使用目地地合理、善意引用,不得对本文内容原意进行曲解、修改,并自负版权等法律责任.Reproduction or quotation of the content of this article must be reasonable and good-faith citation for the use of news or informative public free information. It shall not misinterpret or modify the original intention of the content of this article, and shall bear legal liability such as copyright.。
ansysworkbench瞬态动力学实例

在本文中,我将为您撰写一篇关于ANSYS Workbench瞬态动力学实例的文章。
我们将深入探讨ANSYS Workbench在瞬态动力学仿真方面的应用,从简单到复杂、由浅入深地讨论其原理和实践操作,并共享个人观点和理解。
第一部分:介绍ANSYS Workbench瞬态动力学仿真ANSYS Workbench是一种用于工程仿真的全面评台,包含了结构、流体、热传递、多物理场等多种仿真工具。
瞬态动力学仿真是ANSYS Workbench的重要应用之一,它能够模拟在时间和空间上随机变化的动力学过程,并对结构在外部力作用下的动力响应进行分析。
在瞬态动力学仿真中,ANSYS Workbench可以模拟诸如碰撞、冲击、振动等动态载荷下的结构响应,用于评估零部件的耐久性、振动特性、动态稳定性等重要工程问题。
通过对这些现象的模拟和分析,工程师可以更好地了解结构在实际工况下的性能,进而进行有效的设计优化和改进。
第二部分:实例分析为了更直观地展示ANSYS Workbench瞬态动力学仿真的应用,我们以汽车碰撞仿真为例进行分析。
假设我们需要评估汽车前部结构在碰撞事故中的动态响应,我们可以通过ANSYS Workbench建立汽车前部结构的有限元模型,并对其进行碰撞载荷下的瞬态动力学仿真。
我们需要构建汽车前部结构的有限元模型,包括车身、前保险杠、引擎盖等部件,并设定材料属性、连接方式等。
接下来,我们可以在仿真中引入具体的碰撞载荷,如40km/h车速下的正面碰撞载荷,并进行瞬态动力学仿真分析。
通过仿真结果,我们可以获取汽车前部结构在碰撞中的应力、应变分布,以及变形情况,从而评估其在碰撞事故中的性能表现。
第三部分:个人观点与总结通过以上实例分析,我们可以看到ANSYS Workbench瞬态动力学仿真在工程实践中的重要应用价值。
瞬态动力学仿真不仅能够帮助工程师分析结构在动态载荷下的响应,还可以为设计优化、安全评估等工程问题提供重要参考。
ANSYS Workbench 17·0有限元分析:第12章-热分析

第12章 热分析 热力学分析(简称热分析)用于计算一个系统或部件的温度分布及其他各种热物理参数,如热量的获取与损失、热梯度、热流密度(热通量)等。
热分析在许多工程应用中扮演着非常重要的角色,如内燃机、涡轮机、换热器、电子元件等。
★ 了解传热的基础知识。
12.1 传热概述传热分析(Steady-State Thermal Analysis )遵循热力学第一定律,即能量守恒定律。
对于一个封闭的系统(没有质量的流入或流出),则:PE KE U W Q Δ+Δ+Δ=−式中Q 为热量,W 为所做的功,ΔU 为系统的内能,KE Δ为系统的动能,PE Δ为系统的势能。
对于大多数工程传热问题:0==PE KE ΔΔ若不考虑做功,即0=W ,则U Q Δ=;对于稳态热分析:0=Δ=U Q即流入系统的热量等于流出的热量;对于瞬态热分析:q dU dt =即流入或流出的热传递速率q 等于系统内能的变化。
12.1.1 传热方式热分析包括热传导、热对流、热辐射三种传热方式。
ANSYS Workbench 17.0有限元分析从入门到精通1.热传导热传导可以定义为完全接触的两个物体之间,或一个物体的不同部分之间由于温度梯度而引起的内能交换。
热传导遵循傅里叶定律:dxdT k q −=′′ 式中q ′′为热流密度(W/m 2),k 为导热系数。
2.热对流热对流是指固体的表面与它周围接触的流体之间,由于温差的存在引起的热量交换。
热对流可以分为两类:自然对流和强制对流。
热对流用牛顿冷却方程来描述:)(B T S T h q −=′′ 式中h 为对流换热系数(或称膜传热系数、给热系数、膜系数等),S T 为固体表面的温度,B T 为周围流体的温度。
3.热辐射热辐射是指物体发射电磁能,并被其他物体吸收转变为热的热量交换过程。
物体温度越高,单位时间内辐射的热量就越多。
热传导和热对流都需要有传热介质,而热辐射无须任何介质。
实质上,在真空中的热辐射效率最高。
ANSYS_WORKBENCH热分析

B. 装配体 – 实体接触
• 当导入实体零件组成的装配体时,实体间的接触区将 会被自动创建。
• 面-面接触允许实体零件间的边界上不匹配的网格。
– 接触实现了装配体中零件间的传热。
ANSYS License DesignSpace Entra De signS pa ce P rofe ssi ona l S tru ctu ra l M e cha nica l/M ulti physics
ANSYS License DesignSpace Entra De si gnS pa ce P rofe ssiona l S tructura l M e cha ni ca l /M ultiphysics
Availability
x x
x
… 装配体 – 接触区
• 正如前面的幻灯片所提到的,热量在接触区内沿着接 触法向流动
Initially Touching
Inside Pinball Region Outside Pinball Region
Yes
Yes
No
Yes
Yes
No
Yes
No
No
Yes
No
No
– 接触的 pinball 区域由程序自动定义并被设置一个相对较小的 值,以调和模型中可能出现的小间隙。pinball 区域将在下一
不接触,零件间将不会互相传热。
– 基于不同的接触类型,将热量是否会在接触面和目标面间传 递总结如下:
Contact Type
Bonded No Separation Rough F ric t ionles s
Heat Transfer Betw een Parts in Contact Region?
ansysworkbench热分析教程

. . -•本章练习稳态热分析的模拟,包括:A.几何模型B.组件-实体接触C.热载荷D.求解选项E.结果和后处理F. 作业6.1•本节描述的应用一般都能在ANSYS DesignSpaceEntra或更高版本中使用,除了ANSYSStructural•提示:在ANSYS 热分析的培训中包含了包括热瞬态分析的高级分析•对于一个稳态热分析的模拟,温度矩阵{T}通过下面的矩阵方程解得:[K(T)]{T}={Q(T)}•假设:–在稳态分析中不考虑瞬态影响–[K]可以是一个常量或是温度的函数–{Q}可以是一个常量或是温度的函数•上述方程基于傅里叶定律:•固体内部的热流(Fourier’s Law)是[K]的基础;•热通量、热流率、以及对流在{Q}为边界条件;•对流被处理成边界条件,虽然对流换热系数可能与温度相关•在模拟时,记住这些假设对热分析是很重要的。
•热分析里所有实体类都被约束:–体、面、线•线实体的截面和轴向在DesignModeler中定义•热分析里不可以使用点质量(PointMass)的特性•壳体和线体假设:–壳体:没有厚度方向上的温度梯度–线体:没有厚度变化,假设在截面上是一个常量温度• 但在线实体的轴向仍有温度变化• 唯一需要的材料特性是导热性(ThermalConductivity )• Thermal Conductivity 在Engineering Data 中输入•温度相关的导热性以表格形式输入若存在任何的温度相关的材料特性,就将导致非线性求解。
… 材料特性Training ManualB. 组件-实体接触Training Manual•对于结构分析,接触域是自动生成的,用于激活各部件间的热传导–如果部件间初始就已经接触,那么就会出现热传导。
–如果部件间初始就没有接触,那么就不会发生热传导(见下面对pinball的解释)。
–总结:–Pinball区域决定了什么时候发生接触,并且是自动定义的,同时还给了一个相对较小的值来适应模型里的小间距。
workbench瞬态动力分析

Dx IT20 L 波长方向的长度 c 弹性波速 E 杨氏模量 E
质量密度
非线性响应
• 非线性响应
–全瞬态分析可包括任何非线性类型. – 更小的 ITS 通常有助于平衡迭代收敛. – 塑性、蠕变及摩擦等非线性本质上是非保守的,需 要精确地遵循载荷加载历程.小的 ITS 通常有助于精 确跟踪载荷历程. – 小的ITS可跟踪接触状态的变化.
– 模态叠加法 – 直接积分法
• 运动方程可以直接对时间按步积分。在每个时间点(time = 0, Dt , 2Dt, 3Dt,….) ,需求解一组联立的静态平衡方程 (F=ma);
– 需假定位移、速度和加速度是如何随时间而变化的, (积分方案选择) – 有多种不同的积分方案,如中心差分法,平均加速度 法, Houbolt, WilsonQ, Newmark 等.
积分时间步长
• 如何选择 ITS? • 推荐打开自动时间步长选项 (AUTOTS), 并设置 初始时间步长Dtinitial和最小时间步长Dtmin 、最 大时间步长Dtmax. ANSYS 会利用自动时间步长 功能来自动决定最佳时间步长Dt. • 例如: 如果AUTOTS 是打开的, 并且Dtinitial= 1 sec, Dtmin= 0.01 sec, and Dtmax= 10 sec; 那 ANSYS 起始采用 ITS= 1 sec ,并依据结构的响 应允许其在0.01 和 10 之间变动.
缩减/完整结构矩阵
• 求解时既可用缩减结构矩阵,也可用完整结构矩阵; • 缩减矩阵:
– 用于快速求解; – 不允许非线性因素存在 – 根据主自由度写出[K]、[C]和[M]等矩阵,主自由度是完全自由度 的子集; – 缩减的 [K] 是精确的,但缩减的 [C] 和 [M] 是近似的。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
Workbench瞬态热分析
问题描述:将一个温度为900摄氏度的钢球放在空气中冷却,分别查看钢球和外部空气的温度变化。
分析类型:瞬态热分析分析平台:ANSYS Workbench 17.0分析人:技术邻一无所有就是打拼的理由研究模型:自定义
一、引言结构热分析主要包括热传导、热对流、热辐射,热分析遵循热力学第一定律,即能量守恒。
传热即是热量传递,凡是有温差存在的地方,必然有热量的传递。
传热现象在现实生活中普遍存在,比如食物的加热,冷却,有相变存在的蒸发冷凝换热等。
热分析类型主要有稳态热分析和瞬态热分析。
稳态热分析中,我们只关心物体达到热平衡状态时的热力条件,而不关心达到这种状态所用的时间。
在稳态热分析中,任意节点的温度不随时间的变化而变化。
一般来说,在稳态热分析中所需要的唯一材料属性是热导率。
在瞬态热分析中,我们只关心模型的热力状态与时间的函数关系,比如对水的加热过程。
在瞬态热分析中,需要对材料赋予热导率,密度,比热容等材料属性及初始温度,求解时间和时间增量这些边界条件。
在装配体的热分析中,我们还要考虑到接触区域传热,由于接触面可能存在表面粗糙度,接触压力等情况存在,导致存在接触热阻。
接触面存在两种传热方式,一种是附体间的热传递,另一种是通过空
隙层的热传导,但因为气体的热导率比较低,所以接触热阻不利于传热。
由于钢球散热与时间有关,我们选择瞬态热分析进行钢球的散热分析。
二、分析思路及流程在分析中,我们忽略空气的流动。
先进行稳态热分析,获得瞬态热分析的初始条件,然后将其传递到瞬态热分析中;在瞬态热分析中添加空气对流换热,来求解随时间变化的温度场。
分析流程如下图所示:
三、模型建立及网格划分:由于选取模型比较简单,我们在DM中建立一个钢球,选择钢球的半径为30mm,然后在外侧包络一层空气,包络厚度选择30mm,由于模型是对称的,为了节省计算时间,减少计算量,选取1/4模型进行研究(也可以选取1/8)。
由于模型较为简单,网格采用自动划分,模型及网格如下图所示:四、边界条件施加及结果分析:因为该问题为瞬态热分析,我们需要先进行稳态热分析获得瞬态热分析所需要的初始
条件,对钢球设置初始温度为900摄氏度,空气初始温度为22摄氏度,将稳态热分析的结果作为瞬态分析的初始条件,对空气对流换热系数为10W/m2K。
对瞬态热分析分为2个时间步,两个时间步分别设置为60s,因此钢球散热共计120s。
钢球在散热120s后的温度场如下图所示,从图中可以看出,钢球向空气散热120s后,钢球的最高温度为895.91摄氏度,靠近钢球侧的空气温度上升较为明显,基
本接近钢球温度。
离钢球越远处空气温度越低,最外侧空气最低温度为55.811摄氏度。
t=120s时温度云图t=120s时空气温度分布
钢球散热效果图
从下图中可以看出,钢球在120s中,在每个两个时间步里散热量基本一致,钢球散热处于平稳状态,即每秒中的散热量基本相同。
从空气外侧空气温度上升图表中可以看出,空气在前25秒内温度急剧上升到55.905摄氏度。
由于对流的存在,空气温度基本保持不变,钢球温度随时间的增加而降低。
钢球温度随时间下降过程空气最外侧温度随时间上升
五、总结本例通过ANSYS workbench 讲解了钢球瞬态散热问题的方法和具体应用。
在现实生活中,关于热的问题无处不在,比如我们烧热水,对事物的加热和冷却,保温等,复杂的热分析问题比如沸腾及冷凝等,涉及到相变,更为复杂。
在这个例子中,我们需要注意的是:在进行瞬态热分析时,一般情况下要先进行稳态热分析,以获取瞬态分析的初始天剑,然后通过瞬态热分析的分析设置及边界条件设置,进行多载荷步的求解。