membrane.hpp

Go to the documentation of this file.

#pragma once
 
#include "utils_math.hpp"
#include <string>
 
namespace netdem {
 
class Membrane {
public:
  double radius{0.25};
 
  double height{1.0};
 
  double mesh_size{0.1};
 
  Vec3d center{0, 0, 0};
 
 
  double neo_k{6.94e5};
 
  double neo_mu{5.21e5};
 
  double density{500.0};
 
  double thickness{0.3e-3};
 
  double damp_coef{0.7};
 
  double timestep{1.0e-4};
 
  VecXT<Vec3d> nodes;
 
  VecXT<Vec3i> elements;
 
  VecXT<Vec3d> elemental_stress;
 
  VecXT<double> nodal_vols;
 
  VecXT<Vec3d> nodal_vels;
 
  VecXT<double> bc_facet_pressure;
 
  VecXT<Vec3d> bc_facet_forces;
 
  VecXT<VecNT<double, 6>> bc_nodal_vels;
 
  Membrane(double radius, double height);
 
  Membrane(double radius, double height, double mesh_size);
 
  Membrane(double radius, double height, double mesh_size, double center_x,
           double center_y, double center_z);
 
  void Remesh(double ele_size);
 
  void Init();
 
  void SetBCNodalVelocity(int nid, double vx, double vy, double vz,
                          bool use_prescribed_vx, bool use_prescribed_vy,
                          bool use_prescribed_vz);
 
  void Solve(double dt);
 
  void SaveAsVTK(std::string const &file_name);
 
  ~Membrane();
 
protected:
  double ref_ele_width{0.1}, ref_ele_height{0.1};
 
  double ref_ele_area{0.05};
 
  VecXT<Vec3d> nodal_forces_int;
 
  VecXT<Vec3d> nodal_forces_ext;
 
  void Advance(double dt);
 
  void InitMesh();
 
  void InitInitialCondition();
 
  Mat2d GetDeformationGradient(Vec3d const &v0, Vec3d const &v1,
                               Vec3d const &v2);
 
  Mat2d GetCauchyStress(Mat2d const &def_grad);
 
  std::tuple<Mat3d, Mat3d> GetGlobalForces(Mat2d const &cauchy_stress,
                                           double pressure, Vec3d const &v0,
                                           Vec3d const &v1, Vec3d const &v2);
};
 
} // namespace netdem