【精品】计算机在材料科学与工程中的应用97.ppt

计算机在材料科学与工程中的应用,叶卫平,第四章材料学物理场数值模拟分析,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,本章要点,材料学主要物理场,传热学基本概念,有限差分求解,简单传热学计算机分析,材料学中的应力场简介,材料学中的浓度场简介,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.1材料学主要物理场,温度场,浓度场,应力场,temperaturefieldsimulateanalysis,材料热处理；焊接温度场计算；铸造凝固过程分析；陶瓷烧结和水泥处理；新的加热方式,材料学主要物理场简介,1,位错应力场，热处理相变应力场；焊接应力场计算；铸造凝固应力场计算,浓度扩散分析应力场；渗碳浓度场计算；,上机选题：刃型位错应力场分量的计算机模拟,上机选题：渗碳过程浓度场的计算机模拟,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.1材料学主要物理场,加热设备,冷却过程,加热方式,电阻炉、燃料炉、浴炉、流化床炉、真空炉等,感应加热、电子束加热、激光表面处理、离子轰击加热等各种冷却介质的冷却性能和各种冷却方式高效率、节能的新加热方法。,各种冷却介质冷却，各种冷却方式。,材料学中传热学重要课题,2,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.2传热学基本概念,稳态过程,稳态/非稳态热传递,1,非稳态过程,凡物体中各点温度不随时间变的热传递过程,凡物体中各点温度随时间变化而变化的热传递过程,instability-heat-ransfer,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.2传热学基本概念,fourier定律,导热基本定律,2,在导热现象中，单位时间内通过给定截面所传递的热量，正比例于垂直于该截方向上的温度变化率，而热量传递的方向与温度升高的方向相反。,qx热流密度；材料热传导系数（导热率）w/(mk)；负号表示传热方向与温度梯度方向相反,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.2传热学基本概念,三维非稳态导热微分方程,导热微分方程,3,为材料的密度(kg/m3)，c为材料的比热j/(kgk)=常数，t为时间(s)，热传导系数w/(mk)；q=q(x,y,z,t)是物体内部的热源密度(w/kg),物理意义：反映了物体的温度随时间和空间的变化关系,讨论：1)q为常数；2)稳态,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.2传热学基本概念,三维非稳态导热微分方程,导热微分方程,3,直角坐标系,柱坐标系,球坐标系,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.2传热学基本概念,定解条件,4,第一类,初始时间温度的分布条件,导热物体边界上温度或换热情况,规定了边界上温度值，即t|s=tw(x,y,z,t),第二类,边界上的热流密度q已知，即,讨论：绝热边界,第三类,指物体与其周围环境介质间的对流换热系数h和介质的温度t已知。（tw为物体边界上的温度，tc为介质温度）,讨论：非稳态,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.3有限差分求解,二维稳态导热,1,网格划分,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.3有限差分求解,二维稳态导热,1,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.3有限差分求解,二维稳态导热,1,中间节点用方程(3)计算,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.3有限差分求解,二维稳态导热,1,四周节点用方程(4)计算,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.3有限差分求解,二维稳态导热,1,由方程(3)与边界方程(4)组成定解问题的方程组,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.3有限差分求解,二维稳态导热,1,构成矩阵形式,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.3有限差分求解,有限差分法求解步骤,2,确定导热微分方程式、初始条件和边界条件;对区域进行离散化,确定计算节点s建立离散方程,对每一个节点写出表达式;求解线性方程组;对解进行分析,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,二维稳态问题求解,1,eg1.边界点温度已知，且区域内无内热源。利用有限差分方法来计算节点1、2、3的温度。,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,二维稳态问题求解,1,eg1.边界点温度已知，且区域内无内热源。利用有限差分方法来计算节点1、2、3的温度。,解方程得,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,二维不稳态焊接热传导问题求解,2,eg2.薄板焊接中移动热源,二维焊接离散化模型,(qm-最大热源密度；h-板厚；r-离开热源中心距离),计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,二维不稳态焊接热传导问题求解,2,eg2.薄板焊接中移动热源,二维焊接离散化模型,边界条件,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,二维不稳态焊接热传导问题求解,2,根据二维不稳定导热方程(1)以及初始条件、边界条件,可建立差分方程(2),根据差分方程以及各边界节点差分方程,在满足稳定条件下就可求出不同时刻的温度分布。,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,二维不稳态焊接热传导问题求解,2,时间步长ds=0.01s；热源作用时间s1=0sse=5s=0.0008cal/cm2s；qm=2500cal/cm；板厚h=1cm；焊速v=0.05cm/s,节点数n1=n2=10,t0=ta=20,比热c=0.16cal/g;密度=7.82g/cm2;导热系数k=0.1cal/cms,节点间距dx=0.2cm,输入参数,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,二维不稳态焊接热传导问题求解,2,计算框图,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,二维不稳态焊接热传导问题求解,2,0.5s温度场,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,二维不稳态焊接热传导问题求解,2,3.0s温度场,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,matlabpde二维不稳态焊接热传导求解,3,thepartialdifferentialequation(pde)toolboxcontainstoolsforthestudyandsolutionofpdesintwospacedimensions(2-d)andtime,usingthefiniteelementmethod(fem).itscommandlinefunctionsandgraphicaluserinterfacecanbeusedformathematicalmodelingofpdesinabroadrangeofengineeringandscienceapplications,includingstructuralmechanics,electromagnetics,heattransfer,anddiffusion.,inthematlabcommandwindow:demo(toolbox,partial),计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,matlabpde二维不稳态焊接热传导求解,3,区域设置,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,matlabpde二维不稳态焊接热传导求解,3,边界条件设置,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,matlabpde二维不稳态焊接热传导求解,3,方程类型设置,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,matlabpde二维不稳态焊接热传导求解,3,网格设置,初值和误差设置,解方程,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,matlabpde二维不稳态焊接热传导求解,3,1.5s时温度场,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.4简单传热学计算机分析,matlabpde二维不稳态焊接热传导求解,3,3s时温度场,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.5材料学中的应力场简介,材料学中的应力场,1,热应力分析,残余应力分析,薄膜应力分析,工程材料应力分析,涂层残余应力分析,相变应力分析,应力分析,位错应力场分析,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.5材料学中的应力场简介,轴承滚子组合凹模预应力分析与凹模设计,2,叶卫平等，轴承组合凹模失效分析及提高寿命的研究,武汉工学院学报,1994(1)40-45,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.5材料学中的应力场简介,轴承滚子组合凹模预应力分析与凹模设计,2,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.5材料学中的应力场简介,轴承滚子组合凹模预应力分析与凹模设计,2,r和模具径向和周向应力；d1-凹模内径(mm)；d-凹模外径(mm)；d2-凹模外套外径(mm)；d*-计算预应力处直径。q为过盈量()和材料弹性模量参数有关的函数,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.5材料学中的应力场简介,轴承滚子组合凹模预应力分析与凹模设计,2,凹模预应力计算框图,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.5材料学中的应力场简介,轴承滚子组合凹模预应力分析与凹模设计,2,凹模预应力分布计算,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.5材料学中的应力场简介,轴承滚子组合凹模预应力分析与凹模设计,2,注:三组合凹模过盈量为凹模与增强圈和增加圈与外套的过盈量。,结论,1)结构改进后的三组合凹模；2）用65nb钢代替原crl2mov钢使用寿命明显提高,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.6材料学中的浓度场简介,浓度场简介,1,固体材料扩散是物质传输的唯一方式，通常是指材料内部原子迁移的微观过程，以及由于大量原子迁移而引起物质的宏观流动。扩散与材料在生产和使用过程中的许多重要的物理化学过程密切相关，因此对扩散的浓度场的计算具有重要的意义。,fick第一定律;fick第二定律,浓度场数值模拟,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.6材料学中的浓度场简介,浓度场简介,1,fick第一定律;,浓度场数值模拟,fick第二定律,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.6材料学中的浓度场简介,用pde工具箱求解浓度分布,2,对一10cm10cm5cm的碳含量为wc0.8%长方体钢材进行渗硼，渗硼温度为850。已知850时硼在t8中的扩散系数d=1.1010-13m2/s，分析渗层厚度随时间（09000s）变化的情况。,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.6材料学中的浓度场简介,用pde工具箱求解浓度分布,2,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.6材料学中的浓度场简介,用pde求解扩散方程,3,tosolvethediffusionequationusingtheforward-time/centered-space(ftcs)scheme.,ref:alejandrogarcia,numericalmethodsforphysics,2e,2002.,run:dftcs.mentertimestep:0.0001enterthenumberofgridpoints:51solutionisexpectedtobestable,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.6材料学中的浓度场简介,用pde求解扩散方程,3,entertimestep:0.0001enterthenumberofgridpoints:51solutionisexpectedtobestable,计算机在材料科学与工程中的应用第四章材料学物理场数值模拟分析,4.6材料学中的浓度场简介,用pde求解扩散方程,3,tosolvetheneutrondiffusionequation(中子扩散方程)usingtheforward-time/centered-space(ftcs)scheme.,ref:alejandrogarcia,numericalmethodsforphysics,2e,2002.,run:neutrn.mentertimestep:0.0005enterthenumberofgridpoints:61entersystemlength:2=systemlengthissubcriticalsolutionisexpectedtobestabl