热应力耦合分析

TrainingManual热结构耦合场分析INTRODUCTIONTOANSYS11.0TrainingManual热应力产生•结构受热或变冷时，由于热胀冷缩产生变形。•若变形受到某些限制—如位移约束或相反的压力—则在结构中产生热应力。•产生热应力的另一个原因，是由于材料不同而形成的不均匀变形（如，不同的热膨胀系数）。由约束产生热应力由不同材料产生热应力–耦合场分析…概述INTRODUCTIONTOANSYS11.0TrainingManual•在ANSYS中求解热-应力问题有两种方法。这两种方法各有所长。–顺序耦合-传统方法是使用两种单元类型，将热分析的结果作为结构的温度荷载。+当热瞬态分析时间点很多，但结构时间点很少时效率较高。+很容易用输入文件实现自动处理。–直接耦合+比较新的方法，用一种单元类型就能求解两种物理场问题。+热和结构之间可实现真正的耦合。-在某些分析中可能耗费过多开销。–耦合场分析…概述INTRODUCTIONTOANSYS11.0TrainingManual•顺序耦合方法涉及两种分析：1.先作稳态（或瞬态）热分析。•建立热单元模型。•施加热荷载。•求解并检查结果。2.然后作静力结构分析。•把单元类型转换成结构单元。•定义包括热膨胀系数在内的结构材料特性。•施加包括从热分析得到的温度在内的结构荷载。•求解并检查结果。热分析结构分析jobname.rthjobname.rst温度–耦合场分析顺序耦合方法INTRODUCTIONTOANSYS11.0TrainingManual1.热分析–该过程在第11章中描述。2.结构分析a)进入前处理，把热单元类型转换成结构单元。•MainMenuPreprocessorElementTypeSwitchElemType注意：转换单元类型时，所有单元选项重新设置回原来缺省状态。例如，若用户在热分析中使用的2-D轴对称单元，则需要在转换后重新指定轴对称选项，因此，一定要确保设置正确的单元选项：•MainMenuPreprocessorElementTypeAdd/Edit/Delete[Options]•或用ETLIST和KEYOPT命令–耦合场分析…顺序耦合方法INTRODUCTIONTOANSYS11.0TrainingManualb)定义结构的材料性质(EX等)，包括热膨胀系数(ALPX)(若使用的是ANSYS提供的材料库，材料的热特性和结构特性均已定义，该项可以省略).注意：如果没有定义ALPX，或将该项设置为零，则不计算热应变。可以用该项技巧“关闭”温度的影响。c)指定静力分析类型，这一步只在热分析是瞬态分析时用。MainMenuSolutionAnalysisTypeNewAnalysis•或用ANTYPE命令–耦合场分析…顺序耦合方法INTRODUCTIONTOANSYS11.0TrainingManuald)施加结构荷载，而把温度作为荷载的一部分。•MainMenuSolutionDefineLoadsApplyStructuralTemperatureFromThermAnaly•或用LDREAD命令。e)求解。f)观察应力结果。–耦合场分析…顺序耦合方法INTRODUCTIONTOANSYS11.0TrainingManual•直接耦合方法，通常只涉及用耦合单元的分析，单元必须包括热、结构的自由度。1.首先用以下耦合单元之一建立模型并划分网格。•PLANE13(平面实体单元)。•SOLID5(六面体单元)。•SOLID98(四面体单元)。–在模型上施加结构荷载、热荷载及约束。–求解并查看热和结构结果。Combined热分析结构分析jobname.rst–耦合场分析直接耦合方法INTRODUCTIONTOANSYS11.0TrainingManual•顺序方法–对非高度非线性耦合情况，顺序方法更有效、灵活，因为它可以独立执行两种分析。–在顺序方法热-应力分析中，例如，在非线性瞬态分析之后可以紧接着进行线性静力分析，然后可以把热分析中任意荷载步或时间点的节点温度作为应力分析的荷载。•直接方法–对耦合场是高度非线性情况，直接方法更好，并且该方法用耦合公式单一求解时是最好的。–直接耦合的例子，包括压电分析，有流体流动的共轭传热分析及电路电磁分析。–耦合场分析顺序耦合方法与直接耦合方法比较TrainingManual带散热片的轴对称管热应力耦合分析---顺序耦合INTRODUCTIONTOANSYS11.0TrainingManual13A.热应力分析–顺序耦合带散热片的轴对称管说明•接上面练习12的轴对称散热片问题做一个热应力分析。如下所示，管内部有压力。顶部的线(在Y=1.0处)代表一个对称边界,我们将该线的所有节点的UY自由度耦合起来。INTRODUCTIONTOANSYS11.0TrainingManual13A.热应力分析–顺序耦合带散热片的轴对称管载荷和边界条件从练习12的热分析获取节点温度施加内压力10,000psi在Y=1的节点处耦合自由度UY在线上施加对称边界条件INTRODUCTIONTOANSYS11.0TrainingManual1.按教师指定的工作目录，用“pipe-th-str”作为作业名，进入ANSYS。2.从练习12中恢复数据库文件(或pipe-th.db1):–UtilityMenuFileResumefrom…•选择“pipe-th.db”,然后按[OK]–或用命令:RESUME,pipe-th,db3.在GUI优选框中选择结构:–MainMenuPreferences•选择“Structural”而不选“Thermal”,然后按[OK]4.改变标题:–UtilityMenuFileChangeTitle...•标题为“2DAXI-SYMMTHERMAL-STRESSANALYSISW/INT.PRESS-ESIZE=0.125”•[OK]13A.热应力分析–顺序耦合带散热片的轴对称管INTRODUCTIONTOANSYS11.0TrainingManual5.删除实体模型边上的对流载荷:–MainMenuPreprocessorLoadsDefineLoadsDeleteAllLoadDataAllSolidModLds–[OK]–或用命令:/PREP7LSCLEAR,SOLID6.将温度单元改为相应的结构单元:–MainMenuPreprocessorElementTypeSwitchElemType•选择“ThermaltoStruc”,然后按[OK]•检查警告信息窗,然后按[Close]–或用命令:ETCHG,TTS13A.热应力分析–顺序耦合带散热片的轴对称管INTRODUCTIONTOANSYS11.0TrainingManual13A.热应力分析–顺序耦合带散热片的轴对称管7.将单元选项设为轴对称:–MainMenuPreprocessorElementTypeAdd/Edit/Delete•[Options...]–将K3设为轴对称,然后按[OK]•[Close]8.在练习12的热应力分析中施加温度荷载:–MainMenuPreprocessorLoadsDefineLoadsApplyStructuralTemperatureFromThermAnaly•选择结果文件“pipe-th.rth”,然后按[OK]9.在线（Y=0）上施加对称边界条件:–MainMenuPreprocessorLoadsDefineLoadsApplyStructuralDisplacementSymmetryB.C.OnLines•选择线3,5,11,然后按[OK]INTRODUCTIONTOANSYS11.0TrainingManual13A.热应力分析–顺序耦合带散热片的轴对称管10.将节点自由度（在Y=1）UY耦合:10a.选择节点（在Y=1）:–UtilityMenuSelectEntities...•通过“ByLocation”选择“Nodes”•选择“Ycoordinates”•将Min,Max设为1,然后按[OK]–或用命令:NSEL,S,LOC,Y,110b.将选择的一组节点在UY方向上耦合起来:–MainMenuPreprocessorCoupling/CeqnCoupleDOFs•[PickAll]•将NSET设为1•将Lab设为UY,然后按[OK]–UtilityMenuSelectEverything–或用命令:CP,1,UY,ALLALLSEL,ALLINTRODUCTIONTOANSYS11.0TrainingManual13A.热应力分析–顺序耦合带散热片的轴对称管11.将内部压力施加在线上:–MainMenuPreprocessorLoadsDefineLoadsApplyStructuralPressureOnLines•拾取线9和13,然后按[OK]•将VALUE设为1000,然后按[OK]–或用命令:SFL,9,PRES,1000SFL,13,PRES,100012.通过显示体荷载检查温度载荷:–UtilityMenuPlotCtrlsSymbols•将体载荷符号设为“Structuraltemps”,然后按[OK]–UtilityMenuPlotElements–或用命令:/PBF,TEMP,,1EPLOTINTRODUCTIONTOANSYS11.0TrainingManual13A.热应力分析–顺序耦合带散热片的轴对称管13.存储数据库并计算结果:–Pickthe“SAVE_DB”buttonintheToolbar(orselect:UtilityMenuFileSaveasJobname.db)–MainMenuSolutionSolveCurrentLS•观察“/STATUSCommand”窗口并关闭•[OK]•[Close]–在求解结束时关闭黄色信息框–或用命令::SAVE/SOLUSOLVE14.进入后处理并观察结果:–MainMenuGeneralPostproc–或用命令:/POST1INTRODUCTIONTOANSYS11.0TrainingManual13A.热应力分析–顺序耦合带散热片的轴对称管14a.显示位移:–MainMenuGeneralPostprocPlotResultsContourPlotNodalSolu•拾取“DOFsolution”和“TranslationUSUM”,选择“Def+undefedge”,然后按[OK]–或用命令:PLNSOL,U,SUM,2,1INTRODUCTIONTOANSYS11.0TrainingManual13A.热应力分析–顺序耦合带散热片的轴对称管14b.画vonMises应力:–MainMenuGeneralPostprocPlotResultsContourPlotNodalSolu•拾取“Stress”和“vonMisesSEQV”,选择“Defshapeonly”,然后按[OK]–或用命令:PLNSOL,S,EQVINTRODUCTIONTOANSYS11.0TrainingManual13A.热应力分析–顺序耦合带散热片的轴对称管14c.画径向应力等值线图:–MainMenuGeneralPostprocPlotResultsContourPlotNodalSolu•拾取“Stress”和“X-directionSX”,然后按[OK]–或用命令:PLNSOL,S,XINTRODUCTIONTOANSYS11.0TrainingManual13A.热应力分析–顺序耦合带散热片的轴对称管14d.绕Y轴将轴对称的径向应力扩展90度，然后关于x-z平面做镜面反射:–UtilityMenuPlotCtrlsStyleSymmet