simtk.org

Code: Select all

  CMCmodel = Model(RRAUpdatedModel);
    CMCmodel.initSystem(); 
    
    %     % cant make a call to the current probe, no function getProbe?
    %     Probe = CMCmodel.getProbe()
    

    %  % adding a probe adds as metabolics_0, but does not add probes to
    %  each muscle that the previous model already contains. 
%     activationMaintenanceRateOn = 'True';
%     shorteningRateOn = 'True';
%     basalRateOn = 'False';
%     mechanicalWorkRateOn = 'True';
%     reportTotalMetabolicsOnly = 'False';
    wholeBodyProbe = Umberger2010MuscleMetabolicsProbe(); %activationMaintenanceRateOn, shorteningRateOn, basalRateOn, mechanicalWorkRateOn);
    wholeBodyProbe.set_enabled(1)
    wholeBodyProbe.isEnabled() % prints enable status in command window
    wholeBodyProbe.setName("metabolics");
    CMCmodel.addProbe(wholeBodyProbe)
    
    % save model
    CMCModel = strcat(CMCFolder, '\', Subjects(subj).name, '_', Subjects(subj).Trials(trial).name, '_CMC.osim');
    CMCmodel.print(CMCModel);

Code: Select all

No constructor 'org.opensim.modeling.Umberger2010MuscleMetabolicsProbe' with matching signature found.

Error in ABL_Batch_OpenSim (line 498)
    wholeBodyProbe = Umberger2010MuscleMetabolicsProbe(activationMaintenanceRateOn, shorteningRateOn,
    basalRateOn, mechanicalWorkRateOn);

Code: Select all

No method 'set_report_total_metabolics_only' with matching signature found for class
'org.opensim.modeling.Umberger2010MuscleMetabolicsProbe'.

Error in ABL_Batch_OpenSim (line 504)
    wholeBodyProbe.set_report_total_metabolics_only(reportTotalMetabolicsOnly)

Code: Select all

wholeBodyProbe.set_report_total_metabolics_only(false);

Code: Select all

 clc;
        import org.opensim.modeling.*
        import java.io.*
        
        % define model for CMC
        CMCmodel = Model(RRAUpdatedModel);
        CMCmodel.initSystem();
        
        % get slow twitch data from text file
        fn = 'C:\Users\richa\Documents\OpenSim 4.0\Metabolix\Scripts\metabolicsSlowTwitchRatios_Gait2392.txt';
        fdata = importdata(fn);
        muscleName = fdata.textdata;
        twitchRatio = fdata.data;
        
        % The following booleans are constructor arguments for the Umberger probe.
        % These settings are used for all probes.
        activationMaintenanceRateOn = true;
        shorteningRateOn = true;
        basalRateOn = false;
        mechanicalWorkRateOn = true;
        reportTotalMetabolicsOnly = false;
        
        % The mass of each muscle will be calculated using data from the model:
        % muscleMass = (maxIsometricForce / sigma) * rho * optimalFiberLength
        % where sigma = 0.25e6 is the specific tension of mammalian muscle (in
        % Pascals) and rho = 1059.7 is the density of mammalian muscle (in kg/m^3).
        
        % The slow-twitch ratio used for muscles that either do not appear in the
        % file, or appear but whose proportion of slow-twitch fibers is unknown.
        defaultTwitchRatio = 0.5;
        % Get the slow-twitch ratio as the default value.
        slowTwitchRatio = ones(1,length(twitchRatio)) * defaultTwitchRatio;
        
        % Define a whole-body probe that will report the total metabolic energy
        % consumption over the simulation.
        wholeBodyProbe = Umberger2010MuscleMetabolicsProbe(...
            activationMaintenanceRateOn,...
            shorteningRateOn,...
            basalRateOn,...
            mechanicalWorkRateOn);
        wholeBodyProbe.setOperation("value");
        wholeBodyProbe.set_report_total_metabolics_only(reportTotalMetabolicsOnly);
        
        % Add the probe to the model and provide a name.
        CMCmodel.addProbe(wholeBodyProbe);
        wholeBodyProbe.setName("metabolics");
        
        % Loop through all muscles, adding parameters for each into the whole-body probe.
        for j = 1:CMCmodel.getMuscles().getSize()
            
            % offset for python/C++ looping
            iMuscle = j -1;
            % save muscle name in text file via diary function
            diary('temp.txt');
            thisMuscle = CMCmodel.getMuscles().get(iMuscle)
            diary off
            
            % get muscle name and save
            FID = fopen('temp.txt');
            TXT = textscan(FID, '%s');
            if strcmp(TXT{1}(end), 'off')
                TXT{1}(end) = [];
            end
            if strcmp(TXT{1}(end), 'diary')
                TXT{1}(end) = [];
            end
            txt = TXT{1}(end);
            Muscle = txt{1}(1:end-2);
            
            % Set the slow-twitch ratio to the physiological value, if it is known.
            %     for i = 1:length(twitchRatio) % in twitchRatios.items():
            ind = find(contains({muscleName{:}}, Muscle));
            if strcmp(Muscle, muscleName{ind}) && twitchRatio(ind) ~= -1
                slowTwitchRatio(ind) = twitchRatio(ind);
            end
            
            % Add this muscle to the whole-body probe. The arguments are muscle
            % name, slow-twitch ratio, and muscle mass. Note that the muscle mass
            % is ignored unless we set useProvidedMass to True.
            wholeBodyProbe.addMuscle(thisMuscle.getName(),...
                slowTwitchRatio(ind));
            
        end
        
        % Save the new model to a file      
        CMCModelPath = strcat(CMCFolder, '\', Subjects(subj).name, '_', Subjects(subj).Trials(trial).name, '_CMC.osim');
        CMCmodel.print(CMCModelPath);