内弹道解算程序c++.doc

附录程序#define_CRT_SECURE_NO_DEPRECATE#include<>#include<iostream>#include<>#include<>#include<>#include<iomanip>//Cheadfile#include<>#include<>#include<string>#include<fstream>usingnamespacestd;/*初始参数*///////////燃烧室参数//***燃烧的参数#//燃速系#//燃速指数#//装药质量#//弹丸质量#definef980000//***力#//阻力系数#//(k-1)#defineDP1600//***装填密度#-3//***燃气余容doublePG=30e6;//挤进压力//***参数(单基9/7)#-3//第二种***半弧厚#-3//第二种***孔直径#defineHD2C12e-2//第二种******长//***形状特征量Π1、Q1、βdoubleXHT,XHQ,XHB;//doubleXHB;//碎粒燃烧结束时的Zb,***形状参数χ、λ和μdoubleZb0,CI,ALM,AMU;//碎粒燃烧阶段的形状参数χs、λsdoubleCIS,ALMS;//////////火炮参数#-3#//第一种药质量#definef2980000//第一种药***力#defineDP21600//第一种药密度#defineALF21e-3//第一种药余容//***参数(选用杆状药)#-3//***半弧厚#-3//孔直径#-3//***长//#definef2980000doubleZ2,PSI2,dZ2,Zb2,Zm2;doubleDLT2;//***形状参数χ、λ和μdoubleCI2,ALM2,AMU2;doubleXHB2;#//圆周率#//大气压#defineh1e-7//步长inti,i0,j,k,ii;//次要功系数、身管截面积,燃烧室当量长度doubleFI,s,L0;//***的装填密度doubleDLT;doubleaaa;doubleNST,HIB;//***燃烧时间doubleTIME,TIMm,TIMb;//PTDD:膛底压力doublePTDD;//弹后空间平均压力、弹底压力doubleP,PD;doublePSI,ZJSN;//doublePSI2;doublePV,LPSI;longNN;//弹丸运动距离(不大于身管长度)doubleL,Lb,Lm,dL;//弹丸速度doublev,vm,vb,dv;intnn;//////////////微分过程参数doublea1[6];doubleZ,dZ,Zb,Zm;voidranshaoshi();voidchus