Estabilidad de una presa de gravedad - Historial de revisiones

Carlos Castro: /* Programa en MATLAB */

2017-03-08T12:26:12Z

‎Programa en MATLAB

Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/mat/public_html/w/includes/diff/DairikiDiff.php on line 434

Carlos Castro: /* Programa en MATLAB */

2017-03-08T11:08:54Z

‎Programa en MATLAB

Carlos Castro en 09:16 8 mar 2017

2017-03-08T09:16:56Z

Carlos Castro: /* Programa en MATLAB */

2017-03-07T22:48:42Z

‎Programa en MATLAB

Carlos Castro: /* Programa en MATLAB */

2017-03-07T22:41:09Z

‎Programa en MATLAB

Carlos Castro: /* Observaciones: */

2017-03-05T16:08:49Z

‎Observaciones:

Carlos Castro: /* Observaciones: */

2017-03-03T22:06:07Z

‎Observaciones:

Carlos Castro: /* Observaciones: */

2017-03-03T22:02:48Z

‎Observaciones:

Carlos Castro: /* Observaciones: */

2017-03-03T21:41:04Z

‎Observaciones:

Carlos Castro: /* Observaciones: */

2017-03-03T21:38:54Z

‎Observaciones:

@@ Línea 69: / Línea 69: @@
 % la estabilidad requiere g_rest3(n) < 0
 % De hecho, si S_A<S_Amin, y asumimos que la distribución de tensiones crece linealmente cuando
-% nos movemos hacia el pie de aguas abajo, el momento de las fuerzas normales (para un mismo peso y subpresión)
+% nos movemos hacia el pie de aguas abajo, el momento de las fuerzas normales
-% sería menor que el del caso S_A=S_Amin y la restricción no se verificaría
+% (para un mismo peso y subpresión) sería menor que el del caso S_A=S_Amin
 % %% Planteamiento matemático del problema de optimización:

@@ Línea 59: / Línea 59: @@
 % + momento subpresion S_2
 % + momento subpresion S_1
-% + momento de la normal N_2
+% + momento de la normal N_2 (cuando suponemos S_A=S_Amin)
-% + momento de la normal N_1 (basado en la cota superior de S_B calculada anteriormente)
+% + momento de la normal N_1 (cuando suponemos S_A=S_Amin y por tanto S_B=Cota_S_B)
 g_rest3 = @(n) -P(n)*base(n)*2/3 ...
      + E*h_NMN/3 ...
@@ Línea 68: / Línea 68: @@
      + 0.5*base(n)*(Cota_S_B(n)-S_Amin*10)*(base(n)/3);
 % la estabilidad requiere g_rest3(n) < 0
 % %% Planteamiento matemático del problema de optimización:

@@ Línea 45: / Línea 45: @@
 % %% 2. Compresión en el pie de aguas abajo menor que el máximo admitido por el cimiento
 % N=P-S=1/2*(S_A+S_B)*base(n)
-% Suponiendo que al menos tenemos la compresión mínima, es decir S_A >= S_Amin,
+% Suponiendo que al menos tenemos la compresión mínima en el pie de aguas arriba, es decir S_A >= S_Amin,
-% S_B<=(P-S)/(base(n)/2)-S_Amin que en este caso no depende de n pero puede
+% podemos calcular una cota superior para la tensión aguas abajo
 % depender en general, así que lo definimos como función
-S_B = @(n) (P(n)-S(n))/(0.5*base(n))-S_Amin*10;
+Cota_S_B = @(n) (P(n)-S(n))/(0.5*base(n))-S_Amin*10;
-g_rest2 = @(n) S_B(n)-S_Bmax*10;
+g_rest2 = @(n) Cota_S_B(n)-S_Bmax*10;
 % la estabilidad requiere g_rest2(n) < 0
@@ Línea 59: / Línea 60: @@
 % + momento subpresion S_1
 % + momento de la normal N_2
-% + momento de la normal N_1
+% + momento de la normal N_1 (basado en la cota superior de S_B calculada anteriormente)
 g_rest3 = @(n) -P(n)*base(n)*2/3 ...
      + E*h_NMN/3 ...
@@ Línea 65: / Línea 66: @@
      + 0.5*base(n)*(h_NMN-h_CE)*2/3*base(n)...
      + S_Amin*10*base(n)*(base(n)*0.5) ...
-     + 0.5*base(n)*(S_B(n)-S_Amin*10)*(base(n)/3);
+     + 0.5*base(n)*(Cota_S_B(n)-S_Amin*10)*(base(n)/3);
 % la estabilidad requiere g_rest3(n) < 0

@@ Línea 45: / Línea 45: @@
 % %% 2. Compresión en el pie de aguas abajo menor que el máximo admitido por el cimiento
 % N=P-S=1/2*(S_A+S_B)*base(n)
-% Suponiendo la compresión mínima S_Amin dada,
+% Suponiendo que al menos tenemos la compresión mínima, es decir S_A >= S_Amin,
-% S_B=(P-S)/(base(n)/2)-S_Amin que en este caso no depende de n pero puede
+% S_B<=(P-S)/(base(n)/2)-S_Amin que en este caso no depende de n pero puede
 % depender en general, así que lo definimos como función
 S_B = @(n) (P(n)-S(n))/(0.5*base(n))-S_Amin*10;

@@ Línea 64: / Línea 64: @@
      + h_CE*base(n)*(base(n)/2)*mu_a ...
      + 0.5*base(n)*(h_NMN-h_CE)*2/3*base(n)...
-     + S_Amin*base(n)*(base(n)*0.5) ...
+     + S_Amin*10*base(n)*(base(n)*0.5) ...
-     + 0.5*base(n)*(S_B(n)-S_Amin)*(base(n)/3);
+     + 0.5*base(n)*(S_B(n)-S_Amin*10)*(base(n)/3);
 % la estabilidad requiere g_rest3(n) < 0

@@ Línea 176: / Línea 176: @@
 {{matlab|codigo=
 % Definimos una matriz con las coordenadas del polígono
 c=[-1 0 -1 1 1;-1 0 1 1 -1];

@@ Línea 149: / Línea 149: @@
 % Función para calcular el área de un polígono
 function S=ar_pol(C)
-% Se introduce una matriz C con los vértices del polígono
+% Se introduce una matriz C con las coordenadas de los vértices del polígono
-% ordenados en sentido antihorario.
+% por columnas y ordenados en sentido antihorario.
 % Se devuelve el área S
 % Aplicamos la fórmula que nos envió Rafael Morán
@@ Línea 161: / Línea 161: @@
 nv=length(C);
 % Ahora definimos los vectores con las coordenadas x e y de los vértices
-x=C(1,:);
+x=C(1,:);    % la primera fila contiene las coordenadas x de los vértices
-y=C(2,:);
+y=C(2,:);    % la segunda fila contiene las coordenadas y de los vértices
 % Calculamos S pero separamos el primer y último término de la suma
 % para poder tener en cuenta la regla

@@ Línea 153: / Línea 153: @@
 % Se devuelve el área S
 % Aplicamos la fórmula que nos envió Rafael Morán
-% S = 1/2 abs(sum yi (x(i-1)-x(i+1)))
+% S = 1/2 abs(sum yi (x_(i-1)-x_(i+1)))
 % teniendo en cuenta la siguiente regla: (x_0,y_0)=(x_n,y_n) y
 % (x_(n+1),y_(n+1))=(x_1,y_1)