word材料成形计算机模拟基于MSC Marc有限元实例仿真 专 业：机械设计制造与其自动化班 级：机 械 1212 班姓 名：黄河学 号：0806111406指导教师：湛利华教师摘 要通过三次课外上机，在教师的指导下练习使用Marc2010，意在熟悉有限元建模、求解与结果分析步骤和方法；能利用Marc软件对实例结构进展静力有限元分析；加深有限元理论关于网格划分概念、划分原如此等的理解 Marc仿真分析过程如下： 〔1〕前处理：有限元分析的最终目的是复原一个实际工程系统的数学行为特征，即分析必须针对一个物理原型准确的数学模型Marc分析的前处理主要就是用来进展建模与网格划分，前处理包括单元类型实常数、材料、属性、建模和划分网格 〔2〕加载求解：有限元模型建好后，就可以进入求解器进展加载求解当施加载荷和边界条件的面、节点或单元比拟多时，应该用实体选择命令把这些对象选出来，然后在其上施加载荷或边界条件，以保证所施加的载荷或边界条件的正确性 〔3〕后处理：Marc的后处理过程即为采集求解器处理分析的结果，提取用户所需的信息，了解计算机结果的过程。
后处理处理器可以用来查看整个模型在某一时间段的计算结果和用来查看模型的某一局部在整个时间段上的计算结果，还可以生成动画，Marc提供的现成的动画制作功能已经非常丰富，各种计算结果的变形动画、时间历程动画、切片动画、粒子轨迹动画等 目 录一、 MSC Marc简介................................................................................3二、 软件应用介绍...................................................................................3三、 实例容...........................................................................................5四、 求解步骤1、网格定义……………………………………………………62、几何定义……………………………………………………63、边界条件……………………………………………………74、定义初始化………………………………………………105、材料性质…………………………………………………106、定义接触…………………………………………………147、定义加载历程……………………………………………158、定义JOBS…………………………………………………169、结果分析显示……………………………………………19五、 总结………………………………………………………… 20六、参考文献.......................................................................................21 一、 MSC Marc简介 MSC.Marc是功能齐全的高级非线性有限元软件的求解器，表现了30年来有限元分析的理论方法和软件实践的完美结合，它具有极强的结构分析能力。
可以处理各种线性和非线性结构分析包括：线性/非线性静力分析、模态分析、简谐响应分析、频谱分析、随机振动分析、动力响应分析、自动的静/动力接触、屈曲/失稳、失效和破坏分析等它提供了丰富的结构单元、连续单元和特殊单元的单元库，几乎每种单元都具有处理大变形几何非线性,材料非线性和包括接触在的边界条件非线性以与组合的高度非线性的超强能力 MARC的结构分析材料库提供了模拟金属、非金属、聚合物、岩土、复合材料等多种线性和非线复杂材料行为的材料模型分析采用具有高数值稳定性、高精度和快速收敛的高度非线性问题求解技术为了进一步提高计算精度和分析效率，MARC软件提供了多种功能强大的加载步长自适应控制技术，自动确定分析曲屈、蠕变、热弹塑性和动力响应的加载步长MARC卓越的网格自适应技术，以多种误差准如此自动调节网格疏密，不仅可提高大型线性结构分析精度，而且能对局部非线性应变集中、移动边界或接触分析提供优化的网格密度，既保证计算精度，同时也使非线性分析的计算效率大大提高 此外，MARC支持全自动二维网格和三维网格重划，用以纠正过渡变形后产生的网格畸变，确保大变形分析的继续进展对非结构的场问题如包含对流、辐射、相变潜热等复杂边界条件的非线性传热问题的温度场，以与流场、电场、磁场，也提供了相应的分析求解能力；并具有模拟流－热－固、土壤渗流、声－结构、耦合电－磁、电－热、电－热－结构以与热－结构等多种耦合场的分析能力。
为了满足高级用户的特殊需要和进展二次开发，MSC.Marc提供了方便的开放式用户环境这些用户子程序入口几乎覆盖了MARC有限元分析的所有环节，从几何建模、网格划分、边界定义、材料选择到分析求解、结果输出、用户都能够访问并修改程序的缺省设置在MSC.Marc软件的原有功能的框架下，用户能够极扩展MARC有限元软件的分析能力二、软件应用介绍有限元分析是对于结构力学分析迅速开展起来的一种现在计算方法它是50年代首先在连续体力学领域—飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应于求解热传导、电磁场、流体力学等连连续行问题想要解答，必须先简化结构，采用数值模拟方法分析〔一〕有限元软件开展特点1.单一场计算向多物理耦合场问题的求解开展2. 由求解线性问题开展到求解非线性问题3. 与CAD/CAM等软件的集成4. 提高自动化的网格处理能力5．软件面向专业用户的开放性6. 软件开发强强联合(二)、分析研究过程1、前处理〔1〕建模有限元分析的最终目的是复原一个实际工程系统的数学行为特征，即分析必须针对一个物理原型准确的数学模型Marc分析的前处理主要就是用来进展建模与网格划分，前处理包括单元类型实常数、材料、属性、建模和划分网格。
2)单元选择有限元模型可分为2D和3D两种,可以由点单元,线单元,面单元或实体单元组成,也可将不同类型的单元混合使用〔3〕网格划分网格划分是其中一个重要的步骤，网格划分的好坏，直接影响到计算的精度和速度，网格划分方法主要有自由网格划分、映射网格划分和体扫掠网格划分三种2、加载求解有限元模型建好后，就可以进入求解器进展加载求解当施加载荷和边界条件的面、节点或单元比拟多时，应该用实体选择命令把这些对象选出来，然后在其上施加载荷或边界条件，以保证所施加的载荷或边界条件的正确性Marc的后处理过程即为采集求解器处理分析的结果，提取用户所需的信息，了解计算机结果的过程后处理处理器可以用来查看整个模型在某一时间段的计算结果和用来查看模型的某一局部在整个时间段上的计算结果，还可以生成动画，Marc提供的现成的动画制作功能已经非常丰富，各种计算结果的变形动画、时间历程动画、切片动画、粒子轨迹动画等三、实例容用Marc/Mentat 分析包括大变形、自动热接触的热-机耦合问题的实例结构是一受钢块挤压的铝环，二者变形均予考虑，铝环变形过程中的变形生热和热接触也包括在四、求解步骤 1、网格定义MAINFILESOKMESH GENERATION (创建网格模型)SRFS: ADD 〔生成曲面〕POINT( 0, 13.5, 0)POINT(9, 13.5, 0)POINT(9, 27, 0)POINT(0, 27, 0)SRFS: ADDPOINT(9, 0 , 0)POINT(15, 0, 0)POINT(15, 42, 0)POINT(9, 42, 0)CONVERT (将生成的面转化为网格)DIVISIONS 9 15 〔指定在两个方向上单元的划分数目〕SURFACES TO ELEMENTS 1 〔将曲面1 转化为板状9*15 单元〕END LIST(#)DIVISIONS 6 25SURFACES TO ELEMENTS 2 〔将曲面2 转化为板状6*25 单元〕END LIST(#)RETURNSELECT 〔选择功能〕STORE:ELEMENTS 〔键盘输入集合名domain 1〕domain 1 〔从单元1 到单元135，如下图中的A 所指〕OKEND LIST(#)STORE:ELEMENTS 〔键盘输入集合名domain 2〕domain 2 〔从单元136 到单元285，如下图中的B 指向〕OKEND LIST(#)2、几何定义MAINGEOMETRIC PROPERTIESNEWAXISYMMETRICSOLIDCONSTANT DILATATION(on) (常体积)OKELEMENTS:ADDdomain 1(on)OK3、边界条件MAINBOUNDARY CONDITIONSNEWTHERMALFIXED TEMPERATUREON TEMERATURE(TOP) 20OKNODES:ADD (选择钢块右端边上的所有节点，指定钢块右端的温度为20℃)END LIST(#)NEWMACHANICALFIXED DISPLACEMENTON X DISPLACEMENTOKNODES:ADD (选择铝块左段的所有节点，指定其左端无位移)END LIST(#)NEWFIXED DISPLACEMENTON Y DISPLACEMENTOKNODES:ADD (钢块最下端的边上的所有节点，指定钢块环向无位移)END LIST(#)NEW 〔每个增量步钢块右端向左前进，用表格控制〕FIXED DISPLACEMENTOKTABLESNAME disp.TABLE TYPE timeOKADD POINT0 00.03 1FITRETURNFIXED DISPLACEMENTON X DISPLACEMENT -4.5(mm)TABLE disp.OKOKNODES:ADD(选择钢块右端边上的所有节点，指定其轴向位移)END LIST(#) 〔下示左图〕RETURNNEWTHERMALPLASTIC HEAT GENERATIONON(on)OKELEMENTS:ADD 〔下示右图〕SET domain 1OK4、定义初始条件铝环温度为427℃，钢块温度为20℃。
MAININITIAL CONDITIONSNEWTHERMALTEMPERATUREON TEMPERATURE(TOP)427OKNODES:ADD (指定铝块的初始温度为427℃，节点1-160)END LIST(#)NEWTEMPERATUREON TEMPERATURE(TOP)20OKNODES:ADD (指定钢块的初始温度为20℃，节点161-342)END LIST(#)5、材料性质铝环为热弹塑性E=1000N/mm2，，3，α=1.3*10-5，初始屈服应力2，相对于200℃的温度参考值200℃材料的硬化规律是线性的，且服从使屈服。