Ansys 第33例瞬态热分析实例一水箱Word文档下载推荐.docx

Ansys 第33例瞬态热分析实例一水箱Word文档下载推荐.docx

《Ansys 第33例瞬态热分析实例一水箱Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《Ansys 第33例瞬态热分析实例一水箱Word文档下载推荐.docx（31页珍藏版）》请在冰点文库上搜索。

如果非均匀的初始温度场是已知的，可以用MainMenu→Solution→DefineLoads→Apply→InitialCondit'

n→Define即IC命令施加。

非均匀的初始温度场一般是未知的，此时必须先进行行稳态分析确定该温度场。

该稳态分析与一般的稳态分析相同。

要设定载荷（如已知的温度、热对流等），将时间积分关闭，选择MainMenu→Solution→LoadStepOpts→Time/Frequenc→TimeIntegration→AmplitudeDecay;

设定只有一个子步，时间很短（如（0.01s）的载荷步，MainMenu→Solution→LoadStepOpts→Time/Frequenc→Time→TimeStep。

（3）设置载荷步选项。

普通选项包括每一载荷步结束的时间、每一载荷步的子步数、阶跃选项等，选择MainMenu→Solution→LoadStepOpts→Time/Frequenc→Time-TimeStep.

非线性选项包括：

迭代次数（默认25），选择MainMenu→Solution→LoadStepOpts→Nonlinear→EquilibriumIter;

打开自动时间步长，选择MainMenu→Solution→LoadStepOpts→Time/Frequenc→Time→TimeStep：

将时间积分打开，选择MainMenu→Solution→LoadStepOpts→Time/Frequenc→TimeIntegration→AmplitudeDecay.

输出选项包括：

控制打印的输出，选择MainMenu→Solution→LoadStepOpts→OutputCtrls→SoluPrintout;

结果文件的输出，选择MainMenu→Solution→LoadStepOpts→OutputCtrls→DB/ResultsFile.

（4）如果均匀的初始温度场是通过稳态分析施加的，则须删除稳态分析时施加的温度载荷，选择MainMenu→Solution→DefineLoads→Delete→Thermal→Temperature。

（5）求解。

3.查看结果

可以用POST26或POST1查看结果。

33.2问题描述

图33-1所示为一个温度为500℃的铁块和一个温度为400℃的铜块，突然放入温度为20℃的完全绝热的水箱中。

忽略水的流动，试分析1h后铜块和铁块的最高温度，以及铜块和铁块的温度变化情况。

材料热物理性能参数如表33-1所示。

图33-1水箱示意图

表33-1材料热物理性能参数

热性能参数

铜

铁

水

导热系数/W（m.℃）

383

70

2

密度/kg/m3

8889

7833

966

比热/J/（kg.℃）

390

448

4185

33.3分析步骤

33.3.1改变任务名

拾取菜单UtilityMenu→File→ChangeJobname，在所弹出对话框的“[/FILNAM]”文本框中输入EXAMPLE33，单击“OK”按钮。

图33-2改变任务名对话框

33.3.2选择单元类型

拾取菜单MainMenu→Preprocessor→ElementType→Add/Edit/Delete，弹出如图33-3所示的对话框，单击“Add…”按钮，弹出如图33-4所示的对话框，在左侧列表中选“ThermalSolid”，在右侧列表中选“8node77”，单击“OK”按钮，返回到如图33-3所示的对话框，单击“Close”按钮关闭对话框。

图33-3单元类型对话框

图33-4单元类型库对话框

33.3.3定义村料模型

拾取菜单MainMenu→Preprocessor→MaterialProps→MaterialModels，弹出如图33-5所示的对话框，在右侧列表中依次拾取“Thermal”、“Conductivity”、“Isotropic”，弹出如图33-6所示的对话框，在“KXX”文本框中输入383（导热系数），单击“OK”按钮；

再次拾取如图33-5所示对话框中右侧列表的“SpecificHeat”项，弹出如图33-7所示的对话框，在“C”文本框中输入390（比热），单击“OK”按钮；

再次拾取如图33-5所示对话框中右侧列表的“Density”项，弹出如图33-8所示的对话框，在“DENS”文本框中输入8889（密度），单击“OK”按钮，于是定义好了材料模型1（铜）。

单击如图33-5所示对话框的菜单项Material--NewModel，单击弹出的“DefineMaterialTD”对话框中的“OK”按钮，然后重复定义材料模型1时的各步骤，定义材料模型2（铁）的导热系数为70，比热为448，密度为7833。

重复定义材料模型2时的各步骤，定义材料模型3（水）的导热系数为2，比热为4185，密度为996。

最后关闭如图33-5所示的对话框。

略

图33-5材料模型对话框

图33-6定义导热系数对话框

图33-7定义比热对话框

33.3.4创建矩形面，模拟铜块、铁块和水箱

拾取菜单MainMenu→Preprocessor→Modeling→Create→Areas→Rectangle→ByDimension，弹出如图33-9所示的对话框，在“X1，X2”文本框中分别输入0，0.6，在“Y1，Y2”文本框中分别输入0，0.5，单击“Apply”按钮，再次弹出如图33-9所示的对话框，在“Xl,X2”文本框中分别输入0.15,0.225，在“Y1,Y2”文本框中分别输入（0.225，0.27），单击“Apply”按钮，再次弹出如图33-9所示的对话框，在“Xl，X2”文本框中分别输入0.342，0.42），在“Y1，Y2”文本框中分别输入0.225，0.27，单击“OK”按钮。

33.3.5显示面号

拾取菜单UtilityMenu→PlotCtrls→Numbering，在弹出的"

PlotNumberingControls"

对话框中，将Areanumbers（面号）打开，单击“OK”按钮。

图33-8定义密度对话框

图33-9创建矩形对话框

33.3.6交叠面

拾取菜单MainMenu→Preprocessor→Modeling→Operate→Booleans→Overlap→Areas,弹出拾取窗口，单击“PickAll”按钮。

33.3.7划分单元

拾取菜单MainMenu→Preprocessor→Meshing→MeshTool，弹出如图33-10所示的对话框。

选择“ElementAttributes”的下拉列表框为“Areas”，单击下拉列表框后面的“Set”按钮，弹出拾取窗口，选择面2，单击拾取窗口中的“OK”按钮，弹出.”AreasAttributes"

对话框，选择“MAT”下拉列表框为1，单击“Apply”按钮，再次弹出拾取窗口，选择面3，单击拾取窗口中的“OK”按钮，选择“AreasAttributes”对话框的“MAT”下拉列表框为2，单击“Apply”按钮，再次窗口，择面4，单击拾取窗口中的“OK”按钮，选择“AreasAttributes”对话框的“MAT”下拉列表框为3，单击"

OK"

按钮。

单击“SizeControls”区域中“Global”后面的“Set”按钮，弹出如图33-11所示的对话框，在“SIZE”文本中输入0.01，单击“OK”按钮；

在如图33-10所示对话框的"

Mesh"

区域，选择单元形状为“Quad”（四边形），选择划分单元的方法为“Mapped”（映射），单击“Mesh”按钮，弹出拾取窗口，拾取面2和面3，单击“OK”按钮。

图33-10划分单元工具对话框

图33-11单元尺寸对话框

单击“SizeControls”区域中“Global”后面的“Set”按钮，弹出如图33-11所示的对话框，在"

SIZE"

文本框中输入0.03，单击“OK”按钮；

在如图33-10所示对话框的“Mesh”区域，选择单元形状为"

Quad"

（四边形），选择划分单元的方法为“Free”（自由），单击“Mesh”按钮，弹出拾取窗口，拾取面4，单击“OK”按钮，关闭如图33-10所示的对话框。

33.3.8指定分析类型

拾取菜单MainMenu→Solution→AnalysisType→NewAnalysis，在弹出的对话框中选择“TypeofAnalysis”为“Transient”，单击“OK”按钮，在随后弹出的“TransientAnalysis”对话框中，单击“OK”按钮。

以下步骤进行稳态分析，以得到初始温度场。

33.3.9设置时间积分为关闭，进行稳态分析，以得到瞬态分析的初始温度场

拾取菜单MainMenu→Solution→LoadStepOpts→Time/Frequenc→TimeIntegration→AmplitudeDecay，弹出如图33-12所示的对话框，将“TIMINT”关闭，单击“OK”按钮。

图33-12时间积分控制对话框．

33.3.10确定稳态分析的载荷步时间和吋间步长

拾取菜单MainMenu→Solution→LoadStepOpts→Time/Frequenc→Time→TimeStep,弹出如图33-13所示的对话框，在“TIME”文本框中输入0.01，在“DELTIMTimestepsize”文本框中输入0.01，单击“OK”按钮。

图33-13确定载荷步时间和时间步长对话框

33.3.11施加温度载荷

拾取菜单MainMenu→Solution→DefineLoads→Apply→Thermal→Temperature→OnAreas，弹出拾取窗口，拾取面2，单击“OK”按钮，弹出如图33-14所示的对话框，在“Lab2”列表中选择“TEMP”，在“VALUE”文本框中输入400，单击“Apply”按钮。

再次弹出拾取窗口，拾取面3，单击“OK”按钮，弹出如图33-14所示的对话框，在“VALUE”文本框中输入500，单击“Apply”按钮。

再次弹出拾取窗口，拾取面4，单击“OK”按钮，弹出如图33-14所示的对话框，在“VALUE”文本框中输入20，单击"

OK“按钮。

33.3.12求解

拾取菜单MainMenu→Solution→Solve→CurrentLS，单击“SolveCurrentLoadStep"

对话框中的"

OK”按钮。

图33-14在面上施加温度载荷对话框

以下步骤进行瞬态分析。

33.3.13确定瞬态分析的载荷步时间和时间步长

拾取菜单MainMenu→Solution→LoadStepOpts→Time/Frequenc→Time→TimeStep，弹出如图33-13所示的对话框，在“TIME”文本框中输入3600，在“DELTIMTimestepsize”文本框中输入50，选择“AUTOTS”为“ON”，在“DELTIMMinimumtimestepsize”文本框中输入10，在“DELTIMMaximumtimestepsize”文本框中输入200，单击“OK”按钮。

33.3.14打开吋间积分，进行瞬态分析

拾取菜单MainMenu→Solution→LoadStepOpts→Time/Frequenc→TimeIntegration-→AmplitudeDecay，弹出如图33-12所示的对话框，将“TIMJNT”打开，单击“OK”按钮。

33.3.15删除稳态分析吋施加的温度载荷

拾取菜单MainMenu→Solution→DefineLoads→Delete→Thermal→Temperature→OnAreas，弹出拾取窗口，单击“PickAll”按钮。

33.3.16确定数据库和结果文件中所包含的内容

拾取菜单MainMenu→Solution→LoadStepOpts→OutputCtrls→DB/ResultsFile，弹出“ControlsforDatabaseandResultsFileWriting”对话框，选择下拉列表框“Item”为“AllItems”，选中“Everysubstep”，单击“OK”按钮。

33.3.17求解

菜单MainMenu→Solution→Solve→CurrentLS，单击”SolveCurrentLoadStep”对话框中的“OK”按钮。

以下步骤查看结果。

33.3.18查看温度

拾取菜单MainMenu→GeneralPostproc→PlotResults→ContourPlot→NodalSolu，在列表中依次选择“NodalSolution→DOFSolution→NodalTemperature”，单击“OK”按钮，结果如图33-15所示。

图33-.153600s后的温度分布情况

33.3.19定义变量

拾取菜MainMenu→TimeHistPostpro→DefineVariables,弹出“DefinedTime-→HistoryVariables“对话框，单击“Add…”按钮，弹出“AddTime-HistoryVariables”对话框，选择“TypeofVariable”为“NodalDOFresult”，单击“OK”按钮，弹出拾取窗口，在拾取窗口的文本窗口中输入106（位于铜块中心），单击“OK”按钮，弹出如图33-16所示的对话框，在“Name”文本框中输入TEMP_106，单击“OK”按钮。

关闭“DefinedTime-HistoryVariables”对话框。

图33-16定义节点数据对话框

33.3.20用曲线图显示温度变化情况

拾取菜单MainMenu→TimeHistPostpro→GraphVariables，弹出“GraphTime-HistoryVariables”对话框，在“NVAR1”文本框中输入2，单击“OK”按钮，结果如图33-17所示。

图33-17节点106处的温度变化情况

第34例在结构上直接施加温度载荷进行热应力分析实例

------双金属簧片

本例介绍了在结构单元上直接施加温度载荷进行热应力分析的方法和步骤，并使用解析解对有限元分析结果进行了验证。

34.1概述

当一个结构发生温度变化时，它会发生膨胀或收缩。

如果其膨胀或收缩受到限制，或者由于温度分布不均而导致膨胀和收缩的程度不同，就会产生热应力。

ANSYS提供了三种计算热应力的方法，即在结构上直接施加温度载荷、间接法和直接法。

本例介绍第一种方法，该方法主要用于结构温度场为已知且均匀分布的场合。

分析类型为静力学分析，温度场的施加采用MainMenu→Solution→DefineLoads→Apply→Structural→Temperature命令。

此时，节点温度作为静力学分析的体载荷而不是节点自由度。

34.2问题描述及解析解

温度调节器所使用的双金属簧片的计算简图如图34-1所示。

两簧片等厚度、等宽度，材料不同，二者紧密连成一体，其左端可视为固定端。

第1个金属片为钢制，其弹性模量E1=2x1011Pa，线膨胀系数，α1=10x10-6/℃；

第2个金属片为铜制，其弹性模量E2=1.1x1011Pa，线膨胀系数a2=16xl0-'

/℃。

簧片各部分尺寸为长度L=40mm，高度h=0.5mm。

试分析当温度升高△t=l00℃时，簧片自由端的挠度。

图34-1双金属簧片

根据材料力学的知识，容易得出簧片自由端的挠度为

=7.038×

10-4m

34.3分析步骤

34.3.1改变任务名

拾取菜单UtilityMenu→File→ChangeJobname，在所弹出对话框的“[/FILNAM]”文本框中输入EXAMPLE34，单击“OK”按钮。

34.3.2选择单元类型

拾取菜单MainMenu→Preprocessor→ElementType→Add/Edit/Delete，弹出如图34-2所示的对话框，单击“Add…”按钮，弹出如图34-3所示的对话框，在左侧列表中选“StructuralSolid”，在右侧列表中选“20node186”，单击“OK”按钮，返回到如图34-2所示的对话框，单击“Close”按钮。

图34-2单元类型对话框

图34-3单元类型库对话框

34.3.3定义村料模型

拾取菜单MainMenu→Preprocessor→MaterialProps→MaterialModels，弹出如图34-4所示的对话框，在右侧列表中依次拾取“Structural”、“Linear”、“Elastic”、”Isotropic”，弹出如图34-5所示的对话框，在“EX”文本框中输入2e11（弹性模量），在“PRXY”文本框中输入0.3（泊松比），单击“OK”按钮；

再在图34-4所示对话框的右侧列表中依次拾取，Structural下面的“ThenualExpansion”、“SecantCoefficient”、“Isotropic”，弹出如图34-6所示的对话框，在“ALPX”文本框中输入10e-6（线膨胀系数），单击“OK”按钮。

单击如图34-4所示对话框的菜单项Material→NewModel，单击所弹出“DefineMaterialID”对话框中的“OK”按钮，然后重复定义材料模型1时的各步骤，定义材料模型2的弹性模量为1.1e11，泊松比为0.34，线膨胀系数为16e-6，最后关闭如图34-4所示的对话框。

图34-4材料模型对话框

图34-5材料特性对话框

（一）

图34-6材料特性对话框

（二）

34.3.4创建块

拾取菜单MainMenu→Preprocessor→Modeling→Create→Volumes→Block→ByDimension，弹出如图34-7所示的对话框，在“X1，X2”文本框中分别输入0，0.04，在“Y1，Y2”文本框中分别输入0，0.0005，在“21，22”文本框中分别输入0，0.005，单击“Apply”按钮。

再次弹出如图34-7所示的对话框，在“X1，X2”文本框中分别输入0，0.04，在“Y1，Y2”文本框中分别输入0.0005，0.00l，在“Z1，22”文本框中分别输入0，0.005，单击“OK”按钮。

图34-7创建块对话框

34.35黏结体

拾取菜单MainMenu→Preprocessor→Modeling→Operate→Booleans→Glue→Volumes.弹出拾取窗口，单击“PickAll”按钮。

34.3.6划分单元

拾取菜单MainMenu→Pteprocessor→Meshing→MeshTool，弹出如图34-8所示的对话框。

选择“ElementAttributes”的下拉列表框为“Volumes”，单击列表框后面的“Set”按钮，弹出拾取窗口．选择上面的六面体，单击拾取窗口中的“OK”按钮，弹出“VolumeAttributes”对话框，选择“MAT”下拉列表框为1，单击“Apply”按钮，再次弹出拾取窗口，选择下面的六面体，单击拾取窗口中的“OK”按钮，选择“VolumeAttributes”对话框的“MAT”下拉列表框为2，单击“OK”按钮。

单击“SizeControls”区域中“Global”后面的“Set”按钮，弹出如图34-9所示的对话框，在“SIZE”文本框中输入0.0005，单击“OK”按钮。

图34-8划分单元工具对话框

图34-9单元尺寸对话框

在图34-8所示对话框的“Mesh”区域，选择单元形状为“Hex”（六面体），选择划分单元的方法为“Mapped”（映射），单击“Mesh”按钮，弹出拾取窗口，单击“PickAll”按钮，最后单击如图34-8所示对话框中的“Close”按钮。

34.3.7改变视点

拾取菜单UtilityMenu→PlotCtrls→PanZoomRotate，在弹出的对话框中，依次单击“Iso”、“Fit”按钮，或者单击图形窗口右侧显示控制工具条上的剧按钮。

34.3.8施加约束

拾取菜单MainMenu→Solution→DefineLoads→Apply→Structural→Displacement→OnAreas，弹出拾取窗口，拾取x=0平面，单击“OK”按钮，弹出如图34-10所示的对话框，在列表中选择“AllDOF”，单击“OK”按钮。

图34-10在面上施加约束对话框

34.3.9施加温度载荷

拾取菜单MainMenu→Solution→DefineLoads→Apply→Structural→Temperature→OnVolumes，弹出拾取窗口，单击“PickAll”按钮，弹出如图34-11所示的对话框，在“VAL1”文本框中输入100，单击“OK”按钮。

图34-11在体上施加温度载荷对话框

34.3.10定义参考温度

拾取菜单MainMenu→Solution→DefineLoads→Settings-→ReferenceTemp，弹出如图34-12所示的对话框，在“TREF”文本框中输入0，单击“OK”按钮。

图34-12定义参考温度对话框

计算热膨胀大小，温度差等于步骤34.3.9所定义的节点温度减去参考温度，参考温度的默认值为0。

34.3.11求解

拾取菜单MainMenu→Solution→Solve→CurrentL