版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
覆盖件冲压成形数值模拟与回弹控制的深度解析与实践一、引言1.1研究背景与意义在现代工业生产中,覆盖件冲压成形技术扮演着举足轻重的角色,被广泛应用于汽车、航空航天、电子等众多领域。以汽车行业为例,汽车覆盖件如车身外板、车门、发动机罩等,其重量占到汽车用钢材总量的50%以上,这些覆盖件不仅构成了汽车的外观和基本结构,还对汽车的安全性、舒适性和美观性起着关键作用。航空航天领域的飞行器外壳、内饰件等覆盖件,同样对飞行器的性能和安全性至关重要。在电子产品领域,如手机、电脑等设备的外壳,也多采用冲压成形的覆盖件,其质量直接影响产品的外观和用户体验。冲压成形作为覆盖件制造的主要方式,具有生产效率高、成本低、能够制造形状复杂零件等显著优点。然而,在实际冲压成形过程中,板料的变形行为极为复杂,涉及到几何非线性、材料非线性以及接触非线性等多重非线性问题。这使得冲压件的质量控制面临诸多挑战,其中回弹问题尤为突出。回弹是指在冲压成形过程结束后,当载荷卸去时,冲压件由于弹性恢复而产生的形状和尺寸与模具型面不一致的现象。回弹现象的存在严重影响了冲压件的尺寸精度和形状精度,导致产品质量下降,增加了模具调试和产品后续加工的难度,甚至可能造成零件报废,给企业带来巨大的经济损失。据统计,仅在美国汽车行业,每年由于成形回弹而造成的损失就超过五千万美元。为弥补回弹造成的成形偏差,需调整工艺或修正模具,这项工作会增加制造准备时间,约占制造准备时间的20%-30%。传统的冲压工艺设计和模具开发主要依赖于经验和反复试模,这种方式不仅周期长、成本高,而且难以准确预测和控制冲压件的回弹。随着计算机技术和数值计算方法的飞速发展,数值模拟技术应运而生,为解决冲压成形中的难题提供了新的途径。通过数值模拟,可以在计算机上对冲压成形过程进行虚拟仿真,深入分析板料在冲压过程中的应力、应变分布以及材料的流动规律,准确预测可能出现的缺陷,如起皱、破裂和回弹等。这使得工程师能够在模具制造之前对冲压工艺进行优化,减少试模次数,降低生产成本,提高产品质量和生产效率。在当前制造业竞争日益激烈的背景下,对覆盖件冲压成形进行数值模拟和回弹研究具有极其重要的现实意义。通过深入研究板料冲压成形的数值模拟方法和回弹规律,可以为工业生产提供更加科学、准确的理论依据和技术支持,有效提升冲压件的质量和生产效率,增强企业的市场竞争力,推动相关产业的高质量发展。1.2国内外研究现状随着计算机技术和有限元理论的不断发展,国内外学者在覆盖件冲压成形数值模拟和回弹研究方面取得了丰硕的成果。在数值模拟方面,国外起步较早,发展较为成熟。自20世纪70年代起,有限元法开始应用于板料成形过程的模拟,随后各种成熟的CAE仿真软件纷纷应用于模具成形产业中。美国ETA公司开发的Dynaform软件,是目前板料冲压成形模拟领域应用最广泛的软件之一,它能够对复杂的冲压成形过程进行精确模拟,分析板料的变形、应力应变分布等情况。德国的AutoForm软件在汽车覆盖件冲压模拟中也具有重要地位,该软件具有强大的前后处理功能和高精度的求解器,能够快速准确地预测冲压成形过程中的各种问题。此外,法国的PAM-STAMP软件在航空航天等领域的覆盖件冲压模拟中表现出色,其独特的算法和材料模型能够满足不同材料和复杂形状覆盖件的模拟需求。国内在数值模拟领域的研究虽然起步相对较晚,但发展迅速。近年来,许多高校和科研机构开展了相关研究工作,并取得了一系列重要成果。华中科技大学在基于特征的覆盖件模型上,通过特征映射与分解产生一组彼此相关联的工艺特征,再对这些工艺特征进行交互的组合和排序,从而得到工艺方案,为计算机辅助覆盖件工艺方案设计提供了新的思路。上海交通大学针对汽车覆盖件冲压成形过程中的复杂问题,开发了一系列数值模拟算法和软件模块,提高了模拟的精度和效率。同时,国内企业也逐渐认识到数值模拟技术的重要性,开始加大在这方面的投入和应用,如一汽、上汽等汽车制造企业,通过使用数值模拟技术优化冲压工艺,降低了生产成本,提高了产品质量。在回弹研究方面,国内外学者进行了大量的理论、实验和数值模拟研究。国外学者在回弹预测模型和控制方法方面取得了许多重要进展。例如,采用静态隐式算法和动-静态联合算法来仿真汽车覆盖件成形和切边时的回弹,能够更准确地预测回弹量。通过优化模具型面、调整工艺参数等方法来控制回弹,取得了较好的效果。国内学者也在回弹研究领域取得了显著成果。山东大学使用有限元模拟软件DYNAFORM,结合弹塑性有限元方法对覆盖件的回弹规律进行了研究,分析了压边力、摩擦状态、拉延筋、模具间隙等因素对覆盖件成形及回弹的影响,研究结果表明回弹量随凸凹模间隙的增大而增大,随压边力、摩擦系数、拉延筋阻力的增大而减小。南京工业大学以某型汽车底座横梁外板为研究对象,采用板料成形仿真软件Dynaform,在计算机上模拟冲压成形及回弹的过程,预测实际生产过程中可能出现的回弹量,并总结出控制回弹的措施。然而,目前的研究仍存在一些不足之处。在数值模拟方面,虽然各种模拟软件能够对冲压成形过程进行较为准确的模拟,但对于一些复杂的冲压工艺和材料模型,模拟精度仍有待提高。不同软件之间的数据交换和协同工作能力也有待加强,以满足多学科、多领域的联合仿真需求。在回弹研究方面,虽然已经提出了多种回弹预测模型和控制方法,但由于回弹影响因素众多,且各因素之间相互耦合,导致回弹的准确预测和有效控制仍然是一个难题。此外,现有的回弹控制方法往往需要增加模具制造和调试成本,如何在保证控制效果的前提下降低成本,也是需要进一步研究的问题。1.3研究内容与方法本研究旨在深入探究覆盖件冲压成形的数值模拟方法,精准分析回弹规律,并提出切实可行的控制策略,主要内容如下:覆盖件冲压成形数值模拟方法研究:全面剖析冲压成形过程中的多重非线性问题,包括几何非线性、材料非线性和接触非线性的具体表现和作用机制。深入研究有限元理论在冲压成形数值模拟中的应用,包括单元类型选择、网格划分技术、求解算法优化等,以提高模拟的精度和效率。对常用的冲压成形模拟软件,如Dynaform、AutoForm、PAM-STAMP等进行对比分析,明确各软件的特点、优势和适用范围,为后续模拟工作的软件选择提供依据。覆盖件冲压成形回弹规律研究:以典型覆盖件为研究对象,运用数值模拟和实验相结合的方法,系统研究压边力、摩擦状态、拉延筋、模具间隙等工艺参数对回弹的影响规律。通过改变压边力大小,观察回弹量的变化情况,分析压边力与回弹之间的定量关系;研究不同摩擦系数下板料与模具之间的摩擦力变化,以及对回弹的影响机制;探讨拉延筋的设置位置、形状和尺寸对板料流动和回弹的影响;分析模具间隙的大小对板料成形和回弹的影响规律。基于大量的模拟和实验数据,建立考虑多因素影响的回弹预测模型,提高回弹预测的准确性。利用机器学习算法,如神经网络、支持向量机等,对回弹数据进行训练和学习,构建高精度的回弹预测模型。覆盖件冲压成形回弹控制策略研究:根据回弹规律研究结果,提出优化工艺参数、改进模具结构等回弹控制策略。通过数值模拟和实验验证,评估各控制策略的有效性。优化工艺参数方面,采用正交试验设计或响应面法,对压边力、拉延筋阻力、冲压速度等工艺参数进行优化组合,以减小回弹量;改进模具结构方面,研究采用变压边力模具、数控柔性拉边系统等新型模具结构,实现对回弹的有效控制。探索模具型面补偿法在覆盖件回弹控制中的应用,通过对模具型面进行反向修正,补偿回弹变形,提高冲压件的尺寸精度。利用数值模拟技术,对模具型面补偿量进行优化计算,确定最佳的补偿方案。本研究将综合运用多种研究方法,以确保研究的全面性和深入性:有限元分析方法:利用成熟的有限元分析软件,如Dynaform、AutoForm等,对覆盖件冲压成形过程进行数值模拟。通过建立合理的有限元模型,模拟板料在冲压过程中的应力、应变分布,预测可能出现的起皱、破裂和回弹等缺陷,为工艺参数优化和模具设计提供理论依据。在模拟过程中,精确设置材料参数、接触条件、边界条件等,确保模拟结果的准确性。实验研究方法:设计并进行覆盖件冲压成形实验,采用与实际生产相同或相似的材料、模具和设备,获取真实的冲压件数据。通过对实验结果的分析,验证数值模拟的准确性,深入研究冲压成形过程中的物理现象和回弹规律。实验过程中,严格控制实验条件,多次重复实验,以确保实验数据的可靠性。理论分析方法:基于塑性力学、材料力学等相关理论,对覆盖件冲压成形过程中的力学行为进行理论分析。推导应力应变关系、变形协调方程等,深入理解板料的变形机理和回弹产生的原因,为数值模拟和实验研究提供理论支持。结合理论分析,对模拟结果和实验数据进行深入解读,揭示冲压成形过程中的内在规律。二、覆盖件冲压成形数值模拟理论基础2.1弹塑性有限元法2.1.1基本原理弹塑性有限元法是覆盖件冲压成形数值模拟的核心方法,其基本原理基于连续介质力学和变分原理。在冲压成形过程中,板料经历大变形,涉及几何非线性、材料非线性和接触非线性,传统解析方法难以求解,而弹塑性有限元法通过离散化处理,将连续的板料划分成有限个单元,单元之间通过节点连接,从而将复杂的连续体问题转化为有限个单元的集合问题进行求解。在板材成形数值模拟领域,弹塑性有限元法采用Lagrange描述,即以初始时刻物体内各质点的坐标作为描述物体变形的参考坐标。在冲压过程中,随着板料的变形,各质点的位置发生变化,但始终以初始坐标为基准进行描述。这种描述方式能方便跟踪材料质点的运动轨迹,精确计算工件的变形情况。例如,在汽车覆盖件冲压成形过程中,通过Lagrange描述可以清晰看到板料从初始平板状态逐渐变形为复杂曲面形状的过程,准确捕捉材料在各个区域的流动和变形情况。弹塑性有限元法不仅能计算工件的变形,还能有效计算应力、应变分布。在冲压成形中,板料各部分受到不同的外力作用,产生复杂的应力应变状态。通过弹塑性有限元法,基于材料的本构关系(描述材料应力与应变之间关系的方程),可以计算出每个单元在不同时刻的应力和应变值。如在分析汽车发动机罩冲压成形时,能够得到板料在拉延、翻边等工序中的应力应变分布云图,直观展示危险区域,为工艺优化提供依据。此外,该方法还能计算工件回弹和残余应力。回弹是覆盖件冲压成形中的关键问题,直接影响产品尺寸精度。弹塑性有限元法通过考虑材料的弹性卸载过程,模拟冲压结束后板料的回弹行为,预测回弹量。残余应力是冲压过程中材料内部残留的应力,可能导致零件变形、开裂等问题,利用弹塑性有限元法可以计算残余应力的大小和分布,评估其对零件质量的影响。2.1.2关键方程与算法在弹塑性有限元法中,虚功原理方程是重要的基础方程。在汽车覆盖件冲压成形过程中,坯料与模具所构成的大变形接触冲压系统都须满足弱形式的虚功原理。对于一个处于平衡状态的力学系统,虚功原理可表述为:外力在虚位移上所做的虚功等于内力在虚位移上所做的虚功。用数学表达式表示为:\int_{\Omega}\sigma_{ij}\delta\epsilon_{ij}d\Omega=\int_{\Gamma_{t}}\bar{t}_{i}\deltau_{i}d\Gamma+\int_{\Omega}\bar{b}_{i}\deltau_{i}d\Omega其中,\sigma_{ij}是应力张量,\delta\epsilon_{ij}是虚应变张量,\Omega是物体的体积域,\Gamma_{t}是给定面力的边界,\bar{t}_{i}是给定的面力,\deltau_{i}是虚位移,\bar{b}_{i}是体积力。该方程将力学系统的平衡条件与虚位移联系起来,通过求解此方程可以得到物体的位移、应力等物理量。在覆盖件冲压分析中,常用的算法之一是动态显式算法。在显式求解法中,用中心差分法进行动态问题的时域积分。其基本思想是将时间域离散化,通过逐步求解每个时间步的运动方程来获得系统的响应。以某汽车覆盖件冲压成形模拟为例,在每个时间步,根据当前时刻板料的位移、速度和所受外力,利用中心差分法计算下一时刻板料的位移和速度。动态显式算法的优点是无条件收敛,不需要求解联立方程组,计算过程相对简单,对于处理大变形、接触碰撞等复杂问题具有优势。在汽车覆盖件冲压过程中,板料与模具之间存在复杂的接触和摩擦,动态显式算法能够较好地处理这些非线性因素,准确模拟板料的变形过程。然而,该算法的时间步长受到数值积分稳定性的限制,通常需要取较小的时间步长,导致计算时间较长。为了提高计算效率,在实际应用中,一般会采取加大冲压速度和板料质量密度的方法,但这需要在计算精度和计算效率之间进行权衡。2.2数值模拟中的接触与摩擦处理2.2.1接触算法在覆盖件冲压成形的数值模拟中,坯料与模具之间的接触行为是一个关键问题,直接影响模拟结果的准确性。目前,常用的接触算法主要有罚函数法、拉格朗日乘子法和增广拉格朗日法等,其中罚函数法应用较为广泛。罚函数法的基本原理是通过引入一个罚因子,将接触约束条件转化为惩罚项添加到系统的能量泛函中。在冲压过程中,当坯料与模具之间的距离小于设定的接触容差时,认为两者发生接触,罚函数会产生一个较大的接触力,以阻止坯料穿透模具表面。假设接触力为F_c,罚因子为\alpha,穿透深度为\delta,则接触力可表示为F_c=\alpha\delta。罚因子\alpha的大小决定了接触力的大小,若罚因子过大,可能会导致计算结果不稳定;若罚因子过小,则可能无法有效防止穿透。在实际应用中,需要根据具体问题进行合理取值。以汽车覆盖件冲压成形为例,在Dynaform软件中,常采用罚函数法来处理坯料与模具的接触。在模拟汽车发动机罩的拉延过程时,通过设置罚因子,能够准确模拟坯料与模具之间的接触状态,得到坯料在不同时刻的应力应变分布情况。罚函数法的优点是算法简单,计算效率高,易于实现,在大多数情况下能够满足工程计算的精度要求。然而,它也存在一些局限性,如无法精确满足接触约束条件,可能会导致一定的穿透误差。拉格朗日乘子法是通过引入拉格朗日乘子来强制满足接触约束条件,能够精确处理接触问题,但计算过程较为复杂,需要求解额外的拉格朗日乘子方程,计算效率相对较低。增广拉格朗日法结合了罚函数法和拉格朗日乘子法的优点,通过在罚函数法的基础上增加拉格朗日乘子项,既能提高计算精度,又能在一定程度上保持计算效率,在一些对接触精度要求较高的场合得到应用。在实际数值模拟中,应根据具体的冲压工艺和计算精度要求,选择合适的接触算法。2.2.2摩擦模型摩擦在覆盖件冲压成形过程中起着重要作用,它影响着板料的流动、应力应变分布以及成形质量。在数值模拟中,准确描述摩擦行为对于提高模拟精度至关重要。库仑摩擦模型是覆盖件冲压成形数值模拟中应用最广泛的摩擦模型之一。库仑摩擦模型基于库仑定律,认为摩擦力与接触面上的正压力成正比,其表达式为:F_f=\muF_n其中,F_f是摩擦力,\mu是摩擦系数,F_n是接触面上的正压力。在冲压成形过程中,板料与模具之间的摩擦系数\mu受到多种因素的影响,如板料和模具的表面粗糙度、润滑条件、材料特性等。不同的冲压工艺和材料组合,其摩擦系数会有较大差异。例如,在汽车覆盖件冲压中,对于表面光滑且经过良好润滑处理的模具和板料,摩擦系数可能在0.05-0.15之间;而对于表面粗糙度较大且润滑不足的情况,摩擦系数可能会达到0.2-0.3。在实际应用库仑摩擦模型时,准确设置摩擦系数是关键。如果摩擦系数设置过小,板料在模具表面的滑动阻力较小,可能会导致板料流动过快,出现起皱等缺陷;反之,如果摩擦系数设置过大,板料的流动受到过度限制,可能会导致局部应力集中,增加破裂的风险。为了确定合理的摩擦系数,通常需要结合实验测试和经验数据。通过进行板料与模具材料的摩擦试验,测量不同工况下的摩擦力,从而得到准确的摩擦系数。在一些研究中,利用圆环镦粗试验来测量板材与模具之间的摩擦系数,为数值模拟提供可靠的参数依据。然而,库仑摩擦模型也存在一定的局限性。它假设摩擦系数是常数,不考虑摩擦系数随接触状态、相对滑动速度等因素的变化。在实际冲压过程中,摩擦系数往往是动态变化的,尤其是在高速冲压或复杂变形情况下,库仑摩擦模型的准确性会受到影响。为了更准确地描述摩擦行为,一些改进的摩擦模型,如考虑速度效应的摩擦模型、基于微凸体接触理论的摩擦模型等也逐渐被提出并应用于数值模拟中。这些模型能够更好地反映实际冲压过程中的摩擦现象,但计算复杂度相对较高,在实际应用中需要根据具体情况选择合适的摩擦模型。三、覆盖件冲压成形数值模拟方法与应用3.1数值模拟流程与前处理技术3.1.1模型建立在覆盖件冲压成形数值模拟中,模型建立是首要且关键的环节。以汽车覆盖件为例,通常使用专业的CAD(计算机辅助设计)软件,如UG、CATIA等进行三维模型的构建。这些软件具备强大的曲面造型和实体建模功能,能够精确地描绘出汽车覆盖件复杂的几何形状。在构建汽车发动机罩的三维模型时,设计师需要充分考虑发动机罩的外形轮廓、曲率变化、加强筋的布局以及与车身其他部件的装配关系等因素。通过CAD软件的各种建模工具,如拉伸、旋转、扫描、布尔运算等,将设计理念转化为精确的三维数字模型。完成三维模型构建后,需将其导入到CAE(计算机辅助工程)软件中,如Dynaform、AutoForm等,以便进行后续的数值模拟分析。在导入过程中,由于不同软件之间的数据格式存在差异,可能会出现数据丢失、模型变形等问题。因此,需要选择合适的数据接口和转换方式,确保模型的完整性和准确性。一般来说,常用的CAD数据格式有IGES、STEP等,这些格式具有较好的通用性和兼容性,能够在不同软件之间实现数据的有效传输。在将UG创建的汽车覆盖件三维模型导入Dynaform软件时,可先将模型保存为IGES格式,然后在Dynaform中通过相应的导入功能读取该文件。导入后,还需对模型进行检查和修复,如检查模型是否存在破面、重叠面、缝隙等缺陷,对发现的问题及时进行修补,以保证后续模拟分析的顺利进行。3.1.2网格划分网格划分是将连续的覆盖件模型离散化为有限个单元的过程,其质量直接影响数值模拟的精度和计算效率。对于覆盖件这种形状复杂的零部件,传统的规则网格划分方法难以满足需求,因此常采用自适应网格划分技术。自适应网格划分技术能够根据模型的几何特征和变形情况,自动调整网格的疏密程度。在覆盖件冲压成形过程中,板料的变形通常是不均匀的,某些区域如圆角、拐角处以及拉深深度较大的部位,变形较为剧烈,应力应变梯度大;而其他相对平坦的区域,变形则较为平缓。自适应网格划分技术能够在变形剧烈的区域自动加密网格,提高计算精度;在变形平缓的区域适当降低网格密度,减少计算量,从而在保证计算精度的前提下,提高计算效率。以汽车覆盖件中的车门内板冲压成形模拟为例,在初始状态下,对车门内板模型进行初步的网格划分,设置一个相对均匀的初始网格尺寸。在模拟计算过程中,软件会根据预先设定的误差准则,实时监测板料各部分的变形情况。当发现某些区域的计算误差超过允许范围时,判断该区域变形剧烈,需要细化网格。软件会自动对这些区域的单元进行细分,如将一个大的四边形单元分割为四个小的四边形单元,使网格更加细密,从而更准确地捕捉材料的变形行为。自适应网格划分技术对模拟精度的提升具有显著作用。通过在关键区域加密网格,可以更精确地计算应力、应变分布,减少数值计算误差,从而更准确地预测冲压件可能出现的起皱、破裂等缺陷。在分析汽车覆盖件的拉深成形过程时,采用自适应网格划分技术,能够清晰地显示出板料在圆角处的应力集中情况,以及由于材料流动不均匀导致的潜在起皱区域,为工艺优化提供更可靠的依据。然而,自适应网格划分技术也会增加一定的计算时间和内存消耗,因为在模拟过程中需要不断地对网格进行调整和更新。因此,在实际应用中,需要根据具体问题的复杂程度和计算资源的限制,合理设置自适应网格划分的参数,以平衡计算精度和计算效率之间的关系。3.1.3工艺参数设置工艺参数的合理设置是覆盖件冲压成形数值模拟的重要环节,直接影响模拟结果的准确性和可靠性。冲压速度和压边力是两个关键的工艺参数。冲压速度对覆盖件冲压成形过程有着多方面的影响。在实际冲压生产中,冲压速度通常在一定范围内变化,一般汽车覆盖件冲压速度在0.1-10m/s之间。冲压速度过快,会使板料在短时间内受到较大的冲击力,导致材料内部的应力应变分布不均匀,增加破裂的风险。当冲压速度过高时,板料与模具之间的摩擦生热来不及散失,会使局部温度升高,材料的力学性能发生变化,进一步影响成形质量。冲压速度过慢,则会降低生产效率,增加生产成本。在数值模拟中,设置冲压速度时需要综合考虑材料特性、模具结构以及生产效率等因素。对于塑性较好的材料,可以适当提高冲压速度;而对于塑性较差、对变形速率敏感的材料,则需要选择较低的冲压速度。压边力是控制板料在冲压过程中流动的重要参数。其主要作用是防止板料在拉深过程中起皱,通过施加一定的压力,限制板料的径向流动,使板料在凸模的作用下均匀变形。压边力的大小需要根据覆盖件的形状、尺寸、材料性能以及模具结构等因素进行合理确定。如果压边力过小,无法有效抑制板料的起皱,导致冲压件表面出现皱纹,影响产品质量;压边力过大,则会使板料受到过度的挤压,增加破裂的可能性,同时也会增大模具的磨损和设备的负荷。在模拟汽车覆盖件的拉深成形时,通常需要通过多次试算,结合经验公式和实际生产数据,确定合适的压边力范围。对于形状复杂、深度较大的覆盖件,可能需要采用变压边力的方式,即在冲压过程中根据板料的变形情况实时调整压边力的大小,以实现更好的成形效果。除了冲压速度和压边力外,还有其他一些工艺参数,如摩擦系数、模具间隙、拉延筋的设置等,也会对覆盖件冲压成形过程产生重要影响。在数值模拟中,需要准确设置这些工艺参数,以确保模拟结果能够真实反映实际冲压过程,为工艺优化和模具设计提供可靠的依据。3.2数值模拟软件及应用案例3.2.1常用模拟软件介绍DYNAFORM是一款广泛应用于板料冲压成形模拟的专业软件,由ETA公司开发。该软件具有强大的前后处理功能和高效的求解器,能够对复杂的冲压成形过程进行精确模拟。在前后处理方面,DYNAFORM提供了丰富的工具和选项,方便用户进行模型导入、网格划分、工艺参数设置等操作。它支持多种CAD数据格式的导入,如IGES、VDA、UG和CATIA等,能够实现与其他设计软件的无缝对接,确保模型的完整性和准确性。在网格划分方面,DYNAFORM拥有自动网格划分和自适应网格划分功能,能够根据模型的几何特征和变形情况,自动生成高质量的网格,提高计算精度和效率。DYNAFORM的求解器基于著名的LS-DYNA求解器开发,采用动态显式算法,能够快速准确地模拟板料在冲压过程中的大变形、接触和摩擦等复杂问题。在汽车覆盖件冲压成形模拟中,DYNAFORM可以准确预测板料的破裂、起皱、减薄、划痕和回弹等缺陷,为工艺设计和模具优化提供重要依据。通过模拟分析,工程师可以提前发现潜在的问题,并通过调整工艺参数或模具结构来解决这些问题,从而减少试模次数,降低生产成本,提高产品质量。ABAQUS是一款功能强大的通用有限元分析软件,在覆盖件冲压成形模拟中也有广泛应用。与DYNAFORM相比,ABAQUS具有更广泛的材料模型库和更灵活的求解器选项,能够处理各种复杂的工程问题。ABAQUS的材料模型库涵盖了金属、塑料、橡胶等多种材料,并且支持用户自定义材料模型,能够满足不同材料在冲压成形过程中的模拟需求。在求解器方面,ABAQUS提供了显式求解器和隐式求解器,用户可以根据问题的特点和计算效率的要求选择合适的求解器。显式求解器适用于处理大变形、高速碰撞等动态问题,而隐式求解器则适用于处理静态和准静态问题,具有较高的计算精度。在覆盖件冲压成形模拟中,ABAQUS能够精确模拟板料的应力应变分布、材料流动规律以及回弹行为。通过对这些模拟结果的分析,工程师可以深入了解冲压成形过程中的物理现象,为工艺优化和模具设计提供科学依据。在分析汽车覆盖件的冲压回弹问题时,ABAQUS可以考虑材料的各向异性、加工硬化等因素,建立准确的回弹预测模型,从而更准确地预测回弹量,为模具型面补偿提供参考。此外,ABAQUS还具有强大的二次开发功能,用户可以根据自己的需求编写脚本程序,扩展软件的功能,满足特殊的工程应用需求。3.2.2案例分析:汽车侧围外板冲压成形模拟以汽车侧围外板为例,展示运用数值模拟软件进行冲压成形模拟的过程。汽车侧围外板是汽车车身的重要覆盖件之一,其形状复杂,尺寸较大,冲压成形过程中容易出现起皱、破裂和回弹等问题。在本次模拟中,选用DYNAFORM软件进行分析。首先,利用CAD软件构建汽车侧围外板的三维模型,并将其导入DYNAFORM软件中。在导入过程中,仔细检查模型是否存在破面、重叠面等问题,确保模型的质量。接着,对模型进行网格划分,采用自适应网格划分技术,根据侧围外板的几何形状和变形特点,在关键区域如圆角、拐角处加密网格,以提高计算精度;在相对平坦的区域适当降低网格密度,减少计算量。完成网格划分后,设置工艺参数,包括冲压速度、压边力、摩擦系数等。根据实际生产经验和相关资料,将冲压速度设置为5m/s,压边力设置为300kN,摩擦系数设置为0.12。同时,定义材料属性,汽车侧围外板材料选用高强度冷轧钢板,其弹性模量为206GPa,泊松比为0.3,屈服强度为340MPa,抗拉强度为450MPa。模拟结果分析表明,在当前工艺参数下,汽车侧围外板冲压成形过程中存在一些问题。通过查看成形极限图(FLD),发现侧围外板的局部区域处于破裂危险区,如A柱附近和车门开口处,这表明这些区域在冲压过程中可能会发生破裂。分析厚度变化云图,发现部分区域的材料减薄率超过了15%,尤其是在拉深深度较大的部位,减薄现象较为严重,这可能会影响侧围外板的强度和刚度。此外,通过回弹分析,发现侧围外板的回弹量较大,最大回弹量达到了3mm,这将严重影响侧围外板的尺寸精度和装配精度。基于模拟结果,对工艺方案进行改进。为了解决破裂问题,在破裂危险区域增加拉延筋,通过调整拉延筋的位置和高度,改变材料的流动方向和速度,使材料更加均匀地流入模具型腔,从而降低局部区域的应力集中,避免破裂的发生。针对材料减薄问题,适当增加压边力,从300kN提高到350kN,以更好地控制材料的流动,减少材料在拉深过程中的减薄。为减小回弹量,优化模具型面,根据回弹模拟结果,对模具型面进行反向补偿,使冲压件在回弹后能够达到设计要求的尺寸精度。通过再次模拟验证,改进后的工艺方案取得了良好的效果。侧围外板的破裂和减薄问题得到了有效解决,成形极限图显示所有区域均处于安全范围内,材料减薄率控制在了10%以内。回弹量也显著减小,最大回弹量降低到了1.5mm,满足了生产要求。通过本次案例分析,充分展示了数值模拟软件在汽车覆盖件冲压成形工艺设计中的重要作用,能够为实际生产提供可靠的技术支持,有效提高冲压件的质量和生产效率。四、覆盖件冲压成形回弹研究4.1回弹的形成机制与影响因素4.1.1回弹的力学原理在覆盖件冲压成形过程中,板料经历了复杂的加载和卸载过程,回弹现象的产生与这一过程中的应力应变变化密切相关。从材料力学的角度来看,冲压成形过程可分为弹性变形阶段、弹塑性变形阶段和塑性变形阶段。在弹性变形阶段,当外力作用于板料时,板料内部产生应力,应力与应变之间遵循胡克定律,即\sigma=E\varepsilon,其中\sigma为应力,E为弹性模量,\varepsilon为应变。此时,板料的变形是可逆的,一旦外力去除,板料将完全恢复到原来的形状和尺寸。当外力继续增加,应力超过材料的屈服强度时,板料进入弹塑性变形阶段。在这个阶段,板料的变形由弹性变形和塑性变形两部分组成。塑性变形是不可逆的,即使外力去除,这部分变形也会保留下来;而弹性变形部分仍遵循胡克定律,外力去除后会发生弹性恢复。随着外力的进一步增大,板料进入塑性变形阶段,此时塑性变形占主导地位,材料的变形主要是不可逆的塑性变形,但弹性变形依然存在。当冲压成形过程结束,外力卸去时,板料进入卸载阶段。在卸载过程中,板料内部的应力逐渐减小,弹性变形部分开始恢复,导致板料的形状和尺寸发生变化,这就是回弹现象的本质原因。由于板料在冲压过程中各部分的受力情况不同,应力应变分布不均匀,因此卸载后的回弹量在不同部位也会有所差异。以简单的弯曲冲压为例,在弯曲过程中,板料的外层受到拉应力作用,发生拉伸变形;内层受到压应力作用,发生压缩变形。当卸载时,外层的弹性拉伸变形恢复,使板料有向外张开的趋势;内层的弹性压缩变形恢复,使板料有向内收缩的趋势。这两种相反的趋势相互作用,导致弯曲件的角度和曲率发生变化,产生回弹现象。在汽车覆盖件的冲压成形中,由于覆盖件形状复杂,涉及到多种变形方式的组合,如拉深、翻边、胀形等,板料的应力应变分布更加复杂,回弹问题也更加突出。4.1.2影响回弹的因素分析覆盖件冲压成形回弹受到多种因素的综合影响,深入研究这些因素对于准确预测和有效控制回弹至关重要。材料性能是影响回弹的关键因素之一。屈服强度和弹性模量对回弹有着显著的影响。屈服强度是材料开始产生塑性变形的临界应力,屈服强度越高,材料在冲压过程中越不容易发生塑性变形,弹性变形所占比例相对较大。当卸载时,弹性恢复的程度就越大,从而导致回弹量增加。高强度钢在汽车覆盖件中的应用越来越广泛,由于其屈服强度较高,冲压后的回弹问题也更为突出。弹性模量是衡量材料抵抗弹性变形能力的指标,弹性模量越大,材料的刚度越大,在相同外力作用下的弹性变形越小。因此,弹性模量较大的材料,其回弹量相对较小。铝合金材料的弹性模量约为钢材的三分之一,在相同冲压条件下,铝合金覆盖件的回弹量通常比钢质覆盖件大。加工硬化指数也会对回弹产生影响。加工硬化指数表示材料在塑性变形过程中强度提高的程度。加工硬化指数较大的材料,在冲压过程中随着塑性变形的增加,材料的强度不断提高,抵抗变形的能力增强。这使得材料在卸载时的弹性恢复相对较小,从而减小了回弹量。模具结构对回弹也有重要影响。模具间隙是指凸模与凹模之间的间隙,模具间隙的大小直接影响板料在冲压过程中的受力状态和变形方式。当模具间隙过大时,板料在冲压过程中受到的约束较小,容易产生较大的弹性变形。卸载后,弹性恢复量增大,导致回弹量增加。相反,模具间隙过小,板料与模具之间的摩擦力增大,会使板料的变形不均匀,局部应力集中,也可能增加回弹的风险。在汽车覆盖件冲压模具设计中,合理控制模具间隙是减小回弹的重要措施之一。拉延筋的设置是控制板料流动的重要手段,拉延筋通过增加板料的流动阻力,使板料在冲压过程中更加均匀地变形。合理设置拉延筋的位置、形状和高度,可以有效地调整板料的应力应变分布,减小回弹量。在汽车侧围外板的冲压成形中,在容易出现回弹的区域设置合适的拉延筋,能够改善板料的流动,使材料在成形过程中受力更加均匀,从而降低回弹量。工艺参数对回弹的影响也不容忽视。压边力是控制板料在冲压过程中流动的重要工艺参数。适当增大压边力,可以增加板料与压边圈之间的摩擦力,限制板料的径向流动,使板料在凸模的作用下更加均匀地变形。这有助于减小板料的局部应力集中,降低回弹量。但压边力过大,会使板料受到过度的挤压,增加破裂的风险。在实际生产中,需要根据覆盖件的形状、尺寸和材料性能等因素,合理调整压边力的大小。冲压速度对回弹也有一定的影响。冲压速度过快,会使板料在短时间内受到较大的冲击力,导致材料内部的应力应变分布不均匀,增加回弹的可能性。高速冲压时,板料与模具之间的摩擦生热来不及散失,会使局部温度升高,材料的力学性能发生变化,进一步影响回弹。因此,在冲压成形过程中,需要选择合适的冲压速度,以保证板料的变形均匀,减小回弹。4.2回弹研究的实验方法与测量技术4.2.1实验方案设计为深入研究覆盖件冲压成形的回弹特性,精心设计了一套全面且严谨的实验方案。在实验材料选择上,充分考虑覆盖件在实际生产中的常用材料,选用了DC04冷轧钢板。DC04是一种具有良好冲压性能的低碳钢,其化学成分和力学性能稳定,在汽车、家电等行业的覆盖件制造中广泛应用。该材料的屈服强度为140-200MPa,抗拉强度为270-370MPa,断后伸长率不小于34%,这些性能参数使其成为研究冲压成形回弹的理想材料。模具设计是实验方案的关键环节,需确保模具能够准确模拟实际冲压过程。采用CAD/CAM技术设计并制造了一套专用模具,包括凸模、凹模和压边圈。模具的型面设计严格按照目标覆盖件的形状进行,尺寸精度控制在±0.05mm以内,以保证冲压件的尺寸精度和形状精度。凸模和凹模的材料选用Cr12MoV,该材料具有高硬度、高强度和良好的耐磨性,能够满足多次冲压实验的要求。在模具制造过程中,运用数控加工技术进行精密加工,确保模具表面的粗糙度达到Ra0.8μm以下,减少模具表面粗糙度对板料冲压成形和回弹的影响。实验步骤安排如下:首先,将DC04冷轧钢板裁剪成合适尺寸的坯料,坯料尺寸根据覆盖件的形状和冲压工艺要求确定,确保坯料在冲压过程中有足够的材料流动和变形空间。对坯料进行表面处理,去除表面的油污和氧化皮,以保证板料与模具之间的良好接触和润滑效果。在坯料表面均匀涂抹润滑剂,选用专用的冲压润滑剂,其能够在板料与模具之间形成一层均匀的润滑膜,有效降低摩擦系数,减少摩擦力对冲压过程的影响。将坯料放置在模具上,调整好位置后,启动冲压设备进行冲压成形实验。在冲压过程中,通过传感器实时监测冲压速度、压边力等工艺参数,并将数据记录下来,以便后续分析。冲压完成后,小心取出冲压件,避免对冲压件造成额外的变形。对冲压件进行回弹测量,采用先进的测量技术和设备,如三坐标测量仪、数字图像相关法等,准确测量冲压件的回弹量和回弹分布情况。为了确保实验结果的可靠性和准确性,每个实验条件下均进行多次重复实验,取平均值作为实验结果。在不同的压边力条件下,分别进行5次冲压实验,对每次实验得到的冲压件进行回弹测量,然后计算平均值和标准差,以评估实验结果的稳定性和离散性。通过多次重复实验,可以有效减少实验误差,提高实验数据的可信度,为后续的回弹分析和研究提供坚实的数据基础。4.2.2回弹测量技术在覆盖件冲压成形回弹研究中,准确测量回弹量是关键环节,三坐标测量仪和数字图像相关法是两种常用且有效的测量技术。三坐标测量仪是一种高精度的测量设备,其测量原理基于空间坐标测量。该仪器主要由主机、测头和电气系统三大部分构成。主机部分包括框架结构、标尺系统、导轨、驱动装置、平衡部件及转台与附件,为测量提供稳定的物理支撑。测头作为核心部件,能够在三个方向上感受瞄准信号和微小位移,实现瞄准与测微功能。电气系统则负责控制测量过程、处理测量数据并输出结果。在测量覆盖件回弹量时,首先将冲压件放置在三坐标测量仪的工作台上,通过测头在冲压件表面逐点测量,获取冲压件上各测点的三维坐标数据。测量软件根据预先设定的测量路径和测量策略,自动控制测头的运动,确保测量的准确性和高效性。通过测量软件对获取的坐标数据进行处理和分析,将测量得到的冲压件实际形状与模具型面的理论形状进行对比,计算出各测点的回弹量。根据这些回弹量数据,可以绘制出冲压件的回弹分布云图,直观地展示回弹量在冲压件上的分布情况。三坐标测量仪具有测量精度高、测量范围广、可测量复杂形状等优点,其测量精度可达±0.001mm,能够满足对覆盖件回弹量高精度测量的需求。然而,该方法也存在一些局限性,如测量速度相对较慢,对于大型覆盖件的测量需要较长时间;测量过程中测头与冲压件表面接触,可能会对冲压件表面造成轻微损伤。数字图像相关法是一种基于光学原理的非接触式测量技术,近年来在材料力学性能测试和结构变形测量领域得到了广泛应用。其基本原理是利用数字图像处理技术,对物体表面变形前后的数字图像进行相关分析,从而计算出物体表面各点的位移和应变。在覆盖件回弹测量中,首先在冲压件表面制作随机散斑图案,散斑图案的质量对测量精度有重要影响,要求散斑分布均匀、对比度高。利用高分辨率相机从不同角度拍摄冲压件变形前后的图像,获取两组图像数据。将拍摄的图像导入计算机,通过数字图像相关算法对两组图像进行处理和分析。算法通过寻找变形前后图像中散斑的对应关系,计算出冲压件表面各点在x、y、z方向上的位移,进而得到回弹量。数字图像相关法具有非接触、全场测量、测量速度快、精度较高等优点,能够实时测量冲压件在卸载过程中的回弹变形过程,获取回弹的动态变化信息。其测量精度可达亚像素级别,对于一般的覆盖件回弹测量,能够满足工程应用的要求。该方法也受到一些因素的影响,如光照条件、散斑制作质量、图像噪声等,在实际应用中需要对这些因素进行严格控制,以保证测量结果的准确性。五、覆盖件冲压成形数值模拟与回弹研究的关联5.1数值模拟对回弹预测的作用在覆盖件冲压成形领域,准确预测回弹量对于保证产品质量和提高生产效率至关重要。数值模拟技术为回弹预测提供了一种高效、准确的手段,通过在计算机上模拟冲压成形过程,可以深入分析板料的变形行为和应力应变分布,从而预测回弹量的大小和分布情况。在数值模拟中,通常采用弹塑性有限元法来模拟板料的冲压成形过程。首先,将板料离散化为有限个单元,通过定义材料的本构关系和接触条件,建立板料与模具之间的力学模型。在模拟汽车覆盖件的冲压成形时,使用四节点壳单元对板料进行离散化,选择合适的材料模型,如Hill屈服准则来描述材料的塑性行为。在模拟过程中,考虑板料与模具之间的接触和摩擦,采用罚函数法来处理接触问题,设置合理的摩擦系数来模拟实际的摩擦状态。通过数值模拟,可以得到板料在冲压过程中的应力应变分布情况,这些信息是预测回弹量的关键。在冲压成形结束后,板料内部存在残余应力,当卸载时,这些残余应力会导致板料发生弹性恢复,从而产生回弹。通过模拟得到的应力应变分布,可以计算出板料在卸载过程中的弹性恢复量,进而预测回弹量。利用模拟结果中的等效应力和等效应变分布,结合材料的弹性模量和泊松比,通过弹性力学理论计算出板料各部分的弹性恢复变形,从而得到回弹量。为了验证数值模拟对回弹预测的准确性,进行了对比实验。以某汽车覆盖件为例,首先利用DYNAFORM软件进行数值模拟,预测其回弹量。在模拟过程中,严格按照实际冲压工艺参数进行设置,包括冲压速度、压边力、摩擦系数等。模拟完成后,得到了该覆盖件的回弹量预测结果,如在某些关键部位的回弹量大小和回弹方向。随后,进行实际冲压实验,使用相同的材料和模具,在相同的工艺条件下进行冲压成形,并采用三坐标测量仪对冲压件的实际回弹量进行测量。对比模拟结果与实际测量结果发现,在大多数部位,模拟预测的回弹量与实际测量值较为接近,误差在可接受范围内。在覆盖件的边缘部位,模拟预测的回弹量为2.5mm,实际测量值为2.8mm,误差约为10.7%。这表明数值模拟能够较为准确地预测覆盖件冲压成形的回弹量,为后续的工艺优化和模具设计提供了可靠的依据。然而,在一些局部区域,由于实际冲压过程中存在一些难以精确模拟的因素,如模具表面的微观粗糙度、材料性能的微小差异等,导致模拟结果与实际测量值存在一定偏差。在后续的研究中,可以进一步优化数值模拟模型,考虑更多的实际因素,以提高回弹预测的精度。五、覆盖件冲压成形数值模拟与回弹研究的关联5.2基于数值模拟的回弹控制策略5.2.1工艺参数优化工艺参数的优化是控制覆盖件冲压成形回弹的关键手段之一,通过数值模拟分析,可以深入研究冲压速度、压边力等工艺参数对回弹的影响规律,从而实现工艺参数的优化,有效减小回弹量。冲压速度对回弹有着显著的影响。在数值模拟中,通过改变冲压速度进行多组模拟实验,分析冲压速度与回弹之间的关系。当冲压速度较低时,板料在冲压过程中有足够的时间与模具接触并发生变形,应力应变分布相对均匀,回弹量较小。随着冲压速度的增加,板料在短时间内受到较大的冲击力,材料内部的应力应变分布不均匀,局部应力集中现象加剧,导致回弹量增大。在模拟汽车覆盖件的冲压成形时,将冲压速度从1m/s提高到5m/s,发现回弹量增加了约30%。这是因为高速冲压时,板料与模具之间的摩擦生热来不及散失,使局部温度升高,材料的力学性能发生变化,进一步影响了回弹。因此,在实际生产中,应根据覆盖件的材料特性和形状复杂程度,合理选择冲压速度,以减小回弹。对于塑性较好、形状相对简单的覆盖件,可以适当提高冲压速度以提高生产效率;而对于塑性较差、形状复杂的覆盖件,则应选择较低的冲压速度,以保证成形质量。压边力是控制板料流动和变形的重要工艺参数,对回弹也有重要影响。在数值模拟中,通过调整压边力的大小,观察板料的变形情况和回弹量的变化。适当增大压边力,可以增加板料与压边圈之间的摩擦力,限制板料的径向流动,使板料在凸模的作用下更加均匀地变形,从而减小回弹量。当压边力从100kN增大到150kN时,覆盖件的回弹量明显减小。但压边力过大,会使板料受到过度的挤压,增加破裂的风险。在实际应用中,需要通过数值模拟结合实验验证,找到压边力的最佳取值范围。可以采用正交试验设计或响应面法等优化方法,对压边力以及其他工艺参数进行多因素优化,以实现最小的回弹量和最佳的成形质量。除了冲压速度和压边力外,摩擦系数、模具间隙、拉延筋等工艺参数也会对回弹产生影响。摩擦系数的大小影响板料与模具之间的摩擦力,进而影响板料的流动和应力应变分布。模具间隙的大小直接影响板料在冲压过程中的受力状态和变形方式,拉延筋的设置可以调整板料的流动阻力和应力应变分布。通过数值模拟,可以全面研究这些工艺参数对回弹的影响,并进行综合优化,以实现对覆盖件冲压成形回弹的有效控制。5.2.2模具补偿设计模具补偿设计是减少覆盖件回弹的重要方法之一,其核心思想是根据数值模拟结果,对模具型面进行反向修正,使冲压件在回弹后能够达到设计要求的尺寸精度。在覆盖件冲压成形过程中,由于回弹的存在,冲压件的实际形状与模具型面存在偏差。通过模具补偿设计,可以预先在模具型面上增加或减少一定的量,以补偿回弹变形,从而提高冲压件的尺寸精度。基于数值模拟结果进行模具型面补偿设计的过程通常包括以下步骤:首先,利用数值模拟软件对覆盖件冲压成形过程进行模拟,得到冲压件在回弹后的形状和尺寸数据。在模拟过程中,要确保模拟参数与实际生产条件尽可能一致,以保证模拟结果的准确性。通过模拟得到汽车覆盖件回弹后的形状和各部位的回弹量分布情况。然后,根据回弹模拟结果,计算出模具型面需要补偿的量。这需要运用数学算法和几何处理技术,将回弹量转化为模具型面的修正量。根据回弹量的大小和分布,在模具型面上相应部位增加或减少一定的材料厚度,以实现对回弹的补偿。利用专业的模具设计软件,根据计算得到的补偿量,对模具型面进行修改设计。在修改过程中,要考虑模具的加工工艺性和可行性,确保模具能够顺利制造。对修改后的模具进行再次模拟验证,检查补偿效果。如果补偿后的冲压件仍然存在较大的回弹偏差,则需要进一步调整模具型面的补偿量,直到满足设计要求为止。模具补偿设计在实际生产中已经得到了广泛应用,并取得了良好的效果。在某汽车制造企业的汽车覆盖件生产中,通过对模具型面进行补偿设计,成功将覆盖件的回弹量控制在
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 教案嵌入式应用基础(项目2第4周)定时器
- 2026陕西西安雁塔恒通村镇银行社会招聘笔试备考试题及答案详解
- 2026江苏南京工业职业技术大学招聘高层次人才68人笔试备考题库及答案详解
- 2026年宁波余姚市阳明街道办事处公开招聘编外工作人员4人笔试参考题库及答案详解
- 2026中国平安人寿保险股份有限公司嘉兴中心支公司校园招聘笔试备考题库及答案详解
- 2026浙江宁波报业传媒集团有限公司招聘11人笔试备考题库及答案详解
- 2025年华夏银行(苏州分行)校园招聘笔试考试试题及答案详解
- 2025年宁波东海银行校园招聘笔试考试题库及答案详解
- 2026广东江门市五邑人力资源有限公司招聘4人笔试模拟试题及答案详解
- 2025年招商银行(大连分行)人员招聘笔试考试试题及答案详解
- 医院病房深夜起火:如何安全转移卧床病人?医护团队疏散推演
- XX九年一贯制学校2026年春季学期特殊群体学生关爱帮扶工作台账
- 2026浙江广播电视集团社会招聘3人笔试模拟试题及答案解析
- 2025年广东高考生物试卷及答案
- 企业近三年安全生产(施工)记录
- GB/T 9791-2025金属及其他无机覆盖层锌、镉、铝-锌合金和锌-铝合金的铬酸盐转化膜试验方法
- 地源热泵打井协议书
- 征兵智力测试题目集与详解
- 2025ACG临床指南:成人溃疡性结肠炎(更新版)课件
- 口腔器械清洗消毒培训
- 代扣代缴个税协议书
评论
0/150
提交评论