内弹道解算程序C++.docVIP

附录程序 #define _CRT_SECURE_NO_DEPRECATE #include stdio.h #include iostream #include conio.h #include math.h #include stdlib.h #include iomanip//C head file #include time.h #include malloc.h #include string #include fstream using namespace std; /*初始参数*/ //////////燃烧室内参数 //火药燃烧的参数 #define U1 0.2 //燃速系 #define an 0.82 //燃速指数 #define omga 1.76 //装药质量 #define am 21.76 //弹丸质量 #define f 980000 //火药力 #define F1 1.1 //阻力系数 #define STA 0.2 //（k-1） #define DP 1600 //火药装填密度 #define ALF 1.0e-3 //火药燃气余容 double PG=30e6 ; //挤进压力 //火药参数（单基9/7） #define HE1 0.5e-3 //第二种火药半弧厚 #define HD0 0.5e-3 //第二种火药内孔直径 #define HD2C 12e-2 //第二种火药火药长 //火药形状特征量Π1、Q1、β double XHT,XHQ,XHB; //double XHB; //碎粒燃烧结束时的Zb，火药形状参数χ、λ和μ double Zb0,CI,ALM,AMU; //碎粒燃烧阶段的形状参数χs、λs double CIS,ALMS; //////////火炮参数 #define V0 4.2e-3 #define omga2 0.34 //第一种药质量 #define f2 980000 //第一种药火药力 #define DP2 1600 //第一种药密度 #define ALF2 1e-3 //第一种药余容 //火药参数（选用杆状药） #define HE12 0.24e-3 //火药半弧厚 #define HD02 0.3e-3 //内孔直径 #define HD2C2 6.5e-3 //火药长 //#define f2 980000 double Z2,PSI2,dZ2,Zb2,Zm2; double DLT2; //火药形状参数χ、λ和μ double CI2,ALM2,AMU2; double XHB2; #define pi 3.1415926 //圆周率 #define atm 1.013e5 //大气压 #define h 1e-7 //步长 int i,i0,j,k,ii; //次要功系数、身管截面积，燃烧室当量长度 double FI,s,L0; //火药的装填密度 double DLT; double aaa; double NST,HIB; //火药燃烧时间 double TIME,TIMm,TIMb; //PTDD：膛底压力 double PTDD; //弹后空间平均压力、弹底压力 double P,PD; double PSI,ZJSN; //double PSI2; double PV,LPSI; long NN; //弹丸运动距离（不大于身管长度） double L,Lb,Lm,dL; //弹丸速度 double v,vm,vb,dv; int nn; //////////////微分过程参数 double a1[6]; double Z,dZ,Zb,Zm; void ranshaoshi(); void chushitiaojian(); void ranshao_0(); void ranshao_1(); void thy(); void main() { FILE *outhy; ranshaoshi(); NST=5; HIB=NST*h; outhy=fopen(hy0.txt,w); fprintf(outhy,omga(kg)\tFI\tDLT\n); fprintf(outhy,%f\t%f\t%f\n, omga,FI,DLT*1e-3)