ssts__tether_mmet_iii_discrete_number_force.m

空间绳系卫星的仿真软件

function [ delta_force_list ]= SSTS__tether_MMET_III_discrete_number_force(N_number,ro,sectional_area,length,Mp,velocity_psi)

% /*M-FILE Function SSTS__tether_MMET_III_discrete_number_force MMM SSTSLAB */
% /*==================================================================================================
%  Simple Space Tether Simulation Laboratory Toolbox for Matlab 7.x
%
%  Copyright 2007 The SxLAB Family - Yi Chen - [e-mail]
% ====================================================================================================
%File description:
%
% to calculate the force for each point mass within the tether for
% MMET-III, based on SSTS__tether_MMET_II_discrete_number.m
%
% Input:
%     N_number       = 20; in MMET-III, number of point masses
%     ro             = 970;           % tether density - kg/(m^3)
%     sectional_area = 62.83*10^(-6); % tether cross-section - m^2
%     length         = 5000;          % tether length - m
%     Mp             = 500;           % kg
%     velocity_psi   = 65000;         % angular velocity - rad/s
% 
% Output:
% force_list -- force list matrix
% 
% Example:
%  [ force_list ]= SSTS__tether_MMET_III_discrete_number_force( 5 )
% 
% force_list =
% 
%   1.0e+016 *
% 
%          0       NaN       NaN       NaN       NaN
%     0.6086    0.6086       NaN       NaN       NaN
%     0.8472         0    0.8472       NaN       NaN
%     0.9732    0.2842    0.2842    0.9732       NaN
%     1.0510    0.4740         0    0.4740    1.0510
% 
% 
%===================================================================================================
%  See Also:   SSTS__tether_MMET_II_discrete_number
% 
%===================================================================================================
%
%===================================================================================================
%Revision -
%Date        Name    Description of Change email                 Location
%22-Jan-2007 Yi Chen Initial version       [e-mail] Glasgow
%HISTORY$
%==================================================================================================*/

% SSTS__tether_MMET_III_discrete_number_force Begin
if ( nargin == 0 )
    N_number       = 20;
    ro             = 970;           % tether density - kg/(m^3)
    sectional_area = 62.83*10^(-6); % tether cross-section - m^2
    length         = 5000;          % tether length - m
    Mp             = 500;           % kg
    velocity_psi   = 65000;         % angular velocity - rad/s
end

if ( nargin == 1 )
    ro             = 970;           % tether density - kg/(m^3)
    sectional_area = 62.83*10^(-6); % tether cross-section - m^2
    length         = 5000;          % tether length - m
    Mp             = 500;           % kg
    velocity_psi   = 65000;         % angular velocity - rad/s
end

% close('all');
% initialize force_list and get location_list
% put location_list input force_list
[ location_list ] = SSTS__tether_MMET_II_discrete_number(N_number)*(length*2);

% 
location_list = abs(location_list - length);

force_list = (velocity_psi^2*location_list).*( Mp + ro*sectional_area*location_list/2 );

force_end_point = velocity_psi^2*length*( Mp + ro*sectional_area*length/2 );
% get the last line 
% force_list = [force_end_point force_list( N_number , : ) force_end_point];
force_list = [ force_list( N_number , : ) ];

%delta_force_list
% force(1)~force(r)
for idx = 1 : 1 : ( N_number/2-1 )
    
    delta_force_list(idx) = force_list(idx) - force_list(idx+1);

end

%force(r+1)
%first half tether
delta_force_list =[delta_force_list force_list(N_number/2)];
% % save('SSTS__tether_MMET_III_Point_Mass_Number.mat', 'N_number'); 
% % SSTS__tether_MMET_III_Point_Mass_Number = load('SSTS__tether_MMET_III_Point_Mass_Number.mat');
% NN = N_number(2,1);
% %full length tether
% delta_force_list =[delta_force_list delta_force_list((N_number+2)/2:1:1)];

%  figure
%  stem (force_list);
%  xlabel('Number of point masses - Number');
%  ylabel('Force on point masses - Force(N) ');

% SSTS__tether_MMET_III_discrete_number_force End