1. 引言
随着国家对海洋的不断开发,修建在海洋的建筑物或构筑物越来越多,例如港口、灯塔和机场等。修建在海洋的建(构)筑物的地基除了要承受自重荷载外,还要会受到波浪的作用,波浪荷载是一种循环荷载,在此循环荷载作用下,土颗粒会发生位移及重新排列,从而造成土的软化及地基的失稳。部分海洋土在循环荷载下并不是处于正常固结状态,而是以微超固结或超固结为主 [1] 。
目前,国内外的学者针对超固结土和动强度做出了以下的研究:孙磊 [2] 通过对恒定围压与变围压的对比,分析了循环围压与循环偏应力耦合下超固结软土的动力特性;Mayne等 [3] 通过对116组不同地区土样三轴剪切试验和直剪试验的数据进行分析总结,提出利用土体强度试验确定超固结比的计算公式;Azzouz,AS [4] 给出了正常固结和超固结土在双向对称不排水循环剪切下的试验结果;Imai G [5] 提出了超固结土在循环荷载下的本构关系;M Jiang [6] 研究了循环加载频率对饱和海相软土动力特性的影响研究,结果发现,循环频率对累积轴向应变和孔隙水压力有明显的影响,而频率对累积轴向应变的归一化平均孔隙压力影响不显著,频率是影响退化指数的一个重要因素。
以往的研究基本是针对海洋软土在正常固结状态下的动强度,对超固结土的动强度的研究较少,因此本文通过室内动三轴试验研究循环应力比和频率对不同超固结比的海洋软土动强度的影响。
2. 试验设备及方案
DDS-700微机控制电磁式振动三轴仪是完成本次试验的主要设备,是水利行业规范《土工试验规程》(SL237-2009)中振动三轴试验选定的设备,可用于测定各种土的动力特性指标。
本试验采用土样为重塑土样,试验制备根据《土工试验规程》(SL237-2009)具体规定执行,其基本物理性质指标见表1。取样深度为10~12 m,通过施加围压模拟其固结程度。正常固结时,土样在围压为100 kPa下等压固结24 h;超固结时让土样在大于100 kPa的围压下等压固结24 h,固结完毕后先卸载至100 kPa,再进行不排水分级加载动弹性模量试验,试验的波形采用正弦波。表2为本试验的分组情况。
3. 动应变-振次曲线分析
3.1. 循环应力比的影响
图1为不同循环应力比下动应变与振次的关系曲线。图中
为循环应力比,本文定义
,
,其中
为对应超固结比下土体不排水抗剪强度。本文中不排水抗剪强度定义为总应变达到5%时的抗剪强度。从图1中可以看出,不论是正常固结土还是超固结土,随着循环次数的增加,动应变
不断增大,相同循环次数下,循环应力比越大,动应变越大。同时也能看出,振次-动应变曲线随循环应力比的不同有较大的差异。当循环应力比小于0.5时,在加载初期,动应变增长较慢,随着循环次数的增加,当到达某一振次后,曲线出现转折点,此时动应变迅速增大,当动应变增长到某一值后趋于稳定,此时土样已经破坏;当循环应力比大于0.5时,即使在循环初期,动应变也呈较大幅度增长,在较小的循环次数下曲线就达到转折点,稳定后动应变值也较大。循环应力比越大,达到转折点所需的振次越小,即在较小的振次下土样就发生破坏且破坏时的动应变越大。这是因为,Larew和Leonards [7] 认为饱和软土存在临界循环应力比,其定义为土体破坏的最大循环应力比。通过图1可以看出,循环应力比大于0.5
![](Images/Table_Tmp.jpg)
Table 1. Physical-mechanical properties of soft soil
表1. 土样基本物理性质指标
(a) OCR = 1, f = 1 Hz
(b) OCR = 2, f = 1 Hz
(c) OCR = 4, f = 1 Hz
Figure 1. Relationship curves between dynamic strain and cyclic number at different cyclic stress ratio
图1. 不同循环应力比下动应变与振动周次关系曲线
和小于0.5的曲线有较大的不同,可以认为珠江海洋土的临界循环应力比在0.5左右。
3.2. 频率的影响
图2所示为
时,不同频率比下振动周次与动变关系曲线。从图中可以看出,曲线可以分为三个阶段,三个阶段的分隔点在两个应变转折处。第一个阶段,在振动初期,随着振次的增加,三条曲线的应变均有增长,但增长比较缓慢,三条曲线几乎重叠在一起,这表明循环荷载初期,频率对动应变影响不大。第二阶段,当应变到达第一个转折点时,三条曲线的应变开始快速增加,并开始出现明显的不同,主要表现在:相同振次下,频率越大,动应变越小,即频率越小的土样应变发展越大。这是因为频率越大,应力作用在试样上的时间越短,饱和软土的孔压来不及上升,应力还来不及传递给土颗粒就进
图2. 不同频率比下动应变与振动周次关系曲线
入下一轮加载,这会导致试样变形不完全,从而表现为动应变越小,土样在较大的循环次数下才发生破坏。第三阶段,当应变达到第二个转折点时,此时三条曲线趋于平缓,应变不再随着振次的增加而增加,土样已经破坏,且频率越大,破坏时的应变越大。
3.3. 超固结比的影响
图3为循环应力比为0.5时,不同超固结比下振动周次与动变关系曲线。从图中可以看出,不同超固结比的曲线有明显差异。当超固结比8时,动应变增长速度较慢,基本呈线性增长,且没有明显转折点,土样没有破坏;当超固结比为4和2时,土样的应变随循环次数的增加略大于超固结为8的土样,并且出现明显的转折点,且OCR = 2的曲线转折点后的动应变增长速度大于OCR = 4的曲线;当超固结比为1时,土样在初始阶段应变就最大,其随振次增长速度也最快,其破坏所需的振次也越少。相同振次下,超固结比越大,动应变越小,并且其达到破坏所需的振次越多,这是因为初期固结压力越大的土样越密实,颗粒间的空隙越小,变形所需轴力越大,即弹性模量越大。
4. 最大动应变分析
图4为相同频率下,试样在不同循环应力比下循环加载至1500次时的动应变随超固结比的变化曲线,从图中可以看出,在相同的频率和振次下,试样在不同循环应力比下产生的最大动应变随超固结比的增大而下降,采用幂函数拟合,公式如下:
(1)
表3为参数
和
的拟合结果。
从表中可以看出,参数
和
随着循环应力比的不同有很大差异。图5为参数
与循环应力比关系曲线,图中可以看出,参数
与循环应力比呈线性关系,可以用线性函数来表示,图6为参数a2与循环
![](//html.hanspub.org/file/4-2750335x34_hanspub.png)
Figure 3. Relationship curves between dynamic strain and cyclic number at different over consolidation ratio (f = 1 Hz,
)
图3. 不同超固结比下动应变与振动周次关系曲线(f = 1 Hz,
)
![](//html.hanspub.org/file/4-2750335x37_hanspub.png)
Figure 4. Relationship curves between dynamic axial strain (the 1500 th) and over-consolidation ratio at different cyclic stress ratio
图4. 不同循环应力比下第1500圈轴向动应变与超固结比关系曲线
![](Images/Table_Tmp.jpg)
Table 3. Parameter a1 and a2 fitting results
表3. 参数a1、a2拟合结果
![](//html.hanspub.org/file/4-2750335x41_hanspub.png)
Figure 5. Relationship curves of a1 vs. cyclic stress ratio
图5. 参数a1与循环应力比关系曲线
![](//html.hanspub.org/file/4-2750335x42_hanspub.png)
Figure 6. Relationship curves of a2 vs. cyclic stress ratio
图6. 参数a2与循环应力比关系曲线
应力比关系曲线,其关系可以用二次函数表示:
(2)
(3)
把式(2)和式(3)代入公式(1)可以得出最大动应力和超固结公式:
(4)
5. 振次-动应变模型建立与参数分析
为了明确珠江海洋土的振次与动应变的规律,通过观察图1~图3,即不同循环应力比、频率和超固结比下振次-动应变关系曲线,可以发现在第二个转折点前,振次-动应变曲线可以采用指数函数来表示:
(5)
表4为参数a和b的拟合结果。
从上表中可以看出,基本各
参数均在0.9以上,表明拟合度较高,该公式能较好的表示振次与动应变的关系。但是循环应力比为0.3和频率为3 Hz的曲线的
较低,为了找出其原因,观察其曲线发现这些曲线都近似呈线性关系,因而不再适用对数函数来表示,可以用线性函数来表示。
通过观察表4可以发现,参数a和b的值随循环应力比、超固结比和频率的不同而改变,对参数a和b进行数学分析,可得到相应的回归公式。
(1) 关于频率的回归公式
(6) (7)
将上式代入式(5),得到关于频率的振次-动应变关系方程:
![](Images/Table_Tmp.jpg)
Table 4. Parameter a and b fitting results
表4. 参数a、b拟合结果
(8)
(2)关于超固结比的回归公式
(9)
(10)
将上式代入式(5),得到关于超固结比的振次-动应变关系方程:
(11)
(3)关于循环应力比的回归公式
(12)
(13)
将上式代入式5),得到关于频率的振次-动应变关系方程:
(14)
6. 软化特性分析
动三轴试验中常用软化指数来表示试样动弹模的衰减规律。图7所示为应力控制的循环动三轴仪试验得到的珠江海洋土第1次与第N次循环的滞回曲线。根据Idriss [8] 的研究,定义软化指数
,其定义为第N次循环的土动弹性模量与第1次循环的土动弹性模量比值,考虑初始固结剪应力的影响,重新定义软化指数
如下:
(15)
上式中:
和
分别为第N次循环最大与最小偏应力;
和
分别为第N次循环最大与最小动应变;
和
分别为第1次循环最大与最小偏应力;
和
分别为第1次循环最大与最小动应变。
6.1. 循环应力比的影响
因为循环应力比为0.3时,应变缓慢增长,试样基本处于弹性状态,软化不明显,因此采用循环应力比大于0.3试样来研究。图8为不同循环应力比时软化指数与振次关系曲线。从图中可以发现,软化指数随着振次的增加而降低,这是因为随着振次的增加,应变累积不断增加,导致动弹性模量不断降低,从而使土体出现软化。在振动初期,软化指数随振次的增加明显降低,当软化指数降低到某一值后趋于稳定,此时土体动弹模不再降低,土样处于破坏状态。相同振次下,循环应力比越大,软化指数越低,即动弹模衰减越快,这是因为循环荷载越大,土体颗粒破坏越严重,土体破坏加剧。
6.2. 频率的影响
图9为不同频率时软化指数与振次关系曲线。从图中可以看出,无论是正常固结土还是超固结土,频率对软化指数的降低影响较小。但仍可看出,在相同振次下,频率越大,软化指数越大,即动弹性模量衰减越小,这是因为,这是因为频率越大,应力作用在试样上的时间越短,应力还来不及传递给土颗
![](//html.hanspub.org/file/4-2750335x70_hanspub.png)
Figure 7. Relationship curves between stresses and strains (the first and nth) under cyclic loading
图7. 循环荷载下第1次与第N次应力-应变曲线
(a) OCR = 1, f = 1 Hz
(b) OCR = 2, f = 1 Hz
(c) OCR = 4, f = 1 Hz
(d) OCR = 8, f = 1 Hz
Figure 8. Relationship curves between softening index and cyclic number at different cyclic stress ratio
图8. 不同循环应力比时软化指数与振次关系曲线
图9. 不同频率时软化指数与振次关系曲线
粒就进入下一轮加载,这会导致试样变形不完全,从而提高土的刚度。
图中还可看出,在振动初期,三条曲线间隔较小,此时频率对软化指数影响较小,随着振次的增加,三条曲线间隔逐渐增大,表明频率的影响逐渐增强,当振次达到一定数值时,三条曲线间隔逐渐减小,最后趋于相同,此时软化指数不再随着振次的增加而降低,趋于稳定。在循环振次达到1000次时,频率1、2和3的曲线软化指数都在0.2左右,此时试样动弹模衰减为原来的20%。
6.3. 超固结比的影响
图10为不同超固结比时软化指数与振次关系曲线。图中可以看出,相同振次下,超固结比越大,软化指数越高,一方面说明了超固结减缓了试样的软化;另一方面说明了超固结比越大的土样强度越高。这是因为初期固结压力越大的土样越密实,颗粒间的空隙越小,前期同样轴向压力下的变形越小,所以动模量越大。
7. 软化模型建立及参数分析
目前主要有两种经验软化模型,Idriss [8] 提出的指数表达式
以及Kazuya Yasuhara [9] 的半对数表达式
,其中d,D均为软化参数,但是通过观察图8~图10可以发现,珠江海洋土的软化指数与
呈曲线关系,难以用这两种经验模型来描述其关系。因此,本文通过研究发现可以采用王军等 [10] 提出的二次修正曲线描述软化指数与循环次数的关系:
(16)
其中A、B为软化参数,通过对试验结果得出的曲线进行回归分析得出。表5为参数A、B拟合结果。由于频率影响不大,忽略其作用。
由上表可以发现参数A、B随循环应力比的变化而变化,图11和图12分别为不同超固结比下循环应力比与参数A、B的变化曲线。从图中看出,参数A、B与循环应力比可以用二次曲线来描述其关系:
(17)
(18)
其中:A1,A2,A3,B1,B2,B3为曲线参数,通过试验来获得。表6为不同超固结比下试验参数拟合结果。
![](//html.hanspub.org/file/4-2750335x87_hanspub.png)
f = 1 Hz,
Figure 10.Relationship curves between softening index and cyclic number at over consolidation ratio
图10. 不同超固结比时软化指数与振次关系曲线
![](Images/Table_Tmp.jpg)
Table 5. Parameter A and B fitting results
表5. 参数A、B拟合结果
![](//html.hanspub.org/file/4-2750335x90_hanspub.png)
Figure 11. Relationship curves of A vs. cyclic stress ratio
图11. 参数A与循环应力比关系曲线
![](//html.hanspub.org/file/4-2750335x91_hanspub.png)
Figure 12. Relationship curves of B vs. cyclic stress ratio
图12. 参数B与循环应力比关系曲线
将公式(17)和(18)带入式(16),得到不同超固结比的珠江海洋土在循环荷载下的软化模型:
(19)
![](Images/Table_Tmp.jpg)
Table 6. Parameter A1, A2, A3, B1, B2 and B3 fitting results
表6. 参数A1,A2,A3,B1,B2,B3拟合结果
8. 动强度分析
动强度是指在一定动荷载往返作用次数N下土产生某一指定破坏应变
或满足某一破坏标准所需的动应力。如果破坏标准不同,相应的动强度就会不同,其数值会根据破坏标准的不同而不同,一般情况下取5%或10%应变。
8.1. 破坏准则分析
由前文振次-动应变曲线可以发现,当振次达到一定数值后,曲线会出现明显的转折点,此时的应变增幅急剧增大,并且在不同循环应力比下到达转折点所需的振次也不同,针对这种情况,选取5%应变和该应变转折点作为土样破坏标准,绘出其动强度曲线进行分析。
图13为不同频率下超固结珠江海洋土采用转折点和5%应变为破坏标准的动强度曲线。从图中可以看出,不论采用哪一种破坏标准,破坏振次越大,动强度越低。当破坏振次较小时,两条曲线比较接近,此时两种标准下的动强度相同,随着破坏振次的增加,两条曲线间隔逐渐增大,此时相同破坏振次下,5%标准的动强度明显大于转折点的动强度。此时若采用5%应变为破坏标准的话,的得到的动强度大于采用转折点的动强度,是不安全的,在此之前土体已经出现破坏。因此在岩土工程设计中,采用转折点为破坏标准是偏于安全的。
8.2. 超固结比的影响
图14为不同超固结比的动强度曲线。图中可以看出,超固结比越大,相同破坏振次下动强度越大。这是因为超固结比越大,土样越密实,土体强度越大,越不容易破坏,达到相同破坏振次所需的动强度越大。超固结比越大,曲线越平滑。f = 1时,三条曲线的斜率基本相同,动强度以同样的速率降低;f = 2时,破坏振次较低时,三条曲线间隔较大,此时超固结比对动强度的影响比较明显,随着破坏振次的增加,三条曲线间隔逐渐减小。
8.3. 频率的影响
图15为不同频率的动强度曲线。从图中可以看出频率越大,相同振次下动强度越大,这是因为频率越大,应力作用在试样上的时间越短,导致土样变形不完全,从而提高土体强度。但从图中也能看出,不论是正常固结土还是超固结土,三条曲线间隔都不大,表明频率对动强度有一定的影响,但影响不大。
9. 动强度模型的建立及参数分析
观察图13和图14可以发现,动强度-振次曲线呈线性关系,可以采用线性公式来表示其内在联系:
(20)
式中a和b为拟合参数。
表7所示为a和b参数拟合结果。表中可以发现,参数a和b与超固结比有关,与频率影响关系不大,可以忽略。
图16和图17分别为参数a和b与超固结比关系曲线。观察图可以发现,其参数与超固结比可以用二次线性描述其关系:
(21)
(22)
![](//html.hanspub.org/file/4-2750335x104_hanspub.png)
(a) f = 1 Hz(b) f = 2 Hz
Figure 14. Dynamic strength curve at different over consolidation ratio
图14. 不同超固结比的动强度曲线
![](//html.hanspub.org/file/4-2750335x106_hanspub.png)
(a) OCR = 1(b) OCR = 2
(c) OCR = 4
Figure 15. Dynamic strength curve at different frequency
图15. 不同频率的动强度曲线
得出参数
、
、
、
、
和
。将这
![](Images/Table_Tmp.jpg)
Table 7. Parameter a and b fitting results
表7. 参数a和b拟合结果
![](//html.hanspub.org/file/4-2750335x114_hanspub.png)
Figure 16. Relationship curves of a vs. over consolidation ratio
图16. 参数a和超固结比关系曲线
![](//html.hanspub.org/file/4-2750335x115_hanspub.png)
Figure 17. Relationship curves of b vs. over consolidation ratio
图17. 参数b和超固结比关系曲线
些参数带入式(20)得出:
(23)
10. 结论
对不同超固结比的珠江海洋软土进行不同循环应力比和不同频率的动强度试验,可得到以下的结论:
(1) 在循环荷载作用下,随着振动次数的增加,土体动应变不断增大,同时变化曲线存在两个应变转折点,加载初期,动应变缓慢增大,达到第一个转折点之后,总应变迅速增大,达到第二个转折点,动应变增长速度减缓,最后趋于稳定。第二个转折点前的曲线可以用指数函数进行拟合。
(2) 在同一振次下,循环应力比越大、动应变越大;超固结比越大、动应变越小;频率越大、动应变也越小。在相同的频率下,试样在不同循环应力比下产生的最大动应变随超固结比的增大而降低,且下降规律相同,可以采用幂函数进行拟合。
(3) 分析软化特性,发现软化指数随着振次的增加而降低,循环应力比越大,软化越快,超固结比越大,软化越慢,但频率对软化指数影响不大。建立软化模型并分析其参数,即可得到其关于振次的公式。
(4) 以第一个应变转折点和总应变5%为破坏标准进行分析,发现用转折点作为破坏标准更安全。以转折点应变作为破坏强度标准,通过绘出不同超固结比和频率的动强度曲线,可以发现,破坏振次越大,动强度越低,超固结比越大,动强度越大,频率越大,动强度也越大,但不明显。建立动强度模型并分析其参数,即可以得到关于超固结比的动强度经验公式。