Commit a0a5013d authored by Gilles MARCKMANN's avatar Gilles MARCKMANN

Merge branch 'master' of git.gem.ec-nantes.fr:TTK/TLSDuctile

parents 3b464648 5da71693
Pipeline #137 skipped
This diff is collapsed.
......@@ -250,6 +250,11 @@ void AsymElastoDam::sensitivityTo(const std::string& phys_token, double& sensiti
Observer::tell("load_factor", load_factor);
sensitivity=tmp.contract(eps_ref)*load_factor;
}
else if (phys_token == "delta_dissipated_energy") {
double d=curr->scalar("d");
double delta_d=d-old->scalar("d");
sensitivity=Y_c*hardening_function->getVal(d)*delta_d;
}
else if(phys_token == "Y") {
sensitivityTo("a", sensitivity);
}
......
......@@ -4,7 +4,7 @@ disp_space_dim = V2Dxy
disp_space_order = 1
disp_bc_integ_order = 1
mean_field_space_type = PolyLagrange
damage_shape = Exp2
damage_shape = Poly2Revert
lc = 6
delta_phi_max_ratio = 0.99
do_forced_delocalization = 0
......@@ -17,5 +17,5 @@ curve_plot_file_ext = none
export_manager = {
disp 10
damage 10 }
export_sensors_label = { disp_max 1 force_y 1 load_factor 1 }
export_sensors_label = { load_factor 1 disp_y 1 }
export_sensors_point = { }
......@@ -8,4 +8,4 @@ ZONE 302 = { MAT_CLASS = AsymElastoDam MAT_PARAM = mate_far.mat }
BC_POINT 101 = { DISPLACEMENT_X FIX = 0.
DISPLACEMENT_Y FIX = 0. }
BC_POINT 102 = { DISPLACEMENT_Y FIX = 0. }
BC_LINE 201 = { TRACTION_Y FIX = -1.e3 }
BC_LINE 201 = { TRACTION_Y FIX_AND_MEASURE = -1.e3 }
// can be changed
r=5.;
R=1.3*r;
rho=r/10.;
costheta=0.70710678118;
sintheta=0.70710678118;
Nn=30;
a=1.;
lc=1.;
lcmin=lc/Nn;
lcmax=10.*lcmin;
distmin=1.5*lc;
distmax=3.*lc;
Point(1) = {0, 0, 0, a};
Point(2) = {-r*costheta, r*sintheta, 0, a};
Point(3) = {r, 0, 0, a};
Point(4) = {-r*costheta, -r*sintheta, 0, a};
Point(5) = {-R*costheta, R*sintheta, 0, a};
Point(6) = {R, 0, 0, a};
Point(7) = {-R*costheta, -R*sintheta, 0, a};
Point(8) = {r-distmax, 0, 0, a};
Point(9) = {-rho/sintheta, 0, 0, a};
Point(10) = {-rho*costheta*costheta/sintheta, rho*costheta, 0, a};
Point(11) = {-rho*costheta*costheta/sintheta, -rho*costheta, 0, a};
Point(12) = {-rho/sintheta+rho, 0, 0, a};
Line(1) = {2, 10};
Line(2) = {2, 5};
Line(3) = {4, 11};
Line(4) = {4, 7};
Circle(5) = {2, 1, 3};
Circle(6) = {3, 1, 4};
Circle(7) = {5, 1, 6};
Circle(8) = {6, 1, 7};
Circle(9) = {10, 9, 11};
Line(10) = {12, 8};
Line Loop(10) = {6, 3, -9, -1, 5};
Plane Surface(11) = {10};
Line Loop(12) = {7, 8, -4, -6, -5, 2};
Plane Surface(13) = {12};
Physical Point(1) = {1};
Physical Line(11) = {7, 8};
Physical Surface(101) = {11};
Physical Surface(102) = {13};
Field[1]=Attractor;
Field[1].NNodesByEdge=Nn;
Field[1].EdgesList={10};
Field[2]=Threshold;
Field[2].IField=1;
Field[2].LcMin=lcmin;
Field[2].LcMax=lcmax;
Field[2].DistMin=distmin;
Field[2].DistMax=distmax;
Background Field=2;
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