聊聊動態強度和損傷演化*

王禮立,胡時勝,楊黎明,董新龍,王 暉

(1.寧波大學省部共建教育部沖擊與安全工程重點實驗室,浙江寧波315211; 2.中國科學技術大學中國科學院材料力學行為和設計重點實驗室,安徽 合肥230027; 3.寧波市中醫醫院王暉工作室,浙江寧波315000)

聊聊動態強度和損傷演化*

王禮立1,胡時勝2,楊黎明1,董新龍1,王 暉3

(1.寧波大學省部共建教育部沖擊與安全工程重點實驗室,浙江寧波315211; 2.中國科學技術大學中國科學院材料力學行為和設計重點實驗室,安徽 合肥230027; 3.寧波市中醫醫院王暉工作室,浙江寧波315000)

材料強度在傳統上常理解為材料在外載荷下抵抗流動/變形和破斷的能力。由流變階段到貌似突發的破斷,其實源于一個隱含的應變率/時間相關的損傷演化過程。動態損傷演化研究的難點在于損傷與流變總是耦合在一起發展的。研究發現,熱激活損傷演化模型可成功描述材料宏觀損傷的動態演化。在此基礎上,從實測的含損傷演化的表觀應力應變曲線,可將兩者解耦分開,并可確定各自相關的材料參數。這一思路可推廣到中醫脈診的客觀化研究,通過脈搏波信息定量反演脈搏波系統的正常及病態本構關系,可診斷生命體偏離正常狀態的程度(病情),這可視為一種廣義的損傷演化和強度問題。上述思路還可推廣到地震預報研究,即“對地球把脈”。與加卸載響應比理論相結合,通過相鄰3處的地震波信息來反演地球相關板塊含損傷演化的非線性載荷-響應曲線,再區分出損傷演化程度,將有利于改進地震預報,這可視為另一種廣義的損傷演化和強度問題。

動態強度;損傷演化;材料動力學;脈搏波;地震預報

固體的強度和流體的湍流常常并列為力學研究中的兩大經典難題。

根據研究對象的不同,對強度往往有不同的理解。按中國大百科全書《力學卷》[1],強度定義為“材料或結構抵抗外力的能力”。在外力作用下材料或結構的力學響應主要表現為流動/變形過程、直至破壞,所以強度可以更具體地理解為材料或結構在外力作用下抵抗廣義破壞(流動/變形和斷裂)的能力。

材料強度和結構強度又有什么區別?是什么關系呢?早期對它們常常不加區分。例如,大學課程《材料力學》教材可以追溯到20世紀30年代出版的鐵木辛柯的經典著作,原名是《Strength of materials》[2],直譯是材料強度。但實際上只有小部分內容講授材料的力學性質和強度,大部分內容主要講授代表性結構元件(桿、軸、梁等)受力時的應力應變分析,不如稱為“結構元件強度分析”或“構件力學”更合適。課程中真正涉及材料強度的是四大經典強度準則,即最大拉應力強度準則(第一強度理論)、最大伸長線應變強度準則(第二強度理論)、最大切應力理論(第三強度理論)和最大畸變比能理論(第四強度理論)。更一般地,一個可供定量分析的廣義強度準則可概括為:

上式將力學特征量Σ和材料特征量Σc相聯系,其中Σ屬于結構響應,要依靠力學家對結構力學場進行分析得到(但在力學分析中要依賴材料本構關系);而Σc屬于材料響應,要依靠材料學家對材料進行實驗研究確定(但在材料實驗中要依賴力學分析)。強度準則建立了結構響應和材料響應之間的聯系,兩者又是互相耦合的。一旦滿足不等式(1),即結構力學場的最大特征值達到/超過材料強度特征值,就判定為結構破壞/失效了。

與準靜載荷下的力學問題相比較,動載荷下的力學問題必須考慮2個主要的動態效應:慣性效應和應變率效應。前者導致了有關應力波傳播的研究,后者則導致了有關材料動態力學性能(材料動力學)的研究。這樣,沖擊動態下的廣義強度準則應將不等式(1)改寫為:

式中:Σ(t)是計及應力波傳播的非定常動態力學場特征量,現在是時間t的函數;則是計及應變率效應的表征材料動態強度特征的臨界參量,現在是應變率的函數。材料動態強度是《材料動力學》的核心問題,一般地可歸結為兩方面,即:表征材料流動/變形規律的動態本構關系和表征材料破斷的動態破壞規律的研究。這可以用圖1來示意[3],圖中σeff和εeff分別為基于應力偏量第二不變量J2所定義的所謂“有效應力”及相對應的“有效應變”,T*=Tln(ε·0/ε·)是所謂的率-溫等效參數,用以計及應變率效應和溫度效應之間等效性。圖1在不同T*值下的應力應變曲線表征材料率/溫相關的流變特性,而以右側的破壞極限包絡線為界;正是這一右界表征了材料率/溫相關的破壞特征量Σc。人們會對圖1提出一個問題:由動態本構關系所控制的流變階段怎么變為貌似突發的動態破壞呢?

圖1 材料的率-溫相關的動態力學行為之示意Fig.1 Rate-temperature-dependent behaviors of materials

其實,材料的動態破壞,不論其具體機制如何,實際上都有一個發生-發展的時間過程,本質上是時間/速率相關的過程,而并非瞬時事件。設以Tf來表征材料的內稟破壞特征時間,當Tf與沖擊/爆炸載荷特征時間Tl可比較時,就應計及材料動態破壞過程的時間(速率)相關性;只有當Tf?Tl時才可以按準靜態處理,把材料破壞近似看作瞬態事件。

從細觀角度看,材料的動態破壞是一個不同形式的微損傷(微裂紋、微空洞、微剪切帶等)以有限速率演化(成核-生長-連通)的時間過程。因此,從機理上說,材料動態破壞的研究已經離不開對損傷動態演化規律的研究。

問題的復雜性在于:一方面,損傷是隨流變過程而發展的,損傷的演化依賴于材料所經受的應力、應變、應變率等材料本構參量;另一方面,損傷演化及其引發的弱化效應,必將反過來影響材料外在表現的力學行為,體現為,材料表觀的本構關系其實同時包含著材料內在的本構行為,以及宏觀上不易直接觀察到的損傷演化。因此,損傷動態演化的研究往往只能與表觀率型動態本構關系的研究相耦合地進行,再設法將兩者加以區分。計及動態損傷演化的率型本構關系的研究,以及動態破壞準則的研究,已成為力學家和材料科學家們共同的前沿研究課題。其核心問題之一乃是下一節將要討論的宏觀連續損傷的率型動態演化律。

1 基于熱激活機制的連續損傷的率型動態演化律和動態破壞準則

研究者們從微觀和細觀角度對于損傷演化已經開展了大量的研究,提供了損傷演化的物理機制。但為了工程應用,最后還是要回到宏觀層次。通常在如下定義的宏觀連續損傷D[4]的基礎上:

建立可供工程應用的宏觀損傷演化模型,式中σ0為無損傷材料的應力,σ為含損傷材料的表觀應力。

一般容易認識到損傷隨流變過程而演化,因而D顯然是應變ε的函數,D=D(ε)。而大量動態實驗還進一步表明,各種形式的損傷在沖擊載荷下的演化都同時依賴于應變和應變率,即D=D(ε,ε·)。這方面研究的細節可參考文獻[5-8],不在此處詳述。下面將集中討論宏觀連續損傷D=D(ε,ε·)基于熱激活機制的動態損傷演化模型[9-11]。

聯系到材料的率型本構關系,在微觀機制上常用晶體缺陷(微觀損傷)—位錯的熱激活運動來解釋:

式中:US為與作用應力σ相關的熱激活能為頻率因子,k為Boltzman常數,T為絕對溫度。

再聯系到材料細觀損傷研究中,眾所熟知的Curran-Shockey-Seaman的細觀損傷成核與生長(nucleation and growth,NAG)模型中[12-15],其成核過程也參照S.N.Zhukov等的研究成果[16],建立在細觀損傷熱激活成核的機制上,即如果以N表示單位體積內各種尺寸裂紋的數目,則成核速率:

式中為頻率因子,σn0為成核應力閾值,σ1為表征成核的應力敏感性的材料常數。注意式(5)已經假定成核的熱激活能與外加應力成正比。

對于式(6),問題的關鍵在于如何確定熱激活能UD與作用應力σ之間的關系。其實,不論就流變過程中位錯演化的熱激活方程(4)而言,還是就損傷演化的熱激活方程(6)而言,難點和關鍵都是如何確定熱激活能與作用應力之間的關系。

對于同一材料,設想UD和US應該有內稟的聯系。暫不具體考慮它們如何分別依賴于作用應力σ的具體函數形式,而作為一級近似,可設US和UD之間有正比關系:

式中:λ為材料參數。把式(4)和(6)代入式(7),經演算后可得:

積分后可得應變率顯式相關的率型損傷演化律:

對于恒定應變率的過程,且設損傷演化存在某個應變閾值εth,則對式(8c)積分后可得:

式(9a)顯式地刻畫了連續損傷的演化同時依賴于應變和應變率。在更一般的情況下,D與應變之間可能有非線性關系,則式(9a)可推廣為如下更一般的形式:

式中:κ≥1為材料常數。在式(9)的基礎上,動態破壞準則D≥Dc可以表達為如下的率相關形式:

式中:Dc為與材料動態破壞的臨界狀態相對應的臨界損傷值。式(10)表明,對于每一給定的臨界損傷Dc,有一條臨界曲線,不同的應變率載荷下將有不同的破壞應變(雙變量破壞準則),這與圖1給出的右界曲線一致。顯然,由于λ的不同將會出現以下3種情況:

(1)如果λ＞1,則隨著應變率的提高,破壞應變減小,呈現所謂的“沖擊脆化”。

(2)如果λ＜1,則隨著應變率的提高,破壞應變增大,呈現所謂的“沖擊韌化”。

(3)如果λ=1,則式(10)簡化為臨界應變準則,εc=εth+Dc/KD。

綜上所述,按照連續損傷演化的TADE模型,損傷演化律歸結為如何通過實驗確定εth、KD、λ和κ;而動態破壞準則歸結為如何確定Dc。

2 計及動態損傷演化的率型本構關系

考慮到損傷演化與本構流變演化是耦合的,各種動態實驗迄今又無法直接測量式(9)和式(10)所包含的損傷演化參量,實際上必須耦合地開展“計及動態損傷演化的率型本構關系”的研究,再對兩者加以區分。對此,下面通過兩個代表性實例來進一步加以討論。

(1)材料率型本構方程為已知時

以有機玻璃PMMA為例[17],已知其無損傷時的非線性黏彈性本構行為滿足朱-王-唐方程(簡稱Eq.ZWT),當計及損傷演化時,則按式(3)有:

式中:σe(ε)描述非線性彈性平衡響應,σm、m和正整數n均為材料參量;E1和θ1分別為描述低頻響應的Maxwell單元的彈性常數和松弛時間;E2和θ2則分別為描述高頻響應的Maxwell單元的彈性常數和松弛時間。PMMA的ZWT方程的材料參數可在無損傷條件下測知,在應變率10-4～103s-1的范圍內,實驗所確定的參數值為:σm=91.8 MPa,m=22.3,n=4,E1=0.897 GPa,θ1=15.3 s,E2=3.07 GPa,θ2=95.4μs。如圖2所示,理論曲線與實驗數據能令人滿意地符合,除了在高應變率下隨著應變在6%附近進一步增大,未損傷PMMA的理論預示偏離實驗曲線,并出現應變軟化(dσ/dε＜0);相應地則在透明的PMMA試樣中觀察到隨應變增大而增多的微裂紋。因而可把這種應變軟化歸結為損傷弱化所致的本構失穩。未損傷ZWT方程的理論預示與實驗曲線的差值正代表損傷弱化所致的應力差(σ0-σ)。由此可按式(3)確定不同應變率和應變下的D,并進而通過數值擬合方法來確定TADE損傷演化方程的材料參數,得到KD=1.82,λ=1.17,κ=1,εth=0.06。再由實測的試樣破壞應變值與式(10),即可確定臨界損傷Dc=10.2%。

圖2 PMMA不同應變率下的理論曲線與實驗數據的對比Fig.2 Comparisons between the experimental data and theoretical predictions for PMMA

不同應變率下的理論預示與實驗數據的對比,如圖2所示?？梢娫趶V泛應變率范圍內,不僅加載曲線部分吻合得很好,反映了高應變率下出現的損傷弱化本構失穩,而且卸載曲線部分也吻合得很好;還顯示雖然不同應變率下的破壞應變值不同,但破壞點的臨界損傷值Dc幾乎接近恒值。由此充分支持了TADE損傷演化律和動態破壞準則的有效性。

(2)材料率型本構方程為隱函數時

上一個實例適用于材料的本構方程類型為已知的情況。當材料率型本構方程類型未知時,材料在損傷演化前后(以門檻應變εth為界),其本構關系一般可以分別表達為:

以下分別簡稱為“情況1”和“情況2”。

這時如何確定材料本構響應諸變量與損傷演化諸變量的過程,相當于系統辨識的過程,即根據已知的輸入、輸出數據,選擇一個系統模型,采用優化方法使系統模型逼近真實系統。在系統辨識的各種方法中,人工神經網絡(artificial neural network,ANN),尤其是反向傳播神經網絡(back-propagation neural networks,BP神經網絡),特別適合于處理式(12)和式(13)這類非線性問題,已被用來確定損傷演化參數等,而不必預先假定材料本構關系和損傷演化律的函數形式[18-22]。

以聚丙烯-尼龍共混高聚物(polypropylene-polyamide(PP/PA)polymer blend)為例[22],其損傷演化門檻應變εth值可由“損傷凍結法”預先測得,即對經歷過不同沖擊加載條件的試樣,事后在靜載下測量其損傷程度來初估,εth≈6%。BP神經網絡分別按“情況1”和“情況2”兩種情況來訓練。在“情況1”時,取SHPB(split Hopkinson pressure bar)實驗測得的作為輸入,而取相應的σ(t)作為輸出。而在“情況2”時,本應取和D(t)作為輸入,但在SHPB實驗過程中迄今尚未解決D(t)的直接測量技術,因此改為取可測的的反函數,即取SHPB實驗過程中可測的和t-1(D)作為輸入,而取相應σ(t)作為輸出。

在“情況1”下所得代表性結果(應變率為1.22×103s-1)如圖3(a)中點線所示。作為比較,圖中還給出了實驗曲線(實線)。圖中的縱坐標以歸一化應力表示,此處σmax是全部實驗中出現的最大應力值,因而有

由圖3(a)可見,在ε＜εth的范圍內,BP神經網絡預示曲線與實驗曲線吻合得很好。在ε＞εth的范圍內,兩曲線的偏差則隨應變的增大而增大,顯示損傷演化導致的弱化效應。事實上,當計及損傷演化后,采用BP神經網絡在“情況2”下所得代表性結果(應變率為1.22×103s-1)如圖3(b)中點線所示。這時,BP神經網絡預示曲線與實驗曲線(實線)在整個應變范圍內都吻合得很好。

圖3 PP/PA共混高聚物BP神經網絡預示曲線與實驗曲線對比Fig.3 Comparisons between the BP neural networks predictions and experimental data for PP/PA polymer blend

由上述的識別結果,按照式(3)可最后確定連續損傷作為應變率和應變的函數以不同的恒應變率下D隨應變而演化的形式作圖,如圖4所示,與TADE損傷演化模型演化規律一致[22]。

由以上討論可知,材料動態強度更具體地應理解為材料在強動載荷下抵抗動態流變、動態損傷演化和動態破斷的能力。相應地,圖1應修正為圖5所示的計及損傷演化的材料率-溫相關的動態力學行為。圖中增加了一條表征損傷開始演化的門檻等效應變εth曲線,而右界曲線則理解為D=Dc曲線。從εth曲線開始,材料進入與損傷演化相耦合的流變過程,直到損傷達到臨界值(D=Dc)而最終導致破壞。

圖4 PP/PA共混高聚物由BP神經網絡確定的D=(ε)Fig.4 TheD=(ε,)determined by the BP neural networks for PP/PA polymer blend

3 由脈搏波來診斷病情(一種廣義損傷)

中醫脈診是通過脈搏波來診斷病人生命體偏離正常運行的狀態—病情,這可以看作一種廣義的損傷演化和強度問題。

在材料動力學研究中,研究者們已發展了一套“解反問題”的方法,即通過測量一系列應力波來反演材料含損傷演化的本構關系,然后如上節所述進一步區分材料本構關系和損傷演化律。依據同樣的思路,能不能對中醫脈搏波進行定量分析來反演和判定人體的損傷演化程度呢?用力學研究的語言說,中醫切脈通過脈搏波來診斷病情,都是在解人體健康狀態及其損傷的反問題。

中醫的整體觀點認為,脈搏波反映的是生命體整體狀況的綜合信息。據此我們把脈搏波系統理解為生命能量以波的形式在血液中傳播的整個系統,而不只是心血管等循環系統本身局部實體器官的信息。按照這一思路,為脈搏波系統建立一個等價的一維力學模型[23-24]?？刂品匠探M由3個守恒方程及系統的本構方程共同組成。其中,3個守恒方程即動量守恒、質量守恒和能量守恒方程分別為:

圖5 計及損傷演化的材料率-溫相關的動態力學行為的示意圖Fig.5 Rate-temperature dependent behavior of materials taking account of damage evolution

式中:X為物質坐標,t為時間;P為壓力,V為比容(與密度ρ互為倒數,ρ=1/V),v為質點速度,E為比內能,均為(X,t)函數。如果把式(15)代入式(16a),可以得到能量守恒方程的另一種形式:

意味著比內能E等于P-V曲線下方包圍的面積。

至于“本構方程”一詞,在醫學界不大熟悉。在力學界把反映材料力學性能的各力學量之間的內在關系稱為本構關系,而本構關系的數學表達形式稱為本構方程。在上述脈搏波系統的3個守恒方程中,有3個力學量與本構性能相關,即P、V和E(v是運動參量,而不是物質的本構參量),當考慮到率效應時,率相關本構方程為:

當忽略率效應并考慮E可以由P和V的關系確定(式16(b)),則有簡化的率無關本構關系:

式(14)～(17)共同組成脈搏波系統的控制方程組,其中的3個守恒方程是普適的,反映了各種脈搏波的共性,而系統本構方程則反映了不同脈搏波的特性,不同脈搏波的傳播特性主要取決于本構方程。

對于脈搏波的控制方程組(式(14)～(17)),由應力波理論可知[25]:

(1)脈搏波包含著壓力波P(X,t)、質點速度波v(X,t)、比容波V(X,t)和內能波E(X,t)等的傳播(對非線性波,它們間不存在線性比例關系)。脈搏波波速由本構方程P=P(V)的局部斜率決定,不是恒值(對于大擾動有可能形成沖擊波,將另行討論)?？梢娒}搏波是這些不同形式波動的總成。目前有些研究者只憑壓力傳感器單點測得的脈搏壓力波P(X,t)來分析,看來難以完整地模擬中醫脈診。

(2)脈搏現象包含著以v流動的血液和在血液中以C傳播的攜帶能量的脈搏波。前者是人眼易見的實體血液的“物質流”,后者是人眼不易見的以波的形式傳播的“能量流”,可以分別解釋為中醫的“血”和“氣”。這為中醫的氣提供了一個新的理解。血是氣傳播的物質載體(媒介),而氣是以波的形式傳播的推動血運行的能量。

(3)要注意區分血液質點的運動速度v(10-1m/s量級)和脈搏波的傳播速度C(100～101m/s量級)。脈搏波以比血液流速v快得多的傳播速度C,把心臟施加給循環血液的力學擾動(脈動載荷,即血液的壓力、容積、流速和能量脈動等)由近及遠地傳遞到生命體的各部分,不到一秒已傳遍全身?？梢?心臟并非像常識理解那樣猶如把血液直接泵向全身的泵,血液狀態的變化是通過脈搏波傳播而實現的。

關于脈搏波系統的本構方程,迄今未見研究報導,甚至還沒有人提出這一問題,而這正是我們要進一步研究的。中醫“切脈象”相當于根據脈診獲得的脈搏波信息,去解“反問題”,即反求系統本構方程。進而根據系統本構方程是否對健康狀態有偏離(病態)以及偏離的原因,是什么類型的病態等等,對病情(廣義損傷)做出診斷,對癥下藥,予以診治。在這個意義上,“切脈象”的診斷對象和根本目標,在于診斷各個具體生命體的脈搏波系統的本構方程的變化狀態,及各種類型脈搏波的變化狀態。推而廣之,不論中醫還是西醫,實際上大量工作都在捕捉廣義的血和氣的各種時空變化信息,然后去解“反問題”。

解一維脈搏波反問題的出發點是守恒方程。中醫“切脈象”時,手指主要感受的脈搏壓力波,我們就來討論這類情況下的反演。設用n個壓力計在n個物質坐標Xi(i=1,2,…,n)處測知了n個壓力波剖面P(Xi,t)。既然不難確定?P/?X,由動量守恒方程(14)就可求得?v/?t。注意到v的零初始條件:t=0時v=0,積分可得v(Xi,t)。既然不難求得?v/?X,由質量守恒方程(15)就可求得?V/?t。注意到V的零初始條件t=0時V=0,積分可得V(Xi,t)。由P(Xi,t)和V(Xi,t)消去t,最終可建立系統的本構關系P=P(V)。由不同相鄰Xi處的不同的P=P(V)可以確定率相關性,得出計及黏性效應的率相

·關本構關系P=P(V,V)。

那么測點數n最少應該是多少呢?注意到由v(Xi,t)求V(Xi,t)時,是通過?v/?X進行的,因而至少應該知道兩個物質質點處的v(Xi,t)波形(i=1,2);然而,2個v(Xi,t)波形是通過?P/?X進行的,因而至少應該知道3個相近物質質點處的P(Xi,t)波形(i=1,2,3)。所以,由實測P(Xi,t)波形反演系統本構關系時至少要測知3個相近物質質點處的。

令人驚訝的是,至今仍指導著中醫臨床實踐的“寸口三部九候診法”,一般均采用北宋醫家丁徳用的所謂“密排三指”之法來診斷寸、關、尺三處脈象。這與上述至少要測知3個相近物質質點處P(Xi,t)波形(i=1,2,3)的要求完全一致(三點律)。既證明了“寸口三部九候診法”的合理性,又令人深感我國古代中醫名家之高明。

由此可知,通過相鄰3點處的脈搏壓力波信息,可以反演脈搏波系統本構關系。健康人的本構關系相當于無損傷的正常本構關系,病人的本構關系則相當于含損傷的病態本構關系,對比兩者的差別,就可以定量地描述病態。

傳統中醫學中雖無脈搏波系統的本構方程之說,甚至于迄今還沒有人提出這一問題,但聯系到近年來引人注目的中醫體質學,就體質(body constitution)的重要性和地位而言,相當于前述量化的脈搏波系統中的本構關系。近代中醫體質學的開創人王琦[26-27]把人體生命過程中,在先天稟賦和后天獲得的基礎上所形成的形態結構、生理功能和心理狀態方面綜合的、相對穩定的固有特質稱之為“體質”,并提出9種基本類型體質:平和質、氣虛質、陽虛質、陰虛質、痰濕質、濕熱質、瘀血質、氣郁質和特稟質等。名老中醫王暉則將體質學說、陰陽五行、易理洛書等引入五行體質,形成了五行體質觀:木型體質、火型體質、土型體質、金型體質和水型體質。但還沒有定量化的描述。如果能夠把中醫的體質與本文量化的脈搏波系統本構方程之間建立定量的內在聯系,顯然是很有意義的,是值得進一步研究的。

4 由地震波來診斷震情(一種廣義損傷)

中醫把脈是通過脈搏波來診斷病人生命體偏離正常運行的狀態—病情(一種廣義的損傷演化和強度問題)。與此相類似,地震預報相當于給地球把脈,通過地震波來診斷地球板塊偏離安全而可能發生地震的狀態—震情,這是另一種廣義的損傷演化和強度問題。

尹祥礎正確地指出[30-31]:“地震的孕育過程就是震源區介質的損傷、演化并最終導致破壞的過程,這一過程主要是力學過程”;并相應地提出“加卸載響應比地震預測”新思路,試圖通過“加卸載響應比”(load-unload response ratio,LURR)的測量來把握震源區介質的損傷演化過程及預測其與破壞相對應的臨界狀態。

圖6 加卸載響應比概念之示意Fig.6 Schematics of the concept of load-unload response ratio

LURR理論的出發點是:損傷演化過程表現在應力應變曲線上,如圖6所示,加載時的變形模量小于卸載時的變形模量(或加載時的柔度大于卸載時的柔度),其差異反映了材料損傷或力學性質劣化的程度。

引入如下定義的響應率X(廣義柔度)和加卸載響應比Y:

式中:P和R分別為載荷和響應,如圖6所示;X+和X-分別為加載的響應率和卸載的響應率。當介質處于線彈性階段時,X+=X-,Y=1;到了損傷階段,由于X+＞X-,就有Y＞1。隨著損傷的增加,Y增大,直至臨近破壞,Y達到其峰值。因此,Y可定量刻畫震源區介質的損傷程度,其變化可以作為地震發生的前兆現象。

如何對尺度為上百公里甚至是上千公里的地殼板塊進行加載和卸載?方法之一就是利用由月亮和太陽引力產生的固體潮對地殼加、卸載。方法之二是利用在地震學中常用的地震能量(在實驗中用聲發射能量)來定義加卸載響應比:

式中:E為地震時輻射的地震波能量,N+和N-分別為加載和卸載過程中發生的地震次數,m取0～1的任意值。

用已發生的上百例地震資料進行檢驗,證明加卸載響應比方法效果良好,80%以上的震例表明:在地震孕育初期,Y在1附近波動;在強震來臨的時候,Y會升高;然后,Y升高到最大值(明顯大于1);在主震發生的前夕,Y又會迅速下降[30]。

與此同時,尹祥礎也指出,雖然“取得了一些新認識、新成果,但距離地震預測的宏偉目標,還有很長的路要走…”。LURR理論的主體思路和本文前幾節的思路是一致的,結合上述有關材料損傷演化和強度以及脈搏波反演的討論,至少可以從以下3方面對LURR理論作進一步的發展:

(1)LURR理論是建立在線彈性本構關系基礎上的,表觀加載曲線的非線性都歸為損傷演化。其實,更一般情況下,表觀應力應變曲線的非線性乃是材料本構非線性與損傷演化響應之耦合。如果要把LURR理論推廣到非線性本構關系,就要區分本構非線性響應和損傷演化所致的非線性響應。

(2)不要光停留在地震波本身的分析上,要通過地震波信息去反演更為本質的本構關系及相耦合的損傷演化,這樣至少要在相鄰3點處對地球板塊“把脈”(三點律)。

(3)從相鄰3點處實測的地震波反演出含損傷演化的本構曲線后,還要進一步區分相耦合的損傷演化信息及規律。

期望上述幾點考慮有助于提升LURR理論,通過地震波來給地球把脈,以實現更有效的地震預報。

5 結 論

(1)材料動態強度是材料在動載荷作用下抵抗流變、損傷演化直至破斷的能力,以計及應變率效應為特征。

(2)破斷主要是損傷演化過程的后果,而損傷與流變是耦合在一起發展的;從含損傷的表觀本構關系如何把兩者解耦是研究的難點和關鍵。

(3)基于損傷演化微觀/細觀研究的熱激活機理及實驗,發現宏觀損傷的動態演化也遵循熱激活演化律(TADE模型),并為實驗研究所證實。

(4)為中醫脈診客觀化研究提出了一個新途徑,即通過相鄰3處的脈搏波信息來定量反演脈搏波系統的正常及病態本構關系(與中醫的體質相當),由此反映生命體偏離正常狀態的程度—病情(廣義的損傷演化和強度)。

(5)地震預報是通過地震波信息來診斷地球偏離正常運行的狀態(震情)—另一種廣義的損傷演化和強度。與LURR理論相結合,建議至少在相鄰3點處對地球板塊“把脈”,從同時測得的地震波來反演含損傷演化的載荷-響應曲線,進而從中將兩者加以區分。這將有利于提高地震預報的有效性。

[1] 中國大百科全書總委員會.中國大百科全書:力學[M].北京:中國大百科全書出版社,1985:397.

[2] Timoshenko S P.Strength of materials[M].New York:Van Nostrand Company,1930.

[3] Lindholm U S.Review of dynamic testing techniques and material behavior[C]∥Proceedings of Conference on Mechanical Properties of Materials at High Rates of Strain.London:Institute of Physics,1974:3-21.

[4] Kachanov L M.Time of the rupture process under creep conditions[J].Izv AN SSSR Otd Tekhn Nauk,1958,8: 26-31.

[5] 王禮立.絕熱剪切:材料在沖擊載荷下的本構失穩[M].王禮立,余同希,李永池.沖擊動力學進展.合肥:中國科學技術大學出版社,1992:3-33.

[6] 王禮立,蔣昭鑣,陳江瑛.材料微損傷在高速變形過程中的演化及其對率型本構關系的影響[J].寧波大學學報(理工版),1996,9(3):47-55. Wang Lili,Jiang Zhaobiao,Chen Jiangying.Micro-damage evolution in high velocity deformation and its influence on rate-dependent constitutive relation of materials[J].Journal of Ningbo University(Natural Science&Engineering Edition),1996,9(3):47-55.

[7] Wang Lili,Jiang Zhaobiao,Chen Jiangying.Studies on rheological relation of materials by taking account of ratedependent evolution of internal defects at high strain rates[M]∥Wang Ren.Rheology of Bodies with Defects.Dordrecht:Kluwer Academic Publishers,1999:167-178.

[8] Huang Dejin,Shi Shaoqiu,Wang Lili.Studies on rate-dependent evolution of damage and its effects on dynamic constitutive response by using a random fuse network model[M]∥Chiba A,Tanimura S,Hokamoto K.Impact Engineering and Application.Tokyo:Elsevier Science Ltd.,2001:743-748.

[9] Wang Lili,Zhou Fenghua,Sun Zijian,et al.Studies on rate-dependent macro-damage evolution of materials at high strain rates[J].International Journal of Damage Mechanics,2010,19(7):805-820.

[10] Wang Lili,Hu Shisheng,Yang Liming,et al.Development of experimental methods for impact testing by combining Hopkinson pressure bar with other techniques[J].Acta Mechanica Solida Sinica,2014,27(4):331-344.

[11] 王禮立,胡時勝,楊黎明,等.材料動力學[M].合肥:中國科學技術大學出版社,2016.

[12] Curran D R,Shockey D A,Seaman L.Dynamic fracture criteria for a polycarbonate[J].Journal of Applied Physics,1973,44(9):4025-4038.

[13] Seaman L,Curran D R,Shockey D A.Computational models for ductile and brittle fracture[J].Journal of Applied Physics,1976,47(11):4814-4826.

[14] Curran D R,Seaman L,Shockey D A.Dynamic failure of solid[J].Physics Report,1987,147(5/6):253-388.

[15] Curran D R,Seaman L.Simplified models of fracture and fragmentation[C]∥High-pressure Shock Compression of Solids II:Dynamic Fracture and Fragmentation.New York:Springer-Verlag,1996:340-365.

[16] Zhurkov S N,Sanfirova T P.The temperature and time dependence of the strength of pure metals[J].Dokl Akad-Nauk SSSR,1955,101:237.

[17] 周風華,王禮立,胡時勝.有機玻璃在高應變率下的損傷型非線性粘彈性本構關系及破壞準則[J].爆炸與沖擊, 1992,12(4):333-342. Zhou Fenghua,Wang Lili,Hu Shisheng.A damage-modified nonlinear visco-elastic constitutive relation and failure criterion of PMMA at high strain-rates[J].Explosion and Shock Waves,1992,12(4):333-342.

[18] Wang Lili,Xu Mingqiao,Shi Shaoqiu.Application of BP neural network to SHPB technique for the investigation of impact response of polymers[C]∥Proceedings of 2003 SEM Annual Conference on Experimental and Applied Mechanics.Charlotte,North Carolina,USA,2003:167-174.

[19] 王禮立,賴華偉,孫紫建,等.高聚物計及損傷演化的動態變形和斷裂[J].寧波大學學報(理工版),2003,16(6): 372-380. Wang Lili,Lai Huawei,Sun Zijian,et al.Dynamic deformation and fracture of polymers taking account of damage evolution[J].Journal of Ningbo University(Natural Science&Engineering Edition),2003,16(6):372-380.

[20] Xu Mingqiao,Wang Lili.A new method for studying the dynamic response and damage evolution of polymers at high strain rates[J].Mechanics of Materials,2006,38(1/2):68-75.

[21] Wang L,Xu M,Zhu J,et al.A method of combined SHPB technique and BP neural network to study impact response of materials[J].Strain,2006,42(3):149-158.

[22] Sun Zijian,Wang Lili.Studies on impact constitutive behavior and dmage evolutionfor PP/PA polymer blends at large deformation[J].Journal de Physique IV,2006,134:117-124.

[23] Wang Lili,Wang Hui.Mechanics modeling and inverse analyses of pulse wave system from the view-point of traditional Chinese medicine[C]∥Proceedings of ASME 2016 35th International Conference on Ocean,Offshore and Arctic Engineering.Busan,South Korea,2016:OMAE2016-55106.

[24] 王禮立,王暉.脈搏波系統的力學模型及反演兼對若干中醫學問題的討論[J].力學學報,2016,48(6):1416-1424. Wang Lili,Wang Hui.Mechanical modeling and inverse analyses of pulsewaves system with discussions on some concepts in the traditional Chinese medicine[J].Chinese Journal of Theoretical and Applied Mechanics,2016,48(6):1416-1424.

[25] Wang Lili.Foundations of stress waves[M].Amsterdam:Elsevier,2007.

[26] 王琦.中醫體質學[M].北京:人民衛生出版社,2009.

[27] 王琦.九種體質使用手冊[M].長春:北方婦女兒童出版社,2010.

[28] 王暉.體質的中醫保健[M].寧波:寧波出版社,2009.

[29] 王暉.全國名老中醫王暉學術經驗擷英[M].北京:中國中醫藥出版社,2014.

[30] 尹祥礎.地震預測新途徑的探索[J].中國地震,1987,3(1):1-7. Yin Xiangchu.The new approach of earchquake prediction[J].Earthquake Research in China,1987,3(1):1-7.

[31] 尹祥礎.加卸載響應比理論及其應用[M].北京:科學出版社,2016.

Talk about dynamic strength and damage evolution

Wang Lili1,Hu Shisheng2,Yang Liming1,Dong Xinlong1,Wang Hui3
(1.MOE Key Laboratory of Impact and Safety Engineering,Ningbo University,Ningbo315211,Zhejiang,China; 2.CAS Key Laboratory for Mechanical Behavior and Design of Materials,University of Science and Technology of China,Hefei230027,Anhui,China; 3.Wang Hui Workroom,Ningbo Hospital of Traditional Chinese Medicine,Ningbo315000,Zhejiang,China)

The strength of a material is traditionally understood as its ability to resist flow/deformation and breakage(brittle fracture or ductile rupture)under applied load.The breakage,though looking like an abrupt occurrence,actually results from a strain-rate/time dependent process of damage evolution.The difficulty in studying dynamic damage evolution lies in that the damage evolution and the flow/deformation process are coupled and influence each other.It was found that the dynamic evolution of macro-damage could be successfully described by the thermo-activated damage evolution model.Based on this model,the damage evolution and the flow/deformation can be decoupled from the experimentally measured apparent stress-strain curve with damage evolution,and the related material parameters can be determined.Such an approach is then generalized to the objective study of the TCM pulse.The normal and pathological constitutive relations of the pulse wave system can be inversely determined by pulse wave signals,and then the degree of deviation from the normal condition (illness state)of patients can be diagnosed,and their illness state can be regarded as a kind of generalized damage.The same approach is further generalized to the study of earthquake prediction through the‘pulse-taking for the earth’.Combined with the load-unload response ratio theory,by measuring the seismic wave signals on three adjacent positions,the nonlinear constitutive load-response curve with damage evolution can be inversely determined for the tectonic plates of the earth concerned,then the degree of the damage can be finally discriminated,which is the key information for an improved earthquake prediction.

dynamic strength;damage evolution;dynamics of materials;pulse wave;earthquake prediction

O347國標學科代碼:13015

:A

10.11883/1001-1455(2017)02-0169-11

(責任編輯 張凌云)

2016-12-23;

:2017-01-19

王禮立(1934— ),男,教授,博士生導師,wanglili@nbu.edu.cn。本文根據作者在第十一屆全國爆炸力學學術會議的大會特約報告整理而成。