内弹道解算程序C.doc

附录程序
#define _CRT_SECURE_NO_DEPRECATE
#include <>
#include <iostream>
#include <>
#include <>
#include <>
#include <iomanip>//C head file
#include <>
#include <>
#include <string>
#include <fstream>
using namespace std;
/*初始参数*/
//////////燃烧室内参数
//***燃烧的参数
#define U1 //燃速系
#define an //燃速指数
#define omga //装药质量
#define am //弹丸质量
#define f 980000 //***力
#define F1 //阻力系数
#define STA //(k-1)
#define DP 1600 //***装填密度
#define ALF -3 //***燃气余容
double PG=30e6 ; //挤进压力
//***参数(单基9/7)
#define HE1 -3 //第二种***半弧厚
#define HD0 -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 -3
#define omga2 //第一种药质量
#define f2 980000 //第一种药***力
#define DP2 1600 //第一种药密度
#define ALF2 1e-3 //第一种药余容
//***参数(选用杆状药)
#define HE12 -3 //***半弧厚
#define HD02 -3 //内孔直径
#define HD2C2 -3 //***长
//#define f2 980000
double Z2,PSI2,dZ2,Zb2,Zm2;
double DLT2;
//***形状参数χ、λ和μ
double CI2,ALM2,AMU2;
double XHB2;
#define pi //圆周率
#define atm //大气压
#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;