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

Changes to Sme

parent da95f9bb

Java/Lagoudas/bin/Lag.class

+252 B (17 KiB)

File changed.

No diff preview for this file type.

View original file
View changed file

Java/Lagoudas/bin/Lag_sme.class

+1.25 KiB (18.1 KiB)

File changed.

No diff preview for this file type.

View original file
View changed file

Java/Lagoudas/src/Lag.java

+6 −0
Original line number Diff line number Diff line
@@ -119,6 +119,12 @@ public class Lag { 
          model.component("comp1").physics("ht").feature("sma1").selection().set(1);

          model.component("comp1").physics("ht").feature("sma1").set("EffectiveConductivity", "PowerLaw");



          model.component("comp1").physics("ht").feature("sma1").set("MartensiteMaterial", "mat2");

          model.component("comp1").physics("ht").feature("sma1").set("AusteniteMaterial", "mat1");

          model.component("comp1").physics("ht").feature("sma1").set("xi", "0.0");

          model.component("comp1").physics("ht").feature("sma1").set("rho_mat", "userdef");

          model.component("comp1").physics("ht").feature("sma1").set("rho", "6500");          



          // temp1

          model.component("comp1").physics("ht").create("temp1", "TemperatureBoundary", 1);

          model.component("comp1").physics("ht").feature("temp1").selection().set(1, 4);

Java/Lagoudas/src/Lag_sme.java

+129 −105
Original line number Diff line number Diff line
@@ -17,7 +17,7 @@ public class Lag_sme{ 
          model.param().set("t5", "0.56[s]");

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



          model.param().set("T_ref", "260[K]", "Applied Temperature");

          model.param().set("T_ref", "305[K]", "Applied Temperature");

          model.param().set("para", "0", "Continuation paramater");

          model.param().set("width", "15.1[mm]", "Length");

          model.param().set("height", "1.67[mm]", "width");
@@ -25,6 +25,26 @@ public class Lag_sme{ 
          model.param().set("d", "0.034[mm]", "thickness");

          model.param().set("maxdisp", "0.85 [mm]", "Strain dependent");



          model.param().set("p_ratio", "0.33", "Poissonsratio");

          model.param().set("rho_", "6500", "Density");

          model.param().set("Ymod_A", "41.2e9", "Youngs Modulus Austenite");

          model.param().set("Ymod_M", "32.9e9", "Youngs Modulus Martensite");

          model.param().set("Hcap", "400", "Heat capacity");

          model.param().set("kappa_A", "18.0", "Thermal conductivity Austenite");

          model.param().set("kappa_M", "8.6", "Thermal conductivity Martensite");

          model.param().set("hflux", "18.0", "Heat Flux coefficient");



          model.param().set("TMs_1", "291", "Martensite start temperature");

          model.param().set("TMf_1", "271", "Martensite finish temperature");

          model.param().set("TAs_1", "295", "Austenite start temperature");

          model.param().set("TAf_1", "315", "Austenite finish temperature");



          model.param().set("C_M", "4e6", "Martensite slope");

          model.param().set("C_A", "4e6", "Austenite slope");

          model.param().set("etr_max", "0.04", "max transformation strain");

          model.param().set("xi_0", "0.0", "Initial martensite fraction");





// Displacement Function Piecewise



          model.func().create("pw1", "Piecewise");
@@ -86,20 +106,20 @@ public class Lag_sme{ 
 //Material block Austenite

          model.component("comp1").material().create("mat1", "Common");

          model.component("comp1").material("mat1").label("Austenite");

          model.component("comp1").material("mat1").propertyGroup("def").set("poissonsratio", "0.33");

          model.component("comp1").material("mat1").propertyGroup("def").set("density", "6500");

          model.component("comp1").material("mat1").propertyGroup("def").set("youngsmodulus", "41.2e9");

          model.component("comp1").material("mat1").propertyGroup("def").set("heatcapacity", "400");

          model.component("comp1").material("mat1").propertyGroup("def").set("thermalconductivity", "18.00");

          model.component("comp1").material("mat1").propertyGroup("def").set("poissonsratio", "p_ratio");

          model.component("comp1").material("mat1").propertyGroup("def").set("density", "rho_");

          model.component("comp1").material("mat1").propertyGroup("def").set("youngsmodulus", "Ymod_A");

          model.component("comp1").material("mat1").propertyGroup("def").set("heatcapacity", "Hcap");

          model.component("comp1").material("mat1").propertyGroup("def").set("thermalconductivity", "kappa_A");



//Material block Martensite

          model.component("comp1").material().create("mat2", "Common");

          model.component("comp1").material("mat2").label("Martensite");

          model.component("comp1").material("mat2").propertyGroup("def").set("poissonsratio", "0.33");

          model.component("comp1").material("mat2").propertyGroup("def").set("density", "6500");

          model.component("comp1").material("mat2").propertyGroup("def").set("youngsmodulus", "32.9e9");

          model.component("comp1").material("mat2").propertyGroup("def").set("heatcapacity", "400");

          model.component("comp1").material("mat2").propertyGroup("def").set("thermalconductivity","8.6");

          model.component("comp1").material("mat2").propertyGroup("def").set("poissonsratio", "p_ratio");

          model.component("comp1").material("mat2").propertyGroup("def").set("density", "rho_");

          model.component("comp1").material("mat2").propertyGroup("def").set("youngsmodulus", "Ymod_M");

          model.component("comp1").material("mat2").propertyGroup("def").set("heatcapacity", "Hcap");

          model.component("comp1").material("mat2").propertyGroup("def").set("thermalconductivity","kappa_M");

          model.component("comp1").material("mat2").selection().all();



// Heat transfer
@@ -115,7 +135,13 @@ public class Lag_sme{ 
          // sma1

          model.component("comp1").physics("ht").create("sma1", "ShapeMemoryAlloy", 2);

          model.component("comp1").physics("ht").feature("sma1").selection().set(1);

          

          model.component("comp1").physics("ht").feature("sma1").set("EffectiveConductivity", "PowerLaw");

          model.component("comp1").physics("ht").feature("sma1").set("MartensiteMaterial", "mat2");

          model.component("comp1").physics("ht").feature("sma1").set("AusteniteMaterial", "mat1");

          model.component("comp1").physics("ht").feature("sma1").set("xi", "xi_0");

          model.component("comp1").physics("ht").feature("sma1").set("rho_mat", "userdef");

          model.component("comp1").physics("ht").feature("sma1").set("rho", "rho_");          





          // temp1
@@ -132,7 +158,7 @@ public class Lag_sme{ 
          model.component("comp1").physics("ht").feature("hf1").selection().set(2, 3);

          model.component("comp1").physics("ht").feature("hf1").set("materialType", "solid");

          model.component("comp1").physics("ht").feature("hf1").set("HeatFluxType", "ConvectiveHeatFlux");

          model.component("comp1").physics("ht").feature("hf1").set("h", 40);

          model.component("comp1").physics("ht").feature("hf1").set("h", "hflux");

          model.component("comp1").physics("ht").feature("hf1").set("Text", "T_ref");



          //ophf1
@@ -140,11 +166,11 @@ public class Lag_sme{ 
          model.component("comp1").physics("ht").feature("ophf1").selection().set(1);



          model.component("comp1").physics("ht").feature("ophf1").set("UpsideHeatFluxType", "ConvectiveHeatFlux");

          model.component("comp1").physics("ht").feature("ophf1").set("h_u", 40);

          model.component("comp1").physics("ht").feature("ophf1").set("h_u", "hflux");

          model.component("comp1").physics("ht").feature("ophf1").set("Text_u", "T_ref");



          model.component("comp1").physics("ht").feature("ophf1").set("DownsideHeatFluxType", "ConvectiveHeatFlux");

          model.component("comp1").physics("ht").feature("ophf1").set("h_d", 40);

          model.component("comp1").physics("ht").feature("ophf1").set("h_d", "hflux");

          model.component("comp1").physics("ht").feature("ophf1").set("Text_d", "T_ref");



          // hs1
@@ -200,23 +226,21 @@ public class Lag_sme{ 
          model.component("comp1").physics("solid").feature("sma1").set("AusteniteMaterial", "mat1");

          model.component("comp1").physics("solid").feature("sma1").set("MartensiteMaterial", "mat2");



          model.component("comp1").physics("solid").feature("sma1").set("TMs", "291+rn(x,y)");

          model.component("comp1").physics("solid").feature("sma1").set("TMf", "271+rn(x,y)");

          model.component("comp1").physics("solid").feature("sma1").set("CM", "4e6");

          model.component("comp1").physics("solid").feature("sma1").set("TMs", "TMs_1+rn(x,y)");

          model.component("comp1").physics("solid").feature("sma1").set("TMf", "TMf_1+rn(x,y)");

          model.component("comp1").physics("solid").feature("sma1").set("CM", "C_M");



          model.component("comp1").physics("solid").feature("sma1").set("TAs", "295+rn(x,y)");

          model.component("comp1").physics("solid").feature("sma1").set("TAf", "315+rn(x,y)");

          model.component("comp1").physics("solid").feature("sma1").set("CA", "4e6");

          model.component("comp1").physics("solid").feature("sma1").set("TAs", "TAs_1+rn(x,y)");

          model.component("comp1").physics("solid").feature("sma1").set("TAf", "TAf_1+rn(x,y)");

          model.component("comp1").physics("solid").feature("sma1").set("CA", "C_A");



          model.component("comp1").physics("solid").feature("sma1").set("etrmax", 0.04);

          model.component("comp1").physics("solid").feature("sma1").set("xi0", 1);

          model.component("comp1").physics("solid").feature("sma1").set("etrmax", "etr_max");

          model.component("comp1").physics("solid").feature("sma1").set("xi0", "xi_0");



          model.component("comp1").physics("solid").feature("sma1").set("TransformationFunction", "Quadratic");

          model.component("comp1").physics("solid").feature("sma1").set("TransformationDirectionType", "UserDefined");

          model.component("comp1").physics("solid").feature("sma1").set("trdir", "if(t<t3,1,-1)");



          

          

          model.component("comp1").physics("solid").feature("sma1").set("CalculateDissipatedEnergy", true);



          // Initial values
Impressum   Privacy Policy   WW8 Homepage   Explore   Help