鞏巖博,鄭大勇,王維彬
(北京航天動(dòng)力研究所,北京 100076)
大推力液體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)復(fù)雜,涉及到流體、燃燒、傳熱、力學(xué)、機(jī)械、控制等專業(yè),是一個(gè)復(fù)雜的非線性系統(tǒng),其性能受到自身部件特性和外部邊界條件的影響。在研制過(guò)程中,掌握對(duì)發(fā)動(dòng)機(jī)性能影響較大的參數(shù)并在組件研制階段對(duì)其加以關(guān)注和控制,對(duì)于保證發(fā)動(dòng)機(jī)穩(wěn)定可靠地工作是非常必要的[1]。
工程中常使用敏感性來(lái)評(píng)價(jià)影響因素與關(guān)鍵參數(shù)的相關(guān)性和對(duì)關(guān)鍵參數(shù)的影響程度[2]。將發(fā)動(dòng)機(jī)性能參數(shù)隨內(nèi)外影響因素的變化而產(chǎn)生的相對(duì)變化量叫做發(fā)動(dòng)機(jī)性能對(duì)影響因素的敏感性。相對(duì)變化量越大,敏感性就越高,表明該影響因素對(duì)發(fā)動(dòng)機(jī)性能參數(shù)的影響越大[3]??梢?jiàn),通過(guò)敏感性分析可以定量地確定各影響因素對(duì)發(fā)動(dòng)機(jī)性能參數(shù)影響程度的量級(jí)和大小,篩選出對(duì)性能影響較大的關(guān)鍵因素,忽略對(duì)性能影響微弱的次要因素,讓設(shè)計(jì)人員在發(fā)動(dòng)機(jī)研制過(guò)程中能夠做到主次分明,有的放矢,一方面可以提高發(fā)動(dòng)機(jī)的設(shè)計(jì)水平,另一方面還可以減輕設(shè)計(jì)人員的工作負(fù)擔(dān)。
設(shè)計(jì)人員對(duì)發(fā)動(dòng)機(jī)組合件的參數(shù)敏感性[4-6]和整機(jī)的性能敏感性[1,3,7-8]均進(jìn)行過(guò)相關(guān)研究,但到目前為止,在發(fā)動(dòng)機(jī)整機(jī)性能的多因素敏感性分析時(shí)使用的均是極差分析法,不能很好地反映影響因素對(duì)性能參數(shù)的影響顯著性[9]。為了提高發(fā)動(dòng)機(jī)性能敏感性分析的準(zhǔn)確度,本文基于發(fā)動(dòng)機(jī)非線性靜態(tài)特性模型,使用方差分析法對(duì)某型液氧/甲烷發(fā)動(dòng)機(jī)進(jìn)行性能敏感性分析,得到各內(nèi)外影響因素對(duì)發(fā)動(dòng)機(jī)性能影響的顯著性大小,從中篩選出對(duì)發(fā)動(dòng)機(jī)性能影響較大的關(guān)鍵因素,并對(duì)關(guān)鍵因素進(jìn)行分析,為該型發(fā)動(dòng)機(jī)的設(shè)計(jì)和試驗(yàn)提供理論依據(jù)。
敏感性分析可以確定發(fā)動(dòng)機(jī)內(nèi)部組件特性和外部邊界條件對(duì)發(fā)動(dòng)機(jī)性能影響的量級(jí)與大小。本文通過(guò)發(fā)動(dòng)機(jī)靜態(tài)特性的非線性模型建立目標(biāo)參數(shù)與各影響因素的關(guān)系。由于影響因素眾多,因此可以對(duì)單個(gè)因素作用時(shí)的影響進(jìn)行分析,也可以對(duì)多個(gè)因素同時(shí)作用時(shí)的影響進(jìn)行分析。前者叫做單因素敏感性分析,后者叫做多因素敏感性分析[10]。
根據(jù)發(fā)動(dòng)機(jī)部件的特點(diǎn)建立各部件的靜態(tài)特性方程,然后根據(jù)發(fā)動(dòng)機(jī)的系統(tǒng)組成對(duì)發(fā)動(dòng)機(jī)各流路建立壓力平衡、流量平衡和功率平衡的關(guān)系式,最終建立如下的發(fā)動(dòng)機(jī)系統(tǒng)非線性模型
Fi(D,X)=0,i=1,2,…,n
(1)
式中:D(d1,d2,…,dn)為關(guān)注的目標(biāo)參數(shù);X(x1,x2,…,xn)為發(fā)動(dòng)機(jī)推力和混合比,為包含隨機(jī)信息的發(fā)動(dòng)機(jī)各部件參數(shù)。
(2)
(3)
(4)
其中
當(dāng)系數(shù)行列式不等于零時(shí),上述方程組有唯一解。于是由上式得到的解為
(5)
當(dāng)?shù)玫讲钪敌∮谠O(shè)定值時(shí),計(jì)算收斂。
單因素敏感性分析即研究發(fā)動(dòng)機(jī)性能與某一個(gè)影響因素的相關(guān)性。
若發(fā)動(dòng)機(jī)某項(xiàng)性能參數(shù)Fi會(huì)受到m個(gè)因素的影響,分別記為a1,a2,a3,…,am,則可以根據(jù)發(fā)動(dòng)機(jī)靜態(tài)特性方程建立數(shù)學(xué)模型
Fi=f(a1,a2,a3,…,am) (i=1,2,…,n)
(6)
(7)
可以看出,敏感性越高,發(fā)動(dòng)機(jī)性能參數(shù)隨影響因素變化的變化幅度就越大,即影響因素對(duì)發(fā)動(dòng)機(jī)性能的影響越明顯。單因素敏感性分析方法簡(jiǎn)單,能夠快速直觀地反映發(fā)動(dòng)機(jī)性能隨某一影響因素變化的變化梯度。但是單因素分析方法每次只能對(duì)影響因素空間的某一點(diǎn)進(jìn)行分析,通過(guò)對(duì)某一方向的局部梯度判斷影響因素的作用效果,存在一定的片面性[4]。
單因素敏感性分析操作簡(jiǎn)單,結(jié)果直觀,但只能就單一因素變化時(shí)發(fā)動(dòng)機(jī)性能受到的影響進(jìn)行研究,無(wú)法準(zhǔn)確分析多因素共同作用對(duì)發(fā)動(dòng)機(jī)性能產(chǎn)生的影響。液體火箭發(fā)動(dòng)機(jī)結(jié)構(gòu)復(fù)雜,部件的非線性和參數(shù)的耦合性都很強(qiáng),在工作時(shí)往往會(huì)有很多參數(shù)偏離設(shè)計(jì)值,使發(fā)動(dòng)機(jī)性能受到多個(gè)參數(shù)變化的共同作用。在這種情況下,單因素敏感性分析就不能滿足需求了,需要進(jìn)行多因素敏感性分析。
正交試驗(yàn)法即從所有的排列組合中選取部分有特點(diǎn)、有代表性的子樣,使各影響因素、各因素水平的分布具有均勻散布和整齊可比的特點(diǎn),即各個(gè)因素的每個(gè)不同水平在試驗(yàn)中出現(xiàn)的次數(shù)均相同;任意兩個(gè)因素所有不同水平的搭配均在試驗(yàn)中出現(xiàn),且出現(xiàn)的次數(shù)都相同。這樣就可以通過(guò)部分試驗(yàn)得到與全部試驗(yàn)同等可信的結(jié)果,極大減少試驗(yàn)的次數(shù)。一般使用正交表來(lái)安排正交試驗(yàn)[12]。
在以前的研究中,常采用極差分析法對(duì)試驗(yàn)結(jié)果進(jìn)行統(tǒng)計(jì)和研究。
假設(shè)多因素敏感性分析的試驗(yàn)結(jié)果為yI(i=1,2,…,n),將第j個(gè)因素同水平的結(jié)果加起來(lái),記作Xk=(k=1,2,…,l),表示第j個(gè)因素中水平為k的結(jié)果之和。若每種因素水平為k的試驗(yàn)有z次,則將該因素水平為k的平均結(jié)果記作xk=Xk/z,計(jì)算該因素所有水平的平均結(jié)果的極差Rj=max(xk)-min(xk),就可以反映該因素對(duì)目標(biāo)參數(shù)的影響程度。極差越大,說(shuō)明該因素水平變化對(duì)目標(biāo)參數(shù)的影響越大。
極差法原理簡(jiǎn)單,計(jì)算速度快,結(jié)果直觀易懂,但是只能分析出不同因素對(duì)目標(biāo)參數(shù)影響程度的相對(duì)大小,無(wú)法判斷每個(gè)因素對(duì)目標(biāo)參數(shù)影響程度的顯著性,而且也不能區(qū)分影響因素水平變化引起的目標(biāo)參數(shù)變化和試驗(yàn)誤差引起的目標(biāo)參數(shù)變化。要排除試驗(yàn)誤差帶來(lái)的影響,并分析出不同因素對(duì)目標(biāo)參數(shù)影響程度的顯著性,就需要使用方差分析法[13-17]。
若使用正交表Ln(qm)進(jìn)行正交試驗(yàn),將第k號(hào)試驗(yàn)結(jié)果記為yk,k=1,2,…,n。并記
(8)
(9)
(10)
(11)
將二次型ST的秩稱為ST的自由度,記作fT;同理定義Sj的自由度,記為fj。則有
fT=n-1
(12)
fj=q-1,j=1,2,…,m
(13)
(14)
若正交表中有未安排影響因素的空列,則該列的離差平方和作為誤差平方和,記作Se,它的自由度記作fe。
如果試驗(yàn)的結(jié)果相互獨(dú)立,且服從等方差為σ2的正態(tài)分布,那么S1,S2,Sm互相獨(dú)立。如果在某一列j中,T1j,T2j,…,Tqj相互獨(dú)立,且同樣服從方差為σ2的正態(tài)分布,則認(rèn)為該列因素對(duì)試驗(yàn)結(jié)果的影響不顯著,則有
Sj/σ2~x2(fj),j=1,2,…,m
(15)
選取統(tǒng)計(jì)量
(16)
如果影響因素的作用不顯著,則有
Fj~F(fj,fe),j=1,2,…,m
(17)
給定顯著性水平α,可以通過(guò)查F分布表可以得到Fj(fj,fe)的值,然后由公式F1-α(fj,fe)=1/Fα(fj,fe)求出F1-α(fj,fe)的值,如果Fj>F1-α(fj,fe),則認(rèn)為該因素的影響顯著,反之則認(rèn)為影響不顯著。當(dāng)判斷因素的影響顯著時(shí),顯著性水平α越小,則認(rèn)為影響的顯著性越高。
在進(jìn)行方差分析時(shí),可以首先求出每一列的平均離差平方和,即
(18)
(19)
可以看出,方差分析能夠有效地把影響因素的離散程度和與誤差的離散程度區(qū)分開(kāi),并通過(guò)F檢驗(yàn)法分析出每個(gè)因素對(duì)試驗(yàn)結(jié)果影響的顯著性大小,進(jìn)而找出對(duì)試驗(yàn)結(jié)果造成影響的主要因素。與極差法相比,方差分析法更準(zhǔn)確、更可信,分析得出的信息量也更大。
某型在研65 t級(jí)液氧/甲烷發(fā)動(dòng)機(jī)采用燃?xì)獍l(fā)生器循環(huán),同軸渦輪泵布局,發(fā)動(dòng)機(jī)系統(tǒng)如圖1 所示。在地面級(jí)發(fā)動(dòng)機(jī)中,推力Fe和混合比re是設(shè)計(jì)人員最關(guān)心的性能參數(shù),因此選取推力和混合比作為研究的目標(biāo)參數(shù)。
圖1 液氧/甲烷發(fā)動(dòng)機(jī)系統(tǒng)示意圖Fig.1 Schematic diagram of liquid oxygen/ methane engine system
能夠?qū)Πl(fā)動(dòng)機(jī)性能產(chǎn)生影響的干擾因素分為內(nèi)部和外部干擾因素。外部干擾因素是由發(fā)動(dòng)機(jī)工質(zhì)和外部環(huán)境帶來(lái)的干擾,可以精確測(cè)量并采取措施加以控制,一般認(rèn)為是非隨機(jī)變量;內(nèi)部干擾因素則是由發(fā)動(dòng)機(jī)本身帶來(lái)的干擾,即在制造過(guò)程中,因加工誤差,零件安裝差異以及部件液流試驗(yàn)測(cè)量誤差等造成發(fā)動(dòng)機(jī)部件特性具有一定的隨機(jī)散布性,一般認(rèn)為是隨機(jī)變量[2]。
在這些隨機(jī)變量中,有兩類對(duì)發(fā)動(dòng)機(jī)工況的影響最直接,即流路的流阻特性和部件的性能參數(shù)。這兩類參數(shù)也恰好最容易受到加工誤差和工藝穩(wěn)定性的影響,是需要重點(diǎn)關(guān)注的。此外,邊界條件也會(huì)影響發(fā)動(dòng)機(jī)的工作狀態(tài),是需要關(guān)注的。因此,我們從外部邊界條件、流路的流阻特性和渦輪泵的性能特性中選出19個(gè)比較有代表性的參數(shù)[18-19],分別為甲烷泵入口壓力(pipf)、氧泵入口壓力(pipo)、甲烷主管路流阻系數(shù)(Cqpmf)、氧主管路流阻系數(shù)(Cqpmo)、甲烷副管路流阻系數(shù)(Cqpsf)、氧副管路流阻系數(shù)(Cqpso)、推力室甲烷噴嘴壓降(Cqcf)、推力室氧噴嘴壓降(Cqco)、發(fā)生器甲烷噴嘴壓降(Cqgf)、發(fā)生器氧噴嘴壓降(Cqgo)、渦輪噴嘴流量系數(shù)(μt)、渦輪效率系數(shù)(βt)、甲烷泵效率系數(shù)(ηpf)、氧泵效率系數(shù)(ηpo)、甲烷主閥流阻系數(shù)(Cqmfv)、氧主閥流阻系數(shù)(Cqmov)、甲烷副閥流阻系數(shù)(Cqgfv)、氧副閥流阻系數(shù)(Cqgov)、燃?xì)夤苈妨髯柘禂?shù)(Cqpg)。
選取影響因素后,需要確定影響因素的水平數(shù)。因素的水平數(shù)是指該因素取值的個(gè)數(shù)。一般情況下,選擇兩水平數(shù)進(jìn)行分析,得到對(duì)目標(biāo)參數(shù)的影響一般是線性的,選擇三水平數(shù)或者更多的水平數(shù),得到的結(jié)果一般呈二次曲線,甚至是更高次曲線,可以更好地趨近真實(shí)的情況。因此,選擇三水平數(shù)或更多的水平數(shù),能夠提高分析的準(zhǔn)確性。本文從仿真準(zhǔn)確性和計(jì)算規(guī)模兩個(gè)方面綜合考慮,在單因素敏感性分析中,每個(gè)因素均選擇五水平數(shù),即基準(zhǔn)值和分別偏離基準(zhǔn)值-20%,-10%,10%,20%;在多因素敏感性分析中,每個(gè)因素均選擇三水平數(shù),即基準(zhǔn)值和分別偏離基準(zhǔn)值-10%,10%。
使用2.1節(jié)中介紹的原理,對(duì)選出的19個(gè)參數(shù)進(jìn)行單因素敏感性分析,得到發(fā)動(dòng)機(jī)推力敏感度如表1所示;發(fā)動(dòng)機(jī)混合比敏感度如表2所示;推力和混合比對(duì)不同影響因素的平均敏感度直方圖如圖2所示,其中,橫坐標(biāo)1~19代表的影響因素分別為甲烷泵入口壓力(pipf)、氧泵入口壓力(pipo)、甲烷主管路流阻系數(shù)(Cqpmf)、氧主管路流阻系數(shù)(Cqpmo)、甲烷副管路流阻系數(shù)(Cqpsf)、氧副管路流阻系數(shù)(Cqpso)、推力室甲烷噴嘴壓降(Cqcf)、推力室氧噴嘴壓降(Cqco)、發(fā)生器甲烷噴嘴壓降(Cqgf)、發(fā)生器氧噴嘴壓降(Cqgo)、渦輪噴嘴流量系數(shù)(μt)、渦輪效率系數(shù)(βt)、甲烷泵效率系數(shù)(ηpf)、氧泵效率系數(shù)(ηpo)、甲烷主閥流阻系數(shù)(Cqmfv)、氧主閥流阻系數(shù)(Cqmov)、甲烷副閥流阻系數(shù)(Cqgfv)、氧副閥流阻系數(shù)(Cqgov)、燃?xì)夤苈妨髯柘禂?shù)(Cqpg)。
圖2 平均敏感度直方圖Fig.2 Average sensitivity histogram
由表1和表2可以看出,在選取的水平下,發(fā)動(dòng)機(jī)性能參數(shù)與各影響因素的變化呈正相關(guān),即二者具有相同的變化趨勢(shì)。其中,有部分影響因素的變化與性能參數(shù)的變化呈線性關(guān)系。對(duì)于發(fā)動(dòng)機(jī)推力,影響最大的10個(gè)因素對(duì)其影響程度依次為:βt>ηpf>ηpo>Cqgov>Cqgfv>Cqgo>μt>Cqgf>Cqmov>Cqco,其中渦輪效率、甲烷泵效率和氧泵效率對(duì)發(fā)動(dòng)機(jī)推力的影響最大,未參與排序的9個(gè)因素對(duì)其影響非常小;對(duì)于發(fā)動(dòng)機(jī)混合比,影響最大的10個(gè)因素對(duì)其影響程度依次為:βt>ηpf>ηpo>Cqmov>Cqco>Cqcf>Cqgfv>μt>Cqgf>Cqgov,其中對(duì)發(fā)動(dòng)機(jī)混合比影響最大的同樣是渦輪效率、甲烷泵效率和氧泵效率,未參與排序的9個(gè)因素對(duì)混合比幾乎沒(méi)有影響。
表1 發(fā)動(dòng)機(jī)推力敏感度Tab.1 Thrust sensitivity of engine
表2 發(fā)動(dòng)機(jī)混合比敏感度Tab.2 Mixing ratio sensitivity of engines
雖然對(duì)發(fā)動(dòng)機(jī)推力和混合比影響最大的3個(gè)因素是相同的,但是有一些因素對(duì)推力和混合比的影響差別較大。例如副系統(tǒng)氧閥流阻系數(shù)Cqgov,對(duì)推力的影響程度排在第4位,敏感度為0.132 3,對(duì)混合比的影響則排在第10位,敏感度僅為0.058 1,二者相差2.3倍;又如推力室氧噴嘴流阻系數(shù)Cqco,對(duì)推力的影響程度排在第10位,敏感度僅為0.039 8,但對(duì)混合比的影響程度卻排在第5位,敏感度為0.164 0,二者相差4.1倍。
根據(jù)單因素敏感性分析,發(fā)現(xiàn)部分因素對(duì)發(fā)動(dòng)機(jī)推力和混合比的影響非常小,為減小不必要的工作量,提高計(jì)算效率,可以在進(jìn)行多因素敏感性分析的時(shí)候?qū)⑵浜雎?。通過(guò)分析,從19個(gè)因素中分別選出對(duì)推力和混合比影響較為明顯的因素進(jìn)行多因素敏感性分析。
由于發(fā)動(dòng)機(jī)靜態(tài)特性方程組的非線性很強(qiáng)[20],因此當(dāng)多因素同時(shí)變化時(shí),求解方程組的難度會(huì)增大,甚至可能出現(xiàn)計(jì)算不收斂的情況。為了保證方程收斂并提高計(jì)算效率,僅選取3個(gè)水平進(jìn)行仿真。根據(jù)選擇的因素?cái)?shù)和水平數(shù),選擇L27(313)正交表設(shè)計(jì)仿真試驗(yàn)。發(fā)動(dòng)機(jī)推力敏感性試驗(yàn)的組合方式及試驗(yàn)結(jié)果如表3所示,發(fā)動(dòng)機(jī)混合比敏感性試驗(yàn)的組合方式與試驗(yàn)結(jié)果如表4所示。
表3 推力敏感度試驗(yàn)及結(jié)果Tab.3 Thrust sensitivity test and results
表4 混合比敏感度試驗(yàn)及結(jié)果Tab.4 Mixed ratio sensitivity tests and results
使用方差分析法對(duì)正交試驗(yàn)結(jié)果進(jìn)行分析。結(jié)果如表5和表6所示,已知F0.95(2,2)=19,F0.99(2,2)=99,F0.999(2,2)=999,對(duì)比得到各因素的影響顯著性,其中*的數(shù)量越多,表示影響越顯著。
由表5和表6分析可以看出:
1)傳統(tǒng)的極差分析法只能得到各因素對(duì)性能參數(shù)影響程度的相對(duì)大小,而方差分析法可以獲得各因素對(duì)發(fā)動(dòng)機(jī)性能影響的顯著性大小,結(jié)果直觀,與極差分析法相比,既能夠提高發(fā)動(dòng)機(jī)性能敏感性的分析準(zhǔn)確度,又可以加深設(shè)計(jì)人員對(duì)各因素實(shí)際影響程度大小的認(rèn)識(shí)。
2)對(duì)于發(fā)動(dòng)機(jī)推力和混合比,渦輪效率對(duì)性能的影響均高度顯著,即渦輪效率與發(fā)動(dòng)機(jī)性能參數(shù)有著極強(qiáng)的相關(guān)性,渦輪效率的改變會(huì)引起發(fā)動(dòng)機(jī)性能參數(shù)較大范圍的改變。
3)甲烷泵效率和氧泵效率對(duì)發(fā)動(dòng)機(jī)推力的影響均為高度顯著,對(duì)發(fā)動(dòng)機(jī)混合比的影響較顯著,但其顯著性明顯高于其他影響較顯著的因素。通過(guò)查表得F0.995(2,2)=199,即當(dāng)選擇的顯著性水平α=0.005時(shí),甲烷泵效率和氧泵效率對(duì)發(fā)動(dòng)機(jī)混合比的影響顯著性比其他因素要高出一個(gè)等級(jí)。
綜上,渦輪泵的效率對(duì)發(fā)動(dòng)機(jī)性能有著十分顯著的影響,其影響程度遠(yuǎn)大于其他因素,因此在設(shè)計(jì)與加工過(guò)程中,對(duì)渦輪泵的關(guān)鍵幾何參數(shù)要提高精度要求,降低加工散差,同時(shí)還需要通過(guò)渦輪氣流試驗(yàn)和泵水力試驗(yàn)或通過(guò)實(shí)際熱試車(chē)獲取渦輪泵的實(shí)際效率,并以此為基礎(chǔ),使用調(diào)節(jié)元件對(duì)發(fā)動(dòng)機(jī)進(jìn)行參數(shù)調(diào)節(jié),以保證發(fā)動(dòng)機(jī)性能參數(shù)滿足總體要求。
表5 推力敏感度方差分析Tab.5 Variance analysis of thrust sensitivity
表6 混合比敏感度方差分析Tab.6 Variance analysis of mixed ratio sensitivity
為了提高發(fā)動(dòng)機(jī)性能敏感性分析的準(zhǔn)確度,本文基于方差分析法對(duì)發(fā)動(dòng)機(jī)性能開(kāi)展了敏感性分析,分別得到了各因素對(duì)發(fā)動(dòng)機(jī)推力和混合比影響的顯著性大小,提高了發(fā)動(dòng)機(jī)性能敏感性分析的準(zhǔn)確度。
通過(guò)分析,發(fā)現(xiàn)渦輪泵的效率對(duì)發(fā)動(dòng)機(jī)性能的影響高度顯著。在發(fā)動(dòng)機(jī)設(shè)計(jì)與生產(chǎn)過(guò)程中,需要提高渦輪泵關(guān)鍵幾何參數(shù)的加工公差要求,并提高產(chǎn)品放行標(biāo)準(zhǔn),還要通過(guò)渦輪氣流試驗(yàn)和泵水力試驗(yàn)摸清每臺(tái)渦輪泵的實(shí)際特性,并使用調(diào)節(jié)元件對(duì)發(fā)動(dòng)機(jī)參數(shù)進(jìn)行調(diào)節(jié),以減小由于加工和裝配導(dǎo)致的渦輪泵性能散差對(duì)發(fā)動(dòng)機(jī)性能造成的影響。對(duì)于影響顯著的閥門(mén)流阻系數(shù),要通過(guò)水力試驗(yàn)獲得具體的參數(shù),若與設(shè)計(jì)值偏差較大,則需要在發(fā)動(dòng)機(jī)裝配過(guò)程中放置節(jié)流圈進(jìn)行相應(yīng)的調(diào)整,以保證發(fā)動(dòng)機(jī)性能的穩(wěn)定性。
本文使用的方差分析法普遍適用于各類液體火箭發(fā)動(dòng)機(jī),在之前的發(fā)動(dòng)機(jī)敏感性分析中很少用到,對(duì)將來(lái)其他發(fā)動(dòng)機(jī)的敏感性分析具有借鑒意義。