亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        DQFEM求解正交各向異性薄板單元靜力問題

        2020-07-09 03:28:40鄧小環(huán)許華南陳水梅黃清云
        武夷學院學報 2020年3期
        關鍵詞:薄板微分結點

        鄧小環(huán),許華南,陳水梅,黃清云

        (龍巖學院 資源工程學院,福建 龍巖 364000)

        1 理論基礎

        通過采用DQFEM(微分求積有限元法)[1],分析正交各向異性薄板單元[2],并將所求得與有限元軟件Abaqus所得結果進行比較,進一步驗證微分求積有限元法的精確性。

        本文采用薄板理論[3],即不考慮橫向剪切變形。在彈性力學中,將兩個平行的面和垂直于這兩個平行面的棱柱或柱體所圍成的物體稱為板,如圖1所示即為矩形薄板示意圖。

        圖1 矩形薄板示意圖Fig.1 Schematic diagram of a rectangular plate

        利用薄板理論在分析平板彎曲時假設[4],將平板彎曲問題簡化為二維問題,且全部應力和應變可以用板中面的撓度w表示,即:

        2 DQFEM求解正交各向異性材料薄板問題

        2.1 二維函數(shù)的微分求解法則

        對于二維函數(shù)f(x,y),通??捎脙蓚€一維函數(shù)的積表示[5],即:

        經典的微分求積方法通常適合用于矩形區(qū)域[6]。根據(jù)微分求積法則,一維函數(shù)p(x)在結點xi處對坐標x的r階導數(shù)和q(y)在結點yi處對坐標y的s階導數(shù)可由結點函數(shù)值分別表示為:

        為方便推導二維微分求積有限單元矩陣,引入以下矩陣和向量:

        由下面的遞推關系可得到權系數(shù)為:

        2.2 DQFEM求解薄板單元

        對于線彈性問題,總勢能泛函包括應變能和外力勢[7],即:

        式中:ε和D分別是應變向量和材料矩陣;u是位移向量;q是荷載向量。

        結構的動能泛函為:

        對于薄板彎曲問題,撓度函數(shù)也可用拉格朗日函數(shù)表示[8]為:

        薄板的本構關系:

        因為εz=0,γzx=0,γzy=0,可得:

        其應變-位移關系[9]為:

        定義如下向量:

        應用微分求積法則,可得結點應變向量為:

        將式(25)和式(22)帶入式(17)和式(18),離散薄板勢能泛函為:

        C=diag(Ck),其中分別為x和y方向的高斯-洛巴托積分系數(shù),

        q為結點荷載列向量,其形式同式(24)。

        為了滿足薄板單元之間的C1連續(xù)條件,施加固支邊界條件需要對位移向量做一些修正。邊界結點參數(shù)應該既包括結點撓度,也包括轉角,于是有:

        式中:wkx=(?w/?x)k,wky=(?w/?y)k,wkxy=(?2w/?x?y)k。

        為方便書寫轉換矩陣T,此處設板x方向結點數(shù)M=2,y方向結點數(shù)N=2,則:

        則轉換矩陣T為:

        將式(30)帶入式(28),并進行變分,得到與結點位移向量w相對應的微分求積有限薄板單元的質量矩陣、剛度矩陣和荷載列向量,即:

        綜上可得:

        2.3 算例

        薄板如圖1所示,板長寬均為1,板厚h=0.01,板密度為1,彈性模量E=1,泊松比v=0.3,x方向結點數(shù)為M,y方向結點數(shù)為N。薄板對側兩邊界約束條件為簡支約束,在結點上施加均勻力為1N的力,計算板固定結點的位移值,并與有限元軟件的計算結果進行對比。如下分別為網格密度為10×10,20×20,30×30,40×40,50×50,100×100薄板示意圖:

        圖2 薄板有限元軟件網格密度示意圖Fig.2 Grid density diagram of thin plate finite element software

        通過DQFEM對上述4種網格密度進行編程分析,取薄板相同位置的結點位移值,得到如下表1至表3所表示的位移值。

        表1 有限元軟件與DQFEM計算薄板結點位移值及誤差對比Tab.1 Comparison of finite element software and DQFEM for calculating plate node displacement values and errors

        表2 DQFEM計算薄板結點位移值及誤差對比Tab.2 The calculation of displacement value and error comparison of thin plate node by DQFEM

        表3 有限元軟件計算薄板結點位移值及誤差對比Tab.3 The displacement value and error comparison of thin plate node calculated by finite element software

        由表2可得采用DQFEM計算得到的薄板結點位移值,網格大小對于位移值取值影響不大,即該方法可采用較少的結點數(shù)解決問題。由表3可知采用有限元軟件得到的位移值隨網格大小的不同而相差較大,表1則說明有限元軟件需要更密的網格才能接近采用DQFEM計算得到的位移值,進一步說明DQFEM方法的高效性及精確性。

        3 小結

        對于正交各向異性薄板單元,采用微分求積有限元法離散最小勢能原理的泛函,通過相應的Matlab程序得到所求問題的剛度矩陣,并最終得到正交各向異性薄板單元在受壓工況下的位移值,并采用有限元軟件ABAQUS中各向異性殼單元模擬了相應的工況,數(shù)值模擬結果和上述所求計算結果進行系統(tǒng)對比分析,進一步驗證了微分求積有限元法的精確性。

        猜你喜歡
        薄板微分結點
        一角點支撐另一對邊固支正交各向異性矩形薄板彎曲的辛疊加解
        擬微分算子在Hp(ω)上的有界性
        10MN鋁合金薄板拉伸機組的研制
        重型機械(2020年2期)2020-07-24 08:16:12
        上下解反向的脈沖微分包含解的存在性
        Ladyzhenskaya流體力學方程組的確定模與確定結點個數(shù)估計
        借助微分探求連續(xù)函數(shù)的極值點
        鋁薄板高速DP-GMAW焊接性能的研究
        焊接(2016年5期)2016-02-27 13:04:42
        基于PLC的薄板激光焊機控制系統(tǒng)
        焊接(2015年10期)2015-07-18 11:04:46
        對不定積分湊微分解法的再認識
        基于Raspberry PI為結點的天氣云測量網絡實現(xiàn)
        亚洲毛片网| 亚洲av无码乱码国产一区二区| 欧美最大胆的西西人体44| 人妻夜夜爽天天爽三区麻豆AV网站| 久久精品国产久精国产69| 日本老熟妇五十路一区二区三区 | 国产精品三级在线专区1| 国产视频导航| 91色婷婷成人精品亚洲| 东京热东京道日韩av| 97cp在线视频免费观看| 国产精品毛片完整版视频| 在线免费观看国产精品| 亚洲av黄片一区二区| 亚洲中文字幕第一页免费| 无码a级毛片免费视频内谢5j| 先锋影音最新色资源站| 亚洲熟妇少妇任你躁在线观看| 中文字幕无码免费久久| 玖玖资源站亚洲最大的网站| 无码字幕av一区二区三区| 国产成人综合久久精品推| 久久久久久无中无码| 91国产精品自拍视频| 亚洲av中文无码乱人伦下载| 激情久久av一区av二区av三区| 国产精品电影久久久久电影网| 国产激情免费观看视频| 国产精品亚洲专区无码不卡| 国产成人精品久久一区二区三区| 国产免费资源高清小视频在线观看 | 亚洲精品国产美女久久久| 69av视频在线观看| 精品国产车一区二区三区| 国产大屁股视频免费区| 国产精品伦一区二区三级视频| 伊人久久综在合线亚洲不卡| 国产一级做a爱视频在线| 性感女教师在线免费观看| www插插插无码视频网站| 麻豆AⅤ精品无码一区二区|