動三軸試驗滯回曲線的橢圓擬合分析1

毛遠(yuǎn)鳳 葉慶東 沈宇鵬

?

動三軸試驗滯回曲線的橢圓擬合分析1

毛遠(yuǎn)鳳1）葉慶東1）沈宇鵬2）

1）中國地震局第一監(jiān)測中心，天津 3001802） 2）北京交通大學(xué)土木工程學(xué)院，北京 100044

采用正弦波形循環(huán)荷載動三軸試驗獲取土的動剪切模量和阻尼比的數(shù)據(jù)處理過程中通常存在兩個困難：一是因為試驗中存在各種噪聲使滯回曲線橢圓形態(tài)不明顯；二是橢圓擬合會因方法不當(dāng)出現(xiàn)不收斂或者誤差大等情況。為了在一定程度上克服第一個困難，本文將濾波技術(shù)引入到對應(yīng)力應(yīng)變時間序列的處理中來；為了更好地擬合滯回曲線橢圓，本文結(jié)合了主成份分析技術(shù)與橢圓的幾何擬合方法，通過計算橢圓焦點位置、長半軸長度等來確定擬合橢圓，結(jié)果表明，該方法容易線性化且相對穩(wěn)定，值得在動三軸數(shù)據(jù)處理中進(jìn)行推廣。

動三軸試驗 滯回曲線 橢圓擬合 幾何方法 濾波技術(shù)

引言

土的動剪切模量比和阻尼比是土動力學(xué)特征的重要參數(shù)（Hardin等，1968），在工程場地土層地震反應(yīng)分析和地震安全性評價工作中不可缺少（郭婷婷等，2016）。通過描述土在動力荷載作用下的動應(yīng)力應(yīng)變關(guān)系（動本構(gòu)關(guān)系）模型，將動剪切模量和阻尼比與動剪應(yīng)變幅值的函數(shù)關(guān)系具體化。Seed等（1970）首先給出了砂土和黏性土的動剪切模量和阻尼比與動剪應(yīng)變幅值的關(guān)系曲線。Hardin等（1972）給出了動剪切模量與動剪應(yīng)變的關(guān)系（即Hardin-Drnevich雙曲線模型）。Martin等（1982）在Hardin-Drnevich模型的基礎(chǔ)上進(jìn)行了改進(jìn)，提出了三參數(shù)的Davidenkov模型和具有冪次形式的阻尼比擬合公式。陳國興等（2005）基于Davidenkov模型，采用上限剪應(yīng)變幅值作為分界點，推導(dǎo)出了修正的Davidenkov模型阻尼比計算公式。在所有模型中，Hardin-Drnevich模型需要求解的參數(shù)最少，只需要2個參數(shù)，因而較為簡便，目前應(yīng)用較為廣泛。下文中討論動剪切模量和阻尼比與動剪應(yīng)變幅值關(guān)系就基于此模型。

土的動剪切模量和阻尼特性測定是地震等動荷載作用下土工建筑物設(shè)計和計算的基本依據(jù)（趙紅芬等，2009）。通常，為獲得土體的動剪切模量和阻尼比，最直接的方法是進(jìn)行土動力學(xué)試驗。在土的動力學(xué)特性與巖土地震工程試驗研究中最常用的方法主要是動三軸和自振柱試驗兩種方法。其中，動三軸試驗主要用于在室內(nèi)測定粘性土的動強度、各類土的動模量和阻尼比、砂土抗液化強度等相關(guān)參數(shù)（吳世明等，2001；張濤，2004；申權(quán)等，2013；郭可骍，2015）?；谇叭藢ν羷恿μ匦缘姆治?，同時考慮到試驗測試及數(shù)據(jù)采集系統(tǒng)的特征，試驗荷載一般采用正弦波形循環(huán)荷載（南京水利科學(xué)研究院，2003）。在動三軸試驗中，數(shù)據(jù)處理是保證試驗成果有效產(chǎn)出的關(guān)鍵，而應(yīng)力應(yīng)變滯回曲線橢圓擬合則是判斷試驗成果準(zhǔn)確性的決定因素。

1 循環(huán)應(yīng)力動三軸試驗

圖1 確定動彈性模量和阻尼比示意圖

2 擬合滯回曲線橢圓的方法分析

試驗采用正弦波循環(huán)荷載，即動應(yīng)力可以表示為：

對應(yīng)的動應(yīng)變可以表示為：

橢圓擬合的方法有很多種，如蒙特卡羅法、非線性最小二乘法（許正文等，2008；閆蓓等，2008）、參數(shù)擬合法（陳基偉，2007a；2007b）和幾何法（彭青玉，2003）等。蒙特卡羅法也稱統(tǒng)計模擬法，需要借助計算機生成大量樣本并對其特征進(jìn)行統(tǒng)計分析，通常較為費時；非線性最小二乘法直觀明了，但是需要求解5個參數(shù)，同時還要求二元二次函數(shù)滿足橢圓的約束條件，事實上是一個二次規(guī)劃問題，容易因陷入局部極小而不能很好地擬合觀測數(shù)據(jù)；參數(shù)擬合法處理已知初相的問題較為方便，但對本文涉及的問題無法線性化，求解存在困難。彭青玉（2003）根據(jù)橢圓到兩定點距離之和為常數(shù)的幾何定義出發(fā)，提出了橢圓擬合的幾何方法。與傳統(tǒng)方法相比，該方法直接基于幾何定義，因此相對穩(wěn)定，并且很容易線性化，可操作性強，下文中將對該方法進(jìn)行詳細(xì)介紹。

圖2 應(yīng)力、應(yīng)變初始（（a）、（b））及濾波后（（c）、（d））的時間曲線

橢圓擬合的幾何方法中，需要給定半長、短軸以及長軸與橫坐標(biāo)的夾角。我們的做法是首先將濾波后的應(yīng)變-應(yīng)力曲線畫在一張圖上，然后手動刪去明顯異常的點，如圖3。

理想情況下橢圓上的點到兩焦點距離之和為常數(shù)2，實際中由于誤差的存在，橢圓擬合方程表示為：

圖4 橢圓擬合迭代前（a）后（b）變化示意圖

3 試驗結(jié)果

表1 試樣動三軸試驗結(jié)果，不同動剪應(yīng)變下的及

圖5 不同動剪應(yīng)變下的及動阻尼比曲線

從圖5可以看出，本文中橢圓擬合的幾何方法能夠很好地擬合曲線，是能夠應(yīng)用于工程實踐的實用方法。

4 結(jié)論

通過上述試驗和分析，得到了不同動剪應(yīng)變下該粘土樣的動剪切模量和阻尼比隨剪應(yīng)變變化的關(guān)系曲線?？梢钥闯觯疚倪M(jìn)行橢圓擬合的幾何方法是一種相對穩(wěn)定、可操作性強的方法。在實際工程中，采用人機交互可有效地控制數(shù)據(jù)質(zhì)量。將該方法應(yīng)用于同一批其他試驗土樣的滯回曲線擬合，均得到了較好的結(jié)果，進(jìn)一步證明了該方法可以確保土動三軸試驗數(shù)據(jù)結(jié)果的可靠性以及各土動力學(xué)特性參數(shù)取值的準(zhǔn)確性。

陳國興，莊海洋，2005．基于Davidenkov骨架曲線的土體動力本構(gòu)關(guān)系及其參數(shù)研究．巖土工程學(xué)報，27（8）：860—864．

陳基偉，2007a．橢圓直接擬合算法研究．工程勘察，（6）：49—51．

陳基偉，2007b．工程測量中一類參數(shù)曲線的擬合．大地測量與地球動力學(xué)，27（1）：100—103．

郭可骍，2015．非飽和黃土的動三軸試驗研究．西安：長安大學(xué)．

郭婷婷，秦梅梅，2016．土動三軸試驗參數(shù)選取的理論分析與計算．工程抗震與加固改造，38（2）：144—149．

李松林，1990．動三軸試驗的原理與方法．北京：地質(zhì)出版社．

南京水利科學(xué)研究院土工研究所，2003．土工試驗技術(shù)手冊．北京：人民交通出版社．

彭青玉，2003．木星土星邊緣的橢圓擬合．云南天文臺臺刊，（4）：43—48．

申權(quán)，李明俊，蔣文明等，2013．動三軸試驗測試土阻尼的影響因素與不足．江西科學(xué)，31（1）：84—89．

孫靜，袁曉銘，孫銳，2004．土動剪切模量和阻尼比的推薦值和規(guī)范值的合理性比較．地震工程與工程振動，24（2）：125—133．

吳世明，周健，楊挺，2001．土動力學(xué)理論與計算．北京：中國建筑工業(yè)出版社．

許正文，姚連璧，2008．基于穩(wěn)健估計的直接最小二乘橢圓擬合．大地測量與地球動力學(xué)，28（1）：77—80．

閆蓓，王斌，李媛，2008．基于最小二乘法的橢圓擬合改進(jìn)算法．北京航空航天大學(xué)學(xué)報，34（3）：295—298．

楊淑瑩，張樺，2015．模式識別與智能計算-MATLAB技術(shù)實現(xiàn)．3版．北京：電子工業(yè)出版社．

張茹，何昌榮，費文平等，2006．固結(jié)應(yīng)力比對土樣動強度和動孔壓發(fā)展規(guī)律的影響．巖土工程學(xué)報，28（1）：101—105．

張濤，2004．中型動三軸儀的研制及城市垃圾土動力特性試驗研究．大連：大連理工大學(xué)．

趙紅芬，何昌榮，王莉娜，2009．動模量阻尼的動三軸試驗研究．路基工程，（4）：158—160．

Hardin B. O., Black W. L., 1968. Vibration modulus of normally consolidated clay. American Society of Civil Engineers, 94(2): 353—370.

Hardin B. O., Drnevich V. P., 1972. Shear modulus and damping in soils. Journal of the Soil Mechanics and Foundations Division, 98(7): 667—692.

Martin P. P., Seed H. B., 1982. One-dimensional dynamic ground response analyses. Journal of the Geotechnical Engineering Division, 108(7): 935—952.

Seed H. B., Idriss I. M., 1970. Soil moduli and damping factors for dynamic response analyses. Report No EERC70—10. Berkeley: Earthquake Engineering Research Center, University of California Berkeley.

Analysis on Ellipse Fitting of Hysteresis Curve inDynamic Rriaxial Test

Mao Yuanfeng1), Ye Qingdong1)and Shen Yupeng2)

1) The First Monitoring and Application Center, China Earthquake Administration, Tianjin 300180, China 2) School of Civil Engineering, Beijing Jiaotong University, Beijing 100044, China

There exist two difficulties in the data processing of dynamic tri-axial tests when we use sinusoidal cyclic loading to obtain dynamic shear modulus and damping ratio. One difficulty is that the elliptic hysteresis curve is not clear enough due to various noises, and another one is that the procedure of fitting ellipse is divergent or the results of fitting with larger uncertainty because selected ellipse fitting method is improper. In this study, in order to overcome these difficulties to some extents, we processed time series of dynamical stress and stain with filtering technology. To fit the elliptic hysteresis curve well, we calculated the fitted ellipse by determining the two focuses and the length of the long axis of ellipse which we obtained by combining principal component analysis with ellipse geometrical fitting method. The result shows that this method is easy to linearize the nonlinear ellipse fitted problem and relatively stable, and it is worth popularizing in the dynamical tri-axial test data processing.

Dynamic tri-axial test; Hysteresis curve; Ellipse fitting; Geometric method; Filtering technology

10.11899/zzfy20170402

中國地震局第一監(jiān)測中心科技創(chuàng)新主任基金（FMC2017002）與中國地震局震情跟蹤項目（2017010124）共同資助

2017-03-06

毛遠(yuǎn)鳳，女，生于1987年。碩士。主要研究領(lǐng)域：巖土工程、地震科研。E-mail：myf4440@163.com

毛遠(yuǎn)鳳，葉慶東，沈宇鵬，2017．動三軸試驗滯回曲線的橢圓擬合分析．震災(zāi)防御技術(shù)，12（4）：743—750．