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

        Variant;

        
    public:

        SprengerUnitCube8Elem(varargin) {
            varargin(end+1:end+2) = [" FL ",2];
            this = this@models.muscle.AMuscleConfig(varargin[:]);
            this.init;
        }


        function  configureModel(m) {
            configureModel@models.muscle.AMuscleConfig(this, m);
            m.T= 1;
            m.dt= .01;
            mu = m.DefaultMu;
            mu([1 2]) = [1 .1];
            m.DefaultMu= mu;

            /* % Material configuration from CMISS/3Elem_sprenger.xml
             * c3M = 3.05907e-10 MPa */
            mu(5) = 3.05907e-10; /*  [MPa]
             * c4M  */

            mu(6) = 47.270456264135881; /*  [-] */


            /*  c01 ^= c1M = 3.56463903963e-02 MPa */
            mu(9) = 3.6463903963e-2; /*  [MPa]
             * c01 ^= c2M = 3.859558659683e-3 MPa */

            mu(10) = 3.859558659683e-3; /*  [MPa] */


            /*  sigma_max_calculation */
            mu(13) = .300; /*  [MPa]
             * Exponential FL-fun */

            mu(14) = .57; /*  [-] */


            m.DefaultMu= mu;
        }

    
    protected:


        function geo = getGeometry() {
            s = load(fullfile(fileparts(which(mfilename))," Sprenger8Elem.mat "));
            geo = s.geo27;
        }


        function displ_dir = setPositionDirichletBC(displ_dir) {
            geo = this.FEM.Geometry;
            displ_dir(:,geo.Elements([1 3 5 7],geo.MasterFaces(1,:))) = true;
            displ_dir(:,geo.Elements([2 4 6 8],geo.MasterFaces(2,:))) = true;
        }
        function anull = seta0(anull) {
            anull(1,:,:) = 1;
        }
    public: /* ( Static ) */

        static function  test_SprengerUnitCube8Elem_ForceComparison() {
            c = models.muscle.examples.SprengerUnitCube8Elem;
            m = c.createModel;
            dfpos = m.getPositionDirichletBCFaceIdx([1 3 5 7],1);
            dfpos2 = m.getPositionDirichletBCFaceIdx([2 4 6 8],2);

            alphas = [.1 .2];/*  .3 .4]; */

            vis = [1 1];/*  4 5]; */

            forces = zeros(2,length(alphas));
            pi = ProcessIndicator(" Running tests for %d alpha values ",length(alphas),false,length(alphas));
            mu = m.DefaultMu;
            for k = 1:length(alphas)
                mu(1) = vis(k);
                c.ActivationRampMax= alphas(k);
                [t,y] = m.simulate(mu);
                df = m.getResidualForces(t,y);
                forces(1,k) = sum(df(dfpos,end));
                forces(2,k) = sum(df(dfpos2,end));
                pi.step;
            end
            pi.stop;
            disp(" Units are mN (milliNewton) ");
            disp(forces);
        }
};
}
}
}