焊接数值模拟

1、哈尔滨工业大学材料加工过程数值模拟基础实验课程25号钢平板对接焊温度场数值模拟实验报告姓名学号班级王保全14SD09004焊接一班材料科学与工程学院ANSYS及其应用ANSYS软件是美国ANSYS公司研制的大型通用有限元分析（FEA）软件，能够进行包括结构、热、声、流体、电磁热场等科学的研究。在核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、造船、生物医学、轻工、地矿、水利、日用家电等领域有着广泛的应用。ANSYS的功能ANSYS的功能如下： 建立计算模型或者输入结构、产品、组件或系统的CAD模型。 应用施加载荷或者其他设计条件。 研究模型的物理相应，如应

2、力水平、温度分布或者电磁场等。 对产品进行优化设计，以降低产品的费用。 做数值模拟实验。ANSYS包括100多个单元，提供了对各种物理场的分析功能，可以将其应用于如下学科：结构分析、热分析、高度非线性瞬态动力分析（ANSYS/LS-DYNA）、流体静力学和动力学分析（FLOTRAN）、电磁场分析、声学分析、压电分析、多长耦合分析、设计灵敏度及优化分析。ANSYS的设计优化功能允许优化任何方面的设计变量和约束变量，如形状、应力、自然频率、质量、费用、温度、次势、压力、速度或离散量等，可进行参数、形状、拓扑优化。ANSYS的特点ANSYS作为一个功能强大、应用广泛的有限元分析软件，其技术特点主要表

3、现在以下几个方面： 能实现多场及多场耦合分析。用户不但可用其进行诸如结构、热、流体流动、电磁等的单独研究，还可以进行这些问题的相互影响研究。例如：热-结构耦合，磁-结构耦合以及电-磁-流体-热耦合等。 实现前后处理、求解和多场分析统一数据库的一体化大型FEA软件。 具有多物理场优化功能的FEA软件。 强大的并行计算功能支持分布式并行和共享内存式并行。 多种自动网格划分技术。 良好的用户开发环境。ANAYS综合应用菜单、对话框、工具条、命令并行输入、图形化输出等多种方式，从而使应用更加方便。 方便的二次开发功能。应用宏、参数设计语言、用户可编程特性、用户自定义界面语言、外部命令等功能，可以开发个

4、性化的应用程序。ANSYS的结构ANSYS主要包括三个部分：前处理模块、分析计算模块和后处理模块。前处理模块提供了一个强大的实体建模及网格划分工具，用户可以方便的构建有限元模型；分析计算模块包括结构分析（可进行线性分析、非线性分析和高度非线性分析）、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦和分析，可模拟多种物理介质的的相互作用，具有灵敏度分析及优化分析能力；后处理模块可将计算结果以彩色等值线显示、梯度显示、矢量显示、粒子流迹显示、立体切片显示、透明及半透明显示等图形方式显示出来，也可将计算结果以图表、曲线形式显示或输出。ANSYS有限元分析下面以平板对接焊的温度场分析为例

5、详细解析ANSYS有限元分析。焊接材料为25号钢，平板尺寸为120mmX60mmX6mm，施加热源为高斯热源，热源的有效半径为6mm。焊接工艺参数如表1所示。表1焊接工艺参数工艺参数电压/V电流/A焊接速度/(m/s)焊接效率/(%)数值251800.0175各种有限元软件都是由前处理、计算过程和后处理三个模块组成，但具体到每个软件其具体界面以及操作方式并不一样。下面详细解析利用ANSYS分析该问题的操作过程。(1)定义分析路径点击ANSYSProductLauncher，如图1所示。在SimulationEnvironment中选ANSYS，在License中选ANSYSMultiphysi

6、cs，定义ANSYS的分析环境。然后定义工作目录，在WorkingDirectory中通过点击Browse定义分析，然后定义文件名，可在JobName中更改文件名，本例中选用file为文件名。然后点击run进行分析。从实菜单选择File|Changetitle命令，将弹出Changetitle对话框，在enternewtitle中输入标题，点击OK按钮完成标题的定义，如图2所示。图2定义标题3）选择菜单过滤参数从主菜单中选择Preference命令，将弹出PreferenceforGUIFiltering（菜单过滤参数选择）对话框，如图3所示。选中Thermal复选框，点击OK按钮确定。/Pr

7、eferencesforGUTFilterinqKEYWPrsfertncforGUIFilteringIndividualdisciplinefs）toshowintheGUI厂呂电u匚tural再rhermal厂ANSY5-FluidElectromagne-ti口0Magnetic-INqdal厂Magnetic-Edge尸HighFrr-que-nEyJEle-cfricNotKIfnoindividualdisciplinesarete-lectedthwilIaIshew.Disciplineoptonef*hrMethodCi風|匚hue|Hdp图3菜单过滤参数选择4）定义单元类

8、型进行有限元分析时，首先应根据问题的几何结构、分析类型和所分析问题的精度要求等，选定适合具体分析的单元类型。本例中选用热分析单元类型S0LID70,为了划分网格的需要还选用单元类型Mesh200。Mesh200单元类型划分网格生成的单元和节点在求解中是无效的，即不会对这些单元和节点求解。从主菜单中选择preprocessor|elementtype|add/edit/delete,将弹出elementtype（单元类型）对话框。点击add按钮，将弹出Libraryofelementtype（单元类型库）对话框，如图4所示。然后在左边的列表中选择Solid选项，选择实体单元类型。在右边的列表中选

9、择Brick8node70选项，选择把节点实体单元类型Solid70。Elementtypereferencenumber输入1，即定义单元类型号为1。LibraryofElementTypes图4单元类型对话框ALibraryofElementTypesElsmenc&ypsreferencenumberK点击OK按钮，将Solid70单元添加，如图5所示。然后再点击Add按钮添加Mesh200单元类型。再次进入单元类型库，如图6所示。图12各项同性材料的热传导定义图6定义Mesh200单元类型Solved，在右边的列表中选择Mesh200单元类型。图5单元类型在左边的列表中选择Not点击O

10、K按钮进入单元类型对话框，如图7所示。选择Mesh200单元类型，点击Options按钮，进入MESH200elementtypeoptions，如图8所示。在Elenmentshapeand#ofnodeKI中选择QUAD4-NODE,即定义Mesh200为平面四节点单元类型，点击OK即可。返回ElementType对话框，然后点击Close按钮，关闭单元类型对话框。MESH200elefnenttypeop:ionsSjtionsforMESHZM.ElementTypeRef.No.2ElefTiertshapeand#ofnodsK1Eemr:besting图7单元类型OKCanedH

11、elp图8定义Mesh200单元类型5)定义单位在本例的有限元分析工程中，分析过程中选用国际单位制，点击主菜单钮，完成对单位的定义。Preprocessor|Materialsprops|Selectunits,弹出对话框，如图9所示，选择SI(MKS)，定义国际单位制，然后点击OK按钮确定。定义角度单位为度，点击菜单Parameters|Angularunits，弹出对话框，如图10所示，选择角度单元DegreeDEG，然后点击OK按ASetectFilteringUniMforMaterialLibuarybUnirESttecrittLNfTS绸i曲naidiirnpouFinodd1冋

12、曲Te站rvuribet吕emed”略曲jcfier用曾环NOcoiwewiaIs3亡射百蹴也be诃出1AngularUnitsforParametricFunctionsPAPJMLrih;margular-saramEticfunncmHMpI（6）定义材料参数从主菜单中选则Preprocessor|Materialprops|Materialmodels命令，将弹出DefineMaterialModelBehavior（定义材料模型）属性窗口，如图11所示。丄IGraphHdp图11定义材料模型属性窗口依次双击Thernal|Conductivity|lsotropic，展开材料属性的树

13、形结构，将弹出材料热传导定义对话框。因为本例中要输入不同温度下的材料参数，所以点击AddTemperature按钮，增加不同温度下的热传导数值，如图12所示。在KXX文本框中输入不同温度下的热传导参数后，点击OK按钮，返回定义材料模型属性窗口。ConductivityforMaterialNumber1：!-500：.i|iii.rii1-II2500474077|i35140142CdnductiTity(Iedtrcipic)forJfat&rialNunber1T2T3T4T5T6T7T8TerojeraturesXXXAddTemperatureDeleteTemperatureMXT

14、E.XA1.SKXXKMXPREVIEWANSYSR15.0DEC3dZOlLSI06DBOOloco1BOO200026002：307BQ17SCI口TEME-CanEel依次双击ThermalSpecificheat,展开材料属性的树形结构，将弹出材料比热定义对话框。因为本例中要输入不同温度下的材料参数，所以点击AddTemperature按钮,增加不同温度下的热传导数值,如图13所示。在C文本框中输入不同温度下的比热参数后，点击OK按钮，返回定义材料模型属性窗口。依次双击Thermal|Density，展开材料属性的树形结构，将弹出材料密度定义对话框。因为本例中要输入不同温度下的材料参数

15、，所以点击AddTemperature按钮，增加不同温度下的密度数值，如图14所示。在DENS文本框中输入不同温度下的密度参数后，点击OK按钮，返回定义材料模型属性窗口。250EOO7501O0DII1EOOI1TE0BOO|4SQ63Q&70|560|7SO|s20图13各向同性材料的比热定义DensityforMaterialNumber1TS二T?T6TrTenperatures250500福01X150017X郊7700&N5DEJJ5PREVXHWANSYSR15.0DEC30ED20L27257640BETTSTEMPTTtSJ76(1V62CIT-a-fla7s-ea7260mc

16、i712(170401260图14各向同性材料的密度定义在DefineMaterialModelBehavior窗口中，从菜单中选择Material|Exit命令，或者单击窗口右上角的X按钮，完成材料模型的定义。（7）建立实体模型在建立模型的过程中，应考虑到后面的有限元网格划分。而且在焊接模拟中，焊缝及其附近的网格采用较密的网格划分，而离焊缝较远的位置采用较疏的网格首先定义模型的标量参数。单击菜单Parameter|ScalarParameter,在Selction中输入LENGTH=0.12，点击Acecept按钮，完成对长度参数的定义。然后依次输入WIDTH=0.06，HEIGHT=0.0

17、06，RB=0.002，U=25，I=180，V=0.01,EFF=0.75,R=0.006,q二U*1*EFF,x_center=O,y_center=0,如图15所示。最后点击Close按钮，完成对各个参数的定义。建立模型，点击主菜单Preprocessor|Modeling|Volumes|Block|Bydimension，弹出创建实体对话框，在其中输入各个尺寸参数，如图16所示，建立模型如图17所示。点击主菜单Preprocessor|Modeling|Volumes|Block|Bydimension建立两个过渡实体，如图18所示，命令流为：block,0,2*rb,3*rb,5*

18、rb,0,-heightblock,2*rb,4*rb,3*rb,5*rb,0,-height图15定义标量参数图16定义实体对话框图17有限元模型ANSVSR&0DK卄n：lriMt3?ZIL4LlrHULANSYSE15.D图18有限元过渡实体模型图19有限元网格划分模型指定各个过渡实体的线的网格划分。然后通过APDL语言声称其他的实体如图19所示。命令如下：vsel,s,2,3!此实体为所建立的过渡实体cm,vl,voluallsel,all*do,i,1,14,1vgen,2,v1,i*4*rb,1,0,*enddo依照上面的方法，建立随后的实体模型，如图20所示。点击主菜单Prepr

19、ocessor|Numberingctrls|Mergeitems，弹出关键点压缩菜单，在Label选择关键点，然后点击OK按钮完成压缩，如图21所示。PCfljipaiTbiiMergeColncldwcexEqpubdlentllyDefinedAN5YSR1S-IIU.CIII-3-111l*13-1QTCXJEFoSdmoddiotr-arwe-yriLlHMribmior-Mkct?|TJUMMftG|Me*geBktlde曲orEq卫出即!4,DefinedHehLaba!TypeilmibobrrMrgvTDLERRmgr-nci4ncGWTTChRetainI充HqhesP|l

20、OWe5inimbis*t-arKsl图20实体模型图21关键点压缩点击Moddeling|Booleans|Glue|Volumes，弹出GlueVolumes对话框，如图22所示，点击Pickall，对所建立的模型进行布尔运算，即把所有的模型粘在一起。（8）网格划分对焊缝实体进行网格划分，首先选择焊缝区域的实体，点击实菜单Select|Entities，选择实体Volumes-Bylocation-Ycoordinate，分别在最小值、最大值文本框中输入最大值和最小值0,3*rb，点击OK按钮完成对实体的选择，如图23所示。GlueVolumes仔口丄U1C厂UHlZPlUSX：nsr丄a

21、厂Eqjl广C&1广U-uul.a-u.n_tOSTbLUTl=dmXZbLUTl=二Vo-1LLI-Tca=活工丄：3匸&工匸GirgLHxlpMnjl.Inc：/VSelectEntities|VolumesIByl_ocstionOK口七PickA1J-OKRlot匚劉n匸图22体的布尔运算图23选择实体然后对实体进行网格划分。随后对过渡单元进行网格划分。依次对其他的实体进行网格划分，最后实体的网格划分的情况如图24所示。p-1ftnrurici1ngANSYSR15.0图24模型的网格划分（9）有限元求解点击主菜单Solution|Analysistype|Newanalysis,弹出

22、分析类型对话框，如图25所示。选择Transistent复选框，点击OK按钮，弹出瞬态分析对话框，如图26所示，点击OK即可。打开自动时间步长，点击Solution|Loadstepopts|Time/frequent|Time-timestep,弹出对话框，如图27所示，在AUT0TS中选择ON,然后点击OK按钮。定义热源函数，点击菜单Parameter|Function|Define/Edit,弹出FunctionEditor对话框，如图28所示。在Result文本框中输入高斯热源的函数表达式，然后点击File|Save保存文件。随后点击Parameter|Functions|Readfr

23、omfile，选取刚刚命名的函数表达式。在Tableparametername中命名arc，在ConstantValues中输入各个参数，点击OK按钮，完成函数的定义。NewAnalysisANTVPECSteady-GtaireTr旨nsient厂Subs-trurturingCancelHelp图25分析类型图26瞬态分析选项TimeandTinneSeepOptionsTimesndTimeSt&pOptionsTIMETime-Atendoflod-oDELTI：MTimetepsizeKBC5icpp?dorrampodb.c.金RampedL5tflppedAUTOTSAuSomatictimesteppingnON厂OFFLProgChos-snD-ELTTMMinimumti弓trpsInlMsKimumtime34zv-U-sepireviou-sstepsizt?VY巴舌TSRTime-stepresetba-sedonspe-ciF匚time-pointsTirriepo-iHtsfro-mifNoireselLEsistingarrayLarmyNote;T5RE5command13-validforthc-rmalpJemcntB.,thcrm-sl-