void CosseratLinearElasticMaterial::computeQpStress() { //Calculation and apply stress free strain RankTwoTensor stress_free_strain = computeStressFreeStrain(); _elastic_strain[_qp] += stress_free_strain; // stress = E * e _stress[_qp] = _elasticity_tensor[_qp] * _elastic_strain[_qp]; _stress_couple[_qp] = _elastic_flexural_rigidity_tensor[_qp] * _curvature[_qp]; }
void LinearElasticMaterial::computeQpStress() { //Calculation and Apply stress free strain RankTwoTensor stress_free_strain = computeStressFreeStrain(); // add the stress free strain on here // ther derivatives of elastic_strain w.r.t. c are built down in EigenstrainBaseMaterial _elastic_strain[_qp] += stress_free_strain; _total_strain[_qp] = _elastic_strain[_qp]; // stress = C * e _stress[_qp] = _elasticity_tensor[_qp] * _elastic_strain[_qp]; }