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

    基于模態(tài)疊加法的面壓載荷與基礎(chǔ)激勵等效方法

    2017-04-12 03:31:46費(fèi)慶國吳邵慶李彥斌
    關(guān)鍵詞:頻響方根軸向

    張 鵬 費(fèi)慶國 吳邵慶 李彥斌

    (東南大學(xué)土木工程學(xué)院, 南京 210096)

    基于模態(tài)疊加法的面壓載荷與基礎(chǔ)激勵等效方法

    張 鵬 費(fèi)慶國 吳邵慶 李彥斌

    (東南大學(xué)土木工程學(xué)院, 南京 210096)

    基于模態(tài)疊加法,將隨機(jī)面壓載荷等效為隨機(jī)基礎(chǔ)激勵,為薄壁構(gòu)件的振動臺試驗替代成本較高的噪聲試驗提供依據(jù).該方法要求等效前后2種載荷作用下主模態(tài)位移一致,將單源載荷等效這一超靜定問題合理簡化為靜定問題.以炭/炭復(fù)合材料懸臂板為研究對象,來驗證載荷等效方法的有效性.分析結(jié)果顯示,在隨機(jī)面壓載荷及等效隨機(jī)基礎(chǔ)激勵作用下,板上典型部位隨機(jī)正應(yīng)力響應(yīng)的功率譜密度高度吻合,在分析頻段內(nèi)危險點(diǎn)處x軸向正應(yīng)力響應(yīng)均方根分析誤差低于1%,分布在板上的5個響應(yīng)點(diǎn)處x軸向正應(yīng)力響應(yīng)均方根分析誤差低于2%.該方法能有效地將隨機(jī)面壓載荷等效為基礎(chǔ)激勵,保證等效前后2種載荷作用下薄壁構(gòu)件的隨機(jī)動應(yīng)力響應(yīng)在頻域及空間域上的分布一致性.

    面壓載荷; 基礎(chǔ)激勵; 載荷等效; 替代試驗; 模態(tài)疊加法

    高強(qiáng)隨機(jī)脈動壓力作用下先進(jìn)復(fù)合材料結(jié)構(gòu)的力學(xué)行為是工程上較為關(guān)注的問題之一[1-3].在高強(qiáng)隨機(jī)脈動壓力作用下,復(fù)合材料薄壁結(jié)構(gòu)可能會發(fā)生疲勞破壞[4].在設(shè)計階段,必須通過地面試驗分析校核結(jié)構(gòu)的可靠性.然而,噪聲試驗對設(shè)備的要求較高.在構(gòu)件級設(shè)計階段,復(fù)合材料的性能具有較大離散性,需開展批量試驗分析,從而導(dǎo)致噪聲試驗的實(shí)施成本較高.因此,尋求一種節(jié)約成本的替代試驗分析方法具有現(xiàn)實(shí)意義.振動臺基礎(chǔ)激勵對平板類構(gòu)件的作用效果類似于面壓載荷,等同于一種分布載荷.振動臺可以同時對多個較小尺寸構(gòu)件施加基礎(chǔ)激勵,振動臺試驗的實(shí)施成本較低,因而采用振動臺基礎(chǔ)激勵試驗代替噪聲試驗具有較高的可行性.開展替代試驗的關(guān)鍵是進(jìn)行載荷等效,即將隨機(jī)面壓載荷等效為振動臺基礎(chǔ)激勵的強(qiáng)迫位移載荷.

    載荷等效在工程中得到了廣泛的應(yīng)用,主要包括將動載荷等效為靜載荷[5-10]以及將動載荷等效為動載荷[11]兩種情況.將動載荷等效為靜載荷可有效簡化結(jié)構(gòu)在動載荷作用下的響應(yīng)分析過程,提高結(jié)構(gòu)強(qiáng)度、剛度設(shè)計以及結(jié)構(gòu)優(yōu)化設(shè)計的效率.該類方法主要包括加速度等效靜載荷法[5]、位移等效靜載荷法[6-7]、應(yīng)力等效靜載荷法[8-9]和能量等效靜載荷法[10],其等效原則為,在等效前后2種載荷作用下,結(jié)構(gòu)特定響應(yīng)量相等.將動載荷等效為動載荷一般服務(wù)于疲勞壽命評估,其目的主要是簡化疲勞試驗加載、加速疲勞試驗過程,等效原則為,等效前后結(jié)構(gòu)關(guān)鍵部位的損傷量相同或呈特定比例[11].

    本文以C/C復(fù)合材料懸臂板為研究對象,基于模態(tài)疊加法,將板面受到的隨機(jī)面壓載荷等效為基礎(chǔ)激勵的隨機(jī)強(qiáng)迫位移載荷.該載荷等效方法將動載荷等效為動載荷.將動應(yīng)力響應(yīng)作為疲勞壽命評估的輸入條件.為服務(wù)于疲勞壽命評估,要求在隨機(jī)面壓載荷和等效基礎(chǔ)激勵的隨機(jī)強(qiáng)迫位移載荷這2種載荷作用下,復(fù)合材料板上的隨機(jī)動應(yīng)力響應(yīng)一致.特定頻率點(diǎn)處結(jié)構(gòu)的最大動應(yīng)力響應(yīng)多由單個主模態(tài)主導(dǎo),因此上述載荷等效的直接目標(biāo)為,載荷等效前后分析頻率點(diǎn)處主模態(tài)的隨機(jī)位移響應(yīng)一致.分析結(jié)果表明,以主模態(tài)位移為直接等效目標(biāo),可使等效前后2種載荷作用下,結(jié)構(gòu)上的隨機(jī)動應(yīng)力響應(yīng)在頻域及空間域上分布一致.

    1 基礎(chǔ)理論

    1.1 模態(tài)疊加法

    對于一個具有N自由度的線性平板結(jié)構(gòu),在板一側(cè)的面壓載荷p(ω)作用下(見圖1(a)),頻域內(nèi)的動力學(xué)方程為

    (-ω2M+jωC+K)X(ω)=Lp(ω)

    (1)

    式中,M,C,K分別為N×N階質(zhì)量矩陣、阻尼矩陣和剛度矩陣;X(ω)為結(jié)構(gòu)在頻域內(nèi)的N×1階位移響應(yīng)向量,其中ω為角頻率;L為N×1階面壓載荷影響向量,取決于面壓載荷的有限元離散.式(1)由時域變換而來,該變換的初始條件為,初始位移與初始速度均為零.

    (a) 面壓載荷 (b) 基礎(chǔ)激勵

    模態(tài)坐標(biāo)系下頻域內(nèi)結(jié)構(gòu)動力學(xué)方程為

    (-ω2Mq+jωCq+Kq)Xq(ω)=ΦTLp(ω)

    (2)

    式中,Φ=[φ1φ2…φq]為N×q階滿足質(zhì)量歸一化條件的正則模態(tài)振型矩陣,其中,φi為第i階模態(tài)振型,q為模態(tài)截斷后的模態(tài)數(shù);Mq,Cq,Kq分別為q×q階模態(tài)質(zhì)量矩陣、模態(tài)阻尼矩陣和模態(tài)剛度矩陣,且Mq=ΦTMΦ=I;Xq(ω)為結(jié)構(gòu)在頻域內(nèi)的q×1階模態(tài)位移響應(yīng)向量,且X(ω)≈ΦXq(ω).

    求解方程(2),可得第i階模態(tài)的位移響應(yīng)xi(ω)為

    xi(ω)=Hip(ω)p(ω)

    (3)

    式中,Hip(ω)為第i階模態(tài)位移響應(yīng)對面壓載荷p(ω)的頻響,即結(jié)構(gòu)在單位面壓載荷作用下第i階模態(tài)的位移響應(yīng),且

    (4)

    式中,ωi為第i階模態(tài)的固有頻率.

    進(jìn)一步可求得隨機(jī)面壓作用下結(jié)構(gòu)在任一自由度k上隨機(jī)位移響應(yīng)的自功率譜密度為

    (5)

    式中,Spp(ω)為隨機(jī)面壓載荷的自功率譜密度.

    1.2 等效基礎(chǔ)激勵

    令結(jié)構(gòu)在固支或簡支邊界上的節(jié)點(diǎn)為約束節(jié)點(diǎn),其余節(jié)點(diǎn)為非約束節(jié)點(diǎn).當(dāng)結(jié)構(gòu)受到基礎(chǔ)作用時,約束節(jié)點(diǎn)的位移等于基礎(chǔ)激勵的強(qiáng)迫位移載荷xe(ω) (見圖1(b)),結(jié)構(gòu)在頻域內(nèi)的動力學(xué)方程為

    (-ω2M+jωC+K)X(ω)=0

    (6)

    式中,結(jié)構(gòu)在頻域內(nèi)的位移響應(yīng)X(ω)可表示為基礎(chǔ)激勵強(qiáng)迫位移xe(ω)與非約束節(jié)點(diǎn)相對約束節(jié)點(diǎn)的相對位移Xm(ω)之和,即

    (7)

    式中,Xm(ω)為m×1階相對位移向量,其中,m+n=N,m為非約束自由度個數(shù),n為約束自由度個數(shù);E為N×1階強(qiáng)迫位移載荷的影響向量,對應(yīng)xe(ω)方向位移自由度的元素為1,其余為0.

    (8)

    將式(7)和(8)代入式(6)可得

    (-ω2M+jωC+K)Y(ω)=Q(ω)xe(ω)

    (9)

    式中,Q(ω)為基礎(chǔ)激勵作用下等效節(jié)點(diǎn)力激勵的N×1階影響向量,且

    Q(ω)=-(-ω2M+jωC+K)E

    (10)

    式(9)與式(1)具有相同的形式,等同于在分布載荷作用下的結(jié)構(gòu)動力學(xué)方程.可通過模態(tài)疊加法求解式(9),得到第i階模態(tài)的位移響應(yīng)為

    yi(ω)=Hie(ω)xe(ω)

    (11)

    式中,Hie(ω)為第i階模態(tài)位移響應(yīng)對基礎(chǔ)強(qiáng)迫位移載荷xe(ω)的頻響,即結(jié)構(gòu)在單位基礎(chǔ)強(qiáng)迫位移載荷作用下第i階模態(tài)的位移響應(yīng),且

    (12)

    進(jìn)一步可求得隨機(jī)基礎(chǔ)激勵作用下結(jié)構(gòu)在任一自由度k上隨機(jī)相對位移響應(yīng)的自功率譜密度為

    (13)

    式中,See(ω)為基礎(chǔ)激勵隨機(jī)強(qiáng)迫位移的自功率譜密度.

    按位移等效原則,當(dāng)基礎(chǔ)激勵作用下任一自由度k上的相對位移響應(yīng)與面壓載荷作用下該自由度上的絕對位移響應(yīng)一致時,2種載荷作用下結(jié)構(gòu)上的應(yīng)力分布也一致,此時可認(rèn)為該基礎(chǔ)激勵為面壓載荷的等效載荷,有

    (14)

    面壓載荷及基礎(chǔ)激勵均為單源載荷,因此式(14)為超靜定問題.結(jié)構(gòu)低頻響應(yīng)一般由頻率點(diǎn)處的主模態(tài)主導(dǎo),忽略其余模態(tài)對響應(yīng)的貢獻(xiàn),式(14)可簡化為

    (15)

    式中,Hdp(ω)為第d階模態(tài)位移響應(yīng)對面壓載荷p(ω) 的頻響;Hde(ω)為第d階模態(tài)位移響應(yīng)對面壓載荷xe(ω)的頻響.

    基于式(15),可求得等效基礎(chǔ)激勵的強(qiáng)迫位移功率譜為

    (16)

    根據(jù)式(16),將隨機(jī)面壓載荷等效為基礎(chǔ)激勵的隨機(jī)強(qiáng)迫位移載荷,其等效原則為,等效前后主模態(tài)的位移響應(yīng)一致.由式(16)獲得的等效載荷可最大限度地滿足式(14)所示的超靜定方程.

    2 算例驗證

    為驗證所提載荷等效方法的有效性,以一開孔C/C復(fù)合材料懸臂板為研究對象,開展仿真研究.該懸臂板的有限元模型如圖2所示.板尺寸為400 mm×300 mm×3 mm,板材料參數(shù)如下:x,y,z方向上的軸向拉伸模量分別為Ex=60 GPa,Ey=55 GPa,Ez=40 GPa;xy,xz,yz平面上的剪切模量分別為Gxy=32 GPa,Gxz=30 GPa,Gyz=29 GPa;xy,xz,yz平面上的泊松比分別為μxy=0.05,μxz=0.09,μyz=0.08;密度ρ=1 800 kg/m3.在一短邊處30 mm范圍內(nèi)采用夾具夾持以實(shí)現(xiàn)固定約束.

    圖2 C/C復(fù)合材料懸臂板有限元模型(單位:mm)

    2.1 模態(tài)位移頻響分析

    在懸臂板上表面施加單位面壓載荷,分析獲得板前39階模態(tài)在0~2 000Hz內(nèi)的位移,即板模態(tài)位移響應(yīng)對面壓載荷的頻響Hip(ω),結(jié)果見圖3(a).類似地,可分析獲得板模態(tài)位移響應(yīng)對基礎(chǔ)強(qiáng)迫位移載荷的頻響Hie(ω),結(jié)果見圖3(b).

    (a) 面壓載荷

    (b) 基礎(chǔ)強(qiáng)迫位移載荷

    圖3中結(jié)果顯示,在分析頻段內(nèi),懸臂板第1,3,6,10,17,24階模態(tài)對2種載荷的頻響,在其固有頻率附近頻段內(nèi)的取值遠(yuǎn)大于其余模態(tài).這些模態(tài)在相應(yīng)頻段內(nèi)即為主模態(tài),板振動響應(yīng)由主模態(tài)主導(dǎo).圖4給出懸臂板在0~2 000 Hz內(nèi)各主模態(tài)振型及固有頻率fi.由圖可知,各主模態(tài)具有一個共同的特征:模態(tài)振型在短邊方向(y軸向)的波數(shù)為零.這表明相對于其他模態(tài),該類模態(tài)更易被面壓載荷及基礎(chǔ)強(qiáng)迫位移載荷激發(fā).

    2.2 載荷等效

    噪聲試驗廣泛采用中華人民共和國軍用標(biāo)準(zhǔn)(GJB 150.17—86)推薦的用于噪聲試驗的1/3倍頻程譜.本研究將該噪聲1/3倍頻程譜的上限轉(zhuǎn)換成面壓功率譜密度Spp,如圖5所示(雙對數(shù)坐標(biāo)),其中基準(zhǔn)功率譜密度S0為

    S0=10(Lp-L0)/10

    (17)

    式中,Lp為噪聲總聲壓級;L0=123.1dB.

    采用式(16)計算等效載荷前,需確定各頻率點(diǎn)處的主模態(tài).圖3中結(jié)果顯示,基于模態(tài)位移對2種載荷的頻響所確定的主模態(tài)存在差異.例如,由圖3(a)可確定第1階模態(tài)在多頻段內(nèi)為主模態(tài), 但由圖3(b)可判定第1階模態(tài)僅在起始頻段內(nèi)為主模態(tài).經(jīng)比較,根據(jù)圖3(b)中結(jié)果確定主模態(tài)階次更合理.載荷等效過程中,各頻率點(diǎn)處主模態(tài)階次d如圖6所示.

    (a) 1階,f1=20 Hz

    (b) 3階,f3=125 Hz

    (c) 6階,f6=350 Hz

    (d) 10階,f10=681 Hz

    (e) 17階,f17=1 132 Hz

    (f) 24階,f24=1 687 Hz

    圖5 雙對數(shù)坐標(biāo)下的隨機(jī)面壓載荷功率譜密度

    在確定主模態(tài)的基礎(chǔ)上,采用式(16)計算得到不同量級面壓載荷下等效基礎(chǔ)激勵強(qiáng)迫位移的功率譜密度(見圖7).由圖可知,在主模態(tài)發(fā)生更替的臨界頻率處,等效基礎(chǔ)激勵強(qiáng)迫位移功率譜的大小會發(fā)生突變.為考查圖7中等效載荷的有效性,對載荷等效前后懸臂板上的應(yīng)力響應(yīng)進(jìn)行了對比.

    圖6 不同頻段內(nèi)主模態(tài)的階次

    圖7 等效基礎(chǔ)激勵強(qiáng)迫位移的功率譜密度

    圖8給出了等效前后2種載荷作用下板上A1,A2兩個典型位置處的應(yīng)力響應(yīng)功率譜密度,隨機(jī)面壓載荷總聲壓級為140 dB.由圖8可知,載荷等效前后板上典型位置處正應(yīng)力響應(yīng)的功率譜密度吻合較好,僅在遠(yuǎn)離峰值頻率的頻段內(nèi)出現(xiàn)差異,該差異對應(yīng)力響應(yīng)均方根的影響甚微.在總聲壓級為140 dB的面壓載荷作用下,載荷等效前后,板上A1~A5五個響應(yīng)點(diǎn)處x軸向正應(yīng)力響應(yīng)在0~2 000 Hz范圍內(nèi)的均方根如表1所示.

    (a) A1

    (b) A2

    圖8 140 dB面壓載荷等效前后板上x軸向正應(yīng)力響應(yīng)功率譜密度

    表1 x軸向正應(yīng)力響應(yīng)均方根

    表1中結(jié)果顯示,最大應(yīng)力響應(yīng)出現(xiàn)在固支端,即危險點(diǎn)出現(xiàn)在固支端.載荷等效引起的A1~A5五個響應(yīng)點(diǎn)處x軸向正應(yīng)力響應(yīng)均方根計算誤差不超過2%,危險點(diǎn)處x軸向正應(yīng)力響應(yīng)均方根誤差不超過1%,這對結(jié)構(gòu)的強(qiáng)度設(shè)計及疲勞評估所產(chǎn)生的影響可忽略不計.圖8及表1中分析結(jié)果驗證了所提載荷等效方法的有效性.

    3 結(jié)語

    本文借助仿真分析手段,基于模態(tài)疊加法,將隨機(jī)面壓載荷等效為隨機(jī)基礎(chǔ)激勵.該方法以主模態(tài)的位移為等效目標(biāo),忽略了其余模態(tài)對結(jié)構(gòu)位移響應(yīng)的貢獻(xiàn),將單源載荷等效這一超靜定問題合理地轉(zhuǎn)化為靜定問題,簡單易行.研究結(jié)果顯示,在等效基礎(chǔ)激勵作用下,結(jié)構(gòu)典型部位正應(yīng)力響應(yīng)功率譜密度與等效前所得結(jié)果吻合較好,所選5個響應(yīng)點(diǎn)處x軸向正應(yīng)力響應(yīng)均方根誤差不超過2%,危險點(diǎn)處x軸向正應(yīng)力響應(yīng)均方根誤差不超過1%,從而驗證了所提載荷等效方法的有效性.

    References)

    [1]孟光, 周徐斌, 苗軍. 航天重大工程中的力學(xué)問題[J]. 力學(xué)進(jìn)展, 2016, 46(1): 267-322. DOI:10.6052/1000-0992-15-018. Meng Guang, Zhou Xubin, Miao Jun. Mechanical problems in momentous projects of aerospace engineering[J].AdvancesInMechanics, 2016, 46(1): 267-322. DOI:10.6052/1000-0992-15-018.(in Chinese)

    [2]Ciappi E, de Rosa S, Franco F, et al. On the dynamic behavior of composite panels under turbulent boundary layer excitations [J].JournalofSoundandVibration, 2016, 364: 77-109. DOI:10.1016/j.jsv.2015.11.024.

    [3]李彥斌, 張鵬, 吳邵慶, 等. 復(fù)合材料加筋板計及熱效應(yīng)的聲-固耦合分析[J]. 振動工程學(xué)報, 2015, 28(4): 531-540. DOI:10.16385/j.cnki.issn.1004-4523.2015.04.005. Li Yanbin, Zhang Peng, Wu Shaoqing, et al. Structural-acoustic coupling analysis of a composite stiffened panel in a thermal environment[J].JournalofVibrationEngineering, 2015, 28(4): 531-540. DOI:10.16385/j.cnki.issn.1004-4523.2015.04.005.(in Chinese)

    [4]Liu L, Guo Q, He T. Thermal-acoustic fatigue of a multilayer thermal protection system in combined extreme environments [J].AdvancesinMechanicalEngineering, 2014, 6: 176891. DOI:10.1155/2014/176891.

    [5]Leung K, Foist B L. Prediction of acoustically induced random vibration loads for shuttle payloads [C]//Proceedingsof36thStructures,StructuralDynamics,andMaterialsConferenceandExhibit. Washington DC, USA,1995: 1200-1208. DOI:10.2514/6.1995-1200.

    [6]Lee H M. Testing for random limit load versus static limit load [R]. Huntsville, Alabama,USA: NASA, 1997.

    [7]馮若強(qiáng), 丁靜鵠, 李慶祥, 等. 肋環(huán)型單層曲板網(wǎng)殼結(jié)構(gòu)的風(fēng)振響應(yīng)及等效靜風(fēng)荷載[J]. 東南大學(xué)學(xué)報(自然科學(xué)版), 2012, 42(2): 328-333. DOI:10.3969/j.issn.1001-0505.2012.02.026. Feng Ruoqiang, Ding Jinghu, Li Qingxiang, et al. Wind-induced response and equivalent static wind load of single-layer geiger lattice shell with curved plate[J].JournalofSoutheastUniversity(NaturalScienceEdition), 2012, 42(2): 328-333. DOI:10.3969/j.issn.1001-0505.2012.02.026.(in Chinese)

    [8]柯世堂, 王浩. 考慮地形效應(yīng)的大型電視塔等效靜風(fēng)荷載分析[J]. 東南大學(xué)學(xué)報(自然科學(xué)版), 2016, 46(3): 545-551. Ke Shitang, Wang Hao. Analysis on equivalent static wind loads of large television tower considering topographic effect [J].JournalofSoutheastUniversity(NaturalScienceEdition), 2016, 46(3): 545-551. (in Chinese)

    [9]Ke S T, Ge Y J, Zhao L, et al. A new methodology for analysis of equivalent static wind loads on super-large cooling towers [J].JournalofWindEngineeringandIndustrialAerodynamics, 2012, 111: 30-39. DOI:10.1016/j.jweia.2012.08.001.

    [10]張艷崗, 毛虎平, 蘇鐵熊, 等. 基于關(guān)鍵時間點(diǎn)的能量等效靜態(tài)載荷法及結(jié)構(gòu)動態(tài)響應(yīng)優(yōu)化[J]. 機(jī)械工程學(xué)報, 2016, 52(9): 151-157. DOI:10.3901/JME.2016.09.151. Zhang Yangang, Mao Huping, Su Tiexiong, et al. Energy equivalent static load method based on critical time point and structure dynamic response optimization[J].JournalofMechanicalEngineering, 2016, 52(9): 151-157. DOI:10.3901/JME.2016.09.151.(in Chinese)

    [11]Rajaneesh A, Satrio W, Chai G B, et al. Long-term life prediction of woven CFRP laminates under three point flexural fatigue[J].CompositesPartB:Engineering, 2016, 91: 539-547. DOI:10.1016/j.compositesb.2016.01.028.

    Equivalence method of pressure load and foundation excitation based on mode superposition technique

    Zhang Peng Fei Qingguo Wu Shaoqing Li Yanbin

    (School of Civil Engineering, Southeast University, Nanjing 210096, China)

    The random pressure load is equivalent to the random foundation excitation based on the modal superposition method to provide evidence for the replacement of the expensive acoustic tests with the shaking table tests for thin-walled structures. The displacement of the leading modes under the two loads are required to be identical. The statically indeterminate problem of the single source load equivalence is properly simplified to a determinate problem. A carbon/carbon (C/C) composite cantilever plate is used as a research subject to investigate the validity of the load equivalence method. The results show that the power spectral density of the random stress responses at typical points of the plate under the random pressure load agrees well with that under the equivalent random foundation excitation. In the analysis of frequency band, the error of the root-mean-square (RMS) of thex-axial normal stress response at the high-risk point on the plate is less than 1%, and the errors of the RMSs of thex-axial normal stress responses at five points distributed on the plate are less than 2%. By using the proposed method, the random pressure load can be effectively equivalent to the foundation excitation, ensuring the consistency of the random stress responses on the thin-walled structures in both the frequency domain and the spacial domain when the structures are subjected to the two loads.

    pressure load; foundation excitation; load equivalence; alternative test; mode superposition method

    10.3969/j.issn.1001-0505.2017.02.022

    2016-08-08. 作者簡介:張鵬(1987—),男,博士生;費(fèi)慶國(聯(lián)系人),男,博士,教授,博士生導(dǎo)師,qgFei@seu.edu.cn.

    國家自然科學(xué)基金資助項目(11402052,11572086)、教育部新世紀(jì)優(yōu)秀人才支持計劃資助項目(NCET-11-0086)、江蘇省自然科學(xué)基金資助項目(BK20140616).

    張鵬,費(fèi)慶國,吳邵慶,等. 基于模態(tài)疊加法的面壓載荷與基礎(chǔ)激勵等效方法[J].東南大學(xué)學(xué)報(自然科學(xué)版),2017,47(2):331-336.

    10.3969/j.issn.1001-0505.2017.02.022.

    TH113

    A

    1001-0505(2017)02-0331-06

    猜你喜歡
    頻響方根軸向
    方根拓展探究
    大型立式單級引黃離心泵軸向力平衡的研究
    基于分塊化頻響函數(shù)曲率比的砌體房屋模型損傷識別研究
    美團(tuán)外賣哥
    荒銑加工軸向切深識別方法
    均方根嵌入式容積粒子PHD 多目標(biāo)跟蹤方法
    頻響函數(shù)殘差法在有限元模型修正中的應(yīng)用
    揭開心算方根之謎
    頻響阻抗法診斷變壓器繞組變形
    微小型薄底零件的軸向車銑實(shí)驗研究
    免费在线观看完整版高清| 午夜精品久久久久久毛片777| 亚洲欧美精品综合久久99| 怎么达到女性高潮| 国产精品99久久99久久久不卡| av中文乱码字幕在线| 哪里可以看免费的av片| 制服诱惑二区| 欧美久久黑人一区二区| 12—13女人毛片做爰片一| 国产私拍福利视频在线观看| 免费人成视频x8x8入口观看| 夜夜夜夜夜久久久久| 国产成人av激情在线播放| 亚洲精品一卡2卡三卡4卡5卡| av中文乱码字幕在线| 色综合站精品国产| av福利片在线| 亚洲av中文字字幕乱码综合| 成人av一区二区三区在线看| 亚洲熟妇中文字幕五十中出| 两性夫妻黄色片| 正在播放国产对白刺激| 一边摸一边做爽爽视频免费| 免费av毛片视频| 午夜a级毛片| 精品免费久久久久久久清纯| 午夜福利高清视频| 一二三四社区在线视频社区8| 一边摸一边抽搐一进一小说| 国产精品永久免费网站| 久久香蕉精品热| 久久婷婷人人爽人人干人人爱| 欧美av亚洲av综合av国产av| 舔av片在线| bbb黄色大片| 脱女人内裤的视频| 国产99白浆流出| 精品乱码久久久久久99久播| 男人的好看免费观看在线视频 | 超碰成人久久| 色噜噜av男人的天堂激情| 国产亚洲精品久久久久5区| 黄色片一级片一级黄色片| 啦啦啦观看免费观看视频高清| 亚洲电影在线观看av| 一进一出好大好爽视频| 欧美成狂野欧美在线观看| 亚洲天堂国产精品一区在线| 亚洲成av人片在线播放无| 亚洲专区国产一区二区| 99热只有精品国产| 一级作爱视频免费观看| 亚洲欧美日韩高清在线视频| 亚洲精品中文字幕在线视频| 久久精品影院6| 成人手机av| 亚洲天堂国产精品一区在线| 日韩中文字幕欧美一区二区| 美女黄网站色视频| 两性夫妻黄色片| 午夜精品一区二区三区免费看| 亚洲 欧美 日韩 在线 免费| 99在线人妻在线中文字幕| 欧美色视频一区免费| 欧美久久黑人一区二区| 久久精品人妻少妇| 欧美日韩乱码在线| 国内精品久久久久久久电影| 中文在线观看免费www的网站 | 日韩欧美国产一区二区入口| 亚洲精品久久成人aⅴ小说| 免费电影在线观看免费观看| 欧美成人性av电影在线观看| 男插女下体视频免费在线播放| 国产区一区二久久| 伦理电影免费视频| 91大片在线观看| 看片在线看免费视频| 午夜视频精品福利| 亚洲狠狠婷婷综合久久图片| 亚洲人与动物交配视频| 熟妇人妻久久中文字幕3abv| 亚洲国产精品sss在线观看| 午夜亚洲福利在线播放| 欧洲精品卡2卡3卡4卡5卡区| 在线观看66精品国产| 男女做爰动态图高潮gif福利片| 老鸭窝网址在线观看| 久久伊人香网站| 国产精品免费视频内射| 亚洲精品美女久久av网站| 久9热在线精品视频| 人妻丰满熟妇av一区二区三区| 12—13女人毛片做爰片一| 国产黄a三级三级三级人| 麻豆国产av国片精品| 怎么达到女性高潮| 国内少妇人妻偷人精品xxx网站 | 国产日本99.免费观看| 国产熟女xx| 久久久国产欧美日韩av| 男女视频在线观看网站免费 | 九色成人免费人妻av| av福利片在线| 91在线观看av| 午夜福利欧美成人| 99riav亚洲国产免费| xxx96com| 天天躁夜夜躁狠狠躁躁| 极品教师在线免费播放| 日韩欧美在线二视频| av福利片在线观看| 中文字幕人妻丝袜一区二区| 女人高潮潮喷娇喘18禁视频| 亚洲av熟女| 蜜桃久久精品国产亚洲av| 激情在线观看视频在线高清| 午夜免费观看网址| 国产成人aa在线观看| 真人做人爱边吃奶动态| 91字幕亚洲| 女生性感内裤真人,穿戴方法视频| 国产成人aa在线观看| 一本综合久久免费| 国产精品久久久人人做人人爽| 成年人黄色毛片网站| 国产探花在线观看一区二区| 亚洲真实伦在线观看| 一夜夜www| 亚洲国产中文字幕在线视频| 精品国产超薄肉色丝袜足j| 国产av不卡久久| 久久久精品国产亚洲av高清涩受| 日本熟妇午夜| 欧美黑人欧美精品刺激| xxxwww97欧美| 国产精品一区二区精品视频观看| 国产av一区二区精品久久| 久久久久九九精品影院| 中文字幕最新亚洲高清| 日韩欧美国产一区二区入口| 夜夜爽天天搞| 日日摸夜夜添夜夜添小说| 午夜免费成人在线视频| 午夜福利在线观看吧| 亚洲在线自拍视频| 亚洲精品在线观看二区| 久久天躁狠狠躁夜夜2o2o| 久久热在线av| 国产亚洲精品第一综合不卡| 老熟妇乱子伦视频在线观看| 欧美黄色片欧美黄色片| 国产成人一区二区三区免费视频网站| 两性午夜刺激爽爽歪歪视频在线观看 | 黄片大片在线免费观看| 在线观看舔阴道视频| 国产又色又爽无遮挡免费看| 国产精华一区二区三区| 国产伦在线观看视频一区| 在线观看一区二区三区| 操出白浆在线播放| 岛国在线观看网站| 老汉色av国产亚洲站长工具| 国产成人欧美在线观看| 国产视频一区二区在线看| 亚洲熟妇熟女久久| 女人被狂操c到高潮| 777久久人妻少妇嫩草av网站| 日韩欧美三级三区| 久久久国产成人精品二区| 精品福利观看| 九九热线精品视视频播放| 精品欧美一区二区三区在线| 欧美成人午夜精品| 国产黄色小视频在线观看| 欧洲精品卡2卡3卡4卡5卡区| www.自偷自拍.com| 中国美女看黄片| 两个人看的免费小视频| 亚洲自拍偷在线| 亚洲精华国产精华精| 国内精品久久久久精免费| 亚洲乱码一区二区免费版| 一级片免费观看大全| 人妻丰满熟妇av一区二区三区| 亚洲乱码一区二区免费版| 久久精品91无色码中文字幕| 精品高清国产在线一区| 美女黄网站色视频| 精品日产1卡2卡| 身体一侧抽搐| 欧美日韩亚洲综合一区二区三区_| 午夜成年电影在线免费观看| 好男人电影高清在线观看| 亚洲七黄色美女视频| 国产又黄又爽又无遮挡在线| 女人高潮潮喷娇喘18禁视频| 给我免费播放毛片高清在线观看| 国产午夜精品论理片| 国产av一区二区精品久久| 亚洲成人久久爱视频| 曰老女人黄片| 丝袜美腿诱惑在线| 国内精品久久久久精免费| 亚洲国产精品sss在线观看| 国产精华一区二区三区| 亚洲av五月六月丁香网| 亚洲av电影在线进入| 亚洲avbb在线观看| 久久欧美精品欧美久久欧美| 国产精品爽爽va在线观看网站| 99精品久久久久人妻精品| 好看av亚洲va欧美ⅴa在| 久久精品人妻少妇| 欧美zozozo另类| 人妻丰满熟妇av一区二区三区| 18禁黄网站禁片免费观看直播| 嫁个100分男人电影在线观看| 欧美高清成人免费视频www| 久久精品夜夜夜夜夜久久蜜豆 | 热99re8久久精品国产| 日韩精品中文字幕看吧| 久久香蕉国产精品| 日本撒尿小便嘘嘘汇集6| 欧美最黄视频在线播放免费| 熟女少妇亚洲综合色aaa.| 久久久精品欧美日韩精品| 欧美成人一区二区免费高清观看 | 男人的好看免费观看在线视频 | 国产私拍福利视频在线观看| tocl精华| 看免费av毛片| 精品人妻1区二区| 日韩av在线大香蕉| 91老司机精品| 久久久国产欧美日韩av| 亚洲七黄色美女视频| 欧美黄色淫秽网站| 亚洲av五月六月丁香网| 欧美黄色片欧美黄色片| 无人区码免费观看不卡| 超碰成人久久| 1024香蕉在线观看| 美女 人体艺术 gogo| 少妇熟女aⅴ在线视频| 男女床上黄色一级片免费看| 三级毛片av免费| 国产亚洲精品综合一区在线观看 | 国产成人一区二区三区免费视频网站| 一本综合久久免费| 久久精品aⅴ一区二区三区四区| 香蕉丝袜av| 精品欧美一区二区三区在线| 精品国产亚洲在线| videosex国产| 老司机午夜福利在线观看视频| av免费在线观看网站| 深夜精品福利| 在线看三级毛片| 欧美不卡视频在线免费观看 | 成人18禁高潮啪啪吃奶动态图| 999精品在线视频| 最近视频中文字幕2019在线8| 狂野欧美白嫩少妇大欣赏| 久久久久精品国产欧美久久久| 国产精品久久久久久久电影 | 少妇熟女aⅴ在线视频| 男女床上黄色一级片免费看| 久久久精品国产亚洲av高清涩受| 欧美乱码精品一区二区三区| 无遮挡黄片免费观看| 免费搜索国产男女视频| 亚洲人成网站高清观看| 伊人久久大香线蕉亚洲五| 久久亚洲真实| 熟女少妇亚洲综合色aaa.| 国产精品综合久久久久久久免费| 亚洲精品国产精品久久久不卡| 精品午夜福利视频在线观看一区| 久久香蕉精品热| 亚洲无线在线观看| 久久草成人影院| av免费在线观看网站| 床上黄色一级片| 国产精品亚洲av一区麻豆| 欧美日韩瑟瑟在线播放| 免费在线观看完整版高清| 制服诱惑二区| 久久久久免费精品人妻一区二区| 最近在线观看免费完整版| 国产精品 欧美亚洲| 三级男女做爰猛烈吃奶摸视频| 欧美在线黄色| 日本a在线网址| 国产激情久久老熟女| 757午夜福利合集在线观看| 国产成人系列免费观看| 性色av乱码一区二区三区2| 老司机午夜十八禁免费视频| 亚洲自偷自拍图片 自拍| 婷婷精品国产亚洲av在线| АⅤ资源中文在线天堂| tocl精华| 50天的宝宝边吃奶边哭怎么回事| 黄色 视频免费看| АⅤ资源中文在线天堂| 国产亚洲欧美98| 麻豆成人av在线观看| 露出奶头的视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲一区二区三区色噜噜| 色在线成人网| 丰满人妻一区二区三区视频av | 99国产精品99久久久久| ponron亚洲| 最近视频中文字幕2019在线8| 久久精品国产99精品国产亚洲性色| 婷婷精品国产亚洲av| av有码第一页| 久久久久国产一级毛片高清牌| 两个人的视频大全免费| 69av精品久久久久久| 成人一区二区视频在线观看| 免费av毛片视频| 欧美丝袜亚洲另类 | 国产日本99.免费观看| 在线观看66精品国产| 淫秽高清视频在线观看| www.熟女人妻精品国产| 天天躁夜夜躁狠狠躁躁| 我要搜黄色片| 伦理电影免费视频| 成人国产一区最新在线观看| 91av网站免费观看| 级片在线观看| 欧美日韩黄片免| 欧美日韩乱码在线| 久久久国产成人精品二区| 久久精品91蜜桃| 中文字幕熟女人妻在线| 法律面前人人平等表现在哪些方面| 午夜视频精品福利| 最近最新中文字幕大全免费视频| 久久人妻av系列| 国产av在哪里看| 久久人人精品亚洲av| 一二三四社区在线视频社区8| 国产单亲对白刺激| 老汉色∧v一级毛片| 狠狠狠狠99中文字幕| 成人一区二区视频在线观看| 欧美av亚洲av综合av国产av| 日本 欧美在线| 丰满人妻熟妇乱又伦精品不卡| 熟女电影av网| 床上黄色一级片| 久久香蕉激情| netflix在线观看网站| 18禁观看日本| 欧美黑人欧美精品刺激| 成人av在线播放网站| 在线观看舔阴道视频| 国产主播在线观看一区二区| 日本一二三区视频观看| 两个人看的免费小视频| 999久久久国产精品视频| 99在线人妻在线中文字幕| 成年免费大片在线观看| 国产精品免费视频内射| 国产精品九九99| 色av中文字幕| 国内少妇人妻偷人精品xxx网站 | www.熟女人妻精品国产| 美女黄网站色视频| 免费观看精品视频网站| 黄色视频不卡| 99精品在免费线老司机午夜| 露出奶头的视频| 亚洲av片天天在线观看| av在线天堂中文字幕| 黑人操中国人逼视频| a级毛片a级免费在线| 99久久99久久久精品蜜桃| 嫩草影院精品99| 中文资源天堂在线| 欧美日韩乱码在线| 亚洲av五月六月丁香网| 99riav亚洲国产免费| 日本黄色视频三级网站网址| 成人av在线播放网站| 日本免费一区二区三区高清不卡| 亚洲人成伊人成综合网2020| 午夜福利18| 国产精品亚洲av一区麻豆| 搡老岳熟女国产| 成年女人毛片免费观看观看9| 精品久久久久久,| 久久久久国产一级毛片高清牌| 国产又色又爽无遮挡免费看| 日韩欧美国产在线观看| 日本精品一区二区三区蜜桃| 国产激情欧美一区二区| 90打野战视频偷拍视频| 不卡av一区二区三区| а√天堂www在线а√下载| 久久久久国产精品人妻aⅴ院| 亚洲一码二码三码区别大吗| 婷婷精品国产亚洲av| 18禁裸乳无遮挡免费网站照片| 国产亚洲av高清不卡| 亚洲五月婷婷丁香| 日本免费一区二区三区高清不卡| 亚洲免费av在线视频| 色在线成人网| 国产久久久一区二区三区| 91麻豆精品激情在线观看国产| 日本一本二区三区精品| 九色成人免费人妻av| 久久久久久久久中文| 舔av片在线| 天天一区二区日本电影三级| 丝袜人妻中文字幕| 国产高清视频在线播放一区| 桃色一区二区三区在线观看| 中文字幕人成人乱码亚洲影| 成人18禁高潮啪啪吃奶动态图| 亚洲av成人av| 人人妻人人看人人澡| 制服丝袜大香蕉在线| www日本在线高清视频| 最近最新免费中文字幕在线| 精品国产超薄肉色丝袜足j| 成人av在线播放网站| 人人妻,人人澡人人爽秒播| 久久草成人影院| 国产精品av视频在线免费观看| 久久午夜综合久久蜜桃| 亚洲精品美女久久av网站| 亚洲av片天天在线观看| 大型av网站在线播放| 欧美精品亚洲一区二区| 成人手机av| 亚洲专区中文字幕在线| 欧美丝袜亚洲另类 | 午夜视频精品福利| 成年免费大片在线观看| 日韩 欧美 亚洲 中文字幕| 精品久久久久久久久久久久久| 天堂动漫精品| 国产熟女xx| 亚洲七黄色美女视频| 久热爱精品视频在线9| 国产av又大| 亚洲精品在线观看二区| 亚洲男人的天堂狠狠| 一进一出好大好爽视频| 中文字幕最新亚洲高清| 国产精品日韩av在线免费观看| 中文资源天堂在线| 人成视频在线观看免费观看| 亚洲熟女毛片儿| 操出白浆在线播放| 亚洲成人中文字幕在线播放| 日本免费一区二区三区高清不卡| 国产亚洲欧美98| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲美女视频黄频| 黑人欧美特级aaaaaa片| 国产精品自产拍在线观看55亚洲| 色老头精品视频在线观看| 亚洲五月婷婷丁香| 国产精品影院久久| 蜜桃久久精品国产亚洲av| 国产精品免费一区二区三区在线| 99国产精品99久久久久| 日韩av在线大香蕉| 亚洲精品中文字幕在线视频| 国产69精品久久久久777片 | 男人的好看免费观看在线视频 | 两人在一起打扑克的视频| 亚洲专区中文字幕在线| 免费无遮挡裸体视频| 亚洲狠狠婷婷综合久久图片| 精品久久久久久久人妻蜜臀av| 亚洲18禁久久av| 日本在线视频免费播放| 一级a爱片免费观看的视频| 国产视频内射| av天堂在线播放| 三级毛片av免费| 欧美+亚洲+日韩+国产| 日本一区二区免费在线视频| ponron亚洲| 欧美日韩福利视频一区二区| 国产99白浆流出| 欧美乱妇无乱码| 亚洲一区二区三区色噜噜| 亚洲成av人片免费观看| 精品国产超薄肉色丝袜足j| 午夜福利视频1000在线观看| 黄色成人免费大全| 最近最新免费中文字幕在线| 在线a可以看的网站| 最近最新免费中文字幕在线| 亚洲精品一卡2卡三卡4卡5卡| av超薄肉色丝袜交足视频| 搞女人的毛片| 日本 欧美在线| 亚洲专区中文字幕在线| 亚洲国产中文字幕在线视频| 老汉色av国产亚洲站长工具| 成人高潮视频无遮挡免费网站| 国产精华一区二区三区| 亚洲乱码一区二区免费版| 精品久久久久久久久久久久久| 琪琪午夜伦伦电影理论片6080| 婷婷亚洲欧美| 操出白浆在线播放| 日韩精品免费视频一区二区三区| 日日爽夜夜爽网站| 狠狠狠狠99中文字幕| 999久久久精品免费观看国产| 欧美丝袜亚洲另类 | 黄片大片在线免费观看| 亚洲真实伦在线观看| 精品久久久久久久久久久久久| 91成年电影在线观看| 国产精品99久久99久久久不卡| 国产熟女午夜一区二区三区| 美女大奶头视频| 两个人看的免费小视频| 在线观看美女被高潮喷水网站 | 夜夜爽天天搞| 少妇粗大呻吟视频| 99热这里只有精品一区 | 制服丝袜大香蕉在线| 久久精品91蜜桃| 久久久国产成人精品二区| 日韩欧美三级三区| 国产免费男女视频| 久久午夜亚洲精品久久| 亚洲国产精品sss在线观看| 观看免费一级毛片| 免费观看人在逋| 日本黄大片高清| 婷婷丁香在线五月| 国产午夜精品久久久久久| 老熟妇乱子伦视频在线观看| 精品熟女少妇八av免费久了| 久久精品aⅴ一区二区三区四区| 99国产精品99久久久久| 18美女黄网站色大片免费观看| 人人妻,人人澡人人爽秒播| 午夜久久久久精精品| 中文字幕熟女人妻在线| 一区二区三区激情视频| 国产区一区二久久| 日本精品一区二区三区蜜桃| 在线十欧美十亚洲十日本专区| 国产真实乱freesex| 国产一区在线观看成人免费| 国产1区2区3区精品| 精品久久久久久久末码| 女人高潮潮喷娇喘18禁视频| 精品久久久久久久久久免费视频| 欧美成人性av电影在线观看| 999久久久国产精品视频| 丰满人妻熟妇乱又伦精品不卡| 听说在线观看完整版免费高清| 欧美成人午夜精品| 久久久国产成人精品二区| 国产欧美日韩精品亚洲av| 老熟妇乱子伦视频在线观看| 午夜激情福利司机影院| 国产精品免费一区二区三区在线| 亚洲自偷自拍图片 自拍| 久久香蕉激情| 99久久精品热视频| 香蕉久久夜色| 国产欧美日韩一区二区精品| 国产主播在线观看一区二区| 亚洲国产日韩欧美精品在线观看 | 一级片免费观看大全| 午夜福利在线在线| 久久久久久久久免费视频了| cao死你这个sao货| 亚洲全国av大片| 欧美大码av| 亚洲七黄色美女视频| 亚洲专区中文字幕在线| 日韩欧美国产在线观看| 一夜夜www| 一区二区三区高清视频在线| 法律面前人人平等表现在哪些方面| 最近视频中文字幕2019在线8| 国产一级毛片七仙女欲春2| 婷婷丁香在线五月| 视频区欧美日本亚洲| 成人特级黄色片久久久久久久| 精品午夜福利视频在线观看一区| 欧美性猛交黑人性爽| 在线观看免费视频日本深夜| 久久久久国产一级毛片高清牌| 色综合站精品国产| 成人特级黄色片久久久久久久| 色av中文字幕| 嫩草影院精品99| 人妻夜夜爽99麻豆av| 精品久久久久久成人av| 听说在线观看完整版免费高清| 免费在线观看影片大全网站| 在线播放国产精品三级| 一夜夜www| 神马国产精品三级电影在线观看 |