Mejora en título, incluye sección

Author: ppizarror

mcurv/section/SectionAnalysis.m

@@ -339,14 +339,16 @@ function plot_lastIter(obj)
             movegui(plt, 'center');
             set(gcf, 'name', 'Momento curvatura');
             hold on;
+            secName = obj.lastsole0p{8}.getName();
             niter = obj.lastsole0p{9};
 
-            plot(1:length(niter), niter, 'k-', 'Linewidth', 1.5);
+            plot(1:length(niter), niter, '-', 'Linewidth', 1.5);
             grid on;
             grid minor;
 
             xlabel('Numero de paso');
             ylabel('Numero de iteraciones');
+            title({'Variacion numero iteraciones', secName});
 
         end % plot_lastIter function
 
@@ -691,6 +693,7 @@ function plot_e0M_mcurv(obj, phi, mxInt, myInt, r, curvAxis, secName, angle) %#o
             mmin = min(ylm);
 
             % Calcula las interpolaciones
+            kphi = 1; % Numero de puntos agregados
             if (strcmp(r.m, 'x') || strcmp(r.m, 'y') || strcmp(r.m, 'T'))
                 for i = 1:length(r.vecphi)
 
@@ -712,14 +715,15 @@ function plot_e0M_mcurv(obj, phi, mxInt, myInt, r, curvAxis, secName, angle) %#o
                         fprintf('\tphi %e: Momento %f %s\n', phiobj, mi, r.unitloadM);
                         plot([phiobj, phiobj], [mmin, mi], '--', ...
                             'Color', r.vecphiColor{i}, 'LineWidth', r.vecphiLw, ...
-                            'DisplayName', sprintf('\\phi=%.3e', phiobj));
+                            'DisplayName', sprintf('\\phi_%d=%.3e', kphi, phiobj));
                         pl = plot([min(phi), phiobj], [mi, mi], '--', ...
                             'Color', r.vecphiColor{i}, 'LineWidth', r.vecphiLw);
                         set(get(get(pl, 'Annotation'), 'LegendInformation'), 'IconDisplayStyle', 'off');
                         pl = plot(phiobj, mi, '.', ...
                             'Color', r.vecphiColor{i}, 'LineWidth', r.vecphiLw, ...
                             'MarkerSize', r.vecphiSize);
                         set(get(get(pl, 'Annotation'), 'LegendInformation'), 'IconDisplayStyle', 'off');
+                        kphi = kphi + 1;
                     end
 
                 end
@@ -795,13 +799,13 @@ function plot_e0M_pcurv(obj, phi, pInt, r, curvAxis, secName, angle) %#ok<INUSL>
             plot(phi, pInt, '-', 'LineWidth', r.linewidth);
             grid on;
             grid minor;
-            ylabel(sprintf('Carga axial P (%s)', r.unitloadP));
+            ylabel(sprintf('Carga axial interna P_{int} (%s)', r.unitloadP));
 
             if ~strcmp(curvAxis, 'a')
-                title({sprintf('Carga axial vs curvatura \\phi_%s', curvAxis), secName});
+                title({sprintf('Carga axial interna vs curvatura \\phi_%s', curvAxis), secName});
                 xlabel(sprintf('Curvatura \\phi_%s (%s)', curvAxis, r.unitlength));
             else
-                title({sprintf('Carga axial vs curvatura \\phi - Angulo %.1f', angle), secName});
+                title({sprintf('Carga axial interna vs curvatura \\phi - Angulo %.1f', angle), secName});
                 xlabel(sprintf('Curvatura \\phi (%s)', r.unitlength));
             end
             if r.limPos

test/section/wallT.m

@@ -27,7 +27,7 @@
 % steel.plot('emin', -2.5e-1, 'emax', 3.1e-1);
 
 % Agrega los elementos
-caseNum = '1';
+caseNum = '1'; % Numero caso analisis
 h = 5000;
 b = 4000;
 bw = 300;
@@ -148,17 +148,6 @@
 analysis = SectionAnalysis('Analisis', 1000, 0.01, 'showprogress', true);
 analysis.calc_e0M_angle(wallt, p, phi, curvang, 'ppos', [0, 0]);
 
-% Grafica, N*mm -> kN*m
-if showSap
-    analysis.plot_e0M('plot', 'mphi', 'factorM', 1e-6, 'm', 'T', ...
-        'sapfile', sprintf('test/section/mcurv-sap2000/wallT%s_%d.txt', caseNum, curvang), ...
-        'sapcolumnPhi', 10, 'sapcolumnM', 11, 'sapfactorM', 1e-6, ...
-        'sapdiff', true, 'vecphi', [phiDef3, phiDef8] .* ~showSap, 'vecphiColor', {'r', 'k'});
-else
-    analysis.plot_e0M('plot', 'mphi', 'factorM', 1e-6, 'm', 'T', ...
-        'vecphi', [phiDef3, phiDef8], 'vecphiColor', {'r', 'k'});
-end
-
 % Grafica tension y deformacion para deformacion de 0.003 y 0.008 si es que
 % aplica
 if showST && ~showSap
@@ -170,4 +159,16 @@
     end
 end
 
-% analysis.plot_lastIter();
+% Grafica momento curvatura, N*mm -> kN*m
+if showSap
+    analysis.plot_e0M('plot', 'mphi', 'factorM', 1e-6, 'm', 'T', ...
+        'sapfile', sprintf('test/section/mcurv-sap2000/wallT%s_%d.txt', caseNum, curvang), ...
+        'sapcolumnPhi', 10, 'sapcolumnM', 11, 'sapfactorM', 1e-6, ...
+        'sapdiff', true, 'vecphi', [phiDef3, phiDef8] .* ~showSap, 'vecphiColor', {'r', 'k'});
+else
+    analysis.plot_e0M('plot', 'mphi', 'factorM', 1e-6, 'm', 'T', ...
+        'vecphi', [phiDef3, phiDef8], 'vecphiColor', {'r', 'k'});
+    analysis.plot_e0M('plot', 'pphi');
+end
+
+analysis.plot_lastIter();