Skip to content
Snippets Groups Projects
Commit b2494a14 authored by Ujwal Kishore Jinaga's avatar Ujwal Kishore Jinaga :clown:
Browse files

fixed several bugs in extraBranch tangents TVP. dPdF fixed. model converges

parent 48a67f5b
Branches
Tags
1 merge request!390AddedNumericalTangentstoCheck_dGamma_dDgamma_dCepr_dT_derivatives___PLEASEcomm...
...@@ -125,14 +125,7 @@ void mlawHyperViscoElastic::evaluatePhiPCorrection(double tr, const STensor3 &de ...@@ -125,14 +125,7 @@ void mlawHyperViscoElastic::evaluatePhiPCorrection(double tr, const STensor3 &de
else if(method == 1) else if(method == 1)
{ {
A_v = getVolumeCorrection()*(tanh(getXiVolumeCorrection()/3.*tr*tr-getZetaVolumeCorrection())+tanh(getZetaVolumeCorrection())); A_v = getVolumeCorrection()*(tanh(getXiVolumeCorrection()/3.*tr*tr-getZetaVolumeCorrection())+tanh(getZetaVolumeCorrection()));
// if(A_v > -1.){
dA_vdE = getVolumeCorrection()*getXiVolumeCorrection()*2./3.*(1.-tanh(getXiVolumeCorrection()/3.*tr*tr-getZetaVolumeCorrection())*tanh(getXiVolumeCorrection()/3.*tr*tr-getZetaVolumeCorrection())); dA_vdE = getVolumeCorrection()*getXiVolumeCorrection()*2./3.*(1.-tanh(getXiVolumeCorrection()/3.*tr*tr-getZetaVolumeCorrection())*tanh(getXiVolumeCorrection()/3.*tr*tr-getZetaVolumeCorrection()));
// }
/*else{
A_v = -0.999;
dA_vdE = 0.;
}*/
intA=0.; intA=0.;
if (getXiVolumeCorrection()!=0.) if (getXiVolumeCorrection()!=0.)
...@@ -142,21 +135,11 @@ void mlawHyperViscoElastic::evaluatePhiPCorrection(double tr, const STensor3 &de ...@@ -142,21 +135,11 @@ void mlawHyperViscoElastic::evaluatePhiPCorrection(double tr, const STensor3 &de
} }
B_d = getDevCorrection()*(tanh(getThetaDevCorrection()*dev.dotprod()-getPiDevCorrection())+tanh(getPiDevCorrection())); // 1. - B_d = getDevCorrection()*(tanh(getThetaDevCorrection()*dev.dotprod()-getPiDevCorrection())+tanh(getPiDevCorrection()));
STensorOperation::zero(dB_vddev); STensorOperation::zero(dB_vddev);
// if(B_d > -1.){
dB_vddev=dev; dB_vddev=dev;
dB_vddev*=2.*getPiDevCorrection()*getDevCorrection()*(1.-tanh(getThetaDevCorrection()*dev.dotprod()-getPiDevCorrection())*tanh(getThetaDevCorrection()*dev.dotprod()-getPiDevCorrection()) ); dB_vddev*=2.*getPiDevCorrection()*getDevCorrection()*(1.-tanh(getThetaDevCorrection()*dev.dotprod()-getPiDevCorrection())*tanh(getThetaDevCorrection()*dev.dotprod()-getPiDevCorrection()) );
// }
/*else{
B_d = -0.999;
STensorOperation::zero(dB_vddev);
}*/
// Msg::Error(" Inside evaluatePhiPCorrection, A_v = %e, B_d = %e !!", A_v, B_d);
intB=0.; intB=0.;
if(getThetaDevCorrection()!=0) if(getThetaDevCorrection()!=0)
{ {
......
...@@ -1150,26 +1150,6 @@ void mlawNonLinearTVE::extraBranchLaw(const STensor3& Ee, const double& T, const ...@@ -1150,26 +1150,6 @@ void mlawNonLinearTVE::extraBranchLaw(const STensor3& Ee, const double& T, const
*DsigDEe = _Idev; *DsigDEe = _Idev;
*DsigDEe *= (2.*Gextra*B); *DsigDEe *= (2.*Gextra*B);
STensorOperation::prodAdd(devEe, dB_dDevEe, 2.*Gextra, *DsigDEe); STensorOperation::prodAdd(devEe, dB_dDevEe, 2.*Gextra, *DsigDEe);
/*
static STensor43 temp1, temp2;
STensorOperation::zero(temp1);
STensorOperation::zero(temp2);
for (int i=0; i<3; i++)
for (int j=0; j<3; j++)
for (int k=0; k<3; k++)
for (int l=0; l<3; l++){
temp1(i,j,k,l) = 2.*Gextra*devEe(i,j)*dB_dDevEe(k,l);
}
for (int i=0; i<3; i++)
for (int j=0; j<3; j++)
for (int k=0; k<3; k++)
for (int l=0; l<3; l++)
for (int p=0; p<3; p++)
for (int q=0; q<3; q++){
temp2(i,j,k,l) += temp1(i,j,p,q)*_Idev(p,q,k,l);
}
*DsigDEe += temp2;*/
STensorOperation::prodAdd(_I,_I,Kextra*(A+trEe*dA_dTrEe),*DsigDEe); STensorOperation::prodAdd(_I,_I,Kextra*(A+trEe*dA_dTrEe),*DsigDEe);
} }
......
This diff is collapsed.
This diff is collapsed.
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment