• ISSN 1673-5722
  • CN 11-5429/P

应用非线性高阶谐波衰减和尾波干涉监测高温作用后水泥制品的损伤演化

张尧 马强 肖武军

任争争, 梅雨辰, 李鸿晶. 微分求积法在结构动力分析中的应用[J]. 震灾防御技术, 2018, 13(4): 829-838. doi: 10.11899/zzfy20180410
引用本文: 张尧, 马强, 肖武军. 应用非线性高阶谐波衰减和尾波干涉监测高温作用后水泥制品的损伤演化[J]. 震灾防御技术, 2018, 13(1): 52-64. doi: 10.11899/zzfy20180105
Ren Zhengzheng, Mei Yuchen, Li Hongjing. The Application of Differential Quadrature Method in Structural Dynamic Analysis[J]. Technology for Earthquake Disaster Prevention, 2018, 13(4): 829-838. doi: 10.11899/zzfy20180410
Citation: Zhang Yao, Ma Qiang, Xiao Wujun. Application Nonlinear High Order Harmonics and Coda Wave Interferometry on Monitoring Damage Evolution of Cement Specimens Subject to Elevated Temperature[J]. Technology for Earthquake Disaster Prevention, 2018, 13(1): 52-64. doi: 10.11899/zzfy20180105

应用非线性高阶谐波衰减和尾波干涉监测高温作用后水泥制品的损伤演化

doi: 10.11899/zzfy20180105
详细信息
    作者简介:

    张尧, 男, 生于1986年。工程师。主要从事重大项目管理, 抗震性能模拟分析等研究。E-mail:zhangyao@seis.ac.cn

Application Nonlinear High Order Harmonics and Coda Wave Interferometry on Monitoring Damage Evolution of Cement Specimens Subject to Elevated Temperature

  • 摘要: 非线性高阶谐波和尾波波速变化均能够反映水泥材料内部微结构的应力变化。利用高阶谐波和尾波干涉实验测量系统,对引入高温作用后的3类不同粒径共6块水泥试样进行单轴加载的损伤演化实验,并与无高温作用的完整试样的实验结果进行对比。结果表明,从初始状态到25%抗压强度的过程中,高温作用后的试样的谐波幅值和尾波波速变化出现明显增强的现象(谐波幅值最大增幅约20%),而无高温作用的完整试样的谐波幅值和尾波波速变化较为平稳(谐波幅值最大增幅约5%);在达到65%抗压强度的过程中,高温作用后的试样的谐波幅值和尾波波速变化急剧增大(谐波幅值最大增幅约100%),且粒径较大的试样的增幅高于粒径较小的试样,而无高温作用的完整试样的谐波幅值和尾波波速变化的增幅较小(谐波幅值最大增幅约10%);当抗压强度超过75%以后,高温作用后的试样的谐波幅值和尾波波速变化急剧衰减(谐波幅值最大衰减幅度约140%),而无高温作用的完整试样的谐波幅值和尾波波速的最大衰减幅值在40%以内。基于以上观测结果对高温作用后水泥制品损伤演化的物理机制以及这两类监测方法的适用性进行了讨论。
  • 地震等灾害环境的作用使结构产生复杂的动态响应,可能导致结构失效、破坏甚至倒塌,从而形成灾害、造成损失。对新建结构进行抗灾设计及对已建结构进行抗灾加固是防御和减轻工程灾害及其损失的有效途径,这就要求认识结构在这些灾害作用下的性态和响应行为,而结构动力分析则是实现这一目标的基本手段。

    结构动力分析的本质是实现对动力荷载激励下的结构运动微分方程(组)的求解,属于常微分方程的初值问题。目前用于分析结构动力的方法大致可以分为2类:一类是变换方法,如振型叠加法利用振型的正交性和完备性将结构动态响应向各阶振型分解,再通过叠加各阶振型响应以获得结构动态响应的结果,这类方法采用叠加原理,一般只用于线弹性结构动力分析,且要求结构具有经典阻尼特性;另一类为直接方法,即直接对结构运动微分方程进行求解而不必引入任何假定,这类方法既可用于线弹性结构,也可用于非线性结构的分析,以Newmark-β法( Newmark,1959)等逐步积分法为代表,一般采用数值求解手段。

    由于现代结构不断向大型化、复杂化发展,加之结构精细化模型的采用,导致结构动力分析的计算需求呈爆发式增长,因而需要寻求高效率的分析方法。微分求积法(Differential Quadrature Method,DQM)是由Bellman等(19711972)发展起来的1种求解微分方程的数值方法,可通过较小的计算工作量获得较高的计算精度。DQM要求求解域比较规则,在工程分析中一般用于时间无关问题的求解。如在结构力学领域,多用来求解静力问题和固有振动问题等(Bert等,1988, 1993, 1996, 1997Jang等,1989Kang等,1995, 1996Liew等,1996a, 1996bSherbourne等,1991Striz等,1988Wang等,1993, 1994Zeng等,2001)。这类问题都属于边值问题,即在DQM中计算的是解函数对空间坐标的导数。而结构动力分析属于初值问题,使用DQM对其求解的研究工作相对较少,Fung(2001a2001b)Liu等(2008)李鸿晶等(2011a, 2011b)和廖旭等(2013)开展过相关的研究工作。本文在此基础上发展1种基于DQM的结构动力分析的高精度方法。不同于静力边值问题,实际结构的动力响应问题有其特殊性,许多情况下时间跨度很长,像边值问题一样一次性对所有时间区域进行离散求解,将出现病态问题而产生错误。本文借鉴单元法的思想,以期提高结构动力分析的计算效率。

    微分求积法是1种用于求解微分方程的数值方法。它的实质是将函数在某一离散节点处的各阶导数值,近似表示成计算域内所有节点处离散函数值的线性加权和,从而将复杂的微分方程化为关于离散点的线性方程(组)。由于本文讨论的是结构动力反应的计算,求解的运动微分方程仅是关于时间t的常微分方程,因此仅介绍一维区域内的微分求积原理。

    设函数f(x)为在区间[ab]上k阶连续可微,将区间[ab]划分为m段,共(m+1)个互不相同的节点,分别记为x0x1,……,xm-1xm,其中x0=axm=b

    根据计算数学的函数逼近理论,函数fx)可做如下逼近:

    $$ f(x) \approx \sum\limits_{j = 0}^m {{q_j}(x)f({x_j})} $$ (1)

    其中,$ {q_j}(x)$为函数空间中各线性无关的基函数。

    对式(1)求k阶导数,然后将所有节点代入,得到:

    $$ {f^{(k)}}({x_i}) \approx \sum\limits_{j = 0}^m {q_j^{(k)}({x_i})f({x_j}){\rm{ }}}, \;\;i = 0, \;1, \cdots \cdots, m $$ (2)

    记$a_{ij}^{(k)} = q_j^{(k)}({x_i}), {f^{(k)}}({x_i}) = f_i^{(k)}, f({x_j}) = {f_j} $,则:

    $$ f_i^{(k)} \approx \sum\limits_{j = 1}^N {{a^{(k)}}_{ij}{f_j}{\rm{ }}}, \;\;i = 0, \;1, \cdots \cdots, m $$ (3)

    式(3)即为一维区域微分求积的基本公式。

    微分求积法的关键问题是确定式(3)中的系数$ {a^{(k)}}_{ij}$。该系数称为微分求积权系数,它与网格节点的分布及选择的基函数空间类型有关,当节点的位置确定后,将基函数代入式(3),即可得到相应的权系数。

    常用的函数空间是(m+1)维多项式空间,选择该空间中的所有基函数都能得到相同的权系数。最常见的多项式基函数是幂指数插值函数和拉格朗日插值函数2种,分别如式(4)、(5)所示:

    $$ f(x) = {x^r}, r = 0, 1, \cdots \cdots, m $$ (4)
    $$ f(x) = \prod\limits_{i = 0\atop i \ne j}^m {\frac{{x - {x_i}}}{{{x_j} - {x_i}}}} \;\;{\rm{, }}\;\;j = 0, 1, \cdots \cdots, m $$ (5)

    选用式(4)的幂指数函数计算权系数,得到权系数的隐式表达式,需要求解范德蒙矩阵的逆矩阵,但当m较大时,不但计算量大,而且矩阵将容易出现病态。为了解决这些问题,目前大多采用式(5)给出的拉格朗日插值基函数,因为拉格朗日插值的各项系数为各节点的函数值,形式与式(1)完全相同,直接求k阶导数,便可得到相应的k阶权系数的显示表达式,计算效率大幅度提高。

    实际计算权系数时往往只需计算1阶导数的权系数,其它各阶导数的权系数可由1阶权系数以矩阵的形式方便地表示:

    $$ {{\bf{A}}^{k}}={({\bf{A}})^k} $$ (6)

    其中,A表示1阶权系数矩阵,Ak表示k阶权系数矩阵。$ {\bf{A}} = \left[ {\begin{array}{*{20}{c}} {{a_{00}}}&{{a_{01}}}& \cdots &{{a_{0m}}}\\ {{a_{10}}}&{{a_{11}}}& \cdots &{{a_{1m}}}\\ \vdots&\vdots&\ddots&\vdots \\ {{a_{m0}}}&{{a_{m1}}}& \cdots &{{a_{mm}}} \end{array}} \right]$,$ {{\bf{A}}^{\left(k \right)}} = \left[ {\begin{array}{*{20}{c}} {{a^{(k)}}_{00}}&{{a^{(k)}}_{01}}& \cdots &{{a^{(k)}}_{0m}}\\ {{a^{(k)}}_{10}}&{{a^{(k)}}_{11}}& \cdots &{{a^{(k)}}_{1m}}\\ \vdots&\vdots&\ddots&\vdots \\ {{a^{(k)}}_{m0}}&{{a^{(k)}}_{m1}}& \cdots &{{a^{(k)}}_{mm}} \end{array}} \right]$, 1阶权系数${a_{ij}} = {a^{\left( 1 \right)}}_{ij} $,可通过选定拉格朗日基函数,直接得到如下的表达式:

    $$ {a_{ij}} = \left\{ \begin{array}{l} \frac{{M({x_i})}}{{({x_i} - {x_j})M({x_j})}}\;\;\;, \;\;\;i \ne j\\ - \sum\limits_{k = 0\atop k \ne i}^m {{a_{ik}}} \;\;\;, \;\;\;i = j \end{array} \right. $$ (7)

    其中,$ M\left({{x_i}} \right){\rm{ = }}\prod\limits_{k = 0\atop k \ne i}^m {({x_i} - {x_k})} {\rm{ }}$,$M({x_j})= = \prod\limits_{k = 0\atop k \ne j}^m {({x_j} - {x_k})} $。

    应用微分求积法时还需确定采样网格节点的位置,大致分为均匀网格点和非均匀网格点2大类。虽然均匀网格点的精度总体上没有非均匀点高,但是其在处理离散荷载,如地震荷载或风荷载时,可直接使用原始采样点作为节点,不需要对荷载进行额外的插值,计算效率比不均匀网格点更高。

    用上述微分求积法求解结构的动力反应,并以线弹性单自由度体系为例进行讨论。虽然实际结构大都为多自由度体系,且为非线性体系,但其动力反应都可以通过线性迭代和振型分解转化为线弹性单自由度体系。因而,讨论线弹性单自由度体系的微分求积法更具普遍意义,且简单直观、易于理解。

    线弹性单自由度体系的动力反应运动方程为:

    $$ \frac{{{{\rm{d}}^2}u}}{{{\rm{d}}{t^2}}} + 2\xi \omega \frac{{{\rm{d}}u}}{{{\rm{d}}t}} + {\omega ^2}u = p(t) $$ (8)

    其中,ωξ分别表示体系的自振频率和阻尼比;u表示体系的位移;$ \frac{{{\rm{d}}u}}{{dt}}$、$ \frac{{{{\rm{d}}^2}u}}{{d{t^2}}}$分别表示体系的速度和加速度;pt)为结构所受到的随时间变化的荷载。

    在实际工程中,动荷载随类型的不同,作用时间差异很大,短则瞬间(如冲击荷载),长则几分钟甚至几十分钟(如风荷载)。对于一些长时间作用的荷载,在整个荷载作用时域内实施微分求积法,要保证计算结果的精确性,需要成千上万的时间节点,这将导致求解时系数矩阵的阶数过于庞大,引起矩阵的条件数大,病态效果严重,无法得到可靠的结果。为了解决此问题,可以借鉴有限单元法的思想,将整个荷载持时划分为多个等间距的时间单元,称为1个时步,在每个时步内,使用微分求积法求解。

    考虑动力荷载持时内长度$ \Delta t$的时步[tjtk]内的反应,tjtk一般与荷载pt)的采样时刻点重合。在时步内定义一局部坐标$tt \in {\rm{[0, }}\;\Delta t{\rm{]}} $,坐标起点为tj,方向同时间坐标。为了方便处理,正则化局部坐标,作$\tau= tt/\Delta t $,则定义域被正则化为[0, 1],时步内关于$\tau $的运动方程可写为:

    $$ \ddot u(\tau) + 2\xi \omega \Delta t\dot u(\tau) + {(\Delta t)^2}{\omega ^2}u(\tau) = {(\Delta t)^2}p(\tau \Delta t) $$ (9)

    其中,$\dot u(\tau) $和$\ddot u(\tau)$分别表示$ u(\tau)$对局部坐标τ的1阶和2阶导数。

    将时步离散为m段,记节点为τ1τ2,……,τm,将$ u(\tau= {\tau _i})$,$ \dot u(\tau = {\tau _i})$,$ \ddot u(\tau = {\tau _i})$,$p(\tau = {\tau _i}) $分别简记为ui,${\dot u_i} $,$ {\ddot u_i}$,pi,则在各离散节点处的方程为:

    $$ {\ddot u_i} + 2\xi \omega \Delta t{\dot u_i} + {(\Delta t)^2}{\omega ^2}u{}_i = {(\Delta t)^2}p{}_i, \;\;i = 0, \;1, \cdots \cdots, m $$ (10)

    对这些节点使用微分求积法:

    $$ {\dot u_i} = \sum\limits_{j = 0}^N {{a_{ij}}{u_j}}, \;\;i = 0, \;\;1, \cdots \cdots, \;m $$ (11)
    $$ {\ddot u_i} = \sum\limits_{k = 0}^m {{a_{ik}}{{\dot u}_k} = } \sum\limits_{k = 0}^m {{a_{ik}}\sum\limits_{j = 0}^m {{a_{kj}}{u_j}} }, \;\;i = 0, \;\;1, \cdots \cdots, \;m $$ (12)

    其中,${a_{ij}}$是微分求积的权系数,将式(11)、(12)写成矩阵的形式:

    $$ {\bf{\dot u}}={\bf{Au}} $$ (13)
    $$ {\bf{\ddot u}} = {{\bf{A}}^2}{\bf{u}} $$ (14)

    其中,$ {\bf{u}} = {({u_0}, {u_1}, \ldots \ldots, {u_m})^{\rm{T}}}$。

    已知时步的初始位移和速度分别为u0v0,且$ {\dot u_0} = {v_0}\Delta t$,将式(11)、(12)改写为:

    $$ {{\bf{\dot u}}_{\bf{s}}} = {{\bf{G}}_{\bf{0}}}u{}_0 + {\bf{G}}{{\bf{u}}_{\bf{s}}} $$ (15)
    $$ {{\bf{\ddot u}}_{\bf{s}}} = {{\bf{G}}_{\bf{0}}}v{}_0\Delta t + {{\bf{G}}_{\bf{0}}}{\bf{G}}{u_0} + {{\bf{G}}^{\bf{2}}}{{\bf{u}}_{\bf{s}}} $$ (16)

    其中,$ {{\bf{u}}_{\bf{s}}} = {({u_1}, {u_2}, \ldots \ldots, {u_m})^{\rm{T}}}$,$ {{\bf{G}}_{\bf{0}}} = ({a_{10}}, {a_{{\rm{2}}0}}, \ldots \ldots, {a_{m0}})$,$ {\bf{G}} = \left[ {\begin{array}{*{20}{c}} {{a_{11}}}& \cdots &{{a_{m1}}}\\ \vdots&\ddots&\vdots \\ {{a_{1n}}}& \cdots &{{a_{mm}}} \end{array}} \right]$。

    将式(15)、(16)代入式(10),并将已知量与未知量分离在等式的两侧,整理后得到如下方程:

    $$ {\bf{\hat k}}{{\bf{u}}_{\bf{s}}}={\bf{\hat p}} $$ (17)

    其中,${\bf{\hat k}} = {{\bf{G}}^{\bf{2}}} + 2\xi \omega \Delta t{\bf{G}} + {(\Delta t)^2}{\omega ^2}{\bf{I}} $,${\bf{\hat p}} = {(\Delta t)^2}{\bf{p}} - {u_0}(2\xi \omega \Delta t + {{\bf{G}}_{\bf{s}}}{{\bf{G}}_{\bf{0}}}) - {v_0}\Delta t{{\bf{G}}_{\bf{0}}} $,Im阶单位矩阵,$ {\bf{p}} = {({p_1}, {p_2}, \ldots \ldots, {p_m})^{\rm{T}}}$。

    式(17)为各时步内微分求积的基本方程,求解该方程可得到时步内各点的位移反应,再运用微分求积原理,可进一步求得各节点速度反应:

    $$ {{\bf{v}}_{\bf{s}}} = \frac{1}{{\Delta t}}({{\bf{G}}_{\bf{0}}}{u_0} + {\bf{G}}{{\bf{u}}_{\bf{s}}}) $$ (18)

    其中,${{\bf{v}}_{\bf{s}}} = {\left({{{\left. {\frac{{{\rm{d}}u}}{{{\rm{d}}t}}} \right|}_{\tau = {\tau _1}}}, {{\left. {\frac{{{\rm{d}}u}}{{{\rm{d}}t}}} \right|}_{\tau = {\tau _2}}}, \ldots \ldots, {{\left. {\frac{{{\rm{d}}u}}{{{\rm{d}}t}}} \right|}_{\tau = {\tau _m}}}} \right)^{\rm{T}}} $。

    将本时步末的位移$ {u_m}$和速度${\left. {\frac{{{\rm{d}}u}}{{{\rm{d}}t}}} \right|_{\tau = {\tau _m}}} $作为下一时步的初始位移和速度,从初始0时刻开始逐时步进行求解,可得到荷载作用的所有时刻的位移和速度反应。除了位移和速度外,工程上关心的加速度$ \frac{{{{\rm{d}}^2}u}}{{{\rm{d}}{t^2}}}$可通过运动方程(8)变形得到:

    $$ \frac{{{{\rm{d}}^2}u}}{{{\rm{d}}{t^2}}} = p(t) - 2\xi \omega \frac{{{\rm{d}}u}}{{{\rm{d}}t}} - {\omega ^2}u $$ (19)

    通过具体的算例验证上述微分求积法求解结构动力反应计算结果的精确性和可靠性。选择3种不同自振周期的单自由度体系,其频率范围大致覆盖低频、中频和高频,阻尼比都取为工程中常见的0.05,在上面施加不同频率的正弦荷载,体系的基本参数如表 1所示,荷载信息如表 2所示。

    表 1  体系的基本特性
    Table 1.  Basic characteristics of systems
    体系编号 自振周期Tn/s 阻尼比ξ
    1 0.68 0.05
    2 0.25 0.05
    3 0.08 0.05
    下载: 导出CSV 
    | 显示表格
    表 2  简谐荷载的信息
    Table 2.  Information of simple harmonic load
    荷载函数/N 周期/s 振幅/m 持时/s
    sin2πt 1 1 40
    sin10πt 0.2 1 40
    sin20πt 0.1 1 40
    下载: 导出CSV 
    | 显示表格

    用微分求积法求解3种体系在以上简谐荷载下的动力反应。为了方便计算,时步内采用均匀网格离散方案,时步的长度Δt取为简谐荷载的周期。由于1个周期的长度是简谐荷载的最小重复单元,这样取值不仅能方便地分析时步分段数m对数值稳定性和精度的影响,更能清楚地发现荷载周期与步长Δt取值的规律。

    首先,计算该时步长度下采用不同m所得到的体系的位移和速度反应,然后采用体系在简谐荷载下的解析解作为精确解进行校核。表 35列出了分段数m为20内的偶数时,简谐荷载激励下3种体系动力反应DQM解与精确解的平均相对误差。

    表 3  荷载周期1s的平均相对误差
    Table 3.  Average relative error with load period of 1s
    m 体系编号
    1 2 3
    位移/% 速度/% 位移/% 速度/% 位移/% 速度/%
    2 98.734 98.734 98.734 98.734 98.734 98.734
    4 88.218 81.849 99.415 389.850 53.532 1293.500
    6 42.402 32.547 15.926 4.728 10.979 9.609
    8 6.698 8.020 1.543 7.814 2.042 2.741
    10 1.623 2.687 0.907 2.338 1.036 0.349
    12 0.121 0.138 0.897 3.341 0.819 4.156
    14 0.008 0.009 0.935 5.153 0.716 1.290
    16 0 0 0.930 3.899 0.635 0.825
    18 0 0 1.005 3.292 0.577 2.28
    20 0 0 1.268 4.715 0.527 0.899
    下载: 导出CSV 
    | 显示表格
    表 4  荷载周期0.2s的平均相对误差
    Table 4.  Average relative error with load period of 0.2s
    m 体系编号
    1 2 3
    位移/% 速度/% 位移/% 速度/% 位移/% 速度/%
    2 99.749 99.749 99.749 99.749 99.749 99.749
    4 799.760 2460.400 403.870 411.780 108.540 225.480
    6 90.884 11.229 23.620 19.707 11.436 5.439
    8 6.052 15.762 2.102 2.135 1.068 2.520
    10 0.264 0.046 0.116 0.140 0.765 1.155
    12 0.009 0.022 0.005 0.004 0.807 3.325
    14 3681.600 256.240 251.760 2935.800 972.740 1334.400
    16 0 0 0 0 0.093 0.167
    18 0 0 0 0 0.010 0.013
    20 0 0 0 0 0.001 0.001
    下载: 导出CSV 
    | 显示表格
    表 5  荷载周期为0.1s的平均相对误差
    Table 5.  Average relative error with load period of 0.1 seconds
    m 体系编号
    1 2 3
    位移/% 速度/% 位移/% 速度/% 位移/% 速度/%
    2 99.875 99.875 99.875 99.875 99.875 99.875
    4 1620.200 1370.200 578.690 459.540 51.868 52.298
    6 175.570 9.444 65.345 7.903 6.80×10282 8.63×10282
    8 11.571 9.569 4.413 3.254 1.951 2.655
    10 0.498 0.046 0.204 0.044 0.142 0.266
    12 0.016 0.013 0.007 0.005 0.006 0.007
    14 27683 306.960 1322.400 246.430 234.600 3888.600
    16 0 0 0 0 0 0
    18 0 0 0 0 0 0
    20 0 0 0 0 0 0
    下载: 导出CSV 
    | 显示表格

    表 35可以看出,当时步分段数m很小时,DQM计算结果严重偏离精确解,但随着m的增大,除少部分点外,误差大致呈迅速减小的趋势,个别数据(如m=14时)反应误差巨大是由于算法在该时步分段下,时步长度取为荷载周期时出现了数值不稳定现象,初始误差随着逐时步推进不断放大,最后完全湮灭真实的结果。当m增大到10时,3种体系在不同周期的简谐荷载下的位移和速度反应都减小到了5%以内,对于实际工程已足够精确。为了更形象地描述这种精确程度,图 13给出了0—20s内荷载周期为1s下,m=10时DQM的位移计算结果与精确解,为了更清晰地进行对比,对每张图进行了局部的放大。

    图 1  简谐荷载激励下体系1的位移反应
    Figure 1.  Displacement response of system 1 under simple harmonic load
    图 2  简谐荷载激励下体系2的位移反应
    Figure 2.  Displacement response of system 2 under simple harmonic load
    图 3  简谐荷载激励下体系3的位移反应
    Figure 3.  Displacement response of system 3 under simple harmonic load

    从图中可以看出,DQM的计算结果与精确解几乎完全重合,充分说明了分段数m=10时,DQM计算动力反应足够精确,且数值稳定性能得到很好的保证。以地震荷载为例,中低频地震荷载的大部分周期一般在0.2s以上,偏于保守取为0.2s,分10段后节点间距为0.02s。而目前一般的地震地面加速度采样周期仅为0.005s或0.1s,节点间距取采样周期的几倍以上,都能得到精确的结果,计算精度较高。

    对比表 345可以发现,当m由10继续向上增大时(排除m=14时因失稳误差剧增的情况),误差基本上呈继续减小的趋势,很多情况下误差甚至降至1%或1‰以内。但这对工程实际已无太大的意义,反而会随着m的增大,不断地增大计算量。因此,对于结构在简谐荷载作用下的反应,采用均匀节点方案,即在1个荷载周期内取m=10,时步内插入9个内节点,既能达到土木工程领域内精度的要求,又能最大限度地减少计算工作量,是相对较优的时步节点数量。而且,以往DQM求解工程结构静力问题的大量经验表明,m取8—10可得到满意的计算结果(王鑫伟,1995),这与本文得出的m=10的取值较优也相符,这虽然是DQM解决静力问题的经验,但在1个周期内的动力问题与静力问题存在不少的联系,从而也可间接说明本文将m=10作为较优参数的合理性。

    表 35已在工程可接受的范还可以发现,对于简谐荷载激励下的反应,当m=10时,取时步长度Δt与荷载周期相等时,计算误差已在工程可接受的范围内。因此,计算体系在简谐荷载激励下的反应,在选定m=10的前提下,将时步长度和简谐荷载的周期保持一致即可。

    虽然简谐荷载是1种理想的荷载,但是对计算实际荷载激励下的反应具有重要的意义,这是因为工程中实际的动荷载都是一定持时的暂态荷载,都可以通过周期延拓表示成一系列简谐荷载的叠加。实际计算时,可以首先采用谐波分析,检测出其包含的不同周期的简谐波的幅值,得到该动力荷载的频谱;然后以最大幅值所对应的周期,即卓越周期为基准,确定等效周期;最后,取时步的长度为荷载的等效周期,将时步等分成10段进行计算,即可得到较精确的计算结果。

    本文将微分求积法引入结构动力反应的分析与计算,并通过数值算例得到如下的结论:

    (1)采用微分求积法求解结构在动力荷载激励下的反应合理可行。在较大的节点距离下依然能够得到精确的结果,计算精度高,且具有普适性,对不同自振周期的结构、不同频率的动力荷载都适用。

    (2)用DQM进行动力分析时,能一次性求得多个时刻的反应。相比传统的每次只能求1个时刻的逐步积分法,计算效率得到提高,计算成本也得到了降低。

    (3)在时步长度Δt一定时,DQM求解动力反应的计算精度和数值稳定性与时步分段数m有关。排除一些失稳飘移的情况,一般m越大,计算精度越高,但计算量也越大。综合考虑,对于均匀网格离散方案,实际计算时取m=10相对较优。

    (4)使用DQM进行实际结构动力反应分析时,时间步长Δt可选为动力荷载的等效周期,然后将各时步等分成10段来计算,这样可获得较满意的计算结果。

  • 图  1  穿过1#水泥试样的单频脉冲时域信号

    Figure  1.  Time domain signal through cement specimen 1#

    图  2  1#试样原始信号和使用脉冲反转技术得到的二次谐波频谱图

    Figure  2.  Fourier spectra of original signal and second harmonic signals extracted by pulse-inversion technique of cement specimen 1#

    图  3  4#试样在不同抗压强度下的尾波波形

    Figure  3.  Illustration for coda waves of sample 4# under different percentage of failure force

    图  4  实验装置及工作流程示意图

    Figure  4.  Illustration of experiment system and its workflow

    图  5  1#和4#试样随抗压强度的基频(a)与二次谐波(b)的衰减变化

    Figure  5.  Attenuation in amplitude of fundamental (a) and second order(b) harmonics of specimen 1# and 4# as a function of percentage of failure force

    图  6  2#和5#试样随抗压强度的基频(a)与二次谐波(b)的衰减变化

    Figure  6.  Attenuation in amplitude of fundamental (a) and second order (b) harmonics of specimen 2# and 5# as a function of percentage of failure force

    图  7  3#和6#试样随抗压强度的基频(a)与二次谐波(b)的衰减变化

    Figure  7.  Attenuation in amplitude of fundamental (a) and second order (b) harmonics of specimen 3# and 6# as a function of percentage of failure force

    图  8  1#和4#试样随抗压强度的尾波波速变化(a)与P波波速(b)

    Figure  8.  Velocity variation of tail wave and P wave speed of specimen 1# and 4# as a function of percentage of failure force

    图  9  2#和5#试样随抗压强度的尾波波速变化(a)与P波波速(b)

    Figure  9.  Velocity variation of tail wave and P wave speed of specimen 2# and 5# as a function of percentage of failure force

    图  10  3#和6#试样随抗压强度的尾波波速变化(a)与P波波速(b)

    Figure  10.  Velocity variation of tail wave and P wave speed of specimen 3# and 6# as a function of percentage of failure force

    表  1  水泥试样的物理参数

    Table  1.   Physical parameters of cement specimens

    试样编号 描述状态 尺寸/cm 水:水泥:骨料颗粒(质量) 粒径/cm 传播时间/µs 波速/m·s-1
    1 完整状态 4×4×8 0.4:1:1 0.1—0.2 9.94 4024
    2 完整状态 4×4×8 0.4:1:1 0.3—0.6 9.69 4124
    3 完整状态 4×4×8 0.4:1:1 0.8—1.2 9.46 4224
    4 高温加热 4×4×8 0.4:1:1 0.1—0.2 10.84 3687
    5 高温加热 4×4×8 0.4:1:1 0.3—0.6 10.65 3753
    6 高温加热 4×4×8 0.4:1:1 0.8—1.2 10.30 3883
    下载: 导出CSV
  • 陈小佳, 沈成武, Jacobs L. L., 2007.一种基于非线性超声谐波幅值比的微裂缝探测方法.武汉大学学报(工学版), 40(6):61-65. http://industry.wanfangdata.com.cn/dl/Detail/Periodical?id=Periodical_whsldldxxb200706013
    宋丽莉, 葛洪魁, 郭志伟等, 2012a.利用多次散射波监测介质性质变化的试验研究.岩石力学与工程学报, 31(4):713-722. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=yslxygcxb201204009
    宋丽莉, 葛洪魁, 梁天成等, 2012b.小应力扰动下岩石弹性波速变化的波形检测.中国石油大学学报(自然科学版), 36(4):60-65. http://d.wanfangdata.com.cn/Periodical_sydxxb201204011.aspx
    谢凡, 葛洪魁, 郭志伟, 2013.利用超声二次谐波测试水泥试样的非线性系数.机械工程学报, 49(14):9-15. http://d.wanfangdata.com.cn/Periodical_jxgcxb201314002.aspx
    周正干, 刘斯明, 2011.非线性无损检测技术的研究、应用和发展.机械工程学报, 47(8):2-11. https://www.wenkuxiazai.com/doc/3c3488975ef7ba0d4b733b1a-4.html
    朱金林, 刘晓宙, 周到等, 2009.声波在有裂纹的固体中的非经典非线性传播.声学学报, 34(3):234-241. http://xueshu.baidu.com/s?wd=paperuri%3A%28590f97add7fce8ab1465369066629b88%29&filter=sc_long_sign&tn=SE_xueshusource_2kduw22v&sc_vurl=http%3A%2F%2Fkns.cnki.net%2FKCMS%2Fdetail%2Fdetail.aspx%3Ffilename%3Dxiba200903009%26dbname%3DCJFD%26dbcode%3DCJFQ&ie=utf-8&sc_us=17138570087307439368
    Bruneau M., Potel C., 2009. Materials and acoustics handbook. Hoboken, NJ, USA:Wiley Online Library.
    Brunet T., Jia X. P., Johnson P. A., 2008. Transitional nonlinear elastic behaviour in dense granular media. Geophysical Research Letters, 35(19):L19308. doi: 10.1029/2008GL035264
    Clarke D., Zaccarelli L., Shapiro N. M., et al., 2011. Assessment of resolution and accuracy of the Moving Window Cross Spectral technique for monitoring crustal temporal variations using ambient seismic noise. Geophysical Journal International, 186(2):867-882. doi: 10.1111/j.1365-246X.2011.05074.x
    Frouin J., Sathish S., Matikas T. E., et al., 1999. Ultrasonic linear and nonlinear behavior of fatigued Ti-6Al-4V. Journal of Materials Research, 14(4):1295-1298. doi: 10.1557/JMR.1999.0176
    Hadziioannou C., Larose E., Coutant O., et al., 2009. Stability of monitoring weak changes in multiply scattering media with ambient noise correlation:laboratory experiments. The Journal of the Acoustical Society of America, 125(6):3688-3695. doi: 10.1121/1.3125345
    Jhang K. Y., 2009. Nonlinear ultrasonic techniques for nondestructive assessment of micro damage in material:a review. International Journal of Precision Engineering and Manufacturing, 10(1):123-135. doi: 10.1007/s12541-009-0019-y
    Jia X. P., Laurent J., Khidas Y., et al., 2009. Sound scattering in dense granular media. Chinese Science Bulletin, 54(23):4327-4336. doi: 10.1007/s11434-009-0609-1
    Kim J. Y., 2006. Experimental characterization of fatigue damage in a nickel-base superalloy using nonlinear ultrasonic waves. The Journal of the Acoustical Society of America, 120(3):1266-1273. doi: 10.1121/1.2221557
    Larose E., Hall S., 2009. Monitoring stress related velocity variation in concrete with a 2·10-5 relative resolution using diffuse ultrasound. The Journal of the Acoustical Society of America, 125(4):1853-1856. doi: 10.1121/1.3079771
    Mehta P. K., 1986. Concrete:Structure, properties and materials. Englewood Cliffs, NJ, USA:Prentice-Hall.
    Payan C., Garnier V., Moysan J., et al., 2009. Determination of third order elastic constants in a complex solid applying coda wave interferometry. Applied Physics Letters, 94(1):011904. doi: 10.1063/1.3064129
    Schurr D. P., Kim J. Y., Sabra K. G., et al., 2011. Monitoring damage in concrete using diffuse ultrasonic coda wave interferometry. AIP Conference Proceedings, 1335(1):1283-1290. doi: 10.1063/1.3592081
    Shah A. A., Ribakov Y., 2012. Damage detection in concrete using nonlinear signal attenuation ultrasound. Latin American Journal of Solids and Structures, 9(6):713-730. http://www.lajss.org/index.php/LAJSS/login?source=%2Findex.php%2FLAJSS%2Farticle%2Fview%2F407%2F973
    Shah A. A., Ribakov Y., Zhang C., 2013. Efficiency and sensitivity of linear and non-linear ultrasonics to identifying micro and macro-scale defects in concrete. Materials & Design, 50:905-916. https://www.sciencedirect.com/science/article/pii/S0261306913002872
    Snieder R., Grêt A., Douma H., et al., 2002. Coda wave interferometry for estimating nonlinear behavior in seismic velocity. Science, 295(5563):2253-2255. doi: 10.1126/science.1070015
    Tournat V., Gusev V. E., 2010. Acoustics of unconsolidated "model" granular media:an overview of recent results and several open problems. Acta Acustica united with Acustica, 96(2):208-224. doi: 10.3813/AAA.918271
  • 期刊类型引用(0)

    其他类型引用(1)

  • 加载中
图(10) / 表(1)
计量
  • 文章访问数:  71
  • HTML全文浏览量:  49
  • PDF下载量:  2
  • 被引次数: 1
出版历程
  • 收稿日期:  2017-05-05
  • 刊出日期:  2018-03-01

目录

/

返回文章
返回