材料成形现在模拟技术

1、西 安 建 筑 科 技 大 学研 究 生 课 程 考 试 试 卷考试科目：材料成型现代模拟技术课程编码：任课教师：邹德宁考试时间：学 号：1204240572学生姓名：宋晔题 号成 绩总 成 绩学 分123456阅卷人签字789试题总页数101. 通过考察处于导热过程中各向同性物质的微元体积（xyz）的能量平衡，见下图1，建立导热微分方程式。（15分）xzydQx+xdQxdQzdQydQz+z图1 直角坐标系导热方程式的微元体解：首先要了解傅立叶定律，傅立叶定律其中，分别为x，y，z坐标轴上的单位矢量。为导热率（单位）。其含义表示，单位时间内，通过某单位截面上的热流（单位），与该处的温度梯度

2、成正比，但方向相反。 其次，导热微分方程的推导依据，根据能量守恒定律与傅立叶定律，建立导热物体中的温度场应满足的数学表达式，即导热微分方程；，物体在单位时间内获得的热量；，物体在单位时间内内能的增加；，物体对外界所做的功。对于固体来说，温度改变导致体积变化对环境所做的功可忽略不计，上式变为： 最后，直角坐标系下导热微分方程的推导 考察时间内微元体中：导入与导出净热量 + 内热源发热量 = 热力学能的增加1. 导入与导出微元体的净热量（1）时间内、沿x 轴方向、经垂直于x轴 的热量导入表面导入的热量： (单位)同理，时间内、沿x 轴方向、经垂直于x轴 的热量导出表面导出的热量： (单位)，分别为

3、热量导入面和导出面上的热流密度，单位。请注意，事实上这里有： ，所以导入与导出的热量差为： (单位)同理：（2）时间内、沿y 轴方向、经垂直于y轴 的两表面导入导出的热量差： (单位)（3）时间内、沿z 轴方向、经垂直于z轴 的两表面导入导出的热量差： (单位)2 微元体自身的发热量时间内，微元体自身的发热量：3微元体热力学能的增量（即微元体温度升高耗费的能量）时间内，微元体温度升高耗费的能量：根据前面所述的能量守恒，有：即整理得：又因为傅立叶定律，即，所以：， ， ，带入上式，得直角坐标系下的导热微分方程：2. 右下图所示了绝热、给定热流密度qr、对流、定温、辐射换热等五种边界表面的二维矩形

4、区域网格，请用能量平衡法建立各种换热边界条件下的节点方程。（15分）解：位于平直边界上的节点如图所示边界节点(m,n)只能代表半个元体，若边界上有向该元体传递的热流密度为，据能量守恒定律对该元体有： （4-4a）若时，则： （1）求解区域边界条件如下：给定温度边界条件，y=l2，0xl1，t=tw热流边界条件，y=0，0xl1，-kty=qr绝热边界条件，x=0，0<y<l2，tx=0对流传热边界条件，x=l1，0<y<l2，ktx=h(t-tf) 将上述条件分别代入（1）就可求得各换热边界条件下的节点方程。3 针对给定文献，提交学习报告。（15分）解：1.介绍及提出问

5、题阿维迪公司ISP薄板坯连铸连轧生产线钢水包容量为120吨，结晶器出口处连铸坯厚60毫米，经液芯压下后为43毫米，经三个机架压下至15毫米25毫米。由于采用了液固两相坯料轧机，轧制压力较小，每个机架功率为600千瓦。中间带坯经剪断机、感应加热器后送至克雷莫纳炉，除保温外，还可将带坯尾部变成头部，使其先进入精轧机以减少头、尾部温差。精轧机原设计为4机架，后增设1机架。精轧机后设有冷却轨道和地下卷取机，由连铸机至带钢卷取机的全长仅为168米。 工厂布局及操作。根据设计，克利蒙娜厂计划在1994年产能达到55万t/a，相当于每天生产16炉钢，最大连浇炉数为4炉钢水。在连铸坯芯部为液态的情况

6、下就开始减薄厚度，3个大压下量的轧机机架将连铸坯的厚度减薄到结晶器出口厚度的15%。经过感应加热以后，薄板坯在卷取炉内卷取并保温。经过除鳞后，薄板坯进入4辊精轧机。根据这种布局，钢液可在170m的距离和15min内转变为钢卷。ISP(内联带钢生产)和AST(Arvedi钢技术)过程, 允许在轧机入口具有非常高的温度。即离开铸造机的滚动的钢板,允许利用在内联机连铸存储的热。内联的轧制设备生产过程包括铸带轧制机,液芯还原设备,高还原机、感应加热器,克雷莫纳绕线站,除鳞机,传统精轧机和冷却区。最初的模具带厚度为74毫米, 在液体核心还原过程是减少到55毫米。这个地区的最大带温度高还原机放置在带中心，

7、主要好处技术有:1.非常低的投资成本；2.不需要隧道炉,英语和好的产品质量；3.供暖能耗下降到0时,连铸速度是在约0.14米/秒；4.高达20倍降低水的消耗；5.(相比传统轧制)温度梯度在截面是非常有用的对于轧制过程；6.低水平的安装机功率高还原机(3滚动站在0.5,0.6和0.8 MW)提供还原从55到12.5毫米的带钢宽度1300毫米。AST是最紧凑的热带钢生产过程用振动模技术与优秀的效率和成本。设备的总长度是70到80米包括铸造机在前线和最终冷却器在后部。最大温度带发生在横截面并且随着连铸速度从1340到1420摄氏度变化。表明中间区域经历卷取还是多空隙的。主要提到的新技术的目标是显著地

8、降低轧制力和达到非常有利的温度场相比传统流程。然而某些特定问题针对金属处理提升。中央部分材料仍是多空隙的。这个结果在改变材料密度和发生温度的特点,具有重大影响的塑料的行为材料。一个至关重要的问题也缺乏相关数据材料的热,机械性能和显著的变化的密度。数学和实验建模的糊状钢变形是一个创新主题对于非常高的温度范围内变形过程。对有色金属的变形测试更容易。不断上升的通过热机械应力的能力使钢样品的测试和模拟结果都计算机仿真和改进新的轧制技术,类似于ISP和AST流程。本章揭示了这些问题。它侧重于轴向对称的计算机模型,从而保证了其实验验证的可能性的帮助物理模拟。2半固态钢的热机械变形模型数值分析样品的变形有液

9、相的中央部分显示了极高的应变不均匀性要求应用程序混合分析数值解的问题，该模型致力于建模流程,这需要精度高产生的参数字段。分析解决方案的不可压缩性和质量都保护(糊状区)条件模型的重要部分。分析条件消除问题造成无意的标本量的变化应用数值方法。现有的、物理变化的钢的密度糊状区有影响的变化控制体积。真正另一方面数值错误可以是一个来源的容积损失造成干扰真正的改变。这个效果是非常不受欢迎的行为在造型的钢的热温度范围为转换特征的聚集状态。本文做了以下几个步骤a.热模拟b.控制收缩率的力学模型分析c.密度变化3相关开发的数值模拟系统Def_Semi_Solid和JMatPro（JMatPro是英国SenteS

10、oftware公司开发的，一款功能强大的材料性能模拟软件）计算机辅助测试的力学性能在非常高的钢温度使用Gleeble®3800物理仿真器来避免问题的出现传统的测试程序。第一个模块允许建立新项目或与之前已有的工作。整数部分的每个项目有:输入数据一个特定的压缩/拉伸试验以及结果的测量和优化。在当前版本的程序模块允许应用程序的数量的数据库引擎(在其他标准MSAcces,数据库IV和悖论7 - 8进行实时的系统),并允许实施材料数据库和程序的自动数据验证。第二模块(求解)给用户管理的可能性工作条件的仿真过程。逆分析可以被关闭或使用这个系统的一部分。4.计算机辅助实验过程5.结论 计算机辅助测

11、试的钢样品变形共存液体和固体阶段需要解决许多问题。其中一个问题是如何确定材料热性能和机械性能，如:热传递系数和其他热性能。最主要的问题是如何描述压力测试得到的拉力-压力曲线。该模型与不可压缩性条件解析形式允许模拟材料的变形与糊状区避免容积损失,引起一些密度。提出的程序是一个极好的半固态独特的工具,这可以是非常有用的,并可能使正确的解释结果非常高的吗温度测试。本文表明其预测能力有关:温度、形状的大小和变形区。关注的焦点是力学性能研究钢的具体特性的理论模型应用到分析。一个可以观察到压缩测试解释可能只是由于应用对模型和实现的反分析。4 结合具体实例，采用正交试验，进行实验设计。（15分）解：实例：为

12、了提高某产品的合格率，经分析考虑了温度A，时间（B），用碱量（C）三个因素，再考虑三个水平，采用L934正交表的前三列进行正交试验，合格率为31,54,38,53,49,42,57,62,64；1)用直观分析法和极差分析法分析试验结果，确定最适应的因素水平组合；2)对实验结果进行方差分析后，确定各因素的显著性次序及适宜水平组合；3)计算各因素的贡献率。（1）计算数据数据分析：1、直观法：第9方案 y=64 ，最佳方案为： 2、极差法：，各因素的趋势图：方差分析和波动贡献率计算表，最佳水平组合是 ：，考虑B为不显著因素，取经济方案。5 借助MATLAB软件，对焊接过程二维不稳态焊接热传导问题进行

13、数值模拟。（20分）解：薄板焊接过程温度场分析：取焊件的一半为模型进行离散化，起点为O点，以后以v速度沿y轴运动，如图所示。 根据题意该问题为二维非稳态导热，导热方程为：题目可化为以下微分方程组：(以y轴正方向为上，x轴正方向为右)相关的实验数据：参数 数值 备注 7.82 g/cm3密度v 0.4 cm/s 焊接速度h 1cm 板厚度Qm 4000 cal/cm3热源分布密度 0.0008cal/cm2·s·换热系数Te，T020 周边介质温度，初始温度 k 0.1 cal/cm·s· 导热系数 实验过程： 区域设置启动MATLAB软件，在命令窗口输入

14、pdetool，打开PDE Tool Box工具箱。单击Rectangle/Square矩形选框工具，在窗口中拉出一个矩形。双击矩形区域，在Object Dialog对话框输入Left为0，Bottom为0，Width为2，Height为2。与默认坐标相比，图形显示不在屏幕中间，所以还需要调整坐标显示比例。选择Options Axes Limits，勾选X、Y轴的自动选项。然后在Options Application选择Heat Transfer热传输模型。边界条件设置单击，使边界变红色，然后分别双击每段边界，打开Boundary Conditions对话框，如下图所示。根据边界条件，采用Ne

15、umann条件，设置各边界条件。边界条件输入值见下表。边界条件设置值边界 g (Heat Flux) q (Heat Transfer Coefficient) 左边界0 0 右边界0.0008*20 0.0008 下边界（x轴） -0.0008*20 -0.0008 上边界0.0008*20 0.0008 方程类型设置单击，打开PDE Specification对话框，设置方程类型为Parabolic（抛物型），rho(密度）为7.82，C（比热）为0.16，k（导热系数）为0.1，Q（热源）为4000*exp(-3*(x.2+(y-0.4*t).2)/0.49)，其他参数为0。网格划分单击

16、，或者加密网格，单击。初值和误差的设置 单击Solve菜单中Parameters选项，打开Solve Parameters对话框，输入Time，如Time=5s则输入0:0.5:5，u(t0)为20,其他不变。解方程和整理数据单击，开始解方程。单击Mesh->Export Mesh输出pet的数值，单击Solve->Export Solution输出u。单击，设置所绘制的图形的颜色、维度等，单击Plot即可完成。6 根据专业特点，选择任意模拟软件，对一工程问题进行模拟研究。（20分）解：1.问题描述分析模型如图1所示，梁的横截面为矩形 宽×高 = 1 × 2 m

17、2。受到分布载荷作用。材料的弹性模量200 GPa，泊松比0.3。图1 悬臂梁的计算分析模型2.计算模型2.1有限元模型及网格划分此问题属于平面问题，采用的ANSYS单元信息为：Element Type：Solid Quad 4node 182；K1：Reduced integration；K3：Plane Stress；网格划分信息为：长边NDIV：50；短边NDIV：10；Mesh：Areas；Shape：Quad，Mapped。即将用实体单元进行离散，网格划分为50 × 10的正方形网格。网格划分图见图2。图2 网格划分图2.2边界条件本问题的边界约束条件为给悬臂梁左边施加一固

18、定约束，并给梁的上边施加线性分布载荷。2.3材料模型本问题具体的材料模型和材料参数为：Structural：LinearElasticIsotropic；即材料属各向同性，按线弹性模型处理；材料参数：EX：200e9；PRXY：0.3，即弹性模量为200GPa，泊松比为0.3。2.4特殊处理本问题一般生成步骤是先生成特征点，而后依次连接四个特征点生成面。具体过程为：(1) 生成特征点Main Menu：Preprocessor Modeling Create Key points In Active CS input:1(0,0,0), 2(10,0,0), 3(10,2,0), 4(0,2,0) OK；(2) 生成面Main Menu：Preprocessor Modeling Create Areas Arbitrary Through KPS 依次连接四个特征点，1 2 3 4