基于数值模拟探究木拱廊桥抗洪性能：结构、载荷与防护策略

基于数值模拟探究木拱廊桥抗洪性能：结构、载荷与防护策略一、绪论1.1研究背景木拱廊桥作为中国传统木构桥梁中技术含量最高的品类，不仅是交通设施，更是承载着深厚历史文化价值的瑰宝。它凝聚着古代劳动人民的智慧与创造力，是中国桥梁建筑史上的杰出代表，反映了特定历史时期的社会、经济、文化等多方面的状况。在漫长的岁月中，木拱廊桥横跨山川溪流，见证了无数的历史变迁，成为地域文化的重要象征。2009年，“中国木拱桥传统营造技艺”被列入联合国教科文组织《急需保护的非物质文化遗产名录》，2012年“闽浙木拱廊桥”正式列入中国世界文化遗产预备名单，这一系列荣誉足以彰显其在世界文化遗产中的重要地位。然而，木拱廊桥大多坐落于山区河流之上，这些地区每逢雨季或台风季节，极易遭受洪水侵袭。洪水的强大冲击力对木拱廊桥的结构安全构成了严重威胁，众多木拱廊桥在洪灾中受损甚至被冲毁。如2016年9月15日，超强台风“莫兰蒂”带来的暴雨引发洪水，致使浙江温州泰顺县的文兴桥、薛宅桥和文重桥等国家级重点保护古廊桥被冲毁。这一事件引起了社会各界对木拱廊桥命运的广泛关注与痛心。据不完全统计，近年来因洪水灾害而遭受不同程度破坏的木拱廊桥已达数十座，且随着全球气候变化，极端天气事件增多，洪水发生的频率和强度呈上升趋势，木拱廊桥面临的抗洪形势愈发严峻。倘若不能深入了解木拱廊桥的抗洪性能，采取有效的保护措施，这些珍贵的历史文化遗产很可能在未来的洪灾中遭受更大的损失，甚至永远消失。因此，开展木拱廊桥抗洪性能研究迫在眉睫，具有极其重要的现实意义。1.2研究目的与意义本研究旨在借助数值模拟技术，深入剖析木拱廊桥在洪水作用下的力学响应和结构性能变化，全面揭示其抗洪机理，为木拱廊桥的保护和修复提供科学的理论依据与技术支持。通过建立精准的数值模型，模拟不同洪水工况下木拱廊桥的受力情况，包括结构内力分布、变形模式以及关键节点的应力应变状态等，明确木拱廊桥在洪水中的薄弱环节和潜在风险点。同时，基于模拟结果，评估现有木拱廊桥的抗洪能力，为制定针对性的保护措施和修复方案提供量化的数据参考，提高保护工作的科学性和有效性。木拱廊桥作为珍贵的历史文化遗产，对其开展抗洪性能研究具有多方面的重要意义。从文化传承角度来看，木拱廊桥承载着深厚的历史文化内涵，是古代桥梁建造技艺的杰出代表，保护好木拱廊桥就是守护人类文明的瑰宝，传承中华民族的历史记忆和文化基因。对木拱廊桥抗洪性能的研究，有助于我们更好地了解古人在应对自然灾害时的智慧和创造力，为文化遗产的保护和传承提供理论支撑，让这一独特的文化遗产得以延续和发展，使其在现代社会中继续发挥文化纽带和精神寄托的作用。从工程技术角度而言，木拱廊桥独特的结构体系蕴含着丰富的力学原理和建造技术，通过对其抗洪性能的研究，可以深入挖掘其中的科学价值，为现代桥梁工程的设计和建造提供有益的借鉴。研究木拱廊桥在洪水作用下的力学行为和结构响应，有助于完善桥梁结构的抗灾理论和设计方法，提高现代桥梁在极端水文条件下的安全性和可靠性。数值模拟技术在木拱廊桥抗洪性能研究中的应用，不仅能够弥补传统试验方法的局限性，降低研究成本和时间，还可以实现对复杂洪水工况和结构响应的精细化模拟，为木拱廊桥的保护和工程技术的发展提供了一种高效、准确的研究手段，推动相关领域的技术进步和创新。1.3国内外研究现状1.3.1木拱廊桥结构研究现状木拱廊桥作为一种独特的传统桥梁形式，其结构研究一直是国内外学者关注的焦点。在国内，众多学者从不同角度对木拱廊桥的结构展开了深入探究。例如，有学者对浙江南部山区木拱廊桥的结构设计进行研究，以庆元举水如龙桥为例，详细分析了木拱廊桥的基本结构，指出桥木拱架部分由上下两层系统组成，下层为“三节苗”，上层为“五节苗”，二者相互穿插，并通过“牛头”使拱架相互联系成为一整体，顶部的水平拱木“平苗”在结构上相当于梁的作用。还有学者通过对泰顺木拱廊桥的研究，发现其桥梁由二组拱骨相贯而成，为避免桥柱受洪水冲击，加强梁架的横向稳定性，还设有横向斜撑架，即“剪刀架”，节点横梁起着横向力的均衡分配作用，整座桥力的分布非常合理。这些研究成果深入揭示了木拱廊桥的结构特点和力学原理，为后续的研究奠定了坚实基础。在国外，虽然木拱廊桥主要分布于中国，但随着对中国传统建筑文化研究的深入，一些国外学者也开始关注木拱廊桥的结构。他们运用先进的数字化技术和力学分析方法，对木拱廊桥的结构进行模拟和分析，从国际视角为木拱廊桥结构研究提供了新的思路和方法，促进了木拱廊桥结构研究的国际化交流与合作。1.3.2木拱廊桥抗洪性能研究现状木拱廊桥的抗洪性能研究对于其保护和维护具有重要意义。国内许多学者针对木拱廊桥在洪水作用下的响应和破坏模式开展了研究。有研究通过对历史上木拱廊桥水毁案例的分析，总结出洪水对木拱廊桥的破坏形式主要包括桥墩冲毁、桥身倾斜、构件脱落等，并探讨了水毁的原因，如洪水流量过大、桥梁基础不稳、结构抗冲击能力不足等。还有学者通过实地调研和模型试验，研究了木拱廊桥在不同洪水工况下的受力情况和变形特征，提出了一些提高木拱廊桥抗洪能力的措施，如加强桥墩基础、优化结构连接方式等。然而，目前对于木拱廊桥抗洪性能的研究仍存在一定局限性，缺乏系统全面的理论体系和量化分析方法。国外在桥梁抗洪研究方面具有丰富的经验和先进的技术手段，虽然针对木拱廊桥抗洪性能的专门研究较少，但在桥梁水动力学、结构抗冲击等相关领域的研究成果，为木拱廊桥抗洪性能研究提供了有益的借鉴。例如，国外在桥梁结构水动力荷载计算、洪水与桥梁相互作用的数值模拟等方面的研究，为深入研究木拱廊桥在洪水中的受力和变形提供了理论和技术支持。1.3.3数值模拟在桥梁领域应用现状随着计算机技术的飞速发展，数值模拟在桥梁领域的应用日益广泛。在国内，数值模拟技术已被广泛应用于现代桥梁的设计、分析和评估中。例如，在桥梁结构的力学分析中，通过有限元软件对桥梁结构进行建模，模拟不同荷载工况下桥梁的应力、应变分布，评估桥梁的承载能力和安全性。在桥梁抗风、抗震研究中，利用数值模拟技术可以模拟风场、地震波作用下桥梁的动力响应，为桥梁的抗风、抗震设计提供依据。在桥梁施工过程模拟中，通过数值模拟可以预测施工过程中桥梁结构的变形和内力变化，指导施工方案的制定和优化。同时，数值模拟技术也逐渐应用于古桥保护研究中，为古桥的结构分析和保护修复提供了新的手段。在国外，数值模拟技术在桥梁领域的应用更为成熟和深入。国外学者利用先进的数值模拟软件和高性能计算机，开展了大量关于桥梁结构性能、灾害响应等方面的研究。例如，在桥梁抗洪研究中，运用计算流体动力学（CFD）方法模拟洪水对桥梁的作用，分析桥梁周围的水流场分布和水动力荷载，评估桥梁在洪水中的稳定性。在桥梁疲劳寿命预测、材料性能研究等方面，数值模拟技术也发挥了重要作用，为桥梁工程的发展提供了强大的技术支持。1.4研究内容与方法1.4.1研究内容本研究聚焦于木拱廊桥抗洪性能，借助数值模拟技术展开深入探究，主要涵盖以下几个关键方面：木拱廊桥所在流域的水文计算：全面收集木拱廊桥所在流域的地形、气象、水文等基础数据，运用水文分析方法，精准确定洪水重现期。基于此，详细计算设计暴雨，包括平均点雨量、设计面雨量以及设计暴雨雨型。进而，通过合理的计算方法，准确推算设计洪峰流量和洪水过程线，为后续研究提供可靠的洪水工况数据，明确木拱廊桥在不同洪水条件下可能面临的水文环境。木拱廊桥所受洪水作用的仿真计算：采用计算流体动力学（CFD）方法，借助专业的Fluent软件，构建木拱廊桥与洪水相互作用的数值模型。根据木拱廊桥的实际结构特点和尺寸，对模型进行精细设置，包括附加模型的选定、物理求解方程的确定以及相关参数的设定。通过模拟不同洪水工况下木拱廊桥周围的水流场分布，深入分析木拱廊桥所受到的水动力荷载，如压力、流速、冲击力等，明确洪水对木拱廊桥结构的作用机制和影响程度。洪水作用下木拱廊桥的有限元结构分析：运用有限元分析方法，借助SAP2000软件，建立木拱廊桥的三维有限元力学模型。合理选定单元对象，准确模拟木拱廊桥各构件的力学行为，如梁、柱、拱等。将CFD仿真计算得到的水动力荷载施加到有限元模型上，模拟木拱廊桥在洪水作用下的结构响应，包括应力、应变分布以及变形情况等。通过分析计算结果，确定木拱廊桥在洪水中的薄弱部位和潜在风险点，评估其结构的安全性和稳定性。木拱廊桥抗洪方法研究：基于上述研究结果，深入探讨提高木拱廊桥抗洪能力的有效方法。一方面，研究拆除风雨板对木拱廊桥抗洪能力的影响，通过对比拆除前后模型的结构受力情况和抗洪性能指标，分析拆除风雨板在减轻洪水冲击力、改善结构受力状态方面的作用机制。另一方面，研究桥面加重对木拱廊桥抗洪能力的影响，选取不同的桥面加重工况，模拟加重过程中木拱廊桥的结构响应，确定安全状态的判断标准，评估桥面加重在增强结构稳定性、提高抗洪能力方面的效果。此外，还将结合实际工程经验和研究成果，提出其他可行的木拱廊桥抗洪思路和措施，为木拱廊桥的保护和修复提供技术支持。1.4.2研究方法本研究综合运用多种研究方法，以确保研究的全面性、科学性和准确性，具体包括：文献研究法：广泛查阅国内外关于木拱廊桥结构、抗洪性能以及数值模拟在桥梁领域应用等方面的文献资料，全面了解相关研究现状和发展趋势，总结前人的研究成果和经验，为本次研究提供坚实的理论基础和研究思路。通过对文献的梳理和分析，明确木拱廊桥的结构特点、受力机制以及在洪水中的破坏模式等关键信息，为后续的数值模拟和分析提供参考依据。数值模拟法：这是本研究的核心方法。利用CFD软件对木拱廊桥所受洪水作用进行仿真计算，模拟洪水在木拱廊桥周围的流动特性和水动力荷载分布；运用有限元软件对木拱廊桥的结构进行分析，模拟其在洪水作用下的力学响应和变形情况。通过数值模拟，可以实现对不同洪水工况和结构参数下木拱廊桥抗洪性能的精细化研究，弥补传统试验方法在成本、时间和工况模拟等方面的局限性，深入揭示木拱廊桥的抗洪机理和结构性能变化规律。案例分析法：选取具有代表性的木拱廊桥作为研究案例，如泰顺的文兴桥、薛宅桥等，收集这些桥梁在历史洪灾中的受损情况和相关数据资料。结合数值模拟结果，对案例进行深入分析，验证数值模拟的准确性和可靠性，同时进一步了解木拱廊桥在实际洪灾中的破坏过程和原因，为提出针对性的抗洪保护措施提供实践依据。通过案例分析，还可以总结不同类型木拱廊桥的抗洪特点和共性问题，为其他木拱廊桥的保护和修复提供借鉴。1.5技术路线本研究的技术路线旨在系统、全面地开展基于数值模拟的木拱廊桥抗洪性能研究，其核心流程如下：首先，通过广泛查阅国内外相关文献，收集木拱廊桥所在流域的地形、气象、水文等基础数据，为后续研究奠定坚实的数据基础。在水文计算环节，运用专业的水文分析方法，确定洪水重现期，进而计算设计暴雨和设计洪水，包括平均点雨量、设计面雨量、设计暴雨雨型、设计洪峰流量以及洪水过程线等，明确木拱廊桥面临的洪水工况。随后，借助CFD方法，利用Fluent软件构建木拱廊桥与洪水相互作用的数值模型，模拟不同洪水工况下木拱廊桥周围的水流场分布，分析其所受到的水动力荷载。同时，运用有限元分析方法，通过SAP2000软件建立木拱廊桥的三维有限元力学模型，将CFD仿真得到的水动力荷载施加到该模型上，模拟木拱廊桥在洪水作用下的结构响应，确定其薄弱部位和潜在风险点。最后，基于上述研究结果，深入探讨提高木拱廊桥抗洪能力的方法，如研究拆除风雨板和桥面加重对木拱廊桥抗洪能力的影响，并提出其他可行的抗洪思路和措施。整个技术路线环环相扣，从资料收集、数据计算、模型模拟到结果分析与应用，形成一个完整、科学的研究体系，以实现对木拱廊桥抗洪性能的深入探究，为其保护和修复提供有力的技术支持，具体技术路线图见图1-1。图1-1技术路线图二、木拱廊桥结构特点与抗洪性能关系2.1木拱廊桥结构类型与特点2.1.1结构类型木拱廊桥作为中国传统桥梁中的独特类型，其结构类型丰富多样，主要包括贯木拱廊桥和八字撑木拱廊桥等，每种类型都蕴含着古人的智慧与独特的力学原理。贯木拱廊桥是木拱廊桥中最为典型的一种，其结构精巧复杂，具有极高的技术含量。它以梁木穿插别压形成拱桥，足支撑在两岸的岩石上，底座由数十根粗大圆木纵横拼接对拱而成“八字结构”。这种结构不用钉铆，完全依靠自身的强度、摩擦力以及构件直径的大小、所成的角度和水平距离等因素巧妙搭接。其拱架部分通常由上下两层系统组成，下层为“三节苗”，由3根圆木纵向连成八字形拱架，共9组；上层为“五节苗”，采用5根稍短的圆木纵连成五折边形拱架，共8组，与“三节苗”相互穿插，并通过“牛头”使拱架相互联系成为一个整体。顶部的水平拱木“平苗”，与“牛头”以燕尾榫卯接，在结构上相当于梁的作用，承担着桥面传来的荷载，并将其传递至拱架。这种独特的结构形式使得贯木拱廊桥能够以较短的构造材料实现较大的跨度，在世界桥梁史上独树一帜，如福建寿宁、屏南和浙江泰顺、庆元、景宁等地的许多木拱廊桥都属于贯木拱廊桥类型，它们不仅是交通设施，更是地域文化的象征，承载着丰富的历史文化内涵。八字撑木拱廊桥则具有另一种独特的结构特点。其在拱架的基础上，增设了八字形的斜撑结构，这些斜撑从桥墩或桥台向上斜向支撑拱架，与拱架形成稳定的三角形结构。八字撑的设置增强了拱架的稳定性，使其能够更好地承受竖向荷载和水平力的作用。在一些八字撑木拱廊桥中，斜撑与拱架之间通过榫卯连接，节点构造精巧，既能有效地传递力，又保证了结构的整体性。这种结构类型的木拱廊桥在抵御洪水冲击时，八字撑能够分散洪水对拱架的冲击力，减少拱架的变形和损坏风险。与贯木拱廊桥相比，八字撑木拱廊桥的结构更加稳固，在应对复杂的水流条件和较大的洪水冲击力时具有一定的优势，但其建造工艺也更为复杂，对工匠的技术要求更高，如在一些山区的河流上，八字撑木拱廊桥以其独特的结构屹立于洪水之中，见证了岁月的变迁。2.1.2典型结构剖析以泰顺县的文兴桥为例，对木拱廊桥的典型结构进行深入剖析，有助于更全面地了解木拱廊桥的结构特点及其抗洪性能。文兴桥是一座具有较高历史文化价值和独特结构特征的木拱廊桥，其建造工艺精湛，结构设计巧妙，在木拱廊桥中具有一定的代表性。文兴桥的拱架是其关键的受力结构，采用了贯木拱的形式。由“三节苗”和“五节苗”相互穿插组成，通过“牛头”连接成稳固的整体。“三节苗”和“五节苗”的木材选用粗大的圆木，直径较大，以保证拱架具有足够的强度和刚度来承受桥梁的自重和各种荷载。在实际受力过程中，拱架将桥面传来的竖向荷载转化为轴向压力，并通过拱的曲线形状将力传递至桥台和桥墩。当受到洪水冲击时，拱架能够凭借其独特的结构形式，将洪水的冲击力分散到各个构件上，从而减小单个构件所承受的压力。“牛头”作为连接“三节苗”和“五节苗”的关键节点，承受着较大的弯矩和剪力，其节点构造采用了榫卯连接，这种连接方式不仅具有良好的传力性能，还能在一定程度上适应构件的变形，增强了拱架的整体性和稳定性。在历史上的洪灾中，文兴桥的拱架多次经受住了洪水的考验，尽管部分构件有所损坏，但整体结构依然保持相对稳定，展现了其卓越的抗洪能力。廊屋是文兴桥的重要组成部分，不仅为行人提供了遮风挡雨的场所，还对桥梁的结构稳定性起到了辅助作用。廊屋采用穿斗式木构架，这种构架形式具有结构轻巧、整体性好的特点。其由多根立柱和横梁相互穿插连接而成，立柱垂直于地面，承担着上部结构的重量，并将力传递至基础。横梁则横向布置，与立柱连接，形成稳定的框架结构。廊屋的屋面采用双坡顶或悬山顶，有利于排水，避免雨水积聚对桥梁结构造成损害。在抗洪方面，廊屋能够增加桥梁的自重，使桥梁在洪水冲击下更加稳定。同时，廊屋的存在还能在一定程度上阻挡洪水对拱架和桥面的直接冲击，起到缓冲作用。然而，廊屋也会增加桥梁所受的风荷载和水压力，在设计和建造时需要综合考虑这些因素，以确保廊屋与拱架等结构的协同工作，共同抵御洪水的侵袭。桥台是文兴桥与河岸连接的重要结构，它承担着拱架传递过来的巨大压力，并将其分散到地基中。文兴桥的桥台通常采用石材砌筑，石材具有较高的抗压强度和耐久性，能够保证桥台在长期的使用过程中稳定可靠。桥台的基础深入地下，与地基紧密结合，以增加桥台的稳定性。在洪水作用下，桥台除了要承受拱架传来的竖向力和水平力外，还要抵御洪水的冲刷和浸泡。为了增强桥台的抗洪能力，其表面通常会进行加固处理，如采用水泥砂浆勾缝，防止洪水渗入石材缝隙，导致桥台结构松动。桥台的形状和尺寸也经过精心设计，以减小洪水对桥台的冲击力。例如，桥台的迎水面通常设计成流线型，使洪水能够顺畅地流过，减少局部冲刷。在一些洪水灾害中，文兴桥的桥台虽然受到了洪水的强烈冲击，但由于其坚固的结构和合理的设计，依然能够保持稳定，为桥梁的整体安全提供了有力保障。2.2结构特点对抗洪性能的影响2.2.1拱架结构的抗洪优势木拱廊桥的拱架结构是其抗洪的关键所在，展现出卓越的力学性能和独特的优势，使其在洪水冲击下能够保持相对稳定。以贯木拱廊桥的拱架为例，其由“三节苗”和“五节苗”相互穿插组成，这种结构形式巧妙地利用了木材之间的相互作用和力学原理，实现了力的有效分散和传递。在洪水来临时，强大的冲击力会作用于拱架上，而拱架能够凭借其独特的结构，将冲击力沿着拱的曲线方向分散到各个构件上。“三节苗”和“五节苗”的交叉布置使得拱架形成了多个三角形稳定结构，三角形具有稳定性强的特点，能够有效地抵抗变形和外力。当洪水的冲击力作用于拱架时，这些三角形结构能够协同工作，将力分散到整个拱架体系中，避免了单个构件承受过大的荷载而导致破坏。拱架的结构形式还使得其在承受竖向荷载时，能够将荷载转化为轴向压力，通过拱的推力传递至桥台和桥墩，从而减小了对拱架自身的弯矩和剪力作用。这种力的转化和传递方式使得拱架在洪水作用下能够保持较好的稳定性，不易发生弯曲和断裂等破坏现象。通过数值模拟分析可以进一步验证拱架结构的抗洪优势。在模拟不同洪水工况下，当洪水冲击力作用于拱架时，拱架各构件的应力分布相对均匀，没有出现局部应力集中的现象。在高水位、大流速的洪水模拟中，拱架的最大应力值远低于木材的极限强度，表明拱架结构能够承受较大的洪水冲击力。与其他结构形式的桥梁相比，如梁式桥，在相同的洪水条件下，梁式桥的梁体往往会承受较大的弯矩和剪力，容易在跨中或支座处出现破坏，而木拱廊桥的拱架结构则能够有效地分散这些力，提高了桥梁的抗洪能力。在一些实际的洪水灾害中，木拱廊桥的拱架结构经受住了考验，尽管部分附属构件受到损坏，但拱架主体依然能够保持相对完整，为桥梁的修复和后续使用提供了可能。2.2.2廊屋对桥梁抗洪的作用廊屋作为木拱廊桥的重要组成部分，在抗洪过程中发挥着多重作用，对桥梁的抗洪性能产生着重要影响。廊屋的重量能够增加桥梁的整体稳定性，使其在洪水冲击下更不易被冲走。廊屋采用的穿斗式木构架，由多根立柱和横梁组成，这些构件的重量分布在桥面上，增加了桥梁与基础之间的摩擦力。当洪水的冲击力试图推动桥梁时，廊屋的重量所产生的摩擦力能够有效地抵抗这种推力，使桥梁保持在原位。在一些洪水灾害中，观察到廊屋保存相对完好的木拱廊桥，其整体结构的稳定性明显优于廊屋受损严重的桥梁，这充分说明了廊屋重量对桥梁抗洪稳定性的积极作用。廊屋的防风挡水设计也对桥梁的抗洪性能起到了重要的保护作用。廊屋的屋面通常采用双坡顶或悬山顶，这种设计有利于雨水的快速排除，避免雨水在桥面上积聚。雨水积聚可能会增加桥面的荷载，导致桥梁结构承受额外的压力，甚至可能引发桥面的变形和损坏。而合理的屋面排水设计能够及时将雨水排离桥面，减轻了桥梁的负担。廊屋的侧面通常设置有木板或栏杆等防风挡水设施，这些设施能够在一定程度上阻挡洪水对拱架和桥面的直接冲击。当洪水携带杂物和水流冲击桥梁时，防风挡水设施能够起到缓冲作用，减少了洪水对桥梁结构的冲击力。在洪水模拟实验中，对比设置了防风挡水设施和未设置的桥梁模型，发现设置了防风挡水设施的模型在洪水冲击下，其拱架和桥面的受损程度明显减轻，这表明廊屋的防风挡水设计能够有效地保护桥梁结构，提高其抗洪能力。然而，廊屋在抗洪过程中也存在一些潜在的问题。由于廊屋增加了桥梁的迎风面积，在遭遇强风时，可能会受到较大的风荷载作用。如果风荷载过大，可能会导致廊屋结构的损坏，进而影响桥梁的整体稳定性。廊屋的存在也会增加桥梁所受的水压力，特别是在洪水水位较高时，廊屋下部会承受较大的水压。因此，在设计和建造木拱廊桥时，需要综合考虑廊屋的结构强度、防风性能和防水性能，确保其在抗洪过程中能够发挥积极作用，同时避免因风荷载和水压力过大而导致的损坏。2.2.3桥台设计与抗洪关联桥台是木拱廊桥与河岸连接的重要结构，其设计与抗洪性能密切相关，对抵抗洪水冲刷和推力起着至关重要的作用。桥台的形式对其抗洪能力有着显著影响。常见的桥台形式有重力式桥台和轻型桥台等。重力式桥台通常采用较大体积的石材或混凝土砌筑而成，依靠自身的重量来抵抗洪水的推力和冲刷。这种桥台具有较高的稳定性和抗冲刷能力，能够承受较大的水平力和竖向力。在洪水来临时，重力式桥台能够凭借其巨大的自重，将洪水的冲击力分散到地基中，减少对桥台自身的破坏。而轻型桥台则相对轻巧，采用钢筋混凝土等材料制成，其结构形式较为灵活。轻型桥台虽然在自重方面不如重力式桥台，但通过合理的结构设计和配筋，也能够有效地抵抗洪水的作用。一些轻型桥台采用了悬臂式或框架式结构，能够在一定程度上减小洪水对桥台的冲击力，同时提高了桥台的承载能力。在实际应用中，需要根据木拱廊桥所在河道的水文条件、地质情况等因素，选择合适的桥台形式，以确保桥台在洪水中的稳定性。桥台的基础稳固性是其抗洪的关键因素之一。桥台的基础深入地下，与地基紧密结合，为桥台提供了稳定的支撑。在洪水作用下，桥台的基础不仅要承受桥台自身的重量和拱架传递过来的荷载，还要抵御洪水的冲刷和浸泡。如果基础不稳固，可能会导致桥台下沉、倾斜甚至倒塌，从而危及整个桥梁的安全。为了增强桥台基础的稳固性，通常会采取一系列措施。在基础施工时，会对地基进行加固处理，如采用夯实、换填等方法，提高地基的承载力。会增加基础的埋深，使基础深入到较稳定的土层中，减少洪水对基础的影响。还会在基础周围设置防护设施，如石笼、挡土墙等，防止洪水对基础的冲刷。在一些洪水灾害中，发现桥台基础稳固的木拱廊桥能够较好地抵御洪水的侵袭，而基础存在隐患的桥梁则容易出现桥台损坏的情况，这充分说明了桥台基础稳固性对桥梁抗洪性能的重要性。2.3木拱廊桥水毁案例分析2.3.1案例选取本研究选取了泰顺县的文兴桥和庆元县的兰溪桥作为典型案例进行深入分析。文兴桥位于浙江省泰顺县筱村镇东洋村，横跨玉溪之上，是一座具有重要历史文化价值的木拱廊桥。它始建于明朝永乐年间，清乾隆甲申年（1764年）重建，全长46.2米，净跨31.46米，矢高5.9米，桥屋16间。其独特的不对称结构在木拱廊桥中独树一帜，具有较高的研究价值。兰溪桥位于浙江省庆元县五大堡乡西洋村，始建于明万历二年（1574年），清乾隆五十九年（1794年）重建，桥长48.12米，净跨36.8米，矢高9.25米，桥屋21间。它是庆元县现存单孔跨度最大的木拱廊桥，其结构设计和建造工艺代表了当时的较高水平。选择这两座桥作为案例，主要是因为它们在结构类型、建造年代、地理位置等方面具有一定的代表性，能够更全面地反映木拱廊桥在洪水作用下的不同破坏情况和水毁原因。2.3.2水毁原因分析在2016年9月15日的超强台风“莫兰蒂”侵袭中，文兴桥遭受了严重的破坏。台风带来的暴雨引发了洪水，洪水的冲击力对文兴桥造成了多方面的损害。从洪水冲击力角度来看，强大的洪水携带大量的漂浮物，如树木、杂物等，以极高的速度冲击着桥体。这些漂浮物与桥体的碰撞产生了巨大的冲击力，导致桥体的部分构件松动、脱落。桥身的廊屋木板被冲毁，部分立柱被撞断，使得廊屋的结构完整性受到严重破坏。洪水的强大推力还作用于拱架，使拱架承受了远超设计荷载的水平力。在这种巨大的冲击力作用下，拱架的节点连接部位出现了松动，“三节苗”和“五节苗”之间的榫卯连接部分失效，导致拱架的稳定性下降。从结构破坏方面分析，文兴桥的拱架结构虽然具有一定的抗洪能力，但在此次超强洪水的冲击下，其薄弱环节逐渐显现。由于洪水的冲击力不均匀，导致拱架各部分受力不均，部分构件承受了过大的压力和弯矩。在拱架的顶部和底部，出现了木材的压溃和断裂现象，这使得拱架的承载能力大幅降低。廊屋的破坏也对桥体结构产生了连锁反应。廊屋作为桥体的一部分，不仅起到遮风挡雨的作用，还对拱架起到一定的支撑和稳定作用。廊屋的受损使得拱架失去了部分侧向支撑，在洪水的持续冲击下，更容易发生变形和倒塌。基础冲刷也是文兴桥水毁的重要原因之一。洪水的强大水流对桥的基础产生了强烈的冲刷作用。文兴桥的桥台基础采用的是石材砌筑，在长期的洪水冲刷下，基础周围的土体被逐渐掏空，导致桥台的稳定性下降。桥台的倾斜和下沉使得拱架的支撑条件发生改变，拱架所承受的荷载分布不均，进一步加剧了桥体的破坏。在洪水过后的现场勘查中，发现桥台周围的土体出现了明显的冲刷痕迹，部分石材已经外露，这充分说明了基础冲刷对文兴桥水毁的影响。兰溪桥在历史上也多次遭受洪水侵袭，其水毁原因同样与洪水冲击力、结构破坏和基础冲刷密切相关。在某次洪水灾害中，兰溪桥的桥身受到了洪水的猛烈冲击。洪水的高流速和大流量使得桥体受到了巨大的水动力荷载。桥身的侧面受到洪水的直接冲击，产生了较大的压力，导致桥身的部分木板和构件被冲毁。洪水携带的大量漂浮物对桥体的撞击也造成了严重的破坏，一些关键的连接部位被撞坏，使得桥体的整体性受到影响。兰溪桥的结构在洪水作用下也出现了破坏。拱架作为桥体的主要承重结构，在洪水的冲击下，承受了巨大的压力。由于拱架的木材长期受到自然环境的侵蚀，其强度有所降低，在洪水的冲击下更容易发生损坏。在拱架的节点处，由于榫卯连接的松动，导致节点的传力性能下降，使得拱架的稳定性受到威胁。廊屋的破坏也削弱了桥体的整体结构性能。廊屋的部分屋顶被洪水掀翻，立柱被冲倒，这不仅影响了廊屋的正常使用功能，还使得桥体的抗风、抗震能力下降，在后续的洪水冲击中更容易遭受破坏。基础冲刷同样是兰溪桥水毁的重要因素。兰溪桥位于河流弯道处，水流速度较快，对基础的冲刷作用更为强烈。在洪水期间，河水对桥台基础的冲刷使得基础周围的防护层被破坏，基础逐渐暴露在洪水之中。随着冲刷的加剧，基础的埋深减小，稳定性降低，无法为桥体提供足够的支撑。这导致桥体在洪水的冲击下发生倾斜和位移，最终导致桥体的破坏。在对兰溪桥的修复过程中，发现桥台基础的石材出现了严重的磨损和松动现象，这进一步证实了基础冲刷对桥体水毁的影响。2.3.3案例启示通过对文兴桥和兰溪桥水毁案例的分析，我们可以得到以下对提高木拱廊桥抗洪性能研究和保护的重要启示。在结构设计优化方面，应进一步加强拱架结构的稳定性。通过改进榫卯连接方式，增加连接的强度和可靠性，使拱架各构件之间能够更好地协同工作，共同抵抗洪水的冲击力。可以采用新型的连接材料或加固措施，提高节点的承载能力和抗震性能。合理设计廊屋的结构，增强其与拱架的连接，使廊屋在洪水冲击下能够更好地发挥对拱架的支撑和稳定作用。可以增加廊屋的立柱数量，优化立柱的布置方式，提高廊屋的整体刚度。加强桥台基础的设计和加固，采用抗冲刷能力强的基础形式和材料，如扩大基础、桩基础等，并在基础周围设置有效的防护措施，如石笼、挡土墙等，减少洪水对基础的冲刷。材料选择与维护对于提高木拱廊桥的抗洪性能也至关重要。应选用质量优良、耐久性好的木材作为桥体的主要材料。在木材的选择上，优先考虑抗腐蚀、抗虫蛀能力强的树种，并对木材进行必要的防腐、防虫处理。定期对木拱廊桥进行维护和检查，及时发现和修复木材的损坏、腐朽等问题。建立完善的维护制度，明确维护的周期和内容，确保桥体的结构安全。在洪水来临前，对桥体进行全面的检查和加固，如对松动的构件进行紧固，对受损的部位进行修复，增强桥体的抗洪能力。预警与应急措施的建立也是保护木拱廊桥的重要环节。建立洪水预警系统，及时获取洪水的信息，提前做好木拱廊桥的防护准备。与当地的气象、水利等部门建立密切的合作关系，共享洪水监测数据，提高预警的准确性和及时性。制定应急预案，明确在洪水发生时的应对措施和责任分工。在洪水发生时，能够迅速采取有效的防护措施，如设置防洪挡板、疏散桥上人员等，减少洪水对桥体的破坏。加强对周边环境的管理，减少漂浮物对桥体的冲击。定期清理河道，防止河道堵塞，降低洪水对桥体的冲击力。三、数值模拟技术在木拱廊桥抗洪研究中的应用3.1数值模拟方法与工具3.1.1计算流体动力学（CFD）原理计算流体动力学（CFD）是一种通过数值模拟来研究流体流动和传热现象的方法，其核心在于利用计算机求解流体力学和传热学的基本方程，如纳维-斯托克斯方程（N-S方程）、能量守恒方程等，从而得到流体流动和传热的数值解。在木拱廊桥抗洪研究中，CFD主要用于模拟洪水流动以及洪水与桥梁结构的相互作用。纳维-斯托克斯方程是描述粘性不可压缩流体动量守恒的运动方程，在直角坐标系下，其矢量形式为：\rho(\frac{\partial\vec{v}}{\partialt}+(\vec{v}\cdot\nabla)\vec{v})=-\nablap+\mu\nabla^2\vec{v}+\vec{F}其中，\rho为流体密度，\vec{v}为速度矢量，t为时间，p为压力，\mu为动力粘度，\vec{F}为质量力。该方程反映了流体的惯性力、压力梯度力、粘性力和质量力之间的平衡关系。在实际应用中，由于N-S方程的非线性和复杂性，通常难以获得解析解，因此需要采用数值方法进行求解。CFD常用的数值解法包括有限差分法、有限元法和有限体积法等。有限差分法是将求解域划分为差分网格，用有限个网格节点代替连续的求解域，然后将偏微分方程的导数用差商代替，推导出含有离散点上有限个未知数的差分方程组，通过求解差分方程组得到数值解。有限元法是将求解域离散为有限个单元，通过对每个单元进行插值和加权余量法，建立单元的有限元方程，然后将所有单元的方程组装成总体有限元方程，求解该方程得到数值解。有限体积法是将计算区域划分为一系列不重复的控制体积，并使每个网格点周围有一个控制体积，将待解的微分方程对每一个控制体积积分，从而得到一组离散方程。在木拱廊桥抗洪模拟中，有限体积法因其守恒性好、物理意义明确等优点而被广泛应用。以模拟洪水对木拱廊桥的冲击为例，首先需要根据木拱廊桥的实际尺寸和周围河道的地形，建立包含桥梁和水流区域的计算模型。将该模型离散为大量的网格单元，在每个网格单元上应用N-S方程和连续性方程，并结合合适的边界条件，如入口边界的流速、压力条件，出口边界的自由出流条件，以及桥梁和河道边界的无滑移条件等。通过迭代求解这些方程，逐步计算出每个网格单元内的流体速度、压力等物理量，从而得到整个流场的分布情况。通过CFD模拟，可以清晰地观察到洪水在木拱廊桥周围的流动形态，如水流的绕流、分离、漩涡的产生等，以及木拱廊桥所受到的水动力荷载，如压力分布、冲击力大小和方向等，为后续分析木拱廊桥的结构响应提供重要依据。3.1.2有限元分析（FEA）原理有限元分析（FEA）是一种用于求解工程力学问题的数值方法，其基本思想是将连续的求解域离散为有限个单元，通过对每个单元进行力学分析，建立单元的力学方程，然后将所有单元的方程组装成总体力学方程，求解该方程得到整个求解域的近似解。在木拱廊桥抗洪性能研究中，FEA主要用于分析木拱廊桥在洪水作用下的结构力学响应，包括应力、应变分布以及变形情况等。有限元分析的基本步骤包括：首先是结构离散化，将木拱廊桥的结构按照一定的规则划分为有限个单元，如梁单元、杆单元、壳单元等。每个单元通过节点与其他单元相连，节点是单元之间传递力和位移的关键部位。对于木拱廊桥的拱架结构，可以采用梁单元来模拟其受力行为，因为梁单元能够较好地考虑弯曲、拉伸和剪切等多种受力情况。而对于廊屋的屋面和墙体等结构，可以采用壳单元进行模拟，壳单元能够有效地模拟薄板结构在平面内和平面外的受力特性。材料属性定义是有限元分析的重要环节，需要根据木拱廊桥实际使用的木材特性，确定其弹性模量、泊松比、密度等材料参数。不同种类的木材具有不同的力学性能，在定义材料属性时，应尽可能准确地反映木材的实际情况。对于常见的杉木，其弹性模量一般在10-12GPa之间，泊松比约为0.3。考虑木材的各向异性特性，即在不同方向上力学性能存在差异，在有限元模型中合理设置材料的各向异性参数，以提高分析结果的准确性。边界条件的设定直接影响有限元分析的结果，在木拱廊桥的有限元模型中，需要根据实际情况设置边界条件。桥台与地基的连接部位通常设置为固定约束，限制其在三个方向的位移和转动，以模拟桥台在地基上的固定状态。对于桥梁与桥墩的连接处，可以根据实际的连接方式，设置为铰接或刚接约束，铰接约束只限制垂直于连接面的位移，允许绕节点的转动；刚接约束则同时限制位移和转动，确保节点处的力和力矩能够有效传递。荷载施加是有限元分析中模拟实际工况的关键步骤，将CFD模拟得到的水动力荷载，包括压力、冲击力等，准确地施加到木拱廊桥的有限元模型上。根据水动力荷载在桥面上的分布情况，将其转化为节点荷载或面荷载施加到相应的单元上。考虑洪水作用的动态特性，在有限元分析中采用动态分析方法，如瞬态动力学分析，以更真实地模拟木拱廊桥在洪水冲击下的结构响应过程。通过求解总体力学方程，可以得到木拱廊桥在洪水作用下各单元的应力、应变和位移等结果。分析这些结果，能够确定木拱廊桥结构中的薄弱部位，如拱架的节点处、廊屋与拱架的连接部位等，这些部位在洪水作用下可能承受较大的应力和变形，容易发生损坏。通过有限元分析，还可以评估木拱廊桥在不同洪水工况下的结构安全性，为制定保护措施和修复方案提供科学依据。3.1.3常用模拟软件介绍在木拱廊桥抗洪性能研究中，Fluent和ANSYS是两款常用的模拟软件，它们各自具有独特的特点和强大的功能，为研究工作提供了有力的支持。Fluent是一款广泛应用于计算流体动力学领域的专业软件，由ANSYS公司开发。其具有丰富的物理模型库，涵盖了各种流体流动和传热现象，如湍流、多相流、燃烧等，能够满足不同复杂程度的工程问题模拟需求。在木拱廊桥抗洪研究中，Fluent的主要优势在于其强大的流体模拟能力。它采用有限体积法对控制方程进行离散求解，能够准确地模拟洪水在木拱廊桥周围的流动特性和水动力荷载分布。Fluent提供了多种湍流模型，如标准k-ε模型、RNGk-ε模型、Realizablek-ε模型等，用户可以根据实际情况选择合适的模型来模拟洪水的湍流流动。对于木拱廊桥周围复杂的水流情况，如水流的分离、漩涡的产生等，Fluent能够通过合理选择湍流模型和设置相关参数，准确地捕捉这些流动现象，为分析洪水对木拱廊桥的作用提供详细的流场信息。Fluent的网格划分功能也十分强大，支持多种网格类型，包括结构化网格、非结构化网格以及混合网格。对于木拱廊桥这种复杂的几何结构，非结构化网格能够更好地适应其形状，提高网格质量和计算精度。Fluent还提供了网格自适应功能，能够根据计算结果自动调整网格密度，在流场变化剧烈的区域，如木拱廊桥的桥墩周围、拱架附近等，加密网格，以提高计算的准确性。在模拟洪水与木拱廊桥的相互作用时，Fluent能够通过设置合理的边界条件，如入口边界的流速、压力条件，出口边界的自由出流条件，以及桥梁和河道边界的无滑移条件等，准确地模拟洪水对木拱廊桥的冲击过程，得到木拱廊桥所受到的水动力荷载，如压力分布、冲击力大小和方向等。Fluent还具有强大的后处理功能，能够以直观的方式展示模拟结果，如绘制速度场、压力场的云图和矢量图，生成力随时间变化的曲线等，方便用户对模拟结果进行分析和理解。ANSYS是一款功能全面的大型通用有限元分析软件，除了具备强大的结构分析功能外，还涵盖了流体、热、电磁等多个物理场的分析能力。在木拱廊桥抗洪性能研究中，ANSYS主要用于木拱廊桥的结构力学分析。它拥有丰富的单元库，包括梁单元、杆单元、壳单元、实体单元等，能够满足木拱廊桥复杂结构的建模需求。对于木拱廊桥的拱架结构，可以使用梁单元来模拟其受力行为，梁单元能够准确地考虑弯曲、拉伸和剪切等多种受力情况，为分析拱架在洪水作用下的力学响应提供了有效的手段。对于廊屋的屋面和墙体等结构，壳单元能够很好地模拟其在平面内和平面外的受力特性，确保对廊屋结构的分析准确性。ANSYS的材料模型库也非常丰富，能够定义各种材料的力学性能，包括木材的弹性模量、泊松比、密度等参数。考虑木材的各向异性特性，ANSYS可以通过设置相应的材料参数，准确地模拟木材在不同方向上的力学行为。在边界条件设置方面，ANSYS提供了多种约束类型，如固定约束、铰接约束、弹性约束等，用户可以根据木拱廊桥的实际支撑情况和连接方式，合理地设置边界条件，确保模型能够真实地反映木拱廊桥的受力状态。在荷载施加方面，ANSYS能够方便地将CFD模拟得到的水动力荷载施加到木拱廊桥的有限元模型上，通过瞬态动力学分析等方法，模拟木拱廊桥在洪水冲击下的动态响应过程。ANSYS的后处理功能同样强大，能够对计算结果进行深入分析，如查看结构的应力、应变分布云图，提取关键部位的应力、应变和位移数据，生成结构变形动画等，帮助用户全面了解木拱廊桥在洪水作用下的结构力学响应，为评估木拱廊桥的抗洪性能和制定保护措施提供科学依据。3.2建立木拱廊桥数值模型3.2.1模型简化与假设在建立木拱廊桥数值模型时，为了在保证计算精度的前提下提高计算效率，需对复杂的木拱廊桥结构进行合理的简化与假设。考虑到木拱廊桥的结构特点和实际受力情况，对一些次要结构和细节进行简化处理。例如，对于桥面上的一些附属设施，如栏杆、神龛等，由于它们对木拱廊桥在洪水作用下的整体力学性能影响较小，在模型中可予以忽略。这些附属设施的重量相对较轻，且其与桥体的连接方式较为简单，在洪水冲击时，它们所承受的荷载主要是自身的重力和较小的风力，对桥体的主要受力构件影响不大。而对于木拱廊桥的主体结构，如拱架、廊屋和桥台等，则进行详细建模，以准确模拟其力学行为。在模型假设方面，基于木材的力学性能特点和木拱廊桥的实际工作状态，做出以下假设。假设木材为各向同性材料，尽管实际上木材具有一定的各向异性特性，但其在顺纹方向和横纹方向的力学性能差异在某些情况下对整体分析结果影响较小。在木拱廊桥的主要受力方向，如拱架承受竖向荷载和水平力时，顺纹方向的力学性能起主导作用，因此在一定程度上忽略木材的各向异性，将其视为各向同性材料，可简化计算过程，同时也能满足工程分析的精度要求。假设木拱廊桥各构件之间的连接为刚性连接。虽然在实际结构中，木拱廊桥的构件连接采用榫卯等方式，具有一定的柔性和变形能力，但在数值模拟中，为了便于分析和计算，将连接视为刚性连接。这种假设能够简化模型的力学分析，突出主要结构构件的受力特性。在后续的研究中，可以通过进一步的试验和分析，考虑连接的柔性对木拱廊桥结构性能的影响，对模型进行修正和完善。通过以上模型简化与假设，既能减少计算量，提高计算效率，又能抓住木拱廊桥结构在洪水作用下的主要力学特征，为后续的数值模拟分析奠定基础。在简化和假设过程中，充分考虑了木拱廊桥的实际结构和受力情况，确保模型能够合理地反映木拱廊桥的抗洪性能。3.2.2材料参数设定准确设定材料参数是建立木拱廊桥数值模型的关键环节，它直接影响到模拟结果的准确性和可靠性。木拱廊桥主要由木材、石材等材料构成，不同材料具有不同的力学性能，因此需要根据实际情况合理确定其参数。对于木材，其力学性能参数主要包括弹性模量、泊松比、密度等。在木拱廊桥中，常用的木材种类有杉木、松木等，不同种类木材的力学性能存在一定差异。以杉木为例，其弹性模量一般在10-12GPa之间，泊松比约为0.3。在数值模拟中，可根据实际使用的木材种类和相关标准，选取合适的弹性模量和泊松比。考虑木材的含水率对其力学性能的影响，含水率的变化会导致木材的弹性模量、强度等参数发生改变。一般来说，木材的含水率越高，其弹性模量和强度越低。在设定木材参数时，需参考实际木材的含水率情况，对参数进行适当调整。根据相关研究和试验数据，当木材含水率从12%增加到20%时，其弹性模量可能会降低10%-20%。因此，在实际工程中，可通过测量木材的含水率，结合经验公式或试验数据，对木材的力学性能参数进行修正，以提高模型的准确性。石材作为桥台等结构的主要材料，其参数设定同样重要。石材的主要力学参数包括抗压强度、弹性模量、泊松比等。常见的用于桥台的石材有花岗岩、石灰岩等，花岗岩的抗压强度较高，一般在100-200MPa之间，弹性模量约为50-70GPa，泊松比在0.2-0.3之间。在数值模型中，根据石材的具体种类和质量，合理确定其抗压强度、弹性模量等参数。考虑石材的耐久性和抗风化性能，在长期的使用过程中，石材可能会受到风化、侵蚀等作用，导致其力学性能下降。在设定石材参数时，需考虑这些因素，适当降低石材的力学性能指标，以更真实地反映石材在实际工程中的性能。通过对石材进行长期的风化试验，获取其力学性能随时间变化的数据，根据试验结果对石材参数进行修正，确保模型能够准确模拟桥台在洪水作用下的力学行为。3.2.3边界条件设定边界条件的合理设定是数值模拟中至关重要的环节，它直接影响到模拟结果的准确性和可靠性，对于木拱廊桥抗洪性能的研究具有关键作用。在模拟木拱廊桥与洪水相互作用时，需要准确设置洪水流量、流速、水位等边界条件，以真实反映木拱廊桥在洪水中的受力环境。洪水流量是描述洪水规模的重要参数，其大小直接决定了洪水对木拱廊桥的冲击力和破坏力。根据木拱廊桥所在流域的水文资料，如历史洪水记录、水文监测数据等，确定不同重现期下的洪水流量。通过对历史洪水数据的统计分析，结合水文计算方法，如推理公式法、经验公式法等，推算出不同重现期（如50年一遇、100年一遇等）的洪水流量。对于某一特定的木拱廊桥，经过计算确定其50年一遇的洪水流量为500m³/s，100年一遇的洪水流量为800m³/s。在数值模拟中，将这些计算得到的洪水流量作为入口边界条件，输入到CFD模型中，以模拟不同规模洪水对木拱廊桥的作用。流速是影响洪水冲击力的关键因素之一，不同位置和不同时刻的流速变化会导致木拱廊桥所受的水动力荷载发生改变。在设定流速边界条件时，考虑洪水在河道中的流动特性，如流速分布的不均匀性、流速随时间的变化等。通过实地测量或参考相关研究成果，获取木拱廊桥周围河道在不同洪水工况下的流速分布数据。在河流中心和靠近岸边的流速存在差异，在洪水过程中，流速也会随着水位的上升和下降而发生变化。在数值模拟中，根据这些实际情况，采用合适的流速分布模型，如对数流速分布模型、幂律流速分布模型等，来设定流速边界条件。对于某一河道，根据测量数据和相关研究，确定在洪水期采用对数流速分布模型，在入口边界处，根据不同的洪水工况，设定相应的流速值，以准确模拟洪水在木拱廊桥周围的流动情况和对桥体的冲击作用。水位是洪水的重要特征之一，它决定了木拱廊桥被淹没的程度和所受水压力的大小。通过对木拱廊桥所在流域的水位监测数据进行分析，结合历史洪水记录和地形地貌信息，确定不同洪水工况下的水位变化过程。在模拟中，将水位作为边界条件，设置为随时间变化的函数。在洪水初期，水位逐渐上升，达到峰值后又逐渐下降，在数值模型中，根据实际的水位变化曲线，设置相应的边界条件，以模拟木拱廊桥在不同水位下的受力情况。考虑水位变化对木拱廊桥基础的影响，如基础的冲刷、浸泡等，在设定边界条件时，综合考虑这些因素，确保模型能够全面反映洪水对木拱廊桥的作用。3.3数值模拟结果验证与分析3.3.1与实际监测数据对比为了验证数值模拟结果的准确性和可靠性，将模拟结果与实际洪水监测数据进行对比分析。以泰顺县某木拱廊桥为例，在2019年的一次洪水中，对该桥进行了实时监测，获取了洪水过程中桥体的位移、应力等数据。同时，运用建立的数值模型，模拟相同洪水工况下该木拱廊桥的受力和变形情况。在位移对比方面，实际监测数据显示，在洪水峰值时刻，桥体跨中部位的竖向位移达到了35mm。而数值模拟结果表明，在相同的洪水条件下，桥体跨中部位的竖向位移为33mm。两者之间的误差在合理范围内，相对误差约为5.7%。这说明数值模拟能够较为准确地预测桥体在洪水作用下的竖向位移情况。在水平位移方面，实际监测到桥体桥台部位的水平位移为12mm，数值模拟结果为11mm，相对误差约为8.3%。通过对比可以看出，数值模拟结果与实际监测数据在位移方面具有较好的一致性，验证了数值模型在模拟桥体位移响应上的准确性。在应力对比方面，实际监测数据显示，桥体拱架的某些关键部位在洪水作用下的最大应力达到了8MPa。数值模拟结果表明，这些关键部位的最大应力为7.8MPa，相对误差约为2.5%。这表明数值模拟能够准确地模拟桥体在洪水作用下的应力分布情况，为评估桥体的结构安全性提供了可靠的依据。通过将模拟结果与实际监测数据进行全面对比，验证了数值模型的准确性和可靠性。这不仅为进一步分析木拱廊桥在洪水作用下的力学响应提供了有力支持，也为基于数值模拟的木拱廊桥抗洪性能研究奠定了坚实的基础。3.3.2模拟结果分析通过对数值模拟结果的深入分析，能够全面了解洪水作用下木拱廊桥的应力、应变和位移响应，为评估其抗洪性能提供关键依据。在应力响应方面，当洪水冲击木拱廊桥时，拱架作为主要受力结构，承受着较大的应力。从模拟结果的应力云图可以清晰地看到，拱架的节点处和拱顶部位应力集中较为明显。在节点处，由于“三节苗”和“五节苗”通过“牛头”连接，力的传递较为复杂，导致节点处承受着较大的弯矩和剪力，从而产生较高的应力。在某次模拟中，节点处的最大应力达到了10MPa，接近木材的许用应力。拱顶部位由于直接承受桥面传来的荷载和洪水的冲击力，也承受着较大的压力，应力值相对较高。廊屋的立柱和横梁在洪水作用下也会产生一定的应力，尤其是靠近拱架的立柱，由于受到拱架传来的水平力和自身重力的作用，应力分布较为复杂。应变响应与应力密切相关，在应力较大的部位，应变也相应较大。拱架节点处和拱顶部位的应变明显高于其他部位。在节点处，由于木材的受力复杂，不仅存在轴向应变，还存在剪切应变。通过模拟结果可以计算出节点处的最大剪切应变达到了0.005，这表明节点处的木材在洪水作用下发生了一定程度的剪切变形。拱顶部位的轴向应变也较为显著，最大轴向应变达到了0.003，这说明拱顶部位在压力作用下产生了一定的压缩变形。廊屋的立柱和横梁在应力作用下也会产生相应的应变，立柱的轴向应变和横梁的弯曲应变都在一定程度上影响着廊屋的结构稳定性。位移响应是衡量木拱廊桥在洪水作用下结构变形的重要指标。模拟结果显示，桥体在洪水作用下主要发生竖向位移和水平位移。竖向位移在桥体跨中部位最为明显，这是由于跨中部位承受的荷载最大，在洪水冲击力和自身重力的共同作用下，跨中部位产生了较大的挠度。在不同洪水工况下，跨中部位的竖向位移随着洪水流量和流速的增加而增大。在一次模拟中，当洪水流量为500m³/s时，跨中部位的竖向位移为25mm；当洪水流量增加到800m³/s时，跨中部位的竖向位移增大到40mm。水平位移主要发生在桥台和桥墩部位，这是由于洪水的冲击力对桥台和桥墩产生了水平推力，导致其发生一定程度的水平移动。水平位移的大小与洪水的流速和流向密切相关，在流速较大且流向与桥体垂直的情况下，桥台和桥墩的水平位移会显著增大。通过对模拟结果中应力、应变和位移响应的分析，能够明确木拱廊桥在洪水作用下的薄弱部位和潜在风险点，为制定针对性的抗洪保护措施提供科学依据。四、木拱廊桥抗洪性能影响因素的数值模拟分析4.1洪水参数对木拱廊桥抗洪性能的影响4.1.1流速变化影响为深入探究流速变化对木拱廊桥抗洪性能的影响，运用数值模拟技术，构建了木拱廊桥与洪水相互作用的数值模型，设定不同的流速工况，模拟洪水对木拱廊桥的冲击力和结构响应。在模拟过程中，保持其他洪水参数（如水位、流量等）不变，将流速分别设定为1m/s、2m/s、3m/s、4m/s和5m/s。随着流速的逐渐增大，洪水对木拱廊桥的冲击力呈现出显著的变化。当流速为1m/s时，木拱廊桥所受的冲击力相对较小，桥体各构件的应力和应变也处于较低水平。随着流速增大到2m/s，冲击力明显增大，桥体的拱架和桥墩部位所受应力有所增加，部分关键节点的应变也开始增大。当流速达到3m/s时，冲击力进一步增强，拱架的节点处应力集中现象更加明显，部分节点的应力接近木材的许用应力，存在一定的安全隐患。流速增大到4m/s时，桥体所受冲击力大幅增加，拱架出现明显变形，部分构件的应变超过了允许范围，结构的稳定性受到严重威胁。当流速达到5m/s时，木拱廊桥所受冲击力达到极大值，桥体结构出现严重破坏，部分构件断裂，拱架坍塌，桥梁失去承载能力。通过对模拟结果的详细分析，绘制了冲击力随流速变化的曲线，以及桥体关键部位（如拱架节点、桥墩底部等）的应力和应变随流速变化的曲线。从冲击力曲线可以看出，冲击力与流速的平方成正比，即流速的微小变化会导致冲击力的大幅增加。在桥体关键部位的应力应变曲线中，随着流速的增大，应力和应变均呈现出非线性增长的趋势，当流速超过一定阈值时，应力和应变急剧增大，表明桥体结构已接近破坏极限。这些曲线直观地展示了流速变化对木拱廊桥抗洪性能的影响规律，为评估木拱廊桥在不同流速洪水下的安全性提供了重要依据。4.1.2水位变化影响为了深入了解不同水位高度对木拱廊桥受力和稳定性的影响，通过数值模拟设置了多种水位工况，包括低水位、中水位和高水位，分别模拟不同水位条件下木拱廊桥的受力和变形情况。在低水位工况下，水位高度仅达到桥墩高度的1/3。此时，木拱廊桥所受的水压力相对较小，桥体主要承受自身重力和少量的水流冲击力。通过模拟分析发现，桥体的应力分布较为均匀，拱架和桥墩的应力水平较低，结构变形也较小。在拱架的关键节点处，应力值仅为木材许用应力的20%左右，桥墩底部的应力也处于安全范围内。桥体的竖向位移和水平位移都非常小，对桥梁的正常使用和结构安全影响较小。当中水位工况时，水位高度达到桥墩高度的2/3。随着水位的升高，木拱廊桥所受的水压力显著增加，水流冲击力也相应增大。模拟结果显示，桥体的应力分布发生了明显变化，拱架的节点处和桥墩与桥面连接部位的应力集中现象较为明显。在拱架节点处，应力值达到木材许用应力的40%左右，桥墩与桥面连接部位的应力也有所增加。桥体的竖向位移和水平位移也有所增大，其中竖向位移较之前增加了约50%，水平位移增加了约30%。虽然结构仍处于稳定状态，但已接近安全阈值，需要密切关注。在高水位工况下，水位高度超过桥墩高度，桥体大部分被淹没。此时，木拱廊桥所受的水压力和水流冲击力达到最大值，对桥体结构产生了巨大的威胁。模拟结果表明，桥体的应力分布极不均匀，拱架的多个部位出现应力集中，部分节点的应力超过了木材的许用应力，导致木材出现局部破坏。桥墩底部的应力也急剧增大，接近其极限承载能力。桥体的竖向位移和水平位移大幅增加，竖向位移较中水位时增加了约100%，水平位移增加了约80%。桥体结构出现明显的倾斜和变形，稳定性受到严重破坏，随时可能发生倒塌。通过对不同水位工况下木拱廊桥受力和稳定性的模拟分析，可以得出结论：水位高度的增加会显著增大木拱廊桥所受的水压力和水流冲击力，导致桥体应力分布不均，关键部位应力集中，结构变形增大，稳定性降低。在高水位情况下，木拱廊桥的结构安全面临极大挑战，需要采取有效的防护措施来确保其在洪水中的安全。4.1.3流量变化影响洪水流量的变化对木拱廊桥抗洪性能有着显著影响，通过数值模拟不同流量工况，深入研究其对桥梁抗洪性能的作用机制。在模拟过程中，设定了一系列不同的洪水流量值，分别为100m³/s、200m³/s、300m³/s、400m³/s和500m³/s。随着流量的逐渐增大，木拱廊桥所面临的抗洪挑战也愈发严峻。当流量为100m³/s时，洪水对木拱廊桥的冲击力相对较小，桥体结构基本能够保持稳定。通过模拟分析发现，桥体各构件的应力和应变处于较低水平，拱架的关键节点应力仅为木材许用应力的15%左右，桥墩底部的应力也在安全范围内。桥体的位移变形微小，对桥梁的正常使用和结构安全影响不大。随着流量增大到200m³/s，洪水的冲击力明显增强，桥体所受荷载增大。模拟结果显示，拱架节点处的应力有所增加，达到木材许用应力的25%左右，桥墩与桥面连接部位的应力也相应增大。桥体的位移变形开始显现，竖向位移和水平位移分别增加了约3mm和2mm。虽然结构仍处于稳定状态，但已需要关注应力和位移的变化趋势。当流量达到300m³/s时，洪水对木拱廊桥的作用更加显著，桥体所受冲击力大幅增大。此时，拱架的多个部位出现应力集中现象，部分节点的应力达到木材许用应力的40%左右，接近安全阈值。桥墩底部的应力也显著增加，对桥墩的稳定性构成一定威胁。桥体的竖向位移和水平位移进一步增大，分别达到约6mm和4mm，结构的变形开始影响其正常使用功能。流量增大到400m³/s时，洪水的冲击力使木拱廊桥的结构安全面临较大挑战。模拟结果表明，拱架的部分节点应力超过木材许用应力，木材出现局部破坏的迹象。桥墩底部的应力接近其极限承载能力，桥墩的稳定性受到严重影响。桥体的竖向位移和水平位移急剧增大，分别达到约10mm和7mm，桥体结构出现明显的倾斜和变形，稳定性降低。当流量达到500m³/s时，洪水对木拱廊桥的冲击力达到极大值，桥体结构遭受严重破坏。模拟结果显示，拱架多处构件断裂，节点连接失效，拱架整体坍塌。桥墩也发生严重位移和倾斜，无法支撑桥体重量。桥体的竖向位移和水平位移超出了允许范围，桥梁完全失去承载能力和使用功能。通过对不同流量工况下木拱廊桥抗洪性能的模拟研究，可以清晰地看到，随着洪水流量的增大，木拱廊桥所受的冲击力、应力和变形均显著增加，结构的稳定性逐渐降低，当流量达到一定程度时，桥体结构将遭受严重破坏，失去抗洪能力。这为评估木拱廊桥在不同洪水流量下的安全性提供了重要依据，也为制定有效的抗洪保护措施提供了参考。4.2桥梁结构参数对抗洪性能的影响4.2.1拱架形式变化影响为深入探究拱架形式对木拱廊桥抗洪性能的影响，构建了三种不同拱架形式的木拱廊桥数值模型，分别为传统贯木拱、改进型贯木拱和新型组合拱。传统贯木拱采用常见的“三节苗”和“五节苗”相互穿插的结构形式；改进型贯木拱在传统结构基础上，对“三节苗”和“五节苗”的连接方式进行了优化，增强了节点的连接强度；新型组合拱则引入了斜撑结构，与传统贯木拱相结合，形成了一种新的结构形式。在模拟过程中，设置相同的洪水工况，包括流速为3m/s、水位达到桥墩高度的2/3、流量为300m³/s。模拟结果显示，在相同洪水条件下，三种拱架形式的木拱廊桥抗洪性能存在显著差异。传统贯木拱的木拱廊桥，在洪水冲击下，拱架节点处应力集中明显，部分节点的应力接近木材的许用应力。由于节点连接方式相对传统，在较大的洪水冲击力作用下，节点容易出现松动，导致拱架的整体性受到影响，结构变形较大。在此次模拟中，拱架节点处的最大应力达到了9MPa，桥体跨中的竖向位移达到了8mm。改进型贯木拱的木拱廊桥，由于优化了节点连接方式，节点处的应力分布得到改善，应力集中现象有所缓解。在相同洪水工况下，节点处的最大应力降低到了7MPa，桥体跨中的竖向位移减小到了6mm。这表明改进后的节点连接方式能够有效提高拱架的整体性和稳定性，增强木拱廊桥的抗洪能力。新型组合拱的木拱廊桥，由于斜撑结构的引入，拱架的受力更加均匀，应力集中现象得到了进一步改善。斜撑结构能够分担部分洪水冲击力，减轻拱架的负担，使得桥体在洪水作用下的变形明显减小。在模拟中，节点处的最大应力仅为5MPa，桥体跨中的竖向位移为4mm。新型组合拱在抗洪性能方面表现出明显的优势，能够更好地适应洪水的冲击。通过对不同拱架形式木拱廊桥抗洪性能的对比分析，可以得出结论：拱架形式的变化对木拱廊桥的抗洪性能有着重要影响。改进节点连接方式和引入斜撑结构等优化措施，能够有效提高拱架的整体性和稳定性，降低桥体在洪水中的应力和变形，增强木拱廊桥的抗洪能力。4.2.2廊屋结构变化影响廊屋结构的变化对木拱廊桥抗洪性能有着显著影响，通过构建不同廊屋结构的数值模型，包括有无廊屋、廊屋结构形式变化等，深入研究其在洪水中的作用。首先对比有廊屋和无廊屋的木拱廊桥模型在相同洪水工况下的抗洪性能。设置洪水流速为3m/s、水位达到桥墩高度的2/3、流量为300m³/s。模拟结果显示，有廊屋的木拱廊桥在洪水冲击下，整体稳定性明显优于无廊屋的木拱廊桥。有廊屋时，廊屋的重量增加了桥体的自重，使得桥体与基础之间的摩擦力增大，抵抗洪水推力的能力增强。在此次模拟中，有廊屋的木拱廊桥桥台水平位移为5mm，而无廊屋的木拱廊桥桥台水平位移达到了8mm。廊屋还能在一定程度上阻挡洪水对拱架和桥面的直接冲击，起到缓冲作用，减少了拱架和桥面的受损程度。接着研究廊屋结构形式变化对木拱廊桥抗洪性能的影响。构建了三种不同廊屋结构形式的模型，分别为传统穿斗式廊屋、加强型穿斗式廊屋和框架式廊屋。传统穿斗式廊屋采用常见的穿斗式木构架；加强型穿斗式廊屋在传统结构基础上，增加了立柱和横梁的数量，增强了廊屋的结构强度；框架式廊屋则采用框架结构，具有更高的整体性和稳定性。在相同洪水工况下，三种廊屋结构形式的木拱廊桥抗洪性能存在差异。传统穿斗式廊屋的木拱廊桥，在洪水冲击下，廊屋部分构件出现较大变形，部分立柱和横梁的应力接近木材的许用应力。由于结构强度相对有限，在较大的洪水冲击力作用下，廊屋的稳定性受到一定影响。在此次模拟中，廊屋立柱的最大应力达到了7MPa，部分横梁出现明显弯曲变形。加强型穿斗式廊屋的木拱廊桥，由于增加了立柱和横梁的数量，廊屋的结构强度得到提高，在洪水冲击下的变形明显减小。在相同洪水工况下，廊屋立柱的最大应力降低到了5MPa，横梁的弯曲变形也得到有效控制。这表明加强型穿斗式廊屋能够更好地承受洪水的冲击，增强了木拱廊桥的抗洪性能。框架式廊屋的木拱廊桥，由于其结构整体性和稳定性较高，在洪水冲击下表现出较好的抗洪性能。廊屋各构件的应力分布较为均匀，变形较小。在模拟中，廊屋立柱的最大应力仅为3MPa，桥体的整体位移也相对较小。框架式廊屋在抗洪性能方面具有明显优势，能够为木拱廊桥提供更有效的保护。通过对不同廊屋结构木拱廊桥抗洪性能的研究，可以得出结论：廊屋结构的变化对木拱廊桥的抗洪性能有着重要影响。有廊屋的木拱廊桥抗洪性能优于无廊屋的木拱廊桥，加强廊屋结构强度和采用框架式廊屋等优化措施，能够有效提高廊屋的稳定性和承载能力，增强木拱廊桥的抗洪能力。4.2.3桥台尺寸变化影响桥台尺寸的改变对木拱廊桥抵抗洪水冲刷和推力有着重要影响，通过构建不同桥台尺寸的木拱廊桥数值模型，研究其在洪水中的力学响应。设置三种桥台尺寸工况，分别为标准尺寸桥台、增大尺寸桥台和减小尺寸桥台。标准尺寸桥台按照实际工程中的常见尺寸设计；增大尺寸桥台在标准尺寸的基础上，将桥台的长度、宽度和高度分别增加20%；减小尺寸桥台则在标准尺寸的基础上，将桥台的长度、宽度和高度分别减小20%。在模拟过程中，设置相同的洪水工况，包括流速为3m/s、水位达到桥墩高度的2/3、流量为300m³/s。模拟结果显示，在相同洪水条件下，不同桥台尺寸的木拱廊桥抵抗洪水冲刷和推力的能力存在显著差异。标准尺寸桥台的木拱廊桥，在洪水冲击下，桥台底部受到一定的冲刷力，基础周围土体出现一定程度的松动。桥台能够承受拱架传递过来的水平推力和竖向荷载，桥体整体保持相对稳定。在此次模拟中，桥台底部的最大冲刷深度为0.5m，桥台的水平位移为6mm。增大尺寸桥台的木拱廊桥，由于桥台尺寸的增加，其抵抗洪水冲刷和推力的能力明显增强。更大的桥台尺寸使得其与地基的接触面积增大，能够更好地分散拱架传递过来的荷载，减小地基的应力。在相同洪水工况下，桥台底部的冲刷深度减小到了0.3m，桥台的水平位移减小到了4mm。这表明增大桥台尺寸能够有效提高木拱廊桥抵抗洪水冲刷和推力的能力，增强桥体的稳定性。减小尺寸桥台的木拱廊桥，由于桥台尺寸的减小，其抵抗洪水冲刷和推力的能力明显下降。较小的桥台尺寸使得其与地基的接触面积减小，无法有效分散拱架传递过来的荷载，导致地基应力增大。在洪水冲击下，桥台底部的冲刷深度增大到了0.8m，桥台的水平位移增大到了8mm。桥台出现明显的倾斜和下沉，桥体的整体稳定性受到严重威胁。通过对不同桥台尺寸木拱廊桥抗洪性能的研究，可以得出结论：桥台尺寸的变化对木拱廊桥抵抗洪水冲刷和推力的能力有着重要影响。增大桥台尺寸能够有效提高木拱廊桥的抗洪能力，增强桥体的稳定性；而减小桥台尺寸则会降低木拱廊桥的抗洪能力，增加桥体在洪水中的安全风险。4.3防护措施对木拱廊桥抗洪性能的影响4.3.1增设防护设施模拟为探究防护设施对木拱廊桥抗洪性能的影响，运用数值模拟技术，对增设挡水板、加固桥墩等防护设施后的木拱廊桥进行模拟分析。在模拟增设挡水板的情况时，考虑到挡水板的位置和高度对其防护效果的影响，分别设置了不同的工况。在桥体两侧靠近桥墩的位置设置了高度为1米的挡水板，模拟洪水冲击时，挡水板能够有效地阻挡部分洪水，改变水流方向，减少洪水对桥体的直接冲击力。通过模拟结果可以看到，在相同洪水工况下，增设挡水板后，桥体所受的冲击力明显减小，拱架和桥墩部位的应力也相应降低。拱架节点处的最大应力从增设挡水板前的8MPa降低到了6MPa，桥墩底部的应力也有所下降，这表明挡水板能够在一定程度上提高木拱廊桥的抗洪性能。在模拟加固桥墩的情况时，采用增加桥墩直径和在桥墩周围设置加固材料的方式进行模拟。将桥墩直径增加20%，并在桥墩周围包裹一层钢筋混凝土加固材料。模拟结果显示，加固后的桥墩抵抗洪水冲刷和推力的能力显著增强。在洪水冲击下，桥墩底部的冲刷深度明显减小，从加固前的0.6m减小到了0.3m。桥墩的水平位移也大幅降低，从加固前的8mm减小到了4mm。这说明加固桥墩能够有效地提高木拱廊桥的抗洪稳定性，减少桥墩在洪水中的损坏风险。4.3.2加固措施效果分析通过对不同加固方式下木拱廊桥抗洪性能的模拟分析，深入探讨了加固措施对提高桥梁抗洪性能的作用。对于增设挡水板的加固方式，其主要作用在于改变水流方向，减少洪水对桥体的直接冲击力。当洪水冲击挡水板时，挡水板将一部分水流引导至桥体两侧，避免了洪水直接冲击桥体，从而降低了桥体所受的冲击力。挡水板还能够在一定程度上阻挡洪水携带的杂物，减少杂物对桥体的撞击破坏。从模拟结果来看，增设挡水板后，桥体各构件的应力和应变明显减小，结构的稳定性得到提高。这表明增设挡水板是一种有效的提高木拱廊桥抗洪性能的措施，尤其适用于洪水冲击力较大的情况。加固桥墩的措施则主要通过增强桥墩的强度和稳定性来提高木拱廊桥的抗洪性能。增加桥墩直径能够增大桥墩与地基的接触面积，分散桥墩所承受的荷载，从而减小桥墩底部的应力。在桥墩周围设置钢筋混凝土加固材料，能够提高桥墩的抗冲刷能力和抗压强度，减少洪水对桥墩的侵蚀和破坏。模拟结果显示，加固后的桥墩在洪水作用下，其变形和位移明显减小，能够更好地支撑桥体，保证桥体的整体稳定性。这说明加固桥墩对于提高木拱廊桥的抗洪性能具有重要作用，特别是在洪水冲刷较为严重的情况下，能够有效地保护桥墩，确保桥梁的安全。通过对比不同加固方式的模拟结果，可以发现不同加固措施对木拱廊桥抗洪性能的影响各有侧重。在实际工程中，应