系留氣球平臺動(dòng)力學(xué)仿真研究

第卷第12期計算機仿真2012年12月文章編號:1006-9348(2012)12-0066-05系留氣球平臺動(dòng)力學(xué)仿真研究楊方,張青斌,豐志偉(國防科技大學(xué)航天與材料工程學(xué)院湖南長(cháng)沙410073)摘要:研究系留氣球平臺設計問(wèn)題,隨著(zhù)人類(lèi)探索太空技術(shù)的不斷進(jìn)步,空間任務(wù)也變的越來(lái)越復雜。以空間繩網(wǎng)為柔性部件是最近發(fā)展起來(lái)的空間概念,在空間任務(wù)中發(fā)揮越來(lái)越大的作用,而柔性部件研究是當今世界學(xué)術(shù)研究的重點(diǎn)技術(shù),需要通過(guò)大量的試驗分析繩網(wǎng)的性能。根據繩網(wǎng)地面試驗的需要,建立了系留繩氣球平臺動(dòng)力學(xué)模型,分析了不同構建方式下的系留氣球平臺的穩定性,同時(shí)分析地面環(huán)境下風(fēng)力干擾作用對系留氣球平臺穩定性的影響,提出了構建穩定系留氣球平臺的方法,并將仿真結果與地面試驗結果進(jìn)行對比,結果一致關(guān)鍵詞:系留氣球;仿真;空氣動(dòng)力;地面實(shí)驗中圖分類(lèi)號:0313文獻標識碼:BDynamic Model and Analysis of Tethered Balloon SystemYANG Fang, ZHANG Qing-bin, FENG Zhi-weiCollege of Aerospace and Material Engineering, National Univ. of Defense TechnologyChangsha Hunan 410073, China)ABSTRACT: With the advance of explorating of space technology, space missions become more complex. Repre-sentatives of space webs technology as a flexible rope course is recently developed in the parts of the space concept inthe near future, and plays an increasing role in space missions. The academic research of flexible parts in the worldis a difficult problem that requires a lot of ground tests to analyse the performance of rope. Based on the rope test, weet up a dynamic model of the tethered balloon platform, analysed the stability of the different modes in different tethered balloon platform and the effects on the stability of the platform of the wind power in the ground environment, pro-posed the advise of building a tethered balloon platform approach, and tested the platform with the ground rope emission. The results are no differenceKEYWORDS: Tethered Balloon; Simulation; Aerodynamics; Ground Experiments等參數對高空系留氣球設計的影響得出風(fēng)場(chǎng)對系留氣球系1引言統的影響最大。史獻林,余莉施紅等利用有限元法建立隨著(zhù)空間技術(shù)的發(fā)展,對以空間繩網(wǎng)為代表的柔性在軌繩索模型,分析了系留氣球的繩孔過(guò)程。服務(wù)技術(shù)的需求越來(lái)越迫切。相對剛性結構而言,柔性部件對系留氣球的研究主要就是要建立一個(gè)穩定的空間繩結構小,質(zhì)量輕,設計靈活性高,在太空任務(wù)中能夠更好的完網(wǎng)地面試驗的平臺。本文利用集中質(zhì)量彈簧模型建立了系成對空間非合作目標的捕獲。為了研究繩網(wǎng)在太空的展開(kāi)留繩的仿真模型;主要考慮氣動(dòng)力的影響將氣球簡(jiǎn)化為剛抓捕的動(dòng)力學(xué)過(guò)程,分析繩網(wǎng)空間抓捕性能,同時(shí)降低成本,體建立氣球動(dòng)力學(xué)模型;對比不同構建模型,同時(shí)分析在自需要進(jìn)行合理的地面試驗。本文的目的就是在地面環(huán)境下然風(fēng)力作用下系留氣球系統的穩定性提出構建系留氣球平建立一個(gè)穩定的繩網(wǎng)地面發(fā)射試驗平臺。臺的方法,為繩網(wǎng)拋射試驗提供良好的試驗基礎。A. J. Euler等'對分析了高空系留氣球的可靠性,研究了在20km高空布置系留氣球的可行性;1996年,S.S.Bade.2系留氣球動(dòng)力學(xué)模型sha等分析了阻力系數下沉氣流、系留繩長(cháng)度和大氣環(huán)境2.1試驗模型繩網(wǎng)發(fā)V山中國煤化工氣球、三條固定繩索和有效載收稿日期:2012CNMHGE繩索的一段連接為懸置在空中熱回定在地面上,有效載荷懸掛在氣球的下方。2.2.1繩索的阻尼彈簧模型有限段法是分析繩索動(dòng)力學(xué)問(wèn)題的一個(gè)常用的方法,其主要途徑是采用離散化的方法將繩索離散為若干有限段,然后將每個(gè)繩段的質(zhì)量集中在兩個(gè)端點(diǎn),即繩段的結點(diǎn)處。在熱氣球離散化模型中,如果考慮繩索的彈性可伸長(cháng)的特點(diǎn)則可假A設繩段結點(diǎn)之間由“彈簧”相連,而且該“彈簧”只能承受張固定繩索力,不能承受壓力同時(shí)考慮到繩索的阻尼效應,故可以將繩段處理為質(zhì)量集中在繩段兩端點(diǎn)的集中質(zhì)量阻尼彈簧,這就是所謂的“集中質(zhì)量彈簧模型”。如圖3所示,在模型中將第n個(gè)繩段處理為質(zhì)量集中在端點(diǎn)的“彈簧阻尼”元件,并將其長(cháng)度記為ln,彈性系數記為kn阻尼系數記為Cn。圖1氣球系留繩試驗示意圖繩段如圖1,定義坐標系為:將氣球在地面的投影點(diǎn)O設置為原點(diǎn),A、BC分別為固定繩索在地面的固定點(diǎn),設定OA為X軸,垂直地面向上為Z軸,Y軸由右手坐標系確定。如圖2所示。定義跨度為地面固定點(diǎn)A、B、C到O點(diǎn)的距離圖3繩段的阻尼彈簧模型如圖3所示,建立第n個(gè)繩段單元的ψ=0連體坐標系ox,y2,該坐標系的三個(gè)坐標軸分別對應繩索曲線(xiàn)的主法向圖2繩索地面固定方式示意圖量、次法向量和切向量。連體坐標系與慣性坐標0xy24之間可以由三個(gè)歐拉角dn、6,、。來(lái)確定其轉換關(guān)系。慣性坐2.2繩索動(dòng)力學(xué)模型標系可以依次經(jīng)過(guò)三次旋轉與連體坐標系重合。首先繞坐基于有限元方法將柔性大變形的系留繩離散成一系標軸o轉過(guò)中角,從x到達OMz的位置;然后繞列首尾相連的繩段,各繩段采用集中質(zhì)量彈簧模型進(jìn)行模節線(xiàn)ON轉過(guò).角,從Oy到達OMy2的位置;最后繞擬。系留氣球簡(jiǎn)化為一個(gè)正球體,主要考慮重力,氣動(dòng)力以∞軸轉過(guò)ψ。角,從oNy2,到達0xy。繩段連體坐標系及系留繩牽引力影響的剛體模型。,yn2與慣性坐標系axay4=的轉換矩陣為cos , cosψ- sin cos 0, sinψ。 sin , cosψn+cosφ cos 6. sinψnsin6nsinψnB,=-cos sin v,-sin p, cos 8, cos v, -sin p, sin v,+ cos o, cos os y,n sin 8, os y.sin中nsin6n中sin由于繩索極度柔軟,不能承受法向的彎矩和剪切力僅和即第n個(gè)繩段的張力由下式給出能承受切向的張力。于是可以忽略歐拉角ψn,則兩坐標系間En≤0的變換矩陣可以簡(jiǎn)化為B. tsin d0在上式中pn(n)是第n段傘繩的線(xiàn)性張力函數,B。為B.--ndc0. cos cos 8.smn.(2)繩段的張力阻尼系數,應變n由下式得到o, sin 6o, sin e。cosb由于很難得到繩索準確的動(dòng)態(tài)應力一應變曲線(xiàn),在工程H中國煤化工上通常利用繩索的靜態(tài)應力一應變曲線(xiàn),采用分段線(xiàn)性化的在上后的長(cháng)度,即CNMH方法,將繩段的本構關(guān)系簡(jiǎn)化為線(xiàn)彈性與線(xiàn)性速度阻尼之+(y,-y)+(xn-n-1)2(5)C(8)上式中,F和F,分別表示繩段n法向力和切向力繩段p為周?chē)黧w的密度,d為繩索圓截面的直徑,CN和C分別是法向力系數和切向力系數,根據經(jīng)驗公式得到。2.2.3系統動(dòng)力學(xué)方程系留氣球假設為球體,質(zhì)量集中在球心。系留氣球的受力分析如圖6所示,主要有浮力F、重力G地面風(fēng)場(chǎng)作用力F.和系留繩拉力F1,上述外力可以表示為圖6。Fr圖4繩段坐標系對上式微分可以得到繩段的應變微分En=[(x-xn-1)(xn-xn1)+(y。-y-1)(yn-yn1+(z,,)(2,-Z,)J/(LAOL,)(6)2.2.2繩索氣動(dòng)力模型圖6系留氣球受力分析圖繩索所受到的外力包括重力、流體阻力等。本節主要討論在流體中運動(dòng)的繩段單元所受到的流體作用力。流體作F,=(p.-Pair)Vsg用力又可以細分為附加質(zhì)量力( Added Mass Force)、法向阻G=G+G, =(p.,+ m,)g力( Normal Drag)、切向阻力( Tangential Drag)等。式中p表示標準狀態(tài)空氣密度pm=1.209kg/m3,V為氣球的體積,g為重力加速度,取|g|=9.8m/s32。p4表示氮氣密度,m,表示氣球凈質(zhì)量。地面風(fēng)場(chǎng)作用力為F=-A(10)上式中C4為氣動(dòng)力系數、A為氣球的截面積、a為當地節點(diǎn)風(fēng)速系留氣球的動(dòng)力學(xué)方程可以寫(xiě)為F圖5流體中運動(dòng)的繩段如圖5所示,采用集中質(zhì)量法建立繩段n的氣動(dòng)力模F. +F+F型,并假設繩段為長(cháng)圓柱體。兩個(gè)端點(diǎn)分別記為i和j,質(zhì)心為G。繩段軸向的單位向量記為r,在來(lái)風(fēng)方向與繩段面內與T垂直的單位向量記為n。質(zhì)心G相對于流體的速度v3仿真分析可以表示為根據已經(jīng)建立的繩網(wǎng)發(fā)射平臺模型,模擬在不同條件下發(fā)射平臺的穩定性v=2(V+y)-V(7)3.1主要參數上式中V、V分別為節點(diǎn)i和節點(diǎn)j的速度,V為繩索周?chē)鶕K網(wǎng)拋射試驗的具體要求,選取直徑為7.8m的熱的流體速度。氣球將熱氣球假設為正圓球.豐要考慮其氣動(dòng)力和重力影作用在繩段n上的氣動(dòng)力表示為響。地面上中國煤化工點(diǎn)O,位于正三角F。=F+F形中心。主CNMHG表1系統參數從圖7~9中可以看出,在初始時(shí)刻,固定繩索處于松弛氦氣密度0.1347kg/m3狀態(tài)。在繩索繃緊過(guò)程中,氣球位置變化比較劇烈。在20s熱氣球凈重后,氣球位置開(kāi)始震蕩,最終趨于穩定。地面空氣密度1.29kg/m3氮氣重量33.47kg繩索質(zhì)量2kg/100m熱氣球直徑7.8--+--風(fēng)壓系數0.968氣球高度在仿真過(guò)程中,主要仿真發(fā)射平臺不同構建方式下在不同風(fēng)力情況下的穩定性。系留氣球平臺的第一種建構方式是氣球高度為90m,繩索固定點(diǎn)離O。點(diǎn)的距離為150m;第種是氣球高度為90m,繩索固定點(diǎn)離O點(diǎn)的距離為100m;第三種是氣球高度為90m,繩索固定點(diǎn)離O點(diǎn)的距離為50m地面風(fēng)速選擇4級(風(fēng)速6.7m/s)和六級(風(fēng)速12.8m圖8Y方向位移變化示意圖s)兩種情況進(jìn)行分析。主要可變參數如表2。其中,來(lái)風(fēng)方向定義為,在O-XY面內風(fēng)向與X軸的夾角。表2可變參數表項目范圍氣球高度M體熱氣球浮力200~1700N跨度50m.100m,150m風(fēng)速6.7m/s,12.8m/s來(lái)風(fēng)方向0~120°3.2仿真結果仿真分析中,主要觀(guān)察對象為氣球,通過(guò)觀(guān)察氣球的穩定性就可以判斷整個(gè)發(fā)射平臺的穩定性。在仿真模型中,H圖9Z方向位移變化示意圖點(diǎn)的位移變化情況就能夠完全反應系留氣球平臺的穩定為了研究地面風(fēng)場(chǎng)對系留氣球平臺穩定性的影響,本節情況。分析系留氣球平臺在來(lái)風(fēng)角度分別為0和90°的四級風(fēng)場(chǎng)氣球高度90m,繩索固定點(diǎn)離原點(diǎn)0距離150m。在四下的穩定性,仿真結果見(jiàn)表3和表4。從表中可以看出,四級級風(fēng)作用下仿真結果。其中,來(lái)風(fēng)方向為90°。風(fēng)場(chǎng)下系留氣球平臺的穩定性總體較好。表30°風(fēng)場(chǎng)作用下氣球位移情況平臺穩定性平臺平臺三D,漂移點(diǎn)/mAlAnNA圍3-416-1833-36D,漂移點(diǎn)/m范圍/m0D2漂移點(diǎn)/91.6范圍/m91.1下系面78-81主要分析在極限情況下系留氣球平臺的穩定性。如:北京年平均風(fēng)速在18-3m/s之間。風(fēng)速受地理環(huán)境的影響圖7X方向位移變化示意圖較大。城區中國煤化工霞云嶺為1.8ms;城區為2.CNMHG延慶縣、古北口都為3m/s;海拔100mn以上的佛爺頂年平均風(fēng)速最大為好。隨著(zhù)系留繩底端跨度增大,發(fā)射平臺的穩定性越好。5.7m/s。全年以春季風(fēng)速最大,冬季次之,夏季風(fēng)速最小。以北京氣象臺為例:4月份平均風(fēng)速3.4m/s,8月份平均風(fēng)參考文獻:速1.5m/s,相差一倍以上唯有延慶、昌平縣冬季風(fēng)速大于春1] A J Euler, SSBadesha, L D Schroeder. Very high altitude季風(fēng)速。其平均風(fēng)速都小于四級風(fēng)。hered balloon feasibility study[ R]. AIAA295216122CP, 1995表490°風(fēng)場(chǎng)作用下氣球位移情況表[2]SS Badesha, A J Euler, L D Schroeder. Very high altitude teth-ered balloon parametric sensitivity study R J平臺穩定性平臺一平臺平臺三AIAA29620579,1996漂移點(diǎn)/m-0.05[3] Paul Williams. Optimal Wind Power Extraction with a Tethered范圍/m0.3-0.213~-15Kie[R].AIAA2006-6193,2006D,漂移點(diǎn)/m[4]史獻林,余莉,施紅.系留氣球升空過(guò)程的動(dòng)態(tài)模擬[冂]·航空范圍/m1.5-2.59~1l學(xué)報,2009,30(14)D2漂移點(diǎn)/m[5]張晶敏,陸宇平.平流層氣球軌跡控制系統平衡點(diǎn)估算研究[J].航天控制,2007,25(3):428范圍/m92.5~9383.5~85.574~78[6 X Provot. Def ormat ion con st raint s in a mass spring model to de-進(jìn)行了系留氣球平臺的地面試驗,氣球漂移范圍在仿真scribe rigid cloth behavior[C]. Proceedin gs of Graphics Inter-結果范圍以?xún)?仿真結果與試驗一致。face,,1995,95:147-154[7 T Takagi, T Shimizu. Validity and layout of"NaLA": a net con小結figuration and loading analysis system[ J]. Fisheries Research本文主要是通過(guò)建立系留氣球系統的仿真模型,分析了2004,66(2/3):235-243.在不同情況下繩網(wǎng)發(fā)射平臺的穩定性[8]周雨青葉兆寧,吳宗漢.球類(lèi)運動(dòng)中空氣阻力的計算和分析1)系留氣球平臺的穩定性比較好,在風(fēng)力小于4級得情[冂].物理與工程,202,12(1):55-57況下可以用于繩網(wǎng)拋射試驗。[作者簡(jiǎn)介]2)系留氣球平臺在風(fēng)力小于4級,熱氣球靜載荷在楊方(1984-),男(漢族),山西省長(cháng)治市人,碩1500N以上時(shí),平面內位移在10m范圍以?xún)?。垂直方向穩定性相對比較穩定,其穩定性主要由熱氣球的靜載荷和繩索強士研究生,主要研究領(lǐng)域為動(dòng)力學(xué)與控制。張青斌(1975-),男(漢族),山西省大同市人,副度決定;教授,碩士研究生導師,主要研究領(lǐng)域為動(dòng)力學(xué)與3)不同的來(lái)風(fēng)方向對系留氣球平臺的穩定性有一定影控制。,但影響不大。隨著(zhù)氣球凈載荷的增加,氣球的穩定性越豐志偉(1984-),男(漢族),山西省臨汾市人,博士研究生,主要研究領(lǐng)域為動(dòng)力學(xué)與控制。(上接第52頁(yè))[2]許菁,雷定猷,鄧煜陽(yáng).基于區位理論的物流配送中心車(chē)輛調調度的效率與傳統算法基本一致。通過(guò)對表2中的數據進(jìn)度優(yōu)化算法[J.鐵道貨運,2007-9行整理分析能夠得知,假設軍用車(chē)輛運輸路徑阻塞,那么利3]郎茂祥,胡思繼車(chē)輛路徑問(wèn)題的禁忌搜索算法研究[J].管用本文算法進(jìn)行車(chē)輛調度的效率高于傳統算法。理工程學(xué)報,2004-1[4]李軍,郭強.車(chē)輛調度問(wèn)題的改進(jìn)表上作業(yè)法[冂].西南交通5結束語(yǔ)大學(xué)學(xué)報,2000-5本文提出了一種基于蟻群信息反饋算法的救災物資軍5]宋維堂張淑梅.車(chē)輛調度優(yōu)化模型仿真研究[冂].計算機仿真,201-9:346-349用車(chē)輛調度方法。計算路段阻塞系數,并將其反饋到調度端,利用蟻群極值方法選取最優(yōu)路徑。從而實(shí)現救災物資軍[作者簡(jiǎn)介]用車(chē)輛合理調度。實(shí)驗證明,這種算法提高了軍用車(chē)輛調度李冬梅(1977-),女,(漢族),山東梁山人,博士,的效率。工程師,主要研究方向:后勤裝備建設。嚴代彪(1973-),男(漢族)重慶梁平人,博士,工參考文獻白化建設。[1]李寧陳彬徐凱一種基于道路網(wǎng)絡(luò )拓撲改進(jìn)的格網(wǎng)空間索中國煤化工引算法[J].上海師范大學(xué)學(xué)報(自然科學(xué)版),2008-5CNMHG