BTA=BAO;
VA=AO.*OMGA;
VP3=VA.*sin(BTA)./sin(pi-BTA-(pi./2-XTA));
VZ=VA.*sin(pi./2-XTA)./sin(pi-BTA-(pi./2-XTA));
OMGAZ=VZ./AB;
B1BC=ABC-ABB1;
VC=sqrt(VP3.^2+(CB.*OMGAZ).^2-2.0.*VP3.*(CB.*OMGAZ).*cos(pi-(pi./2-B1BC)));
ALFG=acos(CB.*OMGAZ.*cos(B1BC)./VC);
RJ3=RM3.*(VC./OMGA).^2+RJ3.*(OMGAZ./OMGA).^2;
F3=RM3.*9.8.*VC.*cos(ALFG)./OMGA;
if XTA<=(pi/4.0)
FM3=0.0D0;
else
FM3=RM3*9.8/2.0*U*VP3/OMGA*cos(pi);
end
%弹筒4的运动分析,认为弹筒4在90°周期内做平动:
RM4=68.8398;
RJ4=0.294227;
U=0.15;
RJ4=RM4.*(VP3./OMGA).^2;
F4=0.0;
FM4=-RM4.*9.8.*U.*VP3./OMGA;
%弹筒5的运动分析,认为弹筒5在90°周期内作平动:
RM5=68.8398;
RJ5=0.294227;
U=0.15;
RJ5=RM5.*(VP3./OMGA).^2;
F5=0.0;
FM5=-RM5.*9.8.*U.*VP3./OMGA;
%弹筒6的运动分析,认为弹筒6在90°周期内做平动:
RM6=68.8398;
RJ6=0.294227;
U=0.15;
RJ6=RM6.*(VP3./OMGA).^2;
F6=0.0;
FM6=-RM6.*9.8.*U.*VP3./OMGA;
%弹筒7的运动分析,认为弹筒7在90°周期内做平动:
RM7=68.8398;
RJ7=0.294227;
U=0.15;
RJ7=RM7.*(VP3./OMGA).^2;
F7=0.0;
FM7=-RM7.*9.8.*U.*VP3./OMGA;
%弹筒8的运动分析,认为弹筒8在90°周期内做平动:
RM8=68.8398;
RJ8=0.294227;
U=0.15;
RJ8=RM8.*(VP3./OMGA).^2;
F8=0.0;
FM8=-RM8.*9.8.*U.*VP3./OMGA;
%弹筒9的运动分析,认为弹筒9在90°周期内做平动:
RM9=68.8398;
RJ9=0.294227;
U=0.15;
RJ9=RM9.*(VP3./OMGA).^2;
F9=0.0;
FM9=-RM9.*9.8.*U.*VP3./OMGA;
%弹筒10的运动分析,认为弹筒10在90°周期内做平动:
RM10=68.8398;
RJ10=0.294227;
U=0.15;
RJ10=RM10.*(VP3./OMGA).^2;
F10=0.0;
FM10=-RM10.*9.8.*U.*VP3./OMGA;
%弹筒11的运动分析,认为弹筒11在90°周期内做平动:
RM11=68.8398;
RJ11=0.294227;
U=0.15;
RJ11=RM11.*(VP3./OMGA).^2;
F11=0.0;
FM11=-RM11.*9.8.*U.*VP3./OMGA;
%弹筒12的运动分析,认为弹筒12在90°周期内做平动:
RM12=68.8398;
RJ12=0.294227;
U=0.15;
RJ12=RM12.*(VP3./OMGA).^2;
F12=0.0;
FM12=-RM12.*9.8.*U.*VP3./OMGA;
%弹筒13的运动分析,认为弹筒13在90°周期内做平动:
RM13=68.8398;
RJ13=0.294227;
U=0.15;
RJ13=RM13.*(VP3./OMGA).^2;
F13=0.0;
FM13=-RM13.*9.8.*U.*VP3./OMGA;
%弹筒14的运动分析,认为弹筒14在90°周期内做平动:
RM14=68.8398;
RJ14=0.294227;
U=0.15;
RJ14=RM14.*(VP3./OMGA).^2;
F14=0.0;
FM14=-RM14.*9.8.*U.*VP3./OMGA;
%弹筒15的运动分析,认为弹筒15在90°周期内做平动:
RM15=68.8398;
RJ15=0.294227;
U=0.15;
RJ15=RM15.*(VP3./OMGA).^2;
F15=0.0;
FM15=-RM15.*9.8.*U.*VP3./OMGA;
%弹筒16的运动分析,认为弹筒16在0——45°范围内做完全的平动,在45——90°范围内做平面运动:
%0——45°
RM16=68.8398;
RJ16=0.294227;
U=0.15;
if XTA<=(pi/4.0);
RJ16=RM16.*(VP3./OMGA).^2;
F16=0.0;
FM16=-RM16.*9.8.*U.*VP3./OMGA;
%45——90°
else
OB1=131.0e-3;