#include "tutorials/adv-1_contact/src/TestingContact.hpp"

Inheritance diagram for ContactOps::HertzAxisymmetricPressureNormCheck:

Collaboration diagram for ContactOps::HertzAxisymmetricPressureNormCheck:

Public Member Functions
	ContactPressureNormCheck (MoFEM::Interface &m_field, double tol_norm, std::string check_name)

Public Member Functions inherited from ContactOps::ContactPressureNormCheck
	ContactPressureNormCheck (MoFEM::Interface &m_field, double tol_norm, std::string check_name)

virtual	~ContactPressureNormCheck ()=default

MoFEMErrorCode	operator() (DM dm)

Protected Member Functions
double	analyticalPressure (double x, double y, double) const override

Additional Inherited Members
Protected Attributes inherited from ContactOps::ContactPressureNormCheck
MoFEM::Interface &	mField

double	tolNorm

std::string	checkName

Detailed Description

Test 5: Hertz sphere-on-plane contact, axisymmetric.

Definition at line 157 of file TestingContact.hpp.

Member Function Documentation

◆ analyticalPressure()

double ContactOps::HertzAxisymmetricPressureNormCheck::analyticalPressure	(	double	x,
		double	y,
		double	z
	)		const

inlineoverrideprotectedvirtual

Analytical pressure at a boundary gauss point's physical coordinates.

Implements ContactOps::ContactPressureNormCheck.

Definition at line 161 of file TestingContact.hpp.

                                                                       {
    double E_star = young_modulus / (1 - poisson_ratio * poisson_ratio);
    double R = 100.;
    double d = 0.01;
    double F = (4. / 3.) * E_star * std::sqrt(R) * std::pow(d, 1.5);
    double a = std::pow((3. * F * R) / (4. * E_star), 1. / 3.);
    double p_max = (3. * F) / (2. * M_PI * a * a);
    double r = std::sqrt((x * x) + (y * y));
    if (r > a)
      return 0.;
    return p_max * std::sqrt(1 - ((r * r) / (a * a)));
  }