基于J積分的含裂紋加筋板在單軸拉伸載荷下的極限強(qiáng)度計(jì)算

夏世林，磨季云

(武漢科技大學(xué) 理學(xué)院，湖北 武漢 430065)

0 引 言

加筋板結(jié)構(gòu)是船舶、航空航天領(lǐng)域最為常用的結(jié)構(gòu)類型之一。對(duì)于此類結(jié)構(gòu)，當(dāng)結(jié)構(gòu)受力達(dá)到材料屈服強(qiáng)度后，其仍有一定的承載能力，直至載荷水平達(dá)到其極限強(qiáng)度。在加筋板結(jié)構(gòu)設(shè)計(jì)中，極限狀態(tài)（ULS）設(shè)計(jì)已經(jīng)被公認(rèn)為優(yōu)于傳統(tǒng)的許用應(yīng)力方法，因?yàn)榍罢吣芨鎸?shí)反映結(jié)構(gòu)的安全極限[1]。加筋板結(jié)構(gòu)在服役年限中不可避免遭受不同種類的損傷，裂紋損傷就是影響結(jié)構(gòu)失效行為的重要因素之一。在加筋板結(jié)構(gòu)中，裂紋經(jīng)常在焊接和應(yīng)力集中區(qū)域萌生，例如在加強(qiáng)筋和底板的交接處。在加筋板壽命周期中，這些裂紋可能在不同載荷狀態(tài)下擴(kuò)展，結(jié)果形成不同的尺寸和位置。

按載荷形式的不同，含裂紋加筋板的剩余強(qiáng)度問(wèn)題可分為壓縮載荷下加筋板的屈曲強(qiáng)度、拉伸載荷下加筋板的極限強(qiáng)度和剪切載荷下加筋板的剩余強(qiáng)度。對(duì)于此類問(wèn)題，傳統(tǒng)的有限元方法是將裂紋建立成類似于狹縫的幾何模型，并在裂紋尖端進(jìn)行網(wǎng)格加密，從而將斷裂力學(xué)問(wèn)題簡(jiǎn)化為一般靜力學(xué)問(wèn)題。Sadek等[2]用數(shù)值模擬的方法分析了含裂紋加筋板在壓縮應(yīng)力下的屈曲特性，分析中考慮了裂紋的演化發(fā)展；Cao等[3]根據(jù)歐洲規(guī)范EN1990，EN1991和EN1993，對(duì)環(huán)向和經(jīng)向制造裂紋的圓形鋼筒倉(cāng)的屈曲性能進(jìn)行了評(píng)估；Duc等[4]通過(guò)數(shù)值模擬計(jì)算，建立了基于裂紋有限元分析的靜穩(wěn)定板模型。Brighenti等[5]研究了裂紋損傷對(duì)拉伸、壓縮、剪切載荷下的矩形板屈曲強(qiáng)度的影響，其中裂紋長(zhǎng)度、裂紋方向和材料泊松比是變化的，結(jié)果表明裂紋方向垂直拉伸和壓縮載荷方向時(shí)是最危險(xiǎn)的工況。Margaritis Y等[9]研究了含裂紋加筋板的極限和失效響應(yīng)，詳細(xì)檢查了裂紋閉合對(duì)結(jié)構(gòu)響應(yīng)的作用。Paik等[10]對(duì)完整和含裂紋加筋板的極限剪切強(qiáng)度進(jìn)行了數(shù)值研究，并基于完整加筋板方程發(fā)展了預(yù)報(bào)含裂紋加筋板的極限剪切強(qiáng)度方程。Saad-Eldeen等[11]對(duì)含中心橢圓開口和封閉裂紋的薄壁鋼板剩余強(qiáng)度進(jìn)行了試驗(yàn)研究，試驗(yàn)中的試件是從服役期中的實(shí)際結(jié)構(gòu)構(gòu)件中取得的，而不是新建試件。

與傳統(tǒng)的一般靜力學(xué)分析方法相比，J積分方法具有以下優(yōu)勢(shì)：一是不需要建立裂紋模型，而是從斷裂力學(xué)機(jī)理出發(fā)，直接定義裂尖和裂紋方向，這更加符合裂紋模型的實(shí)際情況；二是J積分方法避開了直接求解裂端塑性應(yīng)力場(chǎng)的困難。Cherepanov[6]和Rice[7]各自獨(dú)立提出了一種與積分路徑無(wú)關(guān)的J積分，用于綜合度量裂端應(yīng)力應(yīng)變場(chǎng)強(qiáng)度。它對(duì)彈塑性斷裂力學(xué)發(fā)展起到了重要的作用。

通過(guò)以上文獻(xiàn)可以看出，以前的研究主要關(guān)注含裂紋結(jié)構(gòu)在典型拉伸載荷下的失效行為或塑性屈曲強(qiáng)度評(píng)估。含裂紋結(jié)構(gòu)在軸向拉伸下的剩余極限強(qiáng)度在近年來(lái)得到了廣泛關(guān)注，但是該領(lǐng)域的研究仍然是有限、不充分的。因此，需要進(jìn)行更多的研究來(lái)更深理解含裂紋缺陷結(jié)構(gòu)構(gòu)件的失效行為。本文的主旨在于發(fā)現(xiàn)含裂紋損傷的加筋板在單調(diào)軸向拉伸下的極限強(qiáng)度行為。假定板橫向裂紋位于加筋板內(nèi)截面，或只位于板單元上，或位于板單元和加強(qiáng)筋上。需要指出的是，本文不考慮裂紋擴(kuò)展和斷裂相關(guān)的臨界裂紋長(zhǎng)度。

加筋板材料為Q345B高強(qiáng)度鋼，為了達(dá)到目的，進(jìn)行了含一系列不同尺寸、位置裂紋加筋板在單軸拉伸載荷下的非線性彈塑性大變形分析。

1 力學(xué)分析方法

1.1 J積分及其有限元列式

J積分理論避開了直接計(jì)算裂紋尖端附近的彈塑性應(yīng)力、應(yīng)變場(chǎng)，而用遠(yuǎn)場(chǎng)J積分作為表示裂紋尖端應(yīng)力應(yīng)變集中特征的平均參量。對(duì)于二維問(wèn)題，J積分可定義為[3]：

其中： Γ 為由裂紋下表面某點(diǎn)到裂紋上表面某點(diǎn)的簡(jiǎn)單的積分路徑；W為彈塑性應(yīng)變能密度；Ti為作用于積分回路單位周長(zhǎng)上的主應(yīng)力；ui為積分回路邊界上的位移；ds為積分回路線的弧長(zhǎng)。

在線彈性斷裂力學(xué)中，J積分與應(yīng)力強(qiáng)度因子KI有如下關(guān)系式：

式中：對(duì)于平面應(yīng)力，彈性模量E′=E，對(duì)于平面應(yīng)變E′=E/1-μ2。

可以解析證明，J積分與積分路線 Γ 的選取無(wú)關(guān)，即J積分滿足路徑守恒。因此，可選取應(yīng)力應(yīng)變場(chǎng)較易求解的線路得到J積分值，而此值與積分路徑非常靠近裂端的結(jié)果是相同的。J積分有著明確的物理意義，它代表著裂紋沿x方向擴(kuò)展單位長(zhǎng)度時(shí)的能量釋放率。

由于J積分與選擇的路徑無(wú)關(guān)，并能描述裂紋尖端的應(yīng)變狀態(tài)。這一特點(diǎn)使J積分計(jì)算可以在離開高應(yīng)力、高應(yīng)變梯度的裂紋尖端區(qū)外的任意回線進(jìn)行。對(duì)于二維體的J積分列式可以寫為[11]：

其中：

將平面上定義的二維J積分沿著裂紋縫線逐點(diǎn)積分可得裂紋線的三維J積分為：

式中：L為裂紋縫線長(zhǎng)度；dA為包繞裂紋縫線的圓柱曲面上的面積微元， dA=dsdΓ；n為微元dA的外法線方向； λ (s)為裂紋縫線s點(diǎn)處的長(zhǎng)度。

J積分的計(jì)算必須采用數(shù)值積分，如果回線經(jīng)過(guò)單元積分點(diǎn)，可以采用有限元中的高斯積分策略。從遠(yuǎn)離裂紋尖端的應(yīng)力場(chǎng)、應(yīng)變場(chǎng)計(jì)算J積分，只能在無(wú)體力、初始應(yīng)變和裂紋表面拉力情況下能用。當(dāng)有體力或初始應(yīng)變時(shí)，J積分要加入封閉表面的積分。此時(shí)路徑無(wú)關(guān)性將喪失，必須從完整的應(yīng)力和應(yīng)變結(jié)果中計(jì)算，包括裂紋尖端的高應(yīng)力區(qū)。

J積分的主要優(yōu)點(diǎn)是不需要用很細(xì)的網(wǎng)格對(duì)裂紋尖端區(qū)進(jìn)行分析。路徑無(wú)關(guān)性的非線彈性材料也存在，真實(shí)的彈塑性材料與非線彈性材料類似。J積分也可以用于描述彈塑性材料的裂紋尖端的應(yīng)力場(chǎng)和應(yīng)變場(chǎng)，材料用增量塑性理論描述并不允許有卸載，以免W的加載歷程相關(guān)性，J值的增量可以從每個(gè)載荷步的應(yīng)力和應(yīng)變?cè)隽坑?jì)算得到。

1.2 精度驗(yàn)證算例

以緊湊拉伸(CT)式樣為例進(jìn)行說(shuō)明，材料選用304不銹鋼，厚度為5 mm，對(duì)其進(jìn)行單項(xiàng)拉伸，CT式樣尺寸如圖1所示。為了效驗(yàn)計(jì)算結(jié)果的準(zhǔn)確性，將各方法的計(jì)算結(jié)果與(2)式用應(yīng)力強(qiáng)度因子法計(jì)算的理論值做比較。根據(jù)上述計(jì)算J積分的方法，對(duì)該實(shí)例進(jìn)行計(jì)算，結(jié)果如圖2所示。

圖1 CT試樣尺寸Fig.1 Dimension of the CT specimen

圖2 J積分計(jì)算結(jié)果圖Fig.2 The scheme of J-integral calculation result

由圖2可知，各種J積分計(jì)算方法的結(jié)果具有相同的增長(zhǎng)趨勢(shì)且相差不大，通過(guò)與理論值的對(duì)比，驗(yàn)證了各種常用J積分計(jì)算方法的準(zhǔn)確性與可行性。

1.3 有限元模型

1.3.1 幾何和材料

本文的分析對(duì)象為含2根加強(qiáng)筋（截面對(duì)稱）和1個(gè)板單元的加筋板，其長(zhǎng)細(xì)比范圍為4～6，如圖3所示。在圖3中，板的尺寸用a×b×t表述，加強(qiáng)筋的尺寸用hw×tw表述。底板和加強(qiáng)筋的詳細(xì)尺寸見(jiàn)表1。

圖3 加筋板的模型尺寸Fig.3 Model dimension of the stiffened-panel

表1 加筋板的幾何特征Tab.1 Geometrical features of the stiffened-panel

其中長(zhǎng)細(xì)比β是控制加筋板屈曲/失效行為的重要參數(shù)，定義如下[13]：

其中：b為板寬；t為板厚； σs為材料屈服強(qiáng)度；E為楊氏模量。

本文使用的材料為Q345B鋼。一般來(lái)說(shuō)，鋼結(jié)構(gòu)具有應(yīng)變硬化正切模量Et，其典型取值范圍為楊氏模量的5%～15%。顯然應(yīng)變硬化對(duì)加筋板的的非線性行為有一定影響。影響的程度是包括板長(zhǎng)細(xì)比在內(nèi)的眾多因素的函數(shù)。在本文中，板的材料行為建模成應(yīng)變硬化率為E/65的雙線性彈塑性行為（見(jiàn)圖4），該應(yīng)變硬化率的值是Khedmati[14]通過(guò)大量彈塑性大變形分析得到的。詳細(xì)材料屬性如表2所示，材料的應(yīng)力-應(yīng)變行為如圖4所示。

1.3.2 裂紋損傷研究工況

為了考慮裂紋的不同位置，本文應(yīng)用了3種不同的模型，如圖5所示。這些模型按裂紋位置分類，如裂紋只在底板上，裂紋在底板和加強(qiáng)筋上，裂紋只在加強(qiáng)筋上。

裂紋在加筋板中的深度取決于多種因素，其中裂紋深度等于板厚時(shí)最為危險(xiǎn)。在工程中，在相對(duì)較薄的加筋板中，裂紋一般被認(rèn)為貫穿板厚。因此，在本文研究中，底板和加強(qiáng)筋上的裂紋都假設(shè)貫穿厚度。在所有的3種模型中，不論裂紋在底板上還是加強(qiáng)筋上，都沿板的橫向。為了考慮不同裂紋長(zhǎng)度，假設(shè)加筋板上的裂紋長(zhǎng)度c是對(duì)應(yīng)板寬的某一特定比例，其比值c/b的變化范圍為0.05～0.15，而圖5（b）和圖5（c）中在加強(qiáng)筋上的裂紋長(zhǎng)度等于加強(qiáng)筋高度。表3描述了本文考慮的不同幾何參數(shù)。

表2 Q345B鋼材料屬性匯總Tab.2 A summary of Q345B steel material properties

圖4 材料的理想彈塑性應(yīng)力-應(yīng)變行為Fig.4 Ideal elastic-plastic stress-strain behavior of the material

圖5 3種裂紋模型Fig.5 3 types of crack models

表3 FE模型幾何參數(shù)值匯總Tab.3 A summary of FE models’ geometrical parameters

1.3.3 載荷和邊界條件

在進(jìn)行FE分析之前，在FE模型上施加合適的邊界條件盡可能合理模擬結(jié)構(gòu)實(shí)際響應(yīng)是十分重要的。一個(gè)常規(guī)的方法是在極限狀態(tài)評(píng)估中在橫向邊緣施加簡(jiǎn)支或夾持邊界條件。本文主要關(guān)注的是加筋板在單軸拉伸載荷下的極限強(qiáng)度。通常認(rèn)為加筋板橫剖面會(huì)受到均勻的拉伸外載荷力，可采用位移控制的方式來(lái)實(shí)現(xiàn)加筋板均勻軸向載荷的逐步加載。圖6描述了含裂紋加筋板模型的載荷和邊界條件，其中t為時(shí)間。

圖6 含裂紋加筋板載荷和邊界條件Fig.6 Loads and BCs of the stiffened-panel with cracks

1.3.4 單元和網(wǎng)格尺寸

結(jié)構(gòu)極限強(qiáng)度計(jì)算是非線性問(wèn)題。在本文的研究中，含裂紋加筋板的失效行為和極限強(qiáng)度評(píng)估使用有限元軟件Abaqus，評(píng)估中同時(shí)考慮幾何和材料非線性。在Abaqus包含可使用的單元庫(kù)中，選擇8節(jié)點(diǎn)S8R單元來(lái)離散加筋板模型。S8R適合屬于加筋板結(jié)構(gòu)范疇的適度薄壁殼結(jié)構(gòu)。其適合線性、大應(yīng)變和大轉(zhuǎn)動(dòng)非線性應(yīng)用。

在FE分析中，計(jì)算精度和速度在很大程度上取決于單元大小或網(wǎng)格密度。網(wǎng)格尺寸必須足夠精細(xì)來(lái)捕捉含裂紋加筋板的屈曲/塑性失效行為。另一方面，更為精確的網(wǎng)格引起更大的計(jì)算代價(jià)。從這一點(diǎn)來(lái)說(shuō)，必須進(jìn)行精度和計(jì)算時(shí)間的平衡。對(duì)于含裂紋加筋板，精確網(wǎng)格必須分配在局部裂紋區(qū)域來(lái)模擬更早的屈服和由于應(yīng)力集中引起的高應(yīng)力梯度。

進(jìn)行不同網(wǎng)格尺寸的收斂性研究決定相對(duì)合適的非線性有限元分析網(wǎng)格模型。最終，沿板寬方向b和板長(zhǎng)方向a的殼單元數(shù)目指定為20。對(duì)于加強(qiáng)筋腹板在高度方向hw單元數(shù)量為1。同時(shí)在裂紋附近進(jìn)行網(wǎng)格細(xì)化。裂紋區(qū)域網(wǎng)格細(xì)化的含裂紋加筋板網(wǎng)格劃分實(shí)例如圖7所示。

圖7 含裂紋加筋板網(wǎng)格實(shí)例Fig.7 Mesh example of the stiffened-panel with cracks

2 分析和討論

使用Abaqus有限元軟件進(jìn)行一系列彈塑性大變形分析研究含裂紋加筋板在單軸拉伸載荷下的失效行為。在分析中，使用牛頓-拉普生方法來(lái)解決非線性控制方程，考慮上文描述的3種裂紋模型。而且改變裂紋長(zhǎng)度和板厚來(lái)研究其對(duì)含裂紋加筋板極限強(qiáng)度特性的影響。

2.1 長(zhǎng)細(xì)比的影響

考慮表3中所列的不同板厚來(lái)研究其對(duì)含裂紋加筋板極限強(qiáng)度特性的影響。選擇裂紋長(zhǎng)度c/b=0.1和裂紋縱向位置固定于加筋板的中橫剖面。將加載試驗(yàn)中采集的軸向荷載和軸向拉伸變形進(jìn)行換算，即將軸向荷載除以實(shí)測(cè)的試件截面積，得到截面的平均應(yīng)力值；將軸向拉伸變形除以試件的實(shí)測(cè)初始長(zhǎng)度，得到試件的平均應(yīng)變值[15]，圖8給出了不同長(zhǎng)細(xì)比的模型（a）應(yīng)力-應(yīng)變曲線（裂紋長(zhǎng)度c/b=0.1）。圖中也包含了相關(guān)無(wú)裂紋損傷完整加筋板作為參考模型。應(yīng)力-應(yīng)變曲線的峰值即為加筋板的極限拉伸強(qiáng)度，當(dāng)載荷水平達(dá)到該值時(shí)，含裂紋加筋板的應(yīng)力云圖（以β=5為例）如圖9所示。可以看出，當(dāng)載荷水平達(dá)到含裂紋加筋板極限強(qiáng)度時(shí)，加筋板中橫剖面幾乎所有單元的等效米塞斯應(yīng)力都達(dá)到材料的抗拉強(qiáng)度，加筋板失去抗拉能力。

為了測(cè)量裂紋損傷對(duì)加筋板極限拉伸強(qiáng)度的影響，將結(jié)果與相關(guān)無(wú)裂紋損傷完整加筋板進(jìn)行對(duì)比是必要的。從這一點(diǎn)出發(fā)，引入極限強(qiáng)度縮減系數(shù)，如下式：

圖8 不同長(zhǎng)細(xì)比的模型（a）應(yīng)力-應(yīng)變曲線（裂紋長(zhǎng)度c/b=0.1）Fig.8 The stress-strain curve of Model (a) with different slenderness (crack length c/b=0.1)

圖9 含裂紋加筋板的應(yīng)力云圖（β=5）Fig.9 The stress distribution of the stiffened-panel with cracks（β=5）

其中： σu為含裂紋加筋板極限強(qiáng)度； σu0為相關(guān)完整加筋板極限強(qiáng)度。

3種長(zhǎng)細(xì)比下極限強(qiáng)度縮減系數(shù)分別為0.914，0.902，0.895，含裂紋加筋板極限拉伸強(qiáng)度隨加筋板長(zhǎng)細(xì)比的增大略有減小，說(shuō)明當(dāng)加筋板厚度越大，裂紋損傷對(duì)加筋板極限強(qiáng)度的影響越明顯。總體來(lái)看，加筋板長(zhǎng)細(xì)比對(duì)極限強(qiáng)度的影響不大，這是因?yàn)樵诒疚牡乃憷?，?dāng)加筋板板厚增大時(shí)，裂紋深度也隨之增大，表明含貫穿型裂紋的加筋板在單軸拉伸載荷下的剩余強(qiáng)度對(duì)加筋板長(zhǎng)細(xì)比不敏感。

2.2 裂紋長(zhǎng)度的影響

為了研究裂紋長(zhǎng)度對(duì)加筋板失效行為的影響，選擇板厚為10 mm，裂紋位于加筋板中橫剖面的工況，改變底板上裂紋的長(zhǎng)度。加強(qiáng)筋上沒(méi)有裂紋。圖10演示了裂紋長(zhǎng)度對(duì)應(yīng)力-應(yīng)變關(guān)系的影響?？芍鸭y長(zhǎng)度的增加影響加筋板結(jié)構(gòu)響應(yīng)，極限強(qiáng)度和剩余強(qiáng)度都隨裂紋長(zhǎng)度的增大而減小。而且，從圖中還可以看出極限應(yīng)變（失效點(diǎn)應(yīng)變）隨裂紋長(zhǎng)度的增大而減?。?種工況下極限應(yīng)變分別為0.158，0.152，0.142）。3種裂紋長(zhǎng)度下極限強(qiáng)度縮減系數(shù)分別為0.939，0.902，0.793，含裂紋加筋板極限拉伸強(qiáng)度隨裂紋長(zhǎng)度的增大而減小，且減小的幅度逐漸增大，說(shuō)明含裂紋加筋板在單軸拉伸載荷下的極限強(qiáng)度對(duì)裂紋長(zhǎng)度較為敏感。

圖10 裂紋長(zhǎng)度對(duì)應(yīng)力-應(yīng)變關(guān)系的影響Fig.10 The influence of crack length on the stress-strain behavior

2.3 裂紋位置的影響

為了檢查裂紋位置對(duì)含裂紋加筋板極限強(qiáng)度行為的影響，設(shè)計(jì)了裂紋只出現(xiàn)在底板上、裂紋只出現(xiàn)在加強(qiáng)筋上、裂紋同時(shí)出現(xiàn)在底板和加強(qiáng)筋上3組仿真試驗(yàn)。圖11和圖12分別為裂紋在加強(qiáng)筋上和裂紋同時(shí)在底板和加強(qiáng)筋上時(shí)的應(yīng)力云圖，圖13為3種工況下含裂紋加筋板應(yīng)力-應(yīng)變曲線。3種工況下極限強(qiáng)度縮減系數(shù)分別為0.939（裂紋只出現(xiàn)在底板上）、0.975（裂紋只出現(xiàn)在加強(qiáng)筋上）、0.924（裂紋同時(shí)出現(xiàn)在底板和加強(qiáng)筋上），說(shuō)明加強(qiáng)筋上的裂紋對(duì)含裂紋加筋板極限強(qiáng)度的影響小于底板上的裂紋，而裂紋同時(shí)出現(xiàn)在底板和加強(qiáng)筋上時(shí)對(duì)含裂紋加筋板極限拉伸強(qiáng)度的影響最大。這是因?yàn)榱鸭y只位于底板和加強(qiáng)筋上時(shí)，裂紋區(qū)域在單軸拉伸載荷下的變形都受到加強(qiáng)筋或底板的限制，而裂紋同時(shí)位于底板和加強(qiáng)筋上時(shí)，裂紋在底板和加強(qiáng)筋相交的位置能夠自由變形。

圖11 裂紋在加強(qiáng)筋上時(shí)的應(yīng)力云圖Fig.11 The stress distribution when the crack locates on the stiffener

圖12 裂紋同時(shí)在底板和加強(qiáng)筋上時(shí)的應(yīng)力云圖Fig.12 The stress distribution when the crack locates on both the stiffener and the plate

圖13 三種工況下含裂紋加筋板應(yīng)力-應(yīng)變曲線Fig.13 The stress-strain curve of the stiffened-panel with cracks under 3 different cases

3 結(jié) 語(yǔ)

本文建立了一系列不同長(zhǎng)細(xì)比、不同裂紋長(zhǎng)度、不同裂紋位置的含裂紋加筋板有限元模型，基于J積分理論對(duì)含裂紋加筋板在單軸拉伸載荷下的極限強(qiáng)度進(jìn)行了計(jì)算分析，得到以下結(jié)論：

1）長(zhǎng)細(xì)比對(duì)含貫穿型裂紋的加筋板在單軸拉伸載荷下的剩余強(qiáng)度影響不大；

2）含裂紋加筋板極限拉伸強(qiáng)度隨裂紋長(zhǎng)度的增大而減小，且減小的幅度逐漸增大，極限拉伸強(qiáng)度對(duì)裂紋長(zhǎng)度較為敏感；

3）加強(qiáng)筋上的裂紋對(duì)含裂紋加筋板極限強(qiáng)度的影響小于底板上的裂紋，而裂紋同時(shí)出現(xiàn)在底板和加強(qiáng)筋上時(shí)對(duì)含裂紋加筋板極限拉伸強(qiáng)度的影響最大。

本文在進(jìn)行含裂紋加筋板極限強(qiáng)度計(jì)算時(shí)沒(méi)有考慮裂紋擴(kuò)展，這或多或少與實(shí)際情況存在一定的差異。裂紋擴(kuò)展對(duì)含裂紋加筋板極限強(qiáng)度的影響，是下一步研究的方向。