基于ANSYS-LS - DYNA的船体外板非对压多压头冲压成形数值模拟与分析

基于ANSYS/LS-DYNA的船体外板非对压多压头冲压成形数值模拟与分析一、绪论1.1研究背景与意义在船舶制造领域，船体外板作为船体结构的重要组成部分，其成形质量直接关系到船舶的性能、安全性以及使用寿命。传统的船体外板冲压成形方法，如整体模具冲压和水火弯板等，在面对日益复杂的船体外板形状和高精度的制造要求时，逐渐暴露出诸多局限性。整体模具冲压虽在汽车覆盖件等大批量生产中成本低且易实现自动化，但对于形状各异、小批量甚至单件生产的船体外板而言，高昂的模具制造成本和漫长的制造周期使其难以适用。而水火弯板作为目前国内外船厂加工双曲度船体外板的常用方法，主要依靠有经验的工人手工操作，通过热源局部加热并浇水冷却使板材热胀冷缩产生变形。然而，这种方法不仅劳动强度大、能耗高、作业环境恶劣，而且对工人技术要求极高，生产效率低下，质量稳定性差，严重制约了现代造船业向自动化和数字化方向发展的进程。随着科技的不断进步，非对压多压头冲压成形技术应运而生，为船体外板成形领域带来了新的突破。该技术基于活络方形压头，能够灵活调整基本体，确保与板料的最大接触面积，使基本体离散分布，且每个上压头由下压头支撑。这种独特的设计理念使得该技术在船体外板成形中展现出显著的优势。它有效克服了传统方法的弊端，具有更高的成形效率和质量，能够更好地适应船体外板复杂多变的形状需求。通过精确控制压头的运动和压力分布，可以实现对板材变形的精准控制，从而提高成形精度，减少废品率。同时，该技术易于实现自动化和数字化控制，为提高造船效率和质量，降低生产成本提供了有力支持，在船体外板加工技术创新中发挥着重要作用。在现代船舶制造中，数值模拟技术已成为不可或缺的重要工具，对船体外板冲压成形技术的发展具有深远意义。通过数值模拟，可以在实际生产之前对船体外板冲压成形过程进行全面、深入的研究和分析。借助专业的有限元软件，如ANSYS/LS-DYNA，能够精确模拟冲压过程中板材的应力、应变分布，预测可能出现的缺陷，如起皱、拉裂和回弹等。这使得工程师可以在虚拟环境中对冲压工艺参数进行优化，提前发现并解决潜在问题，避免在实际生产中进行大量的试错，从而大大缩短产品开发周期，降低生产成本。数值模拟还可以为新型冲压设备的研发和改进提供理论依据，推动船体外板冲压成形技术不断向更高水平发展，提升整个船舶制造业的竞争力。1.2国内外研究现状1.2.1船体外板冲压成形方法研究常见的船体外板冲压成形方法多种多样，各有其特点和适用范围。整体模具冲压在汽车制造等大批量生产领域应用广泛，其模具通常根据零件的形状和尺寸进行专门设计制造，具有较高的生产效率和成形精度。通过模具的精确配合，能够使板材在短时间内冲压成所需形状，并且零件的一致性较好，便于实现机械化和自动化生产。然而，在船舶制造中，由于船体外板形状复杂且多为单件或小批量生产，定制专用模具的成本极高，而且模具制造周期长，无法满足船舶建造的时效性要求，因此整体模具冲压在船体外板加工中受到很大限制。水火弯板是目前国内外船厂加工双曲度船体外板的主要方法之一。它利用板材在局部加热和冷却过程中产生的热胀冷缩效应，使板材发生塑性变形，从而获得所需的曲率。具体操作时，先使用辊板机或油压机对板材进行横向曲度的初步弯曲加工，然后在纵向上通过水火弯板技术，沿着预定的水火线进行加热，并在加热处浇水冷却，实现纵向曲度的弯曲成形。在加工过程中，工人需要依靠样板或样箱进行对比和校验，以确保弯板的形状符合要求。尽管水火弯板在船体外板加工中应用广泛，但它存在诸多缺点。该方法严重依赖工人的经验和技能水平，不同工人操作可能导致产品质量存在较大差异；生产效率低下，每一块外板的加工都需要耗费大量时间；能耗大，对环境造成一定污染；并且由于其成形过程复杂，影响因素众多，实现自动化控制的难度极大。虽然日本、韩国等造船强国以及我国都曾投入大量资源研发数控水火弯板机，但至今相关设备仍未得到广泛应用。除了上述两种方法，还有一些其他的冲压成形方法也在不断发展和探索中。多点成形技术将整体模具离散成有限个高度可调的基本体，通过冲头群代替整体模具来完成曲面成形。这种方法具有“柔性”加工特点，能够在成形前和成形过程中实时改变基本体的位置，从而调整被成形件的变形路径及受力状态，以适应不同的成形需求。然而，多点成形技术在实际应用中也面临一些问题，如下模调形驱动元器件的尺寸限制了可重构达到的最小曲率，导致在加工小曲率半径外板时效果不佳，往往需要与其他工艺混合使用，增加了加工的复杂性和成本。1.2.2船体外板冲压成形数值模拟研究数值模拟技术在船体外板冲压成形领域的应用日益广泛，为该领域的研究和发展提供了强大的支持。通过数值模拟，可以深入了解冲压过程中板材的变形机理，预测成形缺陷，优化冲压工艺参数，从而提高船体外板的成形质量和生产效率。在数值模拟中，常用的软件包括ANSYS/LS-DYNA、AutoForm、DYNAFORM等。ANSYS/LS-DYNA是一款功能强大的显式动力学分析软件，在船体外板冲压成形模拟中，它能够精确模拟板材在冲压过程中的大变形、接触非线性以及材料非线性等复杂力学行为。通过建立合理的有限元模型，可以准确计算板材的应力、应变分布，预测起皱、拉裂等缺陷的发生位置和程度。AutoForm则是一款专门用于金属板料成形模拟的软件，它在冲压工艺分析和模具设计方面具有独特的优势，能够快速准确地预测板材的成形性能，优化压边力、模具间隙等工艺参数，减少试模次数，降低模具开发成本。DYNAFORM也广泛应用于冲压成形模拟，它具备良好的前后处理功能和求解精度，能够对复杂的冲压过程进行高效模拟，为冲压工艺的优化提供可靠依据。国内外学者利用这些软件在船体外板冲压成形数值模拟方面取得了丰硕的研究成果。一些研究通过数值模拟分析了不同冲压工艺参数对船体外板成形质量的影响，如压边力、冲压速度、模具间隙等，并建立了相应的数学模型，为工艺参数的优化提供了理论依据。还有学者对不同材料的船体外板冲压成形过程进行了模拟研究，分析了材料性能对成形结果的影响，为材料的选择和应用提供了参考。此外，在数值模拟方法和技术方面也不断有新的进展，如采用自适应网格划分技术提高模拟精度，结合人工智能算法实现工艺参数的自动优化等。然而，当前的研究仍存在一些不足之处。一方面，数值模拟模型的准确性和可靠性有待进一步提高。虽然现有的软件能够模拟大部分冲压成形过程，但在某些复杂情况下，如材料的微观组织变化、模具与板材之间的复杂接触摩擦等，模拟结果与实际情况仍存在一定偏差。另一方面，对于船体外板冲压成形过程中的多物理场耦合问题，如热-力耦合、力-电耦合等，研究还相对较少。在实际冲压过程中，这些多物理场因素往往相互影响，对成形质量产生重要作用，因此需要进一步加强这方面的研究，以完善数值模拟理论和方法，提高模拟结果的准确性和实用性。1.3研究内容与方法本文基于ANSYS/LS-DYNA软件展开船体外板非对压多压头冲压成形数值模拟研究，具体内容如下：建立精确的有限元模型：深入分析船体外板非对压多压头冲压成形的工作原理，综合考虑板材的材料特性、压头的结构和运动方式以及模具与板材之间的接触条件等因素，运用ANSYS软件强大的前处理功能，构建高精度的有限元模型。对板材和压头进行合理的网格划分，确保网格质量满足计算要求，准确模拟冲压过程中的力学行为。模拟冲压成形过程并分析结果：借助LS-DYNA软件卓越的非线性动力学分析功能，对建立的有限元模型进行求解，模拟船体外板在非对压多压头冲压作用下的成形过程。详细分析模拟结果，包括板材的应力、应变分布情况，观察板材在冲压过程中的变形规律，明确应力集中区域和变形较大部位，为后续的工艺优化提供数据支持。研究工艺参数对成形质量的影响：系统研究压头压力、冲压速度、压头间距等关键工艺参数对船体外板成形质量的影响。通过改变这些工艺参数的值，进行多组数值模拟实验，对比不同参数组合下的模拟结果，分析各参数对板材起皱、拉裂和回弹等缺陷的影响规律，从而确定最佳的工艺参数范围，为实际生产提供科学的工艺指导。验证模拟结果的准确性：为了确保数值模拟结果的可靠性，将模拟结果与实际冲压实验数据进行对比分析。设计并进行船体外板非对压多压头冲压成形实验，测量实验过程中板材的变形情况和应力分布，与数值模拟结果进行详细比对。根据对比结果，对有限元模型和模拟参数进行调整和优化，进一步提高模拟结果的准确性和可信度。在研究方法上，主要采用以下几种：理论分析方法：运用金属塑性成形理论、材料力学、弹性力学等相关学科的基本原理，深入分析船体外板非对压多压头冲压成形过程中的力学行为和变形机理。建立相应的数学模型，推导相关公式，为数值模拟和实验研究提供理论基础。数值模拟方法：利用ANSYS/LS-DYNA软件强大的数值计算功能，对船体外板非对压多压头冲压成形过程进行全面、细致的模拟。通过建立精确的有限元模型，设置合理的材料参数、接触条件和边界条件，模拟不同工艺参数下的冲压过程，获取丰富的模拟数据，并对这些数据进行深入分析，揭示工艺参数与成形质量之间的内在关系。实验研究方法：设计并开展船体外板非对压多压头冲压成形实验，验证数值模拟结果的准确性。根据实际生产条件和模拟研究的需求，制作实验模具和试件，搭建实验平台，选择合适的实验设备进行数据测量和采集。通过实验，获取实际冲压过程中的数据，如板材的变形量、应力分布等，并与数值模拟结果进行对比分析，对模拟模型进行修正和完善。对比分析方法：在研究过程中，对不同工艺参数下的模拟结果和实验数据进行对比分析，找出各参数对船体外板成形质量的影响规律。同时，将非对压多压头冲压成形技术与传统冲压成形方法进行对比，评估非对压多压头冲压成形技术在提高成形质量、生产效率和降低成本等方面的优势，为该技术的推广应用提供有力依据。二、相关理论基础2.1船体外板非对压多压头冲压成形原理船体外板非对压多压头冲压成形技术是一种创新的板材成形方法，它基于活络方形压头的独特设计，旨在实现船体外板复杂曲面的精确成形。在该技术中，压头的排列方式和调整机制是实现高效、精确成形的关键要素。压头通常呈方形，以矩阵形式离散分布，构成了可灵活调整的成形模具。这种离散分布的方式使压头能够适应船体外板复杂多样的形状，克服了传统整体模具的局限性。在冲压过程中，每个上压头由对应的下压头支撑，形成了稳定的支撑结构，确保在施加压力时，压头能够均匀地作用于船体外板，避免局部应力集中导致的板材损伤或成形缺陷。下压头通过一套精密的控制系统实现高度和位置的精确调整。该控制系统借助先进的计算机技术和传感器反馈，能够根据预先设定的成形目标和板材实时变形情况，快速、准确地调整下压头的参数。通过测量板材的实时应变和位移，控制系统可以动态地改变下压头的高度，使压头与板材之间的接触力始终保持在最佳状态，从而精确控制板材的变形路径和程度。这种实时调整机制使得非对压多压头冲压成形技术能够在不同的工况下，都能保证船体外板的高质量成形。当船体外板放置在压头阵列上时，压头开始按照预设的程序对板材施加压力。随着压头的逐渐下压，板材在压力的作用下开始发生塑性变形。由于压头与板材的接触面积大且分布均匀，板材的变形更加均匀，有效减少了起皱、拉裂等缺陷的产生。在整个冲压过程中，压头与板材之间始终保持着紧密的相互作用，压头通过不断调整自身的位置和压力，引导板材逐步达到预期的形状，最终实现船体外板的精确成形。这种非对压多压头冲压成形原理，不仅提高了船体外板的成形质量和效率，还为船舶制造的自动化和数字化发展提供了有力支持。2.2ANSYS/LS-DYNA软件介绍ANSYS/LS-DYNA是一款功能强大的显式动力学分析软件，在众多工程领域中发挥着重要作用。它由ANSYS公司开发，融合了ANSYS软件在结构分析、热分析、电磁场分析等多物理场领域的卓越能力，以及LS-DYNA在显式动力学计算方面的专业优势，为用户提供了一个全面、高效的仿真平台。该软件具备丰富的功能特点，能够处理多种复杂的物理现象和工程问题。它拥有强大的非线性分析能力，能够精确模拟材料在大变形、高速冲击、爆炸等极端条件下的非线性力学行为。在模拟金属板材冲压成形过程中，它可以准确考虑材料的塑性变形、应变硬化、各向异性等特性，以及模具与板材之间的接触、摩擦等非线性因素，从而得到高精度的模拟结果。ANSYS/LS-DYNA还支持多物理场耦合分析，能够处理热-力耦合、流-固耦合等复杂问题，为研究涉及多个物理过程相互作用的工程现象提供了有力工具。ANSYS/LS-DYNA的适用范围广泛，涵盖了机械工程、航空航天、汽车制造、船舶工业、国防军工等多个领域。在机械工程领域，它可用于模拟各种机械零件的加工过程，如锻造、冲压、切削等，优化加工工艺参数，提高零件的加工质量和生产效率。在航空航天领域，该软件能够模拟飞行器在飞行过程中的气动力、热应力、结构响应等，为飞行器的设计和优化提供重要依据。在汽车制造领域，ANSYS/LS-DYNA常用于汽车碰撞模拟、安全气囊展开分析、车身结构优化等，以提高汽车的安全性和性能。在船舶工业中，它可用于船体外板冲压成形模拟、船舶碰撞分析、波浪载荷作用下的船体结构响应分析等，为船舶的设计、制造和性能评估提供技术支持。在金属成形数值模拟领域，ANSYS/LS-DYNA具有显著的优势。它拥有丰富的材料模型库，包含了各种常见金属材料以及新型材料的本构模型，用户可以根据实际需求快速选择合适的材料模型，并准确设置材料参数，从而提高模拟结果的准确性。该软件提供了多种先进的接触算法和摩擦模型，能够精确模拟模具与金属板材之间的接触状态和摩擦行为，有效解决了金属成形过程中接触非线性的难题。ANSYS/LS-DYNA还具备高效的求解器和强大的并行计算能力，能够在较短的时间内完成大规模、复杂模型的计算，大大提高了数值模拟的效率，满足了工程实际中对计算速度的要求。2.3数值模拟关键技术2.3.1有限元模型建立在基于ANSYS/LS-DYNA进行船体外板非对压多压头冲压成形数值模拟时，建立合理的有限元模型是获得准确模拟结果的基础。首先，需根据船体外板冲压成形的实际情况，对几何模型进行简化。由于船体外板形状复杂，包含众多细节特征，若直接对其进行建模，不仅会增加模型的复杂性和计算量，还可能导致计算效率低下甚至无法求解。因此，需要在不影响模拟结果准确性的前提下，对一些次要的几何特征进行适当简化，如去除微小的圆角、倒角、孔洞等。在简化过程中，需充分考虑这些特征对板材整体力学性能和变形行为的影响，确保简化后的模型能够真实反映实际冲压过程。完成几何模型简化后，接下来进行网格划分。网格划分的质量直接影响模拟结果的精度和计算效率。对于船体外板和压头，通常采用四边形或六面体单元进行网格划分，因为这些单元在模拟复杂形状和大变形时具有较好的适应性和计算精度。在划分网格时，需要根据模型的几何形状、应力应变分布以及计算精度要求，合理控制网格尺寸。在应力集中区域和变形较大的部位，如船体外板的边缘、与压头接触的区域等，应采用较小的网格尺寸，以更精确地捕捉应力应变的变化；而在应力应变变化较小的区域，可以适当增大网格尺寸，以减少计算量。还需注意网格的质量，确保网格单元的形状规则、尺寸均匀，避免出现畸形单元，以免影响计算结果的准确性。可以通过网格质量检查工具，对划分好的网格进行质量评估，如检查单元的长宽比、雅克比行列式等指标，对于质量不合格的网格进行局部调整或重新划分。2.3.2材料模型选择船体外板材料的特性对冲压成形过程和成形质量有着重要影响。常见的船体外板材料主要有高强度合金钢、铝合金等。高强度合金钢具有较高的强度和韧性，能够承受较大的外力作用，但其塑性变形能力相对较弱，在冲压过程中容易出现裂纹等缺陷。铝合金则具有密度小、比强度高、塑性好等优点，但其弹性模量较低，成形后容易产生回弹现象。在数值模拟中，需要根据船体外板材料的具体特性选择合适的材料模型。对于金属材料，常用的材料模型有弹性-塑性模型、弹塑性硬化模型等。弹性-塑性模型适用于描述材料在弹性阶段和塑性阶段的基本力学行为，它假设材料在弹性阶段服从胡克定律，在塑性阶段满足屈服准则，如Mises屈服准则或Tresca屈服准则。弹塑性硬化模型则进一步考虑了材料在塑性变形过程中的硬化效应，即随着塑性变形的增加，材料的屈服强度逐渐提高。这种模型更能准确地描述材料在实际冲压过程中的力学行为，对于预测板材的变形和缺陷具有重要意义。确定材料模型后，还需准确确定模型参数。这些参数通常包括弹性模量、泊松比、屈服强度、硬化参数等，它们可以通过材料试验或查阅相关材料手册获得。弹性模量和泊松比反映了材料的弹性性质，屈服强度则是材料进入塑性变形的临界值，硬化参数决定了材料硬化的程度和速率。在实际模拟中，需要根据具体的材料和冲压工艺条件，对这些参数进行合理调整，以确保材料模型能够准确反映材料的真实力学性能。2.3.3接触与摩擦处理在船体外板冲压成形过程中，船体外板与压头之间存在着复杂的接触关系，这种接触关系对板材的变形和成形质量有着关键影响。在数值模拟中，准确处理接触和摩擦问题是提高模拟结果准确性的重要环节。船体外板与压头之间的接触属于面-面接触，在ANSYS/LS-DYNA中，可以采用多种接触算法来模拟这种接触行为，如罚函数法、拉格朗日乘子法等。罚函数法是一种常用的接触算法，它通过在接触面上引入一个罚因子，将接触约束转化为接触力，从而实现对接触问题的求解。这种方法计算效率较高，但罚因子的选择对计算结果的稳定性和准确性有较大影响，若罚因子过大，可能导致计算结果不稳定；若罚因子过小，则可能无法准确模拟接触行为。拉格朗日乘子法则是通过引入拉格朗日乘子来满足接触约束条件，该方法能够更精确地模拟接触行为，但计算量相对较大。在实际模拟中，需要根据具体问题的特点和计算要求，选择合适的接触算法。摩擦是影响船体外板冲压成形的另一个重要因素。船体外板与压头之间的摩擦会改变板材的受力状态和变形路径，进而影响成形质量。在数值模拟中，常用的摩擦模型有库仑摩擦模型和修正的库仑摩擦模型等。库仑摩擦模型假设摩擦力与接触面上的正压力成正比，其比例系数为摩擦系数。这种模型简单直观，计算效率较高，但它没有考虑摩擦系数随相对滑动速度、接触压力等因素的变化。修正的库仑摩擦模型则对库仑摩擦模型进行了改进，考虑了更多的影响因素，能够更准确地描述实际的摩擦行为。在确定摩擦系数时，可以通过查阅相关文献、进行摩擦试验或参考经验数据来获取。由于摩擦系数在实际冲压过程中可能会受到多种因素的影响而发生变化，因此在模拟中可以采用变摩擦系数的方法，以更真实地反映摩擦行为。2.3.4加载与边界条件设置加载方式和边界条件的设置是数值模拟中至关重要的环节，它们直接关系到模拟结果能否真实反映船体外板的实际冲压成形过程。在实际冲压过程中，压头对船体外板施加压力使其发生塑性变形。在数值模拟中，需要根据实际冲压情况设置合适的加载方式。常见的加载方式有位移加载和力加载两种。位移加载是通过指定压头的位移随时间的变化规律来施加加载，这种加载方式能够准确控制压头的运动过程，适用于对冲压过程中板材变形路径有精确要求的情况。力加载则是直接指定压头对板材施加的力的大小和方向，它更能反映实际冲压过程中力的作用情况，适用于研究板材在不同力作用下的力学响应。在选择加载方式时，需要考虑实际冲压工艺、模拟目的以及计算效率等因素。若模拟目的是研究板材的成形极限，则可以采用位移加载方式，以便更准确地控制板材的变形程度；若关注的是板材在不同力作用下的应力应变分布，则力加载方式更为合适。边界条件的设置同样重要，它主要包括对船体外板和压头的约束条件。对于船体外板，通常在其边缘设置固定约束，以模拟实际冲压过程中板材被夹紧的状态。这种固定约束可以限制板材在某些方向上的位移和转动，确保板材在冲压过程中的稳定性。对于压头，需要根据其实际运动情况设置相应的边界条件。下压头在调整高度时，需要在垂直方向上设置可移动的边界条件，以模拟其能够自由升降的功能；而上压头在冲压过程中，通常在水平方向上设置固定约束，以保证其在施加压力时不会发生横向移动，从而准确地对船体外板施加压力。合理设置加载与边界条件，能够使数值模拟更加贴近实际冲压过程，为准确分析船体外板的冲压成形行为提供保障。三、数值模拟方案设计3.1模拟工况设定为全面研究船体外板非对压多压头冲压成形过程，设定多种模拟工况，涵盖不同船体外板形状、板材厚度以及曲率半径，具体参数如下：模拟工况具体参数船体外板形状球形：用于模拟船舶部分曲面较为规则的区域，如船艏或船艉的某些部位，其形状特点是在各个方向上的曲率较为均匀；鞍形：适用于模拟船舶中部等具有复杂双曲率的区域，该形状在不同方向上的曲率变化明显，更能体现非对压多压头冲压成形技术在处理复杂形状时的性能板材厚度5mm：较薄的板材，在冲压过程中可能更容易发生变形，但也更容易出现起皱等缺陷，研究该厚度板材的冲压成形情况，有助于了解薄板在该技术下的成形特性；10mm：中等厚度的板材，综合考虑了材料成本和结构强度要求，分析其冲压过程，对于实际生产中该厚度板材的应用具有指导意义；20mm：相对较厚的板材，在冲压时需要更大的压力，研究其成形过程可以为厚板冲压工艺提供参考；24mm：进一步增加板材厚度，探索在极限厚度情况下非对压多压头冲压成形技术的可行性和存在的问题曲率半径1000mm：较小的曲率半径，代表了船体外板曲率变化较大的区域，研究该工况下的冲压成形，对于处理船体外板的局部复杂形状具有重要意义；1500mm：适中的曲率半径，在船体外板中较为常见，分析该工况有助于掌握一般曲率情况下的冲压工艺；1800mm：较大的曲率半径，模拟船体外板曲率相对平缓的区域，了解在这种情况下非对压多压头冲压成形技术的效果3.2模型参数设置针对不同模拟工况，在ANSYS/LS-DYNA中进行如下模型参数设置：单元类型：对于船体外板和压头，选用适合大变形分析的壳单元，如基于Mindlin理论的壳单元。这种单元能够较好地模拟板材在冲压过程中的弯曲、拉伸等复杂变形行为，同时在计算效率和计算精度之间取得较好的平衡。在划分网格时，根据模拟工况中船体外板形状的复杂程度和板材厚度，合理调整网格尺寸。对于形状复杂、曲率变化大的区域，如鞍形板的边缘和曲率突变处，以及较薄的板材（如5mm厚板材），采用较小的网格尺寸，一般设置为5mm，以提高计算精度，准确捕捉应力应变的变化；对于形状相对简单、曲率变化小的区域，如球形板的大部分区域以及较厚的板材（如24mm厚板材），网格尺寸可适当增大至10mm，以减少计算量，提高计算效率。同时，通过网格质量检查工具，确保网格单元的长宽比、雅克比行列式等指标满足计算要求，避免出现畸形单元，影响计算结果的准确性。材料参数：假设船体外板材料为高强度合金钢，其主要材料参数如下：弹性模量为200GPa，反映了材料在弹性阶段抵抗变形的能力；泊松比为0.3，用于描述材料在受力时横向应变与纵向应变的比值；屈服强度为350MPa，是材料开始发生塑性变形的临界应力值；硬化参数根据材料的硬化特性确定，如采用等向硬化模型时，硬化模量设置为500MPa，用于描述材料在塑性变形过程中屈服强度的提高程度。这些参数通过查阅相关材料手册或进行材料试验获得，在模拟过程中，根据不同的模拟工况，可对材料参数进行适当调整，以更准确地反映材料在实际冲压过程中的力学行为。接触参数：船体外板与压头之间的接触采用面-面接触算法，选择罚函数法进行接触求解。罚因子的取值经过多次调试确定为0.01，该值既能保证计算结果的稳定性，又能较准确地模拟接触行为。摩擦模型选用修正的库仑摩擦模型，考虑了摩擦系数随相对滑动速度和接触压力的变化。通过参考相关文献和试验数据，确定初始摩擦系数为0.15，并在模拟过程中根据接触状态的变化动态调整摩擦系数。在压头与板材开始接触时，由于相对滑动速度较小，摩擦系数取较小值；随着冲压过程的进行，相对滑动速度增大，摩擦系数相应增大，以更真实地反映实际冲压过程中的摩擦行为。3.3模拟流程基于ANSYS/LS-DYNA的船体外板非对压多压头冲压成形数值模拟流程如下：模型建立：利用ANSYS软件的前处理模块，根据设定的模拟工况，创建船体外板和压头的几何模型。对船体外板的几何模型进行必要的简化，去除对冲压成形影响较小的细节特征，如微小的圆角、倒角等，以减少模型的复杂性和计算量。按照前文所述的单元类型和网格划分原则，对船体外板和压头进行网格划分，确保网格质量满足计算要求。参数设置：在ANSYS/LS-DYNA中，根据模拟工况和材料特性，设置模型参数，包括单元类型、材料参数、接触参数等。详细定义材料模型，输入弹性模量、泊松比、屈服强度、硬化参数等材料参数；设置接触算法和摩擦模型，确定罚因子和摩擦系数等接触参数；根据实际冲压情况，设置加载方式为位移加载，指定压头的位移随时间的变化规律，如在0.01s内压头位移达到设计值；同时，设置边界条件，对船体外板边缘施加固定约束，限制其在水平和垂直方向的位移，对压头设置相应的约束和运动条件，确保模拟过程符合实际冲压情况。求解计算：完成模型建立和参数设置后，提交模拟任务进行求解计算。在计算过程中，密切关注计算状态和计算结果，如计算是否收敛、能量是否守恒等。若计算出现不收敛或异常情况，及时检查模型和参数设置，调整相关参数后重新计算。结果分析：计算完成后，利用ANSYS软件的后处理模块对模拟结果进行分析。观察船体外板在冲压过程中的应力、应变分布云图，确定应力集中区域和变形较大的部位；分析板材的起皱、拉裂和回弹等缺陷情况，通过提取相关数据，如起皱高度、拉裂位置和回弹量等，评估冲压成形质量。根据模拟结果，总结不同模拟工况下船体外板的冲压成形规律，为工艺优化提供依据。四、数值模拟结果与分析4.1不同工况下的成形结果通过ANSYS/LS-DYNA软件对多种模拟工况下的船体外板非对压多压头冲压成形过程进行数值模拟，得到了丰富的模拟结果。图1展示了在板材厚度为10mm、曲率半径为1500mm的球形板和鞍形板冲压成形后的变形形状。从图中可以清晰地看出，球形板在冲压后整体呈较为规则的球形，表面光滑，各部分变形较为均匀；而鞍形板由于其复杂的双曲率形状，冲压后的变形更加复杂，在鞍形的中心区域和边缘区域变形程度存在明显差异，中心区域变形相对较大，边缘区域则相对较小。模拟工况板材厚度曲率半径变形形状描述球形板10mm1500mm整体呈规则球形，表面光滑，变形均匀鞍形板10mm1500mm中心区域变形大，边缘区域变形小，变形复杂[此处插入球形板和鞍形板冲压成形后的变形形状云图，图1：不同形状船体外板冲压成形后的变形形状（a）球形板；（b）鞍形板]图2和图3分别给出了球形板和鞍形板在冲压成形后的应力分布云图和应变分布云图。在球形板的应力分布云图中，应力主要集中在与压头接触的区域以及板的边缘部分，这些区域的应力值相对较高，而板的中心部分应力相对较低。在应变分布云图中，应变较大的区域同样集中在与压头接触的部位和边缘，表明这些区域的变形较为剧烈。对于鞍形板，应力分布呈现出明显的不均匀性，在鞍形的凹陷和凸起部位应力集中现象较为突出，且不同区域的应力方向也有所不同。应变分布也与应力分布相对应，在应力集中区域应变值较大，尤其是在鞍形的边缘和曲率变化较大的区域，应变明显高于其他区域。通过对不同工况下的应力应变分布云图的分析，可以深入了解船体外板在冲压成形过程中的受力和变形情况，为后续研究起皱和回弹等缺陷提供重要依据。[此处插入球形板冲压成形后的应力分布云图，图2：球形板冲压成形后的应力分布云图][此处插入球形板冲压成形后的应变分布云图，图3：球形板冲压成形后的应变分布云图][此处插入鞍形板冲压成形后的应力分布云图，图4：鞍形板冲压成形后的应力分布云图][此处插入鞍形板冲压成形后的应变分布云图，图5：鞍形板冲压成形后的应变分布云图]4.2起皱现象分析4.2.1起皱的判定与特征在数值模拟结果中，通过观察板材的变形形态和厚度方向的应变分布来判定船体外板是否发生起皱。当板材表面出现周期性的波纹状起伏，且在厚度方向上存在较大的压应变时，可判定板材发生了起皱现象。起皱的特征主要表现为板材表面的褶皱，褶皱的形状通常为波浪形，其波长和波高会因不同的工况而有所变化。在鞍形板的模拟结果中，起皱现象较为明显，尤其是在板材的中心区域和边缘部分。在中心区域，由于鞍形的特殊形状，板材在冲压过程中受到的压应力分布不均匀，导致中心部分的材料更容易发生失稳起皱，形成的褶皱较为密集，波高相对较大；在边缘部分，由于板材与压头的接触情况较为复杂，摩擦力和应力集中等因素的影响，也容易引发起皱，此处的褶皱相对较为稀疏，但波幅可能较大。4.2.2影响起皱的因素分析为深入分析不同因素对船体外板起皱的影响规律，通过改变板材厚度、曲率半径、压头压力等参数进行多组模拟实验，并对模拟结果进行数据对比和分析。表1展示了不同板材厚度和曲率半径下船体外板的起皱情况。从表中数据可以看出，板材厚度对起皱有显著影响。当曲率半径固定为1500mm时，随着板材厚度从5mm增加到24mm，起皱高度逐渐减小。这是因为较厚的板材具有更高的抗弯刚度，在冲压过程中更能抵抗因压应力导致的失稳起皱，从而使起皱现象得到有效抑制。曲率半径对起皱也有重要影响。当板材厚度固定为10mm时，随着曲率半径从1000mm增大到1800mm，起皱高度呈现逐渐增大的趋势。这是因为曲率半径越大，板材在冲压过程中所受的弯曲程度相对较小，材料更容易在压应力作用下发生失稳起皱，导致起皱现象加剧。板材厚度（mm）曲率半径（mm）起皱高度（mm）515003.51015002.12015001.32415000.81010001.61015002.11018002.7压头压力对起皱的影响也不容忽视。通过改变压头压力进行模拟实验发现，当压头压力较小时，板材变形不充分，起皱现象相对较轻，但可能无法达到预期的成形形状；随着压头压力逐渐增大，板材变形更加充分，但当压力超过一定阈值时，起皱高度会迅速增加。这是因为过大的压头压力会使板材所受的压应力急剧增大，超过材料的临界失稳应力，从而导致起皱加剧。在实际冲压成形过程中，需要合理控制压头压力，在保证板材充分变形的同时，有效抑制起皱现象的发生。4.3回弹现象分析4.3.1回弹的计算与测量在数值模拟中，通过计算冲压成形后板材的最终形状与卸载瞬间形状之间的差异来确定船体外板的回弹量。具体计算方法是，在ANSYS/LS-DYNA软件的后处理模块中，提取板材上特定节点在冲压结束时刻和卸载完成时刻的坐标值，通过坐标差值计算出这些节点在各个方向上的位移变化，从而得到板材的回弹位移。以板材的某一特征点为例，设该点在冲压结束时刻的坐标为(x_1,y_1,z_1)，卸载完成时刻的坐标为(x_2,y_2,z_2)，则该点在x方向上的回弹位移\Deltax=x_2-x_1，同理可计算出y方向和z方向上的回弹位移\Deltay和\Deltaz。通过对多个特征点的回弹位移进行统计和分析，可以得到板材整体的回弹情况。在测量指标方面，除了回弹位移外，还可以计算回弹角度等指标来全面描述板材的回弹现象。对于弯曲成形的船体外板，回弹角度是一个重要的测量指标，它反映了板材在弯曲方向上的回弹程度。通过计算板材弯曲部位在冲压前后的角度变化，即可得到回弹角度。4.3.2影响回弹的因素分析材料性能和冲压工艺参数对船体外板回弹有着重要影响。材料的弹性模量和屈服强度是影响回弹的关键性能参数。弹性模量反映了材料抵抗弹性变形的能力，弹性模量越大，材料在卸载后的回弹越小；屈服强度则决定了材料开始发生塑性变形的难易程度，屈服强度越高，材料在冲压过程中的塑性变形量相对较小，卸载后的回弹也会相应减小。在模拟过程中，通过改变材料的弹性模量和屈服强度进行多组实验，结果表明，当弹性模量从200GPa增加到220GPa时，在相同冲压工艺条件下，船体外板的回弹位移明显减小；当屈服强度从350MPa提高到380MPa时，回弹现象也得到了一定程度的抑制。冲压工艺参数如冲压速度、压边力等对回弹也有显著影响。冲压速度会影响板材在冲压过程中的变形速率和应力分布。当冲压速度过快时，板材内部的应力来不及均匀分布，可能导致局部变形不均匀，从而增加回弹量。通过模拟不同冲压速度下的船体外板冲压成形过程发现，随着冲压速度从0.05m/s增加到0.1m/s，回弹位移有所增大。压边力主要作用于板材的边缘，控制板材在冲压过程中的流动。合适的压边力可以使板材在冲压过程中保持稳定，避免出现起皱和拉裂等缺陷，同时也能对回弹产生影响。当压边力较小时，板材边缘的约束不足，在冲压过程中容易产生较大的塑性变形，卸载后回弹量较大；随着压边力逐渐增大，板材边缘的约束增强，塑性变形得到一定控制，回弹量会相应减小。但压边力过大时，可能会导致板材边缘过度约束，增加拉裂的风险，同时也会对回弹产生不利影响。因此，在实际冲压过程中，需要根据板材的材料性能和具体的冲压工艺要求，合理调整冲压速度和压边力等工艺参数，以有效控制船体外板的回弹，提高成形精度。五、实验验证与对比5.1实验方案设计为验证数值模拟结果的准确性，设计并开展船体外板非对压多压头冲压成形实验。实验设备选用自主研制的非对压多压头冲压成形试验机，该试验机具备精确控制压头运动和压力施加的功能，能够满足实验对不同工艺参数的要求。试验机的压头采用活络方形压头设计，可根据船体外板的形状进行灵活调整，确保与板材的良好接触和均匀施压。实验材料选用与数值模拟相同的高强度合金钢板材，其尺寸为500mm×500mm。在实验前，对板材进行预处理，去除表面的油污、锈迹等杂质，以保证实验结果的准确性。根据数值模拟的工况设置，选择板材厚度为10mm、曲率半径为1500mm的球形板和鞍形板作为实验对象。实验步骤如下：模具安装与调试：将加工好的球形和鞍形模具安装在冲压试验机上，确保模具的安装精度和稳定性。通过试验机的控制系统，对压头的初始位置和运动参数进行调试，使其符合实验要求。板材定位与固定：将预处理后的板材放置在模具上，通过定位装置确保板材的位置准确。使用夹具将板材的边缘固定，以模拟实际冲压过程中板材的约束条件。冲压实验：启动冲压试验机，按照预设的工艺参数进行冲压实验。在冲压过程中，实时监测压头的压力、位移等参数，并记录下来。数据采集与测量：冲压完成后，使用三维激光扫描仪对成形后的船体外板进行扫描，获取其三维形状数据。通过测量软件对扫描数据进行处理，得到板材的实际成形形状、起皱高度、回弹量等参数。同时，使用应变片测量板材在冲压过程中的应变分布，以验证数值模拟中应变计算的准确性。5.2实验结果与数值模拟对比将实验得到的船体外板冲压成形结果与数值模拟结果进行详细对比，对比内容包括成形形状、起皱情况、回弹量等。图6为球形板和鞍形板的实验成形形状与数值模拟成形形状的对比图。从图中可以看出，实验得到的球形板和鞍形板的成形形状与数值模拟结果基本一致，能够较好地达到预期的设计形状。在球形板的实验中，板材表面光滑，整体形状与模拟结果的球形度误差在允许范围内；鞍形板的实验成形形状也准确地反映了鞍形的特征，双曲率变化与模拟结果相符，表明数值模拟能够较为准确地预测船体外板的冲压成形形状。[此处插入球形板和鞍形板的实验成形形状与数值模拟成形形状对比图，图6：球形板和鞍形板的实验与模拟成形形状对比（a）球形板实验；（b）球形板模拟；（c）鞍形板实验；（d）鞍形板模拟]在起皱情况方面，实验观察到的起皱位置和形态与数值模拟结果具有较高的一致性。对于鞍形板，在实验中，板材的中心区域和边缘部分出现了明显的起皱现象，这与数值模拟中该区域的应力应变分布所预测的起皱情况相符。通过测量实验中起皱的高度，并与数值模拟结果进行对比，发现两者之间存在一定的差异。对于厚度为10mm、曲率半径为1500mm的鞍形板，实验测得的起皱高度平均值为2.3mm，而数值模拟结果为2.1mm，误差约为9.5%。这种差异可能是由于实验过程中材料性能的不均匀性、模具与板材之间的实际接触情况与模拟假设存在一定偏差等因素导致的。在回弹量方面，实验测量得到的球形板和鞍形板的回弹量与数值模拟结果对比如表2所示。从表中数据可以看出，实验测量的回弹量与数值模拟结果在趋势上是一致的，但具体数值存在一定差异。对于球形板，实验测量的回弹位移在x方向上为0.8mm，y方向上为0.7mm，z方向上为0.9mm，而数值模拟结果在x方向上为0.7mm，y方向上为0.6mm，z方向上为0.8mm，各方向的误差分别为12.5%、14.3%、11.1%。对于鞍形板，实验测量的回弹位移在x方向上为1.1mm，y方向上为1.3mm，z方向上为1.2mm，数值模拟结果在x方向上为0.9mm，y方向上为1.1mm，z方向上为1.0mm，各方向的误差分别为18.2%、15.4%、16.7%。这些误差可能是由于实验过程中材料的加工硬化、模具的弹性变形以及测量误差等因素引起的。板形方向实验回弹位移（mm）模拟回弹位移（mm）误差（%）球形板x0.80.712.5球形板y0.70.614.3球形板z0.90.811.1鞍形板x1.10.918.2鞍形板y1.31.115.4鞍形板z1.21.016.75.3误差分析实验结果与数值模拟结果之间存在误差的原因主要包括以下几个方面：材料性能差异：在数值模拟中，材料参数通常是基于理想状态下的测试数据确定的，但实际实验中使用的材料可能存在性能上的不均匀性和波动。材料的微观组织结构、化学成分的微小差异等都可能导致材料的弹性模量、屈服强度等性能参数与模拟设定值不完全一致，从而影响冲压成形过程中的应力应变分布和回弹情况，导致实验结果与模拟结果出现偏差。模具与板材接触状态：虽然在数值模拟中对模具与板材之间的接触和摩擦进行了考虑，但实际冲压过程中，模具与板材的接触状态非常复杂，难以完全准确地模拟。模具表面的粗糙度、板材与模具之间的润滑条件等因素都会影响接触摩擦力的大小和分布，进而影响板材的变形行为。在实验中，由于模具表面的加工精度和润滑情况与模拟假设存在差异，可能导致实际接触摩擦力与模拟值不同，从而使实验结果与数值模拟结果产生误差。测量误差：在实验过程中，对船体外板的成形形状、起皱高度和回弹量等参数的测量不可避免地存在一定的误差。三维激光扫描仪的测量精度、测量过程中的操作误差以及测量数据的处理方法等都可能影响测量结果的准确性。在使用三维激光扫描仪测量板材的回弹量时，由于扫描角度、扫描范围以及点云数据处理算法的不同，可能会导致测量结果存在一定的偏差。这些误差对研究结论的影响总体上是可控的。虽然实验结果与数值模拟结果存在一定差异，但在成形形状、起皱和回弹的趋势等方面基本一致，说明数值模拟能够较好地反映船体外板非对压多压头冲压成形的基本规律。为减小误差，可以采取以下方法和建议：优化材料参数测量：在实验前，对使用的材料进行更加全面和精确的性能测试，获取材料性能参数的实际分布范围。在数值模拟中，考虑材料性能的不确定性，采用随机变量或区间变量来描述材料参数，通过多次模拟计算，分析材料性能波动对冲压成形结果的影响，从而提高模拟结果的可靠性。改进接触摩擦模型：开展更多关于模具与板材之间接触摩擦的实验研究，深入了解实际冲压过程中接触摩擦力的变化规律。根据实验结果，对数值模拟中的接触摩擦模型进行改进和优化，考虑更多影响接触摩擦的因素，如模具表面粗糙度、润滑状态随冲压过程的变化等，