Skip to content
Snippets Groups Projects
Commit 4686d5a5 authored by Ludovic Noels's avatar Ludovic Noels
Browse files

comment lines

parent 85aa850f
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
Showing
with 10 additions and 10 deletions
j2plast.c
+ 10
− 10
View file @ 4686d5a5
...@@ -1115,7 +1115,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double ...@@ -1115,7 +1115,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
//POWER LAW HARDENING //POWER LAW HARDENING
case H_POW: case H_POW:
if(p<DP_MIN){ if(p<DP_MIN){
printf("Should not be here: p: %.10e\n",p); //printf("Should not be here: p: %.10e\n",p);
*R=0.; *R=0.;
if(hexp<1.0){ if(hexp<1.0){
*dR = 1e20; *dR = 1e20;
...@@ -1136,7 +1136,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double ...@@ -1136,7 +1136,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
//POWER LAW HARDENING FOR DRUCKER PRAGER //POWER LAW HARDENING FOR DRUCKER PRAGER
case H_POW_DP: case H_POW_DP:
if(p<DP_MIN){ if(p<DP_MIN){
printf("Should not be here: p: %.10e\n",p); //printf("Should not be here: p: %.10e\n",p);
*R=0.; *R=0.;
if(hexp<1.0){ if(hexp<1.0){
*dR = 1e20; *dR = 1e20;
...@@ -1159,7 +1159,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double ...@@ -1159,7 +1159,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
case H_EXP: case H_EXP:
if(p<DP_MIN){ if(p<DP_MIN){
printf("Should not be here: p: %.10e\n",p); //printf("Should not be here: p: %.10e\n",p);
*R=0.; *R=0.;
*dR = hmod1*hexp; *dR = hmod1*hexp;
*ddR = -hexp*(*dR); *ddR = -hexp*(*dR);
...@@ -1178,7 +1178,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double ...@@ -1178,7 +1178,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
case H_EXP_DP: case H_EXP_DP:
if(p<DP_MIN){ if(p<DP_MIN){
printf("Should not be here: p: %.10e\n",p); //printf("Should not be here: p: %.10e\n",p);
*R=0.; *R=0.;
*dR = hmod1*hexp; *dR = hmod1*hexp;
*ddR = -hexp*(*dR); *ddR = -hexp*(*dR);
...@@ -1195,7 +1195,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double ...@@ -1195,7 +1195,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
//SWIFT LAW //SWIFT LAW
case H_SWIFT: case H_SWIFT:
if(p<DP_MIN){ if(p<DP_MIN){
printf("Should not be here: p: %.10e\n",p); //printf("Should not be here: p: %.10e\n",p);
*R=0.; *R=0.;
*dR = sy0*hmod1*hexp; *dR = sy0*hmod1*hexp;
*ddR = (*dR)*(hexp-1.)*hmod1; *ddR = (*dR)*(hexp-1.)*hmod1;
...@@ -1212,7 +1212,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double ...@@ -1212,7 +1212,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
//SWIFT LAW FOR DRUCKER PRAGER //SWIFT LAW FOR DRUCKER PRAGER
case H_SWIFT_DP: case H_SWIFT_DP:
if(p<DP_MIN){ if(p<DP_MIN){
printf("Should not be here: p: %.10e\n",p); //printf("Should not be here: p: %.10e\n",p);
*R=0.; *R=0.;
*dR = sy0*hmod1*hexp; *dR = sy0*hmod1*hexp;
*ddR = (*dR)*(hexp-1.)*hmod1; *ddR = (*dR)*(hexp-1.)*hmod1;
...@@ -1230,7 +1230,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double ...@@ -1230,7 +1230,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
case H_LINEXP: case H_LINEXP:
if(p<DP_MIN){ if(p<DP_MIN){
printf("Should not be here: p: %.10e\n",p); //printf("Should not be here: p: %.10e\n",p);
*R = 0; *R = 0;
*dR = hmod1 + hmod2*hexp; *dR = hmod1 + hmod2*hexp;
*ddR = -hmod2*hexp*hexp; *ddR = -hmod2*hexp*hexp;
...@@ -1247,7 +1247,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double ...@@ -1247,7 +1247,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
case H_LINEXP_DP: case H_LINEXP_DP:
if(p<DP_MIN){ if(p<DP_MIN){
printf("Should not be here: p: %.10e\n",p); //printf("Should not be here: p: %.10e\n",p);
*R = 0; *R = 0;
*dR = hmod1 + hmod2*hexp; *dR = hmod1 + hmod2*hexp;
*ddR = -hmod2*hexp*hexp; *ddR = -hmod2*hexp*hexp;
...@@ -1263,7 +1263,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double ...@@ -1263,7 +1263,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
//POWER LAW HARDENING EXTRAPOLATED AFTER 16% DEFO TO MIMIC DIGIMAT TO ABAQUS //POWER LAW HARDENING EXTRAPOLATED AFTER 16% DEFO TO MIMIC DIGIMAT TO ABAQUS
case H_POW_EXTRAPOL: case H_POW_EXTRAPOL:
if(p<DP_MIN){ if(p<DP_MIN){
printf("Should not be here: p: %.10e\n",p); //printf("Should not be here: p: %.10e\n",p);
*R=0.; *R=0.;
if(hmod1<1){ if(hmod1<1){
*dR = 1e20; *dR = 1e20;
...@@ -1299,7 +1299,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double ...@@ -1299,7 +1299,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
//POWER LAW HARDENING EXTRAPOLATED AFTER 16% DEFO TO MIMIC DIGIMAT TO ABAQUS DRUCKER PRAGER //POWER LAW HARDENING EXTRAPOLATED AFTER 16% DEFO TO MIMIC DIGIMAT TO ABAQUS DRUCKER PRAGER
case H_POW_EXTRAPOL_DP: case H_POW_EXTRAPOL_DP:
if(p<DP_MIN){ if(p<DP_MIN){
printf("Should not be here: p: %.10e\n",p); //printf("Should not be here: p: %.10e\n",p);
*R=0.; *R=0.;
if(hmod1<1){ if(hmod1<1){
*dR = 1e20; *dR = 1e20;
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment