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

    溫度影響下Winkler-Pasternak彈性地基上多孔FGM矩形板的自由振動(dòng)分析

    2022-04-09 01:54:40滕兆春王偉斌鄭文達(dá)
    工程力學(xué) 2022年4期
    關(guān)鍵詞:無量邊界條件固有頻率

    滕兆春,王偉斌,鄭文達(dá)

    (蘭州理工大學(xué)理學(xué)院,蘭州 730050)

    功能梯度材料(functionally graded materials,F(xiàn)GM)具有材料性能平滑、連續(xù)變化的特點(diǎn),在世界范圍內(nèi)得到了迅速的研究、開發(fā)和應(yīng)用[1],F(xiàn)GM矩形板被廣泛應(yīng)用于土木、機(jī)械和海洋工程等[2]。隨著新型材料的興起,F(xiàn)GM板結(jié)構(gòu)具有常見復(fù)合材料板難以比擬的優(yōu)點(diǎn)。關(guān)于FGM常見結(jié)構(gòu)的力學(xué)行為研究已有較多文獻(xiàn)成果,如Wattanasakulpong等[3]采用微分變換法(differential transformation method,DTM)研究了FGM階梯梁的自由振動(dòng),得到了梁的固有頻率和振型。Roshan等[4]在經(jīng)典平板理論的基礎(chǔ)上,對(duì)面內(nèi)均布力作用下的FGM圓板振動(dòng)問題進(jìn)行了分析和數(shù)值計(jì)算,研究了梯度指數(shù)和面內(nèi)力對(duì)前三階振動(dòng)固有頻率的影響。Ghazaryan等[5]利用DTM對(duì)非均勻截面軸向功能梯度Euler-Bernoulli梁的自由振動(dòng)進(jìn)行了分析,并驗(yàn)證了該方法可用于軸向材料非均勻或者幾何非均勻的各種類型梁的計(jì)算。Shishesaz等[6]基于應(yīng)變梯度理論,研究了FGM納米圓盤的熱彈性行為,主要分析了溫度變化對(duì)納米圓盤應(yīng)力和徑向位移的影響。Abdolhossein等[7]考慮流體-結(jié)構(gòu)相互作用研究了彈性地基上部分充液FGM圓柱殼在熱環(huán)境下的自由振動(dòng)問題。Ladislav等[8]基于經(jīng)典彈性理論研究了FGM薄板在熱載荷作用下的彎曲問題,結(jié)合移動(dòng)最小二乘逼近的數(shù)值方法,重點(diǎn)研究了材料相關(guān)參數(shù)對(duì)板撓度的影響。Saini等[9]基于經(jīng)典板理論,考慮材料溫度依賴性質(zhì),分析了熱環(huán)境中FGM圓板在均勻面內(nèi)軸向載荷和沿厚度方向的非線性溫度分布作用下軸對(duì)稱屈曲和振動(dòng)問題。Chakraverty等[10]基于基爾霍夫板理論,利用Rayleigh-Ritz法研究了FGM矩形板的自由振動(dòng)問題。蔣偉男等[11]基于三階剪切理論提出考慮橫向拉伸影響的位移場(chǎng),研究了具有簡(jiǎn)支和固支邊界條件下的FGM矩形夾層板的自由振動(dòng)。何昊南等[12]基于幾何非線性和物理中面,分析了熱對(duì)材料參數(shù)影響下FGM梁的熱后屈曲特性。當(dāng)然,上述研究工作大多是從材料的理想形態(tài)入手建立力學(xué)模型并進(jìn)行分析,忽略了FGM在實(shí)際應(yīng)用中孔隙的存在。

    FGM在各領(lǐng)域的廣泛需求,促進(jìn)了其制備技術(shù)的快速發(fā)展,各種FGM制備方法先后出現(xiàn)。但由于現(xiàn)有某些制備技術(shù)的缺陷,制備過程中在FGM內(nèi)部不可避免的會(huì)出現(xiàn)小的孔隙,則形成多孔FGM。同時(shí),隨著FGM越來越多的應(yīng)用擴(kuò)展,也可以使用陶瓷和泡沫金屬來制備多孔FGM形成工程中可應(yīng)用的一些特殊功能結(jié)構(gòu),如過濾器、催化器、燃燒室和生物材料等。這些孔隙對(duì)FGM的物理性能、功能使用和靜動(dòng)力特性有著較大的影響,因此越來越多的科研人員開始關(guān)注多孔FGM結(jié)構(gòu)的力學(xué)行為研究。Behravan等[13]基于微分求積和狀態(tài)空間矢量技術(shù)的求解方法,對(duì)彈性地基非均勻面力作用下非對(duì)稱變厚度多孔FGM圓板的三維磁彈性特性進(jìn)行了分析。Aria等[14]首次提出了一種有限元模型,基于冪律形式分析了多孔FGM納米梁的熱彈性行為。蒲育等[15?17]基于一種擴(kuò)展的廣義梁理論,分別采用廣義微分求積法和萊維法,考慮材料性質(zhì)隨溫度的相關(guān)性,研究了熱環(huán)境中軸向機(jī)械力作用含均勻孔隙FGM梁的耦合振動(dòng)和耦合屈曲特性以及濕-熱-機(jī)-彈耦合作用下多孔FGM梁的振動(dòng)和屈曲特性。蘇盛開等[18]采用Euler梁理論和高階三角剪切變形理論研究了多孔FGM梁的熱力耦合屈曲行為。Ebrahimi等[19]利用DTM對(duì)轉(zhuǎn)動(dòng)多孔FGM Timoshenko梁進(jìn)行橫向振動(dòng)分析,并在修正混合率公式的基礎(chǔ)上研究了孔隙率、梯度指數(shù)、轉(zhuǎn)速等參數(shù)對(duì)轉(zhuǎn)動(dòng)多孔FGM梁無量綱固有頻率的影響。綜上所述,已有文獻(xiàn)中大多學(xué)者基本都專注于多孔FGM梁的研究,對(duì)多孔FGM板的研究則相對(duì)較少,同時(shí)一些研究結(jié)果也很好地驗(yàn)證了DTM的實(shí)用性和有效性。

    目前關(guān)于溫度影響下Winkler-Pasternak彈性地基上多孔FGM矩形板自由振動(dòng)問題的研究在國(guó)內(nèi)外還未見有文獻(xiàn)報(bào)道,且考慮到彈性地基板在工程中的廣泛應(yīng)用背景,因此在參照文獻(xiàn)[20]的基礎(chǔ)上,本文將一種隨機(jī)均勻分布的孔隙引入力學(xué)模型,采用DTM對(duì)溫度影響下Winkler- Pasternak彈性地基上多孔FGM矩形板自由振動(dòng)展開研究,分析梯度指數(shù)、孔隙率、升溫、地基剛度系數(shù)、邊界條件、長(zhǎng)寬比和邊厚比等對(duì)多孔FGM矩形板無量綱固有頻率的影響,為今后多孔FGM板的研究提供數(shù)據(jù)支撐。

    1 控制微分方程的建立

    1.1 材料物性描述

    如圖1所示:在三維笛卡爾坐標(biāo)系中,溫度影響下Winkler-Pasternak彈性地基上多孔FGM矩形板厚度為h, 長(zhǎng)為a, 寬為b,xoy平面為板的幾何中面,kw和qw分別為Winkler-Pasternak地基的彈性剛度系數(shù)和剪切剛度系數(shù)。

    圖1 彈性地基上多孔FGM矩形板的幾何模型和截面上均勻分布的孔隙Fig.1 Geometric model of a porous FGM rectangular plate on elastic foundation and cross-section with evenly distributed porosity

    考慮到 FGM 的物性參數(shù)與溫度T有關(guān),采用Voigt混合冪律模型[21]。物性參數(shù)滿足:

    式中:P可表示彈性模量E、密度 ρ和熱膨脹系數(shù)α ;T為絕對(duì)溫度;n為 梯度指數(shù);下標(biāo) c 和 m則分別為陶瓷和金屬材料的物性參數(shù)。

    上述金屬和陶瓷的某一物性參數(shù)G隨溫度T的變化可采用 Touloukian 非線性函數(shù),統(tǒng)一表述為[14]:

    式 中,G0、G-1、G1、G2和G3為 溫 度 相 關(guān) 的 系數(shù),一般由實(shí)驗(yàn)給出,可見表1。

    表1 陶瓷和金屬材料物性參數(shù)的溫度相關(guān)系數(shù)表[14]Table 1 Temperature correlation coefficient of physical parameters of ceramic and metal materials[14]

    當(dāng)FGM中存在孔隙時(shí),可以假設(shè)孔隙的大小和分布是任意的,并且介質(zhì)是連續(xù)致密的,在此基礎(chǔ)上得到了孔隙率 θ與彈性模量E、泊松比ν和密度 ρ的關(guān)系[20]。定義孔隙率如下:

    式中: ρ(z) 為 致密密度; ρ0(z)為自然密度。

    有效泊松比和有效彈性模量分別為:

    式中:ν、E(z) 和ν?、E?(z) 分 別為 θ=0 和 0<θ<1時(shí)的泊松比及彈性模量;bc=2-3ν。將式(1)~式(3)代入式(5)即可得到不同孔隙率下 FGM 矩形板的有效彈性模量。

    1.2 Winkler-Pasternak地基模型和升溫

    雙參數(shù)Winkler-Pasternak彈性地基中考慮了剪力的作用,地基-板相互作用力F和矩形板橫向位移w的關(guān)系為[22]:

    式中:kw、qw分別為地基的彈性剛度系數(shù)和剪切剛度系數(shù); ?2為二維Laplace算子。

    本文僅考慮均勻升溫:

    式中:T0為無應(yīng)力狀態(tài)的參考溫度,本文取T0=300 K; ΔT為升溫值。

    1.3 基本方程

    根據(jù)薄板和小撓度理論的基本假設(shè),溫度影響下Winkler-Pasternak彈性地基上多孔FGM矩形板位移場(chǎng)為:

    式中:u、v和w分 別為板內(nèi)一點(diǎn)沿x、y和z方向的位移分量;z0為物理中面。物理中面內(nèi)應(yīng)力分量與應(yīng)變分量均為零,其表達(dá)式如下[23]:

    幾何方程為:

    本構(gòu)方程為:

    下面應(yīng)用Hamilton原理建立溫度影響下Winkler-Pasternak彈性地基上多孔FGM矩形板的運(yùn)動(dòng)微分方程。已知哈密頓原理如下:

    式中: Π 、U和V分別為地基-板系統(tǒng)的動(dòng)能、應(yīng)變能和外力勢(shì)能;δ為變分符號(hào);t1和t2分別為系統(tǒng)運(yùn)動(dòng)的初始時(shí)刻和終止時(shí)刻。這里只考慮橫向位移,令w=w(x,y,t),可得:

    考慮溫度在板內(nèi)均勻變化,且忽略面內(nèi)和轉(zhuǎn)動(dòng)慣性力,可得溫度影響下彈性地基上多孔FGM矩形板橫向運(yùn)動(dòng)的控制微分方程:

    2 無量綱控制微分方程及DTM變換

    2.1 無量綱控制微分方程

    溫度影響下彈性地基上多孔FGM矩形板在y=0 和y=b處 為簡(jiǎn)支(S)邊界,而在x=0 和x=a處可為固定(C)邊界、簡(jiǎn)支(S)邊界或自由(F)邊界。在本文的問題討論中,邊界條件按照x=0、y=b、x=a和y=0的順序給出。

    綜上可得到溫度影響下Winkler-Pasternak彈性地基上多孔FGM矩形板自由振動(dòng)時(shí)的無量綱控制微分方程:

    式中:

    式中:A1、A2、A3分別為四階微分方程各項(xiàng)系數(shù);η為慣性系數(shù)。

    接下來考慮多孔 FGM 矩形板在X=0 和X=1處的無量綱邊界條件,其形式如下:

    簡(jiǎn)支(S):

    固定(C):

    自由(F):

    2.2 DTM變換

    對(duì)于溫度影響下Winkler-Pasternak彈性地基上多孔FGM矩形板自由振動(dòng)的無量綱控制微分方程式(17)的求解,本文采用一種半解析法—微分變換法(DTM)[24]進(jìn)行求解。按照DTM的求解過程及原理,將溫度影響下彈性地基上多孔FGM矩形板自由振動(dòng)的無量綱控制微分方程轉(zhuǎn)換為如下的迭代代數(shù)方程:

    邊界條件的DTM變換如下:

    在X=0處:

    簡(jiǎn)支(S):

    固定(C):

    自由(F):

    在X=1處:

    簡(jiǎn)支(S):

    固定(C):

    自由(F):

    結(jié)合式(21)~式(27),邊界條件為兩邊固定兩邊簡(jiǎn)支(CSCS)、一邊固定三邊簡(jiǎn)支(CSSS)和一邊固支一邊自由兩對(duì)邊簡(jiǎn)支(CSFS)的關(guān)于無量綱固有頻率的特征方程如下:

    對(duì)于邊界條件為四邊簡(jiǎn)支(SSSS)和一邊自由三邊簡(jiǎn)支(SSFS),關(guān)于無量綱固有頻率的特征方程如下:

    要求得無量綱的固有頻率必須滿足:

    對(duì)于邊界條件為兩邊自由兩邊簡(jiǎn)支(FSFS),關(guān)于無量綱固有頻率的特征方程如下:

    同理,必須滿足:

    式中,γ為迭代誤差限,這里取 γ=0.000001。

    3 計(jì)算結(jié)果及分析

    應(yīng)用MATLAB軟件編寫相關(guān)程序,由此可得到由DTM求解溫度影響下Winkler-Pasternak彈性地基上多孔FGM矩形板自由振動(dòng)時(shí)的無量綱固有頻率。為了驗(yàn)證DTM的正確性,將溫度影響下彈性地基上多孔FGM矩形板,經(jīng)算例1退化為孔隙率為零時(shí)的FGM矩形板和Winkler彈性地基上FGM矩形板,并與文獻(xiàn)[25]中的結(jié)果進(jìn)行比較。經(jīng)算例2得到新的結(jié)果。

    算例1. 300 K下陶瓷和金屬的物性參數(shù)分別取為: ρc=3800kg/m3, ρm=2700kg/m3,Ec=380 GPa,Em=70 GPa。

    表2 CSSS和CSFS邊界條件下不同長(zhǎng)寬比、不同梯度指數(shù)對(duì)FGM矩形板無量綱固有頻率的影響Table 2 Influence of aspect ratio and gradient index on dimensionless natural frequencies of FGM rectangular plate with CSSS and CSFS boundary conditions

    表3 CSCS和SSSS邊界條件下長(zhǎng)寬比對(duì)Winkler彈性地基上FGM矩形板無量綱固有頻率的影響(n=1, K=100)Table 3 Influence of aspect ratio on dimensionless natural frequencies of FGM rectangular plate on Winkler elastic foundation with CSCS and SSSS boundary conditions (n = 1, K = 100)

    算例2. 圖2~圖6中各物性參數(shù)均遵循表1中的實(shí)驗(yàn)數(shù)據(jù),其中,泊松比ν取陶瓷材料和金屬材料的平均值即:ν=(νc+νm)/2=0.280。

    圖2反映了當(dāng)θ =0.2 ,H=0.1 ,λ=1 ,K=100,Q=10 ,ΔT=100K,在CSCS、CSSS、SSSS、CSFS、SSFS和FSFS邊界條件下梯度指數(shù)n與前三階的無量綱固有頻率 Ω的關(guān)系曲線。結(jié)果顯示:隨著n的增大,無量綱固有頻率 Ω總體上在減小。當(dāng)n在小值范圍取值時(shí),無量綱固有頻率減小趨勢(shì)很劇烈;當(dāng)n在較大值范圍取值時(shí),無量綱固有頻率變化趨于平緩;當(dāng)n的值一定時(shí),邊界約束越強(qiáng),無量綱固有頻率越大,即CSCS固有頻率>CSSS固有頻率>SSSS固有頻率>CSFS固有頻率>SSFS固有頻率>FSFS固有頻率。

    圖2 不同邊界條件下梯度指數(shù)對(duì)Winkler-Pasternak彈性地基上多孔FGM矩形板無量綱固有頻率的影響Fig.2 Influence of gradient index on the dimensionless natural frequencies of porous FGM rectangular plates on Winkler-Pasternak elastic foundations under different boundary conditions

    圖3給出當(dāng)n=1 ,H=0.1 ,λ=1 ,K=100,Q=10 , ΔT=100 K時(shí),在CSCS、CSSS、SSSS、CSFS、SSFS和FSFS邊界條件下孔隙率 θ與前三階無量綱固有頻率 Ω的關(guān)系。結(jié)果顯示:一階無量綱固有頻率在CSCS、CSSS、SSSS邊界下隨著孔隙率的增大而減小,在CSFS、SSFS和FSFS邊界下,一階無量綱固有頻率起初隨孔隙率增大而減小,后又微弱增大趨勢(shì)。而二階、三階無量綱固有頻率隨著孔隙率的增大在減小,且減小的幅度比一階要大。當(dāng) θ一定時(shí),邊界約束越強(qiáng),無量綱固有頻率越大,即CSCS固有頻率>CSSS固有頻率>SSSS固有頻率>CSFS固有頻率>SSFS固有頻率>FSFS固有頻率。

    圖3 不同邊界條件下孔隙率對(duì)Winkler-Pasternak彈性地基上多孔FGM矩形板無量綱固有頻率的影響Fig.3 Influence of porosity on the dimensionless natural frequency of porous FGM rectangular plate on Winkler-Pasternak elastic foundation under different boundary conditions

    圖4給出當(dāng)n=1 ,H=0.1 ,λ=1 ,K=100,θ=0.2 ,Q=10時(shí),在CSCS、CSSS、SSSS、CSFS、SSFS和FSFS邊界條件下溫升ΔT與前三階無量綱固有頻率 Ω的關(guān)系。結(jié)果顯示:前三階無量綱固有頻率隨著溫升的不斷增大而減小,且二階和三階無量綱固有頻率減小的幅度比一階大。

    圖4 不同邊界條件下溫升對(duì)Winkler-Pasternak彈性地基上多孔FGM矩形板無量綱固有頻率的影響Fig.4 Influence of temperature rise on the dimensionless natural frequency of porous FGM rectangular plate on Winkler-Pasternak elastic foundation under different boundary conditions

    圖5為n=1 ,H=0.1 ,ΔT=100 K,K=100,θ=0.2 ,Q=10時(shí),在CSCS、CSSS、SSSS、CSFS、SSFS和FSFS邊界條件下長(zhǎng)寬比λ與前兩階無量綱固有頻率 Ω的關(guān)系。結(jié)果顯示:前兩階無量綱固有頻率隨著長(zhǎng)寬比的增大而增大,且二階無量綱固有頻率增大的幅度比一階大。當(dāng)長(zhǎng)寬比λ一定時(shí),邊界約束越強(qiáng),無量綱固有頻率越大,即CSCS固有頻率>CSSS固有頻率>SSSS固有頻率>CSFS固有頻率>SSFS固有頻率>FSFS固有頻率。

    圖5 不同邊界條件下長(zhǎng)寬比對(duì)Winkler-Pasternak彈性地基上多孔FGM矩形板無量綱固有頻率的影響Fig.5 Influence of the aspect ratio on the dimensionless natural frequency of the porous FGM rectangular plate on the Winkler-Pasternak elastic foundation under different boundary conditions

    圖6 給出了在CSCS、CSSS、SSSS、CSFS、SSFS和FSFS邊界條件下,當(dāng)n=1 ,H=0.1 , ΔT=100 K,λ=1 , θ=0.2 ,Q=10無量綱彈性剛度系數(shù)K對(duì)無量綱固有頻率 Ω的影響。結(jié)果表明:前兩階無量綱固有頻率隨著無量綱剛度系數(shù)的增大而增大,且二階無量綱固有頻率幾乎呈線性變化,當(dāng)無量綱彈性剛度系數(shù)K一定時(shí),邊界約束越強(qiáng),無量綱固有頻率越大, 即CSCS固有頻率>CSSS固有頻率>SSSS固有頻率>CSFS固有頻率> SSFS固有頻率>FSFS固有頻率。

    圖6 不同邊界條件下無量綱彈性剛度系數(shù)對(duì)Winkler-Pasternak彈性地基上多孔FGM矩形板無量綱固有頻率的影響Fig.6 Influence of dimensionless elastic stiffness coefficient on the dimensionless natural frequencies of porous FGM rectangular plate on Winkler-Pasternak elastic foundation with different boundary conditions

    圖7 反映了在CSCS、CSSS、SSSS、CSFS、SSFS和FSFS邊界條件下,當(dāng)n=1 ,H=0.1 , ΔT=100 K,λ=1 , θ=0.2 ,K=100,無量綱剪切剛度系數(shù)Q對(duì)無量綱固有頻率 Ω的影響,結(jié)果表明:前兩階無量綱固有頻率隨著無量綱剛度系數(shù)Q的增大而增大。且二階無量綱固有頻率比一階無量綱固有頻率增大得到幅度大,當(dāng)Q的值一定時(shí),邊界約束越強(qiáng),無量綱固有頻率越大,即CSCS固有頻率>CSSS固有頻率>SSSS固有頻率>CSFS固有頻率> SSFS固有頻率>FSFS固有頻率。

    圖7 不同邊界條件下無量綱剪切剛度系數(shù)對(duì)Winkler-Pasternak彈性地基上多孔FGM矩形板無量綱固有頻率的影響無量綱固有頻率的影響Fig.7 Influence of dimensionless shear stiffness coefficient on the dimensionless natural frequencies of porous FGM rectangular plate on Winkler-Pasternak elastic foundation with different boundary conditions

    圖8 給出了在CSCS、CSSS、SSSS、CSFS、SSFS和FSFS邊界條件下,當(dāng)n=1 , ΔT=100K ,λ=1 , θ=0.2 ,K=100 ,Q=10 , 邊 厚 比 1/H對(duì)一階無量綱固有頻率 Ω1(即無量綱基頻)的影響。這里邊厚比的取值為8~60之間,結(jié)果表明:一階無量綱固有頻率均隨著邊厚比 1/H的增大而減小。由于隨著邊厚比的增大,多孔FGM矩形板的整體剛度趨向于減小,從而導(dǎo)致這一必然結(jié)果。

    圖8 不同邊界條件下邊厚比對(duì)Winkler-Pasternak彈性地基上多孔FGM矩形板一階無量綱固有頻率的影響Fig.8 Influence of side-to-thickness ratio on first dimensionless natural frequencies of porous FGM rectangular plate on Winkler-Pasternak elastic foundation with different boundary conditions

    4 結(jié)論

    本文基于Hamilton原理和經(jīng)典薄板理論建立了溫度影響下Winkler-Pasternak彈性地基上多孔FGM矩形板自由振動(dòng)的控制微分方程并無量綱化,然后運(yùn)用微分變換法(DTM)對(duì)無量綱控制微分方程及其邊界條件進(jìn)行變換,通過編寫MATLAB程序求解出無量綱固有頻率,并分別給出了不同參數(shù)值的變化對(duì)無量綱固有頻率的影響。得到主要結(jié)論如下:

    (1) 多孔FGM矩形板無量綱固有頻率 Ω 隨著梯度指數(shù)n的增大而減小,隨著梯度指數(shù)n趨向于無窮大,無量綱固有頻率趨于不變,合理的解釋了FGM中陶瓷材料向金屬材料過渡的特性。

    (2) 多孔FGM矩形板無量綱固有頻率 Ω隨著孔隙率 θ變化:一階無量綱固有頻率在CSCS、CSSS、SSSS邊界下隨著孔隙率的增大而減小,在CSFS、SSFS和FSFS邊界下,一階無量綱固有頻率起初隨孔隙率增大而減小,后有微弱增大趨勢(shì)。而二階、三階無量綱固有頻率隨著孔隙率的增大在減小,且減小的幅度比一階要大。

    (3) 多孔FGM矩形板前三階無量綱固有頻率Ω隨著溫升ΔT的增大而減小,且二階和三階無量綱固有頻率減小的幅度比一階大。

    (4) 多孔FGM矩形板無量綱固有頻率 Ω隨著無量綱彈性剛度系數(shù)K的增大而增大,且二階無量綱固有頻率幾乎呈線性變化。隨著無量綱剪切剛度系數(shù)Q的增大而增大。

    (5) 多孔FGM矩形板前兩階無量綱固有頻率隨著長(zhǎng)寬比的增大而增大。一階無量綱固有頻率Ω隨著邊厚比 1/H的增大而減小。

    猜你喜歡
    無量邊界條件固有頻率
    烏雷:無量之物
    現(xiàn)場(chǎng)測(cè)定大型水輪發(fā)電機(jī)組軸系的固有頻率
    一類帶有Stieltjes積分邊界條件的分?jǐn)?shù)階微分方程邊值問題正解
    帶有積分邊界條件的奇異攝動(dòng)邊值問題的漸近解
    劉少白
    藝術(shù)品(2020年8期)2020-10-29 02:50:02
    論書絕句·評(píng)謝無量(1884—1964)
    炳靈寺第70 窟無量壽經(jīng)變辨識(shí)
    西藏研究(2017年3期)2017-09-05 09:45:07
    帶Robin邊界條件的2維隨機(jī)Ginzburg-Landau方程的吸引子
    總溫總壓測(cè)頭模態(tài)振型變化規(guī)律研究
    A novel functional electrical stimulation-control system for restoring motor function of post-stroke hemiplegic patients
    欧美另类亚洲清纯唯美| 在线观看舔阴道视频| 99在线视频只有这里精品首页| 久久久水蜜桃国产精品网| 久久久国产欧美日韩av| 啦啦啦韩国在线观看视频| 久久久国产成人精品二区| 操出白浆在线播放| 如日韩欧美国产精品一区二区三区| 久久精品91无色码中文字幕| 国产不卡一卡二| 麻豆一二三区av精品| 色综合站精品国产| 午夜福利18| 高清在线国产一区| 两人在一起打扑克的视频| 99在线视频只有这里精品首页| 国产色视频综合| 亚洲 欧美 日韩 在线 免费| 自拍欧美九色日韩亚洲蝌蚪91| 日本五十路高清| 91精品国产国语对白视频| 中文字幕人成人乱码亚洲影| av电影中文网址| 啦啦啦韩国在线观看视频| 午夜久久久在线观看| 久久久国产成人免费| 一个人观看的视频www高清免费观看 | 久久人妻福利社区极品人妻图片| 咕卡用的链子| 给我免费播放毛片高清在线观看| 变态另类丝袜制服| 日韩免费av在线播放| 麻豆成人av在线观看| 又紧又爽又黄一区二区| 精品国内亚洲2022精品成人| 免费在线观看亚洲国产| 99国产综合亚洲精品| 999精品在线视频| 如日韩欧美国产精品一区二区三区| 中文字幕人妻丝袜一区二区| 久久人人爽av亚洲精品天堂| 亚洲五月婷婷丁香| 男人舔女人下体高潮全视频| 精品国产亚洲在线| 超碰成人久久| 亚洲午夜理论影院| 日韩大尺度精品在线看网址 | 亚洲avbb在线观看| www.自偷自拍.com| 欧美 亚洲 国产 日韩一| 欧美乱色亚洲激情| 99久久99久久久精品蜜桃| 亚洲一区二区三区不卡视频| 免费在线观看日本一区| 窝窝影院91人妻| 亚洲欧洲精品一区二区精品久久久| x7x7x7水蜜桃| 国产亚洲精品久久久久5区| 免费在线观看完整版高清| 人人妻,人人澡人人爽秒播| 日日摸夜夜添夜夜添小说| 国产麻豆成人av免费视频| 日韩欧美免费精品| 国产又色又爽无遮挡免费看| 老熟妇乱子伦视频在线观看| 日本在线视频免费播放| 日本vs欧美在线观看视频| 天堂√8在线中文| 天天躁狠狠躁夜夜躁狠狠躁| 午夜福利视频1000在线观看 | 制服丝袜大香蕉在线| 欧美av亚洲av综合av国产av| 亚洲一区中文字幕在线| 又黄又爽又免费观看的视频| 9热在线视频观看99| 亚洲 国产 在线| 国语自产精品视频在线第100页| 国产精品av久久久久免费| 露出奶头的视频| 国产蜜桃级精品一区二区三区| 国产精品影院久久| 大型av网站在线播放| av天堂在线播放| 精品久久久精品久久久| 中文字幕最新亚洲高清| 久久精品国产综合久久久| 国产伦人伦偷精品视频| 最好的美女福利视频网| 亚洲精品av麻豆狂野| 美女扒开内裤让男人捅视频| 国产精品久久视频播放| 久久久水蜜桃国产精品网| 亚洲伊人色综图| 91在线观看av| 9色porny在线观看| 午夜精品久久久久久毛片777| 亚洲久久久国产精品| 神马国产精品三级电影在线观看 | 无限看片的www在线观看| 熟女少妇亚洲综合色aaa.| 香蕉国产在线看| 国产av又大| 日韩大码丰满熟妇| 黄色视频,在线免费观看| 国产真人三级小视频在线观看| 久久国产精品影院| 国产成人啪精品午夜网站| 成人免费观看视频高清| 国产又爽黄色视频| 在线观看免费视频日本深夜| 亚洲欧美激情综合另类| 夜夜爽天天搞| 亚洲性夜色夜夜综合| 欧美激情 高清一区二区三区| 男女下面进入的视频免费午夜 | 别揉我奶头~嗯~啊~动态视频| a级毛片在线看网站| 在线观看午夜福利视频| 亚洲,欧美精品.| 亚洲国产中文字幕在线视频| 操出白浆在线播放| 国产精品av久久久久免费| 又紧又爽又黄一区二区| 日韩精品中文字幕看吧| 99精品欧美一区二区三区四区| 麻豆成人av在线观看| 中文字幕久久专区| 三级毛片av免费| 中文字幕人妻丝袜一区二区| 男女午夜视频在线观看| 黑人巨大精品欧美一区二区蜜桃| 亚洲欧美精品综合一区二区三区| 中文字幕av电影在线播放| 成年女人毛片免费观看观看9| 99国产综合亚洲精品| 国产精品久久久人人做人人爽| 精品电影一区二区在线| 大型黄色视频在线免费观看| 午夜福利高清视频| 午夜免费成人在线视频| 国产又爽黄色视频| 亚洲伊人色综图| а√天堂www在线а√下载| 国产99白浆流出| 91在线观看av| 桃色一区二区三区在线观看| 51午夜福利影视在线观看| 国产91精品成人一区二区三区| 精品少妇一区二区三区视频日本电影| 亚洲精品国产一区二区精华液| 啦啦啦韩国在线观看视频| 国产亚洲欧美精品永久| 视频区欧美日本亚洲| 国产真人三级小视频在线观看| 中文字幕最新亚洲高清| 国产区一区二久久| 一区二区三区精品91| 亚洲欧美一区二区三区黑人| 亚洲av日韩精品久久久久久密| 日韩 欧美 亚洲 中文字幕| 91字幕亚洲| 可以在线观看毛片的网站| 国产成人啪精品午夜网站| 国产成人精品在线电影| 精品国产美女av久久久久小说| 在线观看一区二区三区| 在线免费观看的www视频| 亚洲狠狠婷婷综合久久图片| 国产精品自产拍在线观看55亚洲| 国产一区二区在线av高清观看| 精品久久久久久久毛片微露脸| 久久久水蜜桃国产精品网| 久久精品人人爽人人爽视色| 美女高潮到喷水免费观看| avwww免费| 国产精品久久久人人做人人爽| 久久天堂一区二区三区四区| 久久天堂一区二区三区四区| 十八禁网站免费在线| 久久久久久久午夜电影| 51午夜福利影视在线观看| 亚洲aⅴ乱码一区二区在线播放 | 别揉我奶头~嗯~啊~动态视频| 免费在线观看亚洲国产| 校园春色视频在线观看| 久久久精品国产亚洲av高清涩受| 国产精品影院久久| 亚洲情色 制服丝袜| 波多野结衣巨乳人妻| 日韩国内少妇激情av| 国产一区二区三区视频了| 国语自产精品视频在线第100页| 波多野结衣av一区二区av| 91成年电影在线观看| 最近最新中文字幕大全免费视频| 嫁个100分男人电影在线观看| 两个人免费观看高清视频| 亚洲 欧美 日韩 在线 免费| 欧美+亚洲+日韩+国产| 亚洲第一av免费看| 亚洲第一青青草原| 久久婷婷人人爽人人干人人爱 | 美女午夜性视频免费| 欧美日韩精品网址| 制服人妻中文乱码| 一级毛片精品| 757午夜福利合集在线观看| 亚洲午夜精品一区,二区,三区| 成年版毛片免费区| 久久香蕉精品热| 色综合站精品国产| 亚洲成av片中文字幕在线观看| 51午夜福利影视在线观看| 日本vs欧美在线观看视频| 波多野结衣巨乳人妻| 国产亚洲精品久久久久久毛片| 国产真人三级小视频在线观看| 精品久久久久久久毛片微露脸| 久久香蕉精品热| 色在线成人网| 亚洲中文av在线| 一区在线观看完整版| 色尼玛亚洲综合影院| 色播在线永久视频| 亚洲专区国产一区二区| 韩国精品一区二区三区| 久久精品国产清高在天天线| 免费高清在线观看日韩| 亚洲人成电影观看| 久久精品国产亚洲av高清一级| 琪琪午夜伦伦电影理论片6080| 日日干狠狠操夜夜爽| 夜夜爽天天搞| 精品无人区乱码1区二区| 91老司机精品| 国产亚洲精品av在线| 国产精品1区2区在线观看.| 亚洲精品国产区一区二| 久久国产乱子伦精品免费另类| 宅男免费午夜| 丝袜美腿诱惑在线| 1024视频免费在线观看| 人人妻人人澡人人看| 久久草成人影院| 高潮久久久久久久久久久不卡| 亚洲va日本ⅴa欧美va伊人久久| 国产色视频综合| 午夜福利免费观看在线| 黄片小视频在线播放| 日本精品一区二区三区蜜桃| 精品人妻1区二区| 日韩免费av在线播放| 一区二区三区激情视频| 日韩大尺度精品在线看网址 | 午夜福利一区二区在线看| 身体一侧抽搐| 一区二区三区激情视频| 国产精品1区2区在线观看.| 国产精品乱码一区二三区的特点 | 亚洲人成网站在线播放欧美日韩| 在线十欧美十亚洲十日本专区| 91在线观看av| 色综合婷婷激情| 亚洲第一电影网av| 别揉我奶头~嗯~啊~动态视频| 啦啦啦韩国在线观看视频| 国产精华一区二区三区| 可以在线观看的亚洲视频| 久久久精品国产亚洲av高清涩受| 制服人妻中文乱码| 高潮久久久久久久久久久不卡| 一本综合久久免费| 久久婷婷成人综合色麻豆| 精品一品国产午夜福利视频| 亚洲久久久国产精品| 一二三四社区在线视频社区8| 无遮挡黄片免费观看| 亚洲 欧美 日韩 在线 免费| 丰满人妻熟妇乱又伦精品不卡| 色播亚洲综合网| 18禁裸乳无遮挡免费网站照片 | 久久精品国产亚洲av高清一级| 久久人人精品亚洲av| 欧美乱码精品一区二区三区| 国产人伦9x9x在线观看| 久久人妻福利社区极品人妻图片| 1024香蕉在线观看| 精品一区二区三区视频在线观看免费| 欧美 亚洲 国产 日韩一| 免费在线观看视频国产中文字幕亚洲| 一边摸一边抽搐一进一出视频| 免费在线观看视频国产中文字幕亚洲| 好看av亚洲va欧美ⅴa在| 精品国内亚洲2022精品成人| 怎么达到女性高潮| 真人做人爱边吃奶动态| 国产精品永久免费网站| 国产精品电影一区二区三区| 国产精品乱码一区二三区的特点 | 国产高清激情床上av| 国产精品一区二区三区四区久久 | 欧美成狂野欧美在线观看| bbb黄色大片| 精品欧美一区二区三区在线| 国产av一区在线观看免费| 一个人观看的视频www高清免费观看 | 免费在线观看日本一区| 欧洲精品卡2卡3卡4卡5卡区| 欧美色欧美亚洲另类二区 | 精品午夜福利视频在线观看一区| a在线观看视频网站| 制服丝袜大香蕉在线| 欧美丝袜亚洲另类 | 午夜成年电影在线免费观看| 高潮久久久久久久久久久不卡| 老司机深夜福利视频在线观看| 激情视频va一区二区三区| 国产野战对白在线观看| 视频区欧美日本亚洲| 久9热在线精品视频| 91大片在线观看| 一进一出抽搐gif免费好疼| 久久久久亚洲av毛片大全| av视频免费观看在线观看| 91麻豆精品激情在线观看国产| 高清黄色对白视频在线免费看| 日韩欧美在线二视频| 久久香蕉精品热| 久久久久久久午夜电影| 国产精品亚洲av一区麻豆| 欧美国产精品va在线观看不卡| 亚洲国产毛片av蜜桃av| 国产精品av久久久久免费| 91麻豆av在线| 97人妻天天添夜夜摸| 国产成人精品在线电影| 熟女少妇亚洲综合色aaa.| 亚洲专区国产一区二区| 国产av一区二区精品久久| 男女下面插进去视频免费观看| 久久久久久久久免费视频了| 91成人精品电影| 国产色视频综合| www日本在线高清视频| 丰满的人妻完整版| 亚洲精品粉嫩美女一区| 久久久水蜜桃国产精品网| 1024视频免费在线观看| 久久人妻av系列| 日本黄色视频三级网站网址| 麻豆av在线久日| 精品无人区乱码1区二区| 一区二区三区国产精品乱码| 欧美精品亚洲一区二区| 久久久国产欧美日韩av| 精品午夜福利视频在线观看一区| 黄片小视频在线播放| 亚洲aⅴ乱码一区二区在线播放 | 91在线观看av| 变态另类丝袜制服| 欧美成狂野欧美在线观看| 黄片播放在线免费| 无遮挡黄片免费观看| 久久久久国内视频| 国产97色在线日韩免费| 亚洲欧美精品综合一区二区三区| 国产高清激情床上av| 亚洲九九香蕉| 精品福利观看| 亚洲片人在线观看| 色精品久久人妻99蜜桃| 高清在线国产一区| 中文亚洲av片在线观看爽| 亚洲视频免费观看视频| 久久精品91无色码中文字幕| 午夜久久久久精精品| 亚洲一区中文字幕在线| 日韩欧美一区视频在线观看| 国产精品久久久久久人妻精品电影| 午夜成年电影在线免费观看| 麻豆成人av在线观看| 免费在线观看完整版高清| 99riav亚洲国产免费| 亚洲精品国产色婷婷电影| 精品国产一区二区久久| 看片在线看免费视频| 日韩精品青青久久久久久| 俄罗斯特黄特色一大片| 性欧美人与动物交配| 欧美激情久久久久久爽电影 | 欧美色欧美亚洲另类二区 | 嫩草影院精品99| e午夜精品久久久久久久| 亚洲最大成人中文| 精品一品国产午夜福利视频| av在线天堂中文字幕| 自拍欧美九色日韩亚洲蝌蚪91| 九色亚洲精品在线播放| 国产私拍福利视频在线观看| 国产欧美日韩一区二区精品| 老司机深夜福利视频在线观看| 日韩欧美三级三区| av片东京热男人的天堂| 久久国产精品影院| or卡值多少钱| 天堂影院成人在线观看| 99国产极品粉嫩在线观看| 久久天躁狠狠躁夜夜2o2o| 精品久久久久久久人妻蜜臀av | 久久香蕉国产精品| 国产激情久久老熟女| 一区二区三区高清视频在线| 无限看片的www在线观看| 亚洲,欧美精品.| 在线av久久热| 侵犯人妻中文字幕一二三四区| av视频免费观看在线观看| 女警被强在线播放| 久久精品亚洲熟妇少妇任你| 精品国产美女av久久久久小说| 精品欧美国产一区二区三| 久久人妻av系列| 少妇粗大呻吟视频| 搡老妇女老女人老熟妇| 制服人妻中文乱码| 自线自在国产av| 电影成人av| 国产成人免费无遮挡视频| 免费观看精品视频网站| 又紧又爽又黄一区二区| 天堂√8在线中文| 欧美中文综合在线视频| 大型黄色视频在线免费观看| 久久精品国产亚洲av高清一级| www.熟女人妻精品国产| 欧美一区二区精品小视频在线| 国产亚洲精品久久久久5区| 国内久久婷婷六月综合欲色啪| 亚洲国产欧美日韩在线播放| 国产高清有码在线观看视频 | 精品一区二区三区av网在线观看| 男人舔女人的私密视频| 日日干狠狠操夜夜爽| 欧美久久黑人一区二区| 乱人伦中国视频| av超薄肉色丝袜交足视频| 一级毛片精品| 久久亚洲精品不卡| 美国免费a级毛片| 日本 av在线| 日本一区二区免费在线视频| 岛国在线观看网站| 欧美老熟妇乱子伦牲交| 一区在线观看完整版| 少妇熟女aⅴ在线视频| 变态另类丝袜制服| 最新在线观看一区二区三区| 精品国产超薄肉色丝袜足j| 国产成人精品久久二区二区91| 精品久久久久久久久久免费视频| 免费看a级黄色片| 级片在线观看| 久久久久亚洲av毛片大全| 桃红色精品国产亚洲av| 欧美日韩瑟瑟在线播放| 午夜福利免费观看在线| 麻豆一二三区av精品| 久久久久久国产a免费观看| 国产激情欧美一区二区| 午夜福利,免费看| 国产精品亚洲av一区麻豆| 国产99白浆流出| 亚洲国产精品成人综合色| 国产成+人综合+亚洲专区| 亚洲男人天堂网一区| √禁漫天堂资源中文www| 日本免费一区二区三区高清不卡 | 99国产精品99久久久久| 少妇被粗大的猛进出69影院| 老司机午夜十八禁免费视频| 一边摸一边做爽爽视频免费| 亚洲九九香蕉| 黄片播放在线免费| 午夜福利欧美成人| 最近最新免费中文字幕在线| 国产精品秋霞免费鲁丝片| 亚洲五月色婷婷综合| 久久精品人人爽人人爽视色| 亚洲天堂国产精品一区在线| 午夜影院日韩av| 大香蕉久久成人网| 少妇熟女aⅴ在线视频| 国产精品精品国产色婷婷| 欧美大码av| 巨乳人妻的诱惑在线观看| 极品人妻少妇av视频| 国产99白浆流出| av网站免费在线观看视频| 韩国av一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 午夜精品国产一区二区电影| 亚洲男人的天堂狠狠| 欧美一区二区精品小视频在线| 禁无遮挡网站| 亚洲一区中文字幕在线| 国产精品日韩av在线免费观看 | 最好的美女福利视频网| 真人一进一出gif抽搐免费| 精品少妇一区二区三区视频日本电影| 日日干狠狠操夜夜爽| 精品高清国产在线一区| 亚洲男人天堂网一区| 色尼玛亚洲综合影院| 成人三级黄色视频| 国产麻豆成人av免费视频| 亚洲电影在线观看av| 这个男人来自地球电影免费观看| 欧美乱色亚洲激情| 亚洲中文av在线| 精品久久久久久久久久免费视频| 一区福利在线观看| 首页视频小说图片口味搜索| 亚洲一区二区三区色噜噜| 日本 欧美在线| 又黄又爽又免费观看的视频| 嫩草影院精品99| 日日摸夜夜添夜夜添小说| 日韩欧美三级三区| 免费久久久久久久精品成人欧美视频| 最好的美女福利视频网| 久久久久久亚洲精品国产蜜桃av| 999久久久国产精品视频| 99riav亚洲国产免费| 曰老女人黄片| а√天堂www在线а√下载| 亚洲国产欧美网| 日韩三级视频一区二区三区| 精品国内亚洲2022精品成人| 亚洲,欧美精品.| 国产片内射在线| 天堂动漫精品| 一级毛片高清免费大全| 久久久久久大精品| 国产精品久久久久久亚洲av鲁大| 色播亚洲综合网| 又黄又爽又免费观看的视频| 亚洲九九香蕉| 99久久99久久久精品蜜桃| 91成年电影在线观看| 色综合婷婷激情| 一二三四社区在线视频社区8| 狂野欧美激情性xxxx| 日本五十路高清| 久久这里只有精品19| 午夜视频精品福利| 一级,二级,三级黄色视频| 亚洲va日本ⅴa欧美va伊人久久| 国产精品自产拍在线观看55亚洲| av电影中文网址| 中文字幕人妻丝袜一区二区| 一个人观看的视频www高清免费观看 | 久久青草综合色| 成人三级黄色视频| 日韩有码中文字幕| 久久久久久人人人人人| 久9热在线精品视频| 国产高清视频在线播放一区| 真人一进一出gif抽搐免费| 日韩成人在线观看一区二区三区| 99久久国产精品久久久| 亚洲专区中文字幕在线| 久9热在线精品视频| 少妇粗大呻吟视频| 午夜精品在线福利| 在线观看66精品国产| 午夜福利在线观看吧| 一级a爱视频在线免费观看| 国产精华一区二区三区| 亚洲国产高清在线一区二区三 | 啦啦啦免费观看视频1| 色综合欧美亚洲国产小说| 天天添夜夜摸| 精品久久久精品久久久| 男女之事视频高清在线观看| 色婷婷久久久亚洲欧美| 狠狠狠狠99中文字幕| 少妇熟女aⅴ在线视频| 国产高清有码在线观看视频 | 在线视频色国产色| 十分钟在线观看高清视频www| 午夜免费成人在线视频| 亚洲成a人片在线一区二区| 亚洲精品一区av在线观看| 国内精品久久久久久久电影| 波多野结衣高清无吗| 久久久久久久久久久久大奶| 国产精品久久久av美女十八| 成年女人毛片免费观看观看9| 无遮挡黄片免费观看| 精品欧美一区二区三区在线| 巨乳人妻的诱惑在线观看| 日本三级黄在线观看| 91麻豆精品激情在线观看国产| av天堂久久9| 91精品三级在线观看| 国产亚洲精品综合一区在线观看 | 日韩国内少妇激情av| 国产成人精品无人区| 欧美日韩精品网址| 国产精华一区二区三区| 自线自在国产av| 午夜福利成人在线免费观看| 88av欧美| 成熟少妇高潮喷水视频|