基于刚度视角下桥梁船撞力的数值模拟与试验验证研究

基于刚度视角下桥梁船撞力的数值模拟与试验验证研究一、引言1.1研究背景与意义随着经济全球化和内河航运的迅猛发展，桥梁作为水陆交通的重要枢纽，其数量和规模不断增加，在跨越江河、湖泊和海峡等水域时，桥梁面临着船舶碰撞的潜在风险。近年来，桥梁船撞事故频发，给人民生命财产安全带来了巨大威胁，也对社会经济发展造成了严重影响。据不完全统计，仅在过去四年间，我国就至少发生了19起船撞桥事件，这些事故不仅导致了桥梁结构的严重破坏，部分桥梁甚至完全倒塌，还造成了人员伤亡和航道受阻等一系列问题。例如，2024年2月22日5时30分左右，一艘集装箱船空载从佛山南海开往广州南沙途中，航经洪奇沥水道时触碰沥心沙大桥桥墩，致沥心沙大桥桥面断裂，事故造成5人死亡。2007年6月15日，“南桂机035”轮从佛山高明开往顺德途中偏离主航道，触碰九江大桥非通航孔的桥墩，造成九江大桥三个桥墩倒塌，其所承桥面约200米坍塌，“南桂机035”轮沉没，有4辆汽车坠入河中，造成9人失踪。这些惨痛的事故案例表明，桥梁船撞问题已成为一个亟待解决的重要课题。桥梁船撞事故的发生，不仅会导致桥梁结构的局部损坏，如桥墩开裂、变形，还可能引发桥梁整体的失稳和倒塌，严重危及桥上车辆和行人的安全，对过往船舶也会造成不同程度的损害，影响船舶的航行安全和运营效率。船撞事故还会导致航道堵塞，使水上交通陷入瘫痪，给内河航运带来巨大的经济损失。据相关统计，一次严重的船撞桥事故，其直接经济损失可达数千万元甚至更高，间接经济损失更是难以估量。此外，船撞事故还可能对周边环境造成污染，如船舶燃油泄漏导致水体污染，破坏生态平衡。研究桥梁船撞力对于保障桥梁和船舶的安全具有至关重要的意义。准确掌握桥梁船撞力的大小和分布规律，能够为桥梁的抗撞设计提供科学依据，使桥梁在设计阶段就能充分考虑船撞荷载的作用，合理选择结构形式和材料，提高桥梁的抗撞能力，从而有效降低船撞事故发生时桥梁的损坏程度，保障桥上人员和车辆的安全。研究桥梁船撞力还有助于制定合理的船舶航行规则和桥梁保护措施，加强对桥区水域的交通管理，减少船撞事故的发生概率。通过对船撞力的研究，还可以为事故后的桥梁检测、评估和修复提供技术支持，加快桥梁的修复进程，减少事故对交通和经济的影响。1.2国内外研究现状国外对于桥梁船撞力的研究起步较早，20世纪80年代初，国际上就开始对船桥碰撞问题进行系统研究。1983年，IABSE（InternationalAssociationofBridgeandStructuralEngineering）在哥本哈根召开了一次国际会议讨论此问题，这是对船桥碰撞问题的第一次研讨会。此后，欧美等国家和地区在桥梁船撞力的理论研究、数值仿真和试验研究等方面取得了一系列重要成果。在理论研究方面，Minorsky船-船碰撞理论自1975年公开发表后，为船-桥碰撞的分析奠定了基础，该理论将船-船碰撞问题分为动能损失和结构损伤两个相互独立的部分，并通过统计分析方法将它们联系起来。汉斯和德鲁彻教授提出的理论，将船舶碰撞桥墩及其防撞设施等效成一个弹簧质量系统的数学模型，用于计算碰撞中桥墩或防护系统受撞位置处的最大位移、船舶的最大加速度、船舶的最大撞击力以及撞击过程的持续时间。在数值仿真方面，随着计算机技术和有限元方法的发展，国外学者利用非线性有限元软件对船桥碰撞过程进行了大量的模拟研究，能够较为准确地预测船撞力的大小和变化规律，分析船舶和桥梁在碰撞过程中的动力响应和损伤机制。在试验研究方面，国外开展了一些船桥碰撞的模型试验和足尺试验，通过实际测量碰撞力、变形等参数，验证理论和数值计算的结果，为船撞力的研究提供了重要的实验数据支持。国内对桥梁船撞力的研究相对较晚，但近年来也取得了显著的进展。随着我国内河航运的快速发展和桥梁建设的日益增多，船桥碰撞问题受到了越来越多的关注。国内学者在借鉴国外研究成果的基础上，结合我国的实际情况，开展了一系列的研究工作。在理论研究方面，对国外的船撞力理论进行了深入的分析和研究，并提出了一些适合我国国情的船撞力计算方法和公式。在数值仿真方面，利用ANSYS/LS-DYNA、ABAQUS等大型有限元软件，对不同类型的桥梁和船舶进行了碰撞模拟，研究了船撞力的影响因素，如船舶速度、吨位、碰撞角度等，以及桥梁结构的抗撞性能。在试验研究方面，一些高校和科研机构开展了船桥碰撞的缩尺模型试验，通过测量模型在碰撞过程中的应变、加速度等参数，分析船撞力的传递规律和桥梁结构的损伤模式。尽管国内外在桥梁船撞力的研究方面取得了一定的成果，但仍存在一些不足之处。现有研究中对于船桥碰撞过程中复杂的非线性行为，如材料非线性、几何非线性和接触非线性等，考虑还不够全面和深入，导致计算结果与实际情况存在一定的偏差。不同的船撞力计算方法和公式之间存在较大的差异，缺乏统一的标准和规范，使得在实际工程应用中难以选择合适的方法进行计算。在试验研究方面，由于船桥碰撞试验的成本较高、难度较大，试验数据相对较少，难以全面验证理论和数值计算的结果。此外，对于新型桥梁结构和船舶类型的船撞力研究还相对薄弱，不能满足工程实际的需求。1.3研究内容与方法本研究旨在基于刚度的角度，深入开展桥梁船撞力的数值仿真与试验研究，具体研究内容包括以下几个方面：桥梁和船舶的刚度分析：考虑物理性质和力学因素对桥梁和船舶结构的影响，建立合理的数值模型，对桥梁和船舶进行刚度分析，明确其在碰撞过程中的刚度特性，为后续的船撞力计算和分析奠定基础。数值仿真模拟：运用数值仿真方法，基于刚度理论建立桥梁和船舶的碰撞模型，模拟不同条件下桥梁和船舶的相互作用过程，分析其在碰撞过程中的受力情况，测算撞击力大小、撞毁桥梁的危害程度等，研究船撞力的影响因素和变化规律。试验研究：根据数值仿真结果，针对不同条件、桥梁和船舶的特点，设计并开展相应的试验研究，选择合适的实验方法，如落锤试验等，对数值仿真的可靠性和准确性进行验证和修正，通过试验获取实际的碰撞数据，为理论分析和数值模拟提供依据。结果分析与优化建议：根据仿真和实验结果，深入分析桥梁和船舶的刚度特点、结构缺陷以及制约因素等，提出针对性的优化建议和改进方案，为预防桥梁船碰撞事故提供科学的参考依据，提高桥梁和船舶的安全性。本研究采用以下研究方法：数值仿真法：根据刚度理论，利用大型通用有限元软件，如ANSYS/LS-DYNA等，建立桥梁和船舶的三维有限元模型，模拟船桥碰撞的动态过程，考虑材料非线性、几何非线性和接触非线性等因素，计算船撞力的大小、作用时间和分布规律，以及桥梁和船舶结构的应力、应变和变形等响应。试验研究法：设计并进行船桥碰撞的物理模型试验，采用相似理论设计试验模型，通过测量碰撞过程中的力、位移、加速度等物理量，获取实际的碰撞数据，与数值仿真结果进行对比分析，验证数值模型的准确性和可靠性。理论分析法：对船桥碰撞的基本理论进行深入研究，分析船撞力的简化计算公式，结合刚度理论，推导考虑碰撞体刚度的船撞力修正公式，从理论上揭示船撞力的本质和影响因素，为数值仿真和试验研究提供理论支持。二、桥梁船撞力相关理论基础2.1船桥碰撞基本理论2.1.1Minorsky理论Minorsky理论由美国学者V.U.Minorsky于1959年提出，其最初是针对核动力军舰的设计、核电站的防护以及其它海洋结构物的防撞保护问题，通过对26起船船碰撞事件的研究，得出了变形的刚体积和吸收冲击能之间的一种线性关系，即著名的Minorsky曲线。该理论按照能量守恒定律，认为两艘船舶在碰撞时，速度会发生突然改变，同时船舶的结构也会发生变化。从能量转移的角度来看，碰撞时会发生动能的改变，部分动能消耗于碰撞后两艘船舶的运动和碰撞冲击时周围水的运动，剩余部分的动能则由结构的弹性和塑性变形或结构撕裂所吸收，进而提出了碰撞系统能量吸收系数的概念。在船桥碰撞中，Minorsky理论将船桥碰撞视为一个能量转换的过程。船舶的动能在碰撞过程中转化为桥梁和船舶结构的变形能、碰撞后船舶和桥梁的剩余动能以及碰撞冲击时周围水的运动能量。通过对碰撞过程中能量的分析，可以估算船撞力的大小和碰撞对结构造成的损伤程度。该理论为船桥碰撞的分析奠定了基础，其能量分析的方法为后续研究提供了重要的思路和框架，被国际桥梁工程界和各国学者广泛认可，并在船桥碰撞研究中得到了大量的应用和验证。然而，Minorsky理论也存在一定的局限性。该理论是基于船船碰撞的统计数据得出的，对于船桥碰撞这种具有不同结构特点和碰撞模式的情况，其适用性可能受到一定限制。它没有充分考虑桥梁和船舶结构的具体力学特性，如材料的非线性、几何形状对碰撞力分布的影响等。在实际的船桥碰撞中，桥梁和船舶的结构响应往往非常复杂，涉及到多种非线性因素，而Minorsky理论在处理这些复杂的非线性行为时存在不足，导致计算结果与实际情况可能存在一定的偏差。该理论对碰撞过程中周围水的影响考虑相对简单，实际船桥碰撞中，水的流体动力作用对碰撞力和结构响应有重要影响，这也是Minorsky理论需要进一步完善的地方。2.1.2沃辛碰撞理论沃辛碰撞理论由德国学者沃辛（G.Woisin）在1967-1976年间提出，其背景是针对保护核动力船的反应堆不受其它船舶撞击而进行了若干次高能动力船舶撞击模型试验。沃辛在这些试验的基础上，充实和发展了Minorsky理论，并将其应用于船桥碰撞的分析当中。该理论的要点主要是在Minorsky理论的基础上，进一步考虑了船舶的结构特性和碰撞过程中的动力学响应。通过试验数据的分析，沃辛提出了一些修正系数和经验公式，用于更准确地计算船撞力的大小和碰撞过程中的能量转化。在计算船撞力时，考虑了船舶的质量、速度、撞击角度以及船舶结构的刚度等因素，通过对这些因素的综合分析，得出了更符合实际情况的船撞力计算公式。沃辛碰撞理论对船撞力分析的作用主要体现在以下几个方面。它为船撞力的计算提供了更具体的方法和公式，使得在实际工程中能够更方便地估算船撞力的大小，为桥梁的抗撞设计提供了更有针对性的参考。该理论通过对船舶结构特性的考虑，更深入地揭示了船桥碰撞的力学机制，有助于理解碰撞过程中船舶和桥梁的相互作用，为进一步研究船桥碰撞的动态响应和损伤模式提供了理论支持。沃辛碰撞理论的试验研究方法也为后续的船桥碰撞研究提供了借鉴，通过模型试验获取的数据和经验，为数值仿真和理论分析提供了验证和补充。2.1.3汉斯—德鲁彻理论汉斯—德鲁彻理论是汉斯和德鲁彻教授根据CG—71955—A合同研究提出的，主要用于研究公路桥梁预防船舶的撞击。该理论的核心内容是将船舶碰撞桥墩及其防撞设施等效成一个弹簧质量系统的数学模型。在这个模型中，船舶被视为具有一定质量和速度的物体，桥墩及其防撞设施则被等效为一个具有一定刚度和阻尼的弹簧系统。通过这个弹簧质量系统的数学模型，可以计算碰撞中桥墩或防护系统受撞位置处的最大位移、船舶的最大加速度、船舶的最大撞击力以及撞击过程的持续时间。在计算最大位移时，根据弹簧的弹性变形原理，结合船舶的动能和弹簧系统的刚度，推导出受撞位置处的位移计算公式；计算最大加速度时，利用牛顿第二定律，考虑船舶受到的碰撞力和自身质量，得出船舶的加速度；计算最大撞击力时，则通过分析碰撞过程中弹簧系统的受力情况和能量转化关系来确定；而撞击过程的持续时间则可以根据碰撞力和船舶运动状态的变化来估算。在船桥碰撞研究中，汉斯—德鲁彻理论具有重要的应用价值。它为桥梁防撞设计提供了一种量化的分析方法，通过对碰撞过程中各种参数的计算，可以评估桥墩和防撞设施的抗撞性能，为防撞设施的设计和优化提供科学依据。在桥梁的安全评估和风险分析中，该理论可以帮助确定船舶碰撞对桥梁结构的影响程度，预测可能发生的破坏模式，从而制定相应的防护措施和应急预案。汉斯—德鲁彻理论的弹簧质量系统模型相对简单直观，便于理解和应用，在实际工程中得到了广泛的应用。但该理论也存在一些局限性，例如对实际船桥碰撞中复杂的非线性因素考虑不够全面，模型的准确性在一定程度上依赖于对弹簧刚度和阻尼等参数的合理选取。2.2船撞力简化公式2.2.1常见规范公式介绍Woisin公式：由德国学者沃辛提出，公式为P_{t}=0.88\sqrt{DWT}\pm50，其中P_{t}为船撞力（kN），DWT为船舶载重吨位（t）。该公式是通过对高能动力船舶撞击模型试验的数据进行分析得出的，主要考虑了船舶载重吨位对船撞力的影响。Pedersen公式：该公式考虑了船舶的动能和结构变形能，公式形式较为复杂，涉及到船舶的质量、速度、结构刚度以及碰撞持续时间等多个参数。其基本思路是基于能量守恒原理，将船舶碰撞过程中的能量转化进行量化分析，从而计算出船撞力的大小。在实际应用中，需要准确获取船舶和桥梁的相关参数，才能较为准确地使用该公式计算船撞力。AASHTO规范公式：美国公路与运输协会（AASHTO）规范公式为P_{t}=0.122\sqrt{VDWT}\pm50，其中V为船舶撞击速度（m/s），DWT为船舶载重吨位（t）。该公式在Woisin公式的基础上，进一步考虑了船舶撞击速度对船撞力的影响，认为船撞力与船舶载重吨位的平方根以及撞击速度的平方根成正比。IABSE规范公式：即修正的沃辛公式，为P_{max}=0.024V^{\frac{3}{2}}D，其中P_{max}为最大船撞力（kN），V为船舶撞击速度（m/s），D为船舶排水量（t）。该公式对Woisin公式进行了修正，引入了船舶撞击速度的三次方根和排水量的因素，使得计算结果更能反映实际情况。Euro-code1规范公式：对于内陆航道的船舶，公式为P=KV，其中K=5MN/m，V=3m/s；对于远洋船舶，K=15MN/m，V=3m/s，K指的是船舶的刚度。该公式根据船舶的类型（内陆航道船舶和远洋船舶），分别给定了不同的船舶刚度系数K，通过船舶速度与刚度系数的乘积来计算船撞力。我国公路桥涵设计通用规范公式：漂流物和排筏的平均撞击力公式为P=\frac{WV}{gt}，其中W为漂流物或排筏的重力（kN），V为漂流物或排筏的流速（m/s），g为重力加速度（m/s^{2}），t为撞击时间（s）。该公式主要用于计算漂流物和排筏对桥梁的撞击力，考虑了漂流物或排筏的重力、流速以及撞击时间等因素。我国铁路桥涵设计基本规范公式：公式为P=\frac{W\gamma}{C_{1}+C_{2}}\sin\alpha，动能折减系数\gamma正向撞击时采用0.3，斜向撞击时采用0.2；没有实际资料时可以取C_{1}+C_{2}=0.0005kN/m，对于刚度较大的桥墩，可取C_{2}=0，W为船舶重力（kN），\alpha为撞击角度（°）。该公式考虑了船舶的重力、撞击角度、动能折减系数以及桥墩的弹性变形系数等因素来计算船撞力。挪威公共道路管理局规定：公式为P=3.5DWT^{\frac{1}{3}}，其中P为船撞力（kN），DWT为船舶载重吨位（t）。该公式仅考虑了船舶载重吨位对船撞力的影响，通过对船舶载重吨位开三次方并乘以一定系数来计算船撞力。北欧公共道路管理局规定：公式为P=0.5\sqrt{DWT}，其中P为船撞力（kN），DWT为船舶载重吨位（t）。该公式同样只考虑了船舶载重吨位对船撞力的影响，通过对船舶载重吨位开平方根并乘以一定系数来计算船撞力。2.2.2公式对比分析不同公式的计算结果存在较大差异。以一艘载重吨位为5000t，撞击速度为5m/s的船舶为例，使用Woisin公式计算得到的船撞力约为0.88\sqrt{5000}\pm50\approx149.2\pm50kN；使用AASHTO规范公式计算得到的船撞力约为0.122\sqrt{5\times5000}\pm50\approx42.7\pm50kN；使用IABSE规范公式计算得到的最大船撞力约为0.024\times5^{\frac{3}{2}}\times5000\approx424.3kN。各公式在不同工况下的适用性也有所不同。Woisin公式和AASHTO规范公式相对简单，计算方便，但它们主要考虑了船舶载重吨位和撞击速度的影响，对于船舶和桥梁的结构特性、碰撞角度等因素考虑较少，适用于对船撞力进行初步估算的情况。IABSE规范公式在一定程度上考虑了更多因素，计算结果相对更准确，但公式中的参数获取可能较为困难，适用于对计算精度要求较高的工程设计。Euro-code1规范公式根据船舶类型区分了船舶刚度系数，对于特定类型船舶的船撞力计算具有一定的针对性，但该公式对船舶速度的取值较为固定，限制了其在不同速度工况下的应用。我国公路桥涵设计通用规范公式和铁路桥涵设计基本规范公式分别针对漂流物和船舶撞击桥梁的情况，考虑了各自的特点和相关因素，但在应用时需要准确确定公式中的各项参数，否则会影响计算结果的准确性。挪威公共道路管理局规定和北欧公共道路管理局规定的公式过于简单，仅考虑了船舶载重吨位，在实际工程中应用时，计算结果的可靠性较低，一般仅作为参考。在实际工程应用中，应根据具体情况选择合适的公式。对于重要桥梁的设计和评估，应优先选择考虑因素较为全面、计算精度较高的公式，并结合数值仿真和试验研究等方法，对计算结果进行验证和修正，以确保桥梁的抗撞设计安全可靠。三、基于刚度的桥梁与船舶数值模型建立3.1桥梁结构数值模型3.1.1几何模型构建以某实际的多跨连续梁桥为案例，该桥梁位于繁忙的内河航道上，主桥跨径布置为（60+100+60）m，引桥采用30m跨径的简支梁桥，桥墩采用双柱式桥墩，桥台为重力式桥台。运用有限元软件ANSYS建立桥梁的几何模型，该软件以其强大的建模功能和丰富的单元库，能够准确地模拟桥梁结构的复杂几何形状和力学行为，在桥梁工程领域得到了广泛的应用。在构建桥墩几何模型时，首先根据设计图纸确定桥墩的截面尺寸和高度。对于双柱式桥墩，每根立柱的直径为1.5m，柱高根据不同的墩位在8-15m之间变化。利用ANSYS软件中的实体建模功能，通过创建圆柱体来模拟立柱，在圆柱体的顶部和底部，根据盖梁和基础的连接方式，进行适当的倒角和切角处理，以准确模拟实际的结构连接情况。在创建桥台几何模型时，重力式桥台的形状较为复杂，需要仔细根据图纸确定其各部分的尺寸和形状。通过组合使用ANSYS软件中的长方体、楔体等基本实体单元，逐步构建出台身、基础等部分的几何模型，并进行布尔运算，将各个部分组合成一个完整的桥台模型。主梁的几何模型构建则根据梁的截面形状和长度进行。该桥主梁采用预应力混凝土箱梁，箱梁的顶板宽度为12m，底板宽度为6m，梁高在跨中为2.5m，在墩顶处为4m，通过在ANSYS软件中定义箱梁的截面参数，利用扫掠功能沿着桥梁的纵向路径创建出主梁的三维实体模型，精确地模拟了主梁的空间几何形状。在建模过程中，严格按照设计图纸的尺寸进行输入，确保模型的几何准确性。对于一些复杂的细节部分，如桥墩与主梁的连接处、桥台的护坡等，通过适当的简化和等效处理，在不影响整体力学性能的前提下，提高建模效率。同时，对模型进行合理的网格划分，在关键部位如桥墩底部、主梁跨中以及桥跨连接处，采用较密的网格，以提高计算精度；在次要部位，采用相对较疏的网格，以减少计算量。3.1.2材料参数设定桥梁各部分材料的力学参数直接影响到数值模拟的准确性，其取值依据主要来源于材料的实际试验数据和相关的设计规范。对于桥墩和桥台，采用C40混凝土材料，根据《混凝土结构设计规范》（GB50010-2010），C40混凝土的弹性模量取值为3.25×10^4MPa，泊松比取0.2。弹性模量反映了材料抵抗弹性变形的能力，通过大量的混凝土材料试验，确定了C40混凝土在正常使用状态下的弹性模量取值，泊松比则是根据材料在单轴受力状态下横向应变与纵向应变的比值确定。密度取值为2500kg/m³，这是根据C40混凝土的配合比和组成成分，通过理论计算和实际测量得出的。主梁采用预应力混凝土材料，除了考虑混凝土的力学性能外，还需考虑预应力钢筋的作用。混凝土部分同样采用C50混凝土，其弹性模量取值为3.45×10^4MPa，泊松比为0.2，密度为2500kg/m³。预应力钢筋采用高强度低松弛钢绞线，其弹性模量为1.95×10^5MPa，泊松比取0.3，密度为7850kg/m³。这些参数的取值依据是钢绞线的产品标准和相关的力学性能测试报告，钢绞线具有高强度、低松弛的特点，其弹性模量和泊松比的取值反映了其在受力过程中的力学响应特性。在数值模型中，通过定义不同的材料属性，并将其赋予相应的几何模型部分，实现对桥梁结构材料特性的准确模拟。3.1.3边界条件处理在桥梁模型中，边界条件的设置至关重要，它直接影响到模型的计算结果和对实际结构的模拟准确性。桥墩底部的约束方式模拟实际的支撑情况，对于与基础刚性连接的桥墩底部，在ANSYS软件中采用固定约束，即限制桥墩底部在X、Y、Z三个方向的平动自由度和绕X、Y、Z轴的转动自由度，使桥墩底部在受力时不能发生任何位移和转动，准确地模拟了桥墩底部与基础之间的刚性连接。这种约束方式能够有效地传递桥墩所承受的荷载到基础，反映了实际工程中桥墩底部的支撑状态。对于桥梁的桥台，其与地基的接触采用弹性支撑模拟。根据地基的实际情况，通过查阅地质勘察报告，确定地基的基床系数。在ANSYS软件中，利用弹簧单元来模拟桥台与地基之间的弹性支撑，弹簧单元的刚度根据基床系数和接触面积计算得出，通过这种方式，考虑了地基在受力时的弹性变形，使模型更符合实际情况。在设置边界条件时，充分考虑了桥梁结构在实际使用过程中的受力特点和支撑方式，确保边界条件的合理性和准确性。对边界条件进行敏感性分析，通过改变边界条件的设置方式，观察模型计算结果的变化，进一步验证边界条件设置的正确性。3.2船舶数值模型3.2.1船舶外形简化在建立船舶数值模型时，为提高计算效率，同时保证模拟结果的准确性，需要对船舶外形进行合理简化。以一艘内河常见的散货船为例，该船船长为80m，型宽12m，型深5m。船舶的实际外形较为复杂，包含众多细节结构，如上层建筑、甲板设备等，这些细节结构在船桥碰撞过程中对整体力学性能的影响相对较小，若全部考虑会极大地增加计算量和计算时间。在简化过程中，主要保留船舶的主体结构，即船体的外壳和内部的主要支撑结构。将船舶的上层建筑简化为质量块，根据实际上层建筑的重量和重心位置，在模型中相应位置添加集中质量，以考虑其对船舶整体质量分布和惯性的影响。对于甲板上的各种设备，如起重机、通风管等，根据其重量和分布情况，进行适当的等效处理，将其重量均匀分布在甲板上，简化为均布质量。通过这样的简化方式，既减少了模型的复杂程度，提高了计算效率，又能保证船舶的主要力学性能，如质量、重心、惯性矩等在模型中得到准确体现，从而确保模拟结果的可靠性。3.2.2材料本构关系船舶材料主要为钢材，其本构关系考虑材料的非线性特性，如塑性、应变硬化等。钢材在受力过程中，当应力超过屈服强度后，会进入塑性变形阶段，此时材料的应力-应变关系不再遵循胡克定律，呈现出非线性的特征。采用vonMises屈服准则来描述钢材的屈服行为，该准则认为当材料的等效应力达到屈服强度时，材料开始屈服。在屈服后的塑性变形阶段，考虑应变硬化效应，即随着塑性应变的增加，材料的屈服强度会逐渐提高。通过定义合适的材料参数来描述这些特性，在ANSYS软件中，输入钢材的弹性模量、屈服强度、强化模量等参数。弹性模量为2.06×10^5MPa，屈服强度根据钢材的型号和标准取值，如Q345钢材的屈服强度为345MPa，强化模量则根据材料的试验数据或经验公式确定，一般取值在弹性模量的0.01-0.1倍之间。通过这些参数的设置，能够准确地模拟钢材在船桥碰撞过程中的力学响应，包括弹性变形、塑性变形以及应变硬化等行为。3.2.3附连水质量模拟附连水质量是指船舶在运动时带动周围水体一起运动所增加的质量，它对船舶撞击力有着重要影响。在数值模型中，采用莫里森方程来模拟附连水质量。莫里森方程将附连水质量的作用等效为附加在船舶结构上的惯性力，其表达式为F=\rho\pir^2C_ma，其中F为附连水质量产生的惯性力，\rho为水的密度，取1000kg/m³，r为船舶结构的特征半径，对于圆柱形结构，r为圆柱半径，对于其他形状的结构，可通过等效方法确定，C_m为附连水质量系数，一般取值在1.0-2.0之间，根据船舶的形状和运动状态通过经验公式或试验确定，a为船舶的加速度。在ANSYS软件中，通过在船舶结构的表面施加相应的惯性力来模拟附连水质量的影响。在船桥碰撞过程中，附连水质量会使船舶的质量增加，从而增大船舶的惯性，在碰撞瞬间，附连水质量产生的惯性力会对船舶的运动状态和撞击力产生显著影响。通过准确模拟附连水质量，能够更真实地反映船桥碰撞过程中的力学行为，提高数值模拟的准确性。3.3模型验证与校准3.3.1与已有实验数据对比将建立的数值模型结果与已有的船桥碰撞实验数据进行对比，以验证模型的准确性。已有实验数据来自某科研机构进行的船桥碰撞缩尺模型实验，实验中采用的船舶模型和桥梁模型与实际结构按照一定的相似比制作，通过在实验中测量碰撞力、结构变形等参数，得到了较为准确的实验数据。在对比过程中，将数值模型的计算结果与实验数据在相同工况下进行比较，包括船舶的撞击速度、撞击角度、船舶和桥梁的结构参数等。以碰撞力为例，数值模型计算得到的碰撞力时程曲线与实验测量的碰撞力时程曲线进行对比，观察两者在碰撞力峰值、作用时间以及变化趋势等方面的差异。通过对比发现，数值模型计算得到的碰撞力峰值与实验测量值相差在10%以内，碰撞力的作用时间和变化趋势也与实验数据基本吻合，表明数值模型能够较好地模拟船桥碰撞过程中的碰撞力变化情况。对于桥梁结构的变形，数值模型计算得到的桥墩和主梁的变形量与实验测量值进行对比，同样验证了数值模型在模拟结构变形方面的准确性。3.3.2参数敏感性分析通过改变模型中的关键参数，如船舶速度、质量、桥梁刚度等，分析这些参数对船撞力的敏感性。在改变船舶速度时，将船舶速度分别设置为3m/s、5m/s、7m/s，保持其他参数不变，计算不同速度下的船撞力。结果表明，船撞力随着船舶速度的增加而显著增大，当船舶速度从3m/s增加到5m/s时，船撞力峰值提高了约50%，这是因为船舶速度的增加导致其动能增大，在碰撞过程中传递给桥梁的能量也相应增加，从而使船撞力增大。改变船舶质量时，将船舶质量分别设置为1000t、1500t、2000t，计算不同质量下的船撞力。发现船撞力随着船舶质量的增加而增大，但增长幅度相对较小，当船舶质量从1000t增加到1500t时，船撞力峰值提高了约20%，这是因为船舶质量的增加主要影响其惯性，在碰撞过程中，质量较大的船舶更难改变运动状态，从而产生较大的撞击力，但相比速度的影响，质量对船撞力的影响相对较小。改变桥梁刚度时，通过调整桥墩和主梁的截面尺寸来改变桥梁的整体刚度，计算不同刚度下的船撞力。结果显示，桥梁刚度越大，船撞力越大，这是因为刚度较大的桥梁在受到撞击时，变形较小，能够承受更大的撞击力，而刚度较小的桥梁在撞击作用下容易发生较大的变形，从而吸收部分撞击能量，使船撞力相对减小。通过参数敏感性分析，明确了各关键参数对船撞力的影响程度，为进一步优化桥梁的抗撞设计和制定合理的船舶航行规则提供了依据。四、桥梁船撞力数值仿真分析4.1不同工况设置4.1.1船舶速度变化设定船舶撞击速度分别为3m/s、5m/s、7m/s和9m/s。在实际的内河航运中，船舶的航行速度受到多种因素的影响，如航道条件、船舶类型、交通规则等。根据相关统计数据，内河船舶在桥区水域的航行速度一般在3-9m/s之间，因此选取这几个速度值具有一定的代表性。通过数值仿真模拟不同速度下船舶与桥梁的碰撞过程，记录船撞力的大小和作用时间。速度的变化会导致船舶动能的改变，根据动能公式E_{k}=\frac{1}{2}mv^{2}（其中E_{k}为动能，m为船舶质量，v为船舶速度），速度的增加会使船舶动能呈平方倍增长，在碰撞过程中，这些动能将更多地转化为对桥梁的撞击能量，从而影响船撞力的大小和作用时间。4.1.2撞击角度差异考虑船舶以0°、15°、30°和45°的角度撞击桥梁。撞击角度是影响船撞力分布和桥梁受力状态的重要因素之一。0°撞击表示船舶正面撞击桥梁，此时船撞力主要集中在撞击点处，对桥墩的轴向压力较大；15°和30°的斜向撞击会使船撞力在桥墩上产生水平和竖向的分力，导致桥墩同时受到弯曲和剪切作用；45°的撞击角度下，船撞力的水平和竖向分力相对较为均衡，对桥墩的受力状态影响更为复杂。在实际的船桥碰撞事故中，由于船舶的航行姿态、航道条件以及驾驶员的操作失误等原因，船舶撞击桥梁的角度具有很大的随机性，因此研究不同撞击角度下的船撞力分布和桥梁受力状态具有重要的实际意义。4.1.3桥梁刚度改变通过调整桥墩和主梁的截面尺寸来改变桥梁的整体刚度。将桥墩的直径分别设置为1.5m、1.8m和2.0m，主梁的高度分别设置为2.5m、2.8m和3.0m。桥墩和主梁是桥梁结构的主要承重构件，它们的截面尺寸直接影响桥梁的刚度。增大桥墩直径和主梁高度可以提高桥梁的抗弯和抗剪能力，从而改变桥梁在船撞作用下的受力性能。在实际工程中，桥梁的刚度设计需要综合考虑多种因素，如桥梁的跨度、荷载等级、地质条件等。通过改变桥梁刚度进行数值仿真，可以研究桥梁刚度与船撞力之间的关系，为桥梁的抗撞设计提供参考依据。4.2仿真结果分析4.2.1船撞力时程曲线绘制不同工况下的船撞力时程曲线，以船舶速度为5m/s，撞击角度为0°为例，船撞力时程曲线如图1所示。从图中可以看出，船撞力在碰撞瞬间迅速上升，达到峰值后逐渐下降。船撞力峰值约为2500kN，作用时间约为0.2s。不同工况下的船撞力峰值和作用时间存在明显差异。随着船舶速度的增加，船撞力峰值显著增大，作用时间略有缩短。当船舶速度从3m/s增加到9m/s时，船撞力峰值从约1000kN增加到约4500kN，这是因为速度的增加导致船舶动能增大，在碰撞过程中传递给桥梁的能量增多，从而使船撞力峰值增大；而作用时间缩短则是因为更高的速度使得碰撞过程更加迅速。随着撞击角度的增大，船撞力峰值先减小后增大，作用时间逐渐延长。在0°撞击时，船撞力峰值最大，随着撞击角度增大到15°和30°，船撞力峰值有所减小，这是因为斜向撞击时船撞力在桥墩上产生分力，分散了部分撞击能量；当撞击角度增大到45°时，船撞力峰值又有所增大，这是因为此时船撞力的水平和竖向分力共同作用，对桥墩的受力更为复杂，导致船撞力峰值增大。作用时间延长是因为斜向撞击时船舶与桥墩的接触时间相对较长。改变桥梁刚度对船撞力峰值和作用时间也有影响，随着桥梁刚度的增大，船撞力峰值增大，作用时间缩短。这是因为刚度较大的桥梁在受到撞击时，变形较小，能够承受更大的撞击力，而刚度较小的桥梁在撞击作用下容易发生较大的变形，从而吸收部分撞击能量，使船撞力相对减小，同时变形过程也会延长碰撞作用时间。[此处插入船撞力时程曲线图片，图片标题为“船舶速度5m/s，撞击角度0°时船撞力时程曲线”]4.2.2能量转化分析研究船桥碰撞过程中的能量转化情况，在船桥碰撞过程中，船舶的动能主要转化为桥梁的变形能、船舶自身的变形能以及碰撞过程中的能量损耗，如声能、热能等。以船舶速度为7m/s，撞击角度为30°的工况为例，能量转化情况如图2所示。从图中可以看出，在碰撞初期，船舶动能迅速减小，同时桥梁变形能和船舶变形能迅速增加。随着碰撞的进行，船舶动能逐渐转化为桥梁变形能和船舶变形能，当碰撞达到一定程度后，能量转化趋于稳定。在整个碰撞过程中，船舶动能的大部分转化为桥梁变形能，约占70%左右，这表明桥梁在碰撞过程中主要通过自身的变形来吸收撞击能量；船舶变形能约占20%，碰撞过程中的能量损耗约占10%。不同工况下的能量转化比例也有所不同。随着船舶速度的增加，船舶动能增大，相应地，桥梁变形能和船舶变形能也会增大，但能量损耗的比例基本保持不变。随着撞击角度的增大，船舶动能在碰撞过程中的分配会发生变化，斜向撞击时，船舶动能更多地转化为桥梁的水平方向变形能，而竖向变形能相对减少。改变桥梁刚度会影响能量转化的速率和比例，刚度较大的桥梁在碰撞初期能够更快地吸收船舶动能，使船舶动能转化为桥梁变形能的速度加快，而刚度较小的桥梁则需要更长的时间来吸收能量，能量转化相对较为缓慢。[此处插入能量转化图图片，图片标题为“船舶速度7m/s，撞击角度30°时能量转化图”]4.2.3桥梁结构响应分析桥梁在船撞作用下的结构响应，如位移、应力分布等。以船舶速度为5m/s，撞击角度为15°的工况为例，桥墩在船撞作用下的位移云图和应力云图分别如图3和图4所示。从位移云图可以看出，桥墩在撞击点附近的位移最大，随着距离撞击点的增加，位移逐渐减小。桥墩顶部的位移约为0.15m，底部的位移相对较小，约为0.05m。从应力云图可以看出，桥墩在撞击点处的应力最大，超过了混凝土的抗压强度设计值，可能会导致桥墩局部混凝土出现开裂和破损。在桥墩与主梁的连接处，也存在较大的应力集中现象，需要加强该部位的构造措施。不同工况下的桥梁结构响应存在差异。随着船舶速度的增加，桥墩的位移和应力均显著增大，对桥梁结构的破坏作用更加明显。随着撞击角度的增大，桥墩的水平位移和弯曲应力增大，竖向位移和轴向应力相对减小。改变桥梁刚度会影响桥墩的位移和应力分布，刚度较大的桥墩在船撞作用下的位移和应力相对较小，能够更好地抵抗船撞作用；而刚度较小的桥墩则容易产生较大的位移和应力，结构的安全性较低。通过对桥梁结构响应的分析，可以评估桥梁在船撞作用下的损伤程度，为桥梁的修复和加固提供依据。[此处插入桥墩位移云图图片，图片标题为“船舶速度5m/s，撞击角度15°时桥墩位移云图”][此处插入桥墩应力云图图片，图片标题为“船舶速度5m/s，撞击角度15°时桥墩应力云图”]五、桥梁船撞力试验研究设计5.1试验目的与方案5.1.1试验目标确定本试验的主要目标是验证基于刚度的桥梁船撞力数值仿真结果的准确性和可靠性，深入研究桥梁船撞力的实际作用规律，为桥梁的抗撞设计和安全评估提供更为真实、可靠的依据。通过试验，能够直观地观察船舶撞击桥梁的过程，获取实际的碰撞数据，包括船撞力的大小、作用时间、作用位置，以及桥梁结构的变形、加速度等参数。将这些试验数据与数值仿真结果进行对比分析，可以检验数值模型的合理性和计算方法的正确性，发现数值仿真中可能存在的问题和不足之处，从而对数值模型和计算方法进行优化和改进，提高数值仿真的精度和可靠性。试验还有助于深入探究桥梁船撞力的实际作用规律。在实际的船桥碰撞中，存在着许多复杂的因素，如船舶和桥梁结构的非线性行为、材料的应变率效应、碰撞过程中的能量耗散等，这些因素在数值仿真中难以完全准确地模拟。通过试验研究，可以直接观察和分析这些复杂因素对船撞力和桥梁结构响应的影响，揭示船桥碰撞的内在力学机制，为进一步完善桥梁船撞力的理论和计算方法提供实验依据。通过试验研究还可以为桥梁的抗撞设计和安全评估提供更具实际意义的参考。在桥梁的设计阶段，准确掌握船撞力的实际作用规律和大小，能够合理确定桥梁的结构形式、材料选择和尺寸参数，提高桥梁的抗撞能力；在桥梁的运营阶段，通过试验研究建立的船撞力评估方法和指标，可以对桥梁的安全性进行实时监测和评估，及时发现潜在的安全隐患，采取相应的防护措施，保障桥梁的安全运营。5.1.2相似模型设计依据相似理论，设计与实际桥梁和船舶相似的试验模型。相似理论是模型试验的基础，它通过相似准则来保证模型与原型在力学行为上的相似性。在设计桥梁试验模型时，首先根据实际桥梁的结构形式和尺寸，按照一定的相似比进行缩尺。以某实际的多跨连续梁桥为例，该桥主桥跨径布置为（60+100+60）m，引桥采用30m跨径的简支梁桥，桥墩采用双柱式桥墩，桥台为重力式桥台。在试验模型中，将相似比确定为1:20，即模型的尺寸为实际桥梁尺寸的1/20。通过相似比的换算，确定模型中桥墩的直径、高度，主梁的截面尺寸和长度等参数。在确定桥墩直径时，实际桥墩直径为1.5m，按照相似比，模型桥墩直径应为1.5m÷20=0.075m。在确定材料参数时，考虑到模型与原型在材料性能上的相似性，模型材料的弹性模量、泊松比、密度等参数也应按照相似比进行换算。对于桥墩和主梁采用的混凝土材料，通过试验选择合适的模型材料，使其弹性模量和强度等性能与实际混凝土材料在相似比的基础上保持一致。在设计船舶试验模型时，同样根据实际船舶的外形和尺寸，按照相似比进行缩尺。以一艘内河常见的散货船为例，该船船长为80m，型宽12m，型深5m，将相似比确定为1:20后，模型船的船长为80m÷20=4m，型宽为12m÷20=0.6m，型深为5m÷20=0.25m。在确定船舶模型的材料参数时，考虑到船舶在碰撞过程中的力学行为，选择合适的材料来模拟船舶的外壳和内部结构，使其材料性能与实际船舶材料在相似比的基础上相匹配。在模型设计过程中，还需要考虑模型的加工精度和安装精度，确保模型的几何形状和尺寸符合设计要求，模型各部分之间的连接牢固可靠，以保证试验结果的准确性和可靠性。5.2试验设备与装置5.2.1撞击系统搭建用于模拟船舶撞击的试验设备为落锤装置，其工作原理是利用重物自由下落产生的动能来模拟船舶的撞击能量。落锤装置主要由落锤、导向轨道、提升系统和释放装置等部分组成。落锤采用钢材制成，其质量和形状根据试验需求进行设计，通过改变落锤的质量和下落高度，可以调节撞击能量的大小。导向轨道用于引导落锤的下落方向，确保落锤能够准确地撞击到桥梁模型上，导向轨道采用高精度的直线导轨，保证落锤在下落过程中的稳定性和准确性。提升系统用于将落锤提升到一定高度，通常采用电动绞车或液压提升装置，提升系统具有精确的高度控制功能，能够准确地设定落锤的下落高度。释放装置用于控制落锤的释放时机，保证落锤在达到预定高度后能够迅速、准确地释放，释放装置采用电磁控制或机械控制方式，具有较高的可靠性和精度。在本次试验中，落锤的质量设计为500kg，最大下落高度为5m，通过调整落锤的下落高度，可以模拟不同速度下的船舶撞击。当落锤从5m高度自由下落时，根据自由落体运动公式v=\sqrt{2gh}（其中v为落锤撞击速度，g为重力加速度，取9.8m/s²，h为下落高度），可计算出落锤撞击速度约为9.9m/s，能够较好地模拟内河船舶在桥区水域的高速撞击工况。5.2.2测量仪器选择测量船撞力使用高精度的力传感器，力传感器安装在桥梁模型的撞击点处，能够实时测量撞击过程中船撞力的大小和变化。本试验选用的力传感器量程为0-500kN，精度为0.1%FS（满量程），具有较高的灵敏度和准确性，能够满足试验对船撞力测量的要求。测量桥梁变形采用位移计，位移计布置在桥墩和主梁的关键部位，如桥墩顶部、底部和主梁跨中等处，通过测量这些部位在撞击过程中的位移变化，来分析桥梁结构的变形情况。选用的位移计量程为0-200mm，精度为0.01mm，能够精确地测量桥梁结构的微小变形。测量加速度使用加速度计，加速度计安装在桥梁模型的不同位置，用于测量桥梁在撞击过程中的加速度响应，分析桥梁结构的动力特性。选用的加速度计量程为0-100g（g为重力加速度），精度为0.1%FS，能够准确地测量桥梁在撞击瞬间的加速度变化。在试验过程中，力传感器、位移计和加速度计采集的数据通过数据采集系统进行实时采集和处理，数据采集系统具有高速、高精度的数据采集能力，能够同时采集多个传感器的数据，并对数据进行实时显示、存储和分析。通过对采集到的数据进行处理和分析，可以得到船撞力的时程曲线、桥梁结构的变形和加速度响应等信息，为后续的试验结果分析提供数据支持。5.3试验工况与步骤5.3.1工况设定设置与数值仿真相对应的试验工况，包括撞击速度、角度、桥梁刚度等条件。在撞击速度工况设置中，分别设定落锤下落高度为1m、2m、3m，对应模拟的船舶撞击速度约为4.4m/s、6.3m/s、7.7m/s，与数值仿真中设定的船舶速度3m/s、5m/s、7m/s相对应，通过不同的下落高度来模拟不同速度下的船舶撞击，研究撞击速度对船撞力和桥梁结构响应的影响。在撞击角度工况设置中，通过调整桥梁模型的摆放角度，使落锤分别以0°、15°、30°的角度撞击桥梁模型，与数值仿真中设定的撞击角度0°、15°、30°相对应，研究撞击角度对船撞力分布和桥梁受力状态的影响。在桥梁刚度工况设置中，通过改变桥墩和主梁的截面尺寸来调整桥梁模型的刚度，设置三组不同的刚度工况，分别为原模型刚度、增加桥墩直径20%后的刚度、增加主梁高度20%后的刚度，与数值仿真中通过调整桥墩和主梁截面尺寸改变桥梁刚度的工况相对应，研究桥梁刚度与船撞力之间的关系。5.3.2试验流程详细描述试验的具体操作步骤，从模型安装、仪器调试到撞击试验的实施过程。在模型安装阶段，首先将桥梁模型按照设计要求安装在试验台座上，确保桥梁模型的位置准确、固定牢固。使用高精度的测量仪器对桥梁模型的尺寸和位置进行测量和校准，保证模型的安装精度符合试验要求。将船舶模型（落锤）安装在落锤装置上，调整落锤的位置和姿态，使其能够准确地撞击到桥梁模型的预定位置。在仪器调试阶段，对力传感器、位移计、加速度计等测量仪器进行校准和调试，确保仪器的测量精度和性能满足试验要求。使用标准力源对力传感器进行校准，检查力传感器的灵敏度和线性度；使用位移校准装置对位移计进行校准，检查位移计的测量精度和重复性；使用加速度校准装置对加速度计进行校准，检查加速度计的测量精度和频率响应。将测量仪器连接到数据采集系统，进行数据采集和传输的测试，确保数据采集系统能够准确地采集和记录测量仪器的数据。在撞击试验实施阶段，首先根据试验工况设定落锤的下落高度和撞击角度，通过提升系统将落锤提升到预定高度。在落锤提升过程中，密切关注提升系统的运行状态和落锤的位置，确保提升过程安全可靠。当落锤达到预定高度后，通过释放装置迅速释放落锤，使落锤自由下落撞击桥梁模型。在落锤撞击桥梁模型的瞬间，测量仪器开始实时采集船撞力、桥梁变形和加速度等数据，数据采集系统对采集到的数据进行实时显示、存储和分析。一次撞击试验完成后，对桥梁模型和测量仪器进行检查，确认是否存在损坏或异常情况。如果桥梁模型和测量仪器正常，则根据试验计划进行下一次撞击试验；如果存在损坏或异常情况，则及时进行修复和调整，确保试验能够顺利进行。按照设定的试验工况，依次完成不同撞击速度、角度和桥梁刚度工况下的撞击试验，获取大量的试验数据。在试验过程中，对试验数据进行实时分析和处理，及时发现试验中存在的问题和异常情况，并采取相应的措施进行解决。试验结束后，对试验数据进行整理和分析，与数值仿真结果进行对比，验证数值仿真的准确性和可靠性，深入研究桥梁船撞力的实际作用规律。六、桥梁船撞力试验结果与数值仿真对比6.1试验结果分析6.1.1船撞力测量值分析在试验过程中，通过高精度的力传感器，对不同工况下的船撞力进行了精确测量。以船舶速度为6.3m/s（对应落锤下落高度2m），撞击角度为15°，桥梁刚度为原模型刚度的工况为例，测量得到的船撞力时程曲线呈现出典型的变化特征。在撞击瞬间，船撞力急剧上升，迅速达到峰值，随后随着时间的推移逐渐下降。该工况下的船撞力峰值为1200kN，作用时间约为0.18s。对不同工况下的船撞力测量值进行统计和分析，结果如表1所示。从表中可以看出，随着船舶速度的增加，船撞力峰值显著增大。当船舶速度从4.4m/s（对应落锤下落高度1m）增加到7.7m/s（对应落锤下落高度3m）时，船撞力峰值从750kN增加到1600kN，增长幅度超过100%。这是因为船舶速度的增加导致其动能增大，在碰撞过程中传递给桥梁的能量增多，从而使船撞力峰值增大。随着撞击角度的增大，船撞力峰值先减小后增大。在0°撞击时，船撞力峰值相对较大；当撞击角度增大到15°时，船撞力峰值有所减小；而当撞击角度进一步增大到30°时，船撞力峰值又有所增大。这是因为斜向撞击时船撞力在桥墩上产生分力，分散了部分撞击能量，但当撞击角度增大到一定程度时，分力的作用效果发生变化，导致船撞力峰值又有所上升。改变桥梁刚度对船撞力峰值也有明显影响，随着桥梁刚度的增加，船撞力峰值增大。当桥墩直径增加20%后，船撞力峰值相比原模型刚度工况下增大了约15%；当主梁高度增加20%后，船撞力峰值增大了约12%。这是因为刚度较大的桥梁在受到撞击时，变形较小，能够承受更大的撞击力。[此处插入船撞力测量值统计分析表，表标题为“不同工况下船撞力测量值统计分析表”]6.1.2桥梁损伤情况观察试验结束后，对桥梁模型的损伤情况进行了详细观察和记录。以船舶速度为7.7m/s，撞击角度为30°，桥梁刚度为原模型刚度的工况为例，桥墩在撞击点附近出现了明显的裂缝。裂缝宽度最大处达到了0.5mm，深度约为50mm，裂缝沿着桥墩的高度方向向上延伸。在桥墩与主梁的连接处，也出现了局部混凝土剥落的现象，剥落面积约为0.1m²。主梁在靠近撞击点的一侧，出现了轻微的弯曲变形，变形量约为15mm。通过对不同工况下桥梁模型损伤情况的观察和分析，发现损伤产生的原因主要是船撞力的作用。船舶的撞击使桥梁结构受到巨大的冲击力，当冲击力超过桥梁材料的极限强度时，就会导致结构出现裂缝、变形和混凝土剥落等损伤。损伤的规律与船撞力的大小和分布密切相关。船撞力越大，桥梁的损伤越严重；在撞击点附近，船撞力集中，损伤也最为明显。撞击角度的变化会影响船撞力在桥梁结构上的分布，从而导致不同的损伤模式。斜向撞击时，桥墩除了受到轴向压力外，还受到水平方向的分力，容易产生弯曲和剪切破坏，因此在桥墩与主梁的连接处和斜向受力部位，损伤较为突出。6.2与数值仿真结果对比6.2.1船撞力对比将试验测量得到的船撞力与数值仿真结果进行对比，以船舶速度为5m/s，撞击角度为15°的工况为例，对比曲线如图5所示。从图中可以看出，试验测量的船撞力峰值为1050kN，数值仿真计算得到的船撞力峰值为1100kN，两者相差约4.8%。在碰撞力的变化趋势上，试验结果和数值仿真结果基本一致，都呈现出在撞击瞬间迅速上升，达到峰值后逐渐下降的特点。对不同工况下的船撞力试验值和仿真值进行对比分析，结果如表2所示。从表中可以看出，在大多数工况下，试验值与仿真值的相对误差在10%以内，说明数值仿真结果与试验测量结果具有较好的一致性。但在个别工况下，如船舶速度为7m/s，撞击角度为30°时，相对误差达到了12%，这可能是由于试验过程中的一些不确定性因素，如模型制作误差、测量误差以及实际碰撞过程中的复杂因素等导致的。[此处插入船撞力试验值与仿真值对比曲线图片，图片标题为“船舶速度5m/s，撞击角度15°时船撞力试验值与仿真值对比曲线”][此处插入船撞力试验值与仿真值对比分析表，表标题为“不同工况下船撞力试验值与仿真值对比分析表”]6.2.2桥梁响应对比对比试验和仿真中桥梁的位移、应力等响应数据，以船舶速度为5m/s，撞击角度为0°的工况为例，桥墩顶部的位移试验测量值为0.12m，数值仿真计算值为0.13m，相对误差约为8.3%。在应力方面，桥墩底部的最大应力试验测量值为15MPa，数值仿真计算值为16MPa，相对误差约为6.7%。从整体上看，试验和仿真中桥梁的位移和应力响应趋势基本一致，数值仿真能够较好地反映桥梁在船撞作用下的响应情况。但在一些细节方面，如局部应力集中区域的应力分布，试验结果和仿真结果存在一定的差异。这可能是因为数值模型在简化过程中，对一些复杂的结构细节进行了处理，导致与实际情况存在一定的偏差。通过对不同工况下桥梁响应的对比分析，进一步评估了数值模型的准确性，为数值模型的改进和完善提供了依据。6.3误差分析与原因探讨6.3.1误差计算与分析通过计算试验结果与数值仿真结果之间的误差，深入分析误差的大小和分布情况。误差计算公式为：\text{相对误差}=\frac{\vert\text{试验值}-\text{仿真值}\vert}{\text{试验值}}\times100\%。以船撞力为例，在不同工况下的相对误差计算结果如表3所示。从表中可以看出，船撞力的相对误差在2%-12%之间，平均相对误差为7.5%。在位移和应力响应方面，相对误差也呈现出类似的分布情况。在速度较低、撞击角度较小的工况下，误差相对较小；而在速度较高、撞击角度较大的工况下，误差相对较大。这表明在复杂工况下，数值仿真结果与试验结果的偏差可能会增大，需要更加关注。[此处插入试验结果与数值仿真结果误差分析表，表标题为“试验结果与数值仿真结果误差分析表”]6.3.2误差产生原因探讨误差产生的原因是多方面的。首先，试验模型与实际结构存在差异。试验模型是按照相似理论进行缩尺制作的，虽然在几何形状和材料性能上尽量与实际结构相似，但仍然无法完全消除差异。在材料性能方面，模型材料的实际性能可能与理论值存在一定偏差；在结构连接方式上，模型的连接方式可能无法完全模拟实际结构的复杂连接情况。测量误差也是导致误差产生的重要原因之一。力传感器、位移计和加速度计等测量仪器在测量过程中存在一定的精度限制，测量环境的干扰也可能影响测量结果的准确性。数值模型的简化也会带来误差。在建立数值模型时，为了提高计算效率，对船舶和桥梁结构进行了一定的简化，如对一些复杂的结构细节进行了忽略，对材料的非线性行为进行了近似处理等。这些简化措施虽然在一定程度上提高了计算效率，但也导致数值模型与实际情况存在一定的偏差。在模拟船桥碰撞过程中，数值模型对碰撞接触算法、材料本构关系等的选择和设定也会影响计算结果的准确性。如果这些参数选择不当，就会导致数值仿真结果与试验结果之间出现较大的误差。七、基于研究结果的桥梁抗撞设计建议7.1桥梁刚度优化策略7.1.1合理刚度取值范围根据数值仿真和试验结果，在不同通航条件下，桥梁的合理刚度取值范围存在差异。在通航船舶吨位较小、航行速度较低的内河小型航道，船舶撞击能量相对较小，桥梁刚度可适当降低。当船舶载重吨位在500t以下，航行速度在3m/s左右时，桥梁桥墩的纵向抗弯刚度EI可控制在1×10^7-3×10^7kN・m²之间，横向抗弯刚度EI可控制在0.8×10^7-2×10^7kN・m²之间。这样的刚度取值既能满足桥梁在正常使用状态下的受力要求，又能在一定程度上吸收船舶撞击能量，减少对桥梁结构的损伤。在通航船舶吨位较大、航行速度较高的内河大型航道或沿海航道，船舶撞击能量较大，需要提高桥梁的刚度以增强其抗撞能力。当船舶载重吨位在1000-5000t，航行速度在5-7m/s时，桥梁桥墩的纵向抗弯刚度EI宜提高到5×10^7-1×10^8kN・m²之间，横向抗弯刚度EI宜提高到4×10^7-8×10^7kN・m²之间。通过增加桥墩的截面尺寸、采用高强度材料等方式提高桥梁刚度，可以有效减小船舶撞击时桥梁的变形和应力，保障桥梁的安全。在确定桥梁合理刚度取值范围时，还需考虑桥梁的结构形式、跨度、地质条件等因素的影响，综合进行分析和判断。7.1.2结构形式调整建议对于梁桥，可通过增加桥墩数量来减小单墩所承受的船撞力。在跨度较大的梁桥中，适当加密桥墩布置，将较大的跨度划分为多个较小的跨度，使船舶撞击力能够分散到多个桥墩上，降低单个桥墩的受力风险。优化主梁的截面形状，采用箱型截面或T型截面等抗撞性能较好的截面形式，增加主梁的抗弯和抗剪能力，提高梁桥的整体抗撞性能。对于拱桥，加强拱脚的连接构造，采用更牢固的锚固方式和增强的混凝土强度，提高拱脚在船撞作用下的承载能力。合理设计拱圈的矢跨比，在满足通航和桥梁使用功能的前提下，选择合适的矢跨比，使拱桥在船撞作用下的受力更加合理，减小拱圈的应力集中。对于斜拉桥和悬索桥，加强索塔的刚度和强度，增加索塔的截面尺寸、采用高强度钢材或增加配筋等措施，提高索塔抵抗船撞力的能力。优化拉索或悬索的布置方式，合理调整索力分布，使桥梁在船撞作用下能够更好地传递和分散荷载，减小结构的变形和应力。7.2防撞设施设计要点7.2.1现有防撞设施分析常见的桥梁防撞设施包括防撞墩、防撞浮筒等，它们各自具有不同的工作原理和效果。防撞墩通常采用钢筋混凝土或钢结构制成，通过自身的强度和刚度来抵抗船舶的撞击。其工作原理是利用防撞墩与桥墩之间的连接结构，将船舶撞击力传递到桥墩上，同时通过防撞墩自身的变形来吸收部分撞击能量。防撞墩能够有效地减小船舶撞击对桥墩的直接冲击力，降低桥墩受损的风险。在一些内河桥梁中，采用钢筋混凝土防撞墩，当船舶撞击时，防撞墩能够承受较大的撞击力，保护桥墩免受严重破坏。但防撞墩的缺点是对船舶的保护作用相对较小，在撞击过程中可能会对船舶造成一定的损伤，而且防撞墩的安装和维护成本较高