基于时域精细算法与扩展有限元技术的粘弹性位移场分析：理论、方法与应用

基于时域精细算法与扩展有限元技术的粘弹性位移场分析：理论、方法与应用一、引言1.1研究背景与意义在现代工程领域，粘弹性材料的应用极为广泛。从航空航天中飞行器的机翼结构、机身蒙皮，到汽车制造里的轮胎、减震器、密封件，再到土木建筑中的桥梁支座、建筑结构的阻尼器，以及生物医学中用于模拟生物组织力学行为的材料等，粘弹性材料都发挥着关键作用。例如在航空航天领域，飞行器在飞行过程中，机翼受到气流的复杂作用，粘弹性材料制成的机翼结构部件能够更好地适应这种动态载荷，减少疲劳损伤，提高飞行安全性；在汽车行业，轮胎需要具备良好的耐磨性和抗冲击性，粘弹性材料的应用可以有效降低轮胎与路面之间的摩擦系数，减少能量损失和热量产生，同时提高轮胎的抗冲击能力，延长使用寿命。然而，传统的分析方法在处理粘弹性材料的位移场时存在诸多局限性。有限元法（FEM）作为常用的数值方法，当材料的应力应变关系不明确时，其计算过程会面临很大困难。因为有限元法基于离散化的思想，将连续体划分为有限个单元，通过单元的力学特性来求解整体结构的力学响应。对于粘弹性材料，其应力应变关系不仅与当前的应变状态有关，还与应变的历史过程相关，这种复杂的关系使得有限元法在建立单元刚度矩阵和求解方程时变得复杂，计算精度和效率受到影响。边界元法（BEM）虽然在处理边界问题上有一定优势，但它只能处理绝对边界，对于任意边界的分析能力有限，无法满足复杂工程结构中粘弹性材料位移场分析的需求。基于时域精细算法与扩展有限元技术对粘弹性位移场进行研究具有重要意义。时域精细算法通过对时间域进行精细划分，能够更准确地描述粘弹性材料力学行为随时间的变化，提高计算精度。扩展有限元技术则在传统有限元的位移模式中增加了能反映间断问题的改进函数项，并附加节点自由度，采用水平集法描述追踪界面动态变化，使其可以处理裂纹、孔洞、夹杂等间断问题，且计算网格与结构内部无关，无需多次进行网格重划分，大大降低了计算成本。这种结合的研究方法，为粘弹性材料在复杂工程结构中的应用提供了更精确的位移场分析手段，有助于优化工程结构设计，提高结构的安全性和可靠性，在航空航天、汽车、土木等众多领域具有广阔的应用前景和实际价值。1.2国内外研究现状在粘弹性位移场分析领域，时域精细算法和扩展有限元技术都取得了一定的研究进展。时域精细算法的研究不断深入，许多学者致力于算法的改进与应用拓展。一些研究通过优化算法的时间步长选取策略，提高了计算效率和稳定性。例如，有学者提出自适应时间步长控制方法，根据问题的复杂程度和计算精度要求动态调整时间步长，在保证精度的前提下减少了计算量。在应用方面，时域精细算法已被应用于航空航天结构的动力学分析，如对飞行器机翼在复杂气动力作用下的粘弹性响应进行计算，能够准确捕捉结构的瞬态位移和应力变化。在机械工程中，对于高速旋转机械的粘弹性部件，时域精细算法也能有效分析其在动态载荷下的力学行为，为部件的设计和优化提供理论依据。扩展有限元技术在处理间断问题上展现出独特优势，其研究主要集中在理论完善和应用范围扩大。在理论方面，学者们对改进函数项的构造进行了深入研究，提出了多种更能准确反映裂纹、孔洞等间断特征的函数形式，提高了扩展有限元对复杂间断问题的模拟能力。在应用上，扩展有限元技术在土木工程领域得到广泛应用，如用于分析混凝土结构中的裂缝扩展，能够直观地展示裂缝的发展过程和对结构位移场的影响。在生物医学工程中，对于模拟生物组织中的损伤和缺陷，扩展有限元技术也为研究生物组织的力学性能提供了有力工具。然而，当前研究仍存在一些不足。在时域精细算法与扩展有限元技术的结合方面，虽然已有相关研究，但二者的融合还不够完善，在处理复杂粘弹性问题时，计算效率和精度仍有待提高。在粘弹性本构模型的选择和参数确定上，还存在一定的主观性和不确定性，不同的本构模型和参数取值会对位移场分析结果产生较大影响。对于多物理场耦合作用下的粘弹性位移场分析，如热-力、流-固耦合等，现有的研究还相对较少，无法满足实际工程中复杂工况的需求。1.3研究内容与方法1.3.1研究内容基于时域精细算法构建粘弹性位移场分析模型：深入剖析粘弹性材料的力学特性，如应力松弛、蠕变等随时间变化的特性。以这些特性为基础，结合时域精细算法的原理，建立能够准确描述粘弹性材料位移场随时间演变的数学模型。在建立模型过程中，详细推导模型的控制方程，明确方程中各项参数的物理意义和取值方法，确保模型的准确性和可靠性。例如，对于粘弹性材料在动态载荷作用下的位移响应，通过模型能够精确计算不同时刻的位移值，为后续分析提供理论基础。利用扩展有限元技术改进模型以适应复杂边界条件：在已建立的基于时域精细算法的位移场分析模型基础上，引入扩展有限元技术。针对复杂工程结构中粘弹性材料可能存在的裂纹、孔洞、夹杂等间断问题，以及不规则的边界形状，对模型进行改进。通过在传统有限元的位移模式中增加能反映间断问题的改进函数项，并附加节点自由度，采用水平集法描述追踪界面动态变化，使模型能够有效处理这些复杂情况，提高模型对实际工程问题的适应性。例如，对于含有裂纹的粘弹性结构，改进后的模型可以准确模拟裂纹扩展过程中位移场的变化。验证基于两种技术结合的粘弹性位移场分析方法的有效性：运用数值模拟软件，对不同类型的粘弹性材料结构进行模拟分析，如不同几何形状的粘弹性板、梁结构等。通过设置多种工况，包括不同的载荷类型（如集中力、分布力、动态载荷等）、边界条件（固定边界、自由边界、弹性支撑边界等）和材料参数，全面验证所提出方法的准确性和有效性。将数值模拟结果与理论解或已有文献中的结果进行对比分析，评估方法的精度和可靠性。同时，设计并开展相关实验，选取合适的粘弹性材料制作实验试件，采用先进的测量技术（如数字图像相关技术、应变片测量技术等）获取实验数据，将实验结果与数值模拟结果相互印证，进一步验证方法的有效性。1.3.2研究方法理论分析：对时域精细算法和扩展有限元技术的基本原理进行深入研究，推导相关的数学公式和理论模型。分析粘弹性材料的本构关系，结合力学基本原理，建立基于这两种技术的粘弹性位移场分析的理论框架。例如，在推导基于时域精细算法的位移场计算公式时，运用数学分析方法，对时间域进行精细划分，推导不同时间步长下位移的递推公式；在研究扩展有限元技术时，从单位分解法的原理出发，详细推导改进函数项的构造和节点自由度的添加方法，为后续的数值模拟和实验研究提供理论依据。数值模拟：利用专业的有限元分析软件（如ABAQUS、ANSYS等），建立粘弹性材料结构的数值模型。根据理论分析的结果，设置模型的参数，包括材料参数（如弹性模量、泊松比、粘性系数等）、几何参数（结构的形状、尺寸等）和载荷条件（载荷的大小、方向、作用时间等）。通过数值模拟，得到粘弹性结构在不同工况下的位移场分布和变化规律，直观地展示所提出方法的计算结果，为方法的验证和优化提供数据支持。实验验证：设计并开展粘弹性材料结构的实验研究，通过实验测量得到结构的位移响应数据。例如，对于一个粘弹性悬臂梁结构，在其自由端施加动态载荷，利用实验设备（如激光位移传感器、应变片等）测量梁上不同位置的位移和应变。将实验结果与数值模拟结果进行对比分析，验证数值模拟方法的准确性和可靠性，同时也可以发现理论分析和数值模拟中存在的问题，为进一步改进方法提供参考。1.4研究创新点提出全新的分析模型：创新性地构建基于时域精细算法与扩展有限元技术相结合的粘弹性位移场分析模型。该模型突破了传统分析方法的局限，将时域精细算法对时间域的精确描述能力与扩展有限元技术处理间断问题和复杂边界的优势有机融合。通过详细推导模型的控制方程和改进函数项，使模型能够更准确地描述粘弹性材料在复杂工况下的位移场变化，为粘弹性材料的力学分析提供了新的理论框架。改进分析方法：在模型建立过程中，对时域精细算法和扩展有限元技术进行了改进。在时域精细算法方面，优化了时间步长的选择策略，采用自适应时间步长控制方法，根据粘弹性问题的特点和计算精度要求动态调整时间步长，提高了计算效率和稳定性。在扩展有限元技术中，提出了新的改进函数项构造方法，使其能更准确地反映裂纹、孔洞等间断特征，增强了模型对复杂间断问题的模拟能力。多场景验证方法有效性：运用多种方法对所提出的分析方法进行全面验证。通过数值模拟，利用专业有限元分析软件对不同类型、不同工况下的粘弹性材料结构进行模拟分析，设置丰富多样的参数组合，包括不同的几何形状、载荷类型、边界条件和材料参数等，得到大量的数值模拟结果。同时，设计并开展实验研究，制作粘弹性材料实验试件，采用先进的测量技术获取实验数据，将数值模拟结果与实验结果进行对比分析，相互印证方法的准确性和有效性，为方法的实际应用提供了可靠依据。二、相关理论基础2.1粘弹性材料特性粘弹性材料是一种在外力作用下，弹性和粘性两种变形机制同时存在的材料，其力学行为介于理想弹性体和理想粘性体之间。这种材料的特性使其在众多工程领域中具有重要应用，如在航空航天、汽车、土木建筑、生物医学等领域，用于减震、隔振、能量吸收和模拟生物组织力学行为等。粘弹性材料的应力-应变关系具有显著的时间依赖性和速率依赖性。与理想弹性材料不同，其应力不仅取决于当前的应变状态，还与应变的历史过程密切相关。当对粘弹性材料施加恒定应力时，材料会发生蠕变现象，即应变随时间不断增加。在一定温度下，对一块粘弹性橡胶材料施加恒定拉力，随着时间推移，橡胶会逐渐伸长，应变持续增大。而当保持材料的应变恒定时，应力会随时间逐渐减小，这就是应力松弛现象。如在固定拉伸长度的情况下，粘弹性材料所承受的拉力会逐渐降低。此外，在循环加载过程中，粘弹性材料的应力-应变曲线不重合，形成一个封闭的滞回环，这表明材料在加载和卸载过程中存在能量损耗。当对粘弹性材料进行周期性的拉伸和压缩加载时，每次加载和卸载的路径不同，形成的滞回环面积代表了能量的损耗。时间和温度是影响粘弹性材料力学性能的两个关键因素。从时间因素来看，加载速率对材料的应力-应变响应有重要影响。当加载速率较慢时，材料有足够时间发生粘性流动，表现出较大的应变；而加载速率较快时，材料来不及产生充分的粘性变形，更多地表现出弹性行为，应力-应变曲线更接近弹性材料。在低速拉伸测试中，粘弹性材料可能产生较大的塑性变形；而在高速冲击加载下，材料则更易表现出脆性断裂。温度对粘弹性材料性能的影响也十分显著。随着温度的升高，材料的粘性成分增加，弹性成分减小，表现为模量降低，阻尼增大。在高温环境下，粘弹性材料的硬度和刚度会明显下降，变得更加柔软和易于变形。当温度降低时，材料的弹性成分相对增加，粘性成分减小，模量升高，阻尼减小。在低温下，粘弹性材料可能会变得更脆，失去部分阻尼特性。而且，粘弹性材料存在玻璃化转变温度，当温度低于玻璃化转变温度时，材料处于玻璃态，分子链段运动受到限制，表现出较高的模量和较低的阻尼；当温度高于玻璃化转变温度时，材料进入橡胶态，分子链段运动较为自由，模量降低，阻尼增大。在玻璃化转变温度附近，材料的力学性能会发生急剧变化。2.2时域精细算法原理时域精细算法是一种基于时间域离散的高精度数值计算方法，其基本原理是将时间域进行精细划分，通过逐步求解离散时间点上的力学方程，来精确描述物理系统随时间的动态变化过程。该算法源于对传统数值积分方法的改进，旨在克服传统方法在处理复杂动力学问题时的精度和稳定性限制。从数学基础来看，时域精细算法基于常微分方程的数值求解理论。对于一个包含时间变量t的动力学系统，其控制方程通常可以表示为一阶常微分方程组：\frac{d\mathbf{y}(t)}{dt}=\mathbf{f}(t,\mathbf{y}(t))，其中\mathbf{y}(t)是包含系统状态变量（如位移、速度、应力等）的向量，\mathbf{f}(t,\mathbf{y}(t))是与时间和状态变量相关的函数向量。时域精细算法通过将时间区间[t_0,t_n]划分为n个小的时间步长\Deltat，即t_i=t_0+i\Deltat，i=0,1,\cdots,n，然后利用某种数值积分公式来逐步求解在每个时间步长上的\mathbf{y}(t_i)。常见的数值积分公式如Runge-Kutta方法、线性多步法等都可用于时域精细算法，但该算法在具体实现中往往采用一些特殊的技巧来提高精度和效率。例如，通过对时间步长进行自适应调整，根据问题的局部特性动态改变\Deltat的大小，在变化剧烈的区域采用较小的时间步长以保证精度，在变化平缓的区域采用较大的时间步长以提高计算效率。采用高精度的数值积分公式，并结合数值稳定性分析，确保在长时间的计算过程中误差不会累积过大。在粘弹性问题中，时域精细算法具有显著的应用优势。由于粘弹性材料的应力-应变关系与时间密切相关，传统的数值方法在处理这种时间依赖性时往往面临挑战。而时域精细算法通过对时间域的精细划分，能够精确捕捉粘弹性材料在不同时刻的力学响应。在分析粘弹性材料的蠕变过程时，该算法可以准确计算出应变随时间的逐渐增加，以及应力在不同时间点的松弛情况。时域精细算法还能有效处理粘弹性材料在动态载荷作用下的复杂响应。对于周期性变化的载荷，它能够精确计算出材料在每个周期内的应力-应变滞回曲线，准确反映材料的能量损耗特性。这对于研究粘弹性材料在振动、冲击等动态环境下的性能具有重要意义。由于其高精度和稳定性，时域精细算法可以减少计算过程中的误差积累，提高计算结果的可靠性，为粘弹性材料在工程实际中的应用提供更准确的理论依据。2.3扩展有限元技术原理扩展有限元技术（XFEM）是在传统有限元方法基础上发展起来的一种数值计算方法，其核心思想基于单位分解法。单位分解法的基本概念是，对于定义在求解域\Omega上的函数u(x)，可以通过一组具有单位分解性质的函数\varphi_{i}(x)与另一组函数u_{i}(x)的乘积之和来近似表示，即u(x)=\sum_{i=1}^{n}\varphi_{i}(x)u_{i}(x)，其中\sum_{i=1}^{n}\varphi_{i}(x)=1在\Omega上成立。在扩展有限元技术中，将这一原理应用于处理含有间断问题的结构力学分析。对于存在裂纹、孔洞、夹杂等间断情况的结构，传统有限元方法需要对间断区域进行精细的网格划分，并且在裂纹扩展等动态过程中需要不断进行网格重划分，这不仅计算成本高，而且容易引入误差。而扩展有限元技术通过在传统有限元的位移模式中增加能反映间断问题的改进函数项，来描述这些间断特征。以二维裂纹问题为例，扩展有限元的位移逼近表达式为：u_{i}(x)=\sum_{j\inN}N_{j}(x)u_{ij}+\sum_{j\inN_{1}}N_{j}(x)a_{j}H(x)+\sum_{j\inN_{2}}N_{j}(x)\sum_{\alpha=1}^{4}b_{j}^{\alpha}F^{\alpha}(x)其中，\sum_{j\inN}N_{j}(x)u_{ij}是传统有限元的位移插值部分，N_{j}(x)是传统的形函数，u_{ij}是节点j的位移自由度；\sum_{j\inN_{1}}N_{j}(x)a_{j}H(x)用于描述裂纹面的不连续，H(x)是阶跃函数，当x在裂纹的一侧时，H(x)=1，在另一侧时，H(x)=-1，a_{j}是附加的自由度，N_{1}是裂纹贯穿的节点集合；\sum_{j\inN_{2}}N_{j}(x)\sum_{\alpha=1}^{4}b_{j}^{\alpha}F^{\alpha}(x)用于描述裂尖的奇异性，F^{\alpha}(x)是与裂尖奇异性相关的函数，b_{j}^{\alpha}是相应的附加自由度，N_{2}是裂尖所在单元的节点集合。在描述和追踪结构内部的间断界面动态变化时，扩展有限元技术通常采用水平集法（LSM）。水平集法通过定义一个符号距离函数\phi(x)，使得在间断界面上\phi(x)=0，在界面一侧\phi(x)>0，在另一侧\phi(x)<0。随着间断界面的移动，\phi(x)的值也相应变化，通过求解水平集函数的演化方程，就可以追踪界面的动态变化。与传统有限元方法相比，扩展有限元技术具有显著的差异和优势。在网格划分方面，传统有限元方法需要根据结构的几何形状和内部特征进行细致的网格划分，尤其是在间断区域附近，需要加密网格以提高计算精度。而扩展有限元技术所使用的计算网格与结构内部的间断无关，不需要在间断区域进行特殊的网格处理，也无需多次进行网格重划分，大大降低了网格划分的难度和计算成本。在处理复杂间断问题的能力上，传统有限元方法对于裂纹扩展、孔洞生长等动态间断问题的模拟较为困难，需要不断更新网格和重新计算。扩展有限元技术通过改进函数项和水平集法，可以方便地模拟这些复杂的间断问题，能够更准确地描述结构在间断情况下的力学行为。三、基于时域精细算法的粘弹性位移场分析模型构建3.1模型假设与简化在构建基于时域精细算法的粘弹性位移场分析模型时，为了使复杂的问题更易于处理和分析，需要提出一系列合理的假设并进行适当的简化。假设粘弹性材料为均匀连续介质。这意味着在材料内部，其物理性质（如弹性模量、粘性系数等）在空间上是均匀分布的，不存在局部的突变或不均匀性。在实际的粘弹性材料中，虽然微观层面可能存在一定的结构差异，但在宏观分析尺度下，这种假设能够忽略微观结构的细节，简化模型的建立和分析过程。对于常见的橡胶类粘弹性材料，在研究其宏观力学性能时，将其视为均匀连续介质可以方便地应用连续介质力学的理论和方法进行分析。假设材料的变形是小变形。即材料在受力过程中产生的位移和应变都非常小，满足线性化的条件。在小变形假设下，几何方程可以线性化，从而简化了应力-应变关系的推导和计算。当粘弹性结构受到较小的外力作用时，其变形量相对结构尺寸很小，此时小变形假设成立，基于此假设建立的模型能够准确描述结构的力学行为。假设粘弹性材料的本构关系满足线性粘弹性理论。线性粘弹性理论假设应力与应变之间存在线性关系，且这种关系与时间相关。具体来说，材料的应力可以表示为应变及其对时间的各阶导数的线性组合。这种假设在一定程度上简化了本构关系的描述，使得模型能够采用相对简单的数学形式进行表达和求解。在实际应用中，许多粘弹性材料在一定的加载条件和变形范围内，其力学行为可以较好地用线性粘弹性理论来近似描述。在实际工程中，结构的几何形状和边界条件往往非常复杂。为了便于模型的构建和分析，需要对结构进行简化。对于一些复杂的几何形状，可以采用等效的简化几何模型来代替。对于具有复杂外形的粘弹性部件，可以将其简化为规则的几何形状（如矩形、圆形等），只要这种简化不会对结构的主要力学性能产生显著影响。在边界条件方面，对于一些难以精确描述的边界条件，可以采用近似的边界条件来替代。将实际的弹性支撑边界简化为固定边界或简支边界，在一定程度上可以简化计算过程，同时通过合理的参数调整，仍然能够保证模型对实际问题的模拟精度。通过以上假设和简化，我们明确了模型的适用范围。该模型适用于在小变形条件下，表现出线性粘弹性行为的均匀连续介质材料。在实际应用中，当所研究的粘弹性材料和结构满足这些条件时，基于时域精细算法构建的位移场分析模型能够准确地描述材料的力学行为，为工程设计和分析提供可靠的理论依据。然而，当材料的非线性效应、非均匀性或大变形等因素不可忽略时，该模型的准确性可能会受到影响，需要采用更复杂的理论和方法进行分析。3.2控制方程建立在基于上述假设与简化的基础上，依据力学和数学原理推导粘弹性位移场的控制方程。从基本的平衡方程出发，考虑粘弹性材料的特性，结合几何方程和本构方程，逐步推导得到位移场的控制方程。平衡方程是描述物体受力平衡状态的基本方程。对于一个微元体，在笛卡尔坐标系下，其平衡方程可表示为：\frac{\partial\sigma_{ij}}{\partialx_j}+f_i=0，其中\sigma_{ij}是应力张量，x_j是坐标方向，f_i是单位体积的体力。在粘弹性材料中，应力张量不仅与当前的应变状态有关，还与应变的历史过程相关。根据粘弹性材料的线性粘弹性理论，其本构关系可以通过积分形式来表示，即：\sigma_{ij}(t)=\int_{-\infty}^{t}G_{ijkl}(t-\tau)\frac{\partial\epsilon_{kl}(\tau)}{\partial\tau}d\tau，其中G_{ijkl}(t-\tau)是松弛模量张量，\epsilon_{kl}是应变张量，\tau是时间积分变量。松弛模量张量G_{ijkl}(t-\tau)描述了材料在不同时刻的应力松弛特性，它是一个与时间相关的函数，反映了材料的粘弹性性质。几何方程用于描述物体的变形与位移之间的关系。在小变形假设下，几何方程为：\epsilon_{ij}=\frac{1}{2}(\frac{\partialu_i}{\partialx_j}+\frac{\partialu_j}{\partialx_i})，其中u_i是位移分量。将几何方程代入本构方程中，得到：\sigma_{ij}(t)=\int_{-\infty}^{t}G_{ijkl}(t-\tau)\frac{1}{2}(\frac{\partial^2u_{kl}(\tau)}{\partialx_k\partial\tau}+\frac{\partial^2u_{lk}(\tau)}{\partialx_l\partial\tau})d\tau。再将上式代入平衡方程\frac{\partial\sigma_{ij}}{\partialx_j}+f_i=0中，经过一系列的数学推导（包括积分求导、张量运算等），得到粘弹性位移场的控制方程：\int_{-\infty}^{t}\frac{\partial}{\partialx_j}[G_{ijkl}(t-\tau)\frac{1}{2}(\frac{\partial^2u_{kl}(\tau)}{\partialx_k\partial\tau}+\frac{\partial^2u_{lk}(\tau)}{\partialx_l\partial\tau})]d\tau+f_i=0。为了便于求解，通常会对上述控制方程进行一些数学变换。引入Laplace变换，对控制方程中的时间变量t进行Laplace变换，将时域问题转化为频域问题。设U_i(s)是u_i(t)的Laplace变换，F_i(s)是f_i(t)的Laplace变换，\overline{G}_{ijkl}(s)是G_{ijkl}(t)的Laplace变换。根据Laplace变换的性质，对控制方程进行变换后得到：\frac{\partial}{\partialx_j}[\overline{G}_{ijkl}(s)\frac{1}{2}(\frac{\partial^2U_{kl}(s)}{\partialx_k}+\frac{\partial^2U_{lk}(s)}{\partialx_l})]+F_i(s)=0。这个经过Laplace变换后的控制方程在频域中更易于求解。通过求解该方程得到频域内的位移解U_i(s)，再利用Laplace逆变换将其转换回时域，即可得到粘弹性材料在时域内的位移场u_i(t)。在实际计算中，需要根据具体的材料参数确定松弛模量张量G_{ijkl}(t)的具体形式，进而确定其Laplace变换\overline{G}_{ijkl}(s)。对于不同的粘弹性模型（如Kelvin-Voigt模型、Maxwell模型等），松弛模量张量的表达式不同。在Kelvin-Voigt模型中，松弛模量张量G_{ijkl}(t)可以表示为弹性模量和粘性系数的函数，通过对其进行Laplace变换，得到相应的\overline{G}_{ijkl}(s)，代入上述控制方程进行求解。3.3时域精细算法的应用与求解过程将时域精细算法应用于上述建立的粘弹性位移场控制方程，具体求解过程如下：时间域离散：首先，将时间区间[0,T]划分为N个小的时间步长\Deltat=T/N，即t_n=n\Deltat，n=0,1,\cdots,N。在每个时间步长内，假设控制方程中的各项变量（如位移、应力等）随时间的变化是线性的或可以通过某种近似方法进行处理。建立递推关系：根据时域精细算法的原理，利用数值积分公式对控制方程进行离散求解。以常见的四阶Runge-Kutta方法为例，对于控制方程\int_{-\infty}^{t}\frac{\partial}{\partialx_j}[G_{ijkl}(t-\tau)\frac{1}{2}(\frac{\partial^2u_{kl}(\tau)}{\partialx_k\partial\tau}+\frac{\partial^2u_{lk}(\tau)}{\partialx_l\partial\tau})]d\tau+f_i=0，在第n个时间步长[t_n,t_{n+1}]上，通过Runge-Kutta方法建立位移u_{i}^{n+1}与u_{i}^{n}之间的递推关系。设F(u)表示控制方程中除位移对时间的导数项之外的其他项，即F(u)=\int_{-\infty}^{t}\frac{\partial}{\partialx_j}[G_{ijkl}(t-\tau)\frac{1}{2}(\frac{\partial^2u_{kl}(\tau)}{\partialx_k\partial\tau}+\frac{\partial^2u_{lk}(\tau)}{\partialx_l\partial\tau})]d\tau+f_i，则Runge-Kutta方法的递推公式为：k_1=\DeltatF(u^n)k_2=\DeltatF(u^n+\frac{1}{2}k_1)k_3=\DeltatF(u^n+\frac{1}{2}k_2)k_4=\DeltatF(u^n+k_3)u^{n+1}=u^n+\frac{1}{6}(k_1+2k_2+2k_3+k_4)其中k_1,k_2,k_3,k_4是中间计算量，通过逐步计算这些中间量，可以得到下一个时间步长的位移值u^{n+1}。迭代求解：从初始时刻t=0开始，根据给定的初始条件（如初始位移u_i(0)和初始速度\frac{\partialu_i(0)}{\partialt}），利用上述递推关系进行迭代计算。在每一步迭代中，需要计算控制方程中的各项系数，包括松弛模量张量G_{ijkl}(t-\tau)在不同时间点的值。由于松弛模量张量与时间相关，在计算时需要根据具体的粘弹性模型和材料参数进行确定。对于Kelvin-Voigt模型，松弛模量张量G_{ijkl}(t)可以表示为G_{ijkl}(t)=G_{ijkl}^0+\eta_{ijkl}\delta(t)，其中G_{ijkl}^0是弹性模量张量，\eta_{ijkl}是粘性系数张量，\delta(t)是Dirac函数。在时间步长计算中，根据t-\tau的值代入相应的表达式计算G_{ijkl}(t-\tau)。边界条件处理：在迭代求解过程中，需要考虑边界条件的影响。对于不同类型的边界条件（如固定边界、自由边界、弹性支撑边界等），采用相应的处理方法。在固定边界条件下，边界上的位移值为已知常数，将其代入控制方程的离散形式中，作为约束条件参与迭代计算。若结构的一端为固定边界，在该边界节点上的位移u_i在所有时间步长上都等于给定的固定值，通过这种方式将边界条件融入到迭代求解过程中，确保计算结果满足实际的边界情况。计算结果输出：经过N次迭代计算，得到在时间区间[0,T]内各个时间步长上的位移场u_i(t_n)，n=0,1,\cdots,N。这些计算结果可以以数据文件的形式保存，也可以通过绘图软件绘制出位移随时间和空间的变化曲线或云图，直观地展示粘弹性材料在不同时刻的位移分布情况。将位移场数据绘制成云图，颜色的深浅表示位移的大小，通过云图可以清晰地看到结构中位移较大和较小的区域，以及位移分布的变化趋势。四、基于扩展有限元技术的模型改进4.1扩展有限元技术在粘弹性位移场分析中的适应性分析扩展有限元技术在处理粘弹性位移场分析问题时，展现出独特的优势，同时也面临一些挑战。从优势方面来看，其处理间断问题的能力与粘弹性材料的特性高度契合。粘弹性材料在实际应用中，由于复杂的受力情况和环境因素，内部容易出现裂纹、孔洞、夹杂等间断情况。例如在航空发动机的高温部件中，粘弹性密封材料长期处于高温、高压和振动的环境下，容易产生裂纹；在土木工程的混凝土结构中，粘弹性阻尼器可能会出现孔洞或夹杂等缺陷。扩展有限元技术基于单位分解法，在传统有限元的位移模式中增加能反映间断问题的改进函数项，并附加节点自由度，能够有效地描述这些间断特征。对于裂纹问题，通过引入阶跃函数和裂尖渐近位移场函数，能够准确模拟裂纹面的不连续和裂尖的奇异性，从而精确分析裂纹对粘弹性位移场的影响。在模拟含有裂纹的粘弹性材料板的受力情况时，扩展有限元技术可以清晰地展示裂纹扩展过程中位移场的变化，为评估材料的性能和寿命提供重要依据。扩展有限元技术采用水平集法描述追踪界面动态变化，这使得它在处理粘弹性材料内部界面动态变化问题时具有显著优势。粘弹性材料在加载、卸载过程中，其内部的微观结构可能会发生变化，导致界面动态变化。在循环加载下，粘弹性材料内部的分子链结构会发生重排，从而引起界面的移动和变化。水平集法通过定义一个符号距离函数，能够准确追踪这种界面的动态变化，为分析粘弹性位移场提供了更准确的边界条件。通过水平集法可以实时监测粘弹性材料在不同加载阶段的界面变化，进而分析其对位移场的影响。然而，扩展有限元技术在应用于粘弹性位移场分析时也面临一些挑战。计算效率是一个重要问题。由于扩展有限元技术需要在传统有限元的基础上增加改进函数项和附加节点自由度，导致计算量大幅增加。在处理大规模粘弹性结构时，求解方程组的规模增大，计算时间显著增加。在分析大型粘弹性桥梁结构时，计算时间可能会延长数倍甚至数十倍，这对于实际工程应用来说是一个较大的限制。在处理复杂的粘弹性本构关系时，扩展有限元技术也存在一定的困难。粘弹性材料的本构关系通常较为复杂，与时间、温度等因素密切相关。将这些复杂的本构关系融入扩展有限元的计算框架中，需要进行大量的数学推导和数值计算，增加了计算的复杂性和难度。对于一些考虑多物理场耦合（如热-力耦合、流-固耦合）的粘弹性问题，如何准确地将扩展有限元技术与多物理场的控制方程相结合，也是目前面临的一个挑战。4.2改进模型的建立与离散化基于扩展有限元技术对原有的基于时域精细算法的粘弹性位移场分析模型进行改进。在改进模型中，针对粘弹性材料内部可能出现的裂纹、孔洞、夹杂等间断情况，引入相应的改进函数项来描述这些间断特征。对于裂纹问题，在位移模式中增加阶跃函数和裂尖渐近位移场函数。设裂纹将求解域划分为两个子域，阶跃函数H(x)定义为：当x在裂纹的一侧时，H(x)=1；当x在另一侧时，H(x)=-1。裂尖渐近位移场函数F^{\alpha}(x)（\alpha=1,2,3,4）用于描述裂尖的奇异性，其具体形式与裂纹尖端的应力应变场的渐近展开式相关。以二维裂纹问题为例，改进后的位移逼近表达式为：u_{i}(x)=\sum_{j\inN}N_{j}(x)u_{ij}+\sum_{j\inN_{1}}N_{j}(x)a_{j}H(x)+\sum_{j\inN_{2}}N_{j}(x)\sum_{\alpha=1}^{4}b_{j}^{\alpha}F^{\alpha}(x)其中，\sum_{j\inN}N_{j}(x)u_{ij}是传统有限元的位移插值部分，N_{j}(x)是传统的形函数，u_{ij}是节点j的位移自由度；\sum_{j\inN_{1}}N_{j}(x)a_{j}H(x)用于描述裂纹面的不连续，a_{j}是附加的自由度，N_{1}是裂纹贯穿的节点集合；\sum_{j\inN_{2}}N_{j}(x)\sum_{\alpha=1}^{4}b_{j}^{\alpha}F^{\alpha}(x)用于描述裂尖的奇异性，b_{j}^{\alpha}是相应的附加自由度，N_{2}是裂尖所在单元的节点集合。对于孔洞问题，通过在孔洞边界附近的单元中增加特殊的改进函数项来模拟孔洞对位移场的影响。设孔洞边界的方程为\Gamma，可以定义一个与孔洞边界相关的函数G(x)，当x靠近孔洞边界时，G(x)的值发生变化，从而反映孔洞的存在对位移的影响。在含有孔洞的粘弹性材料模型中，在孔洞边界附近的单元位移模式中添加与G(x)相关的项，如\sum_{j\inN_{hole}}N_{j}(x)c_{j}G(x)，其中N_{hole}是孔洞边界附近单元的节点集合，c_{j}是附加的自由度。在夹杂问题中，考虑夹杂与基体材料的界面特性，通过在界面附近的单元中引入反映界面特性的改进函数项。设夹杂与基体的界面为\Sigma，定义一个与界面相关的函数I(x)，当x在界面附近时，I(x)体现出夹杂与基体之间的力学相互作用。在界面附近单元的位移模式中增加与I(x)相关的项，如\sum_{j\inN_{interface}}N_{j}(x)d_{j}I(x)，其中N_{interface}是界面附近单元的节点集合，d_{j}是附加的自由度。在离散化处理方面，采用有限元网格对改进后的模型进行离散。与传统有限元方法不同的是，扩展有限元技术所使用的计算网格与结构内部的间断无关，不需要在间断区域进行特殊的网格处理。在划分网格时，按照常规的有限元网格划分方法进行，无需考虑裂纹、孔洞、夹杂等间断的具体位置和形状。采用三角形单元或四边形单元对求解域进行划分，确定单元的节点坐标和连接关系。在确定节点自由度时，除了传统有限元的位移自由度外，还需要考虑由于引入改进函数项而增加的附加自由度，如a_{j}、b_{j}^{\alpha}、c_{j}、d_{j}等。将整个求解域离散为有限个单元后，通过节点的连接和自由度的协调，建立起离散化的方程组，为后续的求解奠定基础。4.3求解方法与收敛性分析针对改进后的模型，选择合适的求解方法进行求解。常用的求解方法包括直接解法和迭代解法。直接解法如高斯消去法，它通过一系列的初等行变换将线性方程组的系数矩阵化为上三角矩阵，然后通过回代过程求解方程组。对于规模较小的方程组，高斯消去法具有计算速度快、精度高的优点。在处理小型粘弹性结构模型时，使用高斯消去法可以快速得到精确的位移场解。然而，对于大规模的方程组，高斯消去法的计算量和存储量会随着方程组规模的增大而急剧增加，导致计算效率低下。迭代解法如共轭梯度法，它通过迭代的方式逐步逼近方程组的解。共轭梯度法的基本思想是构造一组共轭方向，在这些方向上进行搜索，使得目标函数（通常是残差的平方和）逐步减小。该方法适用于大规模稀疏矩阵方程组的求解，因为它不需要存储整个系数矩阵，只需要存储矩阵与向量的乘积，从而大大减少了存储量。在处理大型粘弹性结构时，共轭梯度法能够有效地求解大规模的方程组，提高计算效率。共轭梯度法的收敛速度与系数矩阵的条件数有关，条件数越大，收敛速度越慢。对改进模型的收敛性进行分析。收敛性是指随着计算过程的进行，解是否逐渐逼近真实解。通过数值实验和理论推导来评估收敛性。在数值实验中，固定其他参数，逐步增加有限元网格的密度，观察计算结果（如位移场）的变化。如果随着网格密度的增加，计算结果逐渐稳定，且与理论解或精确解的误差逐渐减小，则说明模型具有收敛性。当网格单元尺寸减小到原来的一半时，位移场的计算结果与上一次计算结果的相对误差小于一定的阈值（如0.01），则可以认为模型在该工况下收敛。从理论推导的角度，根据扩展有限元技术的理论基础和改进函数项的性质，分析模型的收敛性。对于扩展有限元中增加的反映间断问题的改进函数项，证明其在一定条件下能够保证位移场的逼近精度随着网格细化而提高。通过数学分析，证明改进后的位移逼近表达式在满足一定的数学条件时，能够收敛到真实的位移场。收敛性分析对于确保计算结果的准确性和可靠性至关重要。只有当模型收敛时，计算得到的位移场才能够真实地反映粘弹性材料在复杂间断情况下的力学行为，为工程设计和分析提供可靠的依据。五、数值算例与结果分析5.1算例选取与模型建立为了验证基于时域精细算法与扩展有限元技术的粘弹性位移场分析方法的有效性和准确性，精心选取带中心裂纹的粘弹性平面薄板作为典型算例。此算例在工程实际中具有广泛的代表性，如航空航天领域中飞行器的机翼结构、汽车制造中的车身板件等，在服役过程中都可能出现类似的裂纹缺陷，对其位移场进行准确分析对于保障结构的安全性和可靠性至关重要。利用专业有限元分析软件ABAQUS建立有限元模型。该软件具有强大的建模和求解功能，能够准确模拟复杂结构的力学行为。模型的几何参数设定为：薄板的长度L=100\text{mm}，宽度W=50\text{mm}，厚度t=5\text{mm}，中心裂纹长度a=20\text{mm}。这些参数的选择综合考虑了实际工程中常见的尺寸范围以及计算的可行性。材料参数方面，根据实际应用中常见的粘弹性材料特性进行设定。弹性模量E=2\times10^{9}\text{Pa}，它反映了材料抵抗弹性变形的能力；泊松比\nu=0.3，用于描述材料在受力时横向应变与纵向应变的关系；粘性系数\eta=1\times10^{5}\text{Pa}\cdot\text{s}，体现了材料的粘性特征，对材料的时间相关力学行为有重要影响。这些材料参数是通过对实际材料的实验测试和数据分析得到的，具有一定的工程实际意义。在模型的边界条件设置上，将薄板的左右两侧边界约束为水平方向位移为零，即u_x=0。这模拟了实际结构中在某些支撑条件下，该方向的位移受到限制的情况。上下两侧边界约束为垂直方向位移为零，即u_y=0，以反映实际结构中可能存在的类似约束条件。在薄板的上表面施加均布载荷q=1000\text{N/mm}^2，模拟实际工程中可能受到的压力作用。通过合理设置这些边界条件和载荷，能够更真实地模拟粘弹性平面薄板在实际工况下的受力状态，为后续的位移场分析提供可靠的基础。5.2计算结果与分析通过ABAQUS软件对带中心裂纹的粘弹性平面薄板模型进行求解计算，得到了丰富的计算结果，这些结果从多个角度展示了粘弹性平面薄板在给定工况下的力学行为，对分析模型的性能和验证方法的有效性具有重要意义。5.2.1位移场结果位移场分布云图直观地展示了薄板在均布载荷作用下的位移情况，如图1所示。从云图中可以清晰地看出，在裂纹尖端附近，位移呈现出明显的集中现象，位移值较大。这是因为裂纹的存在破坏了材料的连续性，使得应力在裂纹尖端集中，从而导致该区域的位移显著增大。在远离裂纹的区域，位移分布相对较为均匀，且位移值逐渐减小。这表明裂纹对位移场的影响主要集中在其尖端附近，随着距离裂纹的增加，影响逐渐减弱。为了更深入地分析位移随时间的变化规律，选取了裂纹尖端处的一点进行研究，绘制了该点的位移-时间曲线，如图2所示。从曲线中可以看出，随着时间的增加，位移逐渐增大，呈现出典型的粘弹性材料的蠕变特性。在加载初期，位移增长速度较快，这是由于材料在突然受到载荷作用时，弹性变形迅速发生，同时粘性变形也开始逐渐发展。随着时间的推移，位移增长速度逐渐减缓，这是因为粘性变形的发展受到材料内部结构的阻碍，需要一定的时间来克服这种阻碍。最终，位移趋于稳定，达到一个相对稳定的值，此时材料的粘性变形和弹性变形达到了一种平衡状态。与传统有限元方法的计算结果进行对比，进一步验证本文方法的准确性。在相同的模型参数和工况条件下，采用传统有限元方法进行计算，并将其位移场计算结果与本文基于时域精细算法与扩展有限元技术的方法所得结果进行对比。对比结果表明，本文方法在处理粘弹性材料的位移场分析时，能够更准确地捕捉裂纹尖端的位移集中现象和位移随时间的变化规律。传统有限元方法在处理裂纹问题时，由于其对间断问题的处理能力有限，往往会在裂纹尖端附近产生较大的误差，导致位移计算结果不准确。而本文方法通过引入扩展有限元技术，能够有效地描述裂纹的间断特征，结合时域精细算法对时间域的精确处理，使得位移场的计算结果更加准确可靠。5.2.2应力分布结果应力分布云图展示了薄板内部的应力分布情况，如图3所示。在裂纹尖端附近，应力集中现象非常明显，应力值远高于其他区域。这是由于裂纹尖端的几何形状导致应力在该区域高度集中，容易引发材料的破坏。在薄板的其他部位，应力分布相对较为均匀，但在靠近边界和加载区域，应力也会出现一定程度的变化。在加载区域，由于直接承受均布载荷的作用，应力值较大；在边界约束区域，由于位移受到限制，也会产生一定的应力集中。分析不同时刻的应力分布变化，发现随着时间的推移，应力分布逐渐发生变化。在加载初期，应力主要集中在加载区域和裂纹尖端附近。随着时间的增加，应力逐渐向周围扩散，裂纹尖端的应力集中程度有所降低，但仍然是应力的高值区域。这是因为粘弹性材料的应力松弛特性，使得应力在材料内部逐渐重新分布。在应力松弛过程中，材料内部的分子结构逐渐调整，以适应外力的作用，从而导致应力分布的变化。与理论解进行对比，验证应力计算结果的准确性。根据弹性力学和粘弹性理论，对于带中心裂纹的薄板在均布载荷作用下的应力分布，存在相应的理论解。将本文方法计算得到的应力结果与理论解进行对比，结果显示二者在趋势上基本一致，在数值上也较为接近。在裂纹尖端附近，由于理论解通常基于一些简化假设，与实际计算结果存在一定的差异，但本文方法的计算结果能够更真实地反映裂纹尖端的复杂应力状态。这表明本文基于时域精细算法与扩展有限元技术的方法在计算粘弹性材料的应力分布时具有较高的准确性和可靠性，能够为工程实际提供有效的理论支持。5.3与传统方法对比验证为了更全面地验证基于时域精细算法与扩展有限元技术的粘弹性位移场分析方法的优势，将其与传统的有限元方法和边界元方法进行详细对比。在计算精度方面，针对带中心裂纹的粘弹性平面薄板算例，采用传统有限元方法进行计算时，由于其对裂纹等间断问题的处理能力有限，在裂纹尖端附近，位移和应力的计算结果存在较大误差。传统有限元方法在模拟裂纹尖端的应力集中现象时，无法准确捕捉应力的急剧变化，导致计算得到的应力值与实际值偏差较大。而基于时域精细算法与扩展有限元技术的方法，通过引入扩展有限元技术，能够准确描述裂纹的间断特征，结合时域精细算法对时间域的精确处理，使得位移和应力的计算结果与理论解和实验结果更为接近。在计算位移时，新方法计算得到的裂纹尖端处的位移值与理论解的相对误差在可接受范围内，而传统有限元方法的相对误差则明显较大。从计算效率来看，传统边界元方法在处理粘弹性问题时，由于其只能处理绝对边界，对于复杂边界条件的问题，需要进行大量的边界离散和积分计算，导致计算时间大幅增加。在处理具有不规则边界的粘弹性结构时，边界元方法的计算时间是新方法的数倍。基于时域精细算法与扩展有限元技术的方法所使用的计算网格与结构内部的间断无关，不需要在间断区域进行特殊的网格处理，也无需多次进行网格重划分，大大降低了计算成本，提高了计算效率。在相同的计算条件下，新方法的计算时间明显短于传统边界元方法，能够更快地得到计算结果。在处理复杂边界和间断问题的能力上，传统有限元方法和边界元方法都存在一定的局限性。传统有限元方法需要对间断区域进行精细的网格划分，并且在裂纹扩展等动态过程中需要不断进行网格重划分，这不仅计算成本高，而且容易引入误差。边界元方法则无法处理任意边界，对于含有裂纹、孔洞、夹杂等间断情况的结构，难以准确模拟其力学行为。而基于时域精细算法与扩展有限元技术的方法，能够有效地处理各种复杂边界和间断问题，准确模拟粘弹性材料在这些情况下的位移场和应力场分布。对于含有多个裂纹和孔洞的粘弹性结构，新方法能够清晰地展示裂纹和孔洞对位移场和应力场的影响，而传统方法则难以准确模拟。通过以上对比分析，充分验证了基于时域精细算法与扩展有限元技术的粘弹性位移场分析方法在计算精度、计算效率和处理复杂问题能力等方面的优越性。六、实验验证与分析6.1实验方案设计为了对基于时域精细算法与扩展有限元技术的粘弹性位移场分析方法进行全面验证，精心设计实验方案。在材料选择上，选用丁腈橡胶作为实验材料。丁腈橡胶是一种典型的粘弹性材料，具有良好的耐油性、耐磨性和抗老化性能，在工业领域应用广泛，其粘弹性特性较为明显，适合用于本次实验研究。通过对丁腈橡胶的实验分析，能够更直观地验证所提出方法在实际粘弹性材料中的应用效果。实验设备选用高精度的电子万能试验机，它能够精确控制加载力的大小和加载速率，为实验提供稳定可靠的加载条件。配备高精度的激光位移传感器，用于测量粘弹性材料在加载过程中的位移变化。激光位移传感器具有高精度、非接触式测量的特点，能够准确获取材料表面的位移信息，减少测量误差。测量方法采用非接触式测量与应变片测量相结合的方式。利用激光位移传感器对丁腈橡胶试件表面的位移进行非接触式测量，获取试件在不同位置和不同加载阶段的位移数据。在试件表面关键部位粘贴高精度应变片，测量试件的应变情况。通过应变与位移之间的关系，进一步验证位移测量的准确性。实验步骤如下：首先，将丁腈橡胶加工成尺寸为长50\text{mm}、宽20\text{mm}、厚5\text{mm}的矩形试件，以保证实验的可重复性和数据的准确性。将加工好的试件安装在电子万能试验机的夹具上，确保试件安装牢固，受力均匀。在试件表面的特定位置粘贴应变片，并连接好应变测量设备，确保应变片能够正常工作。在试件周围布置激光位移传感器，调整传感器的位置和角度，使其能够准确测量试件表面的位移。通过电子万能试验机对试件施加拉伸载荷，加载速率设定为1\text{mm/min}，按照预设的加载方案逐步增加载荷，记录每个加载阶段的载荷值、应变值和位移值。在加载过程中，实时采集应变片和激光位移传感器的数据，确保数据的完整性和准确性。加载完成后，对采集到的数据进行整理和分析，与数值模拟结果进行对比，评估基于时域精细算法与扩展有限元技术的粘弹性位移场分析方法的准确性和可靠性。6.2实验过程与数据采集实验过程严格按照既定方案有条不紊地进行。在试件安装环节，将丁腈橡胶试件小心地安装在电子万能试验机的夹具上，确保试件的安装位置准确无误，夹具对试件的夹持力度均匀，避免因安装不当导致试件受力不均，影响实验结果的准确性。在安装应变片时，先对试件表面进行清洁处理，去除表面的杂质和油污，然后使用专用的胶水将应变片精确地粘贴在预定位置，确保应变片与试件表面紧密贴合，连接导线时注意避免导线松动或短路，保证应变片能够稳定地传输信号。正式加载前，再次检查电子万能试验机、激光位移传感器和应变测量设备的工作状态，确保设备正常运行。设定电子万能试验机的加载速率为1\text{mm/min}，按照预设的加载方案，以一定的载荷增量逐步增加载荷。在加载过程中，密切关注电子万能试验机的运行情况，确保加载过程平稳、无异常波动。同时，实时监测激光位移传感器和应变片采集的数据，观察数据的变化趋势，判断数据是否合理。若发现数据异常，立即暂停加载，检查设备和试件的状态，排查问题原因，确保数据采集的可靠性。在整个实验过程中，准确采集不同载荷下的位移和应变数据。每隔一定的载荷增量，记录一次激光位移传感器测量的位移数据和应变片测量的应变数据。为了确保数据的准确性，每个载荷点的数据采集多次，取平均值作为该载荷点的测量值。在采集位移数据时，激光位移传感器的测量精度能够达到微米级，保证了位移测量的高精度。对于应变数据，应变片的测量精度也能满足实验要求，通过对应变片输出信号的精确测量和处理，得到准确的应变值。实验过程中，仔细观察并记录丁腈橡胶试件的变形情况和破坏现象。在加载初期，试件呈现出弹性变形特征，随着载荷的逐渐增加，试件开始出现明显的粘弹性变形，变形量随时间逐渐增大。当载荷达到一定程度时，试件表面开始出现细微的裂纹，随着载荷的进一步增加，裂纹逐渐扩展，