[9] 吴兴波.对武器装备可靠性指标体系的探讨.火炮发射与控制学报,1997(4):17.19.

[10] 侯保林.大口径火炮弹药自动装填系统研究[D].南京理工大学,2003.

[11] 于惠利.传动零部件设计实例精解[M].北京:机械工业出版社,2009.

[12] 孟祥宾等.德国链传动标准译文集[M].北京:中国标准出版社,2003.

[13] 杨义勇,金德闻.机械系统动力学.清华大学出版社,2009.

[14] 刘寅立,王亮剑等.MATLAB数值计算案例分析.北京航空航天出版社,2011.

[15] 张岳.MATLAB程序设计与应用基础教程.清华大学出版社,2011.

[16] 李庆扬,王能超等.数值分析第5版.清华大学出版社,2010.

[17] 黎明安. 动力学系统建模与仿真.北京国防工业出版社,2012.

[18] 常备功,樊智敏等.带传动和链传动设计手册.北京化学工业出版社,2010.

 

 

 

 

 

 

 

 

 

 

 

 

附录:

弹筒1运动微分方程:

%弹筒1的运动方程

function dy=dantong1(t,y)

dy=zeros(2,1);

%求各时刻,各转角下的等效转动惯量及其对转角的导数

%求各时刻,个转角下的等效力矩(所有弹筒的重力产生的等效力矩,所有弹筒摩擦力产生的等效力矩)

%RJ-不同时刻各弹筒的等效转动惯量

%EJ-不同时刻所有弹筒的等效转动惯量

%F-不同时刻各(每个)弹筒重力产生的等效力矩

%EF-不同时刻所有弹筒的重力产生的等效力矩

%FM-不同时刻各(每个)弹筒摩擦力产生的等效力矩

%EFMDT-不同时刻所有弹筒摩擦力产生的等效力矩

%DJ1DF1-各转角下的等效转动惯量对转角的导数

OMGA=100.0/180.0*pi;

DT=0.5e-3;

N=pi/2.0/OMGA/DT;

for IT=1:(N+1);

    T=(IT-1)*DT;

    XTA(IT)=T*OMGA;

    RM1=68.8398; %RM1--弹筒1的质量,由计算机测量获得

RJ1=0.294227; %RJ1--弹筒1的转动惯量,由计算机测量获得

U=0.15; %U--摩擦系数

OB=131e-3; %参数(1)

BA=sqrt(2*131e-3^2); %参数(2)

OB1=131e-3; %参数(3)

CA=sqrt(72e-3^2+92.63e-3^2); %参数(4)

CAB=atan(72e-3/92.63e-3); %参数(5)

R=pi/180.0;

BB1=sqrt(OB^2+OB^2-2.0*OB^2*cos(XTA)); %(1-1)

BB1O=(pi-XTA)/2.0; %(1-2)

BB1A=pi/2+BB1O; %(1-3)

B1AB=asin(BB1./BA.*sin(BB1A)); %(1-4)

ABB1=pi-BB1A-B1AB; %(1-5)

B1A=sqrt(BB1.^2+BA.^2-2.*BB1.*BA.*cos(ABB1)); %(1-6)

OA=sqrt(B1A.^2+OB1.^2); %(1-7)

S1=BA-B1A; %(1-8)

OBA=acos((BA.^2+OB.^2-OA.^2)./(2.0.*BA.*OB)); %(1-9)

ALF=OBA; %(1-10)

VB=OB.*OMGA; %(1-11)

VZ=VB.*sin(XTA)./sin(pi-XTA-ALF); %(1-12)

OMGAZ=VZ./BA; %(1-13)

VP1=VB.*sin(ALF)./sin(pi-XTA-ALF); %(1-14)

CAB1=CAB-B1AB; %(1-15)

BTA=CAB1; %(1-16)

GMA=pi-(pi./2.0-BTA); %(1-17)

VC=sqrt(VP1.^2+(CA.*OMGAZ).^2-2.0.*VP1.*CA.*OMGAZ.*cos(GMA)); %(1-18)

ALFG=acos(CA.*OMGAZ.*cos(BTA)./VC); %(1-19)

RJ1=RM1.*(VC./OMGA).^2+RJ1.*(OMGAZ./OMGA).^2;

F1=RM1.*9.8.*VC.*cos(ALFG)./OMGA;

FM1=RM1.*9.8/2.0.*U.*VP1/OMGA.*cos(pi);

end

DJ1DF1(1)=(RJ1(2)-RJ1(1))./(XTA(2)-XTA(1));

DJ1DF1(N+1)=(RJ1(N+1)-RJ1(N))./(XTA(N+1)-XTA(N));

DJ1DF1(N/2+1)=0;

DJ1DF1(N/2+2)=0;

for i=2:(N/2);

DJ1DF1(i)=(RJ1(i+1)-RJ1(i))./(XTA(i+1)-XTA(i))./2.0+(RJ1(i)-RJ1(i-1))./(XTA(i)-XTA(i-1))./2.0;

上一篇:反末敏弹战斗部设计+Autodyn仿真
下一篇:某拦截系统拦截概率分析+文献综述

密集光波分复用系统的波...

变节距排缆拖曳绞车电液控制设计+CAD图纸

SolidWorks变节距排缆拖曳绞...

FDS船舶受限空间内高温热...

RAE2822翼型大长径比火箭弹气动特性数值模拟

DEFORM-3D大高径比长杆件镦粗防失稳模具设计

AT89S51双波长的差动技术新...

AT89C52单片机的超声波测距...

神经外科重症监护病房患...

医院财务风险因素分析及管理措施【2367字】

国内外图像分割技术研究现状

承德市事业单位档案管理...

C#学校科研管理系统的设计

中国学术生态细节考察《...

志愿者活动的调查问卷表

公寓空调设计任务书

10万元能开儿童乐园吗,我...