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

        CubePull(varargin) {
            this = this@models.muscle.AMuscleConfig(varargin[:]);
            this.init;
        }
        function  configureModel(m) {
            configureModel@models.muscle.AMuscleConfig(this, m);
            m.T= 150;
            m.dt= .5;
            m.DefaultMu(1) = .5;
            m.DefaultMu(3) = .1;
            m.DefaultInput= 1;
        }


        function  configureModelFinal() {
            m = this.Model;
            m.SpaceReducer.Value= m.System.NumStateDofs;
            configureModelFinal@models.muscle.AMuscleConfig(this);
        }
        function P = getBoundaryPressure(elemidx,faceidx) {
            P = [];
            if elemidx == 1 && faceidx == 3
                P = 1;
            end
        }
        function u = getInputs() {
            u = [this.getAlphaRamp(this.Model.T/3,1)];
        }

    
    protected:


        function geo = getGeometry() {
            geo = fem.geometry.RegularHex8Grid([0 5],-1:10,[0 5]);
            geo = geo.toCube27Node;
            /* geo = geo.toCube20Node; */
        }
        function displ_dir = setPositionDirichletBC(displ_dir) {
            geo = this.FEM.Geometry;
            /*  Fix all on left and only the y,z directions of the back right
             * nodes */
            displ_dir(2,geo.Elements(end,geo.MasterFaces(4,:))) = true;
            displ_dir(:,geo.Elements(end,geo.MasterFaces(4,1))) = true;
        }
        function anull = seta0(anull) {
        }
    public: /* ( Static ) */

        static function  test_CubePull() {
            m = models.muscle.Model(models.muscle.examples.CubePull);
            mu = m.getRandomParam;
            mu(3) = .1;
            m.simulateAndPlot(true,mu,1);
        }


};
}
}
}