999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于混合輸入可控機構(gòu)的連鑄結(jié)晶器非正弦驅(qū)動裝備

2015-12-12 08:38:40劉大偉劉孟朝任廷志
燕山大學(xué)學(xué)報 2015年5期
關(guān)鍵詞:振動

劉大偉,劉孟朝,任廷志,?

(1.燕山大學(xué)機械工程學(xué)院,河北秦皇島066004;2.燕山大學(xué)國家冷軋板帶裝備及工藝工程技術(shù)研究中心,河北秦皇島066004)

基于混合輸入可控機構(gòu)的連鑄結(jié)晶器非正弦驅(qū)動裝備

劉大偉1,劉孟朝2,任廷志2,?

(1.燕山大學(xué)機械工程學(xué)院,河北秦皇島066004;2.燕山大學(xué)國家冷軋板帶裝備及工藝工程技術(shù)研究中心,河北秦皇島066004)

以高性能伺服電機作為動力源已逐漸成為連鑄結(jié)晶器非正弦驅(qū)動裝備的發(fā)展趨勢,針對目前電動式驅(qū)動系統(tǒng)的缺陷,提出兩自由度混合動力驅(qū)動結(jié)晶器非正弦振動的新方案。首先以結(jié)晶器3個基本振動參數(shù)為原始變量,構(gòu)建了與曲柄對稱擺動運行模式相匹配的非正弦振動波形,并闡明了伺服電機和普通電機通過差動輪系運動耦合得到所需運動規(guī)律的機構(gòu)原理;以保證伺服電機跟蹤能力為主要目標(biāo),對NGW型差動輪系的6種運動分配模式進行結(jié)構(gòu)優(yōu)化,通過對比找出了最佳組合模式;最后以實現(xiàn)3組振動規(guī)律為例,對最佳組合模式中的電機運行方式及相應(yīng)的非正弦波形進行仿真計算,以驗證所提出的非正弦振動波形和振動裝備的可行性,并為混合驅(qū)動系統(tǒng)的選型提供理論依據(jù)。

結(jié)晶器;非正弦振動;混合機構(gòu);運動分配

0 引言

結(jié)晶器非正弦振動方式被公認(rèn)為是發(fā)展高效連鑄和提高鑄坯質(zhì)量的關(guān)鍵,因此研發(fā)低成本、運行可靠、維護方便的先進非正弦驅(qū)動裝備一直是國內(nèi)外工程界的熱點問題。

非正弦驅(qū)動裝備根據(jù)動力源可分為液壓式和電動式兩大類。液壓式通過液壓缸直接驅(qū)動結(jié)晶器,可靈活地在線調(diào)整振動參數(shù),但其投資成本高昂,而且液壓缸存在停振、不同步等缺陷[1]。相比較而言,電動式驅(qū)動系統(tǒng)投資成本較低,設(shè)備維護方便,無污染。在電動缸驅(qū)動的結(jié)晶器非正弦振動系統(tǒng)中[2],伺服電機頻繁正反轉(zhuǎn)動,運行模式不合理,影響控制精度,而且其核心傳動構(gòu)件滾柱絲杠的機械性能也難以與液壓缸或曲柄機構(gòu)相媲美;文獻(xiàn)[3?4]將伺服電機的單向變速旋轉(zhuǎn)經(jīng)減速器和曲柄連桿機構(gòu)轉(zhuǎn)化成結(jié)晶器的非正弦振動,該驅(qū)動系統(tǒng)的波形和頻率可在線可調(diào),但卻無法在線調(diào)整振幅;省掉減速機構(gòu),通過直驅(qū)電機帶動曲柄連桿機構(gòu)驅(qū)動結(jié)晶器,在結(jié)構(gòu)上最為簡單,但振幅同樣無法在線調(diào)整[5];文獻(xiàn)[6?8]中采用變頻電機,通過機械式非勻速機構(gòu)和曲柄連桿機構(gòu)合成非正弦振動,這類驅(qū)動系統(tǒng)無需伺服控制,降低了設(shè)備的成本,但運行中除頻率外,振幅和波形均無法在線調(diào)整。

綜上所述,在電動式驅(qū)動裝備中,雖然曲柄連桿機構(gòu)的耐磨損和抗沖擊性能均優(yōu)于滾珠絲杠,但以曲柄連桿機構(gòu)為波形發(fā)生器的驅(qū)動裝備無法實時控制振幅,不符合高效連鑄的發(fā)展要求。因此提出通過混合輸入可控機構(gòu)驅(qū)動曲柄連桿機構(gòu)實現(xiàn)結(jié)晶器非正弦振動,不僅可在線調(diào)整所有振動參數(shù),而且能保證伺服電機以合理模式運行。

1 混合輸入可控驅(qū)動裝備的工作原理

1.1 曲柄擺動運行模式

以常見的弧型連鑄機為例,曲柄連桿機構(gòu)激發(fā)結(jié)晶器振動的結(jié)構(gòu)簡圖如圖1所示,結(jié)晶器振動臺在板簧的約束下按弧形連鑄機設(shè)計半徑R做仿弧運動,其中板簧按照四連桿導(dǎo)向機構(gòu)布置,在初始位置時板簧延長線的交點為仿弧軌跡的圓心,在結(jié)晶器小振幅振動時,可實現(xiàn)精確的仿弧約束。連桿AB與振動臺鉸接,當(dāng)曲柄OA旋轉(zhuǎn)時振動臺產(chǎn)生振動。目前現(xiàn)場使用的結(jié)晶器驅(qū)動裝備中,曲柄均為單向連續(xù)旋轉(zhuǎn),曲柄OA的長度決定了振動臺的振幅。又因為運行過程中曲柄的長度不易改變,故以偏心軸為波形發(fā)生器的電動式驅(qū)動裝備均無法實現(xiàn)振幅在線調(diào)整。

圖1 結(jié)晶器波形發(fā)生機構(gòu)簡圖Fig.1 Oscillation mechanism diagram of mold

若改變偏心軸的運轉(zhuǎn)模式,使其關(guān)于x軸對稱擺動運行,擺動角為β。考慮到振動臺的軌跡半徑和連桿AB長度均遠(yuǎn)遠(yuǎn)大于曲柄OA長度,可認(rèn)為點B沿y軸往復(fù)運動,則振動臺的振幅可表示為

式中,e為曲柄的長度,mm,h為結(jié)晶器振幅,mm。從式(1)中可以看出,在曲柄不變的情況下,改變擺動角β可以調(diào)整振幅h。其中結(jié)晶器的最大振幅為e,因此設(shè)計時只需保證曲柄長度e不小于結(jié)晶器的最大振幅,即可實現(xiàn)結(jié)晶器振幅任意調(diào)整。

為便于對結(jié)晶器振動參數(shù)的控制,以振幅h(mm),頻率f(min-1)和波形偏心率α為原始變量,構(gòu)建符合非正弦振動特征的曲柄速度函數(shù)。設(shè)曲柄一個回轉(zhuǎn)周期中的速度函數(shù)為

式中,ωp為曲柄角速度,rad/s,t為時間參數(shù),s,a、 b、c0、c1、c2、c3和c4為速度函數(shù)的未知系數(shù),在區(qū)間[0,tα)中,曲柄順時針旋轉(zhuǎn),結(jié)晶器向下運動,在區(qū)間[tα,T]中,曲柄逆時針旋轉(zhuǎn),結(jié)晶器向上運動。其中tα與偏斜率α的關(guān)系為

式中,T為結(jié)晶器振動周期,s,T=60/f。

根據(jù)曲柄的運轉(zhuǎn)模式可知,在區(qū)間[0,tα),ωp要滿足t=0和t=tα?xí)r,ωp=0,并且速度函數(shù)在區(qū)間[0,tα)上的積分值為-β,由此可得下述關(guān)于a和b的線性方程組

式中,tα,h和e均為已知,故可解出參數(shù)a和b。為了保證結(jié)晶器振動過程中無沖擊,曲柄的角速度曲線應(yīng)連續(xù)且光滑,并且曲柄正、反轉(zhuǎn)的角度應(yīng)相等,則在時間區(qū)間[tα,T]內(nèi),速度函數(shù)滿足如下約束條件

由于曲柄的角速度函數(shù)均為多項式,故式(5)中積分和求導(dǎo)運算非常方便,將參數(shù)a和b代入式(5),可得關(guān)于c0、c1、c2、c3和c4的線性方程組,解線性方程組(5),可求得曲柄逆時針旋轉(zhuǎn)的角速度。

1.2 驅(qū)動裝備的構(gòu)型

文獻(xiàn)[9]中給出通過直驅(qū)伺服電機正反轉(zhuǎn),帶動曲柄實現(xiàn)可變擺動角和可變角速度的運行規(guī)律,但電機低頻正反轉(zhuǎn)不利于電機壽命,而且影響控制精度,對電機功率造成浪費。為此借鑒多自由度機構(gòu),提出恒速電機和可變速電機通過差動輪系合成出任意擺動規(guī)律的技術(shù)方案,如圖2所示。

變速電機采用伺服電機,恒速電機采用三相異步交流電動機,混合驅(qū)動形式可有效地降低系統(tǒng)對伺服功率的需求,使設(shè)備具有低成本,大功率和柔性輸出的技術(shù)優(yōu)勢。差動輪系為兩自由度機構(gòu),伺服電機和普通電機分別與輪系中兩個定軸旋轉(zhuǎn)的構(gòu)件固接。運行時伺服電機單向變速旋轉(zhuǎn),普通電機單向恒速旋轉(zhuǎn),通過差動輪系速度耦合特性,得到偏心軸所需的擺動規(guī)律。控制伺服電機的角速度可以對結(jié)晶器振動參數(shù)實現(xiàn)在線實時調(diào)整。

圖2 混合驅(qū)動裝置機構(gòu)簡圖Fig.2 Mechanism of hybrid drive system

2 混合驅(qū)動裝備的逆運動學(xué)分析及選型

在設(shè)計圖2中的差動輪系時,考慮盡量減少齒輪的數(shù)量和加工難度,因此選用傳動效率高,結(jié)構(gòu)簡單的NGW行星輪系(由內(nèi)嚙合齒輪副N,外嚙合齒輪副W及共用的行星輪G組成)。

2.1 運動分配模式

NGW行星輪系由太陽輪1、3,行星架H和若干行星輪2組成,為便于電機安裝,引入一對外嚙合齒輪4和5,齒輪5可與輪系中的行星架或太陽輪固接,組成3種構(gòu)型,如圖3所示。

圖3 驅(qū)動系統(tǒng)的3種構(gòu)型Fig.3 Three kinds of configurations of driving system

伺服電機,普通電機和輸出軸分別與太陽輪和行星架相連,能產(chǎn)生6種運動分配模式如表1所示。

不同運動分配模式對電機的運行特性具有決定性的影響,因此通過逆運動學(xué)分析,反求出伺服電機在6中運動模式下的速度與加速度公式,為后續(xù)選奠定基礎(chǔ),如表2所示。其中i13為齒輪3和齒輪1的齒數(shù)比,i45為齒輪5和齒輪4的齒數(shù)比,ωc為普通電機的額定轉(zhuǎn)速,r/min,αp為輸出軸的角加速度,rad/s2。

表1 運動分配模式Tab.1 Distribution model of movement

表2 伺服電機運動學(xué)參數(shù)Tab.2 Kinematic parameters of servo motor

2.2 基于伺服電機跟蹤能力的選型原則

對于混合輸入機構(gòu),通常期望常速電機提供主要動力,伺服電機提供輔助動力,因此功率分配的經(jīng)濟性是混合輸入機構(gòu)的首要考慮因素,如混合驅(qū)動壓力機設(shè)計[10]。但本設(shè)備與其它混合驅(qū)動設(shè)備的重要區(qū)別在于對合成角速度的精度要求高,因為非正弦波形的精確程度對鑄坯質(zhì)量具有重要的影響,故在本設(shè)備中伺服電機的速度跟蹤能力是決定設(shè)備性能的最主要因素。

為保證伺服電機的跟蹤能力,應(yīng)盡量減小伺服電機的速度波動范圍及加速度[11]。由表2可知,在ωc和ωp一定的情況下,差動輪系的結(jié)構(gòu)決定了伺服電機的運動參數(shù),故對差動輪系提出如下設(shè)計原則:

1)保證伺服電機單向旋轉(zhuǎn),速度波動范圍盡量小,并且接近額定轉(zhuǎn)速,以避免出現(xiàn)低頻共振和功率浪費;

2)使伺服電機加速度的絕對值盡量小,以保證電機跟蹤能力;

3)令伺服電機額定轉(zhuǎn)速不超過3 000 r/min,普通電機選取同步轉(zhuǎn)速為3 000 r/min和1 500r/min的兩種類型,以保證設(shè)備的通用性和經(jīng)濟性。

2.3 混合驅(qū)動裝備的最佳構(gòu)型

對于表1中的每種組合模式,在理論上都有滿足上述3個原則的最優(yōu)結(jié)構(gòu),為此對6種組合模式分析對比,以實現(xiàn)伺服電機加速度的絕對值最小為標(biāo)準(zhǔn)選取最易保證高精度運行的組合模式。

當(dāng)曲柄角加速度最大時,反求出的伺服電機角加速度也最大,此時結(jié)晶器處于極限運行狀態(tài),即振幅、頻率和偏斜率均為最大值。根據(jù)本文中的波形特點和現(xiàn)場經(jīng)驗,取hmax=4.5 mm,fmax=200 min-1,αmax=0.25為振動參數(shù)的最大值,曲柄長度e=6 mm>hmax,將這些參數(shù)代入式(3)~(6)可求出差動輪系輸出軸轉(zhuǎn)速的最大值ωpmax=112.5 r/min和最小值ωpmin=-216 r/min。

另外,由于普通電機功率不同時,其額定轉(zhuǎn)速略有差異,為統(tǒng)一計算,以電機同步轉(zhuǎn)速代替額定轉(zhuǎn)速,即ωc取1 500 r/min和3 000 r/min兩個值。

i13和i45是決定差動輪系的結(jié)構(gòu)參數(shù),一般設(shè)計中單級圓柱齒輪的傳動比不大于6,另外考慮差動輪系的尺寸,設(shè)定i13∈[1.2,7],i45∈[1,6]。為得到表1中伺服電機角加速度絕對值最小的組合模式,采用流程圖3中的計算方法,分別求取6種組合模式的最佳的結(jié)構(gòu)參數(shù)i13和i45。

圖3 差動輪系結(jié)構(gòu)參數(shù)計算流程圖Fig.3 Calculation flow of parameters of differential gear train

在原始參數(shù)中,普通電機的轉(zhuǎn)速可選取兩種,應(yīng)分別進行計算,以ωc=3 000 r/min時為例,對表1中第Ⅰ種組合模式的最優(yōu)結(jié)構(gòu)參數(shù)進行計算。將i13和i45作為兩個變量,代入表2中第Ⅰ種組合模式的伺服電機轉(zhuǎn)速公式,求出差動輪系不同結(jié)構(gòu)下伺服電機的最大角速度ωsmax和最小角速度ωsmin,如圖4(a)和(b)所示,其中負(fù)值表示與普通電機轉(zhuǎn)速相反。

圖4 不同機構(gòu)參數(shù)下伺服電機角速度極值Fig.4 Extreme values of servo motor angular acceleration with different mechanism parameters

為保證伺服電機的額定轉(zhuǎn)速不超過3 000 r/min,根據(jù) 3個約束條件:|ωsmax|≤3 000,|ωsmin|≤3 000和ωsmaxωsmin>0,選取相應(yīng)的i13和i45組合數(shù)據(jù)。

然后將其代入表2中第Ⅰ種組合模式的伺服電機角加速度公式,由于αp的值是在結(jié)晶器振動參數(shù)選取后就固定,因此只需比較包含i13和i45的多項式值,將該多項式定義為伺服電機角加速度系數(shù)p。對|p|進行比較,最小|p|對應(yīng)的i13、i45為最佳組合。

當(dāng)|p|為最小值時,往往對應(yīng)多組i13、i45,此時需再對其進一步優(yōu)化,選取|ωsmax|和|ωsmin|中最接近伺服電機額定轉(zhuǎn)速時對應(yīng)的i13和i45,以提高伺服電機功率利用率。經(jīng)計算,6種運動分配模式的最優(yōu)結(jié)構(gòu)參數(shù)如表3所示。對比表3中的數(shù)據(jù)可知,第Ⅵ種組合模式下伺服電機的角加速度系數(shù)p和速度波動范圍最小,而第Ⅲ模式下最大,另外第Ⅴ種與第Ⅵ種差動輪系結(jié)構(gòu)相同,伺服電機運動參數(shù)也相差不多,但是當(dāng)常速電機的同步轉(zhuǎn)速ωc=3 000 r/min,相應(yīng)伺服電機的額定轉(zhuǎn)速超過3 000 r/min,因此考慮伺服電機跟蹤能力以及設(shè)備的經(jīng)濟性,選擇ωc=3 000 r/min時的第Ⅵ種為最佳組合模式。

表3 6種運動分配模式下的傳動比優(yōu)化結(jié)果Tab.3 Optimization results of transmission ratio of six motion allocation modes

3 仿真計算

根據(jù)優(yōu)化后的結(jié)構(gòu)參數(shù)對第Ⅵ種組合方案進行齒數(shù)選配,令各齒輪的齒數(shù)分別為Z1=18,Z2=54,Z3=126,Z4=30,Z5=30,行星輪個數(shù)k=3,經(jīng)驗證行星輪滿足均勻安裝和鄰接條件。另外選取普通三相異步交流電動機額定轉(zhuǎn)速為 2 935 r/min,伺服電機額定轉(zhuǎn)速3 000 r/min。以小方坯連鑄機為例,其中小方坯的斷面尺寸為250 mm× 150 mm,則結(jié)晶器隨拉速不同時,分別取表4中的振動參數(shù),結(jié)晶器振動波形和伺服電機速度曲線如圖5所示,其中代號a,b,c的3條曲線分別對應(yīng)表4中序號為a,b,c的數(shù)據(jù)。

表4 結(jié)晶器基本振動參數(shù)Tab.4 Basic oscillation parameters of mould

圖5 不同振動參數(shù)下的非正弦波形及伺服電機轉(zhuǎn)速Fig.5 Non?sinusoidal wave and servo motor angular velocity under different oscillation parameters

從圖5(a)中可以看出,調(diào)節(jié)偏斜率α可以改變波形,隨著α的減小,波形中正速度和負(fù)速度的時間比值相應(yīng)減小,完全符合非正弦波形的特征。當(dāng)α=0時,正、負(fù)速度的時間相等,結(jié)晶器產(chǎn)生一種類似正弦波形的規(guī)律,所以本文中基于偏心軸擺動運行模式構(gòu)建的結(jié)晶器波形,在生產(chǎn)中可以實現(xiàn)正弦規(guī)律與非正弦規(guī)律的切換,更加便于振動規(guī)律的控制。

從圖5(b)中可以看出,在普通電機恒速輸出的情況下,只需調(diào)節(jié)伺服電機的角速度,就能夠?qū)崿F(xiàn)結(jié)晶器所有振動參數(shù)的在線控制。通過合理設(shè)計差動輪系的結(jié)構(gòu)和運動分配模式,不僅能保證伺服電機單向運行,而且可以使伺服電機的速度及加速度在最小范圍內(nèi)波動,有利于提高伺服電機的跟蹤能力,確保驅(qū)動系統(tǒng)的精度。

4 結(jié)論

1)采用曲柄對稱擺動激發(fā)結(jié)晶器非正弦振動的方式,可實現(xiàn)結(jié)晶器全振動參數(shù)在線實時控制。運行中伺服電機只需單向變速旋轉(zhuǎn)即可實現(xiàn)可控擺動輸出,合理的運行模式對提高伺服電機的壽命和驅(qū)動系統(tǒng)控制精度極為有利。

2)混合驅(qū)動裝備與電液伺服驅(qū)動相比,成本低、無泄漏、維護方便;與電動缸驅(qū)動相比,避免電機頻繁正反轉(zhuǎn),具有節(jié)能降耗的特點,并且其機械系統(tǒng)的抗沖擊、耐磨損性能突出,因此具有巨大的應(yīng)用潛力。

3)采用伺服電機驅(qū)動行星架,普通常速電機驅(qū)動大太陽輪,小太陽輪作為輸出的運動分配方案時,伺服電機的角加速度最大值及其角速度波動范圍最小,其跟蹤能力最好。

[1]方濤,向忠輝,李宏磊.連鑄機結(jié)晶器液壓振動常見故障分析與對策[J].機床與液壓,2014,42(20):147?153.

[2]訾福寧.板坯連鑄結(jié)晶器電動缸非正弦振動技術(shù)與應(yīng)用[J].連鑄,2013(1):22?26.

[3]Yao Yunfeng,Li Xiankui,F(xiàn)ang Yiming,et al.Study of non?sinusoidal oscillation of continuous casting mold driven by servomo?tor[J].Journal of Iron and Steel Research,International,2008,15(1):558?562.

[4]方一鳴,李宮胤,李建雄,等.伺服電機驅(qū)動連鑄結(jié)晶器振動系統(tǒng)建模與分析[J].儀器儀表學(xué)報,2014,35(11):2616?2623.

[5]曹鐵軍,王鵬,劉奎奎.直接驅(qū)動結(jié)晶器振動發(fā)生裝置及其振動方法:中國,200910057435.9[P].2009?06?17.

[6]張立平,汪建曉,盧清華,等.新型結(jié)晶器非正弦振動系統(tǒng)彈性振動響應(yīng)的研究[J].機械工程學(xué)報,2011,47(9):91?96.

[7]劉大偉,任廷志.非圓齒輪驅(qū)動結(jié)晶器非正弦振動的研究[J].中國機械工程,2013,24(3):327?330.

[8]Zhang Xingzhong,Zheng Xueran,Liu Qingguo,et al.Investigation and application of non?sinusoidal oscillation technique of mold[J]. Journal of Iron and Steel,International,2013,20(2):19?24.

[9]杜傳得,武青,徐敦石.一種振幅在線可調(diào)的結(jié)晶器振動發(fā)生裝置及其振動方法:中國,201210291325.0[P].2012?08?16.

[10]He Yupeng,Zhao Dunsheng,Zou Jun,et al.Study of utilizing dif?ferential gear train to achieve hybrid mechanism of mechanical press[J].Science in China,Series E:Technological Sciences,2007,50(1):69?80.

[11]孟彩芳,張策,陸永輝,等.新型混合輸入式曲柄壓力機械的優(yōu)化設(shè)計[J].中國機械工程,2003,14(10):834?835.

Non?sinusoidal driving equipment of continuous casting mold based on hybrid input mechanism

LIU Da?wei1,LIU Meng?zhao2,REN Ting?zhi2
(1.School of Mechanical Engineering,Yanshan University,Qinhuangdao,Hebei 066004,China;2.National Engineering Technology Research Center of Cold Rolling Strip
Equipment and Technology,Yanshan University,Qinhuangdao,Hebei 066004,China)

The application of the high performance servo motor in non?sinusoidal driving equipment becomes trend of the develop?ment of continuous casting.A new project containing two?DOF hybrid input machine is presented in order to overcome the defects of existing electric system.Firstly,a non?sinusoidal oscillation wave that matches symmetrical swing style of the crack is constructed with the three basic oscillation parameters of the mold.A differential gear mechanism driven by a servomotor and ordinary motor to implement the given angular velocity is illustrated.The structure optimization of six motion distribution modes of NGW gear train is carried out to establish the optimum mode with the goal to guarantee the traceability of the servomotor.Finally,the operation way of the servomotor and the corresponding non?sinusoidal waves under the optimum transmission pattern are simulated to fulfill three groups of oscillation parameters,which can validate the feasibility of the presented non?sinusoidal oscillation wave and equipment and provide theoretical basis for its model selection.

mold;non?sinusoidal vibration;hybrid mechanism;motion distribution

TH132

A

10.3969/j.issn.1007?791X.2015.05.005

1007?791X(2015)05?0414?06

2015?08?19 基金項目:河北省自然科學(xué)基金鋼鐵聯(lián)合研究基金資助項目(E2011203078);燕山大學(xué)博士基金(B870)

劉大偉(1984?),男,河北滄州人,博士,講師,主要研究方向為高效連鑄裝備及工藝;?通信作者:任廷志(1960?),男,黑龍江大慶人,博士,教授,博士生導(dǎo)師,主要研究方向為機械設(shè)計及智能化,Email:rtz@ysu.edu.cn。

猜你喜歡
振動
振動的思考
某調(diào)相機振動異常診斷分析與處理
振動與頻率
This “Singing Highway”plays music
具非線性中立項的廣義Emden-Fowler微分方程的振動性
中立型Emden-Fowler微分方程的振動性
基于ANSYS的高速艇艉軸架軸系振動響應(yīng)分析
船海工程(2015年4期)2016-01-05 15:53:26
主回路泵致聲振動分析
UF6振動激發(fā)態(tài)分子的振動-振動馳豫
計算物理(2014年2期)2014-03-11 17:01:44
帶有強迫項的高階差分方程解的振動性
主站蜘蛛池模板: 国产福利一区在线| 亚洲av日韩综合一区尤物| 无码一区18禁| 国产午夜福利在线小视频| 国产丝袜无码精品| 国模沟沟一区二区三区| 77777亚洲午夜久久多人| 国产丝袜91| 九九久久99精品| 亚洲中文字幕久久精品无码一区| 55夜色66夜色国产精品视频| 日日摸夜夜爽无码| 久久综合丝袜长腿丝袜| 老司机午夜精品网站在线观看| 欧美成人影院亚洲综合图| 欧美福利在线观看| 国产精品永久久久久| 伊伊人成亚洲综合人网7777| 中文字幕无码av专区久久| 亚洲国产在一区二区三区| 色窝窝免费一区二区三区| 亚洲色欲色欲www网| 亚洲欧美在线综合图区| 国产女人爽到高潮的免费视频| 国产精品网曝门免费视频| 亚洲精品无码久久久久苍井空| 亚洲第一国产综合| 国产区人妖精品人妖精品视频| 成年午夜精品久久精品| 国产黄在线观看| 国产清纯在线一区二区WWW| 国产一区二区三区精品欧美日韩| 自拍偷拍欧美| 一区二区三区精品视频在线观看| 九九视频在线免费观看| 国产免费人成视频网| 国产精品一区二区在线播放| 国产精品久久久久无码网站| 亚洲首页国产精品丝袜| 九九久久99精品| 国产精品嫩草影院av| 久久毛片网| 国产在线视频欧美亚综合| 亚洲无码91视频| 久久久久人妻一区精品色奶水| 国产一级特黄aa级特黄裸毛片| 在线五月婷婷| 亚洲第一页在线观看| 韩国福利一区| 在线观看国产精美视频| 日韩专区欧美| 制服丝袜在线视频香蕉| 国产精品久线在线观看| 五月六月伊人狠狠丁香网| 一本无码在线观看| 欧美综合一区二区三区| 大陆精大陆国产国语精品1024| 久久国产黑丝袜视频| 久久人人妻人人爽人人卡片av| 全午夜免费一级毛片| 在线中文字幕网| 国产成人久久综合777777麻豆| 久久频这里精品99香蕉久网址| 中文字幕无码中文字幕有码在线| 国产激情影院| 国内精品九九久久久精品| 成人看片欧美一区二区| 欧美激情第一区| 亚洲伊人久久精品影院| 欧美午夜网| 国产成人免费手机在线观看视频 | 欧美一级在线| 97精品久久久大香线焦| 亚洲精品人成网线在线| 国产永久免费视频m3u8| 国产精品久久久久久久伊一| 人妻少妇乱子伦精品无码专区毛片| 亚洲黄色成人| 午夜视频www| 狠狠色成人综合首页| 亚洲成人播放| 亚洲精品视频免费|