Commit edb104d4 authored by Guido Giuntoli's avatar Guido Giuntoli

Small fix and configuring benchmark-mic-4

parent 484b4861
...@@ -310,7 +310,8 @@ void micropp<tdim>::calc_ctan_lin_mix_rule_Chamis(double ctan[nvoi * nvoi]) ...@@ -310,7 +310,8 @@ void micropp<tdim>::calc_ctan_lin_mix_rule_Chamis(double ctan[nvoi * nvoi])
const double E22 = Em / (1 - sqrt(Vf) * (1 - Em / Ef)); const double E22 = Em / (1 - sqrt(Vf) * (1 - Em / Ef));
const double nu12 = Vf * nu_f + Vm * nu_m; const double nu12 = Vf * nu_f + Vm * nu_m;
const double G12 = Gm / (1 - sqrt(Vf) * (1 - Gm / Gf)); const double G12 = Gm / (1 - sqrt(Vf) * (1 - Gm / Gf));
const double nu23 = E22 / (2 * G12); //const double nu23 = E22 / (2 * G12);
const double nu23 = nu12;
const double S[3][3] = { const double S[3][3] = {
{ 1 / E11, - nu12 / E11, - nu12 / E11 }, { 1 / E11, - nu12 / E11, - nu12 / E11 },
......
...@@ -81,8 +81,8 @@ int main(int argc, char **argv) ...@@ -81,8 +81,8 @@ int main(int argc, char **argv)
mic_params.geo_params[0] = 0.1; mic_params.geo_params[0] = 0.1;
mic_params.geo_params[1] = 0.02; mic_params.geo_params[1] = 0.02;
mic_params.geo_params[2] = 0.01; mic_params.geo_params[2] = 0.01;
material_set(&mic_params.materials[0], 0, 1.0e7, 0.3, 0.0, 0.0, 0.0); material_set(&mic_params.materials[0], 0, 3.0e7, 0.25, 0.0, 0.0, 0.0);
material_set(&mic_params.materials[1], 0, 3.0e7, 0.3, 0.0, 0.0, 0.0); material_set(&mic_params.materials[1], 0, 3.0e8, 0.25, 0.0, 0.0, 0.0);
material_set(&mic_params.materials[2], 0, 3.0e7, 0.3, 0.0, 0.0, 0.0); material_set(&mic_params.materials[2], 0, 3.0e7, 0.3, 0.0, 0.0, 0.0);
mic_params.lin_stress = true; mic_params.lin_stress = true;
...@@ -152,17 +152,22 @@ int main(int argc, char **argv) ...@@ -152,17 +152,22 @@ int main(int argc, char **argv)
} }
cout << endl; cout << endl;
micro.update_vars(); file << eps[dir] << "\t";
for (int i = 0; i < 6; ++i) {
file << eps[dir] << "\t" file << sig[i] << "\t";
<< sig[dir] << "\t" << endl; }
file << sqrt(pow(sig[0], 2) + pow(sig[1], 2) + pow(sig[2], 2) + \
pow(sig[3], 2) + pow(sig[4], 2) + pow(sig[5], 2)) << "\t"
<< endl;
if (print) { if (print) {
char filename[128]; char filename[128];
snprintf(filename, 128, "micropp_%d", t); snprintf(filename, 128, "micropp_%d", t);
micro.output (0, filename); micro.output(0, filename);
} }
micro.update_vars();
time += dt; time += dt;
} }
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment