Ejemplo n.º 1
0
Real
SolidMechX::computeQpResidual()
  {
    if (!_constant_properties)
      recomputeConstants();

    _strain(0,0) = _grad_u[_qp](0);
    _strain(1,1) = _grad_y[_qp](1);
    _strain(0,1) = _grad_u[_qp](1)+_grad_y[_qp](0);

    _stress(0) = _c1*_strain(0,0)+_c1*_c2*_strain(1,1);
    _stress(1) = _c1*_c3*_strain(0,1);

    if ( 3 == _mesh_dimension){
      _strain(2,2) = _grad_z[_qp](2);
      _strain(0,2) = _grad_u[_qp](2)+_grad_z[_qp](0);

      _stress(0) += _c1*_c2*_strain(2,2);
      _stress(2) = _c1*_c3*_strain(0,2);
    }

    Real value = (_stress*_grad_test[_i][_qp]);

    return value;
  }
Ejemplo n.º 2
0
Real
SolidMechZ::computeQpResidual()
  {
    if (!_constant_properties)
      recomputeConstants();

    _strain(0,0) = _grad_x[_qp](0);
    _strain(1,1) = _grad_y[_qp](1);
    _strain(2,2) = _grad_u[_qp](2);
    _strain(1,2) = _grad_u[_qp](1)+_grad_y[_qp](2);
    _strain(0,2) = _grad_x[_qp](2)+_grad_u[_qp](0);

    _stress(0) = _c1*_c3*_strain(0,2);
    _stress(1) = _c1*_c3*_strain(1,2);
    _stress(2) = _c1*_c2*_strain(0,0)+_c1*_c2*_strain(1,1)+_c1*_strain(2,2);

    return (_stress*_grad_test[_i][_qp]);

  }