针对各向异性材料铸造残余应力的数值模拟方法技术

本发明专利技术公开了一种针对各向异性材料铸造残余应力的数值模拟方法，用于解决现有铸造残余应力的模拟方法适用范围窄的技术问题。技术方案是利用ProCAST软件对铸件凝固及冷却过程温度场进行数值模拟，并将温度场转换成ABAQUS能接受的温度幅值，生成ABAQUS数据文件INP文件；将经过转换的温度幅值复制至生成的ABAQUS数据文件INP文件中，生成新的INP文件；对于生成新的INP文件，利用ABAQUS Command命令完成对各向异性材料铸件的残余应力数值模拟。该方法结合了ProCAST软件以及ABAQUS软件各自的优点，能够准确地模拟各向异性材料铸件在铸造过程中产生的残余应力，准确率达到81％以上。

【技术实现步骤摘要】
针对各向异性材料铸造残余应力的数值模拟方法
本专利技术涉及一种铸造残余应力的模拟方法，特别涉及一种针对各向异性材料铸造残余应力的数值模拟方法。
技术介绍
各向异性材料在现代制造工业中应用日益广泛，这类材料由于在各微观方向上的热物性性能及力学性能不一致，导致该材料在铸造过程中在各方向上的热量传递及收缩变形均会不一致，从而使铸件的铸造应力分布很难确定，例如航空发动机中单晶涡轮叶片的铸造残余应力分布问题是研制高推重比发动机的一个重要障碍。目前针对铸造残余应力的数值模拟主流软件是ProCAST，但这种软件的缺点是材料库不全，在应力场模拟方面只能针对各向同性材料铸件进行数值模拟。文献“FEpredictionofresidualstressesofinvestmentcastinginaBottomCoreVaneunderequiaxedcooling，JournalofManufacturingProcesses,Vol13(1),January2011,Pages30–40”给出了一种利用ABAQUS软件来模拟铸造残余应力的新方法。此方法将ProCAST对导向叶片的温度场编写成ABAQUS能接受的UTEMP子程序，将该温度场对导向叶片的每个网格节点进行逐一加载，进行顺序热力耦合，以求解导向叶片的铸造残余应力，所得结果与实际测量吻合度较好。但是该方法并不适用于模拟各向异性材料的铸造残余应力，由于每种各向异性材料都有其特殊性，为了准确的表达其特性，一般根据其特性来编写用户材料UMAT子程序，而在ABAQUS中每个仿真文件只能附加一个子程序，因此上述方法不能同时对温度场UTEMP和材料子程序UMAT进行数值模拟，具有较大的局限性。此外，上述方法对温度场编写成UTEMP子程序过程非常繁琐，不具备通用批量化处理的能力，也限制了其使用范围。综上所述，现阶段针对铸造残余应力的方法并不适用于各向异性材料，适用范围窄。
技术实现思路
为了克服现有铸造残余应力的模拟方法适用范围窄的不足，本专利技术提供一种针对各向异性材料铸造残余应力的数值模拟方法。该方法利用ProCAST软件对铸件凝固及冷却过程温度场进行数值模拟，并将温度场转换成ABAQUS能接受的温度幅值，在Viewer-CAST中将浇注系统中除炉体外各组件分别以INP文件导出，生成ABAQUS数据文件INP文件；将经过转换的温度幅值复制至生成的ABAQUS数据文件INP文件中，生成新的INP文件；根据铸件材料的特性利用Fortran语言编写各向异性材料子程序UMAT文件；对于生成新的INP文件，结合子程序UMAT文件，利用ABAQUSCommand命令完成对各向异性材料铸件的残余应力数值模拟。该方法结合了ProCAST软件能较准确模拟铸件温度场以及ABAQUS具有丰富材料库两种优点，利用用户材料UMAT子程序以及高可靠性的温度场转换接口，能够准确地模拟各向异性材料铸件在铸造过程中产生的残余应力，准确率达到81％以上。本专利技术解决其技术问题所采用的技术方案：一种针对各向异性材料铸造残余应力的数值模拟方法，其特点是包括以下步骤：步骤一、利用ProCAST软件对铸件凝固及冷却过程温度场进行数值模拟，并利用ProCAST软件的后处理软件Viewer-CAST将铸件温度场导出；步骤二、在MATLAB操作环境中，利用编写好的温度场转换接口，将步骤一所得温度场转换成ABAQUS能接受的温度幅值，并利用转换接口生成铸件的节点集及预定义场；步骤三、在Viewer-CAST中将浇注系统中铸件、模壳组件分别以INP文件导出，并在ABAQUS中重新装配，加入边界条件，修改有限元网格类型，生成ABAQUS数据文件INP文件；步骤四、将步骤二所得转换好的温度幅值、节点集及预定义场复制至步骤三所生成的ABAQUS数据文件INP文件中，生成新的INP文件；步骤五、根据铸件材料的特性利用Fortran语言编写各向异性材料子程序UMAT文件；步骤六、对于步骤四生成新的INP文件，结合步骤五所得UMAT子程序，利用ABAQUSCommand命令完成对各向异性材料铸件的残余应力数值模拟。本专利技术的有益效果是：该方法利用ProCAST软件对铸件凝固及冷却过程温度场进行数值模拟，并将温度场转换成ABAQUS能接受的温度幅值，在Viewer-CAST中将浇注系统中除炉体外各组件分别以INP文件导出，生成ABAQUS数据文件INP文件；将经过转换的温度幅值复制至生成的ABAQUS数据文件INP文件中，生成新的INP文件；根据铸件材料的特性利用Fortran语言编写各向异性材料子程序UMAT文件；对于生成新的INP文件，结合子程序UMAT文件，利用ABAQUSCommand命令完成对各向异性材料铸件的残余应力数值模拟。该方法结合了ProCAST软件能较准确模拟铸件温度场以及ABAQUS具有丰富材料库两种优点，利用用户材料UMAT子程序以及高可靠性的温度场转换接口，能够准确地模拟各向异性材料铸件在铸造过程中产生的残余应力，准确率达到81％以上。下面结合附图和具体实施方式对本专利技术作详细说明。附图说明图1是本专利技术针对各向异性材料铸造残余应力的数值模拟方法的流程图。图2是本专利技术方法实施例中工字梁模型图。图3是本专利技术方法实施例中温度场转换主界面图。具体实施方式参照图1-3。本专利技术针对各向异性材料铸造残余应力的数值模拟方法具体步骤如下：第一步、首先利用三维建模软件对整个浇注系统进行建模，包括铸件和型芯等模型；利用Hypermesh软件对有限元模型进行前处理，采用四面体单元对各组件进行二维网格划分，对铸件关键部位或是结构突变处的网格进行局部加密处理；在ProCAST前处理软件MeshCAST中将各部件进行装配、网格质量检查、生成模壳及体网格划分等，在PreCAST中设置铸件材料热物性及及力学参数、换热系数及换热面设置、固定及对称面设置等边界条件设置，完成对整个浇注过程进行仿真；在ProCAST后处理软件Viewer-CAST中将铸件温度场以*.ntl格式导出，除炉体外，将铸件、模壳以及其他组件分别以*.inp格式导出。第二步、利用ProCAST与ABAQUS温度场转换接口对第一步中所导出的温度场进行转换，得到ABAQUS温度幅值，并生成节点集及预定义场。第三步、将第一步中所导出的铸件、模壳以及其他组件的*.inp文件分别导入至ABAQUS中进行重新装配，设置静力稳态分析步，除换热条件外，其余边界条件设置均与ProCAST保持一致，将得到的CAE文件以*.inp格式导出。第四步、将第二步中ABAQUS温度幅值、生成节点集及预定义场分别复制至第三步所得的*.inp文件的相应位置。第五步、利用Fortran语言针对铸件各向异性材料编写ABAQUS用户材料UMAT子程序*.for文件，并在Visual-studio中调试成功。第六步、利用ABAQUSCommand命令将第四步得到的*.inp文件与第五步得到的*.for文件联合起来对铸件应力场进行数值模拟。应用实施例：本专利技术利用求解工字梁残余应力作为一个具体实施例，该工字梁材料为各向异性材料DD6，铸造方式为定向凝固，下面结合该实施例对本专利技术细述如下：第一步、利用UG软件对工字梁浇注系统进行三本文档来自技高网...

【技术保护点】
一种针对各向异性材料铸造残余应力的数值模拟方法，其特征在于包括以下步骤：步骤一、利用ProCAST软件对铸件凝固及冷却过程温度场进行数值模拟，并利用ProCAST软件的后处理软件Viewer‑CAST将铸件温度场导出；步骤二、在MATLAB操作环境中，利用编写好的温度场转换接口，将步骤一所得温度场转换成ABAQUS能接受的温度幅值，并利用转换接口生成铸件的节点集及预定义场；步骤三、在Viewer‑CAST中将浇注系统中铸件、模壳组件分别以INP文件导出，并在ABAQUS中重新装配，加入边界条件，修改有限元网格类型，生成ABAQUS数据文件INP文件；步骤四、将步骤二所得转换好的温度幅值、节点集及预定义场复制至步骤三所生成的ABAQUS数据文件INP文件中，生成新的INP文件；步骤五、根据铸件材料的特性利用Fortran语言编写各向异性材料子程序UMAT文件；步骤六、对于步骤四生成新的INP文件，结合步骤五所得UMAT子程序，利用ABAQUS Command命令完成对各向异性材料铸件的残余应力数值模拟。

【技术特征摘要】
1.一种针对各向异性材料铸造残余应力的数值模拟方法，其特征在于包括以下步骤：步骤一、利用ProCAST软件对铸件凝固及冷却过程温度场进行数值模拟，并利用ProCAST软件的后处理软件Viewer-CAST将铸件温度场导出；步骤二、在MATLAB操作环境中，利用编写好的温度场转换接口，将步骤一所得温度场转换成ABAQUS能接受的温度幅值，并利用转换接口生成铸件的节点集及预定义场；步骤三、在Viewer-CAST中将浇注系统中铸件、模壳组件分别以INP文件导...

【专利技术属性】
技术研发人员：卜昆，胡杰，邱飞，田国良，田嘉伟，王志红，金宗李，李扬，赵丹青，廖金明，张明，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西;61