ThinTendon.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 ThinTendon
  :public models.muscle.AExperimentModelConfig {
    public:

        ylen = 200;

        movetime = 80;

        configs;

        
    public:

        ThinTendon(varargin) {
            this = this@models.muscle.AExperimentModelConfig(varargin[:]);
            this.addOption(" BC ",1);
            this.init;
            this.VelocityBCTimeFun= general.functions.ConstantUntil(this.movetime);
            this.configs= (-.9:.1:.2)*this.ylen/this.movetime;
            this.NumConfigurations= length(this.configs);
            this.NumOutputs= 1;
        }


        function  configureModel(m) {
            configureModel@models.muscle.AMuscleConfig(this, m);
            m.T= this.movetime+20;
            m.dt= m.T/200;

            mu = m.DefaultMu;
            /*  Small viscosity */
            mu(3) = 0;
            m.DefaultMu= mu;
        }


        function tmr = getTendonMuscleRatio(points) {
            tmr = ones(1,size(points,2));            
        }
        function u = getInputs() {
            u = [];
            if this.Options.BC == 2
                u = [this.getAlphaRamp(this.movetime,1)];
            end
        }
        function P = getBoundaryPressure(elemidx,faceidx) {
            P = [];
            if this.Options.BC == 2 && any(elemidx == 21:24) && faceidx == 4
                P = 1;
            end
        }
        function o = getOutputOfInterest(double t,matrix<double> y) {
            m = this.Model;
            /*  Get the residual dirichlet forces at end time */
            df = m.getResidualForces(t(end),y(:,end));
            idx = m.getPositionDirichletBCFaceIdx(1:4,3);
            o(1) = sum(df(idx));
            /* o(2) = sum(df(idx,passive_pos));    */
        }

    
    protected:


        function geo = getGeometry() {
            geo = fem.geometry.Belly(6,this.ylen," Radius ",0," InnerRadius ",2," Gamma ",2);
            n = geo.Nodes;
            /*  Slightly deviate a node in the center */
            n(1,137) = n(1,137)*1.1;
            geo = fem.geometry.Cube27Node(n,geo.Elements);
            /* geo = belly.scale(20); */
        }


        function displ_dir = setPositionDirichletBC(displ_dir) {
            geo = this.FEM.Geometry;
            /*  Fix ends in xz direction */
            displ_dir(:,geo.Elements(1:4,geo.MasterFaces(3,:))) = true;
            displ_dir([1 3],geo.Elements(21:24,geo.MasterFaces(4,:))) = true;
        }
        function [velo_dir , velo_dir_val ] = setVelocityDirichletBC(velo_dir,velo_dir_val) {
            if this.Options.BC == 1
                geo = this.FEM.Geometry;
                /*  Fix ends in xz direction */
                velo_dir(2,geo.Elements(21:24,geo.MasterFaces(4,:))) = true;
                velo_dir_val(velo_dir) = this.configs(this.CurrentConfigNr);
            end
        }
    public: /* ( Static ) */

        static function  test_ThinTendon() {
            c = models.muscle.examples.ThinTendon;
            m = c.createModel;
            c.CurrentConfigNr= 4;
            m.simulateAndPlot;
        }


        static function  runExperiment() {
            c = models.muscle.examples.ThinTendon;
            m = c.createModel;
            e = models.muscle.ExperimentRunner(m);
            d = e.runExperimentsCached;
            semilogy(d.o);
            title(" Forces at fixed nodes at end time ");
        }

    

};
}
}
}