Form of the lagged tangential adhesion disapative potential and forces. More...

#include <TangentialAdhesionForm.hpp>

Inheritance diagram for polyfem::solver::TangentialAdhesionForm:

Collaboration diagram for polyfem::solver::TangentialAdhesionForm:

Public Member Functions
	TangentialAdhesionForm (const ipc::CollisionMesh &collision_mesh, const std::shared_ptr< time_integrator::ImplicitTimeIntegrator > time_integrator, const double epsa, const double mu, const ipc::BroadPhaseMethod broad_phase_method, const NormalAdhesionForm &normal_adhesion_form, const int n_lagging_iters)
	Construct a new Tangential Adhesion Form object.

std::string	name () const override

void	force_shape_derivative (const Eigen::MatrixXd &prev_solution, const Eigen::MatrixXd &solution, const Eigen::MatrixXd &adjoint, const ipc::TangentialCollisions &tangential_constraints_set, Eigen::VectorXd &term)

void	init_lagging (const Eigen::VectorXd &x) override
	Initialize lagged fields.

void	update_lagging (const Eigen::VectorXd &x, const int iter_num) override
	Update lagged fields.

void	update_lagging (const Eigen::VectorXd &x)
	Update lagged fields.

int	max_lagging_iterations () const override
	Get the maximum number of lagging iteration allowable.

bool	uses_lagging () const override
	Does this form require lagging?

Eigen::MatrixXd	compute_displaced_surface (const Eigen::VectorXd &x) const
	Compute the displaced positions of the surface nodes.

Eigen::MatrixXd	compute_surface_velocities (const Eigen::VectorXd &x) const
	Compute the surface velocities.

double	dv_dx () const
	Compute the derivative of the velocities wrt x.

double	mu () const

double	epsa () const

const ipc::TangentialCollisions &	tangential_collision_set () const

const ipc::TangentialAdhesionPotential &	tangential_adhesion_potential () const

Public Member Functions inherited from polyfem::solver::Form
virtual	~Form ()

virtual void	init (const Eigen::VectorXd &x)
	Initialize the form.

virtual void	finish ()

virtual double	value (const Eigen::VectorXd &x) const
	Compute the value of the form multiplied with the weigth.

Eigen::VectorXd	value_per_element (const Eigen::VectorXd &x) const
	Compute the value of the form multiplied with the weigth.

virtual void	first_derivative (const Eigen::VectorXd &x, Eigen::VectorXd &gradv) const
	Compute the first derivative of the value wrt x multiplied with the weigth.

void	second_derivative (const Eigen::VectorXd &x, StiffnessMatrix &hessian) const
	Compute the second derivative of the value wrt x multiplied with the weigth.

virtual bool	is_step_valid (const Eigen::VectorXd &x0, const Eigen::VectorXd &x1) const
	Determine if a step from solution x0 to solution x1 is allowed.

virtual double	max_step_size (const Eigen::VectorXd &x0, const Eigen::VectorXd &x1) const
	Determine the maximum step size allowable between the current and next solution.

virtual void	line_search_begin (const Eigen::VectorXd &x0, const Eigen::VectorXd &x1)
	Initialize variables used during the line search.

virtual void	line_search_end ()
	Clear variables used during the line search.

virtual void	post_step (const polysolve::nonlinear::PostStepData &data)
	Update fields after a step in the optimization.

virtual void	solution_changed (const Eigen::VectorXd &new_x)
	Update cached fields upon a change in the solution.

virtual void	update_quantities (const double t, const Eigen::VectorXd &x)
	Update time-dependent fields.

void	set_project_to_psd (bool val)
	Set project to psd.

bool	is_project_to_psd () const
	Get if the form's second derivative is projected to psd.

void	enable ()
	Enable the form.

void	disable ()
	Disable the form.

void	set_enabled (const bool enabled)
	Set if the form is enabled.

bool	enabled () const
	Determine if the form is enabled.

virtual double	weight () const
	Get the form's multiplicative constant weight.

void	set_weight (const double weight)
	Set the form's multiplicative constant weight.

virtual bool	is_step_collision_free (const Eigen::VectorXd &x0, const Eigen::VectorXd &x1) const
	Checks if the step is collision free.

void	set_output_dir (const std::string &output_dir)

virtual void	set_scale (const double scale)
	sets the scale for the form

Protected Member Functions
double	value_unweighted (const Eigen::VectorXd &x) const override
	Compute the value of the form.

void	first_derivative_unweighted (const Eigen::VectorXd &x, Eigen::VectorXd &gradv) const override
	Compute the first derivative of the value wrt x.

void	second_derivative_unweighted (const Eigen::VectorXd &x, StiffnessMatrix &hessian) const override
	Compute the second derivative of the value wrt x.

Protected Member Functions inherited from polyfem::solver::Form
std::string	resolve_output_path (const std::string &path) const

virtual Eigen::VectorXd	value_per_element_unweighted (const Eigen::VectorXd &x) const
	Compute the value of the form multiplied per element.

Private Attributes
const ipc::CollisionMesh &	collision_mesh_
	Reference to the collision mesh.

const std::shared_ptr< time_integrator::ImplicitTimeIntegrator >	time_integrator_
	Pointer to the time integrator.

const double	epsa_
	Smoothing factor for turning on/off tangential adhesion.

const double	mu_
	Global coefficient of tangential adhesion.

const ipc::BroadPhaseMethod	broad_phase_method_
	Broad-phase method used for distance computation and collision detection.

const int	n_lagging_iters_
	Number of lagging iterations.

ipc::TangentialCollisions	tangential_collision_set_
	Lagged tangential constraint set.

const NormalAdhesionForm &	normal_adhesion_form_
	necessary to have the barrier stiffnes, maybe clean me

const ipc::TangentialAdhesionPotential	tangential_adhesion_potential_

Additional Inherited Members
Protected Attributes inherited from polyfem::solver::Form
bool	project_to_psd_ = false
	If true, the form's second derivative is projected to be positive semidefinite.

double	weight_ = 1
	weight of the form (e.g., AL penalty weight or Δt²)

bool	enabled_ = true
	If true, the form is enabled.

std::string	output_dir_

Detailed Description

Form of the lagged tangential adhesion disapative potential and forces.

Definition at line 18 of file TangentialAdhesionForm.hpp.

Constructor & Destructor Documentation

◆ TangentialAdhesionForm()

polyfem::solver::TangentialAdhesionForm::TangentialAdhesionForm	(	const ipc::CollisionMesh &	collision_mesh,
		const std::shared_ptr< time_integrator::ImplicitTimeIntegrator >	time_integrator,
		const double	epsa,
		const double	mu,
		const ipc::BroadPhaseMethod	broad_phase_method,
		const NormalAdhesionForm &	normal_adhesion_form,
		const int	n_lagging_iters
	)

Construct a new Tangential Adhesion Form object.

Parameters

collision_mesh	Reference to the collision mesh
time_integrator	Pointer to the time integrator
epsa	Smoothing factor between static and dynamic tangential adhesion
mu	Global coefficient of tangential adhesion
dhat	Barrier activation distance
broad_phase_method	Broad-phase method used for distance computation and collision detection
normal_adhesion_form	Pointer to normal adhesion form; necessary to have the potential, maybe clean me
n_lagging_iters	Number of lagging iterations

Definition at line 9 of file TangentialAdhesionForm.cpp.

References epsa_.

Member Function Documentation

◆ compute_displaced_surface()

Eigen::MatrixXd polyfem::solver::TangentialAdhesionForm::compute_displaced_surface ( const Eigen::VectorXd & x ) const

Compute the displaced positions of the surface nodes.

Definition at line 53 of file TangentialAdhesionForm.cpp.

References polyfem::solver::NormalAdhesionForm::compute_displaced_surface(), normal_adhesion_form_, and x.

Referenced by update_lagging().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ compute_surface_velocities()

Eigen::MatrixXd polyfem::solver::TangentialAdhesionForm::compute_surface_velocities ( const Eigen::VectorXd & x ) const

Compute the surface velocities.

Definition at line 58 of file TangentialAdhesionForm.cpp.

References collision_mesh_, time_integrator_, polyfem::utils::unflatten(), and x.

Referenced by first_derivative_unweighted(), second_derivative_unweighted(), and value_unweighted().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ dv_dx()

double polyfem::solver::TangentialAdhesionForm::dv_dx ( ) const

Compute the derivative of the velocities wrt x.

Definition at line 65 of file TangentialAdhesionForm.cpp.

References time_integrator_.

Referenced by second_derivative_unweighted(), and value_unweighted().

Here is the caller graph for this function:

◆ epsa()

double polyfem::solver::TangentialAdhesionForm::epsa ( ) const

inline

Definition at line 87 of file TangentialAdhesionForm.hpp.

References epsa_.

◆ first_derivative_unweighted()

void polyfem::solver::TangentialAdhesionForm::first_derivative_unweighted	(	const Eigen::VectorXd &	x,
		Eigen::VectorXd &	gradv
	)		const

overrideprotectedvirtual

Compute the first derivative of the value wrt x.

Parameters

[in]	x	Current solution
[out]	gradv	Output gradient of the value wrt x

Implements polyfem::solver::Form.

Definition at line 75 of file TangentialAdhesionForm.cpp.

References collision_mesh_, compute_surface_velocities(), tangential_adhesion_potential_, tangential_collision_set_, and x.

Here is the call graph for this function:

◆ force_shape_derivative()

void polyfem::solver::TangentialAdhesionForm::force_shape_derivative	(	const Eigen::MatrixXd &	prev_solution,
		const Eigen::MatrixXd &	solution,
		const Eigen::MatrixXd &	adjoint,
		const ipc::TangentialCollisions &	tangential_constraints_set,
		Eigen::VectorXd &	term
	)

Definition at line 29 of file TangentialAdhesionForm.cpp.

References collision_mesh_, normal_adhesion_form_, polyfem::solver::NormalAdhesionForm::normal_adhesion_potential(), tangential_adhesion_potential_, time_integrator_, and polyfem::utils::unflatten().

Here is the call graph for this function:

◆ init_lagging()

void polyfem::solver::TangentialAdhesionForm::init_lagging ( const Eigen::VectorXd & x )

inlineoverridevirtual

Initialize lagged fields.

Parameters

x	Current solution

Reimplemented from polyfem::solver::Form.

Definition at line 62 of file TangentialAdhesionForm.hpp.

References update_lagging(), and x.

Here is the call graph for this function:

◆ max_lagging_iterations()

int polyfem::solver::TangentialAdhesionForm::max_lagging_iterations ( ) const

inlineoverridevirtual

Get the maximum number of lagging iteration allowable.

Reimplemented from polyfem::solver::Form.

Definition at line 73 of file TangentialAdhesionForm.hpp.

References n_lagging_iters_.

◆ mu()

double polyfem::solver::TangentialAdhesionForm::mu ( ) const

inline

Definition at line 86 of file TangentialAdhesionForm.hpp.

References mu_.

◆ name()

std::string polyfem::solver::TangentialAdhesionForm::name ( ) const

inlineoverridevirtual

Implements polyfem::solver::Form.

Definition at line 39 of file TangentialAdhesionForm.hpp.

◆ second_derivative_unweighted()

void polyfem::solver::TangentialAdhesionForm::second_derivative_unweighted	(	const Eigen::VectorXd &	x,
		StiffnessMatrix &	hessian
	)		const

overrideprotectedvirtual

Compute the second derivative of the value wrt x.

Parameters

[in]	x	Current solution
[out]	hessian	Output Hessian of the value wrt x

Implements polyfem::solver::Form.

Definition at line 82 of file TangentialAdhesionForm.cpp.

References collision_mesh_, compute_surface_velocities(), dv_dx(), POLYFEM_SCOPED_TIMER, polyfem::solver::Form::project_to_psd_, tangential_adhesion_potential_, tangential_collision_set_, and x.

Here is the call graph for this function:

◆ tangential_adhesion_potential()

const ipc::TangentialAdhesionPotential & polyfem::solver::TangentialAdhesionForm::tangential_adhesion_potential ( ) const

inline

Definition at line 89 of file TangentialAdhesionForm.hpp.

References tangential_adhesion_potential_.

◆ tangential_collision_set()

const ipc::TangentialCollisions & polyfem::solver::TangentialAdhesionForm::tangential_collision_set ( ) const

inline

Definition at line 88 of file TangentialAdhesionForm.hpp.

References tangential_collision_set_.

◆ update_lagging() [1/2]

void polyfem::solver::TangentialAdhesionForm::update_lagging ( const Eigen::VectorXd & x )

inline

Update lagged fields.

Parameters

x	Current solution

Definition at line 70 of file TangentialAdhesionForm.hpp.

References update_lagging(), and x.

Referenced by update_lagging().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ update_lagging() [2/2]

void polyfem::solver::TangentialAdhesionForm::update_lagging	(	const Eigen::VectorXd &	x,
		const int	iter_num
	)

overridevirtual

Update lagged fields.

Parameters

x	Current solution

Reimplemented from polyfem::solver::Form.

Definition at line 100 of file TangentialAdhesionForm.cpp.

References broad_phase_method_, collision_mesh_, compute_displaced_surface(), polyfem::solver::NormalAdhesionForm::dhat_a(), mu_, normal_adhesion_form_, polyfem::solver::NormalAdhesionForm::normal_adhesion_potential(), tangential_collision_set_, and x.

Referenced by init_lagging().

Here is the call graph for this function:

Here is the caller graph for this function:

◆ uses_lagging()

bool polyfem::solver::TangentialAdhesionForm::uses_lagging ( ) const

inlineoverridevirtual

Does this form require lagging?

Returns: True if the form requires lagging

Reimplemented from polyfem::solver::Form.

Definition at line 77 of file TangentialAdhesionForm.hpp.

◆ value_unweighted()

double polyfem::solver::TangentialAdhesionForm::value_unweighted ( const Eigen::VectorXd & x ) const

overrideprotectedvirtual

Compute the value of the form.

Parameters

x	Current solution

Returns: Computed value

Implements polyfem::solver::Form.

Definition at line 70 of file TangentialAdhesionForm.cpp.

References collision_mesh_, compute_surface_velocities(), dv_dx(), tangential_adhesion_potential_, tangential_collision_set_, and x.

Here is the call graph for this function:

Member Data Documentation

◆ broad_phase_method_

const ipc::BroadPhaseMethod polyfem::solver::TangentialAdhesionForm::broad_phase_method_

private

Broad-phase method used for distance computation and collision detection.

Definition at line 100 of file TangentialAdhesionForm.hpp.

Referenced by update_lagging().

◆ collision_mesh_

const ipc::CollisionMesh& polyfem::solver::TangentialAdhesionForm::collision_mesh_

private

Reference to the collision mesh.

Definition at line 93 of file TangentialAdhesionForm.hpp.

Referenced by compute_surface_velocities(), first_derivative_unweighted(), force_shape_derivative(), second_derivative_unweighted(), update_lagging(), and value_unweighted().

◆ epsa_

const double polyfem::solver::TangentialAdhesionForm::epsa_

private

Smoothing factor for turning on/off tangential adhesion.

Definition at line 98 of file TangentialAdhesionForm.hpp.

Referenced by epsa(), and TangentialAdhesionForm().

◆ mu_

const double polyfem::solver::TangentialAdhesionForm::mu_

private

Global coefficient of tangential adhesion.

Definition at line 99 of file TangentialAdhesionForm.hpp.

Referenced by mu(), and update_lagging().

◆ n_lagging_iters_

const int polyfem::solver::TangentialAdhesionForm::n_lagging_iters_

private

Number of lagging iterations.

Definition at line 101 of file TangentialAdhesionForm.hpp.

Referenced by max_lagging_iterations().

◆ normal_adhesion_form_

const NormalAdhesionForm& polyfem::solver::TangentialAdhesionForm::normal_adhesion_form_

private

necessary to have the barrier stiffnes, maybe clean me

Definition at line 105 of file TangentialAdhesionForm.hpp.

Referenced by compute_displaced_surface(), force_shape_derivative(), and update_lagging().

◆ tangential_adhesion_potential_

const ipc::TangentialAdhesionPotential polyfem::solver::TangentialAdhesionForm::tangential_adhesion_potential_

private

Definition at line 107 of file TangentialAdhesionForm.hpp.

Referenced by first_derivative_unweighted(), force_shape_derivative(), second_derivative_unweighted(), tangential_adhesion_potential(), and value_unweighted().

◆ tangential_collision_set_

ipc::TangentialCollisions polyfem::solver::TangentialAdhesionForm::tangential_collision_set_

private

Lagged tangential constraint set.

Definition at line 103 of file TangentialAdhesionForm.hpp.

Referenced by first_derivative_unweighted(), second_derivative_unweighted(), tangential_collision_set(), update_lagging(), and value_unweighted().

◆ time_integrator_

const std::shared_ptr<time_integrator::ImplicitTimeIntegrator> polyfem::solver::TangentialAdhesionForm::time_integrator_

private

Pointer to the time integrator.

Definition at line 96 of file TangentialAdhesionForm.hpp.

Referenced by compute_surface_velocities(), dv_dx(), and force_shape_derivative().

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

src/polyfem/solver/forms/TangentialAdhesionForm.hpp
src/polyfem/solver/forms/TangentialAdhesionForm.cpp

Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ TangentialAdhesionForm()

Member Function Documentation

◆ compute_displaced_surface()

◆ compute_surface_velocities()

◆ dv_dx()

◆ epsa()

◆ first_derivative_unweighted()

◆ force_shape_derivative()

◆ init_lagging()

◆ max_lagging_iterations()

◆ mu()

◆ name()

◆ second_derivative_unweighted()

◆ tangential_adhesion_potential()

◆ tangential_collision_set()

◆ update_lagging() [1/2]

◆ update_lagging() [2/2]

◆ uses_lagging()

◆ value_unweighted()

Member Data Documentation

◆ broad_phase_method_

◆ collision_mesh_

◆ epsa_

◆ mu_

◆ n_lagging_iters_

◆ normal_adhesion_form_

◆ tangential_adhesion_potential_

◆ tangential_collision_set_

◆ time_integrator_