Working example

examples/linelast/linelast.force_cantilever.yml

@@ -26,7 +26,7 @@ visualization:
   enabled: false
   unified_paraview: true
   file_path:
-    prefix: lattice_output
+    prefix: paraview_output
 
 parameterized_problem:
   name: linelast_force_cantilever

include/parameterized_problem.hpp

@@ -85,6 +85,13 @@ extern double rdisp_f;
 void init_disp(const Vector &x, Vector &u);
 }
 
+namespace linelast_force
+{
+extern double rforce_f;
+
+void tip_force(const Vector &x, Vector &u);
+}
+
 }
 
 class ParameterizedProblem
@@ -240,6 +247,12 @@ class LinElastDispLattice : public LinElastProblem
    LinElastDispLattice();
 };
 
+class LinElastForceCantilever : public LinElastProblem
+{
+public:
+   LinElastForceCantilever();
+};
+
 ParameterizedProblem* InitParameterizedProblem();
 
 #endif

src/parameterized_problem.cpp

@@ -187,6 +187,18 @@ void init_disp_lcantilever(const Vector &x, Vector &u)
 
 }  // namespace linelast_disp
 
+namespace linelast_force
+{
+
+double rforce_f;
+
+void tip_force(const Vector &x, Vector &f)
+{
+   f = 0.0;
+   f(f.Size()-1) = -1.0e-2 * rforce_f;
+}
+
+}  // namespace linelast_force
 
 }  // namespace function_factory
 
@@ -287,6 +299,10 @@ ParameterizedProblem* InitParameterizedProblem()
    {
       problem = new LinElastDispLattice();
    }
+   else if (problem_name == "linelast_force_cantilever")
+   {
+      problem = new LinElastForceCantilever();
+   }
    else
    {
       mfem_error("Unknown parameterized problem name!\n");
@@ -702,3 +718,56 @@ if (i<2)
    general_vector_ptr[0] = NULL;
 
 }
+
+
+LinElastForceCantilever::LinElastForceCantilever()
+    : LinElastProblem()
+{
+   // pointer to static function.
+   bdr_type.SetSize(3);
+   battr.SetSize(3);
+   vector_bdr_ptr.SetSize(3);
+   for (size_t i = 0; i < vector_bdr_ptr.Size(); i++)
+   {
+   battr[i] = i+1;
+   if (i==0)
+   {
+      bdr_type[i] = LinElastProblem::DIRICHLET;
+      vector_bdr_ptr[i] = &(function_factory::linelast_disp::init_disp);
+   }
+   else if (i==1)
+   {
+      bdr_type[i] = LinElastProblem::NEUMANN;
+      vector_bdr_ptr[i] = &(function_factory::linelast_force::tip_force);
+
+   }
+   else if (i==2)
+   {
+      bdr_type[i] = LinElastProblem::ZERO;
+      vector_bdr_ptr[i] = NULL;
+   }
+
+   }
+
+   // Set materials
+   general_scalar_ptr.SetSize(2);
+   general_scalar_ptr[0] = function_factory::linelast_problem::lambda;
+   general_scalar_ptr[1] = function_factory::linelast_problem::mu;
+
+   // Default values.
+   function_factory::linelast_force::rforce_f = 1.0;
+   function_factory::linelast_problem::_lambda = 1.0;
+   function_factory::linelast_problem::_mu = 1.0;
+
+   param_map["rforce_f"] = 0;
+   param_map["lambda"] = 1;
+   param_map["mu"] = 2;
+
+   param_ptr.SetSize(3);
+   param_ptr[0] = &(function_factory::linelast_force::rforce_f);
+   param_ptr[1] = &(function_factory::linelast_problem::_lambda);
+   param_ptr[2] = &(function_factory::linelast_problem::_mu);
+
+   general_vector_ptr.SetSize(1);
+   general_vector_ptr[0] = NULL;
+}