AMuscleConfig.m

Go to the documentation of this file.

namespace models{
namespace muscle{


/* (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 AMuscleConfig
  :public fem.AFEMConfig {
    public:

        PressureFEM;

        
    public:

        FibreTypeWeights = "[]";


        Pool;


        char a0CoordinateSystem = "master";
    public:

        double ActivationRampMax = 1;
        double ActivationRampOffset = 0;
    private:

        iP;

        optArgs = {""};

        
    public:

        AMuscleConfig(varargin) {
            this = this@fem.AFEMConfig(varargin[:]);
            /*  Force-length function type */
            this.addOption(" FL ",1);
        }


        function m = createModel() {
            m = models.muscle.Model(this);
        }
        function  configureModelFinal() {
            configureModelFinal@fem.AFEMConfig(this);
            /*  Automatically set the reducable dimensions for u,v
             * separately. */
            m = this.Model;
            s = m.System;
            /* targetd = s.NumStateDofs + (1:s.NumDerivativeDofs); */
            targetd = 1:s.NumStateDofs;
            /*  Dont reduce those dofs subject to velocity BCs! */
            targetd(s.DerivativeDirichletPosInStateDofs) = [];
            m.SpaceReducer.TargetDimensions= targetd;
        }


        function u = getInputs() {
            u = [];
        }
        function x0 = getX0(x0) {
        }
        function  setForceLengthFun(f) {

            /*  Steps to produce derivative of exponential version
             *             lw = sym('lw'); rw = sym('rw'); lexpo = sym('lexpo'); rexpo = sym('rexpo'); ratio = sym('ratio');
             *             f(ratio) = exp(-((1-ratio)/lw).^lexpo);
             *             df = diff(f)
             *             f2(ratio) = exp(-((ratio-1)/rw).^rexpo);
             *             df2 = diff(f2) */

            switch this.Options.FL
                case 1
                    /*  Linear (Gordon 66)
                     * Exported here to separate class for tidyness */
                    g = models.muscle.functions.Gordon66SarcoForceLength(f.mu(14));
                    [fun, dfun] = g.getFunction;
                case 2
                    /*  Exponential (Schmitt) */
                    lexpo = 4; /*  4 */

                    rexpo = 3; /*  3 */

                    lw = f.mu(14); /*  orig .57 */

                    rw = f.mu(14)*1.3; /*  orig .14 */

                    fun = @(ratio)(ratio<=1).*exp(-((1-ratio)/lw).^lexpo) ...
                        + (ratio>1).*exp(-((ratio-1)/rw).^rexpo);
                    dfun = @(ratio)(ratio<=1).*((lexpo*exp(-(-(ratio - 1)/lw)^lexpo)*(-(ratio - 1)/lw)^(lexpo - 1))/lw) ...
                        + (ratio > 1) .* (-(rexpo*exp(-((ratio - 1)/rw)^rexpo)*((ratio - 1)/rw)^(rexpo - 1))/rw);
                case 3
                    /*  Quadratic Polynomial (Heidlauf) */
                    fun = @(ratio)(-6.25*ratio.*ratio + 12.5*ratio - 5.25) .* (ratio >= .6) .* (ratio <= 1.4);
                    dfun = @(ratio)(12.5*ratio.*(1-ratio)) .* (ratio >= .6) .* (ratio <= 1.4);
            end
            f.ForceLengthFun= fun;
            f.ForceLengthFunDeriv= dfun;
        }
        function alpha = getAlphaRamp(ramptime,double alphamax,double starttime) {
            if nargin < 4
                starttime = this.ActivationRampOffset;
                if nargin < 3
                    alphamax = this.ActivationRampMax;
                end
            end
            ramp = general.functions.Ramp(ramptime, alphamax, starttime);
            alpha = ramp.getFunction;
        }
        function tmr = getTendonMuscleRatio(unused1) {
            tmr = []; /*  zeros(1,size(x,2)); */

        }
    protected:


        function  init() {
            init@fem.AFEMConfig(this);

            /* % Get the geometry */
            geo = this.Geometry;
            if isa(geo," fem.geometry.Cube8Node ")
                press_geo = geo;
            elseif isa(geo," fem.geometry.Cube20Node ") || isa(geo," fem.geometry.Cube27Node ")
                press_geo = geo.toCube8Node;
            else
                error(" Scenario not yet implemented for geometry class '%s' ", class(geo));
            end
            this.PressureFEM= fem.HexahedronTrilinear(press_geo);
            /* this.PressureFEM = fem.HexahedronSerendipity(press_geo.toCube20Node);
             *this.PressureFEM = fem.HexahedronTriquadratic(press_geo.toCube27Node); */
        }


        function anull = seta0(anull) {
        }
        function ftw = getFibreTypeWeights() {
            fe = this.FEM;
            geo = fe.Geometry;
            ftw = zeros(fe.GaussPointsPerElem,length(this.FibreTypes),geo.NumElements);
        }
    public: /* ( Sealed ) */

        function anull = geta0() {
            fe = this.FEM;
            g = this.Geometry;
            anull = zeros(3,fe.GaussPointsPerElem,g.NumElements);
            anull = this.seta0(anull);
            /*  Normalize anull vectors */
            for m = 1:g.NumElements
                anull(:,:,m) = anull(:,:,m) ./ ([1;1;1]*Norm.L2(anull(:,:,m)));
            end
        }

    

};
}
}