ssts__tether_mmet_iii_discrete_number_stress.m

空间绳系卫星的仿真软件

function [ stress_list ]= SSTS__tether_MMET_III_discrete_number_stress(N_number,ro,sectional_area,length,Mp,velocity_psi)

% /*M-FILE Function SSTS__tether_MMET_III_discrete_number_stress 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 stress 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_stress( 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_stress 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

[stress_list] = SSTS__tether_MMET_III_discrete_number_force(N_number,ro,sectional_area,length,Mp,velocity_psi)./sectional_area;

%  figure
%  stem (stress_list);
%  xlabel('Number of point masses - Number');
%  ylabel('Stress on point masses - Stress(N/m) ');

% SSTS__tether_MMET_III_discrete_number_stress End