Belly.m

Go to the documentation of this file.

namespace models{
namespace muscle{
namespace examples{


/* (Autoinserted by mtoc++)
 * This source code has been filtered by the mtoc++ executable,
 * which generates code that can be processed by the doxygen documentation tool.
 *
 * On the other hand, it can neither be interpreted by MATLAB, nor can it be compiled with a C++ compiler.
 * Except for the comments, the function bodies of your M-file functions are untouched.
 * Consequently, the FILTER_SOURCE_FILES doxygen switch (default in our Doxyfile.template) will produce
 * attached source files that are highly readable by humans.
 *
 * Additionally, links in the doxygen generated documentation to the source code of functions and class members refer to
 * the correct locations in the source code browser.
 * However, the line numbers most likely do not correspond to the line numbers in the original MATLAB source files.
 */

class Belly
  :public models.muscle.AMuscleConfig {
    private:

        NumParts;

        
    public:

        Belly(varargin) {
            this = this@models.muscle.AMuscleConfig(varargin[:]);
            this.init;

            /* % Muscle fibre weights */
            types = [0 .2 .4 .6 .8 1];
            ftw = zeros(this.FEM.GaussPointsPerElem,length(types),this.FEM.Geometry.NumElements);
            /*  Test: Use only slow-twitch muscles */
            ftw(:,1,:) = .4;
            ftw(:,2,:) = .05;
            ftw(:,3,:) = .05;
            ftw(:,4,:) = .1;
            ftw(:,5,:) = .2;
            ftw(:,6,:) = .2;
            this.FibreTypeWeights= ftw;
            p = models.motorunit.Pool;
            p.FibreTypes= types;
            this.Pool= p;
        }


        function  configureModel(m) {
            configureModel@models.muscle.AMuscleConfig(this, m);
            m.T= 150;
            m.dt= .1;
            m.DefaultMu(4) = 6;
            m.Plotter.DefaultArgs= [" Pool ",true];
        }

    
    protected:


        function geo = getGeometry() {
            np = 4;
            geo = fem.geometry.Belly(np,10," Radius ",1," InnerRadius ",.5," Gamma ",2);
            this.NumParts= np;
        }


        function displ_dir = setPositionDirichletBC(displ_dir) {
            geo = this.FEM.Geometry;
            for k = geo.NumElements-3:geo.NumElements
                displ_dir(:,geo.Elements(k,geo.MasterFaces(4,:))) = true;
            end
            for k = 1:4
                displ_dir(:,geo.Elements(k,geo.MasterFaces(3,:))) = true;
            end
        }
        function anull = seta0(anull) {
            anull(2,:,:) = 1;
        }

    
    public: /* ( Static ) */

        static function res = test_BellyGeometryGeneration() {
            g = fem.geometry.Belly(4,35," InnerRadius ",.2," Gamma ",7);
            g = fem.geometry.Belly(4,35," InnerRadius ",[.2 .6]," Gamma ",5);
            g = fem.geometry.Belly(4,35," InnerRadius ",.2," Gamma ",[10 20]);
            g = fem.geometry.Belly(4,35," InnerRadius ",[.2 .6]," Gamma ",[10 20]);
            g = fem.geometry.Belly(4,35," Radius ",[4 2]," InnerRadius ",.5," Gamma ",[10 20]);
            g = fem.geometry.Belly(4,35," Radius ",@(x)[sqrt(abs(x)); 1./(x-34).^2]," InnerRadius ",.2);
            g.plot;
            g = fem.geometry.Belly(4,35," Radius ",@(x)sqrt(abs(x)));
            g.plot;
/*              g = fem.geometry.Belly(4,35,'Radius',[4 2],'InnerRadius',.5,...
 *                 'Gamma',[10 20],'MinZ',-1.5);
 *             g.plot; */
            g = fem.geometry.Belly(4,35," Radius ",@(x)sqrt(abs(x))," Layers ",[.4 .7 1]);
            g.plot;
            res = 1;
        }


        static function  test_BellyModel() {
            m = models.muscle.Model(models.muscle.examples.Belly);
            m.dt= min(m.T/10,2);
            m.simulateAndPlot;
        }

    

};
}
}
}