Implementación numérica de una ecuación diferencial de movimiento en un grado de libertad con componente estocástica

Abstract

En dinámica, mediante la ecuación diferencial ordinaria de movimiento, es posible determinar la posición en el tiempo de una masa que se desplaza debido a que es perturbada por alguna acción determinística. En este trabajo se propuso aplicar a la masa una perturbación no determinística de origen sísmico en un grado de libertad vertical y dentro del rango lineal. La pregunta de investigación fue: ¿será´ posible migrar la ecuación diferencial ordinaria (EDO) de movimiento hacia una ecuación diferencial estocástica (EDE) de movimiento? Bajo ese marco, se estudiaron los fundamentos de la teoría de la probabilidad y los procesos estocásticos. Utilizando estas ramas de las matemáticas aplicadas se logró obtener una EDE de movimiento. Se estudió también la aproximación de Euler-Maruyama la cual se implementó, luego de verificar su estabilidad estocástica y numérica, para obtener una solución de la EDE de movimiento encontrada. Los resultados obtenidos permitieron confirmar que el uso de una versión no determinística genera resultados satisfactorios. Se recomienda efectuar análisis similares con otras variables, por ejemplo, en sistemas con un grado de libertad diferente, con más de un grado de libertad y/o considerando un comportamiento no lineal.In dynamics, using the ordinary differential equation of motion, it is possible to determine the position in time of a moving mass because it is disturbed by some deterministic action. In this work it was proposed to apply to the mass a non-deterministic disturbance of seismic origin in a vertical degree of freedom and within the linear range. The research question was: Will it be possible to migrate the ordinary differential equation (ODE) of motion to a stochastic differential equation (SDE) of motion? Under this framework, the foundations of probability theory and stochastic processes were studied. Using these branches of applied mathematics, an SDE of movement was obtained. The Euler-Maruyama approximation was also studied, which was implemented, after verifying its stochastic and numerical stability, to obtain a solution of the EDE of the movement found. The results obtained confirmed that the use of a non-deterministic version generates satisfactory results. It is recommended to carry out similar analyzes with other variables, for instance, in systems with a different degree of freedom, with more than one degree of freedom and / or considering non-linear behavior.