粗粒料颗粒破碎变化规律大型三轴试验研究

刘萌成; 孟锋; 王洋洋

doi:10.11779/CJGE202003018

粗粒料颗粒破碎变化规律大型三轴试验研究 English Version

1.
浙江工业大学建筑工程学院,浙江　杭州　310023
2.
浙江省交通运输科学研究院,浙江　杭州　311305

基金项目:

国家自然科学基金面上项目 51378472

国家自然科学基金面上项目 51778584

详细信息

作者简介:
刘萌成（1972— ）,男,博士,教授,主要从事交通岩土与特殊土土力学方面的教学与科研工作。E-mail：mcliu2002@163.com

中图分类号: TU433
计量
- 文章访问数: 346
- HTML全文浏览量: 28
- PDF下载量: 212
出版历程
- 收稿日期: 2019-03-27
- 网络出版日期: 2022-12-07
- 刊出日期: 2020-02-29

Evolution of particle crushing of coarse-grained materials in large-scale triaxial tests

1.
College of Civil Engineering and Architecture, Zhejiang University of Technology, Hangzhou 310023, China
2.
Zhejiang Scientific Research Institute of Transport, Hangzhou 311305, China

摘要

摘要: 基于土的颗粒破碎渐进演化分析,提出了土的颗粒破碎4项基本性质,以及在长期/反复加载过程中土的颗粒破碎3个特征状态。定义了土的全局与局部破碎指标,用以区分破碎的中止和终止状态,同时分析给出了（修正）相对破碎与破碎指数间的数学关系。定性分析了系列大型三轴试验过程中粗粒料级配变化规律,给出了各类型破碎指标的数学表达式,揭示了单调加载条件下粗粒料颗粒破碎完整规律：剪切过程中全局破碎关于归一化偏应力与平均主应力呈负指数增长关系,而临界状态和固结完成时的全局破碎只关于归一化平均主应力呈负指数增长关系,局部破碎只关于应力水平呈双曲线关系且在临界状态时均归为1。
- 粗粒料 /
- 级配 /
- 颗粒破碎 /
- 破碎指标 /
- 三轴试验
Abstract: By analyzing the progressive evolution of particle crushing of crushable soils, four basic properties of particle crushing are presented, and three characteristic states are also given over long-term or repeated loading. Two breakage indexes, namely the global and or breakage index, are defined for the crushable soils to distinguish temporary stop from perpetual termination of particle crushing, and a mathematical relationship is put forward between the revised relative breakage and the grading index. A series of consolidated drained large-scale triaxial tests are conducted for the coarse-grained materials (CGMs). The evolution of particle grading and breakage is analyzed, and some mathematical formulations are presented for various breakage indexes, which exhibits a entire rule of particle crushing of CGMs under the monotonic loads: (1) The global breakage has a negative exponential increase with the non-moralized deviatoric stress and the mean principal stress upon shearing. (2) At the end of consolidation or in the critical state, it increases merely with the non-moralized mean principal stress which meets the above-mentioned negative exponential function. (3) The local breakage varies with the stress level in a way of hyperbolic function and approaches 1 in the critical state.
- coarse-grained material /
- grading /
- particle crushing /
- breakage index /
- triaxial test

HTML全文

0. 引言

基础频域阻抗是土结动力相互作用研究中描述基础的一种重要方法。目前,关于基础阻抗函数解析解已有大量的研究成果,如：崔春义等^[1]利用桩–土耦合条件推导得到管桩基础的动力阻抗解析解；杨林青等^[2]求解了任意形状基础的频域动力阻抗。频域阻抗只适用于结构线性响应的频域求解,无法用于非线性响应分析和混合试验要求的时域逐步积分求解。为了将该频域阻抗用于土–结相互作用系统动力时程分析,需要将其转换为时域阻抗模型。目前常用的拟合方法有离散有理近似^[3]和连续有理近似函数^[4]两种,有理近似函数的稳定性与精度决定了动力时域分析的稳定性与精度。采用离散时间有理近似函数对基础频域阻抗进行准确拟合主要取决于参数识别方法的有效性。Paronesso等^[5]采用平衡逼近的方法推导了相互作用力与位移的差分方程,并将其表示为离散时间的有理近似函数；Wolf等^[6]、Safak等^[7]、杜修力等^[8]分别利用有理近似函数拟合了基础阻抗,并给出了后验的函数稳定性条件为函数极点的模小于1；Laudon等^[9]采用了Safak等^[7]的方法拟合了频域阻抗,指出该方法虽然可以保证拟合精度,但存在失稳问题；Gash等^[10]将离散时间函数与Newmark-β法结合,评价了滤波器–积分器所需的稳定条件；Sung等^[11]采用最小二乘法对阻抗进行了拟合,分析了土结动力相互作用对结构动力响应的影响。以上研究表明,大多频域阻抗的最优有理近似存在失稳问题,需要一个有效的方法保证识别函数的稳定性。赵建锋等^[12]利用罚函数给失稳的参数组合增加惩罚值,保证了识别程序输出结果的稳定性；赵密等^[3]先在S域中利用连续时间有理近似函数与罚函数法拟合基础阻抗,再通过双线性变换得到Z域的离散时间函数,保证了函数的稳定性。

离散时间有理近似函数的参数识别方法主要存在如下问题：部分方法未对函数的稳定性进行判断,可能导致时域分析失稳；另一部分方法采用后验方式判断辨识函数的稳定性,虽然可以较好地确保稳定性和精度,但由于辨识参数在整个实数域以及参数稳定性判别的复杂性,参数辨识效率较低,在实际中（特别是多自由度基础系统）难以应用。本文从控制理论角度给出了完整的稳定约束条件,并据此提出了一种高效稳定的离散时间有理近似函数参数识别方法。

1. 离散时间有理近似的稳定性

采用解析方法或数值方法求解基础频域动力阻抗^[4]是土结动力相互作用研究中基础的重要描述方法。根据输入与输出物理量不同,频域阻抗函数可以分为两种：广义位移作为输入,广义力作为输出的动力刚度阻抗函数S(ω)；广义力作为输入,广义位移作为输出的动力柔度阻抗函数F(ω)。因此,在频域中基础与无限土体之间力与位移的关系可以表示为

$f (ω) = S (ω) u (ω)$ ,

(1)

式中,f(ω)为广义力,u(ω)为广义位移,S(ω)为基础频域阻抗函数,其可以表示为频率的函数：

$S (ω) = S_{0} [K (ω) + i C (ω)]$ 。

(2)

式中 ω为外荷载频率；S₀为静刚度；i= $\sqrt{(- 1)}$ 为虚数单位；K(ω),C(ω)分别为与频率相关的刚度系数和阻尼系数。采用离散时间有理近似作为时域模型时,其对应的基础频域阻抗S(ω)可描述为

$S (ω) \approx S (z) = S_{0} \frac{b_{0} + b_{1} z^{- 1} + \dots + b_{M} z^{- M}}{1 + a_{1} z^{- 1} + \dots + a_{N} z^{- N}}$ 。

(3)

式中 b_j和a_j分别为分母和分子多项式的待定系数,且均为实数；z=e^iωΔt为Z变换的频域符号,Δt为离散时间步长；M,N分别为分子与分母多项式的阶数,通常根据基础频域阻抗复杂程度选择合适的阶次即可,频域阻抗越复杂需要的阶次越高。该模型的参数识别过程中只需用到离散的频域阻抗值,在采用式（3）进行拟合时,无需K(ω)与C(ω)的具体表达形式,而得到式（3）所示模型后,可以用不同时域模型对其进行等效。为了保证时域模型的稳定性,需要首先保证式（3）所示频域模型的稳定性。根据系统稳定性理论^[13],离散时间有理近似函数的稳定性可由其极点确定：当且仅当其极点全部位于单位圆内,即所有极点的模均小于1时,模型稳定。假设形如式（3）的离散时间有理近似函数的极点为z_j,则其稳定条件可表示为

$| z_{j} | < 1$ 。

(4)

以一个3×3群桩模型为例,假设地基土分为3层,模型如图1所示。其中桩直径d=1 m,桩长L=6 m,土层厚度h=2 m,土泊松比 $ν$ =0.35,土阻尼比β=0.05,桩距S=5 m,桩土密度比 $ρ_{p} / ρ_{s}$ =1.3,混凝土弹性模量E_p=30000 MPa,3层土的桩土弹模比E_p/E_s1=1000,E_p/E_s2=2000、E_p/E_s3=500。土体阻抗由文献[14]的方法计算得到,群桩阻抗由文献[15]的方法计算得到。采用时域递归算法^[12]对其进行描述,假设基础–土界面之间的相互作用力为f,则其可采用前M个时刻的界面位移和前N个时刻的相互作用力描述,其时域力学模型如图2,力学方程如下所示：

图 1 3×3群桩基础模型

Figure 1. Model for foundation with 3×3 pile groups

下载: 全尺寸图片幻灯片

图 2 时域递归模型计算流程

Figure 2. Computational process of time-domain recursive model

下载: 全尺寸图片幻灯片

$f (t) = \sum_{i = 0}^{M} b_{i} u (t - i Δ t) - \sum_{i = 1}^{N} a_{i} f (t - i Δ t)$ 。

(5)

并采用文献[12]给出的修正有理近似函数辨识群桩阻抗进行无约束识别,计算结果如图3所示。由图可知,虽然在频域内有理函数可以精确地与基础阻抗吻合,但由于拟合过程中未考虑函数稳定性问题,因此导致了时程分析发散。因此,对于离散时间有理近似函数,仅令其函数值在频域内精确地拟合基础阻抗函数是不够的,还需要确保时域分析的收敛性才具有应用价值。

图 3 3×3群桩动力响应

Figure 3. Dynamic responses of foundation with 3×3 pile groups

下载: 全尺寸图片幻灯片

2. 稳定参数识别方法

为了确保基础时域模型的收敛性,现有识别方法先满足频域响应精度,而后进行稳定性判别,当发现不稳定时重新识别参数,直到获得同时满足精度和稳定性条件的参数后停止识别。该方法的缺陷在于无法识别过程和结果不可控。如果预先确定满足稳定性的参数条件,然后在该范围内去寻找满足精度要求的参数将会更有效。

2.1 稳定参数约束条件

由式（4）可知,离散时间有理近似函数稳定的充分必要条件为所有极点的模均小于1,但由式（3）可知,只有该模型的所有参数已知才能得到其根,这就是现有方法无法提前确定稳定参数范围的原因所在。从控制系统角度来看,式（3）所示的多阶系统由多个一阶和二阶系统组成。因此,形如式（3）的离散时间有理近似函数,可用下式等效：

$S (ω) = \frac{\sum_{j = 0}^{m} b_{j} z^{- j}}{\prod_{j = 1}^{l} (1 + x_{1 j} z^{- 1} + x_{2 j} z^{- 2}) \times \prod_{j = 1}^{n - 2 l} (1 + x_{0 j} z^{- 1})}$ ,

(6)

式中,x_1j和x_2j为二阶系统系数,x_0j为一阶系统系数,l为二阶系统总数。

当式（3）等效为式（6）后,辨识多阶系统参数变为了辨识一阶和二阶系统参数。由一阶、二阶系统的根可得到其解析式,也就可以根据式（4）限制（6）中参数的取值范围了。式（6）每个一阶与二阶项的根为

$z_{j} = {\begin{matrix} \frac{- x_{1 j} \pm \sqrt{x_{1 j}^{2} - 4 x_{2 j}}}{2} \\ - x_{0 j} \end{matrix}$ 。

(7)

式中,一阶项根x_0j的取值范围可以简单地确定：

$- 1 < x_{0 j} < 1$ 。

(8)

而二阶项的根需分为实数根与复数根。由于未得到系数x的值时,不能判断一阶和二阶系统个数。因此,当n为偶数时,假设一阶项的个数为0,即n=2l；当n为奇数时,假设一阶项的个数为1,即n-1=2l。而当求出未知系数x后,若二阶项的根为实数根,则可自动变为两个一阶项。因此,二阶系统的参数范围分为实极点和复极点分别讨论。

（1）复极点

由于根为复数根,因此x_1j²-4x_2j<0,式（7）的二阶项根改写为

$z = \frac{- x_{1 j} \pm (\sqrt{4 x_{2 j} - x_{1 j}^{2}}) i}{2}$ 。

(9)

而函数稳定要求极点的模小于1,因此存在：

$\begin{array}{l} \frac{x_{1 j}^{2}}{4} + \frac{4 x_{2 j} - x_{1 j}^{2}}{4} < 1, \\ x_{1 j}^{2} - 4 x_{2 j} < 0 。 \end{array}}$

(10)

进一步化简得到x的取值范围：

$\begin{array}{l} 0 \leq x_{2 j} < 1, \\ - 2 \sqrt{x_{2 j}} \leq x_{1 j} \leq 2 \sqrt{x_{2 j}} 。 \end{array}}$

(11)

由式（11）可知,x_1j的取值范围受x_2j的影响,需要非线性约束,其取值范围如图4所示阴影范围。

图 4 复数根时稳定参数界限

Figure 4. Stability parameter boundary for complex roots

下载: 全尺寸图片幻灯片

（2）实极点

当极点为实数时,极点的模小于1的条件为

$\begin{array}{l} | \frac{- x_{1 j} \pm \sqrt{x_{1 j}^{2} - 4 x_{2 j}}}{2} | < 1, \\ x_{1 j}^{2} - 4 x_{2 j} \geq 0 。 \end{array}}$

(12)

将式（12）进一步化简可得

$\begin{array}{l} x_{1 j} - 2 < \sqrt{x_{1 j}^{2} - 4 x_{2 j}} < 2 + x_{1 j}, \\ - x_{1 j} - 2 < \sqrt{x_{1 j}^{2} - 4 x_{2 j}} < 2 - x_{1 j} 。 \end{array}}$

(13)

而由于根为实数根,在实数范围内实数的偶次平方根必然为非负数,因此式（13）不等式组的右端需大于0,可以得到x_1j的取值范围：

$- 2 < x_{1 j} < 2$ 。

(14)

进而根据式（14）可以判断出式（13）的左端必然成立,将不等式组的左端去掉并对两侧取平方可得

$\begin{array}{l} x_{1 j} - 1 < x_{2 j}, \\ - x_{1 j} - 1 < x_{2 j}, \\ x_{1 j}^{2} - 4 x_{2 j} \geq 0, \\ - 2 < x_{2 j} < 2 。 \end{array}}$

(15)

式（15）的不等式组需要取交集以保证式（13）必然成立,则取值范围可简化为式（16）,为便于理解,取值范围可见图5,图中深色填充部分为式（16）所示取值范围,两条交叉直线分别为x_2j=x_1j-1和x_2j=-x_1j-1,抛物线为4x_2j=x_1j²。

图 5 实数根时稳定参数界限

Figure 5. Stability parameter boundary for real roots

下载: 全尺寸图片幻灯片

$\begin{matrix} x_{1 j} - 1 < x_{2 j} \leq \frac{x_{1 j}^{2}}{4} & 0 \leq x_{1 j} < 2, \\ - x_{1 j} - 1 < x_{2 j} < \frac{x_{1 j}^{2}}{4} & - 2 < x_{1 j} < 0 。 \end{matrix}}$

(16)

（3）稳定参数界限

观察式（11）,（16）和图4,5可以发现,当共轭根分别为复数与实数时,参数x的取值范围以曲线为4x_2j=x_1j²界线分为上下两部分。因此,可将两部分合并,则式（6）所示模型的稳定参数界限如式（17）,其参数范围如图6所示,图中浅色填充部分为复数根的取值范围,深色填充部分为实数根的取值范围。

图 6 稳定参数界限

Figure 6. Boundary of stable parameters

下载: 全尺寸图片幻灯片

$\begin{array}{l} - 2 < x_{1 j} < 2, \\ - 1 < x_{2 j} < 1, \\ x_{1 j} - 1 < x_{2 j} (0 \leq x_{1 j} < 2), \\ - x_{1 j} - 1 < x_{2 j} (- 2 < x_{1 j} < 0) 。 \end{array}}$

(17)

而当拟合动力柔度频响函数时,则需要保证函数零点的模小于1,此时可将分子提取b₀,然后再采用部分分式的方法得到相同的取值界限。

2.2 参数识别

得到了新的有理近似函数形式及其稳定参数取值边界后,如何通过参数识别确保最优的精度是参数寻优问题。优化算法可以分为全局优化算法和局部优化算法,常见的全局优化算法有进化算法、遗传算法、粒子群优化算法等等,而常见的局部算法有单纯形法、序列二次规划算法等等。

本文以Wang等^[16]将遗传算法与单纯形法结合形成的混合算法为基础,采用序列二次规划算法(SQP算法)代替单纯形法,形成遗传–序列二次规划算法,用于求解有理近似函数的待定系数。遗传算法作为全局优化算法,不需要给出合理的初始值,其从多个任意初始值开始寻优,搜索范围大,利于全局择优。而SQP算法作为局部优化算法可以从遗传算法得到的全局最优解中选择初始值,避免了因初值选择不适当而造成的拟合效果不佳。此外,采用SQP算法代替单纯形法可以将无约束优化问题变为有约束优化问题,通过约束方程人为地规定参数的取值范围,更适用于本文建立的稳定参数边界。具体的计算流程如图7所示,预先设定有理近似函数阶次和预期误差（比如均方根小于10%）,当辨识参数对应阻抗函数精度满足预期误差条件时跳出循环；当辨识参数对应精度无法满足预期误差条件时,设定连续10次循环每两次输出结果精度之差均小于0.1%时跳出循环,增加有理近似函数阶次,重复前面的循环,直至辨识精度满足预期误差条件。

图 7 参数识别流程

Figure 7. Flow chart of parameter identification

下载: 全尺寸图片幻灯片

3. 数值算例验证

为了验证上述方法的有效性,本节分别采用3组复杂程度不同的基础阻抗进行对比研究。

（1）算例1：3×3群桩基础

对图2所示的3×3群桩基础的频域阻抗,本文采用文献[12]提出的修正离散时间有理近似函数进行辨识,辨识所用目标函数如下式所示：

$f = \sum_{Ω} ({| \frac{R (Ω) - S (Ω)}{S_{0}} |}^{2})$ ,

(18)

式中, $Ω$ 为外荷载频率ω或其无量纲频率a₀,视函数形式而定,但计算过程中需保持量纲一致,R( $Ω$ )为离散时间有理近似函数,S( $Ω$ )为频域阻抗函数。

为了综合验证本文方法在稳定性、精度和计算效率,对图3所示无约束识别失稳案例采用罚函数法^[12]和本文方法进行对比分析,分别对有理函数m=n=3/4/5/6/7阶进行了分析。为了对辨识模型精度进行定量评价,取频域函数精确解与识别函数之差绝对值积分面积占精确阻抗绝对值积分面积的比值作为误差评价指标,其计算方法如下式所示,其中虚部和实部误差权重均定义为0.5。

$\begin{matrix} E = (\frac{\int_{0}^{ω_{N}} | R {(ω)}_{Re} - S {(ω)}_{Re} | d ω}{\int_{0}^{ω_{N}} | S {(ω)}_{Re} | d ω} \times 0.5 + \\ \frac{\int_{0}^{ω_{N}} | R {(ω)}_{Im} - S {(ω)}_{Im} | d ω}{\int_{0}^{ω_{N}} | S {(ω)}_{Im} | d ω} \times 0.5) \times 100 % \end{matrix}$ 。

(19)

式中 E表示识别模型误差,由实部与虚部的权重各一半；R(ω)_Re,S(ω)_Re,R(ω)_Im和S(ω)_Im分别表示有理函数与阻抗精确值的实部与虚部；ω_N表示频率最大值。辨识结果如表1和图8所示。

表 1 3×3群桩基础辨识效率

Table 1. Identification efficiency of foundation with 3×3 pile groups

阶次	3	4	5	6	7
本文方法/s	1.82	2.24	2.24	2.86	2.92
罚函数法/s	3.53	4.67	8.41	16.17	33.67

下载: 导出CSV

| 显示表格

图 8 6阶有理函数拟合结果

Figure 8. Identification results of six-order rational function

下载: 全尺寸图片幻灯片

由图8（a）,8（b）可以看到,两种方法辨识所得模型均可较精确地描述频域阻抗,采用式（19）计算所得误差小于1%,证明两种方法拥有很高精度。由图8（c）,8（d）可以看到,将EL-centro波作为地震动输入,两种方法均可保证时程计算的稳定性,且时程曲线一致。由表1可知,在计算效率方面上,本文方法明显优于罚函数法。两种方法分别采用3到7阶函数对基础阻抗进行识别,本文方法所需时间均在3 s以内,且效率损失随函数阶数增加并不明显,而罚函数法所需时间在3～34 s,为本文方法的2倍到10倍,且拟合用时随阶数增加几乎成倍增长。原因在于：第一,本文方法给定了参数取值区间,而罚函数法的参数取值范围为整个实数域,取值范围的增大造成了计算效率的降低；第二,本文方法取值区间内可保证函数稳定,但罚函数法需要求解高阶多项式判断稳定性,进一步降低了计算效率。

（2）算例2：5×5群桩基础

为验证阶数对于辨识精度的影响,分别采用不同阶数的有理近似函数对图9和图10所示5×5群桩基础的阻抗进行拟合,其阻抗由文献[17]的图6所得。图中d为桩基础直径,S为桩间距,桩中心距筏板边缘长0.5S。由于此基础动力阻抗对频率变化剧烈,因此利用离散时间的有理近似函数拟合具有挑战性。本文采用式（3）给出的有理函数,取Δt=0.02 s,阶数m=n=3～7阶进行拟合。函数拟合精度采用式（19）进行判断,结果如图10和表2。

图 9 5×5群桩基础模型

Figure 9. Model for foundation with 5×5 pile groups

下载: 全尺寸图片幻灯片

图 10 5×5群桩基础不同阶次拟合精度

Figure 10. Identification accuracy of foundation with 5×5 pile groups in different orders

下载: 全尺寸图片幻灯片

表 2 5×5群桩基础拟合精度与效率

Table 2. Identification accuracy and efficiency of foundation with 5×5 pile groups

阶次	3	4	5	6	7
误差/%	18.80	18.23	13.02	2.03	1.97
时间/s	9.02	8.15	9.81	13.74	16.46

下载: 导出CSV

| 显示表格

由表2和图10可以看到,采用不同阶次的有理函数辨识时,辨识精度随阶数增加逐渐提高,当采用3阶函数拟合时,效果较差,整体误差约有19%,难以拟合复杂阻抗。但当阶数增加到6阶时,整体误差降低到2%,此时有理函数可较精确地吻合频域阻抗。从辨识效率来看,3阶函数时拟合用时约为9 s,7阶时约为16 s,增加了仅一倍,可以看出拟合效率的损失随阶次增加较低。采用同阶次的罚函数法进行拟合由于计算效率过低未得到有效结果,证明本文方法更适用于采用高阶次函数对复杂阻抗函数进行准确辨识。

（3）算例3：水平运动的圆形基础

对于随频率变化更为复杂的阻抗函数,以文献[7]中图7的圆形基础水平运动动力刚度函数为例。基础模型如图11所示,图中Pe^iωΔt和ue^iωΔt分别为基础受到的外荷载与产生的位移。基础参数为：半径r=10 m,剪切波速V_s=100 m/s,土体泊松比 $ν$ =1/3,土体阻尼比 $ζ$ =0.05,土层厚度H=2 m。拟合阶数取m=n=3/5/7/9/11阶,拟合精度采用式（19）进行判断,结果如图12与表3所示。

图 11 水平运动的圆形基础

Figure 11. Circular foundation under horizontal motion

下载: 全尺寸图片幻灯片

图 12 圆形基础不同阶次拟合精度

Figure 12. Identification accuracy of circular foundation in different orders

下载: 全尺寸图片幻灯片

表 3 圆形基础拟合精度与效率

Table 3. Identification accuracy and efficiency of circular foundation

阶次	3	5	7	9	11
误差/%	16.82	12.29	7.49	6.74	6.62
时间/s	4.46	4.19	5.79	7.12	8.66

下载: 导出CSV

| 显示表格

对于此圆形基础阻抗,罚函数法仍无法辨识出有效结果,而本文方法结果如表3与图12所示。由表3可知,随有理函数阶次增加,拟合用时逐渐增加,而拟合精度也逐渐提高,但总体用时低于10 s,误差可降低到7%,通过图12的不同阶次对比可知：增加阶次有利于提高阻函数拟合精度,在达到7阶后,增加阶次对精度提高不明显,此时已可以较好的满足应用需求。

4. 结论

针对基础频响离散时间有理近似函数的参数识别难以同时保证稳定性、精度及计算效率的问题,本文提出了时域绝对稳定的离散有理近似函数参数识别方法。采用不同基础阻抗函数,对本文方法和罚函数拟合方法就稳定性、精度及计算效率等方面进行了对比分析,验证了本文方法的有效性与优势。

（1）从线性系统理论角度将基础阻抗离散有理近似函数看成一阶与二阶系统的组合,并根据其根的稳定条件建立了被辨识参数的稳定界限。据此,采用遗传–序列二次规划算法建立了时域稳定的参数识别方法,保证了识别函数时域模型的绝对稳定性。

（2）对单自由度基础阻抗对比仿真表明,本文方法与既有方法拥有同等精度,但计算效率明显优于既有算法。对于简单阻抗函数,本文方法拟合用时均低于3 s,为同等精度既有方法的10%～50%；对于不同的复杂函数,本文方法用时也均低于20 s,提高了多自由度阻抗函数识别的适用性。

（3）通过对复杂基础阻抗的辨识仿真表明,增加有理函数阶次有利于提高拟合精度,但过高的阶次对提高精度不明显。根据本文仿真结果,7阶有理近似函数能比较好的描述复杂基础阻抗。

（4）本文基于单输入单输出控制理论建立了单自由度基础阻抗的时域模型参数识别方法,对于多自由度体系,动力分析方法与有理近似函数特性相互耦合,其稳定条件更为复杂,将在后续工作中开展进一步研究。

图 1 粗粒料级配曲线

Figure 1. Grading curve of coarse-grained materials

下载: 全尺寸图片幻灯片

图 2 粗粒料应力应变与体积应变测试结果

Figure 2. Test results of stress-strain-volume change of coarse-grained materials

下载: 全尺寸图片幻灯片

图 3 大型三轴试验各个阶段粗粒料的级配曲线

Figure 3. Grading curves of coarse-grained materials in various phases of large-scale triaxial tests

下载: 全尺寸图片幻灯片

图 4 大型三轴试验各个阶段粗粒料级配曲线形状指数变化规律

Figure 4. Evolution of shape index of grading curves of coarse-grained materials in various phases of large-scale triaxial tests

下载: 全尺寸图片幻灯片

图 5 大型三轴试验各个阶段粗粒料级配参数变化规律

Figure 5. Evolution of grading parameters of coarse-grained materials in various phases of large-scale triaxial tests

下载: 全尺寸图片幻灯片

图 6 大型三轴试验各阶段粗粒料B_r-p/p_a曲线

Figure 6. B_r-p/p_a curves of coarse-grained materials in various phases of large-scale triaxial tests

下载: 全尺寸图片幻灯片

图 7 大型三轴试验各阶段粗粒料B_r-q/p_a曲线

Figure 7. B_r-q/p_a curves of coarse-grained materials in various phases of large-scale triaxial tests

下载: 全尺寸图片幻灯片

图 8 固结完成与临界状态粗粒料B_r0/B_rc-p/p_a拟合曲线

Figure 8. Fitting curves of B_r0/B_rc-p/p_a of coarse-grained materials at end of consolidation and in critical state

下载: 全尺寸图片幻灯片

图 9 大型三轴试验各阶段粗粒料 ${B^{'}}_{r}$ -S_L拟合曲线

Figure 9. Fitting curves of ${B^{'}}_{r}$ -S_L of coarse-grained materials in..various phases of large-scale triaxial tests

下载: 全尺寸图片幻灯片

表 1 粗粒料基本物理力学性质指标

Table 1 Basic physical and mechanical indexes of coarse-grained materials

土粒相对密度G_s	相对密实度D_r/%	制样干密度/(g·cm^-3)	制样孔隙比e₀	母岩抗压强度(饱和)/MPa	压缩模量E_s4-8/MPa	压缩指数C_c	再压缩指数C_s	渗透系数k/(cm·s^-1)
2.72	96	2.10	0.295	86.9	93.5	1.13×10^-1	1.48×10^-3	0.227

下载: 导出CSV

参考文献(20)

[1]	日本土质工学会(日). 堆石料的现场压实[M]. 郭熙灵,等译.北京: 中国水利水电出版社, 1999. Association of Geotechnique of Japan. In-Situ Compaction of Coarse- Grained Materials[M]. GUO Xi-ling, et al trans. Beijing: China Water Power Press, 1999. (in Chinese)
[2]	SEBASTINA L G. Evaluation of Crushing in Granular Materials Using the Discrete Element Method and Fractal Theory[D]. Pittsburgh: University of Pittsburgh, 2006.
[3]	MARSAL R J. Large scale testing of rockfill materials[J]. Journal of the Soil Mechanics and Foundations Division, ASCE, 1967, 93(2): 27-43. doi: 10.1061/JSFEAQ.0000958
[4]	HARDIN B O. Crushing of soil particles[J]. Journal of Geotechnical Engineering, ASCE, 1985, 111(10): 1177-1192. doi: 10.1061/(ASCE)0733-9410(1985)111:10(1177)
[5]	LADE P V, YAMAMURO J A, BOPP P A. Significance of particle crushing in granular materials[J]. Journal of Geotechnical Engineering, ASCE, 1996, 122(4): 309-316. doi: 10.1061/(ASCE)0733-9410(1996)122:4(309)
[6]	MUIR W D, MAEDA K. Changing grading of soil: effect on critical state[J]. Acta Geotechnica, 2008, 3(1): 3-14. doi: 10.1007/s11440-007-0041-0
[7]	XIAO Y, LIU H L. Elastoplastic constitutive model for rockfill materials considering particle breakage[J]. International Journal of Geomechanics, ASCE, 2017, 17(1): 1-13.
[8]	EINAV I. Breakage mechanics-part I: theory[J]. Journal of the Mechanics and Physics of Solids, 2007, 55(6): 1274-1297. doi: 10.1016/j.jmps.2006.11.003
[9]	MCDOWELLL G R, BOLTON M D. On the micromechanics of crushable aggregates[J]. Géotechnique, 1998, 48(5): 667-679. doi: 10.1680/geot.1998.48.5.667
[10]	尹振宇, 许强, 胡伟. 考虑颗粒破碎效应的粒状材料本构研究：进展及发展[J]. 岩土工程学报, 2012, 34(12): 2170-2180. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201212006.htm YIN Zhen-yu, XU Qiang, HU Wei. Constitutive relations for granular materials considering particle crushing: review and development[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(12): 2170-2180. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201212006.htm
[11]	KIKUMOTO M, MUIR WOOD D, RUSSELL A. Particle crushing and deformation behavior[J]. Soils and Foundations, 2010, 50(4): 547-563. doi: 10.3208/sandf.50.547
[12]	KAN E M, TAIEBAT H A. A bounding surface plasticity model for highly crushable[J]. Soils and Foundations, 2014, 54(6): 1188-1201. doi: 10.1016/j.sandf.2014.11.012
[13]	LIU H B, ZOU D G. Associated generalized plasticity framework for modeling gravelly soils considering particle breakage[J]. Journal of Engineering Mechanics, ASCE, 2013, 139(5): 606-615. doi: 10.1061/(ASCE)EM.1943-7889.0000513
[14]	LIU M C, ZHANG Y, ZHU H Z. 3D elastoplastic model for crushable soils with explicit formulation of particle crushing[J]. Journal of Engineering Mechanics, ASCE, 2017, 143(12): 1-18.
[15]	LIU M C, GAO Y F, LIU H L. An elastoplastic constitutive model for rockfills incorporating energy dissipation of nonlinear friction and particle breakage[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2014, 38(9): 935-960. doi: 10.1002/nag.2243
[16]	刘萌成, 高玉峰, 刘汉龙. 模拟堆石料颗粒破碎对强度变形的影响[J]. 岩土工程学报, 2011, 33(11): 1691-1699. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201111010.htm LIU Meng-cheng, GAO Yu-feng, LIU Han-long. Effect of particle breakage on strength and deformation of modeled rockfills[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(11): 1691-1699. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201111010.htm
[17]	TURCOTTE D L. Fractals and fragmentation[J]. Journal of Geophysical Research, 1986, 91(B2): 1921-1926. doi: 10.1029/JB091iB02p01921
[18]	SAMMIS C G, KING G, BIEGEL R. The kinematics of gouge deformations[J]. Pure Applied Geophysics, 1987, 125(5): 777-812.
[19]	刘萌成, 高玉峰, 刘汉龙. 应力路径条件下堆石料剪切特性大型三轴试验研究[J]. 岩石力学与工程学报, 2008, 27(1): 176-186. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200801025.htm LIU Meng-cheng, GAO Yu-feng, LIU Han-long. Study on shear behaviors of rockfill in large-scale triaxial tests under different stress paths[J]. Chinese Journal of Rock Mechanics and Engineering, 2008, 27(1): 176-186. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200801025.htm
[20]	JIA Y F, XU B, CHI S C, et al. Research on the particle breakage of rockfill materials during triaxial tests[J]. International Journal of Geomechanics, ASCE, 2017, 17(10): 1-11.

施引文献(9)

期刊类型引用(3)

1.	徐荣超，王怡博，张圣哲，孟凡震，李震，高梁，阎震. 轴向及环向变形加载控制方式对花岗岩力学及声发射特性影响的试验研究. 工程地质学报. 2025(01): 123-134 . 百度学术
2.	胡玉波，房敬年，徐荣超，郝小红，阎震，周文朋，李震. 饱水作用对白砂岩损伤破裂及声发射特性的影响. 长江科学院院报. 2024(11): 163-171 . 百度学术
3.	李泓颖，刘晓辉，郑钰，肖文根. 深埋锦屏大理岩渐进破坏过程中的特征能量分析. 岩石力学与工程学报. 2022(S2): 3229-3239 . 百度学术