内弹道程序报告姓名:学号:专业:武器发射工程学院:能源与动力工程学院指导教师:1. 程序功能对于一般的内弹道方程组,微分方程是非线性的,编制通过数值解法进行求解的程序。2. 设计思路和整体方案(1)输入已知的数据,包括火炮构造,弹丸诸元,装药条件,起始条件,步长等参数。(2)进行一些常量的计算,方便在程序中进行使用。(3)计算内弹道方程中变量的初值。(4)将内弹道方程组进行量纲化,然后使用龙格库塔法计算微分方程,包括最大压力搜索和特征点判断,得出各个过程的数据。(5)将程序输出成表格和曲线。3. 程序遇到的主要问题和解决方案(1)在程序中首先要编制一个四阶的龙格库塔子函数,同时还要编制一个微分方程右端子函数。在龙格库塔子函数编制过程中借鉴了之前导弹运动弹道程序的编制。voidrk(n,h)intn;doubleh;{externvoiddery();doublea[4],old_y[6],Y1[6],*dy;inti,j;dy=calloc(n,sizeof(double));a[0]=a[1]=h/2;a[2]=a[3]=h;dery(n,dy,Y);for(i=0;i