OneDContactConstraint Class Reference

A OneDContactConstraint forces the value of a variable to be the same on both sides of an interface. More...

#include <OneDContactConstraint.h>

Inheritance diagram for OneDContactConstraint:

Public Member Functions
	OneDContactConstraint (const InputParameters &parameters)
virtual	~OneDContactConstraint ()
virtual void	timestepSetup ()
virtual void	jacobianSetup ()
virtual void	updateContactSet ()
virtual Real	computeQpSlaveValue ()
virtual Real	computeQpResidual (Moose::ConstraintType type)
virtual Real	computeQpJacobian (Moose::ConstraintJacobianType type)
bool	shouldApply ()

Protected Attributes
NumericVector< Number > &	_residual_copy
bool	_jacobian_update

Detailed Description

A OneDContactConstraint forces the value of a variable to be the same on both sides of an interface.

Definition at line 26 of file OneDContactConstraint.h.

Constructor & Destructor Documentation

OneDContactConstraint()

OneDContactConstraint::OneDContactConstraint

(

const InputParameters &

parameters

)

Definition at line 36 of file OneDContactConstraint.C.

  : NodeFaceConstraint(parameters),
    _residual_copy(_sys.residualGhosted()),
    _jacobian_update(getParam<bool>("jacobian_update"))
{
}

~OneDContactConstraint()

virtual OneDContactConstraint::~OneDContactConstraint ( )

inlinevirtual

Definition at line 30 of file OneDContactConstraint.h.

{}

Member Function Documentation

computeQpJacobian()

Real OneDContactConstraint::computeQpJacobian ( Moose::ConstraintJacobianType  type )

virtual

Definition at line 121 of file OneDContactConstraint.C.

{
  double slave_jac = 0;
  switch (type)
  {
    case Moose::SlaveSlave:
      return _phi_slave[_j][_qp] * _test_slave[_i][_qp];
    case Moose::SlaveMaster:
      return -_phi_master[_j][_qp] * _test_slave[_i][_qp];
    case Moose::MasterSlave:
      slave_jac =
          (*_jacobian)(_current_node->dof_number(0, _var.number(), 0), _connected_dof_indices[_j]);
      return slave_jac * _test_master[_i][_qp];
    case Moose::MasterMaster:
      return 0;
  }
  return 0;
}

computeQpResidual()

Real OneDContactConstraint::computeQpResidual ( Moose::ConstraintType  type )

virtual

Definition at line 104 of file OneDContactConstraint.C.

{
  PenetrationInfo * pinfo = _penetration_locator._penetration_info[_current_node->id()];

  switch (type)
  {
    case Moose::Slave:
      // return (_u_slave[_qp] - _u_master[_qp])*_test_slave[_i][_qp];
      return ((*_current_node)(0) - pinfo->_closest_point(0)) * _test_slave[_i][_qp];
    case Moose::Master:
      double slave_resid = _residual_copy(_current_node->dof_number(0, _var.number(), 0));
      return slave_resid * _test_master[_i][_qp];
  }
  return 0;
}

computeQpSlaveValue()

Real OneDContactConstraint::computeQpSlaveValue ( )

virtual

Definition at line 90 of file OneDContactConstraint.C.

{
  PenetrationInfo * pinfo = _penetration_locator._penetration_info[_current_node->id()];
  Moose::err << std::endl
             << "Popping out node: " << _current_node->id() << std::endl
             << "Closest Point x: " << pinfo->_closest_point(0) << std::endl
             << "Current Node x: " << (*_current_node)(0) << std::endl
             << "Current Value: " << _u_slave[_qp] << std::endl
             << std::endl;

  return pinfo->_closest_point(0) - ((*_current_node)(0) - _u_slave[_qp]);
}

jacobianSetup()

void OneDContactConstraint::jacobianSetup ( )

virtual

Definition at line 50 of file OneDContactConstraint.C.

{
  if (_jacobian_update)
    updateContactSet();
}

shouldApply()

bool OneDContactConstraint::shouldApply

(

)

Definition at line 82 of file OneDContactConstraint.C.

{
  std::set<dof_id_type>::iterator hpit =
      _penetration_locator._has_penetrated.find(_current_node->id());
  return (hpit != _penetration_locator._has_penetrated.end());
}

timestepSetup()

void OneDContactConstraint::timestepSetup ( )

virtual

Definition at line 44 of file OneDContactConstraint.C.

{
  updateContactSet();
}

updateContactSet()

void OneDContactConstraint::updateContactSet ( )

virtual

Definition at line 57 of file OneDContactConstraint.C.

Referenced by jacobianSetup(), and timestepSetup().

{
  std::set<dof_id_type> & has_penetrated = _penetration_locator._has_penetrated;

  std::map<dof_id_type, PenetrationInfo *>::iterator
      it = _penetration_locator._penetration_info.begin(),
      end = _penetration_locator._penetration_info.end();

  for (; it != end; ++it)
  {
    PenetrationInfo * pinfo = it->second;

    // Skip this pinfo if there are no DOFs on this node.
    if (!pinfo || pinfo->_node->n_comp(_sys.number(), _var.number()) < 1)
      continue;

    if (pinfo->_distance > 0)
    {
      dof_id_type slave_node_num = it->first;
      has_penetrated.insert(slave_node_num);
    }
  }
}

Member Data Documentation

_jacobian_update

bool OneDContactConstraint::_jacobian_update

protected

Definition at line 48 of file OneDContactConstraint.h.

Referenced by jacobianSetup().

_residual_copy

NumericVector<Number>& OneDContactConstraint::_residual_copy

protected

Definition at line 46 of file OneDContactConstraint.h.

Referenced by computeQpResidual().

---

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

• contact/include/OneDContactConstraint.h
• contact/src/OneDContactConstraint.C