利用ANSYS 模拟大坝的运行期应力变形
利用ANSYS模拟大坝的运行期应力变形1.主要步骤：
S1：在Auto CAD里面按照给出尺寸绘图并导入到ANSYS软件中。
S2：通过命令流或者GUI操作来设置材料并给其赋值。
S3：通过APDL命令或者建模的“Operate”里面的“Extrude”沿Z轴将面拉成体，注意距离为10。
S4：按照题目要求，通过APDL或者“Meshing”里面的“Mesh Attributes”给体赋值材料号。
S5：划分网格，先将线分段再划分网格，注意在给侧面划分网格时划分份数为2，并将“KYNDIV SIZE,NDIV can be changed”勾选去掉，否则划分网格时划分份数会有变化。
S6：通过APDL或者“Solution”里面的“Define Loads”添加边界条件和自重。
S7：模拟初始状态，“杀死”坝体单元，求解。
S8：模拟完建状态，“激活”坝体状态，求解。
S9：运行期，施加水荷载，计算扬压力，求解运行期。
S10：创建荷载工况，查看结果。
2.导出结果：2.1初始状态模拟
2.2完建状态模拟
2.3运行期模拟
2.4运行期的变形图
2.5运行期的第一应力图
2.6运行期的第一应力等值线图
2.7运行期总位移变形图
2.8运行期总位移变形等值线图
ANSYS计算温度场及应力场
基于ANSYS有限元软件实现施工温控仿真的主要技术（1）研究方法和分析流程 本次计算利用ANSYS软件来进行象鼻岭碾压混凝土拱坝全过程温控仿真计算分析。具体分析流程如下： 1）收集资料：包括工程气象水文资料、大坝体型、热力学参数、工程进度、施工措施、防洪度汛和蓄水等。 2）整理分析资料：参数拟合、分析建模方法。 3）建模：采用ANSYS软件进行建模，划分网格。 4）编写计算批处理程序：根据资料结合模型编写计算温度场的ANSYS批处理程序。 5）检查计算批处理程序：首先检查语句，然后导入计算模型检查所加荷载效果。 6）计算温度：使用ANSYS软件温度计算模块进行计算。 7）分析温度结果：主要分析各时刻的温度场分布和典型温度特征值。 8）应力计算建模：模型结构尺寸与温度分析模型相同，需要改变把温度分析材料参数改为应力分析材料参数。 9）计算应力：使用ANSYS软件温度应力计算模块和自编的二次开发软件进行计算。 10）分析应力结果：主要分析应力场分布和典型应力特征值。 11）编写报告：对计算流程和结果实施进行提炼总结，提出可行的温控指标和措施。 （2）前处理 1）建模方法选择。 有限元建模一般有两种方法：一种为通过点线面几何拓扑的方法建模，这种建模方法精确，但是比较费时。对于较大规模的建模任务花费时间太多。另一种为通过其他软件导入，如CAD，通过在其他软件中建模，然后输出为ANSYS 可以识别的文件类型，再导入ANSYS中完成建模过程，这种建模方式精度较直接建模的精度要稍低一些，但是由于要求建模的模型已经在CAD软件中完成了
初步建模，可以直接拿来稍作处理即可应用，时间花费较少。本计算选用从CAD 软件导入ANSYS中来建立模型。 2）建模范围。 建模范围可以分为全坝段建模和单坝段建模，全坝段建模可以全面反映整个坝体的温度和应力情况，但是建模难度高、计算量大；单坝段建模建模难度小，计算量也相对较小，一般情况下单坝段建模即可满足要求。 3）施工模拟层厚。 根据已建碾压混凝土坝经验，碾压层厚一般为0.3m左右，按照0.3m一层建模是最精确的，但是如果按照0.3m一层建模，计算网格数量巨大，计算时间长，对于硬件要求较高，在硬件和时间达不到要求的情况下，按照3m一层以下精度都是可以基本满足要求的。 4）分区模拟。 由于各分区混凝土水化热差别较大，对于温度计算影响较大，因此建模要尽量反映混凝土大坝内部分区变化。基岩由于对混凝土只是导热作用，且影响范围在10m左右，因此在计算时可以认为是均质体，计算热力学参数采用靠近建基面的地层参数。 5）参数选取。 参数一般选择可研阶段的材料试验报告，如果项目部未能提供这些资料，可以在征求同意的前提下，通过查阅相关书籍，尽量采取相似工程的资料。 （3）计算 1）ANSYS计算模块。 ANSYS计算温度场模块由其自带，可以直接进入模块计算。 2）化学产热模拟。 通过ANSYS中产热命令BFE模拟。 3）边界条件模拟。 ①对流边界条件通过命令SFA模拟。 ②接触散热边界条件通过命令D模拟。 4）浇筑模拟。 通过ANSYS中的生死单元功能实现，初始阶段所有单元均为死单元，死单
ANSYS热应力分析经典例题
ANSYS热应力分析例题 实例1圆简内部热应力分折： 有一无限长圆筒，其核截面结构如图13—1所示，简内壁温度为200℃，外壁温度为20℃，圆筒材料参数如表13．1所示，求圆筒内的温度场、应力场分布。 该问题属于轴对称问题。由于圆筒无限长，忽略圆筒端部的热损失。沿圆筒纵截面取宽度为10M的如图13—2所示的矩形截面作为几何模型。在求解过程中采用间接求解法和直接求解法两种方法进行求解。间接法是先选择热分析单元，对圆筒进行热分析，然后将热分析单元转化为相应的结构单元，对圆筒进行结构分析；直接法是采用热应力藕合单元，对圆筒进行热力藕合分析。 /filname,exercise1-jianjie /title,thermal stresses in a long /prep7 $Et,1,plane55 Keyopt,1,3,1 $Mp,kxx,1,70 Rectng,0.1,0.15,0,0.01 $Lsel,s,,,1,3,2 Lesize, all,,,20 $Lsel,s,,,2,4,2 Lesize,all,,,5 $Amesh,1 $Finish /solu $Antype,static Lsel,s,,,4 $Nsll,s,1 $d,all,temp,200 lsel,s,,,2 $nsll,s,1 $d,all,temp,20 allsel $outpr,basic,all solve $finish /post1 $Set,last /plopts,info,on Plnsol,temp $Finish /prep7 $Etchg,tts Keyopt,1,3,1 $Keyopt,1,6,1 Mp,ex,1,220e9 $Mp,alpx,,1,3e-6 $Mp,prxy,1,0.28 Lsel,s,,,4 $Nsll,s,1 $Cp,8,ux,all Lsel,s,,,2 $Nsll,s,1 $Cp,9,ux,all Allsel $Finish /solu $Antype,static D,all,uy,0 $Ldread,temp,,,,,,rth Allsel $Solve $Finish /post1 /title,radial stress contours Plnsol,s,x /title,axial stress contours Plnsol,s,y /title,circular stress contours Plnsol,s,z /title,equvialent stress contours Plnsol,s,eqv $finish
ANSYS热应力分析实例
ANSYS热应力分析实例 当一个结构加热或冷却时，会发生膨胀或收缩。如果结构各部分之间膨胀收缩程度不同，和结构的膨胀、收缩受到限制，就会产生热应力。 7.1热应力分析的分类 ANSYS提供三种进行热应力分析的方法： 在结构应力分析中直接定义节点的温度。如果所以节点的温度已知，贝U可以 通过命令直接定义节点温度。节点温度在应力分析中作为体载荷，而不是节点自由度 间接法。首先进行热分析，然后将求得的节点温度作为体载荷施加在结构应力分析中。 直接法。使用具有温度和位移自由度的耦合单元，同时得到热分析和结构应力分析的结果。 如果节点温度已知，适合第一种方法。但节点温度一般是不知道的。对于大多数问题，推荐使用第二种方法一间接法。因为这种方法可以使用所有热分析的功能和结构分析的功能。如果热分析是瞬态的，只需要找出温度梯度最大的时间点，并将此时间点的节点温度作为荷载施加到结构应力分析中去。如果热和结构的耦合是双向的，即热分析影响结构应力分析，同时结构变形又会影响热分析（如大变形、接触等），则可以使用第三种直接法一使用耦合单元。此外只有第三种方法可以考虑其他分析领域（电磁、流体等）对热和结构的影响。 7.2间接法进行热应力分析的步骤 首先进行热分析。可以使用热分析的所有功能，包括传导、对流、辐射和表面效应单元等，进行稳态或瞬态热分析。但要注意划分单元时要充分考虑结构分析的要求。例如，在有可能有应力集中的地方的网格要密一些。如果进行瞬态分析，在后处理中要找出热梯度最大的时间点或载荷步。 表7-1热单元及相应的结构单元
重新进入前处理，将热单元转换为相应的结构单元，表7-1是热单元与结构 单元的对应表。可以使用菜单进行转换： Mai n Menu>Prep roeessor>Eleme nt Typ e>Switeh Eleme nt Type ，选择Thermal to Struetual 。 但要注意设定相应的单元选项。例如热单元的轴对称不能自动转换到结构单元中，需要手工设置一下。在命令流中，可将原热单元的编号重新定义为结构单元，并设置相应的单元选项。 设置结构分析中的材料属性（包括热膨胀系数）以及前处理细节，如节点耦 合、约束方程等。 读入热分析中的节点温度, GUI: Solution>Load Apply>Temperature>From Thermal Analysis 。输入或选择热分析的结果文件名*.rth。如果热分析是瞬态的，则还需要输入热梯度最大时的时间点或载荷步。节点温度是作为体载荷施加的，可通过Utility Men u>List>Load>Body Load>On all nodes 列表输出。 设置参考温度，Mai n Men u>Solutio n>Load Setti ng>Refere nee Temp 。 进行求解、后处理。 7.3间接法热应力分析实例 7.3.1 问题描述 图7-1冷却栅示意图
ANSYS热应力分析命令流
/FILNAME,Double,1 ！定义工作文件名。 /TITLE,Temperature Analysis ！定义工作标题。 !* /PREP7 !定义单元。 ET,1,SOLID70 !* !定义材料属性。 MPTEMP,,,,,,,, !定义材料1。 MPTEMP,1,0 MPDATA,KXX,1,,238*3.6 !定义材料1的传热系数KXX1。MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,500 !定义材料1的比热C1。MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,3e-6 !定义材料1密度DENS1。 !* MPTEMP,,,,,,,, !定义材料2。 MPTEMP,1,0 MPDATA,KXX,2,,15*3.6 !定义材料2的传热系数KXX2。MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,2,,100 !定义材料2的比热C2。MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,2.2e-6 !定义材料2密度DENS2。 !* !建立几何模型。 BLC4,-80,-10,160,20,700 VOFFST,3,20, , !* !网格划分。 FLST,5,20,4,ORDE,2 FITEM,5,1 FITEM,5,-20 CM,_Y,LINE LSEL, , , ,P51X
CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1,10, , , , , , ,1 !定义网格大小。 !* TYPE, 1 MAT, 1 REAL, ESYS, 0 SECNUM, CM,_Y,VOLU VSEL, , , , 1 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 TYPE, 1 MAT, 2 REAL, ESYS, 0 SECNUM, CM,_Y,VOLU VSEL, , , , 2 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !定义网格大小完成。 !* FINISH /SOL ANTYPE,4 !定义瞬态分析类型。
ANSYS热应力分析实例教学教材
热流体在代有冷却栅的管道里流动，如图为其轴对称截面图。管道及冷却栅的材料均为不锈钢，导热系数为1.25Btu/hr-in-oF，弹性模量为28E6lb/in2泊松比为0.3。管内压力为1000 lb/in2，管内流体温度为450 oF，对流系数为1 Btu/hr-in2-oF,外界流体温度为70 oF，对流系数为0.25 Btu/hr-in2-oF。求温度及应力分布。 7.3.2菜单操作过程 7.3.2.1设置分析标题 1、选择“Utility Menu>File>Change Title”，输入Indirect thermal-stress Analysis of a cooling fin。 2、选择“Utility Menu>File>Change Filename”，输入PIPE_FIN。 7.3.2.2进入热分析，定义热单元和热材料属性 1、选择“Main Menu>Preprocessor>Element Type>Add/Edit/Delete”，选择PLANE55，设定单元选项为轴对称。 2、设定导热系数：选择“Main Menu>Preprocessor>Material Porps>Ma terial Models”，点击Thermal，Conductivity，Isotropic，输入1.25。 7.3.2.3创建模型 1、创建八个关键点，选择“Main Menu>Preprocessor>Creat>Keypoints>On Active CS”，关键点的坐标如下：
3、设定单元尺寸，并划分网格：“Main Menu>Preprocessor>Meshtool”，设定global size为0.125，选择AREA，Mapped，Mesh，点击Pick all。 7.3.2.4施加荷载 1、选择“Utility Menu>Select>Entities>Nodes>By location>X coordinates，From Full”，输入5，点击OK，选择管内壁节点； 2、在管内壁节点上施加对流边界条件：选择“Main Menu>Solution>Apply>Convection>On nodes”，点击Pick，all，输入对流换热系数1，流体环境温度450。 3、选择“Utility Menu>Select>Entities>Nodes>By location>X coordinates，From Full”， 输入6,12，点击Apply； 4、选择“Utility Menu>Select>Entities>Nodes>By location>Y coordinates，Reselect”，输入0.25,1，点击Apply； 5、选择“Utility Menu>Select>Entities>Nodes>By location>Y coordinates，Also select”，输入12，点击OK； 6、在管外边界上施加对流边界条件：选择“Main Menu>Solution>Apply>Convection>On nodes”，点击Pick，all，输入对流换热系数0.25，流体环境温度70。 7.3.2.5求解 1、选择“Utility Menu>Select>Select Everything”。 2、选择“Main Menu>Solution>Solve Current LS”。 7.3.2.6后处理 1、显示温度分布：选择“Main Menu>General Postproc>Plot Result>Nodal Solution> Temperature”。 7.3.2.7重新进入前处理，改变单元，定义结构材料 1、选择“Main Menu>Preprocessor>Element Type>Switch Elem Type”,选择Thermal to Structure。 2、选择“Main Menu>Preprocessor>Element Type>Add/Edit/Delete”，点击Option，将结构单元设置为轴对称。 3、选择“Main Menu>Preprocessor>Material Porps>Material Models”，输入材料的EX为28E6，PRXY为0.3，ALPX为0.9E-5。
ANSYS热应力分析
题目：ANSYS热应力分析 专业：材料成型及控制工程 班级：型职141 学号：14615118 姓名：武学杰 指导教师：张转转 2017年10月17日至11月13日共4周 指导教师(签字) 系主任(签字)
题目： 第一步：更改文件名
第二步：选择单元 第三步：设置材料属性 1、给定材料的导热系数40W(m·℃)。 Main Menu>Preproessor>Material Props>Material Models 第四步：建立实体模型（国际单位制） 1、创建矩形A1:X1，Y1（0，0）、X2,Y2（0.01，0.07）MainMenu>Preprocessor>Modeling>Creaate>Areas>Rectangle>By Dimensions
2、创建矩形A2：X1,Y1(0,0.05)X2,Y2(0.08,0.07) 3、显示面的编号 Utility Menu>PlotCtrls>Numbering
4、对面A1和A2进行overlap操作 Main Menu>Preocessor>Modeling>Operate>Booleans>Overlap>Areas 第五步：划分网格 1、打开Meshtool对话框； Main Menu>Preprocessor>Meshing>MeshTool 2、设定网格尺寸为0.002，网格形状为四面体映射网格； 3、Mesh。
第六步：施加载荷 1、进入Solution处理器。 Main Menu>Solution 2、设定分析类型为“steady-state” Main Menu>Solution>Define Loads>Apply>Thermal>Temperature>On Lines 3、在外面界线上定义温度载荷60 Main Menu>Solution>Define Loads>Apply>Thermal>Temperature>On Lines 4、在内边界线定义温度载荷0。 5、图形窗口显示线。 Utility Menu>Plot>Lines
ANSYS热应力分析--精选实例.docx
ANSYS 热应力分析实例 当一个结构加热或冷却时，会发生膨胀或收缩。如果结构各部分之间膨胀收 缩程度不同，和结构的膨胀、收缩受到限制，就会产生热应力。 热应力分析的分类 ANSYS提供三种进行热应力分析的方法： 在结构应力分析中直接定义节点的温度。如果所以节点的温度已知，则可以 通过命令直接定义节点温度。节点温度在应力分析中作为体载荷，而不是节点自由度 间接法。首先进行热分析，然后将求得的节点温度作为体载荷施加在结构应 力分析中。 直接法。使用具有温度和位移自由度的耦合单元，同时得到热分析和结构应 力分析的结果。 如果节点温度已知，适合第一种方法。但节点温度一般是不知道的。对于大多数问题，推荐使用第二种方法—间接法。因为这种方法可以使用所有热分析的功能和结构分析的功能。如果热分析是瞬态的，只需要找出温度梯度最大的时间点，并将此时间点的节点温度作为荷载施加到结构应力分析中去。如果热和结构的耦合是双向的，即热分析影响结构应力分析，同时结构变形又会影响热分析（如大变形、接触等），则可以使用第三种直接法—使用耦合单元。此外只有第三种方法可以考虑其他分析领域（电磁、流体等）对热和结构的影响。 间接法进行热应力分析的步骤 首先进行热分析。可以使用热分析的所有功能，包括传导、对流、辐射和表 面效应单元等，进行稳态或瞬态热分析。但要注意划分单元时要充分考虑结构分
析的要求。例如，在有可能有应力集中的地方的网格要密一些。如果进行瞬态分析，在后处理中要找出热梯度最大的时间点或载荷步。 热单元结构单元 LINK32LINK1 LINK33LINK8 PLANE35PLANE2 PLANE55PLANE42 SHELL57SHELL63 PLANE67PLANE42 LINK68LINK8 SOLID79SOLID45 MASS71MASS21 PLANE75PLANE25 PLANE77PLANE82 PLANE78PLANE83 PLANE87PLANE92 PLANE90PLANE95 SHELL157SHELL63 重新进入前处理，将热单元转换为相应的结构单元，表7-1 是热单元与结构单元的对应表。可以使用菜单进行转换：
ANSYS热应力分析经典例题
ANSYS热应力分析例题 实例1——圆简内部热应力分折：有一无限长圆筒，其核截面结构如图13—1所示，简内壁温度为200℃，外壁温度为20℃，圆筒材料参数如表13．1所示，求圆筒内的温度场、应力场分布。 该问题属于轴对称问题。由于圆筒无限长，忽略圆筒端部的热损失。沿圆筒纵截面取宽度为10M的如图1 3—2所示的矩形截面作为几何模型。在求解过程中采用间接求解法和直接求解法两种方法进行求解。间接法是先选择热分析单元，对圆筒进行热分析，然后将热分析单元转化为相应的结构单元，对圆筒进行结构分析；直接法是采用热应力藕合单元，对圆筒进行热力藕合分析。 fini clear /filname,exercise1-jianjie /title,thermal stresses in a long /prep7 $Et,1,plane55 Keyopt,1,3,1 $Mp,kxx,1,70 Rectng,0.1,0.15,0,0.01 $Lsel,s,,,1,3,2 Lesize, all,,,20 $Lsel,s,,,2,4,2 Lesize,all,,,5 $Amesh,1 $Finish /solu $Antype,static Lsel,s,,,4 $Nsll,s,1 $d,all,temp,200 lsel,s,,,2 $nsll,s,1 $d,all,temp,20 allsel $outpr,basic,all solve $finish /post1 $Set,last /plopts,info,on Plnsol,temp $Finish /prep7 $Etchg,tts Keyopt,1,3,1 $Keyopt,1,6,1 Mp,ex,1,220e9 $Mp,alpx,,1,3e-6 $Mp,prxy,1,0.28 Lsel,s,,,4 $Nsll,s,1 $Cp,8,ux,all Lsel,s,,,2 $Nsll,s,1 $Cp,9,ux,all Allsel $Finish /solu $Antype,static D,all,uy,0 $Ldread,temp,,,,,,rth Allsel $Solve $Finish /post1 /title,radial stress contours Plnsol,s,x /title,axial stress contours Plnsol,s,y /title,circular stress contours Plnsol,s,z /title,equvialent stress contours
ANSYS热应力分析实例
A N S Y S热应力分析实例-CAL-FENGHAI-(2020YEAR-YICAI)_JINGBIAN
热流体在代有冷却栅的管道里流动，如图为其轴对称截面图。管道及冷却栅的材料均为不锈钢，导热系数为1.25Btu/hr-in-oF，弹性模量为28E6lb/in2泊松比为0.3。管内压力为1000 lb/in2，管内流体温度为450 oF，对流系数为1 Btu/hr-in2-oF,外界流体温度为70 oF，对流系数为0.25 Btu/hr-in2-oF。求温度及应力分布。 7.3.2菜单操作过程 7.3.2.1设置分析标题 1、选择“Utility Menu>File>Change Title”，输入Indirect thermal-stress Analysis of a cooling fin。 2、选择“Utility Menu>File>Change Filename”，输入PIPE_FIN。 7.3.2.2进入热分析，定义热单元和热材料属性 1、选择“Main Menu>Preprocessor>Element Type>Add/Edit/Delete”，选择PLANE55，设定单元选项为轴对称。 2、设定导热系数：选择“Main Menu>Preprocessor>Material Porps>Material Models”，点击Thermal，Conductivity，Isotropic，输入 1.25。 7.3.2.3创建模型 1、创建八个关键点，选择“Main Menu>Preprocessor>Creat>Keypoints>On Active CS”，关键点的坐标如下：编号 1 2 3 4 5 6 7 8 X 5 6 12 12 6 6 5 5 Y 0 0 0 0.25 0.25 1 1 0.25 2、组成三个面：选择“Main Menu>Preprocessor>Creat>Area>Arbitrary>Throuth Kps”，由1,2,5,8组成面1；由2,3,4,5组成面2；由8,5,6,7组成面3。