Shaker.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 Shaker
  :public models.muscle.AMuscleConfig {
    public:

        ylen;

        
    public:

        Shaker(varargin) {
            this = this@models.muscle.AMuscleConfig(varargin[:]);
            this.init;
            this.VelocityBCTimeFun= general.functions.Sinus(50); /*  50Hz */

        }


        function  configureModel(m) {
            configureModel@models.muscle.AMuscleConfig(this, m);
            m.T= 75;
            m.dt= .5;

            mu = m.DefaultMu;
            /*  Small viscosity */
            mu(1) = .005;
            m.DefaultMu= mu;
        }


    
    protected:


        function geo = getGeometry() {
            belly = fem.geometry.Belly(4,10," Radius ",1," InnerRadius ",.4," Gamma ",2);
            geo = belly.scale(20);
            this.ylen= 100;
        }


        function displ_dir = setPositionDirichletBC(displ_dir) {
            geo = this.FEM.Geometry;
            /*  Fix ends in xz direction */
            displ_dir([1 3],geo.Elements(1:4,geo.MasterFaces(3,:))) = true;
            displ_dir([1 3],geo.Elements(13:16,geo.MasterFaces(4,:))) = true;
        }
        function [velo_dir , velo_dir_val ] = setVelocityDirichletBC(velo_dir,velo_dir_val) {
            geo = this.FEM.Geometry;
            /*  Fix ends in xz direction */
            velo_dir(2,geo.Elements(1:4,geo.MasterFaces(3,:))) = true;
            velo_dir(2,geo.Elements(13:16,geo.MasterFaces(4,:))) = true;
/*              for k = [1 2 7 8 5 6 11 12]  
 *                 pos = geo.Elements(k,geo.MasterFaces(3,:));
 *                 velo_dir(1,pos) = true;
 *                 %velo_dir_val(1,pos) = 1;
 *             end */
            velo_dir_val(velo_dir) = 2;
        }
        function anull = seta0(anull) {
            anull(2,:,:) = 1;
        }
    public: /* ( Static ) */

        static function  test_Shaker() {
            m = models.muscle.Model(models.muscle.examples.Shaker);
            m.simulateAndPlot;
        }

    

};
}
}
}