[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;