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

        FibreDirections(version) {
            if nargin < 1
                version = 1;
            end
            this = this@models.muscle.AMuscleConfig();
            this.init;

            if version == 2
                this.a0CoordinateSystem= " reference ";
            else
                this.a0CoordinateSystem= " master ";
            end
        }


        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;
            m.DefaultMu(2) = 30;
        }

    
    protected:


        function geo = getGeometry() {
            geo = fem.geometry.RegularHex8Grid([0 1],[0 1],[0 1]);
            pts = geo.Nodes;
            theta = -.3;
            R = [cos(theta) -sin(theta) 0
                 sin(theta) cos(theta)  0
                 0          0           1];
            pts = circshift(R,[1 1])*R*pts;
            geo = fem.geometry.Cube8Node(pts, geo.Elements);
            geo = geo.toCube27Node;
        }
        function displ_dir = setPositionDirichletBC(displ_dir) {
            geo = this.FEM.Geometry;
            displ_dir(:,geo.Elements(1,geo.MasterFaces(1,:))) = true;
            /* displ_dir(1,geo.Elements(1,geo.MasterFaces(1,:))) = true;
             *displ_dir(:,geo.Elements(1,geo.MasterFaces(1,5))) = true; */
        }


        function anull = seta0(anull) {
            anull(1,:,:) = 1;
        }

    
    public: /* ( Static ) */

        static function  test_FibreDirections() {
            f = models.muscle.examples.FibreDirections(1);
            m=f.createModel;
            m.plotGeometrySetup;
            m.simulateAndPlot;

            f = models.muscle.examples.FibreDirections(2);
            m=f.createModel;
            m.plotGeometrySetup;
            m.simulateAndPlot;
        }

    

};
}
}
}