• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    改進的降階牛頓迭代數(shù)值子結(jié)構(gòu)方法

    2022-08-01 00:56:58孫寶印孫天舒歐進萍
    工程力學 2022年8期
    關(guān)鍵詞:子結(jié)構(gòu)屈服塑性

    孫寶印,申 偉,孫天舒,歐進萍,4

    (1. 河海大學土木與交通學院,江蘇,南京 210098;2. 廣西大學土木建筑工程學院,廣西,南寧 530004;3. 大連理工大學土木工程學院,遼寧,大連 116024;4. 哈爾濱工業(yè)大學土木工程學院,黑龍江,哈爾濱 150090)

    計算效率和模擬精度是高層建筑等重大工程地震彈塑性全過程分析的關(guān)鍵技術(shù)瓶頸。一方面,此類結(jié)構(gòu)的規(guī)模龐大,為了精確刻畫其非線性行為,有必要針對整體結(jié)構(gòu)建立精細化的有限元模型,自由度動輒數(shù)十上百萬,結(jié)構(gòu)剛度矩陣的集成和分解耗時嚴重。更為關(guān)鍵的是,結(jié)構(gòu)一旦進入非線性狀態(tài),傳統(tǒng)的變剛度法需要反復(fù)迭代更新結(jié)構(gòu)剛度,計算效率更是成為棘手難題。另一方面,高層結(jié)構(gòu)地震倒塌破壞主要是由局部關(guān)鍵構(gòu)件嚴重損傷導(dǎo)致,而這部分關(guān)鍵構(gòu)件的非線性變形特征具有顯著的多樣性和復(fù)雜性,對模型精細程度提出了較高的要求。實際工程中為了滿足計算效率要求,無法針對整體結(jié)構(gòu)建立精細模型,進而難以保證局部關(guān)鍵構(gòu)件的模擬精度。由此可見,計算效率和模擬精度之間相互矛盾、制約?,F(xiàn)有的計算方法較難有效兼顧這類大規(guī)模結(jié)構(gòu)非線性分析的效率和精度,因此,發(fā)展更加高效、精確的非線性分析理論和計算方法十分必要。

    高層結(jié)構(gòu)在地震作用下呈現(xiàn)局部非線性特征,如何充分利用這種局部化特征來提高大規(guī)模結(jié)構(gòu)非線性分析的計算效率,是近幾十年來國內(nèi)外學者關(guān)注的重點。Clough 等[1]利用子結(jié)構(gòu)方法對局部非線性的大規(guī)模結(jié)構(gòu)系統(tǒng)進行動力分析。將整體結(jié)構(gòu)劃分為線彈性和塑性區(qū)域:一方面,通過靜力凝聚消去線彈性區(qū)域內(nèi)部自由度,降低求解矩陣規(guī)模來提高計算效率;另一方面,提高塑性區(qū)模型精細程度來有效兼顧模擬精度。Wong 等[2]提出擬力法,定義能夠描述結(jié)構(gòu)構(gòu)件非線性行為的塑性自由度,建立相應(yīng)的塑性機制,來模擬結(jié)構(gòu)的非線性行為。并采用初始剛度對狀態(tài)方程進行迭代求解,利用顯式求解算法,可實現(xiàn)高效、穩(wěn)定的結(jié)構(gòu)非線性分析。隨后,Li 等[3-5]在擬力法的基礎(chǔ)上作了進一步改進,將材料應(yīng)變分解為線彈性應(yīng)變和非線性應(yīng)變兩部分,并在單元內(nèi)構(gòu)造非線性應(yīng)變的分布場模型,提出了單元非線性變形機制的模擬方法。再結(jié)合虛功原理,建立具有隔離非線性特征的結(jié)構(gòu)控制方程,在該方程迭代分析時僅需對一個小規(guī)模的Schur 補矩陣進行迭代更新,從而能夠顯著提高結(jié)構(gòu)非線性分析的計算效率。此外,陸新征等[6]、Sun 等[7-8]以及陳宇等[9]從多尺度建模角度,對整體結(jié)構(gòu)進行差異化建模分析,在保證模擬精度的前提下,可大幅降低結(jié)構(gòu)模型自由度。

    傳統(tǒng)的子結(jié)構(gòu)方法、擬力法以及多尺度方法均需要預(yù)判結(jié)構(gòu)的塑性區(qū)域,并且模型網(wǎng)格一旦確認,分析過程中較難更改。為了規(guī)避這些缺陷,孫寶印等[10-13]提出數(shù)值子結(jié)構(gòu)方法(如圖1所示),將原始高層結(jié)構(gòu)的非線性分析轉(zhuǎn)化為主結(jié)構(gòu)的等效線彈性分析和屈服構(gòu)件的隔離子結(jié)構(gòu)非線性分析。在整個地震彈塑性時程分析過程中,由線彈性主結(jié)構(gòu)模型計算得到整體結(jié)構(gòu)的動力響應(yīng),根據(jù)結(jié)構(gòu)全局位移實時判斷構(gòu)件的屈服狀態(tài)。將構(gòu)件屈服隔離并建立精細化的子結(jié)構(gòu)模型,由子結(jié)構(gòu)非線性分析計算得到屈服構(gòu)件的真實內(nèi)力,并將其傳回主結(jié)構(gòu)。主結(jié)構(gòu)將轉(zhuǎn)化得到的屈服構(gòu)件等效外荷載作為外力項施加在屈服構(gòu)件端部節(jié)點上,最后,再由主結(jié)構(gòu)進行等效線彈性分析得到整體結(jié)構(gòu)的真實響應(yīng)。經(jīng)研究表明:該方法可以有效兼顧大規(guī)模結(jié)構(gòu)非線性分析的計算效率和模擬精度。

    由于上述數(shù)值子結(jié)構(gòu)方法采用初始彈性剛度進行非線性迭代分析,仍存在一定不足。于是,本文擬對該方法作進一步改進,提出一種降階的牛頓迭代數(shù)值子結(jié)構(gòu)方法。首先,簡介數(shù)值子結(jié)構(gòu)方法基本原理及其不足;其次,提出降階的牛頓迭代算法,并介紹屈服構(gòu)件的隔離子結(jié)構(gòu)非線性分析;最后,通過一平面15 層3 跨鋼結(jié)構(gòu)的地震彈塑性時程分析來驗證所提方法的精度和效率。

    1 數(shù)值子結(jié)構(gòu)方法基本原理及不足

    1.1 數(shù)值子結(jié)構(gòu)方法

    根據(jù)數(shù)值子結(jié)構(gòu)方法基本原理[10-13],結(jié)構(gòu)運動方程可以表示為:

    式中:U、和分別為結(jié)構(gòu)的位移、速度和加速度;M、C0、K0、F和Fc分別為結(jié)構(gòu)的質(zhì)量矩陣、初始阻尼矩陣、初始剛度矩陣、外力和非線性修正力,分別表示為:

    式中:Ne和Ns分別為總單元數(shù)和屈服單元數(shù);B(e)為用于檢索結(jié)構(gòu)中單元所在位置的布爾矩陣;為屈服單元的非線性修正力,可表示為:

    式中:ke和ue分別為屈服單元的初始彈性剛度和節(jié)點位移;re為由隔離子結(jié)構(gòu)計算得到的屈服單元內(nèi)力。

    式(1)等式左邊是一個線彈性系統(tǒng),結(jié)構(gòu)非線性是由等式右邊的非線性修正力Fc來考慮。對于局部非線性結(jié)構(gòu)系統(tǒng),屈服單元數(shù)Ns<<Ne遠小于總單元數(shù)Ne,則非線性修正力Fc的計算量較小。

    1.2 數(shù)值子結(jié)構(gòu)方法的不足

    以文獻[14]中15 層3 跨的平面鋼框架結(jié)構(gòu)為研究對象。采用不動點迭代數(shù)值子結(jié)構(gòu)方法[10-13]對其進行地震彈塑性時程分析時,主結(jié)構(gòu)中每根梁柱構(gòu)件由30 個基于歐拉-伯努利梁理論的彈性單元模擬,分析過程中根據(jù)截面最外側(cè)應(yīng)變是否超過材料屈服應(yīng)變來判斷構(gòu)件的彈塑性狀態(tài)[15],屈服單元由隔離子結(jié)構(gòu)進行非線性分析。在隔離子結(jié)構(gòu)中,屈服單元采用基于歐拉-伯努利梁理論的分布式塑性鉸纖維單元模擬,鋼材采用OpenSees中Steel02 本構(gòu)模型模擬,材料彈性模量E=2.0×106MPa,硬化比b=0.1。在不同工況和收斂誤差下,不動點迭代數(shù)值子結(jié)構(gòu)方法的迭代次數(shù)時程曲線如圖2 所示,迭代信息如表1 所示??梢钥闯觯斒諗空`差為1.0×10-5mm 時,8 度和9 度罕遇地震作用下,數(shù)值子結(jié)構(gòu)方法需要最多24 步和26 步收斂,平均迭代次數(shù)分別為4.6 和6.9;而當收斂誤差減小至1.0×10-8mm 時,最多需要41 步和44 步,平均迭代次數(shù)分別為7.1 和11.1。當荷載水平增大或收斂誤差減小,不動點迭代數(shù)值子結(jié)構(gòu)方法所需的迭代收斂次數(shù)均有所增加,尤其是當收斂條件變得苛刻,收斂速率顯著降低。于是,下文將對該方法作進一步改進,以提高其收斂速率。

    圖2 不動點迭代數(shù)值子結(jié)構(gòu)方法迭代次數(shù)時程曲線Fig. 2 Iteration number time-history curves using the numerical substructure method based on fixed-point iteration

    表1 不動點迭代數(shù)值子結(jié)構(gòu)方法迭代信息統(tǒng)計表Table 1 Table of iteration information using the numerical substructure method based on fixed-point iteration

    2 降階的牛頓迭代算法

    經(jīng)過時間離散后,第n時步結(jié)構(gòu)的等效運動方程滿足:

    不失一般性,以Newmark-β 逐步積分法為例,第n時步速度和位移分別表示為:

    式中, γ 和 β為控制參數(shù)。本文采用平均加速度法,則 γ 和 β取值分別為0.5 和0.25。

    由式(5)可以得到第n時步速度和加速度:

    式 中:c1=γ/(βΔt);c2=1/β(Δt)2;a1=1-γ/β;a2=(1-γ/2β)Δt;a3=-1/βΔt;a4=1-1/2β。

    將式(6)代入式(4),得到結(jié)構(gòu)的等效平衡方程:

    式中,K和Pn分別為結(jié)構(gòu)等效初始剛度和外力:

    由式(7)可以得到結(jié)構(gòu)的位移場:

    定義與屈服單元有關(guān)的自由度為塑性自由度,則塑性自由度位移us,n可以表示為:

    式中,Bs,n為檢索塑性自由度位置的布爾矩陣。

    將式(10)代入式(11),塑性自由度位移殘差φ(us,n)可以表示為:

    采用牛頓迭代算法,已知當前塑性自由度位移us,n和殘差φ(us,n),得到更新后的塑性自由度位移:

    式中,Hs,n為與塑性自由度有關(guān)的迭代矩陣:

    3 隔離子結(jié)構(gòu)非線性分析

    考慮到高層結(jié)構(gòu)在地震作用下發(fā)生倒塌破壞,局部關(guān)鍵構(gòu)件非線性程度較高,而大部分屈服構(gòu)件的非線性程度較低(稱之為一般屈服構(gòu)件)。為了滿足模擬精度需求,一般屈服構(gòu)件和少數(shù)關(guān)鍵構(gòu)件所采用的隔離子結(jié)構(gòu)模型的精細程度不同。為了簡化起見,本節(jié)以梁柱構(gòu)件為例說明。

    3.1 一般屈服梁柱構(gòu)件

    針對一般屈服梁柱構(gòu)件,纖維模型可以較為準確模擬其非線性行為。由于主、子結(jié)構(gòu)中梁柱構(gòu)件均為桿系模型,它們的邊界自由度一致,因此,由隔離子結(jié)構(gòu)可直接計算得到屈服構(gòu)件內(nèi)力和切線剛度,這里不作詳細贅述。

    3.2 關(guān)鍵梁柱構(gòu)件

    以如圖3(a)所示的平面框架結(jié)構(gòu)為例,將主結(jié)構(gòu)中關(guān)鍵柱構(gòu)件隔離并建立精細化的子結(jié)構(gòu)模型(例如:實體單元),如圖3(b)所示。為了消除剛體位移的影響,將主結(jié)構(gòu)中屈服構(gòu)件對應(yīng)的彈性單元節(jié)點位移ue轉(zhuǎn)換得到局部坐標位移ve:

    圖3 關(guān)鍵構(gòu)件隔離子結(jié)構(gòu)精細化模型邊界處理Fig. 3 Boundary treatment of isolated refined substructures for key components

    式中:上標i為精細化子結(jié)構(gòu)第i個邊界;轉(zhuǎn)換矩陣T可以表示為:

    式中, α為桿系單元軸向與全局坐標系X軸之間的夾角。

    如圖3(c)所示,在局部坐標系下,根據(jù)平截面假定,精細化子結(jié)構(gòu)中第i個邊界中第j個節(jié)點的位移滿足以下轉(zhuǎn)換關(guān)系:

    將式(16)代入式(18),并集成得到全局坐標系下精細化子結(jié)構(gòu)的邊界位移轉(zhuǎn)換關(guān)系:

    式中,C為全局坐標系下子結(jié)構(gòu)的邊界轉(zhuǎn)換矩陣:

    式中,m1和m2分別為第1 和第2 個邊界的節(jié)點數(shù)。

    根據(jù)主、子結(jié)構(gòu)邊界反力所作虛功相等,圖3(d)中邊界反力之間滿足:

    將式(20)代入式(22),得到屈服構(gòu)件內(nèi)力re:

    再將式(23)代入式(3),得到屈服構(gòu)件的非線性修正力=keue-CTqb。

    隔離子結(jié)構(gòu)的增量平衡方程可以表示為:

    式中:下標i和b分別為隔離子結(jié)構(gòu)內(nèi)部和邊界自由度;ks和q分別為隔離子結(jié)構(gòu)的剛度和外力;Δ為增量。

    將式(20)代入式(24),并將式(24)第二個等式左乘CT,再通過靜力凝聚,則有:

    4 分析流程

    采用數(shù)值子結(jié)構(gòu)方法進行結(jié)構(gòu)地震彈塑性時程分析時,具體分析流程如下:

    1)初始化。令n=0、Fn=0 、=0 、=0、Bs,n=0 以及=0。根據(jù)式(8)集成等效剛度矩陣K,計算并保存其逆矩陣K-1。根據(jù)總時步NS P,得到荷載增量 dF=F/NS P。

    2)令n=n+1、i=0、=Bs,n-1以 及=,計算當前時步外力Fn=Fn-1+dF,并由式(9)計算得到等效外力Pn。根據(jù)式(10)計算得到主結(jié)構(gòu)位移,更新整體結(jié)構(gòu)屈服構(gòu)件單元狀態(tài)。若屈服單元數(shù)為零,則系統(tǒng)收斂,跳到第7)步;否則,更新布爾矩陣,由式(11)得到塑性自由度位移,進入下一步。

    5)根據(jù)下式判斷塑性自由度系統(tǒng)的收斂性。

    5 算例分析

    為了驗證本文提出的改進的降階牛頓迭代數(shù)值子結(jié)構(gòu)方法(下文簡稱本文方法)的計算效率和模擬精度,將其與傳統(tǒng)的牛頓算法(下文簡稱傳統(tǒng)方法)進行對比分析,同樣以上述15 層3 跨的平面鋼框架結(jié)構(gòu)為例。為避免平臺的差異影響對比分析的客觀性,本文算例均在OpenSees 平臺中進行,采用BandGeneral 求解器,計算硬件條件是Interl Xeon W-2155CPU,主頻是3.3 GHz,內(nèi)存為64 GB。

    5.1 屈服單元隔離非線性分析

    進行動力時程分析之前,將豎向重力荷載分10 步逐漸加到結(jié)構(gòu)上。圖4 給出8 度和9 度罕遇地震作用下整體結(jié)構(gòu)中隔離單元占比時程曲線。從圖4 中可以看出,結(jié)構(gòu)最終屈服單元數(shù)占比分別是4.4%和9.6%;8 度罕遇地震作用下,1.62 s時刻結(jié)構(gòu)進入彈塑性狀態(tài),直到12 s 左右屈服單元數(shù)達到最大;9 度罕遇地震作用下,1.46 s 時刻進入屈服狀態(tài),直到12 s 左右塑性區(qū)不再擴展。

    圖4 隔離單元所占比例時程曲線Fig. 4 Isolated element percentage time-history curves

    圖5 給出該平面鋼結(jié)構(gòu)在8 度和9 度罕遇地震作用下的塑性區(qū)發(fā)展歷程及大小,圖中數(shù)值1~6 表示有1 個~6 個屈服單元。由圖5 可知,8 度罕遇地震下,1 層~4 層梁端先進入屈服狀態(tài),隨后,5 層~9 層梁端出現(xiàn)塑性區(qū),其中,1 層~4 層梁端塑性區(qū)面積較大;9 度罕遇地震下,1 層~6 層梁端先進入屈服狀態(tài),然后,7 層~11 層梁端和底層柱端出現(xiàn)塑性區(qū),其中,1 層~5 層梁端塑性區(qū)擴展面積較大,底層柱端塑性區(qū)面積相對較小。

    圖5 塑性區(qū)發(fā)展歷程及大小Fig. 5 Development history and size of plastic zone

    5.2 準確性和可靠性分析

    作為對比,采用傳統(tǒng)方法分析時,該結(jié)構(gòu)梁柱構(gòu)件同樣劃分30 個非線性纖維梁單元。傳統(tǒng)方法和本文方法得到的結(jié)構(gòu)頂點位移和基底剪力時程曲線如圖6~圖9 所示,可以看出,兩種方法的計算結(jié)果幾乎一致。

    圖6 8 度罕遇地震鋼結(jié)構(gòu)頂層位移時程曲線Fig. 6 Time-history curves of top displacement for the steel structure under 8 degree rare earthquakes

    圖7 9 度罕遇地震鋼結(jié)構(gòu)頂層位移時程曲線Fig. 7 Time-history curves of top displacement for the steel structure under 9 degree rare earthquakes

    圖8 8 度罕遇地震鋼結(jié)構(gòu)基底剪力時程曲線Fig. 8 Time-history curves of base reaction for the steelstructure under 8 degree rare earthquakes

    圖9 9 度罕遇地震鋼結(jié)構(gòu)基底剪力時程曲線Fig. 9 Time-history curves of base reaction for the steelstructure under 9 degree rare earthquakes

    進一步,表2 對比兩種方法的最大頂點位移和最大基底剪力以及它們的相對誤差。由表2 可知,8 度和9 度罕遇地震作用下,結(jié)構(gòu)最大頂層位移和基底剪力的相對誤差均不超過0.2%。由式(28)計算得到頂層位移時程曲線的均方差分別為0.14%和0.04%,基底剪力時程曲線的均方差分別為0.02%和0.01%。進一步驗證了本文方法的準確性、可靠性。

    表2 傳統(tǒng)方法與本文方法的誤差Table 2 Errors using the traditional and proposed methods

    式中:N為總數(shù)據(jù)量;Xi和Yi為第i時刻響應(yīng)值。

    5.3 計算效率分析

    8 度和9 度罕遇地震作用下,傳統(tǒng)方法和本文方法的迭代次數(shù)時程曲線分別如圖10 和圖11 所示。此外,傳統(tǒng)方法、本文方法以及不動點迭代數(shù)值子結(jié)構(gòu)方法(簡稱原方法)的計算信息如表3所示。由圖10 和圖11 可知,本文方法的迭代收斂次數(shù)略高于傳統(tǒng)方法,當?shù)卣鸷奢d水平增大或迭代收斂誤差減小,本文方法的迭代收斂次數(shù)變化較小。由表3 可以看出,當收斂誤差為1.0×10-5mm時,8 度和9 度罕遇地震作用下,傳統(tǒng)方法最多分別需要4 步和5 步收斂,平均迭代次數(shù)分別為2.7和3.1,本文方法最多分別需要6 步和7 步收斂,平均迭代次數(shù)分別為3.2 和3.7,而原方法最多分別需要24 步和26 步收斂,平均迭代次數(shù)分別為4.6和6.9;當收斂誤差為1.0×10-8mm 時,8 度和9 度罕遇地震作用下,傳統(tǒng)方法最多分別需要5 步和6 步收斂,平均迭代次數(shù)分別為3.2 和3.5,本文方法最多分別需要6 步和7 步收斂,平均迭代次數(shù)分別為3.9 和4.1,而原方法最多分別需要41 步和44 步收斂,平均迭代次數(shù)分別為7.1 和11.1。由此可見,本文方法所需的荷載步最大收斂次數(shù)和平均迭代次數(shù)略大于傳統(tǒng)方法,但明顯小于文獻[10 - 13]中的不動點迭代數(shù)值子結(jié)構(gòu)方法。

    表3 不同方法迭代信息統(tǒng)計表Table 3 Table of iteration information using various methods

    圖10 8 度罕遇地震迭代步數(shù)時程曲線Fig. 10 Iteration number time-history curves under an 8-degree rare earthquake

    圖11 9 度罕遇地震迭代步數(shù)時程曲線Fig. 11 Iteration number time-history curves under a 9-degree rare earthquake

    當收斂誤差為1.0×10-5mm 時,8 度和9 度罕遇地震下傳統(tǒng)方法和本文方法的計算信息統(tǒng)計表如表4 所示。可以看出:傳統(tǒng)方法分別需要對整體結(jié)構(gòu)剛度K進行7282 次和8267 次集成和分解,計算耗時分別為313 min 和398 min,而本文方法只需在初始化過程中進行一次集成和分解,非線性分析過程中僅需對Hs矩陣進行8409 次和9743 次集成和分解,計算耗時分別為157 min和332 min。

    表4 傳統(tǒng)方法與本文方法的計算信息統(tǒng)計表Table 4 Computation information using the traditional and proposed methods

    圖12 給出了傳統(tǒng)方法和本文方法求解矩陣維數(shù)的時程曲線。從圖12 可以看出:隨著地震荷載水平的增大,塑性自由度增多,迭代更新和分解的Hs矩陣規(guī)模變大;8 度和9 度罕遇地震下,本文方法求解的最大Hs矩陣的維數(shù)分別為666 和1713,而傳統(tǒng)方法求解的K矩陣維數(shù)高達9315。由此可見,結(jié)構(gòu)非線性程度越低,本文方法求解的Hs矩陣規(guī)模越小,計算效率越高。

    圖12 矩陣運算規(guī)模時程曲線Fig. 12 Matrix operation size time-history curves

    基于變剛度的傳統(tǒng)方法需要重復(fù)集成和分解整體結(jié)構(gòu)的剛度矩陣K,以及更新所有構(gòu)件單元的狀態(tài),而采用本文方法進行非線性分析時,僅對剛度矩陣K進行一次運算,迭代過程中只更新規(guī)模較小的Hs矩陣,且只對數(shù)量較少的屈服構(gòu)件單元進行狀態(tài)更新??梢酝茰y,針對局部非線性的大規(guī)模結(jié)構(gòu)系統(tǒng)的彈塑性分析,本文方法較傳統(tǒng)方法具有顯著的計算優(yōu)勢。

    6 結(jié)論

    本文對文獻[10 - 13]中不動點迭代數(shù)值子結(jié)構(gòu)算法作進一步改進,提出了降階的牛頓迭代數(shù)值子結(jié)構(gòu)方法,并利用該方法對一平面15 層3 跨鋼框架結(jié)構(gòu)進行地震彈塑性時程分析,得出以下結(jié)論:

    (1)改進的降階牛頓迭代數(shù)值子結(jié)構(gòu)方法與傳統(tǒng)牛頓算法得到的結(jié)構(gòu)頂層位移和基底剪力時程曲線一致,從而驗證了本文方法的準確性和可靠性。

    (2)本文方法收斂速率遠高于不動點迭代數(shù)值子結(jié)構(gòu)方法,接近傳統(tǒng)牛頓算法的二次收斂。當?shù)卣鸷奢d水平增加或迭代收斂條件變得苛刻,不動點迭代數(shù)值子結(jié)構(gòu)方法的收斂次數(shù)明顯增多,而本文方法卻沒有顯著變化。

    (3)本文方法的計算效率高于傳統(tǒng)牛頓算法。對于局部非線性結(jié)構(gòu)系統(tǒng),結(jié)構(gòu)模型塑性自由度規(guī)模較小,本文方法求解的Hs矩陣維度遠小于傳統(tǒng)方法求解的K矩陣,計算效率較高,且隨著結(jié)構(gòu)規(guī)模的增加,計算效率優(yōu)勢將更加明顯。

    (4)數(shù)值子結(jié)構(gòu)方法可以有效提高局部非線性結(jié)構(gòu)系統(tǒng)的模擬精度。隔離屈服構(gòu)件單元并建立精細化子結(jié)構(gòu)有限元模型,既可以精確捕捉局部屈服構(gòu)件的非線性行為,又可以有效揭示整體結(jié)構(gòu)的失效路徑和塑性分布。

    猜你喜歡
    子結(jié)構(gòu)屈服塑性
    基于應(yīng)變梯度的微尺度金屬塑性行為研究
    牙被拔光也不屈服的史良大律師秘書
    紅巖春秋(2022年1期)2022-04-12 00:37:34
    完全對換網(wǎng)絡(luò)的結(jié)構(gòu)連通度和子結(jié)構(gòu)連通度
    硬脆材料的塑性域加工
    鈹材料塑性域加工可行性研究
    The Classic Lines of A Love so Beautiful
    勇敢
    鋼框架腹板雙角鋼連接梁柱子結(jié)構(gòu)抗倒塌性能分析
    石英玻璃的熱輔助高效塑性域干磨削
    百折不撓
    丝瓜视频免费看黄片| 色播在线永久视频| 亚洲成人免费电影在线观看| 国产精品麻豆人妻色哟哟久久| 亚洲第一青青草原| 久久中文字幕人妻熟女| 三级毛片av免费| 黄色丝袜av网址大全| 香蕉久久夜色| 乱人伦中国视频| 老汉色av国产亚洲站长工具| 成人手机av| 亚洲第一欧美日韩一区二区三区 | 女性被躁到高潮视频| 黑人巨大精品欧美一区二区蜜桃| 性高湖久久久久久久久免费观看| 欧美成狂野欧美在线观看| 午夜成年电影在线免费观看| 国产精品99久久99久久久不卡| 精品亚洲成国产av| 免费观看av网站的网址| 亚洲精品乱久久久久久| 在线十欧美十亚洲十日本专区| 大型黄色视频在线免费观看| 精品少妇黑人巨大在线播放| 国产熟女午夜一区二区三区| 淫妇啪啪啪对白视频| 777久久人妻少妇嫩草av网站| 欧美精品高潮呻吟av久久| 人人妻人人澡人人爽人人夜夜| 波多野结衣一区麻豆| 久久ye,这里只有精品| 精品福利观看| 1024香蕉在线观看| 久久精品亚洲av国产电影网| 两个人看的免费小视频| 欧美日韩黄片免| 亚洲成人国产一区在线观看| 国产精品成人在线| 看免费av毛片| 丝瓜视频免费看黄片| 国产三级黄色录像| 无人区码免费观看不卡 | 91av网站免费观看| 中文字幕色久视频| 欧美日韩av久久| 日韩精品免费视频一区二区三区| 国产黄色免费在线视频| 天堂8中文在线网| 国产视频一区二区在线看| 80岁老熟妇乱子伦牲交| 欧美亚洲 丝袜 人妻 在线| 高清在线国产一区| 另类精品久久| 日韩熟女老妇一区二区性免费视频| 69精品国产乱码久久久| 精品少妇内射三级| 女人高潮潮喷娇喘18禁视频| 狠狠婷婷综合久久久久久88av| 亚洲人成伊人成综合网2020| 成人永久免费在线观看视频 | √禁漫天堂资源中文www| 欧美性长视频在线观看| 国产高清激情床上av| 99国产精品一区二区蜜桃av | 国产老妇伦熟女老妇高清| tube8黄色片| 一夜夜www| 久久久久国产一级毛片高清牌| 久久精品亚洲精品国产色婷小说| 国产真人三级小视频在线观看| 久久人妻熟女aⅴ| 十分钟在线观看高清视频www| 精品卡一卡二卡四卡免费| 自线自在国产av| 国产亚洲欧美在线一区二区| 亚洲免费av在线视频| 午夜91福利影院| 亚洲伊人色综图| 国产精品自产拍在线观看55亚洲 | 99久久99久久久精品蜜桃| 久久精品国产a三级三级三级| 亚洲人成77777在线视频| 亚洲国产av影院在线观看| 少妇猛男粗大的猛烈进出视频| 韩国精品一区二区三区| 精品免费久久久久久久清纯 | www日本在线高清视频| 天天躁狠狠躁夜夜躁狠狠躁| 日韩人妻精品一区2区三区| 亚洲欧洲日产国产| 久久久久久免费高清国产稀缺| 91精品三级在线观看| 亚洲天堂av无毛| 免费在线观看日本一区| 国产深夜福利视频在线观看| 欧美在线一区亚洲| 人人澡人人妻人| 国产精品久久久久久精品古装| 十八禁网站网址无遮挡| e午夜精品久久久久久久| 夜夜骑夜夜射夜夜干| 国产在线观看jvid| 最新的欧美精品一区二区| 国产日韩一区二区三区精品不卡| 999久久久国产精品视频| 国产免费现黄频在线看| 亚洲综合色网址| 9热在线视频观看99| 人妻 亚洲 视频| 亚洲av国产av综合av卡| 国产1区2区3区精品| 午夜视频精品福利| 肉色欧美久久久久久久蜜桃| svipshipincom国产片| 女同久久另类99精品国产91| 日韩熟女老妇一区二区性免费视频| 日本撒尿小便嘘嘘汇集6| 一本久久精品| 久久天堂一区二区三区四区| 精品国产一区二区久久| 国产老妇伦熟女老妇高清| 在线永久观看黄色视频| 99精品久久久久人妻精品| 精品高清国产在线一区| 精品国产一区二区三区四区第35| 99riav亚洲国产免费| 色视频在线一区二区三区| 大型黄色视频在线免费观看| 69av精品久久久久久 | 三上悠亚av全集在线观看| 国产精品九九99| 国产免费福利视频在线观看| 成人18禁在线播放| 成人国产一区最新在线观看| 交换朋友夫妻互换小说| 看免费av毛片| 国产xxxxx性猛交| av超薄肉色丝袜交足视频| 91大片在线观看| 老熟妇仑乱视频hdxx| 一个人免费在线观看的高清视频| 黑人猛操日本美女一级片| 精品国产亚洲在线| 丝袜美腿诱惑在线| 欧美大码av| 欧美日韩亚洲高清精品| 男女免费视频国产| 免费在线观看日本一区| 日韩精品免费视频一区二区三区| 99精品久久久久人妻精品| 精品久久久精品久久久| 亚洲av美国av| 99久久精品国产亚洲精品| 热99国产精品久久久久久7| 欧美黑人精品巨大| 国产色视频综合| 在线观看免费高清a一片| 男女午夜视频在线观看| 欧美日韩成人在线一区二区| 一区二区三区国产精品乱码| 久热这里只有精品99| 免费观看a级毛片全部| 搡老熟女国产l中国老女人| 国产精品香港三级国产av潘金莲| 亚洲人成电影免费在线| 欧美日韩精品网址| 午夜福利免费观看在线| 老司机福利观看| 国产日韩欧美视频二区| 国产伦人伦偷精品视频| 9色porny在线观看| 亚洲精品在线观看二区| a在线观看视频网站| 精品高清国产在线一区| 宅男免费午夜| 人人妻人人澡人人看| 国产精品熟女久久久久浪| 蜜桃在线观看..| 热99re8久久精品国产| 757午夜福利合集在线观看| 1024香蕉在线观看| 成年女人毛片免费观看观看9 | 老熟女久久久| 国产黄频视频在线观看| 欧美性长视频在线观看| 五月开心婷婷网| 亚洲精品国产一区二区精华液| 动漫黄色视频在线观看| 最近最新中文字幕大全免费视频| 啦啦啦视频在线资源免费观看| 久久国产精品男人的天堂亚洲| 最新的欧美精品一区二区| 精品少妇内射三级| 人妻一区二区av| e午夜精品久久久久久久| 又紧又爽又黄一区二区| 少妇裸体淫交视频免费看高清 | 久久精品91无色码中文字幕| 久久精品亚洲精品国产色婷小说| 欧美黄色片欧美黄色片| 成人国语在线视频| 久久久久久人人人人人| 色精品久久人妻99蜜桃| 丁香六月欧美| 日韩欧美一区视频在线观看| 亚洲伊人久久精品综合| 亚洲成人免费电影在线观看| 欧美日韩亚洲综合一区二区三区_| 成人18禁在线播放| 国产精品久久久av美女十八| 午夜老司机福利片| 久久久久久免费高清国产稀缺| 中文字幕另类日韩欧美亚洲嫩草| 一级毛片女人18水好多| 热99久久久久精品小说推荐| 日本vs欧美在线观看视频| 黑人巨大精品欧美一区二区mp4| 国产一区二区激情短视频| 777米奇影视久久| 欧美亚洲日本最大视频资源| 亚洲精品自拍成人| 中文字幕精品免费在线观看视频| 国产三级黄色录像| 人人澡人人妻人| 日韩成人在线观看一区二区三区| 午夜精品久久久久久毛片777| 久久精品亚洲熟妇少妇任你| 亚洲精品美女久久久久99蜜臀| 一本色道久久久久久精品综合| 激情视频va一区二区三区| 日韩视频一区二区在线观看| 色尼玛亚洲综合影院| 亚洲av成人不卡在线观看播放网| 国产黄频视频在线观看| 水蜜桃什么品种好| 男女之事视频高清在线观看| 丝袜美足系列| 在线天堂中文资源库| 亚洲欧洲日产国产| 国产日韩欧美在线精品| 性少妇av在线| 丰满饥渴人妻一区二区三| 性色av乱码一区二区三区2| 青青草视频在线视频观看| 国产在线观看jvid| 国产国语露脸激情在线看| 水蜜桃什么品种好| 亚洲av日韩在线播放| 99精品在免费线老司机午夜| 国产av一区二区精品久久| 性少妇av在线| 18禁观看日本| 18在线观看网站| 超碰成人久久| 黄色视频,在线免费观看| 最近最新中文字幕大全电影3 | 一进一出好大好爽视频| 亚洲国产欧美在线一区| 男女午夜视频在线观看| 国产成人精品无人区| 国产av国产精品国产| av视频免费观看在线观看| 亚洲人成伊人成综合网2020| 三上悠亚av全集在线观看| 又大又爽又粗| 777久久人妻少妇嫩草av网站| 天天添夜夜摸| 国产免费视频播放在线视频| 另类亚洲欧美激情| av网站在线播放免费| 乱人伦中国视频| 日本av免费视频播放| 老熟妇乱子伦视频在线观看| 国产伦人伦偷精品视频| 久久精品国产综合久久久| 天堂动漫精品| 久久国产精品人妻蜜桃| 久久天堂一区二区三区四区| 宅男免费午夜| 久久这里只有精品19| 亚洲av片天天在线观看| 精品乱码久久久久久99久播| av又黄又爽大尺度在线免费看| 国产精品熟女久久久久浪| 国产欧美亚洲国产| 精品一区二区三区av网在线观看 | 日本撒尿小便嘘嘘汇集6| 国产精品久久久久久精品古装| 久久免费观看电影| 精品福利观看| 精品福利永久在线观看| 在线 av 中文字幕| 国产免费福利视频在线观看| 9191精品国产免费久久| 露出奶头的视频| 午夜免费鲁丝| 一进一出好大好爽视频| 久久久国产一区二区| 人人妻,人人澡人人爽秒播| 免费观看a级毛片全部| 黑人巨大精品欧美一区二区mp4| 久久99一区二区三区| 99九九在线精品视频| 亚洲精品国产精品久久久不卡| 国产人伦9x9x在线观看| av视频免费观看在线观看| 成年女人毛片免费观看观看9 | 人妻 亚洲 视频| 777米奇影视久久| 新久久久久国产一级毛片| 精品国产乱码久久久久久小说| 日韩一区二区三区影片| 久久免费观看电影| 久热爱精品视频在线9| 亚洲av第一区精品v没综合| 大片免费播放器 马上看| 久久天堂一区二区三区四区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲伊人色综图| 99热国产这里只有精品6| 欧美日韩黄片免| 99riav亚洲国产免费| 欧美日韩av久久| 亚洲欧美一区二区三区黑人| 免费日韩欧美在线观看| 国产伦人伦偷精品视频| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美国产精品一级二级三级| 亚洲精品久久成人aⅴ小说| 大码成人一级视频| 国产在线观看jvid| 色综合欧美亚洲国产小说| 免费黄频网站在线观看国产| 操美女的视频在线观看| 午夜免费成人在线视频| 国产成人免费观看mmmm| 丁香六月天网| 国产成人欧美在线观看 | 久久中文字幕人妻熟女| 国产成+人综合+亚洲专区| 90打野战视频偷拍视频| av国产精品久久久久影院| 深夜精品福利| 久久99热这里只频精品6学生| 国产伦理片在线播放av一区| 国产成人av激情在线播放| 免费在线观看完整版高清| 久久ye,这里只有精品| 久久久精品免费免费高清| 欧美中文综合在线视频| 在线永久观看黄色视频| 最新美女视频免费是黄的| 久久精品aⅴ一区二区三区四区| av网站在线播放免费| 巨乳人妻的诱惑在线观看| 可以免费在线观看a视频的电影网站| 精品人妻1区二区| 午夜老司机福利片| 男女午夜视频在线观看| 日韩中文字幕视频在线看片| 少妇精品久久久久久久| 亚洲第一欧美日韩一区二区三区 | 99久久国产精品久久久| 黑人巨大精品欧美一区二区蜜桃| 免费不卡黄色视频| 国产成人系列免费观看| 777米奇影视久久| 啦啦啦免费观看视频1| 久久久久国产一级毛片高清牌| 国产精品久久久av美女十八| 亚洲国产中文字幕在线视频| 日韩欧美一区二区三区在线观看 | 免费一级毛片在线播放高清视频 | 欧美性长视频在线观看| 精品福利永久在线观看| 黑人巨大精品欧美一区二区mp4| 亚洲一区二区三区欧美精品| 动漫黄色视频在线观看| 中文字幕人妻熟女乱码| 人人妻,人人澡人人爽秒播| 极品人妻少妇av视频| 老熟女久久久| 久久久久视频综合| 国产又色又爽无遮挡免费看| 高清在线国产一区| 国产精品二区激情视频| 国产1区2区3区精品| 国产免费现黄频在线看| 亚洲专区字幕在线| 男女下面插进去视频免费观看| www日本在线高清视频| 大型av网站在线播放| 亚洲午夜精品一区,二区,三区| av网站免费在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲精品在线观看二区| 中文字幕色久视频| 精品人妻在线不人妻| 国产人伦9x9x在线观看| 精品人妻熟女毛片av久久网站| 色视频在线一区二区三区| 在线天堂中文资源库| 99国产精品99久久久久| 中文字幕人妻丝袜制服| 757午夜福利合集在线观看| 欧美在线一区亚洲| 亚洲成人免费电影在线观看| 99香蕉大伊视频| 国产欧美日韩一区二区三区在线| 99热网站在线观看| 涩涩av久久男人的天堂| 一级片'在线观看视频| 久久狼人影院| 91大片在线观看| www.999成人在线观看| 久久99热这里只频精品6学生| 国产一区二区激情短视频| 热re99久久国产66热| 丝袜美腿诱惑在线| 老司机在亚洲福利影院| 大香蕉久久成人网| 狠狠精品人妻久久久久久综合| 黄色片一级片一级黄色片| 午夜免费成人在线视频| 2018国产大陆天天弄谢| 1024香蕉在线观看| 丁香六月欧美| 免费观看av网站的网址| 少妇 在线观看| 两性夫妻黄色片| 国产不卡av网站在线观看| 日韩欧美国产一区二区入口| 99热国产这里只有精品6| 好男人电影高清在线观看| 免费不卡黄色视频| av欧美777| 制服人妻中文乱码| 女性生殖器流出的白浆| 99国产精品一区二区三区| 人成视频在线观看免费观看| 精品国产超薄肉色丝袜足j| 精品熟女少妇八av免费久了| 无遮挡黄片免费观看| 欧美中文综合在线视频| 精品少妇黑人巨大在线播放| 午夜日韩欧美国产| 亚洲av第一区精品v没综合| 国产精品一区二区精品视频观看| 日韩大码丰满熟妇| 91字幕亚洲| 欧美黄色淫秽网站| 午夜福利一区二区在线看| 人人妻人人添人人爽欧美一区卜| 国产精品自产拍在线观看55亚洲 | 国产不卡一卡二| av天堂久久9| 黄色视频不卡| 超碰成人久久| 亚洲五月婷婷丁香| 人人妻,人人澡人人爽秒播| 精品欧美一区二区三区在线| 久久中文看片网| 久久国产亚洲av麻豆专区| 狠狠婷婷综合久久久久久88av| 777久久人妻少妇嫩草av网站| 欧美乱妇无乱码| 日韩制服丝袜自拍偷拍| 熟女少妇亚洲综合色aaa.| 日韩欧美三级三区| netflix在线观看网站| 男女下面插进去视频免费观看| 亚洲成人免费电影在线观看| 亚洲国产av新网站| 精品国产一区二区三区四区第35| 久久人人爽av亚洲精品天堂| 色94色欧美一区二区| 免费日韩欧美在线观看| 久久精品亚洲av国产电影网| 午夜免费成人在线视频| 黑人猛操日本美女一级片| 国产欧美日韩精品亚洲av| 欧美乱码精品一区二区三区| 久久精品91无色码中文字幕| 18禁观看日本| 亚洲avbb在线观看| 日韩三级视频一区二区三区| 无限看片的www在线观看| 久久久久国内视频| 嫩草影视91久久| 亚洲综合色网址| 欧美黑人欧美精品刺激| 91九色精品人成在线观看| 久久亚洲真实| 狠狠精品人妻久久久久久综合| 老司机午夜福利在线观看视频 | 国产精品免费大片| 成年人免费黄色播放视频| 高清欧美精品videossex| 18在线观看网站| 日韩制服丝袜自拍偷拍| 多毛熟女@视频| 黄色视频不卡| 国产精品国产高清国产av | 一级毛片女人18水好多| 国产伦人伦偷精品视频| 精品一区二区三区视频在线观看免费 | 日韩成人在线观看一区二区三区| 最新美女视频免费是黄的| 免费高清在线观看日韩| 麻豆成人av在线观看| 操出白浆在线播放| 国产淫语在线视频| 中文字幕av电影在线播放| 亚洲av第一区精品v没综合| 日本黄色视频三级网站网址 | 无限看片的www在线观看| 99riav亚洲国产免费| av天堂在线播放| av一本久久久久| 91成年电影在线观看| 国产精品99久久99久久久不卡| 精品国产乱码久久久久久小说| 99re在线观看精品视频| 一级毛片女人18水好多| 满18在线观看网站| 丝袜在线中文字幕| 日韩中文字幕视频在线看片| 一二三四社区在线视频社区8| av不卡在线播放| 中文欧美无线码| 亚洲国产毛片av蜜桃av| 高潮久久久久久久久久久不卡| 黄频高清免费视频| 美女视频免费永久观看网站| 91麻豆av在线| 啦啦啦 在线观看视频| 美女国产高潮福利片在线看| 日韩有码中文字幕| 人妻 亚洲 视频| 丝袜美足系列| 丝袜人妻中文字幕| 国产日韩欧美亚洲二区| 欧美久久黑人一区二区| 精品免费久久久久久久清纯 | 成年人免费黄色播放视频| 岛国在线观看网站| 国产熟女午夜一区二区三区| 日韩精品免费视频一区二区三区| 国产1区2区3区精品| 亚洲国产欧美一区二区综合| 国产深夜福利视频在线观看| 久久ye,这里只有精品| 国产精品av久久久久免费| 成人特级黄色片久久久久久久 | 18禁观看日本| 亚洲精品国产色婷婷电影| 俄罗斯特黄特色一大片| 午夜两性在线视频| 少妇的丰满在线观看| 精品少妇黑人巨大在线播放| 岛国在线观看网站| bbb黄色大片| 久久亚洲精品不卡| 色综合婷婷激情| 高清毛片免费观看视频网站 | 久久精品亚洲精品国产色婷小说| 如日韩欧美国产精品一区二区三区| tube8黄色片| 国产99久久九九免费精品| 天天操日日干夜夜撸| 18禁观看日本| 国产高清国产精品国产三级| 国内毛片毛片毛片毛片毛片| 久久婷婷成人综合色麻豆| 亚洲avbb在线观看| 91成年电影在线观看| 中文亚洲av片在线观看爽 | 在线十欧美十亚洲十日本专区| 中文字幕人妻熟女乱码| 少妇被粗大的猛进出69影院| av电影中文网址| 国产欧美日韩一区二区精品| 国产高清激情床上av| 午夜91福利影院| tocl精华| 久久精品aⅴ一区二区三区四区| 99九九在线精品视频| 久久午夜亚洲精品久久| h视频一区二区三区| 一区二区三区国产精品乱码| 超色免费av| 国产xxxxx性猛交| 中文字幕高清在线视频| 一区二区三区精品91| 美女午夜性视频免费| 99国产精品99久久久久| 在线观看免费日韩欧美大片| 精品国产一区二区三区四区第35| 精品少妇内射三级| 久久性视频一级片| 久久精品熟女亚洲av麻豆精品| 热re99久久国产66热| 青草久久国产| 欧美黑人精品巨大| 美国免费a级毛片| av福利片在线| 久久天堂一区二区三区四区| 热99国产精品久久久久久7| 午夜老司机福利片| 99精品欧美一区二区三区四区| 下体分泌物呈黄色| 亚洲成a人片在线一区二区| 国产精品二区激情视频| av免费在线观看网站|