• <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
    超碰av人人做人人爽久久| 久久国内精品自在自线图片| 免费看日本二区| 久久久精品大字幕| 亚洲国产高清在线一区二区三| 蜜臀久久99精品久久宅男| 免费高清视频大片| 国产片特级美女逼逼视频| 欧美高清性xxxxhd video| av在线亚洲专区| 三级经典国产精品| 99久久精品热视频| 亚洲精品影视一区二区三区av| 成人鲁丝片一二三区免费| 欧美一级a爱片免费观看看| 91狼人影院| 午夜福利高清视频| 婷婷六月久久综合丁香| 亚洲av熟女| 男人和女人高潮做爰伦理| 成年女人毛片免费观看观看9| 九九爱精品视频在线观看| 久久亚洲国产成人精品v| 日本-黄色视频高清免费观看| 精华霜和精华液先用哪个| 精品久久国产蜜桃| 可以在线观看的亚洲视频| 特大巨黑吊av在线直播| 熟妇人妻久久中文字幕3abv| 午夜福利高清视频| 一级毛片电影观看 | 免费在线观看成人毛片| 男女那种视频在线观看| 日本-黄色视频高清免费观看| 精品国产三级普通话版| 免费大片18禁| 午夜激情欧美在线| 波多野结衣巨乳人妻| 全区人妻精品视频| 国产精品,欧美在线| 久久6这里有精品| 国产中年淑女户外野战色| 精品久久久久久久久久免费视频| 在线天堂最新版资源| 色噜噜av男人的天堂激情| 久久久精品欧美日韩精品| 精品欧美国产一区二区三| 内地一区二区视频在线| 亚洲色图av天堂| 深夜a级毛片| 超碰av人人做人人爽久久| 成人av在线播放网站| 久久久久久久久久久丰满| 欧美激情国产日韩精品一区| 蜜臀久久99精品久久宅男| 国产成人福利小说| 国产中年淑女户外野战色| 精品国内亚洲2022精品成人| 成人av一区二区三区在线看| 变态另类丝袜制服| 国产片特级美女逼逼视频| 男女那种视频在线观看| 国产视频一区二区在线看| 看非洲黑人一级黄片| 国产一区二区在线av高清观看| 亚洲精品国产av成人精品 | 成人国产麻豆网| 18禁在线播放成人免费| 美女 人体艺术 gogo| 欧美一级a爱片免费观看看| 欧美日韩综合久久久久久| 久久精品国产自在天天线| 欧美性感艳星| 欧美3d第一页| 欧美绝顶高潮抽搐喷水| 色在线成人网| 欧美性猛交╳xxx乱大交人| 国产人妻一区二区三区在| 国产高清三级在线| 国内精品美女久久久久久| 日日撸夜夜添| 午夜老司机福利剧场| 一卡2卡三卡四卡精品乱码亚洲| av在线亚洲专区| 日本一二三区视频观看| 中文字幕免费在线视频6| 91狼人影院| 亚洲精品一区av在线观看| 色哟哟·www| 一进一出抽搐动态| 在线国产一区二区在线| 成人三级黄色视频| а√天堂www在线а√下载| 性欧美人与动物交配| 亚洲国产日韩欧美精品在线观看| 五月伊人婷婷丁香| 久久天躁狠狠躁夜夜2o2o| 精品人妻一区二区三区麻豆 | 久久久精品欧美日韩精品| 97超级碰碰碰精品色视频在线观看| 中文资源天堂在线| 最新中文字幕久久久久| 久久久久九九精品影院| 人妻丰满熟妇av一区二区三区| 日韩欧美 国产精品| 国国产精品蜜臀av免费| 免费看日本二区| 看十八女毛片水多多多| АⅤ资源中文在线天堂| 97超级碰碰碰精品色视频在线观看| 欧美+日韩+精品| 亚洲欧美日韩无卡精品| 精品久久久久久久末码| 99精品在免费线老司机午夜| 欧美日韩在线观看h| 精品免费久久久久久久清纯| 国产在视频线在精品| 精品国内亚洲2022精品成人| 亚洲自拍偷在线| 亚洲av第一区精品v没综合| 内地一区二区视频在线| 成人毛片a级毛片在线播放| 久久久久久久午夜电影| 人妻久久中文字幕网| 婷婷亚洲欧美| 99热网站在线观看| 美女cb高潮喷水在线观看| 精品久久久噜噜| 日韩人妻高清精品专区| 超碰av人人做人人爽久久| 免费黄网站久久成人精品| 亚洲经典国产精华液单| 国产精品三级大全| 特大巨黑吊av在线直播| 12—13女人毛片做爰片一| 亚洲无线在线观看| 中文字幕免费在线视频6| 你懂的网址亚洲精品在线观看 | 男女做爰动态图高潮gif福利片| 亚洲精品久久国产高清桃花| 国产激情偷乱视频一区二区| 国产av在哪里看| 亚洲国产精品成人综合色| 精品人妻一区二区三区麻豆 | 国产成人影院久久av| 欧美潮喷喷水| 人人妻人人澡欧美一区二区| 成人无遮挡网站| 日韩大尺度精品在线看网址| 无遮挡黄片免费观看| 91久久精品国产一区二区成人| 尤物成人国产欧美一区二区三区| 亚洲av中文av极速乱| 欧美激情在线99| 中文在线观看免费www的网站| 国产精品一区二区性色av| 色综合站精品国产| 亚洲国产精品久久男人天堂| 久久热精品热| 国产麻豆成人av免费视频| 97超碰精品成人国产| 精华霜和精华液先用哪个| 俺也久久电影网| 中文字幕av在线有码专区| 久久国产乱子免费精品| 夜夜爽天天搞| 少妇的逼水好多| 亚洲自拍偷在线| 日韩高清综合在线| 亚洲无线观看免费| 丰满人妻一区二区三区视频av| 午夜激情欧美在线| 成人特级黄色片久久久久久久| 色哟哟·www| 如何舔出高潮| 亚洲第一区二区三区不卡| 久久久精品94久久精品| 成熟少妇高潮喷水视频| av福利片在线观看| 国国产精品蜜臀av免费| 亚洲美女搞黄在线观看 | 干丝袜人妻中文字幕| 亚洲av二区三区四区| 久久精品夜夜夜夜夜久久蜜豆| 毛片一级片免费看久久久久| 精品久久久噜噜| 亚洲美女黄片视频| 国产亚洲精品久久久com| 久久久久免费精品人妻一区二区| 亚州av有码| 99热全是精品| 熟女电影av网| 日本a在线网址| 久久人妻av系列| 22中文网久久字幕| 天天躁日日操中文字幕| 91狼人影院| 国产精品国产三级国产av玫瑰| 在线免费观看不下载黄p国产| 黄色配什么色好看| 色吧在线观看| 久久草成人影院| h日本视频在线播放| 国产精品三级大全| 国产色爽女视频免费观看| 日日啪夜夜撸| 给我免费播放毛片高清在线观看| 尾随美女入室| 美女黄网站色视频| 中文字幕免费在线视频6| 老熟妇仑乱视频hdxx| 色播亚洲综合网| 亚洲在线自拍视频| 人人妻人人澡欧美一区二区| 成人av在线播放网站| 大香蕉久久网| 亚洲最大成人av| 内地一区二区视频在线| 亚洲欧美日韩无卡精品| 国内精品宾馆在线| 国产成人freesex在线 | 不卡视频在线观看欧美| 天堂√8在线中文| 校园春色视频在线观看| 人人妻,人人澡人人爽秒播| 最后的刺客免费高清国语| 尾随美女入室| 久久久久久久午夜电影| 久久久久久久久久久丰满| 热99在线观看视频| 国产69精品久久久久777片| 99视频精品全部免费 在线| 国产精品精品国产色婷婷| 99热这里只有是精品在线观看| 亚洲国产高清在线一区二区三| 少妇裸体淫交视频免费看高清| 男女那种视频在线观看| 国产伦精品一区二区三区视频9| 秋霞在线观看毛片| 在线观看免费视频日本深夜| 综合色丁香网| 国产在视频线在精品| 2021天堂中文幕一二区在线观| 久久久久国内视频| 精品福利观看| 综合色丁香网| 国产aⅴ精品一区二区三区波| 久久热精品热| 麻豆国产av国片精品| 伦精品一区二区三区| 深夜精品福利| 亚洲美女搞黄在线观看 | 在线免费观看不下载黄p国产| 小说图片视频综合网站| 夜夜爽天天搞| 国产伦在线观看视频一区| 亚洲av成人av| 国产大屁股一区二区在线视频| 亚洲美女黄片视频| av女优亚洲男人天堂| 亚洲最大成人中文| 国产精品女同一区二区软件| 免费一级毛片在线播放高清视频| 97超视频在线观看视频| 尤物成人国产欧美一区二区三区| 国产不卡一卡二| 熟女人妻精品中文字幕| 熟女电影av网| 国模一区二区三区四区视频| 国产黄片美女视频| av黄色大香蕉| 亚洲成人中文字幕在线播放| 天天躁日日操中文字幕| 欧美激情国产日韩精品一区| 国产精品精品国产色婷婷| 日韩大尺度精品在线看网址| 黑人高潮一二区| 蜜桃久久精品国产亚洲av| 欧美又色又爽又黄视频| 51国产日韩欧美| 女同久久另类99精品国产91| 卡戴珊不雅视频在线播放| 精品久久久久久久人妻蜜臀av| 日韩制服骚丝袜av| 一卡2卡三卡四卡精品乱码亚洲| 亚洲精品在线观看二区| 国产真实乱freesex| 在线观看免费视频日本深夜| 午夜福利在线观看吧| 不卡视频在线观看欧美| 日韩成人av中文字幕在线观看 | а√天堂www在线а√下载| 午夜视频国产福利| 久久精品综合一区二区三区| 久久久a久久爽久久v久久| 国内揄拍国产精品人妻在线| 免费在线观看成人毛片| 亚洲精品日韩在线中文字幕 | 免费看av在线观看网站| 久久欧美精品欧美久久欧美| 51国产日韩欧美| 欧美最新免费一区二区三区| 日本爱情动作片www.在线观看 | 久久韩国三级中文字幕| 老熟妇乱子伦视频在线观看| 少妇裸体淫交视频免费看高清| 亚洲人成网站在线观看播放| 国产一区二区在线av高清观看| 国产色爽女视频免费观看| 婷婷六月久久综合丁香| 久久精品国产亚洲av天美| av国产免费在线观看| 非洲黑人性xxxx精品又粗又长| 人妻夜夜爽99麻豆av| 亚洲欧美清纯卡通| 免费av观看视频| 天堂√8在线中文| 久久久久久久亚洲中文字幕| 国产成年人精品一区二区| 精品久久久久久成人av| 婷婷色综合大香蕉| 99久国产av精品| 99国产极品粉嫩在线观看| 毛片女人毛片| 欧美日韩一区二区视频在线观看视频在线 | 亚洲国产精品合色在线| 校园人妻丝袜中文字幕| 欧美+日韩+精品| 舔av片在线| 小蜜桃在线观看免费完整版高清| 97在线视频观看| 亚洲av成人精品一区久久| 久久久久久伊人网av| 干丝袜人妻中文字幕| 最好的美女福利视频网| 欧美色欧美亚洲另类二区| 久久久久久九九精品二区国产| 国产精品嫩草影院av在线观看| 在线观看66精品国产| 男女视频在线观看网站免费| 免费av不卡在线播放| 国产高清有码在线观看视频| 天天一区二区日本电影三级| 日日摸夜夜添夜夜爱| 亚洲,欧美,日韩| 一级毛片电影观看 | 国产中年淑女户外野战色| 亚洲第一区二区三区不卡| 一本久久中文字幕| 亚洲国产高清在线一区二区三| 性色avwww在线观看| 在线a可以看的网站| 美女黄网站色视频| 亚洲av中文av极速乱| 欧美日韩精品成人综合77777| 国产成人一区二区在线| 小说图片视频综合网站| 亚洲美女视频黄频| 成人漫画全彩无遮挡| 国产成人a∨麻豆精品| 欧美成人精品欧美一级黄| 联通29元200g的流量卡| 少妇高潮的动态图| 一a级毛片在线观看| 日日干狠狠操夜夜爽| 午夜精品国产一区二区电影 | 亚洲第一电影网av| 国产伦精品一区二区三区四那| 婷婷亚洲欧美| 少妇丰满av| 国产成人a∨麻豆精品| 99热这里只有精品一区| 欧美+日韩+精品| av中文乱码字幕在线| 淫秽高清视频在线观看| 国产日本99.免费观看| 看片在线看免费视频| 日韩av不卡免费在线播放| 久久精品综合一区二区三区| 国产黄色视频一区二区在线观看 | 97人妻精品一区二区三区麻豆| 久久精品综合一区二区三区| 精品无人区乱码1区二区| 亚洲性久久影院| 欧美xxxx性猛交bbbb| 欧美激情国产日韩精品一区| 日本 av在线| 午夜福利视频1000在线观看| 乱码一卡2卡4卡精品| 99久久久亚洲精品蜜臀av| 久久精品国产清高在天天线| 毛片一级片免费看久久久久| 亚洲第一电影网av| 12—13女人毛片做爰片一| 一级毛片电影观看 | 欧洲精品卡2卡3卡4卡5卡区| 国语自产精品视频在线第100页| 成年女人永久免费观看视频| 熟女人妻精品中文字幕| 又黄又爽又免费观看的视频| 大又大粗又爽又黄少妇毛片口| 国产又黄又爽又无遮挡在线| 欧美区成人在线视频| 久久精品国产亚洲av香蕉五月| 俺也久久电影网| 亚洲成av人片在线播放无| 看非洲黑人一级黄片| 亚洲五月天丁香| 国产色婷婷99| 日本免费一区二区三区高清不卡| 99九九线精品视频在线观看视频| 美女高潮的动态| 欧美激情久久久久久爽电影| 人妻夜夜爽99麻豆av| a级毛片a级免费在线| 女人被狂操c到高潮| 亚洲无线观看免费| 少妇人妻一区二区三区视频| 日本爱情动作片www.在线观看 | 婷婷精品国产亚洲av| 久久久久久久久大av| 亚洲欧美日韩高清专用| 久久久欧美国产精品| 亚洲人成网站在线播放欧美日韩| a级毛片a级免费在线| 欧美成人精品欧美一级黄| 人人妻人人澡人人爽人人夜夜 | 亚洲国产高清在线一区二区三| 成人国产麻豆网| videossex国产| 欧美最黄视频在线播放免费| 99视频精品全部免费 在线| 国产黄a三级三级三级人| 国产精品久久久久久久久免| 国产高潮美女av| 欧美性猛交╳xxx乱大交人| 国产久久久一区二区三区| 一级黄色大片毛片| 国产一区二区激情短视频| 国产精品久久久久久久久免| 插逼视频在线观看| 99久久中文字幕三级久久日本| 一本久久中文字幕| 91在线观看av| 国产成人91sexporn| 亚洲欧美日韩高清专用| 国产精品人妻久久久久久| 美女免费视频网站| 一级毛片我不卡| 日韩一区二区视频免费看| 国产午夜精品久久久久久一区二区三区 | 一进一出抽搐gif免费好疼| 91在线精品国自产拍蜜月| 老司机午夜福利在线观看视频| 国产精品一区二区三区四区免费观看 | 97人妻精品一区二区三区麻豆| 超碰av人人做人人爽久久| 免费看a级黄色片| 亚洲av.av天堂| 丝袜美腿在线中文| 亚洲国产精品国产精品| 中文字幕av在线有码专区| 熟女电影av网| 自拍偷自拍亚洲精品老妇| 久久久精品欧美日韩精品| 午夜爱爱视频在线播放| 男女那种视频在线观看| 99久久九九国产精品国产免费| av在线观看视频网站免费| 男女视频在线观看网站免费| 精品99又大又爽又粗少妇毛片| 婷婷精品国产亚洲av| 久久久久久久久大av| 国产精品一区二区三区四区久久| 日本一二三区视频观看| 精品国内亚洲2022精品成人| 午夜视频国产福利| 欧美日韩一区二区视频在线观看视频在线 | 亚洲国产精品国产精品| 身体一侧抽搐| 久久精品国产亚洲av香蕉五月| 欧美性猛交╳xxx乱大交人| 国产精品一二三区在线看| 欧美高清成人免费视频www| 无遮挡黄片免费观看| 日日摸夜夜添夜夜爱| 亚洲图色成人| 高清午夜精品一区二区三区 | 欧美性猛交黑人性爽| 亚洲国产精品国产精品| 国产一级毛片七仙女欲春2| 欧美精品国产亚洲| 午夜激情欧美在线| 国内精品一区二区在线观看| 99久久中文字幕三级久久日本| 一级av片app| 婷婷精品国产亚洲av| 国产亚洲91精品色在线| 在线播放无遮挡| 毛片一级片免费看久久久久| 身体一侧抽搐| 成年女人毛片免费观看观看9| 热99在线观看视频| 啦啦啦观看免费观看视频高清| 我要看日韩黄色一级片| h日本视频在线播放| 午夜精品在线福利| 成人毛片a级毛片在线播放| a级毛色黄片| 菩萨蛮人人尽说江南好唐韦庄 | 午夜爱爱视频在线播放| 日韩中字成人| av在线老鸭窝| 我要看日韩黄色一级片| 精品不卡国产一区二区三区| 亚洲人与动物交配视频| 一a级毛片在线观看| 成人综合一区亚洲| 狂野欧美激情性xxxx在线观看| 黄色配什么色好看| 欧美xxxx黑人xx丫x性爽| 老熟妇仑乱视频hdxx| 成年女人毛片免费观看观看9| 精品乱码久久久久久99久播| 国产午夜精品论理片| 嫩草影院入口| 三级国产精品欧美在线观看| 国产精品伦人一区二区| 国产精品一区二区三区四区久久| 一个人看视频在线观看www免费| 免费观看的影片在线观看| 久久久午夜欧美精品| 久久精品国产鲁丝片午夜精品| 亚洲av成人av| 欧美人与善性xxx| 性色avwww在线观看| 亚洲精品日韩在线中文字幕 | 国产精品一区二区三区四区免费观看 | 亚洲精华国产精华液的使用体验 | 一区二区三区四区激情视频 | 久久久久国产精品人妻aⅴ院| 九九热线精品视视频播放| 一级毛片aaaaaa免费看小| 国产高清视频在线观看网站| 尤物成人国产欧美一区二区三区| 五月伊人婷婷丁香| 色av中文字幕| 人人妻人人澡欧美一区二区| 日韩 亚洲 欧美在线| 天天躁夜夜躁狠狠久久av| 在现免费观看毛片| 小说图片视频综合网站| 最后的刺客免费高清国语| 天堂av国产一区二区熟女人妻| 欧美日本亚洲视频在线播放| 别揉我奶头~嗯~啊~动态视频| 人人妻人人澡人人爽人人夜夜 | 99久久中文字幕三级久久日本| 久久精品国产鲁丝片午夜精品| 日韩精品有码人妻一区| 女人被狂操c到高潮| 最好的美女福利视频网| 在线免费观看不下载黄p国产| 国产高潮美女av| 国产片特级美女逼逼视频| av在线蜜桃| 99在线视频只有这里精品首页| 99久久久亚洲精品蜜臀av| 午夜精品在线福利| 欧美一区二区精品小视频在线| 亚洲aⅴ乱码一区二区在线播放| 男女那种视频在线观看| 99久久中文字幕三级久久日本| 亚洲图色成人| av天堂在线播放| 日本一二三区视频观看| 国产午夜福利久久久久久| 久久韩国三级中文字幕| 日韩精品中文字幕看吧| 丝袜美腿在线中文| 亚洲第一区二区三区不卡| 女人被狂操c到高潮| 久久精品影院6| 大香蕉久久网| 最近2019中文字幕mv第一页| 久久婷婷人人爽人人干人人爱| 99热精品在线国产| 天天躁日日操中文字幕| 免费看a级黄色片| 国内精品久久久久精免费| 午夜亚洲福利在线播放| 精品一区二区三区人妻视频| 国产精品久久电影中文字幕| 狂野欧美激情性xxxx在线观看| 亚洲国产高清在线一区二区三| 成人二区视频| 亚洲内射少妇av| 三级经典国产精品| 国产日本99.免费观看| 丰满人妻一区二区三区视频av| 搞女人的毛片| 国产成人aa在线观看| www.色视频.com| 国产午夜精品论理片| 亚洲国产色片| 99热网站在线观看| 亚洲,欧美,日韩| 99在线视频只有这里精品首页| 99久久精品热视频| 熟妇人妻久久中文字幕3abv| 最近在线观看免费完整版| 亚洲美女视频黄频| 国产精品亚洲美女久久久| 能在线免费观看的黄片| 日韩欧美一区二区三区在线观看|