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

        Shear(varargin) {
            this = this@models.muscle.AMuscleConfig(varargin[:]);
            this.addOption(" BC ",1);
            this.init;
        }


        function  configureModel(m) {
            configureModel@models.muscle.AMuscleConfig(this, m);
            m.T= 50;
            m.dt= 1;
            m.ODESolver.RelTol= 1e-6;
            m.ODESolver.AbsTol= 1e-6;
            if this.Options.BC == 2
                m.DefaultInput= 1;
                m.DefaultMu(3) = 1;
            end
        }


        function P = getBoundaryPressure(elemidx,faceidx) {
            P = [];
            o = this.Options;
            if o.BC == 2
                if faceidx == 2 && ((o.GeoNr == 1 && elemidx == 1) || ...
                        (o.GeoNr == 2 && any(elemidx == 5:8)))
                   P = zeros(3,3);
                   P(3,1) = 1;
                end
            end
        }
        function u = getInputs() {
            u = [];
            if this.Options.BC == 2
                u = [this.getAlphaRamp(this.Model.T/2,.1)];
            end
        }

    
    protected:


        function geo = getGeometry() {
            if this.Options.GeoNr == 1
                geo = fem.geometry.RegularHex8Grid([0 1],[0 1],[0 1]);
            else
                geo = fem.geometry.RegularHex8Grid(0:2,0:2,0:2);
            end
            geo = geo.toCube27Node;
        }
        function displ_dir = setPositionDirichletBC(displ_dir) {
            geo = this.FEM.Geometry;
            o = this.Options;
            if o.GeoNr == 1
                displ_dir(:,geo.Elements(1,geo.MasterFaces(1,:))) = true;
                if o.BC == 1
                    displ_dir(1:2,geo.Elements(1,geo.MasterFaces(2,:))) = true;
                end
            else
                displ_dir(:,geo.Elements(1:4,geo.MasterFaces(1,:))) = true;
                if o.BC == 1
                    displ_dir(1:2,geo.Elements(5:8,geo.MasterFaces(2,:))) = true;
                end
            end
        }


        function [velo_dir , velo_dir_val ] = setVelocityDirichletBC(velo_dir,velo_dir_val) {
            geo = this.FEM.Geometry;
            o = this.Options;
            if o.BC == 1
                if o.GeoNr == 1
                    velo_dir(3,geo.Elements(1,geo.MasterFaces(2,:))) = true;
                else
                    velo_dir(3,geo.Elements(5:8,geo.MasterFaces(2,:))) = true;
                end
                velo_dir_val(velo_dir) = .1;
            end
        }
        function anull = seta0(anull) {
            anull(3,:,:) = 1;
        }

    
    public: /* ( Static ) */

        static function  test_Shear() {
            for g = 1:2
                for b = 1:2
                    f = models.muscle.examples.Shear(" GeoNr ",g," BC ",b);
                    m = f.createModel;
                    m.plotGeometrySetup;
                    [t,y] = m.simulate;
                    m.plot(t,y," Lambdas ",true);
                end
            end
        }

    

};
}
}
}