飞航导弹系统工程与均匀设计法

航天工业总公司三院三部 张建舟

返回

  摘要:本文通过某导弹武器系统的导弹导弹火控系统的数学模型的建模和验证,说明均匀设计法是一种非常适合於在大型复杂系统工程建立数学模型的过程中,是实时求取一些重要关系表达式的有效方法。与逐步回归分析一起,可以给出符合高精度要求的结果。

一、飞航导弹系统工程的特点

(一)无人驾驶有控飞行器

  所有的飞航导弹和无人驾驶飞行器,虽然在发射制导过程中,还要由人来装定一些参数,但在人不控制的飞行状态,基本上是一个自动控制飞行的飞行器。

(二)独立完整的系统──整体性好

  飞航式导弹虽然有各种各样的形式,各种各样的部件,但大体上都由弹体结构系统、控制系统、动力系统、导航系统、引信战斗部系统、火控系统等构成。这里还不提技术阵地等更大的系统构成部分。各系统虽有不同的性能,但它们是根据一定的原理、合乎飞航导弹逻辑统一性要术构成的整体。为了照顾和保证系统的总体性能好(整体性),必须按照一定的顺序充分注意系统内各个层次和各个组成部分的协调与联结才能成功。

(三)密切联结──相关性好

  飞航导弹各个系统虽然相对独立,但却是相互联结、相互制约、相互依赖完成同一目的地统一体。从系统的角度看导弹系统是密切联结──相关性好。不仅在定性地评价导弹总体性能时,就是定量地研究总体性能时,都可以看到体现系统性能的参数与各影响因素之间的密切关系。

(四)目的明确──唯一性好

  以某飞航式导弹为例,由於导弹为进攻性武器,在其基础上改造成功一个被进攻的靶弹,看上去绝大部分部件相同或相似。但仔细比较不难看到那些构造上的微小区别和性能上的微小差异。也就是这些差异(别)成了它们是不同目的不同功能的飞行器。

(五)环境的作用和影响──重要因素的显著性好

  飞航导弹在离开载体刚性联结,若载体是动基座,则载体对飞航导弹的约束效果体现在初始条件上。另外,飞航导弹在空气中飞行,空气对飞航导弹的作用如阻力、升力、操纵力均与空气的状态如密度、压力、温度有关。扇面装定角引导弹拐弯飞行,环境与导弹系统不断进行着物质的、能量的、信息的交换。所有的作用约束,也就是因素的影响,恰好就是概率统计中用到的,本文重点讨论的问题。

二、两种数学模型

  以某种飞航导弹为例。如图1所示

R──导弹助推器的推力

P──导弹主发动机的推力

G──地球对导弹的重力            g──引力加速度

Q──阻力                 ──阻力系数

──升力                ──升力系数

──升力系数导数            ──侧力

──侧力系数               ──侧力系数导数

α──攻角                 β──侧滑角

θ──弹道偏角               ──助推器与弹轴夹角

m──导弹的质量               ──导弹的滚转角

──导弹的速度               ──导弹的偏转角

──相对於弹体的旋转角速度

──导弹的转动惯量

──力距系数与力距系数的偏导数

──力距系数与力距系数导数

──马赫数               xyz──导弹在大地坐标系中坐标

──俯仰通道控制信号

──分别为俯仰、航向、滚转通道控制信号总和

──分另为升降舵、方向舵、滚转舵的舵偏角

(一)导弹运动参数的精确值

  导弹运动弹道(地面坐标系中的xyz值)与恣态值、各控制参数可通过发射导弹时飞行试验的光测和遥测得到。但为控制导弹准确击中目标而必须装定的控制参数。却需要许多光测弹道──发射许多枚导弹才能得到。由於耗资巨大加之许多实际问题难以解决,这是目前还解决不了的课题。

(二)导弹运动参数的基准值

  随着对导弹运动的深入研究和计算技术的飞速发展,在大量验证和不断修正完善的基础上,产生了精确描述导弹运动弹道(地面坐标系中的xyz)飞行性能。控制参数的空间质点系六自由度弹道控制联立方程组。包括方程在内的这些方程共有多个。用这一联立方程组在导弹设计、运动状态分析、靶场试验设计和分析已足够精确。

  为了既说明问题而又醒目,下面仅列出方程组中个别方程来讨论。

  1. 导弹弹体控制飞行的数学模型

 

    

    

    

    

    

    

 

 

 

 

 

 

      

 

      

  

  

  

 

  

 

 

      

 

 

 

      

控制方程

      

 

 

 

      

  导弹的运动方程和控制方程共有120多个。随着研究内容的不同和研究者的习惯不同表达方式有些差异。

  联立求解这些方程,就可以得到包括导弹坐标(xyz)在内的参数。在这些方程中非常清楚的反映了导弹系统和环境系统的物质的、能量的、信息的交换关系。

  例如,质量m是变质量,导弹质量的减少就是和环境进行质量交换的过程; 升力Y、阻力Q、推力RP、等与环境大气温度、空气密度、空气压力密切相关。也就是说与飞行区间的大气状态密切相关,且不断进行着能量的交换;导弹在飞行过程中通过无线电高度表或气压高度表不断测量飞行高度和垂直速度,用这些信号来控制导弹飞行,也就是不断进行着信息的交换。

2导弹在动基座上运动数学模型

  导弹在动基座上的运动,随着发射方式的不同而有区别。以某飞航式导弹为例,在斜式装置上的,导弹随着舰艇一起在海上摇摆,海浪对导弹的运动影响,通过力和加速度的形式体现在方程里。舰艇航速、航向的影响也在描述导弹运动速度、偏航角的方程里有所体现。描述导弹在舰艇上斜式发射的方程共有近50个运动方程和30多个控制方程。

 

  为了说明问题这里只列必要的方程.

(1)前滑块离轨前导弹的质心运动方程

        

        

        

        

      这里的是海平面坐标系中的角度。

(2)前滑块离轨到后滑块离轨前运动方程

        

        

        

      

        

        

  舰艇对导弹的作用,以牵连速度、作用力N、力距、惯性力等的形式影响着导弹的恣态角和空间坐标(xyz),而装定的扇面角在前滑块离轨前由於导轨的约束,在前滑块离轨后由於导轨和空气动力小的原因。在大推力作用下开始偏转但其效应需过一段才会明显看到。

(3)控制方程

  只有非常清楚地了解导弹和周围环境的相互作用,才会更好地解决构造的数学模型的课题。我们从图2可以看出导弹和 怎样进行物质的、能量的、信息的交换的。

三、用系统工程的观点研究系统

  研究系统和系统工程,要求其最优,就必须掌握系统思想。这就是:

(一)总体观点

  由於系统有大有小,研究者首先要搞明白自己研究的是什么内容?特别要弄明白哪些是外部因素?哪些是内部系统?哪些是输入和输入参数?哪些是输出和输出参数?只有总揽全局,才能看清方向。

(二)层次观点

  研究系统工程,必须深入研究系统的构造,以便搞清楚各个层次的关系,对那些隐含在研究范围里的关系,内部的内部、外部的外部的关系式。可以深入了解,多掌握些材料,但最主要的是搞深搞透描述对象。多余的关系式不仅增加了工作量,而且有可能造成捡了芝麻丢了西瓜的结果。

  在研究火控系统的模型时,对於导弹内部的一些系统,如引信、战斗部、弹体结构、液压系统、气压系统......可不涉及,这绝不是说与其没有关系。

(三)反馈观点

  导弹在飞行中对空气进行冲击挤压,空气(含风速、风向)以升力、阻力、侧力及其力距的形式与其它因素一起对导弹反作用(反馈)。

(四)相关性观点

  既然研制的飞航导弹能相互有机地联结成一个统一的系统,构成系统的各要素之间必然互相作用、相互依存、相互制约。各系统有条件地形成一个紧密地有飞航导弹特征的整体,用逐步回归分析方法可以定量地给出其参数与影响因素的相关系数。在回归方程中只有通过F检验相关系数大的影响因素才能保留下来。

(五)适应性观点

  所有的系统都存在一定的环境之中,当环境在规定的范围变化时,只有能灵敏地适应环境变化的系统才有可能存在。只有性能先进的飞航导弹才能在竞争中有市场。

(六)最优化的观点

  所有的飞航导弹系统之所以在一定时期,一定环境条件下研制成功,总有一定的道理。追求研制的系统最优,是所有参与者、设计者、决策者共同的目标。是否最优?能否存在下去?这需要接受检验。

  同所有真理具有相对性一样,最优化是追求的目标之一。设计、研制、生产、决策都受到许多条件的限制和制约。随着科学技术的发展,人类认识的深化,包括人类和自然界在内的系统都在进化,演变。何况人研制的系统。最好的方法、最优的试验、最优的系统都是相对的、有条件的,在一定时期外部内部环境条件下成立的。一个飞航导弹系统若不及时改进就难永远保持是先进系统。

四、火控系统数学模型

(一)温度对导弹射程的影响

  由於计算机技术的飞速发展;由於科学技术的发展和作战的需要;老式短射程的飞航导弹系统及其机电式指挥仪已不能完成新的任务,在追求较远飞航导弹系统的研制过程中,我院梁守磐、王大卫、周方等科研人员发现了温度对导弹射程的影响,并引入了修正公式。这一过程历时几年。

(二)扇面发射角的引入

  1976年,随着一项新任务的立项,张建舟提出了扇面发射角和兼有扇面发射功能的武器系统设想,滕衣虹、宣平等提出采用无线电高度表降高飞行采用红外导引头的方案。在梁守磐、吴宝初、曹柏桢、宣平、潘荣霖、杨新生的领导下,张建舟、潘乃强、申宗鹤等经过一年的计算分析找到了的精确描述导弹自控段坐标(xz)的表达式,并用正交设计布点计算弹道,采用逐步回归分析的方法给出了火控系统数学模型。柯凡提出舰速影响应预考虑,搞另一任务的赵利华、马恒华和张炳辉、张建舟等先后走访了中国科学院数学所方开泰先生,方开泰、王元深入研究计算创造性地提出了均匀设计法,在方开泰先生的指点下,搞同一任务的张建舟、谷巨卿、张炳辉与柯凡分别采用均匀设计法布点计算弹道,与张炳辉、黄树山一起彩逐步回归分析方法给出了回归关系式的各项系数和火控系统的数学模型。由於种种原因,本文无法提及赵利华等人的杰出贡献和结果。

  本文提及这些 历史,主要是想说明:系统工程是许多人参加的复杂的系统,许多好主意无论由哪一级哪个人提出,都需要决策者批准并经许多技术人员和工人的创造性劳动才能实现,这里就不一一列出他们的贡献了。只有完整地、深入地至少经过一个系统工程的周期,才能体会到钱学森的论断:系统工程是组织系统的规划、研究、设计、制造、试验、和使用的方法,是一种对所有系统都具有普遍意义的科学方法。

  扇面角的引入,彻底的改变了初始的设想,控制体制的彻底改变,由一开始的稳定飞行改为导引飞行产生了质的变化,因而在原有的近射程AB常数或变系数以速度为模式得不到所相似的模型。所以要另找新的模型。在长达两年的分析计算过程中。搞明白了原因,也找到了办法。

(三)用均匀设计法布点计算基准弹道

    由前面的分析已经清楚,找寻新模型已具备了条件。

  1. 为求解的关键

    由前置射击原理图知道,摇摆修正瞄准角,自控段终点时间在下列方程中。

  求解的真值是求的关键

(1)瞬时求解100多个联立方程目前不现实

  火控系统需要瞬时求解并连续装定一些控制参数,20多年前算一条弹道要20多分钟,就是目前也难以作到联立求解100多个方程可以瞬时给出当时值。不是瞬时值会贻误战机。

  正交设计、均匀设计与逐步回归分析方法,为工程技术人员提供了可以用多项式拟合的办法来解决这个课题。

(2)均匀设计的优越性

  虽然一般工程技术人员不是搞和统计专业的,但是我们也知道,工程上的许多问题,都是要求在使用的范围内,使用的机会是均匀的、随机的。也就是说要求充分的均匀。在计算机日益普及和性能大大提高的今天,由於设备造价太高,试验经费的昂贵,有优良性能而又能大大节省经费的均匀设计,当然为科技人员所优先选用。

  如前所述,我们可以选用以100多个联立方程按均匀设计原则布点计算的弹道值作为导弹弹道的真值,依据这些值来求取xz......的回归值.

  文献[6]已证明了

下面简述一下如何排表求取。

(3)欲求的变量与其影响因素的关系已明确

  通过上述分析,影响所求xz的因素的关系已明确;

风在坐标系两个轴上的分量;

舰速在坐标系两个轴上的分量:

大气温度

扇面装定角

射程:x

发射装置瞄准角:

2.选用均匀设计表 选定计算条件

  根据文献1所提供的表,按特性分3个区域安排计算情况。给定了按公式可算出温度T()按使用区域等分成22个水平后换算成绝对温度,按区域特性分别取值 ,注意定义规定。如下表所示

W

T()

Vj

Wfx

Wfz

Vjx

Vjz

-4.0

6.6668

196.363

295.618

14.3072

55

-6.3972

-1.87659

8.2071

-11.7191

-2.619

14.286

32.727

259.19

8.942

..

12.0184

7.72307

5.12943

-7.32449

0.143

7.6192

81.818

279.428

18.7782

..

1.08514

7.54152

10.7718

-15.3814

1.524

15.2834

294.544

243.0

13.413

..

6.3248

-13.8638

7.69415

-10.9867

2.905

0.9524

130.908

299.666

8.0478

..

-0.62356

0.71988

4.61649

-6.59204

4.286

8.5716

243.635

263.238

2.6826

..

-3.80886

-7.67885

1.53883

-2.19734

5.667

16.1908

179.999

319.904

17.884

..

-16.1907

0.00385

10.2588

-14.6489

7.048

1.9048

16.3636

283.476

12.5188

..

1.82765

0.5366

7.1812

-10.2542

8.429

9.524

229.09

247.047

7.1536

..

-6.23907

-7.19586

4.10355

-5.85959

9.81

17.1432

65.454

303.714

1.7884

..

7.12289

15.5933

1.02588

-1.46489

11.191

2.8572

278.181

267.285

16.9898

..

0.40558

-2.8286

9.74593

-13.9165

12.572

10.4764

114.545

323.952

11.6246

..

-4.3506

9.53032

6.66826

-9.52184

 

W

T()

Vj

Wfx

Wfz

Vjx

Vjz

12.953

18.0956

327.272

287.523

6.2594

..

15.2187

-9.78981

3.5906

-5.12714

15.334

3.8096

163.636

251.095

0.8942

..

-3.655

1.07408

0.51294

-0.73244

16.715

11.4288

0

307.76

16.0956

..

11.4288

0

9.23298

-13.184

18.096

19.048

212.726

271.323

10.7304

..

-16.027

-10.2936

6.15532

-8.78939

19.477

4.762

49.091

327.999

5.3652

..

1.74042

4.43255

3.07766

-4.39469

20.858

12.3812

261.817

291.571

0

..

2.71242

-12.0804

0

0

22.237

20.0004

98.182

255.142

15.2014

..

-10.1244

17.2485

8.72004

-12.4516

23.62

5.7144

310.908

311.809

9.8362

..

5.15789

-2.45977

5.64238

-8.05694

25.001

13.3336

147.272

275.380

4.471

..

-13.2121

1.79557

2.56471

-3.66224

 

W

T()

Vj

Wfx

Wfz

Vjx

Vjz

0

6.667

196.363

295.618

14.3072

55.00

-6.39723

-1.8766

8.2071

-11.7191

0

14.286

32.727

259.19

8.942

58.095

12.0184

7.72307

4.72652

-7.5907

0

0

245.454

315.856

3.57680

61.19

0

0

1.7239

-3.1339

0

7.6192

81.818

279.428

18.7782

64.285

1.0851

7.5415

8.1492

-16.9177

0

15.2384

294.544

243

13.413

67.38

6.3249

-13.8638

5.1599

-12.3807

0

0.9524

130.951

299.666

8.0478

70.475

-0.6236

0.7199

2.6904

-7.5848

0

8.5716

243.635

263.238

2.6826

73.57

-3.8089

-7.6789

0.7590

-2.5730

0

16.1908

179.999

319.904

17.884

76.665

-16.1907

0.00385

4.1266

-17.4014

0

1.9048

16.364

283.476

12.5188

79.76

1.8277

0.5366

2.2268

-12.3191

0

9.524

229.09

247.047

7.1536

82.855

-6.23907

-7.1959

0.8905

-7.0979

0

17.1432

65.454

303.714

1.7884

85.95

7.12289

15.5933

0.1265

-1.7839

0

2.8572

278.181

267.285

16.9898

89.045

0.40558

-2.8282

0.2851

-16.9874

0

10.4764

114.545

323.952

11.6246

92.14

-4.3506

9.5303

-0.4327

-11.6165

0

18.0956

327.272

287.523

6.2594

95.235

15.2187

-9.7898

-0.5704

-6.2334

0

3.8096

163.636

251.095

0.8942

98.33

-3.65505

1.0741

-0.1294

-0.8848

0

11.4288

0

307.761

16.0956

101.425

11.4288

0

-3.1862

-15.7770

0

19.084

212.726

271.333

10.7304

104.52

-16.027

10.2936

-2.6889

-10.388

0

4.762

49.091

327.999

5.3652

107.615

3.1187

3.5987

-1.6229

-5.1138

0

12.3812

261.817

291.571

0

110.71

-1.7663

-12.2545

0

0

0

20.0004

96.182

255.142

15.2014

113.805

-2.8438

19.7971

-6.1336

-13.909

0

5.7144

310.908

311.809

9.8362

116.9

3.74035

-4.3202

-4.4489

-8.7729

0

13.3336

147.272

275.38

4.471

119.995

-11.2155

7.2109

-2.2345

-3.8725

 

W

T()

Vj

Wfx

Wfz

Vjx

Vjz

4.000

6.667

196.363

295.618

14.3072

120

-6.2508

-2.3182

-7.1517

-12.3915

2.619

14.286

32.727

259.190

8.9420

...

11.653

8.2641

-4.4698

-7.7447

1.238

0

245.454

315.856

3.5768

...

0

0

-1.7879

-3.0979

-0.143

7.619

81.818

279.428

18.7782

...

1.104

7.5388

-9.3866

-16.2638

-1.524

15.238

294.544

243.00

13.413

...

5.9539

-14.0271

-6.7047

-11.617

-2.905

0.952

130.908

299.666

8.0478

...

-0.5863

0.7506

-4.0228

-6.9702

-4.286

8.572

243.635

263.238

2.6826

...

-4.372

-7.3728

-1.3409

-2.3234

-5.667

16.191

179.999

319.904

17.884

...

-16.1112

1.6025

-8.9336

-15.4893

-7.048

1.905

16.364

283.476

12.519

...

1.8797

0.3083

-6.2577

-10.8425

-8.429

9.524

229.09

247.047

7.1536

...

-7.2264

-6.2034

-3.5758

-6.1958

-9.91

17.143

65.454

303.714

1.7884

...

9.6754

-14.1518

-0.8940

-1.5489

-11.191

2.857

278.181

267.285

16.9898

...

-0.151

-2.853

-8.4926

-14.7129

-12.572

10.476

114.545

323.952

11.6246

...

-2.172

10.2487

-5.8107

-10.068

-13.953

18.096

327.272

287.523

6.2594

...

12.409

-13.170

-3.1286

-5.4213

-15.334

3.810

163.636

251.095

0.8942

...

-3.2409

2.0023

-0.4470

-0.7745

-16.715

11.429

0

307.761

16.0956

...

10.9459

-3.2868

-8.0456

-13.9404

-18.096

19.048

212.726

271.333

10.7304

...

-18.4315

-4.8067

-5.3638

-9.2936

-19.477

4.762

49.091

327.999

5.3652

...

4.1401

2.3530

-2.6818

-4.6468

-20.858

12.381

261.817

291.571

0

...

-6.0135

-10.8227

0

0

-22.239

20.000

98.182

255.142

15.2014

...

4.8599

19.4009

-7.5987

-13.1659

-23.620

5.714

310.908

311.809

9.8362

...

1.6962

-5.4569

-4.9169

-8.5195

-25.001

13.334

147.272

275.38

4.471

...

-7.1174

11.275

-2.2349

-3.8723

 

  考虑到实际出现的极少概率情况,故没有采用表中最下边的第23行,使回归结果精度既高又合乎实际情况。

3.弹道参数拟合公式的建立

  弹道参数的计算结果,反映了导弹在空中时在装定情况下的每一瞬时的实际结果, 除去非稳定情况下的值 ,考虑到n次幂时的极性,......采用逐步回归分析中的基本公式便可得到下列表达式。

    

关系式

  以上仅是火控系统数学模型中的几个公式,其它不一一列举。构造的火控模型在打靶中非常成功。经受了长期的检验。

五、几点体会

(一)重示定义域内有效数据区域

  在采集数据时,一定要使数据充分的代表经常使用的情况,数据会改变回归公式,尤其在变加速度的非稳定状态的数据,会破坏回归结果,至少使其精度降低。在采集数据时不是预期的使用区域的数据最好不用。从这一点来看原联立方程组基本上反映了导弹在任一时刻的真实姿态和弹道,而通过均匀设计法和逐步回归分析所得到的某些参数的回归值 只在其定义域内是真实地在数值上反映了其变化量。要较为全面地反映在整个定义域内的值只能采用分段回归的办法。对某些拐点及其剧烈变化处应采取特殊办法。这首先得看必要性和工作量及是否值得那样作。

(二)实践成功与否是试金石

  回归公式是一些影响因素及其交互作用项的多项式代数和,有时在一些气动系数的表达式中出现了马赫数(M)六次方和五次方的项,按一些物理量的定义似难以理解,不必死钻牛角尖,更不必人为的拿掉它。工程实际是唯一的试金石,若回归结果经过试验检验是正确的,说明从另一个角度看,必有深刻的意义等待去认识它。

(三)加权回归看实况

  加权回归,应在出现较大的区域采用.若在一个充分需要均匀的区域里,人为的在某一区域的精度。

(四)不同模型的根源

  同一类型的工程,同一个参数的表达式可能不止一个,我们经历过同一参数表达式不同的两种类型工程都取得圆满成功的系统工程的实践,差别主要是看问题的角度不同而已。出现了同一参数不同模型的表达式别轻易否定它。要分析原因,靠实验验证。

(五)外推要慎重

  回归公式有一定的使用区域,凡外推使用必须慎重。

六、结束语

  本文标题虽然是飞航导弹系统工程与均匀设计,但仅仅是某飞航导弹系统工程中的一例经历,自均匀设计法诞生以来,均匀设计在不同的课题中均取得了成功,更在许多行业中广为流传。本文只是说明均匀设计法在系统工程中使用很有效,很成功。

  关世义先生仔细审阅了本文,对技术部分提了许多宝贵意见,在此表示感谢。

 

 

参考文献

[1] 钱学森等著,论系统工程,湖南技术出版社,13204-67.

[2]