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

    向量式有限元面內(nèi)逆向轉角計算方法

    2019-06-13 09:29:56王雙炎李聰楊玉龍龐志成湯晗青
    中南大學學報(自然科學版) 2019年5期
    關鍵詞:剛體質點轉角

    王雙炎,李聰,楊玉龍,龐志成,湯晗青

    ?

    向量式有限元面內(nèi)逆向轉角計算方法

    王雙炎,李聰,楊玉龍,龐志成,湯晗青

    (浙江大學 建筑工程學院,浙江 杭州,310058)

    為提高向量式有限元方法的計算效率,以三角形平面固體單元為例,提出一種用邊轉角法計算逆向運動的面內(nèi)轉角的方法,并將其與常用的質心轉角方法進行對比,比較2種方法的計算效率;通過自編程序,分別驗證邊轉角法在計算單元剛體運動、小變形、大變形及復雜結構的變形時的準確性和計算效率。研究結果表明:與質心轉角法相比,采用邊轉角法計算面內(nèi)逆向轉角具有較高的準確度和計算效率,該方法對于常見結構面內(nèi)逆向轉角的計算是可行的、有效的。

    向量式有限元;逆向運動;平面三角形單元;計算準確度;計算效率

    向量式有限元[1-3](vector form intrinsic finite element, VFIFE)是基于向量力學和數(shù)值計算所提出的新型有限元數(shù)值計算方法,其核心思想在于將構件離散成若干個質點,通過描述各質點的運動從而對構件整體進行描述。該方法的計算特點在于通過時間積分實現(xiàn)對各點各步的逐步循環(huán)計算。由于不存在單元剛度矩陣和矩陣奇異等問題,又無需求解繁雜的非線性方程組,所以,相比傳統(tǒng)有限元方法,向量式有限元在計算特殊的非線性問題有著獨特的優(yōu)勢。王震等[4-6]從單元形式入手,將向量式有限元理論推廣至薄板、薄殼以及四節(jié)點實體單元。HOU等[7]在此基礎上,推導了八節(jié)點正方體單元的向量式有限元應用。在向量式有限元非線性行為研究中,如有關斷裂[8-9]、碰撞[10?11]等方面的研究,目前也有了一定的進展。但由于目前還沒有商業(yè)的向量式有限元軟件,一般需要通過Matlab程序實現(xiàn)向量式分析計算,而Matlab自身的運行速度具有一定的局限性。而且越是復雜的結構,劃分的單元數(shù)越多,所需的計算總時間也越多,因此,提高計算效率就顯得尤為重要。YANG等[12]在其研究中闡述了向量式有限元的計算效率。傳統(tǒng)的向量式有限元一般采用質心轉角法計算逆向轉動角度。為了提高計算效率,本文以理論體系較為成熟的三角形平面固體單元[13-15]為例,采用邊轉角法計算逆向轉動角度。HOU等[7,16]雖使用了邊轉角法進行計算,但未對該方法在計算過程中的計算準確度和計算效率進行詳細討論。為此,本文作者對單元逆向運動過程中,邊角轉角法與質心轉角法2種面內(nèi)轉角的計算方法的計算準確度和計算效率進行比較分析;通過算例驗證使用邊轉角法計算面內(nèi)逆向轉動角度的準確性及其在計算效率方面所具有的優(yōu)勢。

    1 向量式有限元平面固體單元理論

    1.1 點值描述

    傳統(tǒng)有限元采用1個有限個單元的集合體來描述求解的結構,而向量式有限元則采用點值描述。這些質點通過相互連接形成單元,建立了質點與單元間的拓撲關系。質點是描述構件受力變形、幾何尺寸、空間位置、邊界條件的載體。構件在運動過程中,通過描述各質點的位置和速度向量來描述固體構件的整體運動。

    每1個質點的運動都遵循牛頓第二運動定律。對于質點而言,m為質點的質量,為質點的空間位置向量,則質點的運動平衡方程可表示為

    1.2 途徑單元

    每1個質點的運動過程都是連續(xù)的,這些質點的運動軌跡可以表示為以時間為自變量的空間位置函數(shù)。在質點運動的全運動周期過程中,用一組時間節(jié)點(0,1,2,…,t)將質點的運動軌跡離散成若干微小時間段。在a~b時間段內(nèi),若質點的運動均滿足連續(xù)的時間函數(shù),則稱2個時間節(jié)點間的時間段為途徑單元。

    時間點和途徑單元可根據(jù)結構可能產(chǎn)生的行為進行設定。每1個離散質點連接形成的單元在設定的途徑單元內(nèi)發(fā)生的變形均為小變形。當構件有大的幾何變形時,整個變形的過程被途徑單元離散為許多小變形運動過程的疊加,從而使得描述大變形變得簡單。

    1.3 逆向運動

    為了計算得到單元節(jié)點內(nèi)力,需要從單元節(jié)點的整體位移中扣除剛體位移得到節(jié)點的純變形量,進而由虛功原理獲得單元節(jié)點內(nèi)力。

    在1個途徑單元內(nèi),平面固體單元的剛體位移由剛體平動和剛體轉動2個部分組成。圖1所示為1個三角形平面固體單元在a至時刻的途徑單元內(nèi)的剛體位移變化。a時刻為單元基礎構型,此時單元節(jié)點分別為1a,2a,3a。單元在時刻運動到了圖示位置,此時單元節(jié)點為1,2和3。圖1中,為節(jié)點的剛體平移向量(=1,2,3);2和3為經(jīng)過逆向平移后的節(jié)點;為扣除剛體平移后的節(jié)點位移向量。

    圖1 剛體逆向平移

    在局部坐標系下,和a時刻節(jié)點的位置向量分別為a,選取任意一點作為參考點(本文選取1號節(jié)點),將單元平移,使得2個對應參考點位置重合,則參考點經(jīng)過的逆向平移向量即為單元逆向剛體平移向量:

    由此得到扣除剛體平移后的節(jié)點變形量為

    式中:=2,3。

    對單元進行逆向剛體轉動。假設逆向轉動角度為,則各節(jié)點的剛體轉動位移為

    式中:=2,3。

    1.4 節(jié)點內(nèi)力計算

    最后,經(jīng)過正向運動,將局部坐標系下的單元節(jié)點內(nèi)力轉換至整體坐標系下:

    式中:=1,2,3。獲得單元節(jié)點內(nèi)力后,即可由中央差分公式實現(xiàn)循環(huán)計算。

    2 逆向運動方式探究

    2.1 面內(nèi)轉角的計算方法與計算精度

    由本文第1節(jié)的分析可知:求解逆向運動過程中的逆向轉動角度至關重要,其計算準確度和計算效率將影響整體計算結果和計算效率。文獻[17?18]中,逆向運動過程的面內(nèi)轉角估算值均采用了質心對應轉角取平均值的方法進行計算(以下簡稱質心轉角法),如圖2所示。其中,aO和dO分別為單元在a和時刻的質心位置向量;aO與dO分別為a和時刻單元各節(jié)點相對于質心的方向向量。(=1,2,3)為節(jié)點相對于自身平面單元形心的轉動角度。

    式中:=1,2,3。

    在逆向轉動面內(nèi),取3個節(jié)點轉動角度的平均值作為逆向轉動轉角:

    “面向傳動裝置的ABB AbilityTM 狀態(tài)監(jiān)測解決方案也是本次ACW的亮點之一?!盇BB中國機器人及運動控制事業(yè)部負責人李剛表示,“它是業(yè)內(nèi)首個集成化服務,能夠將每臺設備的關鍵運行參數(shù)集中顯示,借助ABB的技術優(yōu)勢,它還能使客戶提前了解維護需求,確保設備實現(xiàn)理想運行狀態(tài)??偠灾?,借助ABB AbilityTM 狀態(tài)監(jiān)測服務,用戶可以更好地掌握如何優(yōu)化設備運行,減少宕機風險,延長設備壽命,降低成本并且增加收益?!?/p>

    經(jīng)逆向運動扣除了單元節(jié)點位移中的單元剛體平動和剛體轉動后,a基礎構架時刻和d時刻虛擬單元間的節(jié)點位移僅包含單元純變形和殘余轉動變形。這個殘余轉動變形來源于單元變形的不均勻分布。而向量式有限元在設定途徑單元時,保證每一個途徑單元內(nèi)單元變形均為小變形,單元變形接近均勻變形,這就使得單元純變形對節(jié)點內(nèi)力影響為一階量,而殘余變形的影響相對純變形的影響為高階量。

    因此,可采用任意方法對面內(nèi)轉動角度進行估算,只需要滿足以下基本條件:假設單元的純變形為0時(單元僅存在剛體運動),估算方法得到的轉角應為正確的轉動角度。例如,當固體平面單元發(fā)生角度為的剛體轉動且單元變形量為零時,根據(jù)質心轉角估算法,1=2=3=,=成立,是正確的剛體轉角。

    基于上述理論,本文建議使用另1種方法估算面內(nèi)轉角,即邊轉角計算法(見圖2(b))。三角形平面固體單元經(jīng)歷了逆向平動后,參考點1已重合。此時,需進一步扣除逆向剛體轉動。與逆向剛體平動一樣,任取包含作為逆向剛體平移的參考點的單元一邊作為參考邊,以參考邊的轉動角度作為剛體逆向轉動角度。圖2中選取的邊為節(jié)點1和3所形成的單元邊,則有:

    式中:1為單元在a時刻節(jié)點1和節(jié)點3所形成的邊向量;2為單元在時刻經(jīng)逆向剛體平移后,節(jié)點1和節(jié)點3所形成的方向向量。假設單元僅發(fā)生剛體運動而沒有純變形,且單元剛體轉動角度為,此時,根據(jù)邊角計算法,=成立,滿足面內(nèi)轉角估算方法的基本條件。

    圖2 面內(nèi)逆向轉角計算方法對比

    2.2 不同計算方法的計算效率

    向量式有限元方法在求解式(1)的運動方程時,通常采用顯式時間積分法中的中央差分法。

    無初始條件下的中央差分表達式為

    有初始條件下的中央差分表達式為

    式中:+1為第+1步的節(jié)點位置向量;為質點質量;分別表示節(jié)點所受到的外力和節(jié)點內(nèi)力;為時間步長;1=1/(1+/2);2=1/(1?/2)。

    除此之外,在向量式有限元中,還有一項重要的工作即估算臨界步長。在平面三角形固體單元內(nèi),材料軸向應力剛度一般較高,分析單元軸向運動時所需的時間步長較短,因此,一般以單元軸向運動為基準對計算過程中的臨界步長進行考量。將單元簡化成長度為的質點和勁度系數(shù)為的彈簧單自由度體系,材料的彈性模量為,橫截面積為,在外力作用下,構件長度變化量為,則有:

    =(23)

    由此可見,在材料參數(shù)相同的情況下,離散質點數(shù)越多,質點連接而成的單元網(wǎng)格越小,所能取得的臨界步長越小。向量式有限元作為一種在處理碰撞、斷裂、接觸等結構大變形問題時具有獨特優(yōu)勢的新型數(shù)值計算方法,往往需要將構件離散成數(shù)量龐大的質點組來滿足對構件變形精確描述的要求。而在內(nèi)力計算過程中,為避免內(nèi)力誤差迅速積累而出現(xiàn)結果發(fā)散現(xiàn)象,尤其是在使用顯示時間積分求解運動方程的情況下,時間增量很小,迭代計算的循環(huán)步數(shù)一般很龐大,此時,計算效率的提升就顯得尤為重要。

    圖3所示為不同計算方法的計算流程圖。相較于傳統(tǒng)的質心轉角計算方法,邊轉角計算法無需計算單元質心位置,向量計算量也明顯較少。對于1個時間步長內(nèi)的1個單元而言,在面內(nèi)逆向轉角計算步驟中,質心轉角算法需完成10次向量加減運算、8次向量乘法運算、3次向量點乘運算、6次向量求模運算、2次標量加減運算以及1次標量乘法運算;邊轉角法僅需完成2次向量加減運算與1次向量乘法運算。由計算機工作原理可知,邊轉角計算法可以有效提高計算機的運行速度。

    3 算例驗證

    從上述分析可以看出,取質心與節(jié)點的向量轉角平均值作為逆向轉動角度和以邊向量轉角為基準作為逆向轉動角度均是可行的。同時,在計算效率方面,由于邊轉角法的計算量更小,在多質點多單元的情況下更具優(yōu)勢。本節(jié)通過自編Matlab程序進行算例模擬,對上述邊轉角理論的正確性進行驗證。

    3.1 單元剛體運動

    由表1可知:2種方法在處理單元剛體運動問題時,均具有較高的準確性。當單元旋轉10圈時,質心轉角法所用時間為2.7 s,邊轉角法所用時間為2.2 s。

    3.2 正方形薄板受壓問題

    正方形薄板示意圖及單元劃分如圖5所示。由圖5(a)可見:正方形薄板對角受到均布壓力荷載的作用,荷載沿垂直于平面的薄板厚度方向均勻分布,為2 N/m。正方形薄板對角線長=4 m,板厚0.1 m,彈性模量=1×1011Pa,泊松比=0。

    由于結構具有對稱性,只需取1/4部分作為分析對象。由圖5(b)可見:構件的1/4被離散成6個質點,質點間互相連接形成共4個單元。將向量式有限元分析結果與傳統(tǒng)有限元分析結果進行對比,結果如表2所示。由表2可以看出:在小變形情況下,2種面內(nèi)轉角計算方式的向量式有限元分析結果均與傳統(tǒng)有限元的理論值較符合。質心轉角法所用時間為2.2 s,邊轉角法所用時間為1.8 s。

    圖3 不同計算方法流程圖

    圖4 單元剛體轉動

    圖5 正方形薄板示意圖及單元劃分

    3.3 懸臂梁自由端受集中荷載作用

    在Euler?Benoulli梁理論[19]中,懸臂梁端部受集中荷載作用發(fā)生劇烈變形是經(jīng)典的非線性有限元算例,該算例常被用于驗證非線性算法理論的正確 性[20?23]。懸臂梁算例示意圖如圖6所示。懸臂梁長=10 m,梁高=1 m,梁厚為1 m。構件離散成544個質點,以三角形平面單元形式相互連接。材料彈性模量=10 Pa,泊松比=0。當外荷載以荷載因子/2(其中,為懸臂梁橫截面的慣性矩)來衡量時,此時懸臂梁變形為大變形。由于該算例是擬靜力問題,荷載加載形式取斜坡加載方式。同時,施加虛擬阻尼使結果收斂,取虛擬阻尼系數(shù)=10.0。使用質心轉角計算法和邊轉角計算法分別進行計算,驗證大變形情況下理論的正確性。

    懸臂梁荷載施加點豎向位移和橫向位移對比分別如表3和表4。表3和表4中,ABAQUS的單元劃分方式均與算例的相同。由表3和表4可知:對于質心轉角法和邊轉角法2種不同的逆向轉動計算方法,其自由端在豎直和水平方向上的位移均與理論結果較吻合,2種方法得到的結果之間的差異也很小。即使是當荷載因子達到10這樣的高度非線性情況下,依然能取得理想結果,驗證了本文第3節(jié)提出的邊轉角計算方法的計算精確性。在計算效率方面,質點轉角法最終所用時長為878 s,而在相同情況下,邊轉角法所用時長為779 s。

    表1 節(jié)點坐標對比

    表2 節(jié)點位移對比

    圖6 懸臂梁算例示意圖

    表3 懸臂梁荷載施加點豎向位移對比

    表4 懸臂梁荷載施加點水平位移對比

    3.4 半圓弧結構受力大變形

    圖7所示為半圓弧結構示意圖,圓弧底部兩端完全固定。圓弧外徑1=21 m,內(nèi)徑2=20 m,材料彈性模量為10 Pa,泊松比為0.25。將結構分別離散成108,318和860個質點進行計算。半圓弧頂部位置施加1個豎直向上的集中荷載。由于該問題是1個擬靜力問題,荷載加載方式同樣采用斜坡加載方式。取虛擬阻尼系數(shù)=10.0。構件將在荷載作用下發(fā)生 大變形,驗證工程結構條件下邊轉角計算方法的正 確性。

    為驗證除荷載施加點之外,圓弧拱所離散成的其他節(jié)點位移的正確性,分別選取荷載為1.0 N和3.5 N,觀察圓弧拱拱肋上沿部分其他質點產(chǎn)生合位移后連接所成的拱肋變形,如圖8所示。同時,選取荷載施加質點進行定量分析。圖9(a) 所示為圓弧結構離散為不同質點數(shù)情況下荷載施加點的位移與荷載的曲線(由于質心轉角法和邊轉角法的結果差距極小,圖8和圖9(a)僅展示了邊轉角法的結果),并將結果與商業(yè)通用有限元軟件ABAQUS分析結果進行對比(其中,圓弧結構被離散成5 992個三角形單元)。

    圖7 半圓弧結構示意圖

    荷載/N:(a) 1.0; (b) 3.5

    由圖9(a)可知:圓弧頂點將在荷載作用下先發(fā)生非線性位移變化;隨著荷載的增加,頂點兩邊的弧形結構逐漸由曲線變?yōu)橹本€,位移?荷載的變化又重新趨于線性相關;自編向量式有限元程序的分析結果與ABAQUS顯示的位移變化趨勢相吻合。隨著離散質點數(shù)的增加,頂點位移逐漸接近ABAQUS非線性有限元分析結果;且在相同的單元劃分情況下,向量式有限元計算結果也與ABAQUS計算結果相吻合。

    由上述分析可以看出:在保證計算準確度的前提下,邊轉角法較質心轉角法有著更高的計算效率。由于向量式分析的特點,為了獲得更精確的結果,需要將結構離散成更多的質點以達到更精確描述結構的目的。此時,運用邊轉角法進行計算所能節(jié)省的絕對時間會更多。

    圖9 半圓弧結構算例計算結果

    圖9(b)所示為不同質點數(shù)下Matlab程序對于2種方法所需的平均計算時間。由圖9(b)可知:對于3種不同離散情況,邊轉角法相較于質心轉角法都能節(jié)省時間12%~15%;當結構離散成860個質點時,2種方法的絕對計算時差可達約360 s。

    4 結論

    1) 為提高向量式有限元的程序計算效率,基于向量式有限元的基本理論和假定,采用邊轉角法計算單元在逆向面內(nèi)轉動過程中的逆向轉角。

    2) 通過自編的Matlab程序算例對單元剛體運動、小變形、大變形、復雜結構4種情況進行了驗算,計算結果表明邊轉角法在保證較高計算準確度的前提下,相較于傳統(tǒng)的質心轉角法計算效率更高。

    [1] TING E C, SHIH C, WANG Y K. Fundamentals of a vector form intrinsic finite element: Part 1. Basic procedure and a plane frame element[J]. Journal of Mechanics, 2004, 20(2): 113?122.

    [2] TING E C, SHIH C, WANG Y K. Fundamentals of a vector form intrinsic finite element: Part II. Plane solid elements[J]. Journal of Mechanics. 2004, 20(2): 123?132.

    [3] SHIH C, WANG Y K, TING E C. Fundamentals of a vector form intrinsic finite element: Part III. Convected material frame and examples[J]. Journal of Mechanics. 2004, 20(2): 133?143.

    [4] 王震, 趙陽, 胡可. 基于向量式有限元的三角形薄板單元[J]. 工程力學, 2014, 31(1): 37?45. WANG Zhen, ZHAO Yang, HU Ke. Triangular plate element based on vectored finite element[J]. Engineering Mechanics, 2014, 31(1): 37?45.

    [5] 王震. 向量式有限元薄殼單元的理論與應用[D]. 杭州:浙江大學建筑工程學院, 2013: 80?85. WANG Zhen. The theory and application of a vector finite element thin shell element[D]. Hangzhou: Zhejiang University. College of Civil Engineering and Architecture, 2013: 80?85.

    [6] 王震, 趙陽, 楊學林. 基于向量式有限元的實體結構非線性行為分析[J]. 建筑結構學報. 2015(3): 133?140. WANG Zhen, ZHAO Yang, YANG Xuelin. Nonlinear behavior analysis of solid structure based on vector finite element method[J]. Journal of Architectura, 2015(3): 133?140.

    [7] HOU Xiangying, FANG Zongde. Solid structure analysis with large deformation of eight-node hexahedral element using vector form intrinsic finite element[J]. Advances in Structural Engineering, 2018, 21(6): 852?861.

    [8] DUAN Y F, WANG S M, WANG R Z, et al. Vector form intrinsic finite element based approach to simulate crack propagation[J]. Journal of Mechanics, 2017, 33(6): 797?812.

    [9] HOU Xiangying, FANG Zongde, ZHANG Xijin, et al. Static contact analysis of spiral bevel gear based on modified VFIFE (vector form intrinsic finite element) method[J]. Applied Mathematical Modeling, 2018, 60: 192?207.

    [10] XU Leige, LIN Mian. Numerical modeling of the configuration of a long-distance free-spanning submarine pipeline on an uneven seabed[J]. International Journal of Offshore and Polar Engineering, 2017, 27(1): 102?111.

    [11] WU T Y. Dynamic nonlinear analysis of shell structures using a vector form intrinsic finite element[J]. Engineering Structures, 2013, 56: 2028?2040.

    [12] YANG Yulong, CHENG J J R, ZHANG Tuqiao. Vector form intrinsic finite element method for planar multibody systems with multiple clearance joints[J]. Nonlinear Dynamics, 2016, 86(1): 421?440.

    [13] TORNABENE F, FANTUZZI N, BACCIOCCHI M, et al. Mechanical behavior of damaged laminated composites plates and shells: higher-order shear deformation theories[J]. Composite Structures, 2018, 189: 304?329.

    [14] 劉人懷, 薛江紅. 復合材料層合板殼非線性力學的研究進展[J]. 力學學報, 2017, 49(3): 487?506. LIU Renhuai, XUE Jianghong. Research progress in nonlinear mechanics of laminated composite plates and shells[J]. Journal of Mechanics, 2017, 49(3): 487?506

    [15] 韓強.高等板殼理論[M]. 北京: 科學出版社, 2002: 216. HAN Qian. Advanced theory of plates and shells[M]. Beijing: Science Press, 2002: 216.

    [16] XU Leige, LIN Mian. Analysis of buried pipelines subjected to reverse fault motion using the vector form intrinsic finite element method[J]. Soil Dynamics and Earthquake Engineering, 2017, 93: 61?83.

    [17] 王震, 趙陽, 楊學林. 薄殼結構的向量式有限元屈曲行為分析[J]. 中南大學學報(自然科學版), 2016, 47(6): 2058?2064. WANG Zhen, ZHAO Yang, YANG Xuelin. Vector form intrinsic finite element method for buckling analysis of thin-shell structures[J]. Journal of Central South University of Science and Technology, 2016, 47(6): 2058?2064.

    [18] 底大鈞. 三維空間柔性桿件的三角單元向量式有限元模型研究及仿真[D]. 哈爾濱:哈爾濱工程大學理學院, 2017: 24?31. DI Dajun. Research and Simulation of triangular element vector finite element model for three-dimensional flexible members[D]. Harbin:Harbin Engineering University. College of Science, 2017: 24?31.

    [19] MATTIASSON K. Numerical results from large deflection beam and frame problems analyzed by means of elliptic integrals[J]. International Journal for Numerical Methods in Engineering. 1981, 17(1): 145?153.

    [20] 張鵬飛. 結構破壞行為的數(shù)值模擬計算方法研究[D]. 杭州:浙江大學建筑工程學院, 2016: 84?85. ZHANG Pengfei. Study on numerical simulation method for structural failure behavior[D]. Hangzhou: Zhejiang University. College of Civil Engineering and Architecture, 2016: 84?85.

    [21] 羅堯治, 楊超. 求解平面固體幾何大變形問題的有限質點法[J]. 工程力學, 2013, 30(4): 260?268. LUO Yaozhi, YANG Chao. Finite particle method for solving large deformation problems in plane solid geometry[J]. Engineering Mechanics, 2013, 30(4): 260?268.

    [22] El-Abbasi N, Meguid S A. A new shell element accounting for through-thickness deformation[J]. Computer Methods in Applied Mechanics and Engineering, 2000, 189(3): 841?862.

    [23] SZE K Y, CHAN W K, PIAN T. An eight-node hybrid-stress solid-shell element for geometric non-linear analysis of elastic shells[J]. International Journal for Numerical Methods in Engineering, 2002, 55(7): 853?878.

    Research on reverse rotation angle calculation method of vector form intrinsic finite element in plane

    WANG Shuangyan, LI Cong, YANG Yulong, PANG Zhicheng, TANG Hanqing

    (College of Civil Engineering and Architecture, Zhejiang University, Hangzhou 310058, China)

    In order to improve the calculation efficiency of vector form intrinsic finite element (VFIFE), the side-angle-method for calculation of the reverse rotation angle of the triangular plane solid element reverse motion was proposed. By comparing the side-angle-method with the commonly used centroid-angle-method, the calculation efficiencies of these two methods were studied. The accuracy and efficiency of the side-angle-method in calculating the rigid body motion, the small deformation, the large deformation and the deformation of complex structure were discussed and verified by examples using the self compiled program. The results show that compared with conventional centroid- angle-method, the side-angle-method has good computation accuracy and computational efficiency, and it is feasible and effective for calculation of the reverse rotation angle of common structures.

    vector form intrinsic finite element; reverse motion; planar triangular element; calculation accuracy; calculation efficiency

    TU33+9

    A

    1672?7207(2019)05?1135?09

    10.11817/j.issn.1672?7207.2019.05.017

    2018?06?19;

    2018?08?19

    國家水體污染控制與治理科技重大專項(2017ZX07201004)(Project(2017ZX07201004) supported by the Major Science and Technology Program for Water Pollution and Treatment)

    楊玉龍,博士,講師,從事市政工程安全分析研究;E-mail: yulongy@zju.edu.cn

    (編輯 伍錦花)

    猜你喜歡
    剛體質點轉角
    玩轉角的平分線
    巧用“搬運法”解決連續(xù)質點模型的做功問題
    差值法巧求剛體轉動慣量
    車載冷發(fā)射系統(tǒng)多剛體動力學快速仿真研究
    三次“轉角”遇到愛
    解放軍健康(2017年5期)2017-08-01 06:27:42
    質點的直線運動
    質點的直線運動
    永春堂贏在轉角
    下一個轉角:邁出去 開啟“智”造時代
    剛體定點轉動的瞬軸、極面動態(tài)演示教具
    物理實驗(2015年10期)2015-02-28 17:36:56
    免费人成在线观看视频色| 精品人妻一区二区三区麻豆| 干丝袜人妻中文字幕| 蜜桃亚洲精品一区二区三区| 草草在线视频免费看| 亚洲成人精品中文字幕电影| 黄片无遮挡物在线观看| 丰满人妻一区二区三区视频av| 自拍欧美九色日韩亚洲蝌蚪91 | 一区二区三区四区激情视频| 一二三四中文在线观看免费高清| 99九九线精品视频在线观看视频| 一级二级三级毛片免费看| 卡戴珊不雅视频在线播放| 亚洲欧美清纯卡通| 性色av一级| 嘟嘟电影网在线观看| 国产大屁股一区二区在线视频| 乱码一卡2卡4卡精品| 日韩一区二区视频免费看| 国产探花极品一区二区| 我的老师免费观看完整版| 免费观看av网站的网址| 亚洲精品视频女| 欧美日韩视频精品一区| 国产成人91sexporn| 国产成人免费观看mmmm| 爱豆传媒免费全集在线观看| 菩萨蛮人人尽说江南好唐韦庄| 老司机影院成人| 最近中文字幕2019免费版| 人人妻人人爽人人添夜夜欢视频 | 欧美精品人与动牲交sv欧美| 国产精品一区二区在线观看99| 免费av毛片视频| 成人一区二区视频在线观看| 国产成人精品一,二区| 欧美+日韩+精品| 99re6热这里在线精品视频| 大又大粗又爽又黄少妇毛片口| 亚洲va在线va天堂va国产| 深爱激情五月婷婷| 国产精品国产三级国产av玫瑰| 人妻系列 视频| 秋霞在线观看毛片| 大又大粗又爽又黄少妇毛片口| 熟女av电影| 一个人看的www免费观看视频| 欧美日韩综合久久久久久| 熟女电影av网| 建设人人有责人人尽责人人享有的 | 中文字幕人妻熟人妻熟丝袜美| 久久人人爽人人爽人人片va| 国产男人的电影天堂91| 一级二级三级毛片免费看| 国模一区二区三区四区视频| 欧美xxxx性猛交bbbb| 国产一级毛片在线| 男女那种视频在线观看| 别揉我奶头 嗯啊视频| 晚上一个人看的免费电影| 亚洲欧美日韩另类电影网站 | 国产成人精品久久久久久| 亚洲av国产av综合av卡| 精品久久久久久电影网| 免费观看性生交大片5| 一本久久精品| 午夜精品国产一区二区电影 | 不卡视频在线观看欧美| 大话2 男鬼变身卡| 九九久久精品国产亚洲av麻豆| 精品午夜福利在线看| 久久6这里有精品| 熟女人妻精品中文字幕| av在线app专区| 中国美白少妇内射xxxbb| 中国美白少妇内射xxxbb| 高清视频免费观看一区二区| 天天一区二区日本电影三级| 中文在线观看免费www的网站| 又黄又爽又刺激的免费视频.| 人妻 亚洲 视频| 色综合色国产| 国产精品无大码| 日韩av免费高清视频| 女人十人毛片免费观看3o分钟| 免费观看的影片在线观看| 秋霞在线观看毛片| 久热这里只有精品99| 夫妻性生交免费视频一级片| 日韩av在线免费看完整版不卡| av免费观看日本| 一边亲一边摸免费视频| 欧美日韩国产mv在线观看视频 | 天天躁日日操中文字幕| 99热全是精品| 人妻 亚洲 视频| 久久久久久久久久久免费av| 2022亚洲国产成人精品| 精品午夜福利在线看| 一个人观看的视频www高清免费观看| 国产亚洲精品久久久com| 狂野欧美白嫩少妇大欣赏| 婷婷色av中文字幕| 干丝袜人妻中文字幕| 日韩一区二区视频免费看| 成人亚洲精品一区在线观看 | 国产精品人妻久久久久久| 欧美最新免费一区二区三区| 久久99蜜桃精品久久| 777米奇影视久久| 亚洲婷婷狠狠爱综合网| 国产精品蜜桃在线观看| 新久久久久国产一级毛片| 色5月婷婷丁香| 国产午夜福利久久久久久| freevideosex欧美| 大片电影免费在线观看免费| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜视频国产福利| 亚洲精品成人av观看孕妇| 色播亚洲综合网| 亚洲欧洲日产国产| 国产精品久久久久久精品电影小说 | 自拍欧美九色日韩亚洲蝌蚪91 | 国产欧美日韩精品一区二区| 欧美老熟妇乱子伦牲交| 男女啪啪激烈高潮av片| 国产精品精品国产色婷婷| 能在线免费看毛片的网站| 国产精品伦人一区二区| a级一级毛片免费在线观看| 网址你懂的国产日韩在线| 夫妻性生交免费视频一级片| 日本色播在线视频| 日本爱情动作片www.在线观看| 美女xxoo啪啪120秒动态图| 亚洲国产色片| 中文字幕人妻熟人妻熟丝袜美| 搞女人的毛片| 亚洲av.av天堂| 日韩亚洲欧美综合| 免费av观看视频| 日韩伦理黄色片| 免费观看的影片在线观看| 一级爰片在线观看| 九九在线视频观看精品| 国产免费福利视频在线观看| 最后的刺客免费高清国语| 亚洲国产色片| 日韩免费高清中文字幕av| 中文字幕久久专区| 一区二区三区精品91| 亚洲色图综合在线观看| 美女被艹到高潮喷水动态| 婷婷色麻豆天堂久久| 国精品久久久久久国模美| 97在线人人人人妻| 久久久久久九九精品二区国产| 丰满乱子伦码专区| 中文字幕制服av| 精品国产露脸久久av麻豆| 久久久久久国产a免费观看| 99热全是精品| 成年版毛片免费区| 久久久精品欧美日韩精品| 精品99又大又爽又粗少妇毛片| 国产男女超爽视频在线观看| 久久久国产一区二区| 成人鲁丝片一二三区免费| www.av在线官网国产| 中国三级夫妇交换| 久久人人爽人人片av| 一个人看的www免费观看视频| 亚洲av免费在线观看| 国产淫片久久久久久久久| 日日摸夜夜添夜夜添av毛片| 搡女人真爽免费视频火全软件| 午夜老司机福利剧场| 少妇的逼水好多| 国产精品偷伦视频观看了| 久久99热这里只频精品6学生| 亚洲久久久久久中文字幕| 日本-黄色视频高清免费观看| 一边亲一边摸免费视频| 久久久久久久久久久丰满| 欧美精品一区二区大全| 国产成人a∨麻豆精品| 精品视频人人做人人爽| 草草在线视频免费看| 一本一本综合久久| 特级一级黄色大片| 日韩三级伦理在线观看| 亚洲久久久久久中文字幕| 日韩欧美 国产精品| 偷拍熟女少妇极品色| 九九爱精品视频在线观看| 亚洲欧美一区二区三区黑人 | 中文欧美无线码| 97热精品久久久久久| 国产又色又爽无遮挡免| 国产精品国产三级国产av玫瑰| 国产精品99久久99久久久不卡 | 精品久久久精品久久久| 18+在线观看网站| 一区二区av电影网| 建设人人有责人人尽责人人享有的 | 久久久精品免费免费高清| 国产午夜精品久久久久久一区二区三区| 午夜免费鲁丝| 久久这里有精品视频免费| 亚洲精品成人久久久久久| 亚洲高清免费不卡视频| 国产伦理片在线播放av一区| 九九久久精品国产亚洲av麻豆| 日韩亚洲欧美综合| 久久久久久久久久人人人人人人| 一级片'在线观看视频| 女的被弄到高潮叫床怎么办| 亚洲av日韩在线播放| 观看美女的网站| 久久久久久久久久久丰满| 午夜福利在线在线| 2021少妇久久久久久久久久久| 日韩成人伦理影院| 亚洲自偷自拍三级| 日韩一区二区三区影片| 亚洲欧美精品自产自拍| 天天躁日日操中文字幕| 一级爰片在线观看| av免费观看日本| av女优亚洲男人天堂| 亚洲最大成人中文| 爱豆传媒免费全集在线观看| 欧美日韩亚洲高清精品| 国产精品福利在线免费观看| 国产精品嫩草影院av在线观看| 在线播放无遮挡| 我的女老师完整版在线观看| 蜜桃亚洲精品一区二区三区| 1000部很黄的大片| videossex国产| 91精品一卡2卡3卡4卡| 七月丁香在线播放| 中文字幕亚洲精品专区| 亚洲欧美精品自产自拍| 麻豆国产97在线/欧美| 在线观看美女被高潮喷水网站| 国产高清不卡午夜福利| 日韩av在线免费看完整版不卡| 日韩欧美 国产精品| 精品久久国产蜜桃| 欧美xxxx性猛交bbbb| 亚洲在线观看片| 久久久久久久久久成人| freevideosex欧美| 久久97久久精品| 国产免费又黄又爽又色| 亚洲成人精品中文字幕电影| 国产精品一及| 欧美97在线视频| 国产老妇女一区| 欧美少妇被猛烈插入视频| 欧美极品一区二区三区四区| 看免费成人av毛片| 免费播放大片免费观看视频在线观看| 久久韩国三级中文字幕| 成人欧美大片| 国产精品偷伦视频观看了| 国语对白做爰xxxⅹ性视频网站| 国产免费又黄又爽又色| 国产精品偷伦视频观看了| 亚洲va在线va天堂va国产| 亚洲精华国产精华液的使用体验| 久久久午夜欧美精品| 亚洲av电影在线观看一区二区三区 | 精品久久久久久久末码| a级一级毛片免费在线观看| 日日撸夜夜添| 成人亚洲精品av一区二区| 国产爱豆传媒在线观看| 久久国内精品自在自线图片| 欧美激情久久久久久爽电影| 狂野欧美激情性xxxx在线观看| 国产乱来视频区| 精品视频人人做人人爽| 国语对白做爰xxxⅹ性视频网站| 国产人妻一区二区三区在| 久久精品国产a三级三级三级| 日本熟妇午夜| 亚洲一区二区三区欧美精品 | 99热这里只有是精品在线观看| 91精品一卡2卡3卡4卡| 最新中文字幕久久久久| 一区二区三区乱码不卡18| 国产成人精品一,二区| freevideosex欧美| 一二三四中文在线观看免费高清| 亚洲不卡免费看| 免费观看在线日韩| 综合色av麻豆| 麻豆精品久久久久久蜜桃| 人妻一区二区av| 99热国产这里只有精品6| 简卡轻食公司| 国产午夜福利久久久久久| 舔av片在线| 国产 一区 欧美 日韩| 永久免费av网站大全| 一区二区三区免费毛片| 久久久久久久久久人人人人人人| 亚洲精品久久午夜乱码| 欧美+日韩+精品| 狂野欧美激情性xxxx在线观看| 亚洲av成人精品一区久久| 美女主播在线视频| 国产一级毛片在线| av国产免费在线观看| 国产美女午夜福利| 国产女主播在线喷水免费视频网站| 亚洲国产精品成人综合色| 日日摸夜夜添夜夜添av毛片| 91久久精品国产一区二区三区| 中文在线观看免费www的网站| 免费观看a级毛片全部| 中文乱码字字幕精品一区二区三区| 成人亚洲精品一区在线观看 | 日日摸夜夜添夜夜爱| 欧美人与善性xxx| 搡老乐熟女国产| 亚洲丝袜综合中文字幕| 久久久久久国产a免费观看| 国产亚洲av片在线观看秒播厂| 少妇人妻一区二区三区视频| 日韩 亚洲 欧美在线| 亚洲欧美精品自产自拍| 亚洲精品日本国产第一区| 最近中文字幕2019免费版| 美女内射精品一级片tv| 波野结衣二区三区在线| av卡一久久| 人妻制服诱惑在线中文字幕| 最近中文字幕高清免费大全6| 国产又色又爽无遮挡免| 亚洲av电影在线观看一区二区三区 | 成人黄色视频免费在线看| 秋霞在线观看毛片| 国产黄频视频在线观看| 日本黄色片子视频| 自拍偷自拍亚洲精品老妇| 狂野欧美激情性bbbbbb| 久久久亚洲精品成人影院| 人妻少妇偷人精品九色| 男人舔奶头视频| 亚洲国产精品成人久久小说| 永久免费av网站大全| 日韩亚洲欧美综合| 22中文网久久字幕| 在线免费十八禁| 日本黄大片高清| 国产免费一区二区三区四区乱码| 99久久九九国产精品国产免费| 久久影院123| 另类亚洲欧美激情| 欧美一区二区亚洲| 99re6热这里在线精品视频| 超碰av人人做人人爽久久| 亚洲无线观看免费| 一级二级三级毛片免费看| 中文资源天堂在线| 深爱激情五月婷婷| 亚洲国产日韩一区二区| 嫩草影院入口| 久久久久久久久久久免费av| 亚洲婷婷狠狠爱综合网| 免费观看av网站的网址| 在线免费十八禁| 午夜福利高清视频| 简卡轻食公司| 蜜臀久久99精品久久宅男| 亚洲精品第二区| 男女边摸边吃奶| 亚洲成人久久爱视频| 国产免费一级a男人的天堂| 亚洲最大成人手机在线| 国产高潮美女av| 日本一二三区视频观看| 日本-黄色视频高清免费观看| 性色avwww在线观看| av福利片在线观看| 精品久久久精品久久久| 免费观看av网站的网址| 熟女人妻精品中文字幕| 亚洲精品自拍成人| 成人二区视频| 搡女人真爽免费视频火全软件| 国产欧美日韩一区二区三区在线 | videos熟女内射| 亚洲综合色惰| 国产探花在线观看一区二区| 亚洲国产欧美在线一区| 国产欧美日韩精品一区二区| 麻豆国产97在线/欧美| 国产精品麻豆人妻色哟哟久久| 色吧在线观看| 在线观看人妻少妇| 一边亲一边摸免费视频| 国产精品不卡视频一区二区| 亚洲精品久久午夜乱码| .国产精品久久| 最近手机中文字幕大全| 国产精品一区二区性色av| 国产精品99久久99久久久不卡 | 国产久久久一区二区三区| 日韩伦理黄色片| 亚洲精品一区蜜桃| 国产成人精品婷婷| 午夜视频国产福利| 日韩成人av中文字幕在线观看| 特级一级黄色大片| 欧美激情在线99| 少妇高潮的动态图| 最后的刺客免费高清国语| 亚洲第一区二区三区不卡| 免费观看的影片在线观看| 国产精品一区二区在线观看99| 亚洲高清免费不卡视频| 欧美日本视频| 久久久成人免费电影| 高清毛片免费看| 菩萨蛮人人尽说江南好唐韦庄| 日本一二三区视频观看| 九草在线视频观看| 又粗又硬又长又爽又黄的视频| 高清毛片免费看| 只有这里有精品99| 少妇 在线观看| 亚洲av日韩在线播放| av在线播放精品| 男女那种视频在线观看| 看免费成人av毛片| 我要看日韩黄色一级片| 老师上课跳d突然被开到最大视频| 午夜福利在线在线| 美女cb高潮喷水在线观看| 日韩,欧美,国产一区二区三区| 又粗又硬又长又爽又黄的视频| 美女视频免费永久观看网站| 成人毛片a级毛片在线播放| 欧美成人一区二区免费高清观看| 国产伦在线观看视频一区| 国产一区亚洲一区在线观看| 中文字幕亚洲精品专区| 国产极品天堂在线| 国产 精品1| 1000部很黄的大片| 嫩草影院入口| 欧美激情国产日韩精品一区| 国产一区二区亚洲精品在线观看| av线在线观看网站| 国产真实伦视频高清在线观看| 国产 精品1| 国产午夜精品久久久久久一区二区三区| 亚洲精品456在线播放app| 免费播放大片免费观看视频在线观看| 国产大屁股一区二区在线视频| 亚洲av免费在线观看| 国产成人a区在线观看| 91久久精品国产一区二区三区| 91狼人影院| 六月丁香七月| 麻豆久久精品国产亚洲av| 免费观看的影片在线观看| 亚洲欧美精品自产自拍| 国产中年淑女户外野战色| 欧美日韩视频精品一区| 最近最新中文字幕大全电影3| 久久人人爽人人爽人人片va| 久久久久久久久久成人| 欧美xxxx黑人xx丫x性爽| 制服丝袜香蕉在线| 三级国产精品片| 黄片wwwwww| 一级爰片在线观看| 国产精品女同一区二区软件| 亚洲av二区三区四区| 成年人午夜在线观看视频| 午夜福利在线在线| 国产精品久久久久久久久免| 亚洲怡红院男人天堂| 午夜精品一区二区三区免费看| 日本猛色少妇xxxxx猛交久久| 亚洲av中文字字幕乱码综合| 日韩成人av中文字幕在线观看| 亚洲欧美日韩东京热| 久久韩国三级中文字幕| 成人毛片a级毛片在线播放| 男人狂女人下面高潮的视频| 免费不卡的大黄色大毛片视频在线观看| 欧美97在线视频| 在线观看美女被高潮喷水网站| 看十八女毛片水多多多| 人妻少妇偷人精品九色| 欧美bdsm另类| 日韩视频在线欧美| 在线观看美女被高潮喷水网站| 少妇高潮的动态图| 下体分泌物呈黄色| 哪个播放器可以免费观看大片| 亚洲美女搞黄在线观看| 国产成年人精品一区二区| 免费人成在线观看视频色| 亚洲精品久久午夜乱码| 寂寞人妻少妇视频99o| 色婷婷久久久亚洲欧美| 身体一侧抽搐| 日日摸夜夜添夜夜爱| 精品一区二区三区视频在线| 国产毛片a区久久久久| 一本色道久久久久久精品综合| 99热这里只有是精品在线观看| 日韩欧美 国产精品| 一个人观看的视频www高清免费观看| kizo精华| 国产av国产精品国产| 亚洲欧美日韩东京热| 国产综合精华液| 日韩欧美一区视频在线观看 | 国国产精品蜜臀av免费| 99热这里只有是精品50| 免费av毛片视频| 免费黄频网站在线观看国产| 国产成人免费观看mmmm| 久久久久久久亚洲中文字幕| 亚洲综合色惰| 又黄又爽又刺激的免费视频.| 菩萨蛮人人尽说江南好唐韦庄| 熟妇人妻不卡中文字幕| 精品一区二区免费观看| 国产大屁股一区二区在线视频| 亚洲欧美成人综合另类久久久| 干丝袜人妻中文字幕| 久久精品国产自在天天线| 别揉我奶头 嗯啊视频| 国产成人免费观看mmmm| 国产大屁股一区二区在线视频| 少妇被粗大猛烈的视频| 成人欧美大片| 亚洲欧美日韩卡通动漫| 国产精品一区二区性色av| 超碰97精品在线观看| 国产精品国产三级国产av玫瑰| 日韩在线高清观看一区二区三区| 欧美成人a在线观看| 天天躁夜夜躁狠狠久久av| av线在线观看网站| 久久影院123| av.在线天堂| 久热久热在线精品观看| 一级a做视频免费观看| 国产毛片a区久久久久| 久久久久网色| 少妇人妻 视频| 赤兔流量卡办理| 午夜老司机福利剧场| 中文欧美无线码| 国产一区有黄有色的免费视频| 久久ye,这里只有精品| 高清午夜精品一区二区三区| 一级av片app| 国产男女内射视频| av天堂中文字幕网| 久久人人爽人人片av| 日日撸夜夜添| 国产成人freesex在线| 国产老妇伦熟女老妇高清| 久久久久精品久久久久真实原创| 可以在线观看毛片的网站| 色播亚洲综合网| 欧美精品人与动牲交sv欧美| 欧美成人午夜免费资源| 精品一区在线观看国产| 亚洲电影在线观看av| 青春草国产在线视频| 亚洲精品乱久久久久久| 岛国毛片在线播放| 亚洲精品第二区| 亚洲精品色激情综合| 一个人看视频在线观看www免费| 2021少妇久久久久久久久久久| 国产精品无大码| 黄色欧美视频在线观看| 久久韩国三级中文字幕| 国产亚洲5aaaaa淫片| 亚洲人成网站高清观看| 亚洲精品久久午夜乱码| 男女啪啪激烈高潮av片| 真实男女啪啪啪动态图| 午夜免费鲁丝| 91午夜精品亚洲一区二区三区| 亚洲久久久久久中文字幕| 听说在线观看完整版免费高清| 国产精品久久久久久久电影| 女的被弄到高潮叫床怎么办| 亚洲精品国产色婷婷电影| 亚洲欧美一区二区三区黑人 | 少妇熟女欧美另类| 亚洲自偷自拍三级| 国产毛片a区久久久久| 亚洲欧美一区二区三区国产| 免费av毛片视频| 亚洲国产日韩一区二区| 99热全是精品| 又粗又硬又长又爽又黄的视频| 麻豆精品久久久久久蜜桃|