Commit d686d00f authored by Asutosh Padhy's avatar Asutosh Padhy
Browse files

Changes made to results section

parent 2b6dc414
Loading
Loading
Loading
Loading
+983 B (19.4 KiB)

File changed.

No diff preview for this file type.

+1.43 KiB (19.9 KiB)

File changed.

No diff preview for this file type.

+55 −14
Original line number Diff line number Diff line
@@ -14,12 +14,12 @@ public class Lag_sme{
          model.param().set("t2", "0.25 [s]"); // Saturation
          model.param().set("t3", "0.5 [s]"); // Reduction
          model.param().set("t4", "0.75 [s]"); // Temperature increase and load set to zero
          model.param().set("t5", "5.0 [s]");
          model.param().set("t5", "2.0 [s]");

          model.param().set("dt", "0.001 [s]");

          model.param().set("T_ref", "295[K]", "Applied Temperature");
          model.param().set("T_rec", "350[K]", "Recovery Temperature");
          model.param().set("T_rec", "320[K]", "Recovery Temperature");

          model.param().set("width", "15[mm]", "Length");
          model.param().set("height", "1.5[mm]", "width");
@@ -84,6 +84,47 @@ public class Lag_sme{
          model.func("rn1").set("type", "normal");
          model.func("rn1").set("funcname", "rn");
          
// Function for plotting Lagoudas Phase modelPath
          
          model.func().create("Ms", "Analytic");
          model.func("Ms").label("Martensite Start boundary");
          model.func("Ms").set("funcname", "Ms");
          model.func("Ms").set("expr", "TMs_1+C_M*(T-TMs_1)");
          model.func("Ms").set("args", "T");
        
          model.func("Ms").setIndex("plotargs", "TMf_1", 0, 1);
          model.func("Ms").setIndex("plotargs", "TAf_1", 0, 2);

          ////
          model.func().create("Mf", "Analytic");
          model.func("Mf").label("Martensite finish boundary");
          model.func("Mf").set("funcname", "Mf");
          model.func("Mf").set("expr", "TMf_1+C_M*(T-TMf_1)");
          model.func("Mf").set("args", "T");
        
          model.func("Mf").setIndex("plotargs", "TMf_1", 0, 1);
          model.func("Mf").setIndex("plotargs", "TAf_1", 0, 2);
          
          
          ////
          model.func().create("As", "Analytic");
          model.func("As").label("Austenite Start boundary");
          model.func("As").set("funcname", "As");
          model.func("As").set("expr", "TAs_1+C_M*(T-TAs_1)");
          model.func("As").set("args", "T");
        
          model.func("As").setIndex("plotargs", "TMf_1", 0, 1);
          model.func("As").setIndex("plotargs", "TAf_1", 0, 2);
          
          ////
          model.func().create("Af", "Analytic");
          model.func("Af").label("Austenite finish boundary");
          model.func("Af").set("funcname", "Af");
          model.func("Af").set("expr", "TAf1_+C_M*(T-TAf_1)");
          model.func("Af").set("args", "T");
        
          model.func("Af").setIndex("plotargs", "TMf_1", 0, 1);
          model.func("Af").setIndex("plotargs", "TAf_1", 0, 2);

// Component block
          model.component().create("comp1", true);
+80 −22
Original line number Diff line number Diff line
@@ -14,12 +14,12 @@ public class Lag_sme_mart{
          model.param().set("t2", "0.25 [s]"); // Saturation
          model.param().set("t3", "0.5 [s]"); // Reduction
          model.param().set("t4", "0.75 [s]"); // Temperature increase and load set to zero
          model.param().set("t5", "5.0 [s]");
          model.param().set("t5", "2.0 [s]");

          model.param().set("dt", "0.001 [s]");

          model.param().set("T_ref", "280[K]", "Applied Temperature");
          model.param().set("T_rec", "350[K]", "Recovery Temperature");
          model.param().set("T_rec", "320[K]", "Recovery Temperature");

          model.param().set("width", "15[mm]", "Length");
          model.param().set("height", "1.5[mm]", "width");
@@ -76,6 +76,47 @@ public class Lag_sme_mart{
              {"t3", "t4", "T_ref"},
              {"t4", "t5", "((T_rec-T_ref)/(t5-t4))*(t-t4)+T_ref"}});

// Function for plotting Lagoudas Phase modelPath
          model.func().create("Ms", "Analytic");
          model.func("Ms").label("Martensite Start boundary");
          model.func("Ms").set("funcname", "Ms");
          model.func("Ms").set("expr", "TMs_1+C_M*(T-TMs_1)");
          model.func("Ms").set("args", "T");
        
          model.func("Ms").setIndex("plotargs", "TMf_1", 0, 1);
          model.func("Ms").setIndex("plotargs", "TAf_1", 0, 2);

          ////
          model.func().create("Mf", "Analytic");
          model.func("Mf").label("Martensite finish boundary");
          model.func("Mf").set("funcname", "Mf");
          model.func("Mf").set("expr", "TMf_1+C_M*(T-TMf_1)");
          model.func("Mf").set("args", "T");
        
          model.func("Mf").setIndex("plotargs", "TMf_1", 0, 1);
          model.func("Mf").setIndex("plotargs", "TAf_1", 0, 2);
          
          
          ////
          model.func().create("As", "Analytic");
          model.func("As").label("Austenite Start boundary");
          model.func("As").set("funcname", "As");
          model.func("As").set("expr", "TAs_1+C_A*(T-TAs_1)");
          model.func("As").set("args", "T");
        
          model.func("As").setIndex("plotargs", "TMf_1", 0, 1);
          model.func("As").setIndex("plotargs", "TAf_1", 0, 2);
          
          
          model.func().create("Af", "Analytic");
          model.func("Af").label("Austenite finish boundary");
          model.func("Af").set("funcname", "Af");
          model.func("Af").set("expr", "TAf_1+C_A*(T-TAf_1)");
          model.func("Af").set("args", "T");
        
          model.func("Af").setIndex("plotargs", "TMf_1", 0, 1);
          model.func("Af").setIndex("plotargs", "TAf_1", 0, 2);
          
          
// Function Random normal distribution

@@ -380,11 +421,13 @@ public class Lag_sme_mart{
          model.result("pg1").create("surf1", "Surface");
          model.result("pg1").feature("surf1").label("Temperature");
          model.result("pg1").feature("surf1").set("data", "dset1");
          model.result("pg1").feature("surf1").active(false);
          
          model.result("pg1").create("surf2", "Surface");
          model.result("pg1").feature("surf2").label("Stress");
          model.result("pg1").feature("surf2").set("data", "dset1");
          model.result("pg1").feature("surf2").set("expr", "solid.mises");
          model.result("pg1").feature("surf2").active(false);

          model.result("pg1").create("surf3", "Surface");
          model.result("pg1").feature("surf3").label("Martensite fraction");
@@ -394,30 +437,45 @@ public class Lag_sme_mart{
          model.result("pg1").create("surf4", "Surface");
          model.result("pg1").feature("surf4").label("Austenite Fraction");
          model.result("pg1").feature("surf4").set("expr", "solid.xi_A");
          model.result("pg1").feature("surf4").active(false);
          model.result("pg1").run();


          model.result().create("pg2", "PlotGroup1D");
          model.result("pg2").label("1D");
          model.result("pg2").label("Strain based");

          model.result("pg2").create("glob1", "Global");
          model.result("pg2").feature("glob1").setIndex("expr", "aveop1(solid.mises)", 0);
          model.result("pg2").feature("glob1").setIndex("expr", "aveop1(solid.sx)", 0);
          model.result("pg2").feature("glob1").set("xdata", "expr");
          model.result("pg2").feature("glob1").set("xdataexpr", "aveop1(solid.eXX)");
          model.result("pg2").create("glob2", "Global");
          model.result("pg2").feature("glob1").label("Stress/strain");

          model.result("pg2").feature("glob2").label("Stress/T");
          model.result("pg2").feature("glob2").set("data", "dset1");
          model.result("pg2").feature("glob2").setIndex("expr", "aveop1(solid.mises)", 0);
          model.result("pg2").feature("glob2").set("xdata", "expr");
          model.result("pg2").feature("glob2").set("xdataexpr", "aveop1(T)");
          model.result("pg2").run();
          
          model.sol("sol1").feature("t1").set("plotgroup", "pg2");
          model.result("pg2").feature("glob2").active(false);
          model.result("pg2").run();
          
          model.result().create("pg3", "PlotGroup1D");
          model.result("pg3").label("Temperature based");
          
          model.result("pg3").create("glob1", "Global");
          model.result("pg3").feature("glob1").label("Stress/T");
          model.result("pg3").feature("glob1").set("data", "dset1");
          
          model.result("pg3").feature("glob1").setIndex("expr", "aveop1(solid.sx)", 0);
          model.result("pg3").feature("glob1").setIndex("expr", "Ms(aveop1(T))", 1);
          model.result("pg3").feature("glob1").setIndex("expr", "Mf(aveop1(T))", 2);
          model.result("pg3").feature("glob1").setIndex("expr", "As(aveop1(T))", 3);
          model.result("pg3").feature("glob1").setIndex("expr", "Af(aveop1(T))", 4);
          
          model.result("pg3").feature("glob1").set("xdata", "expr");
          model.result("pg3").feature("glob1").set("xdataexpr", "aveop1(T)");
          model.result("pg3").set("axislimits", true);
          model.result("pg3").set("xmin", "TMf_1-50");
          model.result("pg3").set("xmax", "TAf_1+50");
          
          model.result("pg3").run();



          return model;
      }