• <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网| 亚洲国产成人一精品久久久| 国产v大片淫在线免费观看| 国产在视频线在精品| 欧美精品国产亚洲| 国产精品无大码| 国产午夜福利久久久久久| 日韩成人伦理影院| 久久久久久久国产电影| 我要搜黄色片| 日韩人妻高清精品专区| 久久亚洲国产成人精品v| 女人十人毛片免费观看3o分钟| 欧美又色又爽又黄视频| 国产精品久久久久久久久免| 永久免费av网站大全| a级毛片免费高清观看在线播放| 欧美一区二区精品小视频在线| 国产精品人妻久久久久久| 中文字幕人妻熟人妻熟丝袜美| 国产精品嫩草影院av在线观看| 久久久精品94久久精品| 国产高清国产精品国产三级 | 中文字幕人妻丝袜制服| 国产精品99久久99久久久不卡 | 免费大片18禁| 天天影视国产精品| 99久久中文字幕三级久久日本| kizo精华| 精品人妻在线不人妻| 国产欧美日韩一区二区三区在线| 亚洲av综合色区一区| 丰满迷人的少妇在线观看| 亚洲一区二区三区欧美精品| 我要看黄色一级片免费的| 欧美亚洲 丝袜 人妻 在线| 国产xxxxx性猛交| 免费黄色在线免费观看| 亚洲成国产人片在线观看| 久久久久精品人妻al黑| 亚洲第一av免费看| 青春草国产在线视频| av天堂久久9| 如日韩欧美国产精品一区二区三区| 成年av动漫网址| 女性被躁到高潮视频| 午夜老司机福利剧场| 成年人午夜在线观看视频| 老司机亚洲免费影院| 亚洲婷婷狠狠爱综合网| 婷婷色av中文字幕| 国产又色又爽无遮挡免| 久久99一区二区三区| 我的女老师完整版在线观看| 国产成人精品婷婷| 99香蕉大伊视频| 看十八女毛片水多多多| 久久久久久久大尺度免费视频| 侵犯人妻中文字幕一二三四区| 少妇熟女欧美另类| 春色校园在线视频观看| 中国国产av一级| 国产麻豆69| 久久99蜜桃精品久久| 午夜日本视频在线| 大话2 男鬼变身卡| 寂寞人妻少妇视频99o| 寂寞人妻少妇视频99o| 丁香六月天网| 成年人午夜在线观看视频| 男女高潮啪啪啪动态图| 26uuu在线亚洲综合色| 久久久久久久久久人人人人人人| 成人午夜精彩视频在线观看| 国产xxxxx性猛交| 久久精品久久久久久久性| 99久久综合免费| 国产精品欧美亚洲77777| 夫妻性生交免费视频一级片| videos熟女内射| 免费观看无遮挡的男女| 亚洲精品自拍成人| 国产成人欧美| 亚洲综合精品二区| av卡一久久| 国产亚洲最大av| 如日韩欧美国产精品一区二区三区| 国产成人精品久久久久久| 国产日韩欧美亚洲二区| 激情视频va一区二区三区| 亚洲,一卡二卡三卡| 午夜激情久久久久久久| 两性夫妻黄色片 | 欧美bdsm另类| 亚洲美女黄色视频免费看| 久久久久精品久久久久真实原创| 欧美精品一区二区大全| 在线看a的网站| 热99久久久久精品小说推荐| 亚洲国产色片| 高清欧美精品videossex| 天堂8中文在线网| 桃花免费在线播放| 国产黄频视频在线观看| 国产极品粉嫩免费观看在线| 国产 一区精品| 极品人妻少妇av视频| 亚洲国产精品一区二区三区在线| 涩涩av久久男人的天堂| 丁香六月天网| 欧美成人午夜精品| 又大又黄又爽视频免费| 久久久久精品人妻al黑| 免费黄频网站在线观看国产| 亚洲av在线观看美女高潮| 国产成人欧美| 久久久久久久国产电影| 欧美成人精品欧美一级黄| 精品亚洲成a人片在线观看| 丝袜喷水一区| 国产成人欧美| 一级黄片播放器| 精品99又大又爽又粗少妇毛片| 成年女人在线观看亚洲视频| 性高湖久久久久久久久免费观看| 满18在线观看网站| 观看av在线不卡| 久久午夜综合久久蜜桃| 免费日韩欧美在线观看| 18禁观看日本| 国产精品人妻久久久久久| 99视频精品全部免费 在线| 波多野结衣一区麻豆| 国产永久视频网站| 久久人妻熟女aⅴ| 亚洲欧美精品自产自拍| 大陆偷拍与自拍| 岛国毛片在线播放| 新久久久久国产一级毛片| 一区二区三区乱码不卡18| 一二三四中文在线观看免费高清| 亚洲国产av新网站| 国产极品天堂在线| 最新的欧美精品一区二区| 黑丝袜美女国产一区| 大香蕉久久网| 另类亚洲欧美激情| 免费少妇av软件| 久久久久久久亚洲中文字幕| 老熟女久久久| 日韩一区二区视频免费看| 在线观看一区二区三区激情| 少妇被粗大的猛进出69影院 | 热99国产精品久久久久久7| 中文字幕av电影在线播放| 啦啦啦视频在线资源免费观看| 国语对白做爰xxxⅹ性视频网站| 午夜激情久久久久久久| 国产亚洲欧美精品永久| 大香蕉久久成人网| 99国产综合亚洲精品| 熟女电影av网| 日本免费在线观看一区| 精品熟女少妇av免费看| 日本-黄色视频高清免费观看| 免费观看在线日韩| 日韩人妻精品一区2区三区| 免费女性裸体啪啪无遮挡网站| 欧美激情国产日韩精品一区| 精品少妇黑人巨大在线播放| 巨乳人妻的诱惑在线观看| 激情视频va一区二区三区| 久久久久网色| 亚洲av免费高清在线观看| 一区二区三区四区激情视频| 国产成人精品婷婷| 少妇人妻 视频| 欧美人与性动交α欧美精品济南到 | 香蕉国产在线看| www.色视频.com| 国产欧美另类精品又又久久亚洲欧美| 大片电影免费在线观看免费| 九色成人免费人妻av| 精品一区在线观看国产| 视频区图区小说| 黑人高潮一二区| 一本色道久久久久久精品综合| 国产精品欧美亚洲77777| 成人国语在线视频| 制服人妻中文乱码| av卡一久久| 丝袜美足系列| 午夜老司机福利剧场| 一级毛片我不卡| 丰满饥渴人妻一区二区三| 久久久a久久爽久久v久久| av视频免费观看在线观看| 精品少妇黑人巨大在线播放| 免费av不卡在线播放| 在线看a的网站| 日韩av在线免费看完整版不卡| 91国产中文字幕| 国产成人免费观看mmmm| 自拍欧美九色日韩亚洲蝌蚪91| 日韩三级伦理在线观看| 少妇的逼好多水| 啦啦啦视频在线资源免费观看| 亚洲欧美清纯卡通| 国产高清国产精品国产三级| av国产精品久久久久影院| 久久精品久久精品一区二区三区| 日本欧美国产在线视频| 亚洲人与动物交配视频| 欧美日韩亚洲高清精品| 老女人水多毛片| 午夜视频国产福利| 国产白丝娇喘喷水9色精品| 亚洲精品美女久久av网站| 久久久久久久久久成人| 人妻一区二区av| 国产免费现黄频在线看| 在线看a的网站| 美女主播在线视频| 十分钟在线观看高清视频www| 免费不卡的大黄色大毛片视频在线观看| 2018国产大陆天天弄谢| 人人妻人人添人人爽欧美一区卜| 18禁在线无遮挡免费观看视频| 亚洲成人手机| 我要看黄色一级片免费的| 美女福利国产在线| 亚洲色图综合在线观看| 亚洲国产精品一区三区| 伦理电影免费视频| 黄色视频在线播放观看不卡| 校园人妻丝袜中文字幕| 一本久久精品| 亚洲国产成人一精品久久久| 在线观看免费视频网站a站| 国产精品人妻久久久久久| 99国产精品免费福利视频| 天天操日日干夜夜撸| 免费在线观看黄色视频的| 日产精品乱码卡一卡2卡三| 看非洲黑人一级黄片| 日日撸夜夜添| 国产精品女同一区二区软件| 麻豆精品久久久久久蜜桃| 午夜视频国产福利| 精品福利永久在线观看| 99久久中文字幕三级久久日本| 亚洲国产精品一区三区| 欧美国产精品va在线观看不卡| 亚洲国产精品国产精品| 国产成人精品一,二区| 国产又爽黄色视频| 少妇的逼水好多| 全区人妻精品视频| 国产精品无大码| 精品第一国产精品| 亚洲五月色婷婷综合| 国产亚洲av片在线观看秒播厂| 91精品三级在线观看| 亚洲精品自拍成人| 亚洲av福利一区| 夜夜骑夜夜射夜夜干| 久久国产精品男人的天堂亚洲 | 久久人人爽av亚洲精品天堂| 丝袜在线中文字幕| 性高湖久久久久久久久免费观看| 99香蕉大伊视频| 制服人妻中文乱码| 亚洲高清免费不卡视频| 秋霞在线观看毛片| 国产又爽黄色视频| 纵有疾风起免费观看全集完整版| 久久精品国产自在天天线| 精品午夜福利在线看| 视频中文字幕在线观看| 9热在线视频观看99| 在线观看人妻少妇| 免费在线观看完整版高清| 亚洲经典国产精华液单| 亚洲精品视频女| 大片电影免费在线观看免费| 久久人人爽av亚洲精品天堂| 好男人视频免费观看在线| 亚洲av免费高清在线观看| 女性被躁到高潮视频| 久久精品国产亚洲av天美| 亚洲,欧美精品.| 天天操日日干夜夜撸| 色吧在线观看| 国产免费一区二区三区四区乱码| 美女视频免费永久观看网站| 80岁老熟妇乱子伦牲交| 亚洲一码二码三码区别大吗| 国产精品人妻久久久影院| 日韩av免费高清视频| 美女福利国产在线| 秋霞在线观看毛片| 一个人免费看片子| 老司机影院毛片| 五月伊人婷婷丁香| 成人漫画全彩无遮挡| 亚洲四区av| 97精品久久久久久久久久精品| 自线自在国产av| 考比视频在线观看| 日韩,欧美,国产一区二区三区| videosex国产| 国产在线免费精品| 亚洲国产日韩一区二区| 欧美激情国产日韩精品一区| 国产成人精品一,二区| 丁香六月天网| 国产日韩欧美在线精品| 国产av一区二区精品久久| 在线观看一区二区三区激情| 欧美人与善性xxx| 18+在线观看网站| 亚洲精品一二三| 两个人免费观看高清视频| 久久久久久久久久成人| 久久久久久久久久成人| 看免费成人av毛片| 久久国产精品大桥未久av| 美女脱内裤让男人舔精品视频| 日本黄大片高清| 久久久亚洲精品成人影院| 国产精品一国产av| 美女主播在线视频| av在线播放精品| 久久人人97超碰香蕉20202| 不卡视频在线观看欧美| 如日韩欧美国产精品一区二区三区| 免费久久久久久久精品成人欧美视频 | 国产精品一区二区在线观看99| 1024视频免费在线观看| 日韩av免费高清视频| 男人爽女人下面视频在线观看| 涩涩av久久男人的天堂| 国产精品一区二区在线观看99| 国产成人91sexporn| 人人澡人人妻人| 狂野欧美激情性bbbbbb| 午夜免费男女啪啪视频观看| freevideosex欧美| 久久影院123| 三级国产精品片| 免费黄频网站在线观看国产| 满18在线观看网站| 精品久久久久久电影网| 成人18禁高潮啪啪吃奶动态图| 国产精品一二三区在线看| kizo精华| 午夜福利在线观看免费完整高清在| 色吧在线观看| 久久97久久精品| 黄网站色视频无遮挡免费观看| 一区二区日韩欧美中文字幕 | 永久网站在线| 黑人高潮一二区| 久久国内精品自在自线图片| 乱码一卡2卡4卡精品| 久久久久久久久久久免费av| 婷婷色麻豆天堂久久| 国国产精品蜜臀av免费| 侵犯人妻中文字幕一二三四区| 欧美人与善性xxx| 国产一区二区三区av在线| 午夜免费鲁丝| 男女午夜视频在线观看 | av又黄又爽大尺度在线免费看| 性色avwww在线观看| 日日啪夜夜爽| 日本-黄色视频高清免费观看| 亚洲精品色激情综合| 国产黄色视频一区二区在线观看| 亚洲国产成人一精品久久久| av有码第一页| 69精品国产乱码久久久| 国产亚洲午夜精品一区二区久久| 高清欧美精品videossex| 高清在线视频一区二区三区| 日本黄大片高清| 国产男女内射视频| 九色亚洲精品在线播放| 十八禁网站网址无遮挡| 精品国产乱码久久久久久小说| 久久韩国三级中文字幕| 黑人巨大精品欧美一区二区蜜桃 | 丝袜人妻中文字幕| 欧美精品国产亚洲| 秋霞在线观看毛片| 在线精品无人区一区二区三| 日本色播在线视频| 国产综合精华液| 韩国高清视频一区二区三区| 看免费成人av毛片| 街头女战士在线观看网站| 黄片播放在线免费| 妹子高潮喷水视频| 深夜精品福利| 日本vs欧美在线观看视频| 好男人视频免费观看在线| 97在线人人人人妻| 国产在线视频一区二区| 精品国产乱码久久久久久小说| 国产老妇伦熟女老妇高清| 日韩熟女老妇一区二区性免费视频| 国产国拍精品亚洲av在线观看| 免费看av在线观看网站| 欧美xxⅹ黑人| 国产片内射在线| 麻豆精品久久久久久蜜桃| 精品国产一区二区三区久久久樱花| 国产综合精华液| 国产亚洲av片在线观看秒播厂| 男女下面插进去视频免费观看 | 国产欧美日韩一区二区三区在线| 99国产精品免费福利视频| av免费在线看不卡| 亚洲国产精品成人久久小说| 国产不卡av网站在线观看| 亚洲熟女精品中文字幕| 日日爽夜夜爽网站| 国产高清三级在线| 免费高清在线观看视频在线观看| 精品卡一卡二卡四卡免费| 国产又色又爽无遮挡免| 国产精品秋霞免费鲁丝片| 黄色毛片三级朝国网站| 国产精品国产av在线观看| 国产 一区精品| 五月天丁香电影| 人妻一区二区av| 日韩大片免费观看网站| 亚洲精品456在线播放app| 天天躁夜夜躁狠狠久久av| 精品国产国语对白av| 成人综合一区亚洲| av播播在线观看一区| 最后的刺客免费高清国语| 日韩一区二区视频免费看| 26uuu在线亚洲综合色| av不卡在线播放| 久久av网站| 国产一区亚洲一区在线观看| 国产精品久久久久久精品古装| 男女下面插进去视频免费观看 | 妹子高潮喷水视频| www.av在线官网国产| av一本久久久久| 国产精品麻豆人妻色哟哟久久| 制服诱惑二区| 曰老女人黄片| 久久久国产精品麻豆| 亚洲国产精品一区三区| 欧美3d第一页| 全区人妻精品视频| 日日摸夜夜添夜夜爱| 成人综合一区亚洲| 国产精品偷伦视频观看了| 热re99久久国产66热| 两个人免费观看高清视频| 亚洲av成人精品一二三区| 你懂的网址亚洲精品在线观看| a级片在线免费高清观看视频| 亚洲国产精品专区欧美| 日韩一区二区三区影片| 欧美最新免费一区二区三区| 国产男人的电影天堂91| 欧美日韩亚洲高清精品| 亚洲av国产av综合av卡| 中文乱码字字幕精品一区二区三区| 国产探花极品一区二区| 各种免费的搞黄视频| 女人被躁到高潮嗷嗷叫费观| 一级片免费观看大全| 黄色怎么调成土黄色| 成人综合一区亚洲| 美女主播在线视频| 国产片特级美女逼逼视频| 国产欧美亚洲国产| 欧美国产精品一级二级三级| 精品亚洲成a人片在线观看| 视频在线观看一区二区三区| 咕卡用的链子| 人妻一区二区av| 各种免费的搞黄视频| 飞空精品影院首页| 国产永久视频网站| 热re99久久国产66热| 亚洲成av片中文字幕在线观看 | 国产伦理片在线播放av一区| 搡女人真爽免费视频火全软件| 99热6这里只有精品| 人体艺术视频欧美日本| 青青草视频在线视频观看| 欧美变态另类bdsm刘玥| 黑人高潮一二区| 内地一区二区视频在线| 欧美日韩成人在线一区二区| 婷婷色综合www| 免费大片黄手机在线观看| 在线观看国产h片| 久久狼人影院| 亚洲欧洲国产日韩| 老司机影院成人| 免费av不卡在线播放| 国产综合精华液| 尾随美女入室| 国产精品欧美亚洲77777| 精品亚洲成国产av| 亚洲综合色惰| 欧美成人精品欧美一级黄| 大香蕉久久成人网| 一二三四在线观看免费中文在 | 亚洲国产av影院在线观看| 久久久久久久国产电影| 如日韩欧美国产精品一区二区三区| 久久午夜福利片| 熟女av电影| www.熟女人妻精品国产 | 女人精品久久久久毛片| 亚洲av综合色区一区| 这个男人来自地球电影免费观看 | 美女脱内裤让男人舔精品视频| 久久综合国产亚洲精品| 两个人免费观看高清视频| 一边亲一边摸免费视频| 高清毛片免费看| 欧美bdsm另类| 侵犯人妻中文字幕一二三四区| 国产成人精品久久久久久| 亚洲精品av麻豆狂野| av视频免费观看在线观看| av免费在线看不卡| 亚洲精品国产av蜜桃| 欧美人与性动交α欧美精品济南到 | 国产不卡av网站在线观看| 国产欧美另类精品又又久久亚洲欧美| 乱人伦中国视频| 午夜av观看不卡| 侵犯人妻中文字幕一二三四区| 久久人妻熟女aⅴ| 高清欧美精品videossex| 亚洲av电影在线进入| 午夜福利视频精品| 九草在线视频观看| 丰满迷人的少妇在线观看| 少妇的逼水好多| www.色视频.com| 亚洲欧美中文字幕日韩二区| 久久久久久久久久人人人人人人| 制服人妻中文乱码| 久久久欧美国产精品| 多毛熟女@视频| 女的被弄到高潮叫床怎么办| 国产精品久久久久久久电影| 18禁在线无遮挡免费观看视频| 久久狼人影院| 国产熟女午夜一区二区三区| 亚洲五月色婷婷综合| 男女免费视频国产| 国产精品嫩草影院av在线观看| 国产男女超爽视频在线观看| 美女内射精品一级片tv| 国产精品.久久久| 少妇猛男粗大的猛烈进出视频| 伦精品一区二区三区| 综合色丁香网| 宅男免费午夜| 成人影院久久| 伦精品一区二区三区| 国产高清国产精品国产三级| 97超碰精品成人国产| 一级a做视频免费观看| 少妇的逼水好多| 亚洲欧洲精品一区二区精品久久久 | 美女国产高潮福利片在线看| 免费女性裸体啪啪无遮挡网站| 一本久久精品| 精品熟女少妇av免费看| 国产亚洲精品久久久com| 精品国产一区二区三区久久久樱花| 一本大道久久a久久精品| 免费高清在线观看日韩| 在线观看免费高清a一片| 免费大片黄手机在线观看| 久久国产亚洲av麻豆专区| 看免费成人av毛片| 久久99热这里只频精品6学生| 人妻一区二区av| 黑人欧美特级aaaaaa片| 亚洲成av片中文字幕在线观看 | 高清黄色对白视频在线免费看| 中文字幕精品免费在线观看视频 | 日本爱情动作片www.在线观看| 色5月婷婷丁香| 校园人妻丝袜中文字幕| 日韩,欧美,国产一区二区三区| 久久 成人 亚洲| 大香蕉久久网| 亚洲欧美日韩卡通动漫| 性高湖久久久久久久久免费观看| videossex国产| 免费女性裸体啪啪无遮挡网站| 狂野欧美激情性bbbbbb| 日韩中字成人| 久久久久久久亚洲中文字幕| 亚洲av电影在线观看一区二区三区| 九草在线视频观看|