《Ansys-第33例瞬态热分析实例一水箱.pdf》由会员分享，可在线阅读，更多相关《Ansys-第33例瞬态热分析实例一水箱.pdf（30页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、-第第3 3 例例 瞬态热分析实例水箱瞬态热分析实例水箱本例介绍了利用 ANY进行瞬态热分析的方法和步骤、瞬态热分析时材料模型所包含的内容,以及模型边界条件和初始温度的施加方法。3 3概述概述热分析是计算热应力的基础,热分析分为稳态热分析和瞬态热分析，稳态热分析将在后面两个例子中介绍，本例介绍瞬态热分析。3 3.1.1.瞬态热分析的定义瞬态热分析的定义瞬态热分析用于计算系统随时间变化的温度场和其他热参数。一般用瞬态热分析计算温度场,并找到温度梯度最大的时间点，将此时间点的温度场作为热载荷来进行应力计算。33331 12 2嚼态热分析的步骤嚼态热分析的步骤瞬态热分析包括建模、施加载荷和求解、查看
2、结果等几个步骤。1.建模瞬态热分析的建模过程与其他分析相似,包括定义单元类型、定义单元实常数、定义材料特性、建立几何模型和划分网格等。注意:瞬态热分析必须定义材料的导热系数、密度和比热。2.施加载荷和求解(1)指定分析类型,ain MenuSutionayssTypeew alyss,选择 Tansien。(2）获得瞬态热分析的初始条件。定义均匀的初始温度场:Man MnuoutoDfn LoasSigsUnifor Temp,初始温度仅对第一个子步有效，而用 MinMenuSolutioneine LadppyThermalTepraue 命令施加的温度在整个瞬态热分析过程中均不变，应注意二
3、者的区别。定义非均匀的初始温度场：如果非均匀的初始温度场是已知的，可以用Man Meulutonfine Ldpplynital ConditnDefne 即C 命令施加。非均匀的初始温度场一般是未知的,此时必须先进行行稳态分析确定该温度场。该稳态分析与一般的稳态分析相同。注意:要设定载荷(如已知的温度、热对流等）,将时间积分关闭,选择ai Menuliono StpOpsTime/FreqencTime Ieraionmpliu-Day；设定只有一个子步,时间很短（如（01s)的载荷步,Main MenuSoutionLod teOpsTime/reqencTeTm Ste。()设置载荷步选
4、项。普通选项包括每一载荷步结束的时间、每一载荷步的子步数、阶跃选项等,选择 MainMenutinLad Step ptTim/FeqncTime-e Stp非线性选项包括：迭代次数(默认 25)，选择 Mi enuSolutionLoad St OpsNolearEqilibrum Ite;打开自动时间步长,选择Main MeolutionLoadStep OtsTim/FeuecimTime te：将时间积分打开，选择Man MenuSolutooadSteOptsTm/FrequenmeItegrationAmpltude Deay.输出选项包括：控制打印的输出，选择Min uSolut
5、oLo SepOptstput trSouPrntout;结果文件的输出,选择MaiMuluionLoad Stp Optsutput CtrlsDB/Results ie.()如果均匀的初始温度场是通过稳态分析施加的,则须删除稳态分析时施加的温度载荷,选择 Min eSoluionDei LoadsDeeTermalemperatue。(5)求解。3.查看结果可以用ST26 或 POST1 查看结果。33.233.2 问题描述问题描述图 33-1 所示为一个温度为500的铁块和一个温度为0的铜块，突然放入温度为20的完全绝热的水箱中。忽略水的流动，试分析 1h 后铜块和铁块的最高温度,以及铜
6、块和铁块的温度变化情况。材料热物理性能参数如表3-1 所示。-图3-水箱示意图表 331材料热物理性能参数热性能参数铜铁77344水2966485导热系数/W(m)38密度kg/m3比热/J/（kg)8883033.333.3 分析步骤分析步骤3.3.13.3.1 改变任务名改变任务名拾取菜单 Utlity Men Jobe，在所弹出对话框的“/FILNAM”文本框中输入 EXAMPE3，单击“O”按钮。图 33-改变任务名对话框33.333.32 2 选择单元类型选择单元类型拾取菜单 Main MenureprosoElement TypeAdd/dit/Deee,弹出如-图 33-3 所示
7、的对话框,单击“Add”按钮,弹出如图334所示的对话框,在左侧列表中选“Thema Sold”，在右侧列表中选“8ode 77”，单击“OK”按钮，返回到如图3-3 所示的对话框，单击“lse”按钮关闭对话框。图 33-3单元类型对话框图 3-单元类型库对话框3.3.3.3.定义村料模型定义村料模型拾取菜单 MainMeuPrrocsorMaerial ProsMateria del,-弹出如图 3-5 所示的对话框,在右侧列表中依次拾取“el”、“Codctivty”、“Isotroc”,弹出如图 33-6 所示的对话框，在“K”文本框中输入 383(导热系数)，单击“”按钮;再次拾取如图
8、 33-5 所示对话框中右侧列表的“Speii Heat”项，弹出如图 33-7 所示的对话框，在“”文本框中输入39(比热）,单击“K”按钮；再次拾取如图 33-所示对话框中右侧列表的“Density”项,弹出如图3-8 所示的对话框,在“DN”文本框中输入 888(密度），单击“K”按钮,于是定义好了材料模型1(铜）。单击如图33-所示对话框的菜单项Material-Neodel，单击弹出的“eine atea TD”对话框中的“OK”按钮，然后重复定义材料模型时的各步骤,定义材料模型 2(铁）的导热系数为0，比热为48，密度为833。重复定义材料模型 2 时的各步骤,定义材料模型 3（水
9、）的导热系数为2,比热为1，密度为996。最后关闭如图3-5 所示的对话框。略图 35材料模型对话框图 33-定义导热系数对话框-图3定义比热对话框3.43.4 创建矩形面，模拟铜块、铁块和水箱创建矩形面，模拟铜块、铁块和水箱拾取菜单aiMenPrprosr ModlinCate ArasRctangleBy Dimension，弹出如图33-9 所示的对话框，在“1,X2”文本框中分别输入,0.6,在“Y1,Y2”文本框中分别输入 0，05，单击“Apl”按钮，再次弹出如图 339 所示的对话框,在“Xl，X2”文本框中分别输入.15,025，在“Y1,Y2”文本框中分别输入(0 2,.2）
10、,单击“Appy”按钮,再次弹出如图 33-9 所示的对话框,在“Xl，X2”文本框中分别输入 42,0.4）,在“Y1,Y”文本框中分别输入 0.225，0.27,单击“O”按钮。3.3.5 显示面号拾取菜单tilit uPotCtrNumbering，在弹出的PotNumbrg Controls 对话框中,将ra umbes(面号)打开，单击“K”按钮。-图 33-定义密度对话框图 339创建矩形对话框33.3.633.3.6 交叠面交叠面拾取菜单ainenuPeproessoMdelneratBoolensvelapAreas，弹出拾取窗口，单击“PickAll”按钮。33.3.733.
11、3.7 划分单元划分单元拾取菜单 Mi MenuPreprocesrMhingMeshTool，弹出如图 3310 所示的对话框。选择“lemet Attribtes”的下拉列表框为“Aes”,单击下拉列表框后面的“et”按钮，弹出拾取窗口，选择面，单击拾取窗口中的“O”按钮，弹出”Area Atriuts对话框,选择“A”下拉列表框为1,单击“ply”按钮，再次弹出拾取窗口，选择面3,单击拾取窗口中的“OK”按钮,选择“reasAttrbutes”对话框的“AT”下拉列表框为，单击“ppl”按钮，再次窗口，择面 4,单击拾取窗口中的“OK”按钮,选择“Areas Attribtes”对话框的
12、“A”下拉列表框为，单击OK按钮。单击“Sz Control”区域中“Glbal”后面的“Set”按钮，弹出如图 33-11 所示-的对话框，在“SIZ”文本中输入0.01,单击“O”按钮；在如图3-1所示对话框的Ms区域,选择单元形状为“Qad”（四边形），选择划分单元的方法为“Mppd”(映射)，单击“sh”按钮，弹出拾取窗口，拾取面2 和面 3,单击“”按钮。图3-10划分单元工具对话框-图 33-11单元尺寸对话框单击“Siz otros”区域中“Glbal”后面的“Set”按钮，弹出如图3-11 所示的对话框,在 SIZE文本框中输入.03,单击“K”按钮；在如图310所示对话框的“
13、Mes”区域,选择单元形状为Qa（四边形)，选择划分单元的方法为“Free”(自由)，单击“Mesh”按钮,弹出拾取窗口,拾取面，单击“OK”按钮,关闭如图 33-1所示的对话框。3 3.3.3.指定分析类型指定分析类型拾取菜单ainMnuSlutonAalyis Typew Analysis,在弹出的对话框中选择“Tpeo nalys”为“Tranient”,单击“OK”按钮,在随后弹出的“TransienAnalys”对话框中,单击“OK”按钮。以下步骤进行稳态分析,以得到初始温度场。33.3.933.3.9 设置时间积分为关闭，进行稳态分析设置时间积分为关闭，进行稳态分析,以得到瞬态分析
14、的初始温度场以得到瞬态分析的初始温度场拾取菜单Mai MnuSltioLoadStep psTierequenTime InegrationApie ecay，弹出如图33-12 所示的对话框,将“TIINT”关闭，单击“OK”按钮。-图 33-2时间积分控制对话框.3 3.3.1.3.1确定稳态分析的载荷步时间和吋间步长确定稳态分析的载荷步时间和吋间步长拾取菜单 Main MenuSlutoLadSe tTime/FrequncmeTeStp,弹出如图 33-3 所示的对话框，在“TIME”文本框中输入 0.01,在“DLTIMTm tep size”文本框中输入.01,单击“OK”按钮。-
15、图3-1确定载荷步时间和时间步长对话框3 3.1.1施加温度载荷施加温度载荷拾取菜单 MinnSouionDfine asApplyThermalTemeratuen Aras,弹出拾取窗口，拾取面2,单击“OK”按钮，弹出如图3-14 所示的对话框，在“ab2”列表中选择“TEMP”，在“VU”文本框中输入 400，单击“Apl”按钮。再次弹出拾取窗口，拾取面3，单击“OK”按钮，弹出如图314 所示的对话框，在“VALUE”文本框中输入00,单击“Apl”按钮。再次弹出拾取窗口，拾取面 4,单击“K”按钮,弹出如图 314 所示的对话框,在“VAUE”文本框中输入0,单击O“按钮。33.3
16、.1233.3.12 求解求解拾取菜单ainnu SioSolvurretLS,单击“Slve Cur-nt LoSte对话框中的OK”按钮。图 334在面上施加温度载荷对话框以下步骤进行瞬态分析。3 3.3.13.3.13 确定瞬态分析的载荷步时间和时间步长确定瞬态分析的载荷步时间和时间步长拾取菜单 Main MenSoutioLoa Stp ptTim/Frqecimei tp，弹出如图33-13 所示的对话框，在“TIME”文本框中输入300，在“DLTIMTie stpsize”文本框中输入 50,选择“AUS”为“N”，在“DELTM nimumtmestepi”文本框中输入,在“D
17、ELTM Mximumtmstep iz”文本框中输入0，单击“OK”按钮。33.3.33.3.打开吋间积分，进行瞬态分析打开吋间积分，进行瞬态分析拾取菜单 Main MenuSolutonLoad Step OptsTime/FuencTime Intratin-mpiue Dea，弹出如图3-12 所示的对话框,将“TIMJN”打开，单击“OK”按钮。33.333.3删除稳态分析吋施加的温度载荷删除稳态分析吋施加的温度载荷拾取菜单anMenuSoutoDeineoadseeteThermalemperature ea,弹出拾取窗口,单击“ick ll”按钮。-33.333.31616确定数
18、据库和结果文件中所包含的内容确定数据库和结果文件中所包含的内容拾取菜单 Main MenuSlutiooad Ste OptOutpt CtrlsDBsultsFile，弹出“Conrls forDtbsad Rults”对话框,选择下拉列表框“Item”为“All Ies”,选中“Evry suep”,单击“K”按钮。33.3.1733.3.17 求解求解菜单ain MnSolutionSlveent L,单击”SolvCurrent Load Sep”对话框中的“OK”按钮。以下步骤查看结果。3 3.3.31 1查看温度查看温度拾取菜单 MainMeuGnra PotprocPlt esl
19、Contour Plooda Solu，在列表中依次选择“odl SolutinDO SolutoNoal Temperat”,单击“O”按钮,结果如图 33-1所示。图 33-.560 s 后的温度分布情况3 3.3.19.3.19 定义变量定义变量拾取菜 Main enuTimeist stpoDefineVarals,弹出“Defed Tme-Histoyariables“对话框,单击“Add”按钮，弹出“A ime-Histor rables”对话框,选择“Typ fVarable”为“Nodl D esult”，单击“-K”按钮,弹出拾取窗口,在拾取窗口的文本窗口中输入6（位于铜块中
20、心）,单击“O”按钮，弹出如图 33-16 所示的对话框,在“Nam”文本框中输入 TP_10,单击“O”按钮。关闭“Defined Tieisoy Varible”对话框。图 33定义节点数据对话框3 33 32020 用曲线图显示温度变化情况用曲线图显示温度变化情况拾取菜单 Main MuTimHs Pstroraph aribles，弹出“rapTm-ity Vriabls”对话框,在“NVR”文本框中输入2，单击“K”按钮,结果如图 3-17 所示。-图 33-17 节点6 处的温度变化情况第第3434例例在结构上直接施加温度载荷进行热应力分析实例在结构上直接施加温度载荷进行热应力分析
21、实例-双金属簧片双金属簧片本例介绍了在结构单元上直接施加温度载荷进行热应力分析的方法和步骤,并使用解析解对有限元分析结果进行了验证。3 3.1.1 概述概述当一个结构发生温度变化时，它会发生膨胀或收缩。如果其膨胀或收缩受到限制,或者由于温度分布不均而导致膨胀和收缩的程度不同,就会产生热应力。AYS 提供了三种计算热应力的方法，即在结构上直接施加温度载荷、间接法和直接法。本例介绍第一种方法,该方法主要用于结构温度场为已知且均匀分布的场合。分析类型为静力学分析,温度场的施加采用 Main MenuSolutionDeie LoadsAppyStucuralTeperature命令。此时，节点温度作
22、为静力学分析的体载荷而不是节点自由度。4.24.2 问题描述及解析解问题描述及解析解温度调节器所使用的双金属簧片的计算简图如图3-所示。两簧片等厚度、等宽度，-材料不同，二者紧密连成一体,其左端可视为固定端。第个金属片为钢制,其弹性模量1=21011a,线膨胀系数，=10 x1/；第2 个金属片为铜制,其弹性模量 E2=.1 xPa,线膨胀系数 a2=x-/。簧片各部分尺寸为长度L40mm，高度 h=0 5mm。试分析当温度升高=l00时,簧片自由端的挠度。图 34-1双金属簧片根据材料力学的知识，容易得出簧片自由端的挠度为=h（E12+E12+14E1E2）6（21）tE1E2L27.014
23、m3 3.3.3 分析步骤分析步骤3 3.1.1 改变任务名改变任务名拾取菜单tiliy e Jobae,在所弹出对话框的“/ILNAM”文本框中输入 EMPL34,单击“K”按钮。3 33.23.2 选择单元类型选择单元类型拾取菜单ain MenuPrproessoEleent pdEdit/lete，弹出如图-2 所示的对话框,单击“Add”按钮,弹出如图 33 所示的对话框,在左侧列表中选“Sructul old”,在右侧列表中选“20noe 86”，单击“OK”按钮,返回到如图 34-2 所示的对话框,单击“los”按钮。-图-2单元类型对话框图 34-3单元类型库对话框34.3.33
24、4.3.3 定义村料模型定义村料模型拾取菜单 MaienuPreproceoMatrial PpsMaterial Moels,弹出如图 34-所示的对话框，在右侧列表中依次拾取“Structra”、“Liner”、“Elast”、”Iropic”,弹出如图 35 所示的对话框,在“EX”文本框中输入e11（弹性模量），在“PXY”文本框中输入 0.(泊松比），单击“O”按钮;再在图-所示对话框的右侧列表中依次拾取，Structral 下面的“ThnulExpanson”、“ecat Cefcien-”、“topic”,弹出如图 3所示的对话框,在“ALPX”文本框中输入 1-6(线膨胀系数)
25、，单击“O”按钮。单击如图 3-所示对话框的菜单项 MterialN od，单击所弹出“Defin Matial ID”对话框中的“K”按钮，然后重复定义材料模型 1 时的各步骤，定义材料模型 2 的弹性模量为 1.e11,泊松比为 0.3，线膨胀系数为16e-6，最后关闭如图4 所示的对话框。图 34材料模型对话框图 34-材料特性对话框(一)-图46材料特性对话框（二)34.34.4 4 创建块创建块拾取菜单ainMnuPrepocessordelinreateVolumeslky Dimnsion，弹出如图 347 所示的对话框，在“X1，”文本框中分别输入0，0.04，在“，Y2”文本
26、框中分别输入 0，0 00 5，在“2,2”文本框中分别输入,0.005，单击“ppy”按钮。再次弹出如图 347 所示的对话框，在“X，X2”文本框中分别输入 0，.0,在“1，Y2”文本框中分别输入 00005，.0l，在“Z，22”文本框中分别输入 0，.05，单击“OK”按钮。图4-7创建块对话框34.334.3 黏结体黏结体拾取菜单 MaMenuPrerocessorModng OperateBoolean-GluVoles.弹出拾取窗口，单击“ck All”按钮。3 3.3.3划分单元划分单元拾取菜单 Man MenutepocssoMeshingMehTool,弹出如图 34-所
27、示的对话框。选择“Elment ttribues”的下拉列表框为“Vlmes”，单击列表框后面的“Set”按钮,弹出拾取窗口.选择上面的六面体,单击拾取窗口中的“K”按钮,弹出“Volme Atributes”对话框,选择“MT”下拉列表框为,单击“Appy”按钮,再次弹出拾取窗口，选择下面的六面体,单击拾取窗口中的“OK”按钮，选择“oume tribues”对话框的“MAT”下拉列表框为 2,单击“OK”按钮。单击“Size Cntol”区域中“lob”后面的“Se”按钮,弹出如图 3所示的对话框,在“SZE”文本框中输入 0.0005,单击“OK”按钮。-34-8划分单元工具对话框-图-
28、图-单元尺寸对话框在图 348 所示对话框的“Mesh”区域,选择单元形状为“Hex”(六面体）,选择划分单元的方法为“apped”(映射)，单击“Mesh”按钮，弹出拾取窗口，单击“Pi Al”按钮,最后单击如图 34-8 所示对话框中的“Close”按钮。4.34.37 7 改变视点改变视点拾取菜单 Utity MePlttrsZoomae,在弹出的对话框中，依次单击“Iso”、“Ft”按钮,或者单击图形窗口右侧显示控制工具条上的剧按钮。34.3.834.3.8 施加约束施加约束拾取菜单 Main MnSluionDeine LoadsApplStructurllacetOn Areas，
29、弹出拾取窗口,拾取 x=0 平面,单击“K”按钮,弹出如图4-0 所示的对话框,在列表中选择“llOF”，单击“K”按钮。图 34-0在面上施加约束对话框-3 3.施加温度载荷施加温度载荷拾取菜单 Man Menoluionefne oaAStturlmperatureOnVoles,弹出拾取窗口,单击“Pick All”按钮,弹出如图 34所示的对话框，在“VAL1”文本框中输入 10,单击“OK”按钮。图 34-11在体上施加温度载荷对话框34.3.134.3.1定义参考温度定义参考温度拾取菜单 Main MnSoutinDefiLoadsSttins-Rferenc emp，弹出如图 3
30、-1所示的对话框,在“TF”文本框中输入 0，单击“”按钮。图 3-12定义参考温度对话框计算热膨胀大小，温度差等于步骤3439 所定义的节点温度减去参考温度,参考温度的默认值为 0。3 3.3.11.3.11 求解求解拾取菜单ainMeuSoltonolveCrretLS,单击”Solv Crrntoa Step对话框中的“OK”按钮，当出现“oluion is dn!”提示时，求解结束,即可查看结果。4.3.124.3.12 查看结果，显示变形查看结果，显示变形拾取菜单Mn nuGeral tprocPlt ResltDefored-Shap，弹出如图4-13 所示的对话框，选中“Defu
31、ndef ee”(变形+未变形的模型边界），单击“O”按钮，结果如图4-1所示。图 34-13显示变形对话框图 34-14簧片的变形可见,簧片自由端的位移即簧片的最大位移为33 x-4mm，与理论解基本一致，从而,证明有限元解具有较高的精度.第第5 5 例例基于实测温度场进行热应力基于实测温度场进行热应力分析实例分析实例-转炉托圈转炉托圈本例通过实例介绍了利用NSS 参数化设计语言(AL)编制程序,将实测温度场施加到有限元模型上,并对结构的热应力进行有限元分析，得到了结构各部分热应力的精确分布。35.135.1 概述概述托圈是转炉的重要传动和承载部件,工作过程中的热物理过程十分复杂，其温度场分
32、布通过解析法或数值法进行计算几乎是不可能的，通常托圈的温度场只能是通过现场测试得到,但由于托圈结构比较复杂,其热应力计算一般都是用有限元方法进行的，因此存在着如何将实测得到的温度场施加到有限元模型上的问题。-某转炉托圈形状可近似为如图35-l(a)所示的回转体,图 35(b)为其横截面形状，已知 r=4.18、r1=5.155m、h=2.5m，经测量托圈在、0。、8。、270。四个横截面处 A、B、C、四点的溫度如表 35-1 所示,沿托圈圆周和横截面高度方向温度近似线性分布，现计算其热应力。图 3-转炉托圈的横截面表 35托圈温度分布横截面位置90180180270A 点1411711 12
33、12B115909012C2213123123168087 735由于托圈的温度按线性规律分布,所以有限元分析时可以由、9、180、270四个横截面的温度值通过线性插值得到其他截面上A、B、C、四点的温度值,再由该截面上 A、B、C、D 四点的温度值线性插值得到截面上其他点的温度值。于是,由若干个测量点的温度值插值得到整个托圈上所有点的温度值。使用AS对托圈热应力进行分析时,可以利用其参数化设计语言(APL)编制程序来模拟上述线性插值过程,然后把所计算得到的各个节点的温度值作为载荷施加到有限元模型上。程序框图如图52 所示。35352 2 命令流命令流CEAR！清除数据库，新建分析-/FILN
34、AME,EXAMLE35!定义任务名为”AMPLE35“R0=4.185!定义参数,尺寸R1=5.135=.5H=0./PRP7！进入预处理器E,1,PLANE182!选择单元类型ET,2,SOLID1MP，EX,1,2E11!定义材料模型，弹性模量MP，PY，1,.3!泊松比M，ALPX,3.5E6!线膨胀系数RECG,R0,R0+TH，-/2,H/2RETING，R1TH，R1,-H2,/ECTI，R0,/-T,H/RCTIG，R0,R1，-2+H，-H/AOVLAP，ALL!交叠所有面MSHP,0!六面体单元MSH,1！映射网格ESIE，.05!单元边长度MESH,AL!对面划分单元K,
35、100!创建关键点K,101,0,XTOP,ESIZE,10！挤出段数EXTOPT，ACE,1！挤出后清除面单元！旋转挤出!创建矩形面VROTA,ALL,100,101,60INI!退出预处理器/OU！进入求解器NSEL,，LOC,X，!在选择的节点上施加位移约束,AL,UXNSEL,S,L,Y,-D,LL，YSEL,S,，Z,0,L,ULLS！选择所有A_0=!定义参数，实测温度_0=88C_0=12TD_0=8TA_1790115TC90=3D_0=8T_18138180=90TC_80=2TD_8=7T_270=212TB_27=142TC2706D_2035*GET，NOE_X,NOD
36、E,，NUM,MAX!获得最大节点号CSYS,5!切换活跃坐标系为全局圆柱坐标系（轴)*DO，I,1,NOMAX!循环开始,计算每一个节点的温度值*GET，R,NOE，,LOC,!获得节点 I 的 r 坐标*ET，NIT，NOD,I,LOC,Y!获得节点 I 的坐标*GET,Z,OD,I，LOC，Z!获得节点 I 的 z 坐标*I，NSI,E,90,TET=TA_10(T_-TA_1)/0*（NSI-80)!计算节点 I 所在横截面上 A、B、-!D 四点的温度值T=TB_180-（TB_90TB_180)/90*(NSIT-180）TC=TC_180-(T_0-T_10)/9（NSI180)
37、TD=T_10-（D90-TD180）90(ST-180)*EEI,NST,GE,TA=T_0(T_0-TA_0)/0*(SIT-90）TB=TB90-(T_-TB_90)/90*(NSI-90)C=TC_90（T_0TC_90)/90*（IT-90）T=TD0(TD_0TD_90)90*（NIT-0)*EEF,N,G,-90TA=A_-（T270-TA_)/90(NSIT-)TB=TB_-(TB_2-B_0)/90（NSIT-0)TC=TC_0-（_20C_0）90*（NSIT-0）TD=TD0-（D270TD)/0(NSI-0)ELEIF,NIT,GE,180TA=A270(TA_0-T_
38、70)/90(NIT90)TT_27-（TB_18B_270）/90*(NSIT90)TC=TC_27-(T_8_270)90(N90)TD=T_70-(TD1-D20)/90(NI）EDIFTP0TA+(T-TA)*（NR-0）/（RR0)!由 A、B,，D 四点的温度值计算 P0和 P1 点的温度值T1=TB+（TD-TB)*(N-0)/(R1-R)TP=T1-(TP0-TP1)*(NZ-/2)/H!由 P0 和 P点的温度值计算 P 点即节点 P1 的温度值,I,TEM,TP!施加温度载荷ENDDOSLVE!求解AVE!保存数据库FNISH!退出求解器-PO1!进入普通后处理器PLNSO,BE,TEMPPLNOL,，SM!显示温度分布，如图-所示！显示变形，如图 3-4 所示INH!退出普通后处理器图5-3温度分布-图 35-4变形情况-