内弹道算课程论文.doc

1、沈阳理工大学课程设计说明书 摘 要本文简单介绍了内弹道的计算方法，内弹道解法是内弹道理论的核心，主要包括解析法、经验法、表解法和图解法等，随着计算机在弹道学的广泛应用，可以建立更精确的数学模型，通过数值方法可以获得膛内压力、弹丸速度等参量变化规律的精确描述。本文主要讨论了内弹道方程组的建立过程和数值解法，通过计算弹道的特殊点来分析弹丸在膛内的运动状况，通过计算和matlab综合仿真得到相应的弹道曲线。关键词：内弹道；初值计算；数值解法；matlab仿真目 录摘 要I1 引 言11.1 加农炮的定义及其发展史11.2 自行火炮的国内外发展趋势11.3 1983式152毫米加农炮性能22 内弹道定

2、义52.1 内弹道简介52.2 研究对象53 内弹道方程组的建立73.1 基本假设73.2 内弹道方程组84 内弹道方程组的建立过程和解法94.1 内弹道方程组的分析94.2 初始条件计算公式和火药诸元94.2.1 内弹道方程组所用到的初值的求解104.3 对内弹道方程组求解114.3.1热力学第一时期114.3.2热力学第二时期144.4 仿真和计算15结 论17附 录22结束语24参考文献25281 引 言1.1 加农炮的定义及其发展史加农炮是指发射仰角较小，弹道低平，可直瞄射击，炮弹膛口速度高的火炮。常用于前敌部队的攻坚战中。严格上说，坦克炮也属于加农炮。榴弹炮则是发射仰角大，弹道较高而

3、弯曲，不能直瞄射击而炮弹能飞越障碍物攻击目标的火炮，一般射程较加农炮远，常用于第二线的阵地上对最前线的火力支援和对敌阵的火力压制。军事上没有“普通”炮的概念，加农炮和榴弹炮都属于很“普通”的火炮的两大不同种类，现代军事装备上还有将加农炮和榴弹炮的功能合而为一的火炮加农榴弹炮，简称加榴炮。加农炮是一种身管较长、弹道平直低伸的野战炮，它最早起源于14世纪，到16世纪时，欧洲人便开始把这种身管较长的炮称之为加农炮，当时身管长为1622倍口径。18世纪，身管长一般为2226倍口径。二次大战前后，口径在105108毫米之间的加农炮得以迅速发展，炮身长一般为3052倍口径，初速达880米/秒，最大射程30

4、公里。世纪 60年代，炮身长为4070倍口径，初速达950米/秒，最大射程达35公里。60年代以后，加农炮基本没研制新型号，性能仍保持在60年代水平，身管长、弹道低伸的火炮。二次大战前后，口径在105108毫米之间的加农炮得以迅速发展，炮身长一般为3052倍口径，初速达880米/秒，最大射程30公里。世纪60年代，炮身长为4070倍口径，初速达950米/秒，最大射程达35公里。60年代以后，加农炮基本没研制新型号，性能仍保持在60年代水平。加农炮按口径可分为：小口径加农炮，75毫米以下；中口径加农炮，76130毫米；大口径加农炮，130毫米以上。按运动方式可分为：牵引式、自运式、自行式和装载到

5、坦克、飞机、舰艇上载运式4种。反坦克炮、坦克炮、高射炮、航空炮、舰炮、海岸炮均属加农炮之类。加农炮炮管长度一般为4070倍口径，所以射程较其它类型的火炮都远，例如，美国175 毫米自行加农炮，最大射程32.7公里；而口径比它大的203毫米榴弹炮，最大射程却 只有29公里。1.2 自行火炮的国内外发展趋势俄罗斯 上世纪80年代，俄罗斯出口了大量122mm和152mm自行火炮，最近几年出口量寥寥无几。152mm 2S19(MSTA-S)是在俄陆军服役的最新自行火炮，并供出口。专供出口的155mm/52倍口径2S19自行火炮装有ASUNO自动火炮火控系统，目前处于原型阶段，整个炮塔可安装在其它履带式

6、底盘上。 波兰 选用英国皇家军械厂防御分部155mm/52倍口径AS90炮塔，满足未来的要求。头两个炮塔现已安装在OBRUM工厂设计制造的底盘上。整个火炮系统叫“克拉布”，HSW公司是全权主承包商。 斯洛伐克 上世纪70年代后期，ZTS公司利用“塔特拉”装甲88卡车底盘开发了152mm“达纳”自行火炮系统，已为国内和出口建造了大约750辆。以后在“达纳”基础上发展了“祖扎娜”火炮系统，安装155mm/45倍口径炮管。目前正在研究155mm/52倍口径型。该公司开发的A40炮塔系统已安装在T-72M1坦克底盘上进行试验。 以色列 索尔塔姆公司是该国唯一的火炮系统制造商，开发最新型ATMOS 20

7、00自行火炮系统供出口，印度是潜在客户。 AIMOS 2000采用改进型“塔特拉”66轮式底盘，装有155mm/52倍口径火炮系统，也可安装39倍或45倍口径的炮管。并装有先进火控系统，提高打击目标的精度。该公司还建造可空运的“拉斯卡尔”和“斯拉默”履带式155nm/52倍口径自行火炮样机。后者采用“梅卡瓦”主战坦克底盘，配备60发弹药，另有15发备用。 伊朗 已开发生产155mm“拉德”-2和122mm“拉德”-1两种自行火炮系统。前者与美制M109相似，装有155mm/39倍口径炮管，发射底排弹时最大射程24km，战斗重量36t。标准型安装一台840马力的V-84MS柴油机，“拉德”-2M

8、的动力是一台700马力的5TDF柴油机。后者采用改进的“波拉格”履带装甲人员输送车底盘，装有外观与俄罗斯122mm 2S1自行火炮相似的炮塔，发射标准122mm弹药，最大射程15.2km。 新加坡 是部署牵引式155mm/52倍口径FH-2000火炮系统的第一个国家。新加坡技术动力学公司的产品有155mm FH-88和FH-2000。其开发的首件产品155mm/39口径44轮式轻型自行榴弹炮现正进行射击试验。该炮约重7t，可用C-130“大力神”运输机空运，或者用CH-47或CH-53直升机吊运。该公司还在为新加坡武装部队开发新型履带式自行155mm/39倍口径自行火炮，采用美国联合防御公司的

9、改进型M109底盘，已进行射击试验。 韩国 美国联合防御公司的155mm M109A6是韩国陆军标准的自行火炮系统，与三星技术公司合作生产1040辆。三星公司新近开发生产155mm/52倍口径K9“雷声”自行火炮系统，载弹量48发，炮弹自动装填，手工装药。土耳其在德国拒供PzH2000火炮系统后决定选择该型火炮。该公司还利用相同的底盘开发出弹药补给车。 南非 155mm/45倍口径G6 66轮式自行火炮系统是南非军队唯一装备的自行火炮系统，已装备53辆，其它非洲国家也有装备，且出口阿曼24辆、阿联酋78辆。G6的平台和弹药仍在继续开发，以提高射速和改进火炮的布置。该公司的新型155mm/52倍

10、口径T6炮塔已在印度T-72坦克底盘上演示，“康达尔”轮式155mn/52倍口径火炮系统是参与印度轮式自行火炮的竞争对象。1.3 1983式152毫米加农炮性能 装备陆军集团军属炮兵师，用来压制和歼灭敌有生力量，通信联络系统和后方设施，击毁敌主战坦克、装甲车辆和自行火炮，破坏敌野战防御工事。作为海岸炮，可消灭敌水上目标。全炮重10吨左右。身管采用低镍合金钢，进行自紧处理，装单室冲击式炮口制退器；炮闩为半自动立楔式，配有闩柄保险机构；采用炮身回转结构，炮身可回转180度进行牵引，缩短行军长度；采用气胎式车轮；座盘为分离式，行军时平放在大架架尾部平面上，大架并架时自行紧固；采用独立轮转式前车机构及

11、倒车杠杆。主要技术特点表现为：一、初速高，达到了每秒955米，这也是实现远射程的前提。二、射程远，最大射程达到30公里，这已经超出了中国以往火炮的射程。三、射速快，每分钟为34发，在当时的大口径火炮中，这样的射速已是比较快的了。此外，该炮的高低射界为-2.5度+45度，方向射界左24度，右26度，行军战斗转换时间34分钟，运动方式可以用汽车，也可以用履带车牵引。2 内弹道定义2.1 内弹道简介内弹道（internal ballistics）是弹道的一部分，内弹道研究弹丸从点火到离开发射器身管的行为。内弹道学研究对各种身管武器都有重要意义。击发方法：任何类型的身管武器第一步需要击发火药。最早的枪

12、支、大炮由一个一端密封的金属管组成。 内弹道学是研究发射过程中枪炮膛内及火箭发动机内的火药燃烧、物质流动、能量转换、弹体运动和其它有关现象及其规律的弹道学分支学科。燃烧的发射药产生具有很高压力的气体，使弹丸加速穿过炮膛，直到以预定初速离开炮口。初速是具有一定质量和形状的弹丸最终要达到的整个射程的基础。在设计火炮时必须进行计算以保证最正常、最有效地产生所需要的初速。发射装药产生的能量用于完成好几种工作。大部分能量用于赋予弹丸速度。能量还消耗在做下述功上：使弹丸旋转，克服弹丸与膛壁之间的摩擦力，使发射药和发射药气体在膛内运动以及使火炮后坐部分后坐。有些能量还以热能的形式损失在身管、炮尾、弹丸和药筒

13、（如果使用药筒的话）上。 发射过程都是从点火开始，通过机械击发、电热或其他方式将点火药点燃，所产生的高温气体及灼热粒子再点燃火药装药，迅即扩展到整个装药表面，并同时沿着药粒厚度向内层燃烧。燃烧进行在一个封闭的空间中，这个空间前由弹丸的弹带封闭，后有火炮所采用的紧塞装置封闭，紧塞装置用于防止火药气体从后面逸出。在发射药气体的压力达到能使弹丸运动的程度之前，发射药的燃烧速度与膛压增加的速度是成正比例的。所谓“弹丸启动压力”就是指使弹丸开始向前运动的压力。当弹丸沿身管向前运动时，供发射药气体占用的空间增大，因此膛压的增加速度减小。当空间增加所导致的压力的增加相等时，膛压达到最大值。自此以后膛压开始下

14、降，同时弹丸却在继续加速，甚至在发射药全部燃尽后弹丸仍在继续加速，只是加速度逐渐减小，弹丸一出炮口即变为减速。2.2 研究对象内弹道学的研究对象，主要是有关点火药和火药的热化学性质，点火和火药燃烧的机理及规律；有关枪炮膛内火药燃气与固体药粒的混合流动现象，有关弹带嵌进膛线的受力变形现象，弹丸和枪炮身的运动现象；有关能量转换、传递的热力学现象和火药燃气与膛壁或发动机之间的热传导现象等。 弹丸在膛内的运动大约要消耗掉发射药产生的能量的25-35%。其余的能量都在弹丸离开炮口后排入大气。通过增加身管长度以延长发射药气体作用于弹丸时间的方法，还有可能使弹丸初速增加。只是用这种方法增加初速也有其缺点，因

15、为在身管增长超过一定限度后所增加的初速与所带来的缺点相权衡，是得不偿失的。从发射药燃尽点开始，弹丸速度的增加是越来越平缓的。 内弹道学主要从理论和实验上对膛内的各种现象进行研究和分析，揭示发射过程中所存在的各种规律和影响规律的各有关因素；应用已知规律提出合理的内弹道的方案，为武器的设计和发展提供理论依据；有效地利用能源及探索新的发射方式等。 利用所掌握的内弹道规律，改进现有的发射武器和设计出新型的发射武器，这是内弹道设计 的研究内容。它是以内弹道方程组为基础的 ，例如根据战术技术要求所给定的火炮口径，及外弹道设计所给出的初速、弹重等主要起始数据，解出合适的炮膛结构数据、装填条件，以及相应的压力

16、和速度变化规律。 在内弹道设计方案确定之后，方案的数据就是进一步进行炮身、炮架、药筒、弹丸、引信及发动机等部件设计的基本依据。因此，发射武器的性能在很大程度上决定于内弹道设计方案的优化程度。 能源是实现内弹道过程的主要物质基础，如何选择合适的能源，有效地控制能量释放规律，合理地应用释放的能量以达到预期的弹道效果，一直是内弹道学研究的一个主要问题。 火药是最常用的主要能源。早在无烟药开始应用时对于成形药粒的燃烧，就采用了全面着火、平行层燃烧的假设，并以单一药粒的燃烧规律代表整个装药的燃烧规律，称为几何燃烧定律。它是内弹道学的一个重要理论基础。长期以来，应用这个定律指导改进火药的燃烧条件，控制压力

17、变化规律，以达到提高初速和改善弹道性能的目的。 在火炮设计中发射药在膛内的燃尽位置很重要。如果燃尽位置在膛内过于靠前，则很可能会增加耀眼的炮口焰，从而增加被敌人发现的可能性。如果燃尽位置在炮口外，则炮闩在发射药全部燃尽前有被打开的危险。在设计火炮及其装药系统时，必须非常注意这种可能性，特别是对发射后自动开闩的火炮。使燃尽位置适当靠后还有其他一些理由，其中比较重要的是，这样做能减小各发弹之间的初速差异。3 内弹道方程组的建立3.1 基本假设在写出内弹道方程组时采用了以下假设： （1）火药气体的流动是零维的、无黏性的和不可压缩的； （2）火药固体和气体混合物可由Nobel-Abel状态方程描述；

18、（3）火药燃烧服从几何燃烧定律； （4）可以使用药粒的平均尺寸（长度、半径等）来描述药粒的实际几何尺寸，并假定所有药粒具有相同大小和外形，对于多孔火药，认为孔径是均匀分布的，而且孔对应的所有弧厚都是均匀相等的； （5）在t=0时刻，所有药粒同时着火； （6）在任一瞬间，单位质量火药固体分解后所释放的能量数量都是在当前平均气温下进行的；（7）不考虑火药气体混合物主要成分的再分解； （8）通过火炮身管表面的热量损失可以根据火药燃烧所释放的总能量来算，它占火药燃烧总放热量的百分比可以用系数K。来表示；（9）弹带与炮膛形成了一个完全的气体密封；（10）膛内气流边界层效应可以忽略不计； （11）可以利用

19、Lagrange问题的解来建立起平均压力、膛底压力与弹底压力之间的相互关系； （12）火药气固昆合物在膛内是均匀分布的； （13）绝热系数O=k-l是一个常数； （14）火药气体混合物是一种机械混合； （15）火药气体混合物的物理特性，可由机械混合的相应公式计算； （16）弹丸在膛内运动时期，弹丸前端所受的空气阻力可以忽略不计； （17）火药燃气生成物成分保持不变； （18）不考虑火药的侵蚀燃烧。【3】3.2 内弹道方程组根据以上假设，单一装药内弹道学方程组归纳如下：（1）形状函数： （2）燃速方程： （公式 3.2.1）（3）弹丸运动方程： （4）内弹道基本方程： 其中 （公式3.2.2 ）

20、 弹丸速度与行程关系式 （公式3.2.3）次要功计算系数 （公式3.2.4 ）装药密度 （公式3.2.5 ） 药室容积缩颈长 （公式3.2.6 ）根据射击现象的特点将射击过程划分为三个不同的阶段，即前期、第一时期和第二时期。在这三个不同阶段之间又是互相联结的，前期的最终条件就是第一时期的起始条件，而第一时期的最终条件又是第二时期的起始条件。因此，对于这三个阶段就应该根据各阶段的特点，按顺序地作出各阶段的解法。【1】4 内弹道方程组的建立过程和解法4.1 内弹道方程组的分析内弹道方程组中共有P、v、l、t、和Z六个变量，其它各量都是已知常量，有五个独立的方程，如取其中一个变量为自变量，则其余五个

21、变量作为自变量的函数，可以从上述方程组中解出，方程组是封闭的。【1】 在选择自变量时，我们应以自变量是否有已知的边界条件作为选择的主要标准。在第一时期的所有变量中，只有及Z这两个变量的边界条件是已知的，即从0到l，Z从Z0到l。从数学处理来讲，选择Z作为自变量比选择方便。因此，在现有的弹道解法中大多是采用Z作为自变量。不过在具体解方程组时。由于z的起始条件Z0同Z总是以Z-Z0的形式出现，所以令x=Z-Z0。则所解出的各变量都将以x的函数形式来表示。分析解法 ：从弹道方程组利用数学解析的方法，直接或者间接解出 P=P(l)、v=v(l)、P=P(t)和v=v(t) 的函数关系。表解法 ：在一定

22、的条件下预先将弹道解编成数值表，应用时只需要经过简单的运算和查表就可以求得弹道解。计算机解法：通过计算机编程求弹道解。在射击过程中，弹丸向前运动，弹后空间不断增加，因此膛内压力是弹后空间容积的函数。4.2 初始条件计算公式和火药诸元多孔火药然山哦到腰孔彼此想切时的相对气体质量 （公式4.2.1） 火药燃烧结束时相对厚度 （公式4.2.2 ）当多孔火药燃烧到药孔彼此相切的瞬间棱柱的半径（公式4.2.3） 极限速度 （公式4.2.4 ）符号数值单位符号数值单位符号数值单位火药密度1600kg/m3炮膛断面积S0.019m2火药形状特征量0.704装药量21.09kg药室容积V00.03057m3火

23、药形状特征量0.242火药力f950kJ/kg弹丸全行程长lg7.091m火药形状特征量-0.0245火药热力参数0.25燃速系数u15.19210-8分裂后形状特征量1.225火药气体余容0.001m3/kg燃速指数n0.82分裂后形状特征量-0.3阻力系数1.021/2火药厚度e11.1mm弹丸质量m48kg药室实际长度0.969m火药内孔直径d1mm次要功计算系数1.196挤进压力p03107Pa火药外径D011.8mm火药孔数7单一装药表4.2.1 火药诸元4.2.1 内弹道方程组所用到的初值的求解利用matlab编程实现了初值的计算，通过运行程序输入火药诸元得到的初值大小如下 t01

24、 =0P01 =5.7738L0 =1.6089vj =52.8395B =17.4901=5.5913Zk =5.0840e+0030 = 2.0275Z0 =3.44974.3 对内弹道方程组求解4.3.1热力学第一时期弹丸是瞬时挤进膛线，并在压力达到挤进压力P0时才开始运动。这一时期的特点应该是定容燃烧时期，此时在这一时期中，火药在药室容积W0中燃烧，压力则由PB 升高到P0，与P0相应的前期结束的瞬间标志火药形状尺寸的诸元也将相应地为0、0及Z0。这些量既是这一时期的最终条件，又是第一时期的起始条件。所以，这一时期解法的目的，实际上就是根据已知的P0分别解出0、0及Z0这三个前期诸元。

25、首先根据定容的状态方程解出0 ： 忽略PB，则求得了0后，应用所给出的及Z的公式分别计算出0及Z0 火药相对燃烧面积 （公式4.3.1） （公式4.3.2）求出了这三个诸元之后，即可以作为起始条件进行第一时期的弹道解。 （1）速度的函数式将燃速方程和弹丸运动方程联立消去Pdt （公式4.3.3）从起始条件v=0及Z=Z0积分到任一瞬间的v及Z （公式4.3.4）x=Z-Z0， （公式4.3.5） （公式4.3.6）该式表明，在一定装填条件下，弹丸速度与火药的已燃厚度成比例。（2）解火药的已燃部分的函数式 将Z=x+Z0代入形状函数中导出由于 （公式4.3.7） （公式4.3.8） （公式4.3

26、.9）从而导出 （公式4.3.10）(3)解弹丸行程的函数式 （公式4.3.11）将弹丸运动方程和内弹道基本方程联立消去SP得 (公式4.3.12）根据药室自由容积缩径长 （公式4.3.13） （公式4.3.14） 装填参量 （公式4.3.15）从而求得以下弹丸行程函数 （公式4.3.16）（4）压力函数式 从内弹道学基本方程可以得出 （公式4.3.17）将前三式代入有 （公式4.3.18）（5）燃烧结束瞬间的各弹道诸元 燃烧结束点的各弹道诸元既是第一时期的最终条件，又是第二时期的起始条件，所以燃烧结束点的诸元是必须计算出来的。在火药燃烧结束瞬间的条件为 （公式4.3.19）因此可列出、 及

27、各诸元的表达式为： （公式4.3.20）4.3.2热力学第二时期 （1）速度函数式在第二时期中，由于火药已经燃完，不再有火药燃烧的现象，因此这一时期的基本方程组为 （公式4.3.21）在这个方程组中，有v、l及P三个变量。为了解出这些变量的函数关系，必须指定其中一个变量作为自变量。由于这一时期是从燃烧结束点一直到炮口，所以就起始条件而言，这三个变量的起始条件都是已知的。但是就最终条件而言，只有l是已知的，即所谓弹丸全行程长lg。显然，在这种情况下，选择l作为自变量是恰当的，把v和P作为l的函数来表示。将以上两个方程消去SP，得到如下的微分式 （公式4.3.22）从而可以进行如下的积分 （公式4

28、.3.23）极限速度 （公式4.3.24）则速度函数是 （公式4.3.25）（2）压力的函数式 求出了之后，将给定的l所求得的v分别代入内弹道基本方程，即求得相应的压力 （公式4.3.25）再根据燃烧结束点公式炮口点压力【2】我们即求得第二时期的P-l及v-l曲线，再加上第一时期的P-l及v-l曲线，从而求得整个的P-l及v-l曲线。 4.4 仿真和计算最大膛压点代入数据可得=328080(kPa)时，求出对应的时间的膛压和速度即可确定火药燃烧分裂点。时，求出对应的时间的膛压和速度即可确定火药燃烧结束点。时，求出对应的时间的膛压和速度即可确定火炮口点。方程组通过求得的初值可化简为 （公式 4.

29、4.1 化简后的方程组）其中初值利用matlab初值计算程序算出，通过初始值计算出t=0ms、 2 ms、3.5ms、 5.3ms、 7ms、 8.7ms、 10.5ms、 12.7ms、 14ms的压强P、速度v、以及弹丸行程的值的大小,利用matlab将这些点拟合得到内弹道的Pt P vt v曲线通过已知时间得到对应的膛压，速度和行程的值，通过matlab绘制出Pt P vt v曲线【4】计算得到对应点的膛压，行程，速度值如下t=0, 2 ,3.5, 5.3, 7 ,8.7, 10.5, 12.7, 14 ; P=300, 900, 1700, 2780, 3281, 3170, 2361

30、, 1593 ,1251;L=0 ,0.28, 1.2, 4.7 9, 17 ,28, 45, 56;v=0 ,34.4 ,95, 225.46, 380.53, 556.5, 702.5, 825, 902;结 论1983式152mm加农炮的膛压随时间、行程，速度随时间、行程的变化曲线可以通过matlab绘制出来将所计算的点的值通过matlab绘制得到的所要求的曲线如下图图4.1 膛压随时间变化趋势图 4.2 弹丸速度随时间变化趋势图 4.3 膛压随弹丸的运动行程变化曲线图 4.4 速度随弹丸行程变化曲线图 4.5 初值的输入程序图 4.6 程序的执行附 录1 绘图程序如下t=0 2 3.5

31、 5.3 7 8.7 10.5 12.7 14 ; P=300 900 1700 2780 3281 3170 2361 1593 1251;L=0 0.28 1.2 4.7 9 17 28 45 56;v=0 34.4 95 225.46 380.53 556.5 702.5 825 902;%tangya-shijia figure(1) plot(t,P),gridtitle(P-t)xlabel(t/ms )ylabel(P/(100KPa) （膛压时间曲线绘制程序）pause(6) % sudu-shijian figure(2);plot(t,v),gridtitle(v-t)xl

32、abel(t/ms)ylabel(v/(m/s) （速度时间曲线绘制程序）pause(6) % tangya-xingchengfigure(3)plot(L,P),gridtitle(P-L)xlabel(L/dm )ylabel(P/(100KPa) （膛压行程曲线绘制程序）pause(6)% sudu-xingchengfigure(4)plot(L,v),gridtitle(v-L)xlabel(L/dm )ylabel(v/(m/s) （速度行程曲线绘制程序）初值赋予和计算程序如下S=input(please input the S S=);V0=input(please input

33、 the V0 V0=);lg=input(please input the lg lg=);m=input(please input the m m=);f=input(please input the f f=);w=input(please input the w w=);a=input(please input the a a=);Pp=input(please input the Pp Pp=);cta=input(please input the cta cta=);u1=input(please input the u1 u1=);n=input(please input the

34、 n n=);e1=input(please input the e1 e1=);X=input(please input the X X=);LAM=input(please input the LAM LAM=);u=input(please input the u u=);LAMs=input(please input the LAMs LAMs=);P0=input(please input the P0 P0=);fai=input(please input the fai fai=);Xs=input(please input the Xs Xs=);d=input(please

35、input the d d=); （输入火药诸元和初值程序） %初值计算disp(S)disp(V0)disp(lg)disp(m)disp(f)disp(w)disp(a)disp(Pp)disp(cta)disp(u1)disp(n)disp(e1)disp(X)disp(LAM)disp(u)disp(LAMs)disp(P0)disp(fai)det=w/V0,V01=0,t01=0,P01=P0/(f*det),L0=V0/Svj=sqrt(2*f* w/(cta*fai*m)B=S2*e12*(f*det)(2-2*n)/(f* w* fai*m*u12)rou=0.5*(d+2e

36、1)*0.5591Zk=(e1+rou)/e1ksai0=(1/det-1/Pp)/(f/P0+(a-1/Pp)Z0=0.5*sqrt(1+4*LAM*ksai0/X)-1)/LAM （初值计算并输出程序）结束语 首先感谢我的指导老师王健老师!经过三周的忙碌和工作，本次课程设计已经接近尾声，通过这次设计，把自己所学的理论知识运用于实践，受益非浅。对弹道学特别是内弹道有了更深的认识,对我以后独自处理问题有了很大的帮助。在整个课程设计过程中，我的指导老师王健老师一直给予我们很多的关注和指导，勤勉教学，诲人不倦，帮助我们顺利的完成这次设计任务。通过这次课程设计让我明白理论和实践还是有很大的距离，增强了自己查阅资料解决问题的能力。参考文献【1】 金志明.枪炮内弹道学.第一版.北京理工大学出版社，2004年.【2】 王连荣.火炮内弹道计算手册.第一版.国防工业出版社.1987年.【3】 钱林方.火炮弹道学.第一版.北京理工大学出版社.2009年.【4】 周开利.matlab基础及应用教程.第一版.北京大学出版社.2009年.