SiebertTendonFun.m

Go to the documentation of this file.

namespace models{
namespace muscle{
namespace functions{


/* (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 SiebertTendonFun
  :public general.functions.AFunGen {
    public:

        l0;

                lm = .742;
lam_m;

            F1 = 1345;
    public:

        SiebertTendonFun(l0) {
            if nargin < 1
                l0 = [];
            end
            if ~isempty(l0)
                this.l0= l0;
                this.lam_m= 1+this.lm/l0;
            end
        }


        function [fhandle , dfhandle ] = getFunction() {
            k = 4.36;
            F_1 = this.F1;
            c = F_1/(exp(k)-1);
            K = 8030; /*  [mN/mm] */


            lnull = this.l0;
            if ~isempty(lnull)
                lm = this.lam_m;
                fhandle = @(l)(l<=lm)*c.*(exp(k*(l-1)/(lm-1))-1) + (l > lm).*(F_1+K*(l-lm)*lnull);
            else
                lm = .742;
                fhandle = @(l)(l<=lm)*c.*(exp(k*l/lm)-1) + (l > lm).*(F_1+K*(l-lm));
            end
            dfhandle = [];
        }


        function str = getConfigStr() {
            if ~isempty(this.l0)
                str = sprintf(" l0: %g, \\lambda_m: %g ",this.l0,this.lam_m);
            else
                str = sprintf(" Original from paper Siebert 2012 ");
            end
        }


        function  plot(range) {
            point = this.lm;
            if ~isempty(this.l0)
                point = this.lam_m;
                if nargin < 2
                    range = [1 1.2*point];
                end
                lab = " \lambda [-] ";
            else
                if nargin < 2
                    range = [0 4];
                end
                lab = " \Delta stretch [mm] ";
            end
            plot@general.functions.AFunGen(this, range);

/*              ax = get(gcf,'Children');
 *             ax = ax(2); % second one is the left one - hope this is reproducible */

            ax = gca;
            hold(ax," on ");
            plot(ax,point,this.F1," rx "," MarkerSize ",16);
            xlabel(ax,lab);
            ylabel(ax," force [mN] ");
        }

    

};
}
}
}