Diferencia entre revisiones de «Guía de optimización en ingeniería»
(→Escribir la función coste y las restricciones en términos de los parámetros) |
(→Planteamiento matemático) |
||
| Línea 13: | Línea 13: | ||
Vamos a plantear este problema matemáticamente. Los pasos a seguir son: | Vamos a plantear este problema matemáticamente. Los pasos a seguir son: | ||
# Identificar los grados de libertad y asociarle uno o varios parámetros para describir todas las posibles geometrías. Llamaremos a los parámetros <math> (a_1,a_2,...,a_n)</math> | # Identificar los grados de libertad y asociarle uno o varios parámetros para describir todas las posibles geometrías. Llamaremos a los parámetros <math> (a_1,a_2,...,a_n)</math> | ||
| − | # Escribir la función coste, los | + | # Escribir en función de los parámetros anteriores: la función coste <math> f_{coste}(a_1,a_2,...,a_n)</math>, los intervalos <math> [L_i,M_i] </math> en los que se pueden mover los parámetros <math> a_i, \; i=1,...,n,</math> y las restricciones en términos de los parámetros <math> g_{rest}(a_1,a_2,...,a_n)</math>. |
| + | |||
| + | Una vez hecho esto el planteamiento matemático es: Buscamos los valores <math> (a_1,a_2,...,a_n)</math> que minimicen el coste | ||
| + | |||
| + | <math> f_{coste} (a_1,a_2,...,a_n) </math> | ||
| + | |||
| + | cuando los parámetros toman valores en los intervalos | ||
| + | |||
| + | <math> L_i \leq a_i \leq M_i, \quad i=1,2,...,n, </math> | ||
| + | |||
| + | y además verifican las restricciones | ||
<math> \begin{array}{l} | <math> \begin{array}{l} | ||
| − | |||
| − | |||
g_{rest1} (a_1,a_2,...,a_n)\leq 0\\ | g_{rest1} (a_1,a_2,...,a_n)\leq 0\\ | ||
g_{rest2} (a_1,a_2,...,a_n)\leq 0\\ | g_{rest2} (a_1,a_2,...,a_n)\leq 0\\ | ||
| Línea 25: | Línea 33: | ||
</math> | </math> | ||
| − | |||
Vamos a seguir estos dos pasos en nuestro ejemplo: | Vamos a seguir estos dos pasos en nuestro ejemplo: | ||
=== Identificar los grados de libertad y asociarles parámetros === | === Identificar los grados de libertad y asociarles parámetros === | ||
| Línea 71: | Línea 78: | ||
<math> | <math> | ||
| − | + | g_{rest}=\frac{a^2+a+1}{3(a+1)}-a \leq 0 | |
</math> | </math> | ||
Revisión del 17:20 1 feb 2017
Este artículo contiene una guía básica para plantear y resolver un problema de optimización en ingeniería. Para fijar ideas nos centraremos en un ejemplo sencillo e iremos explicando el proceso sobre ese ejemplo.
Contenido
1 Ejemplo: definición del problema
Consideramos un trapecio construido de un material homogéneo del que podemos elegir la inclinación del lado derecho (ver la figura). Buscamos la inclinación de ese lado derecho que proporcione menor área, pero manteniendo la componente horizontal del centro de gravedad dentro de la base para darle estabilidad.
Es importante identificar en este punto dos cosas:
- Coste que queremos optimizar
- Restricciones que tenemos en el problema
En este caso sencillo, el coste consiste en el área total del trapecio mientras que la restricción es que el centro de gravedad caiga en la base.
2 Planteamiento matemático
Vamos a plantear este problema matemáticamente. Los pasos a seguir son:
- Identificar los grados de libertad y asociarle uno o varios parámetros para describir todas las posibles geometrías. Llamaremos a los parámetros [math] (a_1,a_2,...,a_n)[/math]
- Escribir en función de los parámetros anteriores: la función coste [math] f_{coste}(a_1,a_2,...,a_n)[/math], los intervalos [math] [L_i,M_i] [/math] en los que se pueden mover los parámetros [math] a_i, \; i=1,...,n,[/math] y las restricciones en términos de los parámetros [math] g_{rest}(a_1,a_2,...,a_n)[/math].
Una vez hecho esto el planteamiento matemático es: Buscamos los valores [math] (a_1,a_2,...,a_n)[/math] que minimicen el coste
[math] f_{coste} (a_1,a_2,...,a_n) [/math]
cuando los parámetros toman valores en los intervalos
[math] L_i \leq a_i \leq M_i, \quad i=1,2,...,n, [/math]
y además verifican las restricciones
[math] \begin{array}{l} g_{rest1} (a_1,a_2,...,a_n)\leq 0\\ g_{rest2} (a_1,a_2,...,a_n)\leq 0\\ ...\\ g_{restk} (a_1,a_2,...,a_n)\leq 0\\ \end{array} [/math]
Vamos a seguir estos dos pasos en nuestro ejemplo:
2.1 Identificar los grados de libertad y asociarles parámetros
En el ejemplo que tratamos sólo hay un grado de libertad que es la pendiente del lado derecho del trapecio. En general existen varios parámetros que podemos elegir para representar este grado de libertad: la distancia a, la propia pendiente del lado en rojo, el ángulo que forma ese lado con la horizontal, etc. Conviene elegir el que nos haga más fáciles los cálculos.
En este caso elegimos la longitud de la base a. En situaciones más complicadas habrá más parámetros para describir toda la variabilidad del problema. Por ejemplo, podemos incluir la altura h o la longitud del lado superior b. Sin embargo, para fijar ideas, en lo sucesivo estos dos parámetros los tomaremos constantes e iguales a 1.
[math] h=b=1. [/math]
2.2 Escribir la función coste y las restricciones en términos de los parámetros
La función coste es el área del trapecio, que viene dada por
[math] f_{coste}=h\frac{b+a}{2}=\frac{1+a}{2}. [/math]
El parámetro [math]a[/math] supondremos varía en el intervalo [math][0,1][/math] puesto que otros valores bien no tienen sentido (por ejemplo [math]a\lt0[/math]) o bien podemos descartarlos por ser claramente poco eficientes (por ejemplo [math]a\gt1[/math]).
En cuanto a las restricciones, en este caso sólo tenemos una: que la coordenada x del centro de gravedad caiga en la base. Si llamamos [math]S[/math] a la superficie que ocupa el trapecio, la coordenada [math]x_G[/math] del centro de gravedad viene dada por
[math] x_{G}=\frac{\iint_S x dS}{\iint_S dS}=\frac{\int_0^a\int_0^1 x \; dydx+\int_a^1\int_{\frac{x-a}{1-a}}^1x\;dydx}{\frac{1+a}{2}} =\frac{a^2+a+1}{3(a+1)}. [/math]
Por lo tanto, la restricción se escribe [math] x_{G}=\frac{a^2+a+1}{3(a+1)}\lta [/math]
El planteamiento matemático del problema es entonces: Buscamos [math]a[/math] que minimice el coste
[math] f_{coste}=\frac{1+a}{2} \frac{a^2+a+1}{3(a+1)}\leq a [/math]
en el intervalo
[math] 0\leq a \leq 1, [/math]
con la restricción
[math] g_{rest}=\frac{a^2+a+1}{3(a+1)}-a \leq 0 [/math]
