PolyFEM
Loading...
Searching...
No Matches
Public Member Functions | Private Member Functions | Private Attributes | List of all members
polyfem::assembler::NeoHookeanElasticity Class Reference

#include <NeoHookeanElasticity.hpp>

Inheritance diagram for polyfem::assembler::NeoHookeanElasticity:
[legend]
Collaboration diagram for polyfem::assembler::NeoHookeanElasticity:
[legend]

Public Member Functions

 NeoHookeanElasticity ()
 
double compute_energy (const NonLinearAssemblerData &data) const override
 
Eigen::VectorXd assemble_gradient (const NonLinearAssemblerData &data) const override
 
Eigen::MatrixXd assemble_hessian (const NonLinearAssemblerData &data) const override
 
VectorNd compute_rhs (const AutodiffHessianPt &pt) const override
 
void compute_stiffness_value (const double t, const assembler::ElementAssemblyValues &vals, const Eigen::MatrixXd &local_pts, const Eigen::MatrixXd &displacement, Eigen::MatrixXd &tensor) const override
 
void compute_stress_grad_multiply_mat (const OptAssemblerData &data, const Eigen::MatrixXd &mat, Eigen::MatrixXd &stress, Eigen::MatrixXd &result) const override
 
void compute_stress_grad_multiply_stress (const OptAssemblerData &data, Eigen::MatrixXd &stress, Eigen::MatrixXd &result) const override
 
void compute_stress_grad_multiply_vect (const OptAssemblerData &data, const Eigen::MatrixXd &vect, Eigen::MatrixXd &stress, Eigen::MatrixXd &result) const override
 
void compute_dstress_dmu_dlambda (const OptAssemblerData &data, Eigen::MatrixXd &dstress_dmu, Eigen::MatrixXd &dstress_dlambda) const override
 
void add_multimaterial (const int index, const json &params, const Units &units) override
 
void set_params (const LameParameters &params)
 
const LameParameters & lame_params () const
 
void update_lame_params (const Eigen::MatrixXd &lambdas, const Eigen::MatrixXd &mus) override
 
std::string name () const override
 
std::map< std::string, ParamFunc > parameters () const override
 
void assign_stress_tensor (const OutputData &data, const int all_size, const ElasticityTensorType &type, Eigen::MatrixXd &all, const std::function< Eigen::MatrixXd(const Eigen::MatrixXd &)> &fun) const override
 
template<int n_basis, int dim>
void compute_energy_aux_gradient_fast (const NonLinearAssemblerData &data, Eigen::Matrix< double, Eigen::Dynamic, 1 > &G_flattened) const
 
double assemble_energy (const bool is_volume, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev) const override
 
void assemble_gradient (const bool is_volume, const int n_basis, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev, Eigen::MatrixXd &rhs) const override
 
virtual Eigen::VectorXd assemble_gradient (const NonLinearAssemblerData &data) const=0
 
void assemble_hessian (const bool is_volume, const int n_basis, const bool project_to_psd, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev, utils::MatrixCache &mat_cache, StiffnessMatrix &grad) const override
 
virtual Eigen::MatrixXd assemble_hessian (const NonLinearAssemblerData &data) const=0
 
- Public Member Functions inherited from polyfem::assembler::NLAssembler
virtual ~NLAssembler ()=default
 
double assemble_energy (const bool is_volume, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev) const override
 
Eigen::VectorXd assemble_energy_per_element (const bool is_volume, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev) const override
 
void assemble_gradient (const bool is_volume, const int n_basis, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev, Eigen::MatrixXd &rhs) const override
 
void assemble_hessian (const bool is_volume, const int n_basis, const bool project_to_psd, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, const double dt, const Eigen::MatrixXd &displacement, const Eigen::MatrixXd &displacement_prev, utils::MatrixCache &mat_cache, StiffnessMatrix &grad) const override
 
virtual bool is_linear () const override
 
- Public Member Functions inherited from polyfem::assembler::Assembler
virtual ~Assembler ()=default
 
int size () const
 
virtual void set_size (const int size)
 
virtual void assemble (const bool is_volume, const int n_basis, const std::vector< basis::ElementBases > &bases, const std::vector< basis::ElementBases > &gbases, const AssemblyValsCache &cache, const double t, StiffnessMatrix &stiffness, const bool is_mass=false) const
 
virtual void compute_stress_grad (const OptAssemblerData &data, const Eigen::MatrixXd &prev_grad_u_i, Eigen::MatrixXd &stress, Eigen::MatrixXd &result) const
 
virtual void compute_stress_prev_grad (const OptAssemblerData &data, const Eigen::MatrixXd &prev_grad_u_i, Eigen::MatrixXd &result) const
 
virtual Eigen::Matrix< AutodiffScalarGrad, Eigen::Dynamic, 1, 0, 3, 1 > kernel (const int dim, const AutodiffGradPt &rvect, const AutodiffScalarGrad &r) const
 
void set_materials (const std::vector< int > &body_ids, const json &body_params, const Units &units)
 
virtual bool is_fluid () const
 
- Public Member Functions inherited from polyfem::assembler::ElasticityAssembler
 ElasticityAssembler ()
 
virtual ~ElasticityAssembler ()=default
 
void compute_scalar_value (const OutputData &data, std::vector< NamedMatrix > &result) const override
 
void compute_tensor_value (const OutputData &data, std::vector< NamedMatrix > &result) const override
 
void compute_stress_tensor (const OutputData &data, const ElasticityTensorType &type, Eigen::MatrixXd &stresses) const
 
void compute_von_mises_stresses (const OutputData &data, Eigen::MatrixXd &stresses) const
 
bool is_solution_displacement () const override
 
bool is_tensor () const override
 

Private Member Functions

template<typename T >
T compute_energy_aux (const NonLinearAssemblerData &data) const
 
template<int n_basis, int dim>
void compute_energy_hessian_aux_fast (const NonLinearAssemblerData &data, Eigen::MatrixXd &H) const
 
template<int n_basis, int dim>
void compute_energy_aux_gradient_fast (const NonLinearAssemblerData &data, Eigen::VectorXd &G_flattened) const
 

Private Attributes

LameParameters params_
 

Additional Inherited Members

- Public Types inherited from polyfem::assembler::Assembler
typedef std::pair< std::string, Eigen::MatrixXd > NamedMatrix
 
typedef std::function< double(const RowVectorNd &, const RowVectorNd &, double, int)> ParamFunc
 
- Protected Attributes inherited from polyfem::assembler::Assembler
int size_ = -1
 

Detailed Description

Definition at line 11 of file NeoHookeanElasticity.hpp.

Constructor & Destructor Documentation

◆ NeoHookeanElasticity()

polyfem::assembler::NeoHookeanElasticity::NeoHookeanElasticity ( )

Definition at line 15 of file NeoHookeanElasticity.cpp.

Member Function Documentation

◆ add_multimaterial()

void polyfem::assembler::NeoHookeanElasticity::add_multimaterial ( const int  index,
const json &  params,
const Units &  units 
)
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

Definition at line 19 of file NeoHookeanElasticity.cpp.

References polyfem::assembler::LameParameters::add_multimaterial(), params_, polyfem::assembler::Assembler::size(), and polyfem::Units::stress().

Referenced by polyfem::assembler::MultiModel::add_multimaterial().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ assemble_energy()

double polyfem::assembler::NLAssembler::assemble_energy ( const bool  is_volume,
const std::vector< basis::ElementBases > &  bases,
const std::vector< basis::ElementBases > &  gbases,
const AssemblyValsCache &  cache,
const double  t,
const double  dt,
const Eigen::MatrixXd &  displacement,
const Eigen::MatrixXd &  displacement_prev 
) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

Definition at line 244 of file Assembler.cpp.

◆ assemble_gradient() [1/3]

void polyfem::assembler::NLAssembler::assemble_gradient ( const bool  is_volume,
const int  n_basis,
const std::vector< basis::ElementBases > &  bases,
const std::vector< basis::ElementBases > &  gbases,
const AssemblyValsCache &  cache,
const double  t,
const double  dt,
const Eigen::MatrixXd &  displacement,
const Eigen::MatrixXd &  displacement_prev,
Eigen::MatrixXd &  rhs 
) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

◆ assemble_gradient() [2/3]

Eigen::VectorXd polyfem::assembler::NeoHookeanElasticity::assemble_gradient ( const NonLinearAssemblerData &  data) const
overridevirtual

Implements polyfem::assembler::NLAssembler.

Definition at line 47 of file NeoHookeanElasticity.cpp.

References polyfem::assembler::ElementAssemblyValues::basis_values, polyfem::assembler::Assembler::size(), and polyfem::assembler::NonLinearAssemblerData::vals.

Referenced by polyfem::assembler::MultiModel::assemble_gradient().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ assemble_gradient() [3/3]

virtual Eigen::VectorXd polyfem::assembler::NLAssembler::assemble_gradient ( const NonLinearAssemblerData &  data) const
virtual

Implements polyfem::assembler::NLAssembler.

◆ assemble_hessian() [1/3]

void polyfem::assembler::NLAssembler::assemble_hessian ( const bool  is_volume,
const int  n_basis,
const bool  project_to_psd,
const std::vector< basis::ElementBases > &  bases,
const std::vector< basis::ElementBases > &  gbases,
const AssemblyValsCache &  cache,
const double  t,
const double  dt,
const Eigen::MatrixXd &  displacement,
const Eigen::MatrixXd &  displacement_prev,
utils::MatrixCache &  mat_cache,
StiffnessMatrix &  grad 
) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

◆ assemble_hessian() [2/3]

Eigen::MatrixXd polyfem::assembler::NeoHookeanElasticity::assemble_hessian ( const NonLinearAssemblerData &  data) const
overridevirtual

Implements polyfem::assembler::NLAssembler.

Definition at line 174 of file NeoHookeanElasticity.cpp.

References polyfem::assembler::ElementAssemblyValues::basis_values, polyfem::assembler::Assembler::size(), and polyfem::assembler::NonLinearAssemblerData::vals.

Referenced by polyfem::assembler::MultiModel::assemble_hessian().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ assemble_hessian() [3/3]

virtual Eigen::MatrixXd polyfem::assembler::NLAssembler::assemble_hessian ( const NonLinearAssemblerData &  data) const
virtual

Implements polyfem::assembler::NLAssembler.

◆ assign_stress_tensor()

void polyfem::assembler::NeoHookeanElasticity::assign_stress_tensor ( const OutputData &  data,
const int  all_size,
const ElasticityTensorType &  type,
Eigen::MatrixXd &  all,
const std::function< Eigen::MatrixXd(const Eigen::MatrixXd &)> &  fun 
) const
overridevirtual

Implements polyfem::assembler::ElasticityAssembler.

Definition at line 251 of file NeoHookeanElasticity.cpp.

References polyfem::assembler::OutputData::bs, polyfem::assembler::ElementAssemblyValues::compute(), polyfem::compute_diplacement_grad(), polyfem::assembler::OutputData::el_id, polyfem::F, polyfem::assembler::OutputData::fun, polyfem::assembler::OutputData::gbs, polyfem::assembler::LameParameters::lambda_mu(), polyfem::assembler::OutputData::local_pts, params_, polyfem::PK1, polyfem::pk1_from_cauchy(), polyfem::PK2, polyfem::pk2_from_cauchy(), polyfem::assembler::Assembler::size(), polyfem::assembler::OutputData::t, and vals.

Referenced by polyfem::assembler::MultiModel::assign_stress_tensor().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ compute_dstress_dmu_dlambda()

void polyfem::assembler::NeoHookeanElasticity::compute_dstress_dmu_dlambda ( const OptAssemblerData &  data,
Eigen::MatrixXd &  dstress_dmu,
Eigen::MatrixXd &  dstress_dlambda 
) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

Definition at line 660 of file NeoHookeanElasticity.cpp.

References polyfem::assembler::OptAssemblerData::el_id, polyfem::assembler::OptAssemblerData::global_pts, polyfem::assembler::OptAssemblerData::grad_u_i, polyfem::assembler::OptAssemblerData::local_pts, and polyfem::assembler::OptAssemblerData::t.

◆ compute_energy()

double polyfem::assembler::NeoHookeanElasticity::compute_energy ( const NonLinearAssemblerData &  data) const
overridevirtual

Implements polyfem::assembler::NLAssembler.

Definition at line 300 of file NeoHookeanElasticity.cpp.

Referenced by polyfem::assembler::MultiModel::compute_energy().

Here is the caller graph for this function:

◆ compute_energy_aux()

template<typename T >
T polyfem::assembler::NeoHookeanElasticity::compute_energy_aux ( const NonLinearAssemblerData &  data) const
private

Definition at line 307 of file NeoHookeanElasticity.cpp.

References polyfem::compute_disp_grad_at_quad(), polyfem::assembler::NonLinearAssemblerData::da, polyfem::utils::determinant(), polyfem::assembler::ElementAssemblyValues::element_id, polyfem::get_local_disp(), polyfem::assembler::LameParameters::lambda_mu(), params_, polyfem::quadrature::Quadrature::points, polyfem::assembler::ElementAssemblyValues::quadrature, polyfem::assembler::Assembler::size(), polyfem::assembler::NonLinearAssemblerData::t, val, polyfem::assembler::ElementAssemblyValues::val, and polyfem::assembler::NonLinearAssemblerData::vals.

Here is the call graph for this function:

◆ compute_energy_aux_gradient_fast() [1/2]

template<int n_basis, int dim>
void polyfem::assembler::NeoHookeanElasticity::compute_energy_aux_gradient_fast ( const NonLinearAssemblerData &  data,
Eigen::Matrix< double, Eigen::Dynamic, 1 > &  G_flattened 
) const

Definition at line 371 of file NeoHookeanElasticity.cpp.

References polyfem::assembler::ElementAssemblyValues::basis_values, polyfem::assembler::NonLinearAssemblerData::da, polyfem::assembler::ElementAssemblyValues::element_id, polyfem::assembler::ElementAssemblyValues::jac_it, polyfem::assembler::LameParameters::lambda_mu(), params_, polyfem::quadrature::Quadrature::points, polyfem::assembler::ElementAssemblyValues::quadrature, polyfem::assembler::Assembler::size(), polyfem::assembler::NonLinearAssemblerData::t, val, polyfem::assembler::ElementAssemblyValues::val, polyfem::assembler::NonLinearAssemblerData::vals, and polyfem::assembler::NonLinearAssemblerData::x.

Here is the call graph for this function:

◆ compute_energy_aux_gradient_fast() [2/2]

template<int n_basis, int dim>
void polyfem::assembler::NeoHookeanElasticity::compute_energy_aux_gradient_fast ( const NonLinearAssemblerData &  data,
Eigen::VectorXd &  G_flattened 
) const
private

◆ compute_energy_hessian_aux_fast()

template<int n_basis, int dim>
void polyfem::assembler::NeoHookeanElasticity::compute_energy_hessian_aux_fast ( const NonLinearAssemblerData &  data,
Eigen::MatrixXd &  H 
) const
private

Definition at line 462 of file NeoHookeanElasticity.cpp.

References polyfem::assembler::ElementAssemblyValues::basis_values, polyfem::assembler::NonLinearAssemblerData::da, polyfem::assembler::ElementAssemblyValues::element_id, polyfem::assembler::ElementAssemblyValues::jac_it, polyfem::assembler::LameParameters::lambda_mu(), params_, polyfem::quadrature::Quadrature::points, polyfem::assembler::ElementAssemblyValues::quadrature, polyfem::assembler::Assembler::size(), polyfem::assembler::NonLinearAssemblerData::t, val, polyfem::assembler::ElementAssemblyValues::val, polyfem::assembler::NonLinearAssemblerData::vals, and polyfem::assembler::NonLinearAssemblerData::x.

Here is the call graph for this function:

◆ compute_rhs()

Eigen::Matrix< double, Eigen::Dynamic, 1, 0, 3, 1 > polyfem::assembler::NeoHookeanElasticity::compute_rhs ( const AutodiffHessianPt &  pt) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

Definition at line 27 of file NeoHookeanElasticity.cpp.

References polyfem::assembler::LameParameters::lambda_mu(), polyfem::autogen::neo_hookean_2d_function(), polyfem::autogen::neo_hookean_3d_function(), params_, and polyfem::assembler::Assembler::size().

Here is the call graph for this function:

◆ compute_stiffness_value()

void polyfem::assembler::NeoHookeanElasticity::compute_stiffness_value ( const double  t,
const assembler::ElementAssemblyValues &  vals,
const Eigen::MatrixXd &  local_pts,
const Eigen::MatrixXd &  displacement,
Eigen::MatrixXd &  tensor 
) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

Definition at line 116 of file NeoHookeanElasticity.cpp.

References polyfem::compute_diplacement_grad(), polyfem::utils::flatten(), polyfem::assembler::LameParameters::lambda_mu(), params_, polyfem::assembler::Assembler::size(), and vals.

Here is the call graph for this function:

◆ compute_stress_grad_multiply_mat()

void polyfem::assembler::NeoHookeanElasticity::compute_stress_grad_multiply_mat ( const OptAssemblerData &  data,
const Eigen::MatrixXd &  mat,
Eigen::MatrixXd &  stress,
Eigen::MatrixXd &  result 
) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

Definition at line 573 of file NeoHookeanElasticity.cpp.

References polyfem::assembler::OptAssemblerData::el_id, polyfem::assembler::OptAssemblerData::global_pts, polyfem::assembler::OptAssemblerData::grad_u_i, polyfem::assembler::LameParameters::lambda_mu(), polyfem::assembler::OptAssemblerData::local_pts, mat, params_, and polyfem::assembler::OptAssemblerData::t.

Here is the call graph for this function:

◆ compute_stress_grad_multiply_stress()

void polyfem::assembler::NeoHookeanElasticity::compute_stress_grad_multiply_stress ( const OptAssemblerData &  data,
Eigen::MatrixXd &  stress,
Eigen::MatrixXd &  result 
) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

Definition at line 595 of file NeoHookeanElasticity.cpp.

References polyfem::assembler::OptAssemblerData::el_id, polyfem::assembler::OptAssemblerData::global_pts, polyfem::assembler::OptAssemblerData::grad_u_i, polyfem::assembler::LameParameters::lambda_mu(), polyfem::assembler::OptAssemblerData::local_pts, params_, and polyfem::assembler::OptAssemblerData::t.

Here is the call graph for this function:

◆ compute_stress_grad_multiply_vect()

void polyfem::assembler::NeoHookeanElasticity::compute_stress_grad_multiply_vect ( const OptAssemblerData &  data,
const Eigen::MatrixXd &  vect,
Eigen::MatrixXd &  stress,
Eigen::MatrixXd &  result 
) const
overridevirtual

Reimplemented from polyfem::assembler::Assembler.

Definition at line 616 of file NeoHookeanElasticity.cpp.

References polyfem::assembler::OptAssemblerData::el_id, polyfem::assembler::OptAssemblerData::global_pts, polyfem::assembler::OptAssemblerData::grad_u_i, polyfem::assembler::LameParameters::lambda_mu(), polyfem::assembler::OptAssemblerData::local_pts, params_, polyfem::assembler::Assembler::size(), and polyfem::assembler::OptAssemblerData::t.

Here is the call graph for this function:

◆ lame_params()

const LameParameters & polyfem::assembler::NeoHookeanElasticity::lame_params ( ) const
inline

Definition at line 56 of file NeoHookeanElasticity.hpp.

References params_.

◆ name()

std::string polyfem::assembler::NeoHookeanElasticity::name ( ) const
inlineoverridevirtual

Implements polyfem::assembler::Assembler.

Definition at line 64 of file NeoHookeanElasticity.hpp.

◆ parameters()

std::map< std::string, Assembler::ParamFunc > polyfem::assembler::NeoHookeanElasticity::parameters ( ) const
overridevirtual

Implements polyfem::assembler::Assembler.

Definition at line 677 of file NeoHookeanElasticity.cpp.

References params_, and polyfem::assembler::Assembler::size().

Here is the call graph for this function:

◆ set_params()

void polyfem::assembler::NeoHookeanElasticity::set_params ( const LameParameters &  params)
inline

Definition at line 55 of file NeoHookeanElasticity.hpp.

References params_.

◆ update_lame_params()

void polyfem::assembler::NeoHookeanElasticity::update_lame_params ( const Eigen::MatrixXd &  lambdas,
const Eigen::MatrixXd &  mus 
)
inlineoverridevirtual

Reimplemented from polyfem::assembler::Assembler.

Definition at line 58 of file NeoHookeanElasticity.hpp.

References polyfem::assembler::LameParameters::lambda_mat_, polyfem::assembler::LameParameters::mu_mat_, and params_.

Member Data Documentation

◆ params_

LameParameters polyfem::assembler::NeoHookeanElasticity::params_
private

Definition at line 74 of file NeoHookeanElasticity.hpp.

Referenced by add_multimaterial(), assign_stress_tensor(), compute_energy_aux(), compute_energy_aux_gradient_fast(), compute_energy_hessian_aux_fast(), compute_rhs(), compute_stiffness_value(), compute_stress_grad_multiply_mat(), compute_stress_grad_multiply_stress(), compute_stress_grad_multiply_vect(), lame_params(), parameters(), set_params(), and update_lame_params().

The documentation for this class was generated from the following files: