資源描述:
《參數(shù)化設(shè)計(jì)變量》由會(huì)員上傳分享,免費(fèi)在線閱讀,更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫(kù)。
1、海納仿真交互平臺(tái)Tel:4008803423項(xiàng)目任務(wù)書注:請(qǐng)認(rèn)真填寫表1,分析要求將會(huì)作為驗(yàn)收標(biāo)準(zhǔn)。為了有效保障客戶切身利益,所有擔(dān)??铐?xiàng)均通過平臺(tái)監(jiān)管,待項(xiàng)目完成后支付平臺(tái)供應(yīng)商。表1(客戶填寫)*項(xiàng)目名稱耦合內(nèi)彈道方程彈丸運(yùn)動(dòng)動(dòng)網(wǎng)格仿真*項(xiàng)目周期中標(biāo)后10個(gè)工作日內(nèi)*分析要求(驗(yàn)收標(biāo)準(zhǔn))1、將我給的內(nèi)彈道方程耦合進(jìn)udf,且udf能成功導(dǎo)入fluent并實(shí)現(xiàn)彈丸運(yùn)動(dòng)2、儲(chǔ)存彈丸運(yùn)動(dòng)每個(gè)時(shí)刻的數(shù)據(jù),可以看到炮口流場(chǎng)的氣體速度、壓力、密度分布圖;以及彈丸前部壓力圖。*其他要求:*項(xiàng)目備注:① 供應(yīng)商報(bào)價(jià)
2、將根據(jù)以上填寫內(nèi)容評(píng)估。② 以上填寫內(nèi)容將作為最終驗(yàn)收標(biāo)準(zhǔn),請(qǐng)仔細(xì)填寫。③ 雙方溝通不順利等相關(guān)事宜,可找平臺(tái)項(xiàng)目主管對(duì)項(xiàng)目執(zhí)行情況進(jìn)行評(píng)估,協(xié)調(diào)雙方達(dá)成一致意見。④ 為了更好的為客戶服務(wù),待項(xiàng)目完成后,平臺(tái)供應(yīng)商均免費(fèi)提供售后一個(gè)月的技術(shù)支持服務(wù),請(qǐng)您務(wù)必在服務(wù)期滿后為該供應(yīng)商打分,督促我們提供更優(yōu)質(zhì)的服務(wù)。方案情況:采用2d軸對(duì)稱非結(jié)構(gòu)網(wǎng)格,如下圖:藥室、炮管網(wǎng)格尺寸為1x5較小,y100mm以下網(wǎng)格尺寸為1x10稍大,y100mm-y2500mm區(qū)域尺寸為2x10,尺寸較大。海納仿真交互平臺(tái)Tel
3、:4008803423模型采用k-epsilon模型,參數(shù)可以默認(rèn)。動(dòng)網(wǎng)格要求彈丸運(yùn)動(dòng)為x方向平動(dòng),速度由udf確定;藥室部分氣體壓力設(shè)置初始為30Mpa,炮膛和大氣初始為大氣壓。最后結(jié)果類似如下:海納仿真交互平臺(tái)Tel:40088034231、炮口速度、密度圖2、彈前壓力圖海納仿真交互平臺(tái)Tel:4008803423內(nèi)彈道運(yùn)動(dòng)c語言源程序:#include"studio.h"#include"math.h"/**********************************/#defineroup1
4、600//火藥參數(shù)常數(shù)賦值#definew1.16#definef950000#definesita0.25#definealpha0.001#defineu11.6*10^-8#definen10.83#definee10.00055#defined0.00055#definex10.75#definelamada0.12#definemiu0.0#defineXs1.696#definelamadas-0.4104#definem0.389#defines0.00266#defineV00.0007
5、1#definelg1.502#definefai1.168#defineP03*10^7#definededa828#defineL00.567669#defineVj1641.857#defineB2.3930702#definelgpingjun6.384海納仿真交互平臺(tái)Tel:4008803423realYJ[6],YK[4];//數(shù)組表示內(nèi)彈道參數(shù)intmain()//主函數(shù)YJ[0]=0,YJ[1]=0,YJ[2]=0.0274047,YJ[3]=0,YJ[4]=0.014135,YJ[5]
6、=0.0187976;//分別表示相對(duì)時(shí)間,相對(duì)速度,相對(duì)壓強(qiáng),炮彈相對(duì)位移,相對(duì)燃燒量,相對(duì)燃燒厚度YK[0]=0,YK[1]=0,YK[2]=2*pow(10,7),YK[3]=0;//分別表示時(shí)間,速度,壓強(qiáng),炮彈位移voidrk(intn,doubleh);do{rk(6,0.001);//時(shí)間步長(zhǎng)為0.001s}while(YK[3]<1.502);//Z[3]表示炮彈出炮口后停止,截止條件//龍格庫(kù)塔算法(積分求彈丸速度、彈后壓力等參數(shù))voidrk(intn,doubleh){extern
7、voiddery(intn,doubledy[],doubleYJ[]);doublea[4],old_Y[6],*dy,Y1[6];inti,j;dy=(double*)calloc(n,sizeof(double));a[0]=a[1]=h/2;a[2]=a[3]=h;dery(n,dy,Y);for(i=0;i8、[j+1]*dy[i]/3;}dery(n,dy,Y1);}for(i=0;i