《合金板料激光冲击成形的有限元分析》由会员分享，可在线阅读，更多相关《合金板料激光冲击成形的有限元分析（14页珍藏版）》请在金锄头文库上搜索。
1、newmaker1 引言板料激光冲击成形是近年来出现的一种柔性精密无模成形新方法，它能通过高能激光束对 板料的逐点逐次冲击，积小变形成宏观的大变形。激光冲击成形是在激光冲击强化基础上 发展而来的，用于金属板料塑性成形的一种新技术，它是利用激光诱导的冲击波压力作为 板料塑性成形的变形力，从而实现金属板料的宏观塑性变形。与激光冲击强化相比，约束 的解除为板料的成形提供了空间。图 1 为实验装置示意图。L赧光脉沖龙高濾裔压等務子体3.均竦泾收恳E.睡播费的冲击玻鬼工件图1叛料激光冲击战形示意图金属板料激光冲击波成形涉及激光与材料的相互作用、冲击波传播及其对材料的加载机制、 动态塑性成形理论以及冲压成
2、形工艺等众多学科。受到激光器等硬件条件的制约，目前国 内能开展激光冲击成形试验的研究机构尚且不多，并且试验费用也是不小的开支，而冲击 成形又是一个极其复杂的过程，存在很多不确定性因素。以上原因大大限制了对于激光冲 击成形技术的研究应用。在这种情形下，采用计算机数值模拟技术来获得板料变形与激波 加载之间的规律不失为一种行之有效的方法。因此本章将借助于有限元分析软件对金属板 料激光冲击成形过程进行模拟，为将来建立变形量与激光能量等参数之间的关系以及大面 积冲击成形的理论研究和软件编制提供依据。针对本人前期所做的Al-Mg-Sc材料的激光冲击成形实验,利用有限元模拟软件MSC. Marc 对此过程进
3、行模拟，以便于将模拟结果和试验结果进行比较，优化激光参数。2 金属板料激光冲击成形模拟过程激光冲击成形模拟前要处理的是:建立模型,将所求问题图形化:指定单元类型;给出材料数 据;对模型进行加载;设置边界条件;选择适当的求解器并设置输出要求;求解;查看计算结果。2.1 建立模型激光冲击成形所用的模型为约束模型,即激光辐照能量转换体生成等离子体并发生爆炸产 生强冲击波压力作用于金属板料表面,使金属板料产生宏观塑性变形,这一过程比较复杂 涉及能量转换体如何吸收激光能量、等离子体怎样产生和爆炸、冲击压力波的传播以及金 属板料产生宏观塑性变形等等。由于金属板料的成形与作用在其上的冲击波的压力和作用 时间
4、直接相关,其中,冲击波的压力是与激光参数、能量吸收层、约束层等因素决定,因 此在对金属板料激光冲击成形进行模拟时,不再对冲击波产生和传播进行模拟,而是根据 有关冲击波理论和已有的研究,将激光诱导的冲击波简化为作用在冲击区域内随时间变化 压力载荷,然后将其作为载荷直接作用在金属板料表面。图 2 为激光冲击模型示意图。图 2 激光冲击模型示意图2.2 网格划分及单元选取对于激光冲击成形过程的有限元模拟来说，其建模过程不是很复杂，直接使用Marc软件 的前处理器MSC. Mentat进行网格划分。由于在激光冲击成形过程中所用的金属板料为薄 板，其厚度远小于另外两维尺寸，因此在模拟过程中考虑采用壳体单
5、元。另外，考虑到激 光冲击成形本身的特点，非线性（包括几何非线性和材料非线性）、大变形、而且材料的应 变率非常高，而且壳体单元的节点数应该与网格相对应，对这些因素综合考虑后，在选择单元类型时选用S4R。单元大小为0.05mm*0.05mm,单元类型选取140号壳单元，它是 一种四节点厚壳单元，以整体位移和旋转为自由度，坐标值、平移和转动都采用双线性插 值法。这种壳几何形状采用MITC4壳的几何，采用了和物理稳定方法结合的ANS（假设自 然应变）法，构建的单元就是缩减积分单元，这种平面单元没有对翘曲进行任何人工修整。 这些方法提高了一个积分点单元的精度而没有影响计算速度，可进行大的非线性分析。2
6、.3 边界条件定义 激光冲击成形过程中边界条件的处理对整个分析过程来说是比较重要的一环。这在多点， 多次冲击的时候尤为重要。实际在激光冲击成形过程中，边缘板料沿板料面内的自由度没 有限制，压边圈的存在只是限制了板厚方向的自由度，由于单次激光冲击下板料的变形量 在mm级，经过试验研究发现板料变形后的尺寸在径向方向很小，这样在模拟的时候就可 以对边界条件进行简化处理，既可以将压边圈压住的部分自由度全部限制住.这对最后的变 形结果影响不大。板料固定在凹模基座上，板料在其平面内的转动、移动都没有受到限制， 因此在进行模拟过程中对边界条件进行处理时限制除板料平面内的其他自由度如图 3所示 为激光冲击成形
7、模型力与位移的边界条件定义。图 3 激光冲击边界条件定义2.4 载荷工况定义 在建立金属板料激光冲击成形模型时，激光被直接简化成随时间变化的压力载荷。在激光 冲击成形过程中主要是对冲击载荷和加载时间定义。但是在模拟过程中，确定压力载荷值 大小、载荷作用时间以及加载方式对模拟结果都很重要。2.4.1 激光冲击峰值载荷及加载方式选取激光冲击峰值载荷通过在线测量最为准确。通过压力测量方法得到峰值压力，它是利用VISAR系统得到的，其主要原理基于氢气探针探测受压箔片背面产生的多普勒位移。在L.U.L.I实验室的0. 6ns的测量压力是利用PVDF压电传感器和EMV系统得到的，和PVDF 测量不同的是V
8、ISAR和EMV测量包括自由速度测量。根据自由速度值，利用经典的兰金 雨贡纽关系(Rankin eHugo ni ot)(公式和公式(2)可得到等离子体的压力水平：ZHFL 丫二池蛊(G + SITUU+2/3仙(2) 其中HEL为Hugon血弹性极限(冲击屈鍛强度人 -怕d-氏)* cuU为粒于速度，当 Hugoniot曲线的冲啬和松建部分接近对隸时U = 为戦牲极跟.峰值压力的计算公式P街二咏乏訐蕊门式中各个系数参都是已知的，可以估算出峰值压力的大小，然后根据激光脉冲的作用时间 和实验记录的脉冲信号，将诱导的冲击波按照激光脉冲信号的分布，采取分段逐次逼近的 办法来实现激光脉冲的加载。这里需
9、要注意的是冲击波在整个作用时间内，是随着时间的 变化而变化的，先是上升，随后衰减，呈准高斯分布形式。可以借助Fortran语言，来编 写得到的载荷曲线如图5所示，激光脉宽20ns.激光冲击波作用时间为60ns。圈5锻光冲击随时闾翌就的加我脂浚2.4.2 激光冲击波时空分布 从现已发表的关于激光冲击成形有限元仿真的文章来看，都未考虑到激光冲击波的空间分 布，很显然这种简化导致的误差是不可忽视的。美国哥伦比亚大学的 Y. Lawrence Yao 在 R. Fabbro 模型的基础上给出了激光冲击强化中压力沿光斑径向的变化公式(3)，如下P2、=巩应冲(-気式中，R2是激光光斑的半径，r为离激光束
10、中心的距离。从式(3)可以看出，冲击载荷是关于距离冲击中心半径的函数，考虑到激光脉冲呈高斯分布 将整个载荷作用区域(即光斑直径内的区域)离散为多个带状区域，各区域之间的载荷按高 斯分布加载，而带状区域之内将载荷视为均布载荷。显然，带宽越小越接近于实际的加载状况.当带宽区域趋于 0 时，在整个光斑直径内载荷呈 高斯分布，与实际情况完全吻合。综合考虑计算精度、计算效率和加载难易程度，将半径 为 3. 5mm 的加载区域以 0. 5mm 的间隔离散为 7 个环形带，计算每个带内的平均载荷。 计算过程见附录。图6 显示了按高斯分布离散后的加载的带状云图，可以看出，距激光冲 击中心越近，载荷越大，呈高斯
11、分布。團$驚故北的沖击諏荷；S叭已训!冲托”讪4歸協胡述二丄&卜二：&. 75? I :.；在模拟过程中，还要考虑的是激光诱导的冲击波压力的作用时间，在已有的理论和实验研 究中，激光冲击诱导的冲击波压力的作用时间一直是讨论的主题。这是因为激光脉冲作用 的时间很短，只有几十纳秒，而且在整个激光脉冲的作用过程中，上升沿时间和下降沿时 间又不同，激光脉冲诱导产生的冲击波究竟在板料表面上作用多少时间，有许多学者进行 了研究。R.Fabbro等人的研究认为，激光诱导的冲击波的作用时间大约为激光脉宽的2-3 倍由于脉冲宽度为 22ns ，洪听在其研究中也发现冲击波的时间宽度是激光脉宽的3 倍左右 甚至更高
12、。因此仿真时加载时间设置为70ns。在Marc中，可以用一张表(Table)来控制加 载历程。2.5 材料参数定义 利用有限元分析软件进行材料的成形仿真时，本构关系的建立是至关重要的一步。模拟中 材料参数的选择必须能够反映真实情况，在金属板料激光冲击成形过程中，金属板料发生塑性变形，而且应变率极高，同时材料还会出现硬化效应。因此在进行模拟时，材料参数 的设定既要能够反映板料的变形特点，又要考虑到高应变率下材料的动态屈服应力和硬化 效应对变形的影响。在用Marc做弹塑性分析时弹性变形阶段由弹性模量控制，所以本构关系中只考虑塑性变 形的部分。粘塑性本构方程在一定范围内较好的表达了材料的粘塑性本构特
13、性，但是对 于实际问题，最好通过材料的动态性能试验得到具体的粘一塑性本构方程，这样使计算具 有较高的精度。但由干激光冲击成形的应变速率极高(达106107s-1),目前见到所做的试 验都在 l04s-1 量级以下，远达不到要求的应变速率值。受试验条件的限制，通过插值法， 在有限元仿真时，本章使用应变率 103104s-1 的本构模型，一样能得出在激光冲击成形 高应变率下的数值模拟结果趋势。Marc中自带Johson-Cook模型，而且模型中各个参数 我们已经通过静态下拉伸实验以及动态压缩实验求得，它们分别为A=314. 6MPa B=223, C= 0. 0051, n= 0.1220,假设材
14、料为各项同性。2.6 接触考虑到激光冲击的实际情况，板料放在凹模上，再用压板将板料压紧。可以用定义接触体 的方法来很好地模拟实际冲击时的约束情况，并且建立了接触表来控制接触体之间的接触 类型。如图7所示，定义了三个接触体:sheet为板料，定义为变形体;holddown为压板， 定义为刚体，和板料之间的接触关系为 touching;mold 为凹模，定义为刚体，和板料之间 的接触关系为glue(粘合)。3 激光冲击成形数值模拟结果3.1 激光冲击成形过程利用MSC. Marc的后处理器MSC. Men tat进行后处理，将分析生成的结果以图形、曲线 和动画的形式展现出来，便于和实验结果进行比较
15、。以下都是单点单次激光冲击成形数值 模拟结果，图8以带状云图(CONTOUR BAND)形式显示了激光冲击板料变形过程，厚度 为1.2mm的AI-Mg-Sc合金板在激光冲击载荷下变形量逐渐的变化。在Marc中变形过程 设置成一百步，这里选取了四步来说明其过程的变化。图 9 显示了激光冲击板料的最终变 形量为 1.50lmm。Y 向位移沿路径的变化曲线如图 10 所示。3. 2 激光能量与变形量的关系 为研究激光能量与变形量的关系，我们选取一组激光能量数据来模拟变形量随激光能量变化的过程。图11为激光能量分别为40. 030J, 36. 023J,30. 690J,25. 375J的激光对厚度为1.2mm 的 A1-Mg-Sc 进行冲击后的变形模拟图，可以看出随着激光能量的增大，板料变形 量增大。表1激址礎对板料变形的骼响试萍 犀匱 约東孔艳 薇光逾童毆大变形撻模扭值（严1）疑大更形量的第验值5说煨号（TOTB）mm)(J)01 1.23040.03ftL5O1140002 I.Z303& 0231801350&31.23030,6901.232115004 L 230Z5. 3?50.955875表1 为模拟结果和实验结果的对比，从表中数据可以看出模拟结果均比实验数据略有差距
