什么是d道导d

d道导d是指在火箭发动机的推力下，按照预定的程序飞行，关机后再按照自由抛掷物体的轨迹飞行的导d。其飞行轨迹一般分为主动段和被动段:主动段(也称动力飞行段或助推段)是导d在火箭发动机推力和制导系统的作用下，从发射点到火箭发动机关机时的飞行轨迹；被动段包括自由飞行段和再入段，是导d按照主动段末获得的给定速度和d道角度惯性飞行，然后到d头起爆的路径。[1]

2022年4月9日发表声明称，巴基斯坦当天成功试射一枚“沙欣-3”地对地d道导d。声明没有透露试射导d的d着点和具体发射地点。“沙欣-3”d道导d具备携带核d头和常规d头的能力，射程2750公里。[9]

导d包括了d头和推进装置，也就是火箭

火箭在燃烧完就会自行脱落，就像发射卫星一样。比如说三级推进的洲际导d，第一级火箭大概工作2分钟就会脱落，主要是给导d加速。二级火箭主要维持导d的速度，其工作主要是在外层空间，此时的速度会达到第一宇宙速度；而三级火箭为导d的机动提供动力，是火箭重返大气层，准确击中目标。重返大气层时摩擦会产生大量的热，d头会达到几千度高温，当然，第三级火箭早就脱落焚毁了。d头在到达目的地后爆炸，如果打地下目标，会利用d头高速运动具有的动能延迟爆炸；如果摧毁地面目标的核d，则会在离地面200米爆炸，使冲击波，核辐射 有效作用在地面，同时增大毁伤面积

原因

所给链接的代码不完整，缺少ndd_fun函数。

顺便鄙视一下该链接的提供者，这么广为流传的东西下载居然还要积分，简直穷疯了。

代码

帮你好好找了一下，找到了完整的程序，供参考（全部代码保存到一个M文件运行即可，或直接下载附件）：

%59nian130

A=087;         %q(炮)膛横断面积A  dm^2

G=19;%334;           %d重  kg

W0=204;        %药室容积  dm^3

l_g=250;       %身管行程  dm

P_0 =30000;       %起动压力  kpa

fai1=102;        %次要功系数

K=103;        %运动阻力系数φ1

theta =02;       %火药热力系数

%=========================================

f=950000;          %火药力  kgdm/kg 

alpha=1;         %余容  dm^3/kg

delta=16;         %火药重度γ

%==================================

ome=22;%129;      %第一种装药量  kg

u1=5002410^-5;        %第一种装药烧速系数  dm^3/(skg)

n1=082;         %第一种装药的压力指数n1

lambda=-00071;     %第一种装药形状特征量λ1

lambda_s=0;     %第一种装药分裂点形状特征量λ1s

chi=100716;          %第一种装药形状特征量χ1

chi_s=0;          %第一种装药分裂点形状特征量χ1s

mu=0;            %第一种装药形状特征量μ1

et1=11410^-2;            %第一种装药药厚δ01

d1=2510^-2;             %第一种装药火药内径d1

Ro1=0;             %药型系数α1

%=========================================

%常数与初值计算-----------------------------------------------------

l_0=W0/A;

Delta=ome/W0;

phi=K + ome/(3G);

v_j=196fome/(phithetaG);

v_j=sqrt(v_j);

B = 98(et1A)^2/( u1u1fomephiG );

B=B(fDelta)^(2-2n1);

Z_s=1+Ro1(d1/2+et1)/et1;

p_0=P_0/(fDelta);

psi_0=(1/Delta - 1/delta)/(f/P_0 + alpha - 1/delta);

Z_0=(sqrt(1+4psi_0lambda/chi) - 1)/(2lambda);

%解算子------------------------------------------------------------

C = zeros(1,12);

C(1)=chi;C(2)=lambda;C(3)=lambda_s;C(4)=chi_s;C(5)=Z_s;%

C(6)=theta;C(7)=B;C(8)=n1;C(9)=Delta;C(10)=delta;C(11)=alpha;C(12)=mu;

C;

y0=[Z_0;0;0;psi_0];

options = odeset('outputfcn','odeplot');

[tt,y] = ode45(@ndd_fun,0:100,[Z_0;0;0],options,C);

l = y(:,2);

l = ll_0;

fl = find(l>=l_g);

fl = min(fl);

[tt,y] = ode45(@ndd_fun,0:0005:fl,[Z_0;0;0],options,C);

Z = y(:,1);lx = y(:,2); vx = y(:,3); 

psi = (Z>=0&Z<1)( chiZ(1 + lambdaZ + muZ) ) +%%%%%%%%%

      (Z>=1&Z<Z_s)( chi_sZ(1 + lambda_sZ) ) +

      (Z>=Z_s)1;

l_psi = 1 - (Delta/delta)(1-psi) - alphaDeltapsi;

px = ( psi - vxvx )/( lx + l_psi );

p = pxfDelta/100;

v = vxv_j/10;

l = lxl_0;

t = ttl_01000/v_j;

fl = find(l>=l_g);

fl = min(fl)+1;

p(fl:end)=[];v(fl:end)=[];l(fl:end)=[];t(fl:end)=[];

pd=pxfDelta/100/(1+ome/3/fai1/G);

pt=pd(1+ome/2/fai1/G);

aa=max(px);

M=find(px==aa);

Pm=[tt(M)l_01000/v_j lx(M)l_0 vx(M)v_j/10 px(M)fDelta/100 pt(M) pd(M) psi(M) Z(M)];

%ll=length(tt);

ran=find(Z>=1);

ran=min(ran);

Zf=[tt(ran)l_01000/v_j lx(ran)l_0 vx(ran)v_j/10 px(ran)fDelta/100 pt(ran) pd(ran) psi(ran) Z(ran)];

jie=find(psi>=1);

jie=min(jie);

psij=[tt(jie)l_01000/v_j lx(jie)l_0 vx(jie)v_j/10 px(jie)fDelta/100 pt(jie) pd(jie) psi(jie) Z(jie)];

pg=[tt(end)l_01000/v_j lx(end)l_0 vx(end)v_j/10 px(end)fDelta/100 pt(end) pd(end) psi(end) Z(end)];

Ry1=[Zf;psij;pg;Pm];

Ry2=[ttl_01000/v_j lxl_0 vxv_j/10 pxfDelta/100 pt pd psi Z];

subplot(2,2,1);

plot(t,p,'linewidth',2);

grid on;

xlabel('\fontsize{8}\bft  (ms)');

ylabel('\fontsize{8}\bfp  (kg/cm^{2})');

title('\fontsize{8}\bft-p曲线');

subplot(2,2,2)

plot(t,v,'linewidth',2);

grid on;

xlabel('\fontsize{8}\bft  (ms)');

ylabel('\fontsize{8}\bfv  (m/s)');

title('\fontsize{8}\bft-v曲线');

subplot(2,2,3)

plot(l,p,'linewidth',2);

grid on;

xlabel('\fontsize{8}\bfl  (dm)');

ylabel('\fontsize{8}\bfp  (kg/cm^{2})');

title('\fontsize{8}\bfl-p曲线');

subplot(2,2,4)

plot(l,v,'linewidth',2);

grid on;

xlabel('\fontsize{8}\bfl  (dm)');

ylabel('\fontsize{8}\bfv  (m/s)');

title('\fontsize{8}\bfl-v曲线');

tspan = length(t)/20;

tspan = 1:ceil(tspan):length(t);

tspan(end) = length(t);

fprintf('        t(ms)     p(kg/cm^2)     v(m/s)       l(dm)');

format short g;

Result = [t(tspan) p(tspan) v(tspan) l(tspan)]

format;

%------------------------------------------------------------------

function dy = ndd_fun(t,y,C)

chi=C(1);lambda=C(2);lambda_s=C(3);chi_s=C(4);Z_s=C(5);mu=C(12);

theta=C(6);B=C(7);V=C(8);Delta=C(9);delta=C(10);alpha=C(11);

Z = y(1); l = y(2); v = y(3);

psi = (Z>=0&Z<1)( chiZ(1 + lambdaZ + muZ) ) +

(Z>=1&Z<Z_s)( chi_sZ(1 + lambda_sZ) ) +

(Z>=Z_s)1;

l_psi = 1 - (Delta/delta)(1-psi) - alphaDeltapsi;

p = ( psi - vv )/( l + l_psi );

dy(1) = sqrt(theta/(2B))(p^V)(Z>=0&Z<=Z_s);

dy(2) = v;

dy(3) = thetap/2;

dy = [dy(1);dy(2);dy(3)];

结果

火箭发射后的运动轨迹实际上就是火箭的d道。火箭的d道设计在火箭总体设计中是重要的一项内容，关系到火箭总体方案中的运载能力是否能够在d道约束下满足。对这点做下解释，运载火箭要将有效载荷送入轨道，除了要满足入轨的轨道要求外（可以理解为d道终点即为航天器的轨道起点），还要满足运载火箭的残骸落区要求，分离后的子级一方面要避开人口聚集区，也要避免落在其他国家领土领海，从而引起不必要的冲突，所以火箭的飞行轨迹（d道）是在火箭发射前就已经根据目标轨道的参数确定好了。在火箭发射前根据发射场气象条件（地面风速等）对d道进行修正后，修正后的d道被装订到火箭的飞行程序中。这就意味着，火箭的飞行轨迹理论上来说是应该按照装订好的d道飞行程序飞行的。火箭最完美的飞行轨迹应该是理论上计算的d道轨迹。如果不考虑各种偏差及随机干扰，所有的情况都和设计时预想的一样的话，火箭在点火起飞后，控制系统根据装订的飞行程序，按照时序控制发动机的伺服机构，从而控制发动机喷管的摆动，提供火箭进行程序转弯等动作的的控制力，实际上控制力是发动机推力的一个方向的分量。从而完成火箭的d道飞行。但这个过程实际上是理论上的，可以看出这是一个开环控制的过程。但是在实际情况中，火箭飞行时会遇到风切变、各种飞行误差的累计，从而会使得火箭的实际飞行d道偏离计算的d道轨迹。在这种情况下，如果还按照预定的动作执行控制，势必不能准确入轨。所以，在这里就需要根据飞行误差进行d道的修正。要进行误差的修正，需要先对误差进行识别和敏感，这就涉及到火箭的各种姿态、位置的敏感器件，实际上这是导航和制导系统需要完成的工作内容：对火箭的实际飞行位置和姿态进行感知、敏感和判断，火箭的控制系统根据实际飞行情况与预先装订的d道进行对比，根据计算出的偏差，发出控制指令，伺服机构控制喷管摆动，完成d道轨迹的修正。这样，整个火箭的控制变成了一个闭环过程

以上就是关于什么是d道导d全部的内容，包括:什么是d道导d、d道导d是指在火箭发动机推力作用下按预定程序飞行，那最后火箭去哪了、matlab 内d道计算问题等相关内容解答，如果想了解更多相关内容，可以关注我们，你们的支持是我们更新的动力！