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