不同模式激光加热反应毕业设计

《不同模式激光加热反应毕业设计》由会员分享，可在线阅读，更多相关《不同模式激光加热反应毕业设计（54页珍藏版）》请在装配图网上搜索。

1、安徽工业大学毕业设计说明书（论文）安徽工业大学本科毕业设计(论文)任务书课题名称不同模式激光加热反应学 院专业班级姓 名学 号毕业设计(论文)的主要内容及要求：主要内容：对不同模式的激光对材料辐照过程中的温度效应及热力效应进行分析。本论文主要采用了有限元软件ANSYS对Al2O3陶瓷模型的建立并且模拟两种模式的激光（高斯形激光和帽顶形激光）对Al2O3陶瓷的加热过程及相关的温度响应分析以及对Al2O3陶瓷扫描的应力响应分析。 要求：1.具有查阅相关文献的能力； 2.熟练使用ANSYS软件对材料的建模，相关参数的设置及激光的加载； 3. 对激光辐照过程进行模拟，对数据进行分析整理，完成相关图表；

2、 指导教师签字： 不同模式激光加热反应摘 要激光与材料的相互作用一直是备受人们关注的问题，从工业上的激光热处理技术，到军事上的激光破坏机理研究，均与之密切相关，它在激光加工应用等领域中的诱人前景，激励着这方面的研究工作不断前进。本文利用有限元理论，对不同模式激光加热Al2O3陶瓷的热传导过程及热应力分布进行了数值模拟,得到了Al2O3陶瓷的温度场和热应力场分布。根据热传导方程和热应力方程，建立了不同模式的激光辐照下中瞬态温度场和热应力场的物理模型。模拟计算了Al2O3陶瓷在不同模式激光照射下的温度场。数值模拟结果表明：吸收的激光能量主要分布在光斑半径以内，因此随着加热时间增加形成的径向温度梯度

3、也越大；陶瓷内部沿轴向不同深度各点的温度随时间增长曲线在加热初始阶段均为S型。 本文的研究结果可为激光切割陶瓷加工过程的数值模拟研究提供参考。关键词 激光加工 Al2O3陶瓷 数值分析 ANSYS 温度场 热应力场DIFFERENT MODES OF LASER HEATING REACTIONABSTRACTLaser-material interaction has always been the concern of issues, from the industrial laser heat treatment technology to the military mechanism

4、of laser damage are closely related to it，its attractive future in the field of laser processing applications inspired the research work in this area continues to advance.In this paper, the heat conduction process of different modes laser heating Al2O3 ceramics is numerical simulated by using the fi

5、nite element method (FEM).the temperature field and thermal stress distribution of Al2O3 ceramics is obtained.The physical model of the transient temperature field and thermal stress field under different modes laser irradiation is established in accordance with classic thermal conduction equation a

6、nd thermal stress equation. This paper simulated the temperature field and the coupled stress field under laser of different modes. Numerical results indicate: the laser energy absorbed mainly distributed within the spot radius, for this reason, The radial temperature gradient formed became larger w

7、hen the heating time increased; the temperature growth curves over time of each point to a different depth along the axial direction in the ceramic internal are all S-type in the initial stage of heating. The results in this article may provide the research method and theory for the numerical simula

8、tion of laser cutting ceramic processing.Keywords: laser process ;Al2O3 ceramics; numerical analysis ;ANSYS; temperature field; thermal stress field;目录摘 要3ABSTRACT4第一章 绪 论81.1 课题研究背景81.2 激光与材料相互作用的研究现状101.3 本文研究工作简介11第二章 有限元理论 激光辐照材料的温度场和应力场理论122.1 有限元理论和ANSYS简介122.1.1 有限元的基本思想122.1.2 有限元法的分析过程132.1

9、.3 ANSYS软件简介132.2 激光辐照材料的温度场理论132.2.1 热导方程与定解条件的理论公式：142.2.2 激光辐照固体材料的理论模型152.3 激光辐照材料的应力场理论162.3.1 平衡微分方程162.3.2 几何方程172.3.3 物理方程172.3.4 热弹性方程20182.3.5 控制方程、初始条件、边界条件2119第三章 不同模式的激光辐照材料的温度场分析203.1 激光和材料的参数203.2 高斯激光的温度场分析203.2.1 热传导模型213.2.2 高斯激光的载人223.2.3 计算结果温度场分布图及分析233.3 帽顶形激光的温度场分析263.3.1 帽顶形激

10、光的载人263.3.2 计算结果、温度分布图及分析263.4 两种模式的激光产生的温度场的比较29第四章 激光辐照材料的应力场分析304.1 模型创建、载入激光304.2 计算应力结果及分析31第五章 结论与展望335.1 结论335.2 展望33致谢34参考文献34附录 1 英文原文36附录 2 译文46第一章 绪 论1.1 课题研究背景自从第一台红宝石激光器1960年诞生以来，激光技术的发展已经过了五十多年的历程，几乎渗透到自然科学研究的每一个领域。如日常生活中的激光打印机、CD唱盘与光盘、VCD与DVD视盘、光纤通信、激光测距、激光扫描条码。在工业上，激光被用于切割、焊接、标记与热处理等

11、加工领域。中国的工业激光起步并不比西方国家晚，但由于基础工业的薄弱和投资力度的不足，使商品化的激光切割机研究远远落后与发达国家。目前，人们已经研制出各种激光器，如固体激光器、液体激光器、气体激光器、化学激光器、准分子激光器和半导体激光器等。由于激光具有高单色性、高相干性、高方向性、高亮度和高稳定性等一系列突出的优点，使之容易获得很高的光通量密度。一般的固体或气体激光器，可以将激光束绝大部分能量聚焦在激光焦点上，这就是激光作为热加工热源的重要原因之一1。激光与材料的相互作用一直是备受人们关注的问题，从工业上的激光热处理技术，到军事上的激光破坏机理研究，均与之密切相关，它在激光加工应用等领域中的诱

12、人前景，激励着这方面的研究工作不断前进。激光与材料的相互作用，涉及激光物理、传热学、等离子体物理、非线性光学、固体与半导体物理、热力学、连续介质力学等广泛的学科领域， 激光与材料的相互作用既取决于激光特性和材料特性，也与作用的外部环境有关。激光特性包括波长、能量、功率、脉宽、脉冲结构、重复率等，材料特性包括光学材料和光学薄膜的反射率、透过率、吸收系数、热传导率、抗激光损伤强度等光学、热学、力学参数，其中任一种因素对相互作用的过程就是一项重要的研究课题2。当激光作用于材料表面时，入射激光的能量被分解为如下的几个部分:一部分被材料表面反射，一部分被透射，一部分被散射，剩下的一部分则被材料吸收。不同

13、的材料对不同脉宽和强度的激光的吸收机制是有所不同的，大致可以分为:逆韧致吸收、光致电离、多光子吸收、杂质吸收、空穴吸收等，一般而言逆韧致吸收和光致电离这两种机制起主导作用。材料吸收激光能量后，其中的粒子(电子、离子和原子)将获得过剩的能量。这些获得多余能量的粒子由于相互碰撞传递能量，材料的宏观温度将明显升高。当温度升高到熔点后，材料将发生熔融和汽化等现象。当作用激光强度足够高时，还会产生材料的烧蚀、等离子屏蔽和材料的冲击破坏等。因此，可利用激光束对工件进行刻标、切割、钻孔、焊接、热处理、重熔、表面合金等，应用前景远大3-5。热处理在材料技术中有着广泛的应用，诸如金属的软化或硬化，半导体参杂扩散

14、，陶瓷等脆性材料的切割，混合物的复合形成等，所以研究激光与物质的相互作用的热效应意义重大。激光加工指的是激光束作用于物体表面而引起的物体成形或改性的加工过程。以激光作为加工能源，在硬脆性陶瓷加工方面的发展潜力已见端倪6：它可以实现无接触式加工，减少了因接触应力对陶瓷带来的损伤；陶瓷对激光具有较高的吸收率(氧化物陶瓷对10.6m波长激光的最高吸收率可达80以上)，聚焦的高能激光束作用于陶瓷局部区域的能量可超过108 J/cm2，瞬间就可使材料熔化蒸发，实现高效率加工；由于聚焦光斑小，产生的热影响区小，可以达到精密加工的要求。由于陶瓷是由共价键、离子键或两者混合化学键结合的物质，晶体间化学键方向性

15、强，因而具有高硬度和高脆性的本征特性，相对于金属材料，即使是高精密陶瓷，其显微结构均匀度亦较差，严重降低了材料的抗热震性，常温下对剪切应力的变形阻力很大，极易形成裂纹、崩豁甚至于材料碎裂。因此，高效无损伤激光切割陶瓷类高硬脆无机非金属材料一直是一个的且亟待解决的问题7,8。研究激光与物质热作用过程的理论方法主要有解析法和数值计算法两种。由于涉及到材料熔融和汽化等相变过程的问题非常复杂，理论解析本解比较困难，甚至于根本无法求解。因此一般采用数值方法来处理相关问题。随着计算技术的飞速发展，人们逐渐提出了一些求解激光熔融问题控制方程的新方法，如将激光作用视为点热源、作用对象视为半无限大介质、采用一维

16、近似及假定一个简单的速度场而形成的解析或者半解析法等。而其中的数值计算方法可已采用有限差分法(FDM)和有限单元法(FEM)等，有限单元法包含边界元法、有限分析法。对激光与物质相互作用的数值计算就是从特定的物理模型出发，用计算机进行数值计算或模拟，从而揭示激光与物质相互作用的某些性质和运动规律9。激光与物质的热作用研究主要是计算激光辐照下物质温度的变化，以及温度变化引起的一些热效应。通常情况下，由于解析或半解析方法必须作一些与实际情况相距甚远的假设，计算结果可能与实际结果有较大偏差，但可作为进一步研究和纯数值计算的基础。1.2 激光与材料相互作用的研究现状1971年，Ready出版了研究高功率

17、激光产生系列效应的第一本专著。该专著问世之后，无论在强激光与物质相互作用的基础理论研究方面，还是激光加工以及激光的军事应用方面都有专著问世。在有关专著中，常根据经典电磁理论将物质折射率和消光系数作为常数看待，并与物质的电学特性相联系。前苏联学者在1968年引入精细的微观模型研究了物质光学的微观特性，即研究了n,与物质结构特性的关系，指出吸收系数和光波圆频率的关系决定了物质的吸收光谱，和的关系同样表征了吸收光谱的结构。Ready和Eloy等人10研究了物质对激光的吸收和转化效应，结论是固体材料中在趋肤深度内被吸收的激光能量直接转化为自由电子或束缚电子平均动能的增加，其中大部分的能量再通过电子与晶

18、格或离子的相互作用转化为材料表面层的热能，同时也存在一些其它转化机制。 Carslaws首先给出了简单的激光加热下物体的温度场的解析计解或近似解11. J. A. Mckay和J. T. Schriempf讨论了作用于物体表面的激光功率密度不均情况下靶的表面温度分布并给出了计算结果。而与各种实际问题有关的激光加热温度场的理论和数值研究，其中包括积分变换、积分近似法、有限差分和有限元计算等的研究报道已有不少结果。其中Ready和Warren主要考虑了一维问题，包括材料的激光吸收系数与温度有关及发生气化或熔化的情形; Bragger和 Bechtel的工作给出了高斯分布脉冲激光束照射圆板和方板产生

19、的温度场。Mezines考虑了重复频率激光对薄板的加热问题。当然，应用积分变换方法，原则上可得到轴对称分布、任意时间波形激光加热板材的温度场，包括表面换热和温度相关热物理性质的一般情况，以及镀膜板和多层组合板内(板间理想或非理想热接触)的温度场12,13。E1-Niclany等14对均匀、连续激光辐照下半无限体材料的加热、熔化和气化问题进行了解析研究。1986年，Feisch U等建立了一个LGA(Lattice Gas Automata)数值流体力学模型，经过近些年的不断发展和改进形成了今天的格子Boltzmann模型LBE(LatticeBoltzmann Equation Models)

20、。LBE模型在模拟各种复杂系统物理现象等方面得到了广泛的应用,但也存在一些不足。 我国在激光与物质的相互作用机理方面研究起步相对较晚，由于激光器件的限制，上世纪60, 70年代主要限于短脉冲激光束对靶材冲击破坏效应的理论研究，自80年代以后在强激光引起材料的热学和力学效应方面在基础理论研究、大型数值程序计算、实验测试手段、大型的演示实验等领域都取得了长足的进展。上世纪九十年代以来，国内在高功率激光与物质相互作用方面研究发展较快。闰长春等对脉冲激光加热下材料的温度场进行了研究。夏日源用自相似理论计算对脉冲激光加热下材料的温度场进行了研究。魏在福15等人采用数值模拟法分析了激光加热温度场分布，分析

21、了材料表面热辐射，激光作用时间及功率密度对温度场分布的影响。沈中华等16考虑了材料的热物性参数在计算过程中的变化和入射激光的空间分布，采用二维模型，得到Si材料在强激光作用下的轴向和径向温升分布，给出了表面层开始熔化的时间与激光功率密度关系。王伟平等人考虑了激光光强和时空分布对材料加热的影响，得出光强时间分布不同时，靶面温升的历史不同，但最终温度和温度分布都趋于一致。1.3 本文研究工作简介1、描述了Al2O3陶瓷的加热过程，根据经典热传导方程和热应力方程，建立了不同模式激光辐照下中瞬态温度场和热应力场的物理模型，讨论了激光辐射过程中陶瓷温度场及产生的热应力场分布。2、介绍了激光辐照Al2O3

22、陶瓷的数学模型，利用轴对称的热传导方程和伽辽金法得到了轴对称的有限元的热传导方程，并介绍了空间轴对称物体的温度分布基本理论和有限元解法。3、计算得到了Al2O3陶瓷在不同模式的激光加热时的温度场以及在施加约束后相应的热应力场，并进行了比较分析。第二章 有限元理论 激光辐照材料的温度场和应力场理论2.1 有限元理论和ANSYS简介有限单元法是以变分原理为基础，并吸取了有限差分法中离散的思想而发展起来的一种有效的数值解法。有限单元法则可方便地处理任何复杂形状和复杂边界条件的问题，亦可达到较高精度。2.1.1 有限元的基本思想有限元法是在连续体上直接进行近似计算的一种数值分析方法。这种方法首先是将连

23、续的求解区域离散为一组有限多个、且按一定方式连接在一起的单元的组合体，并且认为单元之间只通过有限个点连接起来，这些连接点称为结点(Node)。有限元法利用在每一个单元内假定的近似函数来分片地表示全求解域上待求的未知的场函数(如位移场、应力场和温度场等)。单元内的近似函数常由未知的场函数(或包括其导数)在单元内各个结点上的数值通过函数插值来表示。这样，未知的场函数(或包括其导数)在单元内各个结点的数值就成为新的未知量(也即自由度)，从而使一个连续的无限自由度问题变成为离散的有限自由度问题。一经求解出这些未知量，就可以通过插值函数计算出各个单元内场函数的近似值，从而得到整个求解域上场函数的近似值。

24、显然，随着单元数量的增加，也即单元尺寸的缩小，或随着单元自由度的增加和插值函数精度的提高，解的近似程度将不断得到改进。若单元满足收敛要求，近似解最后将收敛于精确解。由于单元能按不同的连接方式进行组合，并且其本身可以有不同的形状，所以几何形状复杂的问题都可方便地离散化，因此，有限元法可以方便地处理各种复杂因素，如复杂的几何形状、不均匀的材料特性、任意的边界条件、结构中包含不同类型构件等等，它们都能用有限元法灵活地求解。2.1.2 有限元法的分析过程有限元法的分析过程大体分为前处理、分析和后处理三大步骤。 前处理过程是对实际连续体进行离散化建立有限元分析模型的过程。在这一阶段，要构造计算对象的几何

25、模型，划分有限元网络，生成有限元分析的输入数据。这一步是有限元分析的关键。 有限元分析过程主要包括:单元分析、整体分析、载荷移置、引入约束、求解约束方程等过程。有限元分析的后处理主要包括对计算结果的加工处理、编辑组织和图形表示三个方面。它可以把有限元分析得到的数据，进一步转换为设计人员直接需要的信息，如应力分布状况、温度场分布、结构变形状态等，并且绘成直观的图形，从而帮助设计人员迅速地评价和校核设计方案。2.1.3 ANSYS软件简介数值模拟技术通过计算机程序在工程中得到广泛的应用。国际上较大型的面向工程的有限元通用程序达到几百种，其中著名的有：ANSYS，NASTRAN，ASKA，ADINA

26、，SAP等. 以ANSYS为代表的有限元分析软件将有限元分析、计算机图形学和优化技术相结合，已成为解决现代工程学问题的有力工具。本文就将采用ANSYS软件对激光加热Al2O3陶瓷过程进行仿真。ANSYS分析过程中包含三个主要的步骤首先，创建有限元模型(前处理)。(1)创建或读入几何模型；(2)定义材料属性；(3)划分网格。其次，施加载荷并求解。(1)施加载荷及载荷选项、设定约束条件；(2)求解。最后，查看结果。(1)查看分析结果；(2)检验结果(分析结果是否正确)。2.2 激光辐照材料的温度场理论激光照射靶材时，其能量被材料的表层所吸收并转变为热。该热量通过热传导在靶材内扩散，从而形成温度场，

27、该温度导致靶材性质的变化。本节将主要对激光加热过程中的热传导现象进行分析和计算，这对激光与材料相互作用的物理机理分析具有重要的意义。激光加载条件、物体的形状、初始和边界条件以及物体的热物理性质决定了物体温度场的时空变化。首先给定激光加载条件，通常把被吸收的激光能量或光强作为面热源(表层吸收)或体热源(深层吸收)处理。有时边界上还存在其他种类的热源或冷却条件，如表面辐射损失、气动加热或冷却等。不同介质之间界面的热接触条件比较复杂，理想情况下界面两边的温度和热流量应当相等17,18。由于在物体某处温度达到熔点时，温度暂时不再上升，直到热能达到能够提供继续融化的能量，才能继续融化。这个过程是十分复杂

28、的，在本文中不予讨论。本文只讨论在熔融之前为发生相变时的升温过程。在给定热源、边界和界面条件、初始温度分布和物体的热物理特性时，激光加热各向同性物体的问题就可以归结为不定常热传导方程的计算19。2.2.1 热导方程与定解条件的理论公式：三维热传导方程可写为如下通用形式： (2.1)式中K为热导率，为材料密度，c为材料比热容，T为温度，t为时间变量，A(x, y, z, t)为每单位时间、单位体积传递热给固体材料的加热速率。由于材料的热物理系数是温度的函数，所以方程(2.1)是非线性的，其解非常复杂，很难得到解析解。然而事实上大部分材料的热物理参数随温度变化并不明显，故在一定条件下可假定其与温度

29、无关，在一定的温度范围内取其平均值进行计算，这样方程(2.1)才可能得到解析解。若激光作用下材料是均匀各向同性的，则方程(2.1)可简化为: (2.2)式中为材料的热扩散率。由于求解热传导方程非常复杂，许多学者提出了一些热模型，在求解热传导方程时通常的假定条件是2,19.(1)被加热材料是各向同性物质。(2)材料的热物理参数与温度无关或取特定的平均值。(3)忽略热传导中的辐射和对流，只考虑材料表面的热传导。2.2.2 激光辐照固体材料的理论模型假定激光束垂直入射于物体表面(即x=0)，被加热物体位于右半空间()，物体表面对激光的反射率为R，吸收系数为，x=0处入射激光束的功率密度(光强)是，则

30、物体内部的温度场T可由热传导方程进行描述: (2.3)式中:热物理常数与位置及温度有关； Q是其他体热源项。上式右部是激光束深层吸收的体热源，如果或，用表面吸收率A代替(1-R)，此体热源项可改用边界条件中的面热源表示: (2.4)光斑处激光束的强度通常可表示为空间分布和无量纲时间波形B(t)的乘积。典型的B(t)波形有阶梯波、矩形波、瞬时作用(函数)、三角波和梯形波、指数衰减波形、Gauss波形和重复频率波形等。其他面热源以及常用的温度、热流和换热边界条件的提法与普通传热学问题相同。设t=0为激光开始作用的时刻，是给定的初始温度分布。如果初始时物体与环境等温，则热传导方程式计算的是物体的温升

31、，我们考虑常温即的情况。考虑被加热物体厚度大于激光脉冲作用时间内热量的传播深度的情况，物体几何形状为圆柱形，表面光斑上的激光强度分布对于光束轴线是旋转对称的，物体的温度场对于x轴也一定是旋转对称的。激光强度呈高斯分布，是光束的高斯半径，表面吸收，表面反射率和吸收率分别为R和A，时间波形为连续函数B(t),得到温度场分布为： （2.5）其中是光斑中心处光强的最大值，特别是阶梯波形(连续波)时，B=1，光斑中心温度是。2.3 激光辐照材料的应力场理论轴对称物体是指其几何形状中心轴线，则物体内各处的应力外加载荷或温度作用，以及约束情况都对称于某一应变及位移分量也对称于这一轴线。假设物体是连续的线性弹

32、性且是各向均匀同性的，在外力或温度作用下物体的位移和变形是微小的，物体在外加荷载或温度作用之前处于自然状态，内部没有应力存在。2.3.1 平衡微分方程 在物体内点P外取一微元六面体PABC，六面体由两个半径差为dr的圆柱面、两个夹角为且包含z轴的垂直面和两个相距为dz的水平面所围成，作用在圆柱面的径向正应力为。轴向正应力为。圆柱面上的剪应力，环向正应力。 令体力的径向分量为R，体力的轴向分量为Z。将微元体六个面的作用力都投影到半径方向，则 (2.6)经计算得空间轴对称物体的平衡微分方程为： (2.7)2.3.2 几何方程在轴对称的物体内每一点只能有两个方向的位移，即沿z轴方向的轴向位移及沿半径

33、:方向的径向位移，由于对称关系，物体内的任一点没有沿圆周方向(环向)的位移。径向应变为 (2.8)同理可得轴向应变为 (2.9)环向应变为 (2.10)剪应变为 (2.11)其几何方程为 (2.12)式中:为应变列向量 为径向位移 为轴向位移2.3.3 物理方程由广义虎克定律可得 (2.13)将应力分量表示为应变分量的函数为 (2.14)式中:为泊松比 E为弹性模量2.3.4 热弹性方程20 对于三维轴对称模型，在忽略体积力和惯性力的前提下，材料体内的温度场发生非均匀变化时产生的热应变和热应力可用如下的热弹性方程来描述: (2.15) (2.16)其中，,分别是 (r, z)点的位移在r, z

34、方向上的分量;e为体应变;为热膨胀系数;E,分别为杨氏模量和泊松比。而应变、应力之间分别满足如下的弹性力学关系式: (2.17) (2.18)对应材料的简支状态，热弹性方程的自由边界条件可取为: (a为半径，d为高度) （2.19)对应于靶材夹紧状态，热弹性方程的固定边界条件可取为: （2.20）2.3.5 控制方程、初始条件、边界条件21 1.控制方程 假设陶瓷板为各向同性，器热传导方程为22 (2.21) 式中：c,和k分别是陶瓷的质量热容、密度和导热系数；为体热源的生热率。2.初始条件23t=0, (2.22)式中：T0为陶瓷的初始温度，这里视作环境温度300K。3.边界条件在激光加载面

35、。热传导向陶瓷内部传送的热能，对流换热向外界逸出的热能以及表面向外界辐射的热能，三者之和应该等于其吸收的激光能量，即 (2.23)式中：为激光加载面的表面温度；h为对流传热系数；为表面热辐射系数；为波尔兹曼常数；为陶瓷对激光的吸收率；表示光能密度。在陶瓷板的前后左右和下表面，仅通过对流的方式与外界进行热量交换，因此设置为对流的边界条件24即 (2.24)式中：表示温度沿表面外法线方向的偏导数；表示玻璃表面的温度。激光辐射可以看作表面加热源25,因此方程（2.20）中的=0。方程（2.22）中的为面热源的光能密度，也是激光的功率密度，当激光沿着陶瓷基板的正x轴方向运动时，激光在陶瓷板表面的功率密

36、度可以用函数表示为 (2.25)则 (2.26)式中：为激光功率；r为激光束半径。第三章 不同模式的激光辐照材料的温度场分析3.1 激光和材料的参数激光功率: 激光束半径: 表3.1 Al2O3陶瓷的物理参数密度比热热传导系数热膨胀系数杨氏模量泊松比 396080033.54000.22陶瓷参数：熔点：2050K轴对称模型参数： 半径 5mm 厚度 2.5mm长方体模型参数： 长 20mm 宽5mm 高 1mm3.2 高斯激光的温度场分析本文建立了5mm*2.5mm的矩阵模型，为减少计算量，不影响模拟的情况下，划分网格时，为了把更多的计算用于温度变化明显的区域，采用了不等划分方法，在径向分为1

37、00份，深度方向分为50份。如图3.1图3.1 柱坐标下轴截面的网格划分在计算完结果后可以用ANSYS软件的3D还原功能将其还原。如图3.2图3.2柱坐标3D扩展后网格划分效果图3.2.1 热传导模型为了便于分析激光与陶瓷相互作用时的温度场，本文做以下假设：(l) 激光束直接作用与材料表面(在材料内没有新的热源)，这是因为激光的照射时间是短暂的。(2) 视陶瓷样品表面对激光能量的吸收率为1，且不随时间变化。事实上，吸收系数是根据温度变化而变化的，但是对相同表面的的条件下，测量一个平均值是可行的。(3) 激光照射期间功率等参数不变。(4) 厚板的维度的简化，选择适当的维数进行简化计算，长度，厚度

38、，宽度当到达一定尺度时，可以近似认为无穷大。(5) 忽略热辐射与热对流，认为空气在照射期间是绝热。(6) 忽略材料的密度随温度升高的变化。(7) 材料是各向同性的。(8) 激光束是稳态的。 在有了以上几个假定后，就可以用一般热传到方程求解温度场分布情况了。令激光垂直照射在玻璃板表面，图3.3 激光照射示意图3.2.2 高斯激光的载人初始条件 T0=300K 即视陶瓷的初始温度，也是环境温度。高斯激光光能密度：在柱坐标下分析时高斯激光光能密度公式为：其中是激光功率； 是激光束半径；x为径向深度。3.2.3 计算结果温度场分布图及分析图3.4 功率200W 光斑半径0.5mm 时间步长5ms高斯激

39、光照射陶瓷界面温度分布图由上图可以看出在距离加热点一定范围内经过激光辐照过的陶瓷,陶瓷表面(径向)温度传播速度要比陶瓷内部(轴向)要快,超过该范围,温度趋于一致(即环境温度)。中心点最高温度达到1867K。图3.5 径向距离中心原点分别为0、0.5、1.5、2.0mm随时间温度变化图图3.5说明了距离加热中心越近的温度越高且温度上升的速度要比远心点要快。图3.6 径向温度梯度在径向上，距离加热中心点越近的区域，温度梯度越陡，且梯度随距离的增大而减小，也是由于距离中心点的温度上升速度要比远心点要快，如图3.6。在轴向上的情况基本差不多，但是相比径向，轴向的梯度要比径向要小，且梯度变化相对径向要缓

40、和，如图3.7。图3.7 轴向温度梯度在轴向上，温度梯度要比径向要顺畅，在轴向上各点的温度随时间的增加温度趋势如半S形，如图3.8所示图3.8 轴向不同深度随时间变化温度图图3.9 不同方向上的温度变化在一定的范围内径向温度要比轴向温度高，超出该范围，温度基本一致，即环境温度300K，如图3.9所示。3.3 帽顶形激光的温度场分析3.3.1 帽顶形激光的载人图3.10 帽顶形激光示意图激光功率P=200W 光斑半径r=2.5mm3.3.2 计算结果、温度分布图及分析图3.11 功率200w光斑半径0.5mm时间5ms帽顶形激光照射陶瓷界面温度分布图相比于图3.4，帽顶形激光在轴向温度上要比高斯

41、形激光传播要快，中心温度达到2105K（假设此时陶瓷没有熔化）。说明中心点的温度要比高斯形激光产生的最高温度要高得多。图3.12 径向r=0.5、1.0、1.5、2.0、2.5mm温度随时间变化图一定范围内距离加热中心越近的温度越高且温度上升的速度要比远心点要快。超出该范围，温度基本趋于室温恒定不变，故温度基本无变化。图3.13 径向温度梯度在径向上，距离加热中心点越近的区域，温度梯度越陡，且梯度随距离的增大而减小，也是由于距离中心点的温度上升速度要比远心点要快，如图3.12。在轴向上的情况基本差不多，但是相比径向，轴向的梯度要比径向要小，且梯度变化相对径向要缓和，如图3.14。图3.14 轴

42、向温度梯度图3.15轴向不同深度随时间变化温度变化情况由图3.15 和图3.8可以看出在一定深度内，各部位的温度变化近似为半S形，且在深度超过0.6mm后，陶瓷温度基本和环境温度一致（300K）。图3.16 不同方向上的温度趋势3.4 两种模式的激光产生的温度场的比较本文中比较在相同的功率、时间步长、光束半径相同的不同模式的激光对同一种模型加热所产生的温度：1. 径向上的温度比较，如图3.17图3.17 不同模式激光径向温度场比较结论1：在一定的距离内，帽顶形激光的温度传播比高斯形激光传播要快，超过该距离，帽顶形激光要比高斯形激光慢，直到温度不变。不论是何种模式的激光，吸收激光的能量主要分布在

43、光斑半径以内。因此随着时间的增长形成的径向温度梯度也越大。2. 轴向上的温度比较，如图3.18图3.18 不同模式激光在轴向上的温度比较结论2：在轴向上，一定范围内，帽顶形激光要比高斯激光温度传播要快。结论3：在加热中心点（即光斑半径范围以内），帽顶形激光产生的温度要比高斯形激光产生的温度要高。第四章 激光辐照材料的应力场分析4.1 模型创建、载入激光本次模拟建立了20mm*5mm*1mm的长方体模型，本次工作用的是一个三维的耦合场固体元素SOLID5（该元素有八个节点，每个节点有六个自由度）。并对该模型在长度方向上分了100份，其模型及网格图如图4.1图4.1 应力模拟模型这次设计模拟的是移

44、动的激光束对Al2O3陶瓷基板中心轴线的扫描过程其示意图如图4.2所示图4.2 激光扫描陶瓷基本示意图相关参数：扫描速度v=2m/s;时间步长T=10ms;激光功率=200W；激光束半径r=0.5mm。载入激光方程（即式2.26）4.2 计算应力结果及分析热应力总是由于快速的升温或降温所产生的，当热应力超过材料的临界值的时候，材料就会产生断裂。正如本次工作中，材料出现了断裂。如图4.3所示图4.3 热应力是Y方向产生的位移在Y方向上产生的应力随X变化而变化，如图4.4所示，在一定距离内，陶瓷内部产生的是正应力（即应力方向背向陶瓷基板），超过该距离，内部便产生负应力（应力方向指向陶瓷基板）。过程

45、比较复杂，本次工作不做讨论。图4.4 Y方向上的热应力图4.5 X方向的温度温度变化最快的部分如图4.5，也就是产生断裂的区域，该部位剧烈的热对流产生了超过了陶瓷的临界值，故而产生断裂。图4.6 等效应力显示图4.7 陶瓷基板表面的温度分布在陶瓷基板上，产生了断裂，最高温度为1026K，要远低于熔点2050K,之所以断裂，由于陶瓷基板和环境之间的热对流，在陶瓷基板的表面，之前的加热热量也会沿着陶瓷基板的正方方向传播，而导致了最高温度出现在中部，这样会导致陶瓷基本中部很容易出现剧烈的热对流现象，使陶瓷基板的部分区域发生了迅速的升温和降温，从而产生一个极大的应力足够断裂陶瓷基板。如图4.6和图4.

46、7所示。第五章 结论与展望5.1 结论在不同模式的激光加热陶瓷的过程中，在相同的激光功率、加热时间步长、相同的环境、模型以及材料的情况下：1. 在陶瓷基板距加热中心点的一定距离内，在基本表面上，帽顶形激光的温度传播比高斯形激光传播要快，超过该距离，帽顶形激光要比高斯形激光慢，直到温度都接近环境温度。2. 在陶瓷的轴向方向上，一定范围内，帽顶形激光要比高斯激光温度传播要快，超出该范围，温度都趋于环境温度。3. 在加热中心点，帽顶形激光产生的温度要比高斯形激光瞬间产生的温度要高。4. 在激光扫描陶瓷基本的模拟实验中，激光加热可以使陶瓷基板没有达到熔点而出现断裂，断裂一般出现在陶瓷基板的中后部，这是

47、由于热对流的不均匀导致了局部出现了剧烈的升温或降温而产生超过临界值的热应力导致断裂出现。5.2 展望研究高斯形或者帽顶形激光与Al2O3陶瓷的相互作用过程对激光在工业、军事等各个领域的应用都有非常重要的意义，且作用过程尚有许多问题和工作需要探索、分析和解决。在激光切割陶瓷方面，由于陶瓷本征的硬脆特性 ，使得该技术的发展也存在着许多亟待解决的问题，对陶瓷激光加工的作用机理更需要得到理论上的论证和支持，下面列出一些有待解决的问题和研究内容。 1、考虑陶瓷表面热对流、热辐射，陶瓷材料的各向异性，热物参数随温度的变化及材料对激光的吸收系数随温度的变化对陶瓷温度场和热应力场的影响。2、陶瓷对激光能量的吸

48、收过程中是个复杂的过程，因此不同模式的激光对陶瓷作用时要综合考虑界面吸收、杂质吸收和体吸收。3、需要建立更加符合实际的三维模型，根据陶瓷实际情况施加约束条件得到位移分布和应力场分布，并且可以研究应力场随时间的变化，得到激光作用于陶瓷的应力场瞬态变化数据。4、研究帽顶形激光加载时可能造成的熔融，气化反应和应力引起的损伤和断裂情况。5、可以近一步研究移动激光束对陶瓷作用的温度场和热应力场变化，模拟激光焊接陶瓷过程。6、吸收激光的能量主要分布在光斑半径以内，因此随着加热时间增加形成的径向温度梯度也越大致谢感谢刘剑博士在我完成这篇论文的整个过程的帮助，回顾论文形成的论文的每一步，无不包含着刘剑博士的辛

49、勤指导，正是他深入浅出、不厌其烦的讲解，才使我克服一个个困难。借此机会，谨向刘剑博士致以最衷心的谢意和敬意。参考文献1 兰信钜等编著.激光技术M.北京:科学出版社，2000.2 郑启光，辜建辉.激光与物质相互作用M.武汉:华中理工大学出版社，1996.3 孙承纬.激光辐照效应M.北京:国防工业出版社，2002. 4 奥尔曼.激光束与材料相互作用的物理原理及应用M.北京:科学出版社，1994.5 陆建，倪晓武，贺安之.激光与材料相互作用物理学M.北京:机械工业出版社，1996.6 陈洪新，贾天卿，黄敏等飞秒激光的波长对SiC材料烧蚀的影响J光学学报，2006，26(3)：467-470.7 季凌

50、飞,闫胤洲,鲍勇,蒋毅坚.陶瓷激光切割技术的研究现状与思考J. 中国激光,2008,11(35): 1686-1693.8 张姗，康少英激光加工陶瓷强度的实验研究J应用激光，1996，16(5)：218-2209 王助成，邵敏主编.有限单元法基本原理和数值方法M.北京:清华大学出版社，1995.10 Swanson D G. Plasma Waves M.Boston: Academic Press, 1989.11 Glant V E,Zhilinsky A P, Sakharov I E. Fundamentals of Plasma Physics M.New York：wiley,19

51、80.12 彭惠民等.x射线激光M.北京:国防工业出版社，1997.13 Kruer W L. The Physics of Laser Plasma IneractionsM. Addison-Wesley Publishing Company, 1998. 66-90.14 Estabrook K G, Kruer W L, and Lsinski B F. Phys J. Rev. Letters, 1980(45):1399.15 魏在福，查洪连. 激光加热温度场物理分析J. 光学学报.1994，14(4):355.16 沈中华，陆建，倪晓武.硅材料在强激光作用下的二维温升模拟计算J.

52、南京理工大学学报,1998, 22(3):272-275.17 段晓峰,汪岳峰,牛燕雄,张雏. 激光辐照光学材料热力效应的解析计算和损伤评估J.中国激光，2004, 31(12):1455-1459.18 王玺,李化,聂劲松. 脉冲CO2 激光辐照K9 玻璃的热力效应分析J. 光学与光电技术，2008, 6(4):49-52.19 Carslaw H S, Jaeger J C. Conduction of Heat in solids M.London: Oxford Univ. Press, 1959.20 强希文，张建泉，刘峰等.强激光辐照半导体材料的温升及热应力损伤的理论研究J.中国激

53、光，2000, 27 (8) : 709-713.21 焦俊科，王新兵。二氧化碳激光作用下运动石英玻璃温度分布。华中科技大学 光电子科学与工程学院 2007,1001-4322.22 郑启光. 激光先进制造技术M.武汉:华中科技大学出版社,2002:34-36.23 袁红,赵剑衡,谭福利,等.激光辐照下旋转柱壳温度场的数值模拟J.强激光与粒子束,2005,17(5):681-684.24 席明哲,虞钢,连续移动三维瞬态激光熔池温度场数值模拟J,中国激光,2004,31,1527-1532.25 金花漫-在CO2激光辐射下玻璃温度变化的测量和计算问题J国外建材译丛,1997,(1):43-48.

54、附录 1 英文原文Numerical investigation on machining glass with CO2 lasersJunke JIAO1,2, Xinbing WANG ()11 Wuhan National Laboratory for Optoelectronics, Huazhong University of Science and Technology, Wuhan 430074, China2 Institute of Industry Technology, Guangzhou and Chinese Academy of Sciences, Guangzho

55、u 511458, ChinaAbstract When a glass substrate was irradiated by three different temporal shapes of laser sources, namely, line-time-shape laser, triangle-time-shape laser, and parabola-time-shape laser, the mathematical models were proposed,and the temperature distribution and the resulting thermal

56、 stress were calculated by the finite-element-method (FEM) software ANSYS. With these three types of lasers having the same output laser energy, the resulting thermal stress induced in the glass substrate was analyzed. The results showed that, with the same output laser energy, the thermal stress pr

57、oduced in glass heated by line-time-shape laser is higher than that produced in glass heated by the other two shapes of lasers.Keywords laser machining, soda-lime glass, finite-element-method (FEM), ANSYS1 IntroductionWith the development of laser technology, many studies have been carried out on cu

58、tting glass with lasers 117. Li et al. 3 put forward a mathematical model to explain the heat transfer of glass heated by a laser beam. Wei et al. 4 and Tian et al. 5 investigated the thermal behavior of glass heated by a CO2-laser beam numerically, and concluded that the resulting temperature distr

59、ibution was strongly dependent on the speed and the parameters of the laser beam. Tsai et al. 6 studied the thermal stress of alumina ceramic substrates irradiated by a moving laser beam and some experiments were carried out to investigate how the crack propagation was influenced by laser power, cut

60、ting speed, and specimen geometry. Glass can be cut by continuous-wave lasers in two different ways. One is the controlled fracture method and the other is melting means. The former has attracted more attention and lots of research has been reported in literatures 716. In contrast, very few studies

61、have been made in detail to investigate cutting glass with the melting method except by Chui 17, due to the low thermal conductivity and the brittleness of the glass material. How to reduce the thermal stress in the glass manufacturing process is a challenging task. Thermal stress is always generated by rapid heating or cooling. If the glass is heated slo