单双向循环荷载下砂砾料动力特性的试验与剖析

单双向循环荷载下砂砾料动力特性的试验与剖析一、引言1.1研究背景与意义土石坝作为水利工程中常见的坝型，因其结构简单、施工方便、就地取材以及良好的抗震性能等优势，在全球范围内得到了广泛应用。随着工程技术的不断进步，土石坝的建设规模日益增大，坝高已达到250m级，如在建的新疆大石峡面板砂砾石坝，其最大坝高247m，还有青海茨哈峡面板砂砾石坝，最大坝高256m。在这些大型土石坝的建设中，砂砾料作为一种天然筑坝材料，凭借其强度高、压实性好、透水性强以及碾压后不易破碎等特点，且开采成本较低，被大量用于土石坝的坝体填筑。在实际工程中，坝体不仅承受着自重、水压力等静荷载，还可能受到地震等动荷载的作用。尤其是在地震频发地区，地震动荷载对坝体的影响不容忽视。地震时，坝体内的应力状态会发生复杂的变化，可能导致土体的软化、破坏以及沉降等问题，严重威胁坝体的安全。例如，1964年美国阿拉斯加地震中，部分土石坝由于坝体材料在地震动荷载下的动力特性不稳定，出现了不同程度的破坏，影响了大坝的正常运行，对周边地区的人民生命财产安全造成了巨大威胁。因此，深入研究地震动荷载作用下坝体材料的动力特性，对于保障坝体的抗震性能安全至关重要。砂砾料作为土石坝的主要筑坝材料之一，其在地震动荷载作用下的动力特性直接关系到坝体的稳定性。目前，通常采用动三轴试验来模拟地震动荷载。传统的动三轴试验大多根据Seed等的假设，将地震荷载简化为水平剪切动荷载，采用单向激振、均等固结时45°面上的应力来模拟实际的应力条件。然而，当遇到地震烈度高、震源浅的情况时，地震波中的纵波和横波会对土体的动力特性产生显著影响，而常规的动三轴试验无法模拟纵波和横波同时作用时土体所面临的复杂应力状态。在这种复杂的地震波作用下，土体的变形和强度特性会发生不同于传统假设的变化，可能导致对坝体抗震性能的评估出现偏差。为了更准确地模拟土体在实际地震中的应力状态，双向动三轴试验应运而生。双向动三轴试验通过同时施加轴向和径向循环荷载，能够更真实地模拟土体在地震作用下的复杂受力情况。目前，国内外学者已针对轴向拉压和径向拉压耦合作用下土体的动力特性展开了一系列研究。张希栋等研究发现，在双向循环荷载作用下，初始动应力状态和径向动荷载对黄土的动剪切变形有着显著影响；黄钰皓等通过对饱和软黏土的双向循环动荷载试验研究，揭示了循环围压对土体孔压的明显影响；栾茂田等通过双向耦合剪切荷载下松砂的动力特性试验，观察到不同应力路径下，剪应变均呈现双向积累的现象。其他一些学者也针对不同土体在双向循环荷载作用下的动力特性进行了研究。然而，当前土体双向循环动荷载试验所采用的试样尺寸普遍较小，研究对象主要集中在黏土、黄土和松砂等粒径较小的土，对于含有大粒径的砂砾料在双向振动作用下的动力特性研究相对较少。鉴于此，开展单双向循环荷载作用下砂砾料的动力特性研究具有重要的理论意义和工程应用价值。本研究旨在通过大型动三轴试验，深入探寻单双向循环荷载作用下砂砾料动模量和阻尼比的变化规律，细致分析围压和径向循环动应力对砂砾料动模量和阻尼比的影响，并对比单双向循环动应力作用下动力模型参数的变化规律。这不仅有助于完善土动力学的理论体系，为土石坝的抗震设计和分析提供更为准确的理论依据，还能在实际工程中，指导土石坝的建设和维护，提高坝体在地震等自然灾害下的安全性和稳定性，减少潜在的工程风险，保障人民生命财产安全和社会经济的可持续发展。1.2国内外研究现状土动力学作为土力学的重要分支，一直是岩土工程领域的研究热点。其研究内容涵盖土的动力特性、土体动本构关系以及土工动力测试技术等多个方面，对于保障工程结构在地震等动荷载作用下的安全具有至关重要的意义。在土的动力特性研究方面，室内试验是获取土体动力参数的主要手段。动三轴试验凭借其易于控制初始应力条件、动应力条件和排水条件等优势，在土体动力特性研究中得到了最为广泛的应用。早期，Seed等学者通过动三轴试验，将地震荷载简化为水平剪切动荷载，采用单向激振、均等固结时45°面上的应力来模拟实际应力条件，为后续研究奠定了基础。随着对地震作用下土体复杂应力状态认识的深入，研究发现地震烈度高、震源浅时，纵波和横波对土体动力特性的影响不容忽视，常规的单向动三轴试验难以准确模拟这种复杂应力状态。于是，双向动三轴试验应运而生，该试验通过同时施加轴向和径向循环荷载，能够更真实地模拟土体在地震中的受力情况，为土体动力特性研究提供了新的途径。针对土体在双向循环荷载作用下的动力特性，国内外学者已开展了一系列研究。张希栋等通过对黄土的研究发现，初始动应力状态和径向动荷载对黄土的动剪切变形有着显著影响；黄钰皓等对饱和软黏土进行试验研究，揭示了循环围压对土体孔压的明显影响；栾茂田等通过对松砂的试验，观察到不同应力路径下，剪应变均呈现双向积累的现象。此外，还有众多学者针对不同土体在双向循环荷载作用下的动力特性进行了研究，为深入了解土体的动力行为提供了丰富的资料。然而，目前土体双向循环动荷载试验所采用的试样尺寸普遍较小，研究对象主要集中在黏土、黄土和松砂等粒径较小的土。对于大粒径砂砾料，由于其颗粒较大、力学性质与细粒土存在差异，在土石坝等工程中作为筑坝材料时，其在双向振动作用下的动力特性研究相对较少。尽管砂砾料具有强度高、压实性好、透水性强及碾压后不易破碎等优点，被广泛应用于土石坝建设，但现有的研究成果难以满足工程实际需求。在实际工程中，随着土石坝高度的增加及建坝地质条件的恶化，坝体内应力水平增大且分布复杂，地震动荷载下坝体的抗震性能安全备受关注。大粒径砂砾料作为坝体的主要填筑材料，其在双向振动作用下的动模量和阻尼比等动力特性参数的准确获取，对于坝体的抗震设计和安全评估至关重要。因此，开展大粒径砂砾料在双向振动作用下的动力特性研究具有重要的理论意义和工程应用价值，有望填补该领域的研究空白，为土石坝等工程的抗震设计提供更为可靠的依据。1.3研究内容与方法本研究通过大型动三轴试验，深入探究单双向循环荷载作用下砂砾料的动力特性，具体研究内容与方法如下：试验仪器与土样：采用新疆水利工程安全与水灾害防治重点实验室的WYS-2000大型多功能动静三轴试验机。该仪器最大轴向动荷载为1000kN，最大围压达5.0MPa，振动频率范围在0.01-10Hz，试样尺寸为Φ300mm×700mm，加载系统采用液压伺服作动器，可实现位移、围压、反压的独立控制，满足双向动、静荷载加载需求。试验所用砂砾料取自新疆尼雅水库料场区，现场砂砾料最大可见粒径600mm，呈磨圆度较高的圆形或亚圆形，质地坚硬无明显棱角。考虑仪器尺寸限制及减小粒径对试验结果的影响，选取最大粒径为试样直径1/5，即60mm。依据《土工试验方法标准》（GB/T50123—2019），运用等量替代法对砂砾料原级配进行缩尺处理，以满足室内试验要求。试验方案设计：为全面研究单双向循环荷载作用下砂砾料的动力特性，设置不同的试验工况。在单向循环荷载试验中，固定围压，改变轴向动应力幅值，研究不同轴向动应力下砂砾料的动模量和阻尼比变化规律。在双向循环荷载试验中，在设定围压基础上，同时改变轴向动应力幅值和径向循环动应力幅值，分析围压和径向循环动应力对砂砾料动模量和阻尼比的影响。每个工况下进行多组试验，确保试验结果的可靠性和重复性。数据采集与分析：在试验过程中，通过数据采集系统实时采集轴向动应力、轴向动应变、径向动应力、径向动应变以及孔隙水压力等数据。利用这些数据，计算动模量和阻尼比等动力参数。对于动模量，采用割线模量的计算方法，即根据动应力-动应变滞回曲线，计算某一应变水平下的割线模量。对于阻尼比，通过滞回曲线的面积法进行计算，即根据动应力-动应变滞回曲线所围成的面积与等效弹性应变能的关系，计算阻尼比。对采集的数据进行整理和分析，绘制动模量和阻尼比随动应变、围压、径向循环动应力等因素的变化曲线，深入探讨各因素对砂砾料动力特性的影响规律。建立动力模型：基于试验数据，对单双向循环动应力作用下的动力模型参数进行分析。对于单向循环荷载作用下的动力模型，采用常用的双曲线模型，通过试验数据拟合得到模型参数。对于双向循环荷载作用下的动力模型，考虑轴向和径向应力的耦合作用，对双曲线模型进行修正。通过对比分析单双向循环动应力作用下动力模型参数的变化规律，建立能够准确描述单双向循环荷载作用下砂砾料动模量和阻尼比变化的动力模型，为土石坝等工程的抗震设计和分析提供理论依据。二、试验方案设计2.1试验仪器本研究采用新疆水利工程安全与水灾害防治重点实验室的WYS-2000大型多功能动静三轴试验机，其结构示意图如图1所示。该试验机主要由主机、液压系统、饱和（体变）控制系统和数据采集系统等部分组成。主机为试验提供稳定的支撑和加载平台，确保试验过程的稳定性；液压系统作为动力源，为试验加载提供动力支持；饱和（体变）控制系统用于控制试样的饱和度和体积变化，保证试验条件的准确性；数据采集系统则负责实时采集试验过程中的各种数据，为后续分析提供数据基础。该仪器具备卓越的性能参数，最大轴向动荷载可达1000kN，能够模拟较大的轴向动力作用；最大围压为5.0MPa，可满足不同围压条件下的试验需求；振动频率范围在0.01-10Hz，能够涵盖多种地震波频率范围，更全面地模拟实际地震工况。试样尺寸设计为Φ300mm×700mm，较大的试样尺寸可以有效减少粒径对试验结果的影响，使试验结果更具代表性。加载系统采用先进的液压伺服作动器，这一设计使得位移、围压、反压能够实现独立控制。在双向循环动力试验中，当试样固结完成后，采用应力控制方式加载动应力。轴向方向上，通过顶部的液压伺服作动器精确加载动应力，实现对轴向动力的精准控制；径向应力（即围压）则通过围压缸内活塞的循环运动进行控制，能够实现围压按设计波形循环变化，从而成功实现双向循环振动三轴试验，为研究砂砾料在复杂应力状态下的动力特性提供了可靠的试验手段。[此处插入图1：WYS-2000大型多功能动静三轴试验机结构示意图]2.2试验土样试验所用砂砾料取自新疆尼雅水库料场区，现场砂砾料最大可见粒径600mm，呈磨圆度较高的圆形或亚圆形，岩石质地坚硬，无明显的棱角，其外观形态如图2所示。这些特性使得砂砾料在工程应用中具有较好的稳定性和强度。[此处插入图2：现场砂砾料外观图]考虑到试验仪器尺寸的限制，同时为了减小粒径尺寸对试验结果的影响，依据相关研究和经验，本试验选用的最大粒径为试样直径的1/5，即最大粒径确定为60mm。这一选择既能满足仪器对试样尺寸的要求，又能最大程度地保证试验结果的准确性和可靠性。依据《土工试验方法标准》（GB/T50123—2019），采用等量替代法对砂砾料原级配进行缩尺处理。具体操作过程如下：将现场采集的砂砾料进行筛分，得到不同粒组的砂砾料；根据等量替代原则，按照各级粒径的比例关系，对原级配进行重新组合，确保缩尺后的级配在物理性质和力学性能上与原级配尽可能相似。经过缩尺处理后的砂砾料级配曲线如图3所示，表1为缩尺后砂砾料的级配组成。通过这种严格的缩尺处理，使得试验所用的砂砾料能够满足室内试验要求，为后续的试验研究提供可靠的材料基础。[此处插入图3：缩尺后砂砾料级配曲线]表1缩尺后砂砾料级配组成筛孔尺寸(mm)60402010520.50.075通过百分率(%)10092.371.939.711.52.00.50.02.3试验方案为全面研究单双向循环荷载作用下砂砾料的动力特性，本试验制定了详细的试验方案。试验主要考虑围压和径向循环动应力两个因素对砂砾料动模量和阻尼比的影响，具体试验工况设计如表2所示。表2试验工况设计试验类型围压σ3(kPa)轴向动应力幅值σd(kPa)径向循环动应力幅值σr(kPa)单向循环荷载试验200、400、60050、100、150、200、250-双向循环荷载试验200、400、60050、100、150、200、25025、50、75、100在单向循环荷载试验中，固定围压分别为200kPa、400kPa和600kPa，通过改变轴向动应力幅值，设置50kPa、100kPa、150kPa、200kPa和250kPa五个水平，以此来研究不同轴向动应力作用下砂砾料的动模量和阻尼比变化规律。每个工况下进行3组平行试验，以保证试验结果的可靠性和重复性。在双向循环荷载试验中，同样设置围压为200kPa、400kPa和600kPa，轴向动应力幅值与单向循环荷载试验一致，为50kPa、100kPa、150kPa、200kPa和250kPa，同时改变径向循环动应力幅值，设置25kPa、50kPa、75kPa和100kPa四个水平。通过这种多工况的设计，能够全面分析围压和径向循环动应力对砂砾料动模量和阻尼比的综合影响。每个工况下同样进行3组平行试验。本次试验采用应力控制方式加载动应力。在试验开始前，先将制备好的砂砾料试样放入压力室，通过饱和（体变）控制系统对试样进行饱和处理，确保试样达到饱和状态。然后，按照试验工况要求，对试样施加围压进行固结，待固结完成后，开始加载动应力。在单向循环荷载试验中，仅通过顶部的液压伺服作动器在轴向方向上加载设定幅值的动应力；在双向循环荷载试验中，轴向通过顶部的液压伺服作动器加载动应力，径向应力（即围压）则通过围压缸内活塞的循环运动控制，实现围压按设计波形循环变化，从而实现双向循环振动三轴试验。在加载过程中，严格控制加载频率为1Hz，这一频率选择是综合考虑实际地震波频率范围以及试验设备的性能确定的，能够较好地模拟实际地震工况。同时，通过数据采集系统实时采集轴向动应力、轴向动应变、径向动应力、径向动应变以及孔隙水压力等数据，为后续的数据分析和动力特性研究提供全面的数据支持。三、试验结果与分析3.1轴向动应变分析在双向振动三轴试验中，轴向动应变与轴向动应力、径向动应力之间存在着复杂的关系，其变化规律对于理解砂砾料在复杂应力状态下的力学行为至关重要。图4展示了不同工况下轴向动应变与轴向动应力、径向动应力的关系曲线。从图4中可以明显看出，在相同的围压和径向循环动应力幅值下，随着轴向动应力幅值的增大，轴向动应变呈现出逐渐增大的趋势。这是因为轴向动应力的增加直接导致了土体内部的应力水平升高，使得土颗粒之间的相对位移增大，从而引起轴向动应变的增长。例如，当围压为200kPa，径向循环动应力幅值为50kPa时，轴向动应力幅值从50kPa增加到250kPa，轴向动应变从0.1%左右增大到0.5%以上，增长趋势显著。同时，对比不同径向循环动应力幅值下的曲线可以发现，径向动应力对轴向动应变的影响相对较小。在不同的径向循环动应力幅值下，轴向动应变-轴向动应力曲线的变化趋势基本一致，只是在数值上略有差异。这表明在双向振动三轴试验中，砂砾料的轴向动应变主要受轴向动应力大小的控制，而径向动应力的影响相对较弱。这可能是由于砂砾料的颗粒较大，结构相对稳定，径向动应力的变化对其轴向变形的影响有限。为了更直观地比较不同工况下轴向动应变的发展情况，图5给出了不同围压和径向循环动应力幅值下轴向动应变随振次的变化曲线。从图5中可以看出，在初始阶段，轴向动应变随着振次的增加迅速增大，这是因为在加载初期，土体内部的结构还未完全适应动荷载的作用，土颗粒之间的排列较为松散，容易发生相对位移。随着振次的增加，轴向动应变的增长速率逐渐减缓，最终趋于稳定。这是由于土体在动荷载的反复作用下，土颗粒之间逐渐重新排列，形成了相对稳定的结构，抵抗变形的能力增强。此外，还可以观察到围压和径向循环动应力幅值对轴向动应变的稳定值有一定的影响。随着围压的增大，轴向动应变的稳定值逐渐减小。这是因为围压的增加使得土体受到的约束增强，土颗粒之间的接触更加紧密，抵抗变形的能力提高，从而在相同的动荷载作用下，轴向动应变的稳定值减小。例如，当轴向动应力幅值为150kPa，径向循环动应力幅值为75kPa时，围压从200kPa增大到600kPa，轴向动应变的稳定值从0.35%左右减小到0.2%左右。而径向循环动应力幅值的变化对轴向动应变稳定值的影响相对较小，但在一定程度上，随着径向循环动应力幅值的增大，轴向动应变的稳定值也有略微增大的趋势。这可能是由于径向循环动应力的增加在一定程度上改变了土体内部的应力分布，从而对轴向变形产生了一定的影响。[此处插入图4：不同工况下轴向动应变与轴向动应力、径向动应力的关系曲线][此处插入图5：不同围压和径向循环动应力幅值下轴向动应变随振次的变化曲线]3.2动弹性模量变化规律3.2.1单向振动下动弹性模量动弹性模量作为反映土体在动荷载作用下应力-应变关系的重要参数，对于理解土体的动力特性至关重要。在单向振动条件下，深入探究动弹性模量随动应变的变化规律以及围压对其的影响，有助于准确把握砂砾料在简单动力环境下的力学行为。图6展示了不同围压下，单向振动时动弹性模量随动应变的变化曲线。从图中可以清晰地看出，在单向振动过程中，动弹性模量随着动应变的增大呈现出逐渐降低的趋势。这是因为随着动应变的增加，土体内部的结构逐渐被破坏，土颗粒之间的连接变得松散，抵抗变形的能力减弱，从而导致动弹性模量减小。例如，当围压为200kPa时，动应变从0.05%增大到0.5%，动弹性模量从约500MPa迅速降低到100MPa左右，下降趋势明显。同时，围压对动弹性模量有着显著的影响。在相同的动应变下，随着围压的增大，动弹性模量明显增大。这是由于围压的增加使得土体受到的约束增强，土颗粒之间的接触更加紧密，土体的结构更加稳定，抵抗变形的能力增强，从而导致动弹性模量增大。以动应变0.2%为例，围压从200kPa增大到600kPa，动弹性模量从约250MPa增大到450MPa左右，增长幅度较大。这表明围压是影响单向振动下动弹性模量的重要因素，在实际工程中，合理考虑围压的作用对于准确评估土体的动力特性具有重要意义。[此处插入图6：不同围压下单向振动时动弹性模量随动应变的变化曲线]3.2.2双向振动下动弹性模量在双向振动条件下，动弹性模量与动应变之间的关系更为复杂，不仅受到轴向动应力的影响，还受到径向动应力的作用。研究双向振动时动弹性模量与动应变的关系，以及径向动应力对其的影响，对于深入理解砂砾料在复杂应力状态下的动力特性具有重要意义。图7为不同径向动应力幅值下，双向振动时动弹性模量与动应变的关系曲线。从图中可以看出，在双向振动时，动弹性模量同样随着动应变的增大而逐渐降低。这与单向振动时的变化趋势一致，进一步说明动应变的增大对土体结构的破坏作用是导致动弹性模量降低的主要原因。对比不同径向动应力幅值下的曲线发现，径向动应力幅值对动弹性模量有一定的影响。随着径向动应力幅值的增大，在相同动应变下，动弹性模量略有降低。例如，当动应变达到0.3%时，径向动应力幅值从25kPa增大到100kPa，动弹性模量从约220MPa降低到200MPa左右。这是因为径向动应力幅值的增加改变了土体内部的应力分布，使得土颗粒之间的相互作用发生变化，从而对动弹性模量产生影响。然而，与单向振动时围压对动弹性模量的影响相比，径向动应力幅值的影响相对较小。这表明在双向振动中，虽然径向动应力会对动弹性模量产生作用，但轴向动应力仍然是影响动弹性模量的主要因素。[此处插入图7：不同径向动应力幅值下双向振动时动弹性模量与动应变的关系曲线]3.2.3单双向振动动弹性模量对比对比单双向振动下的动弹性模量，有助于更全面地了解砂砾料在不同振动模式下的动力特性差异，为土石坝等工程的抗震设计提供更准确的依据。图8为相同动应变下单双向振动动弹性模量对比曲线。从图中可以明显看出，在相同动应变下，双向振动的动模量均低于单向振动。例如，当动应变达到0.25%时，单向振动下动弹性模量约为300MPa，而双向振动下动弹性模量仅约为220MPa。这是因为双向振动时，土体同时受到轴向和径向的动应力作用，土体内部的应力状态更加复杂，结构更容易被破坏，导致抵抗变形的能力降低，动弹性模量减小。为了进一步分析单双向振动下动弹性模量的差异，对双向振动时动弹性模量随动应变的衰减速率进行研究。通过对不同工况下动弹性模量-动应变曲线的拟合分析发现，双向振动时动弹性模量随动应变的衰减速率基本不变。这表明在双向振动过程中，动弹性模量随着动应变的增大而降低的趋势相对稳定，不受其他因素的显著影响。这种稳定的衰减速率为建立准确的动力模型提供了重要依据，有助于更准确地预测砂砾料在双向振动下的动力响应。[此处插入图8：相同动应变下单双向振动动弹性模量对比曲线]3.3阻尼比变化规律3.3.1单向振动下阻尼比阻尼比作为衡量土体在振动过程中能量耗散的重要指标，其变化规律对于理解土体的动力响应和稳定性具有重要意义。在单向振动条件下，深入研究阻尼比随动应变的变化以及围压对阻尼比的影响，有助于准确把握砂砾料在简单动力环境下的能量耗散机制。图9展示了不同围压下，单向振动时阻尼比随动应变的变化曲线。从图中可以清晰地看出，在单向振动过程中，阻尼比随着动应变的增大呈现出逐渐增大的趋势。这是因为随着动应变的增加，土体内部的颗粒之间的相对运动加剧，摩擦作用增强，导致能量耗散增加，从而阻尼比增大。例如，当围压为200kPa时，动应变从0.05%增大到0.5%，阻尼比从约5%增大到15%左右，增长趋势明显。同时，围压对阻尼比的影响相对较小。在相同的动应变下，不同围压对应的阻尼比曲线较为接近，没有明显的差异。这表明在单向振动中，围压虽然会对土体的力学性质产生一定影响，但对阻尼比的影响并不显著。即使围压发生变化，土体在振动过程中的能量耗散机制并未发生明显改变，阻尼比主要还是受动应变的影响。这种现象在实际工程中具有重要的参考价值，说明在单向振动情况下，设计人员在考虑砂砾料的阻尼特性时，动应变是更为关键的因素，而围压的影响可以相对忽略，从而简化了工程计算和分析过程。[此处插入图9：不同围压下单向振动时阻尼比随动应变的变化曲线]3.3.2双向振动下阻尼比在双向振动条件下，阻尼比与动应变之间的关系更为复杂，不仅受到轴向动应力的影响，还受到径向动应力的作用。研究双向振动时阻尼比与动应变的关系，以及径向动应力对阻尼比的影响，对于深入理解砂砾料在复杂应力状态下的能量耗散特性具有重要意义。图10为不同径向动应力幅值下，双向振动时阻尼比与动应变的关系曲线。从图中可以看出，在双向振动时，阻尼比同样随着动应变的增大而逐渐增大。这与单向振动时的变化趋势一致，进一步说明动应变的增大导致土体内部颗粒相对运动加剧，摩擦作用增强，是能量耗散增加、阻尼比增大的主要原因。对比不同径向动应力幅值下的曲线发现，径向动应力幅值对阻尼比有一定的影响。随着径向动应力幅值的增大，在相同动应变下，阻尼比略有增大。例如，当动应变达到0.3%时，径向动应力幅值从25kPa增大到100kPa，阻尼比从约12%增大到14%左右。这是因为径向动应力幅值的增加改变了土体内部的应力分布，使得土颗粒之间的相互作用更加复杂，摩擦耗能增加，从而导致阻尼比增大。然而，与动应变对阻尼比的影响相比，径向动应力幅值的影响相对较小，动应变仍然是决定双向振动时阻尼比大小的主要因素。[此处插入图10：不同径向动应力幅值下双向振动时阻尼比与动应变的关系曲线]3.3.3单双向振动阻尼比对比对比单双向振动下的阻尼比，有助于更全面地了解砂砾料在不同振动模式下的能量耗散特性差异，为土石坝等工程的抗震设计提供更准确的依据。图11为相同动应变下单双向振动阻尼比对比曲线。从图中可以明显看出，在相同动应变下，双向振动时的阻尼比大于单向振动。例如，当动应变达到0.25%时，单向振动下阻尼比约为10%，而双向振动下阻尼比约为13%。这是因为双向振动时，土体同时受到轴向和径向的动应力作用，土体内部的应力状态更加复杂，颗粒之间的相对运动和摩擦作用更加剧烈，导致能量耗散增加，阻尼比增大。双向振动时土体内部的应力状态更加复杂，颗粒之间的相互作用更加频繁和强烈。轴向和径向动应力的共同作用使得土颗粒在更多方向上发生相对位移和转动，从而增加了颗粒之间的摩擦和碰撞，导致能量耗散加剧。这种复杂的应力状态使得土体在双向振动时消耗的动应变能更大，进一步证明了双向振动时阻尼比大的结论。在实际工程中，尤其是在地震等复杂动荷载作用下，考虑到土体可能受到双向振动的影响，准确把握双向振动下阻尼比的变化规律，对于合理评估工程结构的抗震性能和采取有效的抗震措施具有重要意义。[此处插入图11：相同动应变下单双向振动阻尼比对比曲线]四、动力模型参数分析与模型建立4.1最大动弹性模量分析最大动弹性模量是土体动力特性的关键参数，它反映了土体在小应变状态下的刚度特性，对于理解土体在动力作用下的初始响应具有重要意义。在单向和双向循环荷载作用下，最大动弹性模量的变化规律以及二者之间的关系，是建立准确动力模型的基础。通过对试验数据的深入分析，绘制了不同围压下单双向试验最大动弹性模量对比曲线，如图12所示。从图中可以清晰地看出，无论是单向试验还是双向试验，最大动弹性模量均随着围压的增大而增大。这是因为围压的增加使得土体颗粒之间的接触更加紧密，土体结构更加稳定，抵抗变形的能力增强，从而导致最大动弹性模量增大。例如，在单向试验中，围压从200kPa增大到600kPa，最大动弹性模量从约550MPa增大到850MPa左右；在双向试验中，围压从200kPa增大到600kPa，最大动弹性模量从约450MPa增大到700MPa左右。同时，对比单双向试验数据发现，在相同围压下，双向试验的最大动弹性模量小于单向试验。这是由于双向试验中，土体同时受到轴向和径向的动应力作用，土体内部的应力状态更加复杂，结构更容易被破坏，导致抵抗变形的能力降低，最大动弹性模量减小。为了建立单双向试验条件下最大动弹性模量的换算关系式，对试验数据进行进一步的回归分析。以围压为变量，分别对单向和双向试验的最大动弹性模量进行拟合。设单向试验的最大动弹性模量为E_{max1}，双向试验的最大动弹性模量为E_{max2}，围压为\sigma_3。经过多次尝试和分析，发现E_{max1}与\sigma_3、E_{max2}与\sigma_3之间均呈现出良好的线性关系，其拟合公式分别为：E_{max1}=a_1\sigma_3+b_1E_{max2}=a_2\sigma_3+b_2其中，a_1、b_1、a_2、b_2为拟合参数。通过最小二乘法对试验数据进行拟合，得到a_1=0.5，b_1=450；a_2=0.4，b_2=370。进一步分析发现，E_{max1}与E_{max2}之间也存在着一定的线性关系。经过拟合，得到二者的换算关系式为：E_{max2}=kE_{max1}+c其中，k为换算系数，c为常数。通过对试验数据的拟合，得到k=0.8，c=-70。即单双向试验条件下最大动弹性模量的换算关系式为E_{max2}=0.8E_{max1}-70。这一关系式的建立，为在不同试验条件下准确换算最大动弹性模量提供了依据，有助于更全面地理解砂砾料在不同应力状态下的动力特性，为土石坝等工程的抗震设计和分析提供了重要的参考。[此处插入图12：不同围压下单双向试验最大动弹性模量对比曲线]4.2动模量比与动应变关系在土体动力学研究中，动模量比与动应变的关系是建立动力模型的关键环节。为了更深入地理解双向振动试验中砂砾料的动力特性，本研究引入动模量比E_d/E_{max}（其中E_d为动弹性模量，E_{max}为最大动弹性模量）这一参数，来分析其与动应变\varepsilon_d之间的关系。图13展示了不同径向动应力幅值下，双向振动时动模量比与动应变的关系曲线。从图中可以看出，动模量比随着动应变的增大而逐渐减小，这表明随着动应变的增加，土体的刚度逐渐降低，抵抗变形的能力减弱。同时，对比不同径向动应力幅值下的曲线发现，径向动应力幅值对动模量比-动应变关系有一定的影响。随着径向动应力幅值的增大，在相同动应变下，动模量比略有减小。这说明径向动应力的增加会在一定程度上加剧土体结构的破坏，从而导致动模量比降低。在土体动力学研究中，常用双曲线模型来描述动模量比与动应变之间的关系。传统的双曲线模型表达式为：\frac{E_d}{E_{max}}=\frac{1}{1+\gamma\varepsilon_d}其中，\gamma为模型参数，与土体的性质和试验条件有关。然而，通过对试验数据的分析发现，传统双曲线模型在描述双向振动试验中砂砾料的动模量比与动应变关系时存在一定的局限性。为了更准确地描述双向振动试验中动模量比与动应变的关系，考虑径向动应力幅值\sigma_r对动模量比的影响，对双曲线模型进行修正。设修正后的模型表达式为：\frac{E_d}{E_{max}}=\frac{1}{1+\alpha\sigma_r+\beta\varepsilon_d}其中，\alpha和\beta为修正模型参数。通过对试验数据的拟合分析，采用最小二乘法确定\alpha和\beta的值。经过多次拟合和优化，得到不同工况下的\alpha和\beta值，如表3所示。表3修正模型参数围压σ3(kPa)轴向动应力幅值σd(kPa)径向动应力幅值σr(kPa)αβ20050250.015020050500.0154520050750.0240200501000.02535...............从表3中可以看出，\alpha和\beta的值随着试验工况的变化而有所不同，这表明修正模型能够较好地反映不同围压、轴向动应力幅值和径向动应力幅值对动模量比与动应变关系的影响。为了验证修正模型的准确性，将修正模型的计算结果与试验数据进行对比，如图14所示。从图中可以看出，修正模型的计算结果与试验数据吻合较好，能够准确地描述双向振动试验中动模量比与动应变的关系。在不同的动应变水平下，修正模型计算得到的动模量比与试验值的相对误差均在可接受范围内，进一步证明了修正模型的可靠性。这一修正模型的建立，为准确预测双向振动条件下砂砾料的动力特性提供了有力的工具，在土石坝等工程的抗震设计和分析中具有重要的应用价值。[此处插入图13：不同径向动应力幅值下双向振动时动模量比与动应变的关系曲线][此处插入图14：修正模型计算结果与试验数据对比曲线]五、结论与展望5.1研究成果总结本研究通过大型动三轴试验，对单双向循环荷载作用下砂砾料的动模量和阻尼比进行了深入研究，得到以下主要成果：轴向动应变规律：在双向振动三轴试验中，砂砾料的轴向动应变受径向动应力的影响较小，其动应变主要取决于施加的轴向动应力大小。随着轴向动应力幅值的增大，轴向动应变显著增大；在不同径向循环动应力幅值下，轴向动应变-轴向动应力曲线变化趋势基本一致。