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

        ?

        基于實(shí)時混合試驗(yàn)誤差累積規(guī)律的時滯影響修正方法

        2021-07-06 07:01:38周子豪李忠獻(xiàn)
        工程力學(xué) 2021年6期
        關(guān)鍵詞:物理模型

        周子豪,李 寧,2,李忠獻(xiàn),2

        (1.天津大學(xué)建筑工程學(xué)院/濱海土木工程結(jié)構(gòu)與安全教育部重點(diǎn)實(shí)驗(yàn)室(天津大學(xué)),天津 300350;2.中國地震局地震工程綜合模擬與城鄉(xiāng)抗震韌性重點(diǎn)實(shí)驗(yàn)室(天津大學(xué)),天津 300350)

        1969年日本東京大學(xué)Motohiko課題組實(shí)現(xiàn)了第一個混合試驗(yàn)[1],開啟了混合試驗(yàn)方法的研究歷程。經(jīng)過30余年發(fā)展,混合試驗(yàn)結(jié)構(gòu)劃分的試驗(yàn)理念得到了廣泛認(rèn)同:將試驗(yàn)對象中本構(gòu)關(guān)系明確、易于模擬的部分劃分為數(shù)值子結(jié)構(gòu),而試驗(yàn)對象中難于建模卻又對整體結(jié)構(gòu)的響應(yīng)起到重要作用的部分劃分為物理子結(jié)構(gòu),并在試驗(yàn)室中進(jìn)行試驗(yàn)以獲得其實(shí)際的響應(yīng)。從擬動力混合試驗(yàn)到快速擬動力混合試驗(yàn)再到實(shí)時混合試驗(yàn),混合試驗(yàn)的應(yīng)用愈發(fā)成熟。近十年來,涌現(xiàn)了許多實(shí)時混合試驗(yàn)成果。Karavasilis等[2]利用實(shí)時混合試驗(yàn)對基于性能抗震設(shè)計的裝有阻尼器的鋼框架進(jìn)行了性能評估測試。Chen等[3]利用實(shí)時混合試驗(yàn)對安裝大型磁流變阻尼器的鋼框架進(jìn)行了分析。Shao等[4]對大型結(jié)構(gòu)的振動臺和作動器聯(lián)合實(shí)時混合試驗(yàn)方法進(jìn)行了研究。Sorkhabi等[5]利用實(shí)時混合試驗(yàn)對TLD在地震作用下的性能進(jìn)行了研究。Chen等[6]研發(fā)了實(shí)時混合試驗(yàn)平臺,考慮了作動器、控制架構(gòu)、數(shù)值積分算法以及時滯補(bǔ)償?shù)龋Υ笮蛯?shí)時混合試驗(yàn)的多作動器時滯自適應(yīng)補(bǔ)償[7]進(jìn)行了探索。Chae等[8]利用實(shí)時混合試驗(yàn)對大型MR阻尼器控制策略進(jìn)行了研究。Xu等[9]利用實(shí)時混合試驗(yàn)對結(jié)構(gòu)主動控制性能進(jìn)行了研究。Eem等[10]對MRE智能基礎(chǔ)隔振的性能進(jìn)行了研究和評估。Asai等[11]利用實(shí)時混合試驗(yàn)對高聳結(jié)構(gòu)智能外伸阻尼系統(tǒng)進(jìn)行了研究。Saouma等[12]進(jìn)了非延性混凝土結(jié)構(gòu)的實(shí)時混合試驗(yàn),并將實(shí)時混合試驗(yàn)的結(jié)果與振動臺試驗(yàn)進(jìn)行了對比。Gao等[13]對多自由度實(shí)時混合試驗(yàn)的一般流程進(jìn)行了研究和試驗(yàn)驗(yàn)證。Cha等[14]對比了多種MR阻尼器半主動控制實(shí)時混合試驗(yàn)結(jié)果和純模擬結(jié)果。Jiang等[15]針對復(fù)雜數(shù)值子結(jié)構(gòu)的情況開展了安裝MR阻尼器的橋梁實(shí)時混合試驗(yàn)。

        目前實(shí)時混合試驗(yàn)幾乎都是在少數(shù)幾個設(shè)備精良的高校實(shí)驗(yàn)室中完成的。這些試驗(yàn)室中環(huán)境干擾少,試驗(yàn)人員專業(yè)知識和試驗(yàn)經(jīng)驗(yàn)較豐富,能夠較好地對試驗(yàn)的過程進(jìn)行把控,對試驗(yàn)中的不足可以及時發(fā)現(xiàn)和修正。還有一個重要的原因是,由于實(shí)時混合試驗(yàn)需保證試驗(yàn)過程的實(shí)時性,對實(shí)時機(jī)、控制器和作動器都有較高的要求,造成試驗(yàn)設(shè)備的成本高,這在一定程度上也阻礙了實(shí)時混合方法的發(fā)展和應(yīng)用[16?17]。此外,隨著作動器最大負(fù)載的增加,要達(dá)到實(shí)時的難度也隨之增加。目前大部分實(shí)時混合實(shí)驗(yàn)中物理子結(jié)構(gòu)依然是縮尺結(jié)構(gòu)模型。對大型的結(jié)構(gòu)而言,如橋墩、高層建筑等,體系并不復(fù)雜,僅僅是體形較大,受限于加載設(shè)備能力,致使實(shí)時混合試驗(yàn)難于開展原型試驗(yàn)。近幾年,隨著實(shí)時子結(jié)構(gòu)試驗(yàn)技術(shù)和設(shè)備的發(fā)展,國內(nèi)外也出現(xiàn)了許多針對大型結(jié)構(gòu)的實(shí)時混合試驗(yàn),例如,Liu等[18]基于實(shí)時混合試驗(yàn)探究了鋼筋混凝土薄壁高墩橋在地震作用下的結(jié)構(gòu)響應(yīng)以及提出了該結(jié)構(gòu)失效模型預(yù)測方法;Wu等[19]進(jìn)行了足尺鋼框架的界面模型更新的實(shí)時混合試驗(yàn);Mei等[20]進(jìn)了高墩混凝土橋梁的基于模型更新的實(shí)時混合試驗(yàn)。

        為了保證實(shí)時混合試驗(yàn)的實(shí)時性,問題主要集中在如何有效地對試驗(yàn)中的時滯影響進(jìn)行消除。在實(shí)時混合試驗(yàn)中,相對于數(shù)值計算和數(shù)據(jù)傳輸?shù)臅r滯,作動器時滯是占主導(dǎo)地位的。為了消除時滯的影響,國內(nèi)外的學(xué)者進(jìn)行了非常多的探索。其中代表性的有Horiuchi等[21]提出的多項(xiàng)式外插時滯補(bǔ)償方法、Chen[22]提出的一階反補(bǔ)償時滯補(bǔ)償算法以及Phillips等[23]提出的基于模型的時滯補(bǔ)償方法。筆者所在的科研團(tuán)隊也在時滯補(bǔ)償算法的開發(fā)上做出了自己的貢獻(xiàn)。例如,將在線時滯預(yù)估[24]分別與多項(xiàng)式外插時滯補(bǔ)償[21]和一階反補(bǔ)償[22]相結(jié)合,提出了二階在線迭代時滯補(bǔ)償算法[25]以及自適應(yīng)時滯補(bǔ)償算法[26];引入模型誤差提出了二階的基于模型的時滯補(bǔ)償方法[27]。但是,目前時滯補(bǔ)償算法依舊無法消除時滯,特別是對于高頻信號,補(bǔ)償?shù)男Ч匀徊焕硐搿?/p>

        考慮到時滯普遍存在于實(shí)時混合試驗(yàn)中,針對由于實(shí)驗(yàn)設(shè)備性能限制或結(jié)構(gòu)負(fù)載較大而導(dǎo)致的時滯補(bǔ)償不理想等情況,本文提出了一種基于雙顯式數(shù)值積分算法誤差累積規(guī)律的試驗(yàn)后時滯修正方法??稍谠囼?yàn)后對結(jié)果進(jìn)行修正以消除時滯影響,降低試驗(yàn)成本。首先,對雙顯示數(shù)值積分算法的誤差累積規(guī)律進(jìn)行了探索,指出在實(shí)時混合試驗(yàn)過程中,時滯導(dǎo)致的誤差是按照一定的規(guī)律進(jìn)行累積傳播的;其次,根據(jù)得到的誤差累積規(guī)律提出了一種試驗(yàn)后的時滯補(bǔ)償方法,在試驗(yàn)后對試驗(yàn)結(jié)果進(jìn)行修正,以達(dá)到減小時滯影響的目的;最后對方法的有效性通過4種不同的物理子結(jié)構(gòu)工況進(jìn)行了驗(yàn)證。

        1 數(shù)值積分算法的誤差累積規(guī)律

        數(shù)值積分算法在實(shí)時混合試驗(yàn)中是關(guān)鍵部分。物理子結(jié)構(gòu)的反饋進(jìn)入數(shù)值子結(jié)構(gòu)之后,通過數(shù)值積分算法得到下一步的位移指令,進(jìn)而發(fā)送到控制器驅(qū)動作動器運(yùn)動。由于時滯的存在,在試驗(yàn)過程中作動器的位置相對于位移指令滯后,導(dǎo)致物理子結(jié)構(gòu)反饋滯后,此時積分得到的結(jié)果也就必然與理想無時滯狀態(tài)存在誤差。因此,無時滯理想狀態(tài)下與有時滯狀態(tài)下數(shù)值子結(jié)構(gòu)計算的結(jié)果是不同的。無時滯、有時滯單自由度體系運(yùn)動方程可寫為式(1)和式(2)。

        雙顯式數(shù)值積分算法的通用表達(dá)式為式(4)、式(5)所示,速度和位移均為顯式遞推得出。

        當(dāng)數(shù)值積分算法分別在無時滯理想情況和有時滯情況下,可分別改寫為式(6)、式(7)及式(8)、式(9)。

        式(8)減去式(6),式(9)減去式(7),得式(10)、式(11)。

        式(10)、式(11)即為雙顯式數(shù)值積分的位移和速度誤差累積公式,其中:Ed、Ev和Ea分別表示位移誤差、速度誤差和加速度誤差??梢钥吹?,其表達(dá)式與積分算法的形式是一致的,而且位移和速度的誤差累積公式是顯式的。將運(yùn)動方程式(2)減去式(1),可以得到式(12)、式(13)。

        綜上所述,式(14)、式(15)及式(16)組成的方程組:

        綜合上述的推導(dǎo),將整體的修正過程總結(jié)如圖1所示。

        圖1 試驗(yàn)后修正方法框圖Fig.1 Flow chart of post-correction method

        Gao等[28]對時滯試驗(yàn)結(jié)果的影響進(jìn)行了分析,其時滯微分方程列式如下:

        令x=Asinωt,忽略Gao等[28]引入的幅值誤差,只考慮時滯,則xm=Asinω(t?δt)。其中:M、C、K分別為質(zhì)量、阻尼和剛度;下角標(biāo)N、P分別代表數(shù)值和物理子結(jié)構(gòu);x¨g為地震動加速度;xm為作動器量測位移;A和ω分別為結(jié)構(gòu)振動的幅值和頻率;δt為時滯。假設(shè)δt很小,將xm進(jìn)行泰勒展開得式(22)~式(23):

        將其代入時滯微分方程得:

        觀察式(25),發(fā)現(xiàn)物理子結(jié)構(gòu)的剛度向運(yùn)動方程中引入了負(fù)阻尼CKp=?KPδt,向系統(tǒng)中引入了能量,這與Horiuchi等[21]得到的結(jié)論一致。從式(25)還可以發(fā)現(xiàn),物理子結(jié)構(gòu)的質(zhì)量向運(yùn)動方程中引入了正阻尼CMp=MPω2δt,耗散了系統(tǒng)的能量;物理子結(jié)構(gòu)的阻尼向運(yùn)動方程中引入了負(fù)質(zhì)量MCp=?CPδt,向系統(tǒng)中引入了能量。

        因此,當(dāng)存在時滯時,物理子結(jié)構(gòu)的質(zhì)量會增加系統(tǒng)的能耗、降低系統(tǒng)的能量;而物理子結(jié)構(gòu)的剛度和阻尼則會向系統(tǒng)中引入能量。故此,當(dāng)物理子結(jié)構(gòu)質(zhì)量耗散的能量小于剛度和阻尼引入的能量時,實(shí)時混合試驗(yàn)的結(jié)構(gòu)響應(yīng)要大于無時滯的結(jié)構(gòu)響應(yīng)。因此,在物理子結(jié)構(gòu)的質(zhì)量相對較小的情況下,可以在有時滯的實(shí)時混合試驗(yàn)結(jié)果中找到無時滯條件下的結(jié)構(gòu)響應(yīng)。

        至此得到應(yīng)用本文提出方法的第一個限制條件:物理子結(jié)構(gòu)質(zhì)量配比相對較小,并使得物理子結(jié)構(gòu)質(zhì)量耗散的能量小于物理子結(jié)構(gòu)剛度和阻尼引入的能量。這一點(diǎn)在實(shí)際的應(yīng)用中很容易滿足,其原因在于質(zhì)量帶來的慣性力可以比較準(zhǔn)確的進(jìn)行測量,且質(zhì)量在試驗(yàn)過程中一般是不變的,其物理規(guī)律是比較清晰的。故在實(shí)時混合試驗(yàn)中一般不會將質(zhì)量占總結(jié)構(gòu)比例非常大的構(gòu)件作為物理子結(jié)構(gòu)[29]。

        圖2 搜索方式示意圖Fig.2 Schematic diagram of search method

        3 模擬驗(yàn)證

        本文采用4種典型工況對所提出方法進(jìn)行模擬驗(yàn)證:物理子結(jié)構(gòu)分別為線性剛度、線性阻尼、非線性剛度及非線性阻尼構(gòu)件。其中非線性剛度構(gòu)件采用的是非線性彈簧[30]以及Bouc-Wen模型[31],而非線性阻尼構(gòu)件則是采用MR阻尼器模型[32]。

        首先,分析了時滯對試驗(yàn)結(jié)果的影響。對于單自由度體系,取數(shù)值子結(jié)構(gòu)的MN、CN、KN分別為7 kg、36 (N·s)/m、100 N/m;物理子結(jié)構(gòu)的MP、CP、KP分別為3 kg、4(N·s)/m、900 N/m;引入如圖3所示的時滯量,根據(jù)目前作動器的性能,其時滯數(shù)值的取值在15/1024 s~25/1024 s波動;地震激勵為El-Centro波。有時滯與無時滯的模擬結(jié)果對比如圖4所示。有時滯情況得到的試驗(yàn)結(jié)果與無時滯情況得到的試驗(yàn)結(jié)果誤差如圖5所示。

        圖3 人為引入時滯Fig.3 Artificially introduced time delay signal

        圖4 有/無時滯試驗(yàn)結(jié)果對比Fig.4 Test results with/without time delay

        圖5 有時滯與無時滯試驗(yàn)結(jié)果誤差Fig.5 Errorsbetween testing resultswith/without delay

        3.1 線性剛度、阻尼構(gòu)件

        表1 數(shù)值與物理子結(jié)構(gòu)參數(shù)Table 1 Parametersof numerical and physical substructures

        從圖6~圖7可以看出,在時滯相同、物理子結(jié)構(gòu)為線性剛度的情況下,隨著物理子結(jié)構(gòu)剛度占總體剛度比例的增加,有時滯與無時滯試驗(yàn)結(jié)果的位移、速度和加速度誤差呈近似指數(shù)增長。且通過上述的模擬結(jié)果可以發(fā)現(xiàn),在相同的時滯條件下,隨著物理子結(jié)構(gòu)剛度占比的增加,有時滯與無時滯情況下的準(zhǔn)確值相差越來越大。上述模擬的結(jié)果均為理想模擬得到,在實(shí)際的試驗(yàn)過程中,由于環(huán)境噪聲、系統(tǒng)擾動等因素的干擾下誤差勢必會進(jìn)一步增加。

        圖7 線性阻尼物理子結(jié)構(gòu)修正前后誤差對比Fig.7 RMScomparison for linear damping physical substructure tests before and after correction

        當(dāng)物理子結(jié)構(gòu)為線性阻尼的情況下,隨著物理子結(jié)構(gòu)阻尼占總體阻尼比例的增加,有時滯與無時滯的試驗(yàn)結(jié)果表明,位移、速度和加速度誤差同物理子結(jié)構(gòu)為線性剛度時都呈現(xiàn)出增大趨勢。當(dāng)物理子結(jié)構(gòu)為線性剛度時,物理子結(jié)構(gòu)的剛度占總體結(jié)構(gòu)的比例為20%時,其RMS值已經(jīng)達(dá)到了400%,試驗(yàn)結(jié)果已不具備任何參考價值;而當(dāng)物理子結(jié)構(gòu)為線性阻尼時,物理子結(jié)構(gòu)的阻尼占總體結(jié)構(gòu)的比例為5%~80%,但是其RMS卻遠(yuǎn)遠(yuǎn)小于物理子結(jié)構(gòu)為線性剛度的情況。而且,物理子結(jié)構(gòu)為線性阻尼時,其占總體阻尼比例的增加對位移、速度和加速度的影響是依次遞增的,并不象物理子結(jié)構(gòu)為線性剛度時,物理子結(jié)構(gòu)剛度占比的增加對位移、速度和加速度誤差增速的影響是基本一致的。

        3.2 非線性剛度構(gòu)件

        對于非線性的剛度物理子結(jié)構(gòu),我們采用了兩種不同的模型:1)高階非線性彈簧模型[30],如式(27)所示;2) Bouc-Wen模型[31]如式(28)、式(29)所示。模型中物理子結(jié)構(gòu)和數(shù)值子結(jié)構(gòu)參數(shù)的取值見表2及表3。高階彈簧及Bouc-Wen模型的位移-力曲線如圖8所示,其中Bouc-Wen模型采用6種不同幅值的正弦激勵,A代表正弦激勵的幅值。

        圖8 高階彈簧及Bouc-Wen模型特性Fig.8 Characteristicsof high-order spring and Bouc-Wen model

        表2高階彈簧及Bouc-Wen模型參數(shù)Table2 Parametersof high-order spring and Bouc-Wen model

        表3數(shù)值子結(jié)構(gòu)參數(shù)Table 3 Parametersof numerical substructure

        為了驗(yàn)證算法的普適性,以及保證實(shí)時混合試驗(yàn)的穩(wěn)定性,本文將物理子結(jié)構(gòu)的反饋力乘以一個系數(shù)之后再反饋到數(shù)值子結(jié)構(gòu)中,分別稱為高階彈簧系數(shù)及Bouc-Wen系數(shù)。隨著系數(shù)不斷增大,物理子結(jié)構(gòu)的剛度占總體剛度的比例不斷增加。圖9為有時滯情況下高階彈簧及Bouc-Wen模型修正前后的誤差(RMS值)對比。通過圖9可以發(fā)現(xiàn),類似于線性剛度情況,時滯對于試驗(yàn)結(jié)果的影響隨著物理子結(jié)構(gòu)的剛度占比增加而增加。而且,本文方法在物理子結(jié)構(gòu)為高階彈簧時表現(xiàn)出了與純剛度情況下一樣的準(zhǔn)確性,這是因?yàn)槲锢碜咏Y(jié)構(gòu)為高階彈簧時,根據(jù)位移就可以準(zhǔn)確地對應(yīng)其反饋力。其修正前、后的位移、速度和加速度誤差對比如圖10~圖12所示。

        圖9 高階彈簧及Bouc-Wen模型修正前后誤差對比Fig.9 Error comparison of high-order spring and Bouc-Wen model before and after correction

        圖10 修正前后誤差對比位移時程圖(系數(shù)0.8)Fig.10 Displacement comparison before and after correction(physical substructure:Bouc-Wen model,coefficient:0.8)

        圖12 修正前后誤差對比加速度時程圖(系數(shù)0.8)Fig.12 Acceleration comparison beforeand after correction(physical substructure:Bouc-Wen model,coefficient:0.8)

        圖11 修正前后誤差對比速度時程(系數(shù)0.8)Fig.11 Velocity comparison before and after correction(physical substructure:Bouc-Wen model,coefficient:0.8)

        但物理子結(jié)構(gòu)為Bouc-Wen模型時,本文所提方法的誤差卻隨著系數(shù)的增加而不斷的增加。這主要是由于:1)Bouc-Wen模型的反饋力是由位移和速度共同決定,本文提出方法的適用條件之一是位移相關(guān)力或者速度相關(guān)力占據(jù)反饋力的主導(dǎo)地位,而Bouc-Wen模型反饋力由位移和速度共同決定的特性決定了根據(jù)位移對其反饋力進(jìn)行確定會造成誤差;2)Bouc-Wen模型的滯回特性是造成誤差的另一個重要原因,本文提出的方法是根據(jù)位移和速度其中之一進(jìn)行反饋力的修正,然而在位移曲線的拐點(diǎn)處,Bouc-Wen模型由于滯回的特性,即使位移大小一致,但在拐點(diǎn)前后所對應(yīng)的反饋力是不一致的,也會造成修正反饋力的誤差。但還是可以非常明顯地看到本文提出的算法對試驗(yàn)結(jié)果的修正效果尚可。

        3.3 非線性阻尼構(gòu)件

        非線性阻尼構(gòu)建的模型采用的是Spencer提出的磁流變阻尼器模型[32],其表達(dá)式如式(30)~式(32)所示,其模型如圖13所示,其在正弦激勵下的響應(yīng)如圖14所示,其中I為輸入磁流變阻尼器的電流。

        圖13 磁流變阻尼器模型Fig.13 Magnetorheological damper model

        圖14 磁流變阻尼器性能Fig.14 Performanceof magnetorheological damper

        式中,α、c0、k0、k1、c1、γ、β、n和A為阻尼器模型的參數(shù),如表4所示[34]。

        表4 磁流變阻尼器模型參數(shù)Table 4 Parameters of magnetorheological damper

        當(dāng)物理子結(jié)構(gòu)為磁流變阻尼器時,同樣將其反饋力乘以一個系數(shù),以此來調(diào)節(jié)物理子結(jié)構(gòu)在整體結(jié)構(gòu)中的占比。結(jié)果表明(如圖15所示),本文修正方法可以較好地降低時滯引起的試驗(yàn)誤差。由磁流變阻尼器的模型可知,模型中也存在Bouc-Wen模型,是位移和速度共同決定的。由于本文提出的方法不能同時根據(jù)位移或速度來完全準(zhǔn)確的對反饋力進(jìn)行修正,故修正結(jié)果的誤差隨著物理子結(jié)構(gòu)占比的增加而不斷增加,但總體而言其誤差還是較大幅度地降低了。

        圖15 修正前后誤差對比Fig.15 Error comparison before and after correction

        4 結(jié)論

        實(shí)時混合試驗(yàn)中,在線補(bǔ)償方法存在以下2個主要的問題:第一,在線時滯補(bǔ)償本質(zhì)上是一種預(yù)測方法,使得這種方式無法完全消除時滯,目前的時滯補(bǔ)償方法對于低頻信號能夠較好的補(bǔ)償,但是對高頻信號效果卻不理想,所以實(shí)時混合試驗(yàn)中,時滯普遍存在,并對試驗(yàn)結(jié)果產(chǎn)生影響;第二,實(shí)時混合試驗(yàn)的實(shí)時加載,對加載設(shè)備和控制系統(tǒng)要求過高,對大型構(gòu)件,對加載能力要求較高,現(xiàn)有試驗(yàn)條件很難或者需要付出巨大的經(jīng)濟(jì)成本才能達(dá)到實(shí)時要求。

        針對上述問題,基于數(shù)值子結(jié)構(gòu)積分算法誤差累計規(guī)律提出了一種試驗(yàn)后時滯修正方法,并進(jìn)行線性剛度、阻尼、非線性剛度阻尼驗(yàn)證,結(jié)論如下:

        (1)該方法可同時對試驗(yàn)的位移、速度和加速度結(jié)果進(jìn)行修正;

        (2)試驗(yàn)后時滯修正方法的提出,允許在試驗(yàn)的過程中存在一定的時滯,等于放寬了對于加載設(shè)備和控制系統(tǒng)的要求,一些性能稍顯不足的作動器也可以在保證試驗(yàn)穩(wěn)定的情況下完成實(shí)時混合試驗(yàn),然后通過本文提出的方法對結(jié)果進(jìn)行修正得到較為準(zhǔn)確的結(jié)果。

        (3)對于試驗(yàn)后發(fā)現(xiàn)時滯補(bǔ)償效果并不理想的試驗(yàn),也可以通過本文提出的方法對結(jié)果進(jìn)行修正,從而避免必須重復(fù)試驗(yàn)的窘境,可以節(jié)約試驗(yàn)成本。

        猜你喜歡
        物理模型
        一半模型
        只因是物理
        井岡教育(2022年2期)2022-10-14 03:11:44
        重要模型『一線三等角』
        如何打造高效物理復(fù)習(xí)課——以“壓強(qiáng)”復(fù)習(xí)課為例
        重尾非線性自回歸模型自加權(quán)M-估計的漸近分布
        處處留心皆物理
        我心中的物理
        三腳插頭上的物理知識
        3D打印中的模型分割與打包
        FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
        自拍偷拍亚洲视频一区二区三区| 日本午夜理论片在线观看| 人妻久久久一区二区三区蜜臀| 精品无码av一区二区三区不卡| 日日噜噜夜夜狠狠久久丁香五月| 玩50岁四川熟女大白屁股直播| 粉嫩少妇内射浓精videos| 中国免费一级毛片| 国产一区二区三区特黄| 久久久熟女一区二区三区| 亚洲精品在线国产精品| 欧美性受xxxx白人性爽| 亚洲碰碰人人av熟女天堂| 香蕉成人啪国产精品视频综合网| 在线毛片一区二区不卡视频| 国产一区二区三区在线爱咪咪| 白浆国产精品一区二区| 内地老熟女老少配视频| 精品国产av 无码一区二区三区| 99久久国产亚洲综合精品| 日本免费一区二区在线| 成人国产一区二区三区| 久久99精品久久久久久秒播 | 美女露内裤扒开腿让男人桶无遮挡| 国产乱子伦农村叉叉叉| 一本大道久久精品 东京热| 蜜臀精品一区二区三区| 日本护士xxxxhd少妇| 国产av无码专区亚洲av琪琪| 久久天天躁狠狠躁夜夜2020!| 久久精品国产亚洲AⅤ无码剧情| 亚洲不卡av一区二区三区四区| 色偷偷av一区二区三区| 秋霞午夜无码鲁丝片午夜精品| 久久久久久无中无码| 日产国产精品亚洲高清| 一本色道久久88亚洲精品综合| 欧美z0zo人禽交欧美人禽交| 女人一级特黄大片国产精品| 日本亚洲视频免费在线看| 牛牛在线视频|