Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 2 vistas

    punto 1

    Cargado por

    santiagozarategarzon
    Punto numero 1

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd

    punto 1

    Cargado por

    santiagozarategarzon
    2 vistas4 páginas
    Punto numero 1

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Punto numero 1

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    2 vistas4 páginas

    punto 1

    Cargado por

    santiagozarategarzon
    Punto numero 1

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    de 4

    In[116]:=

    (* Ejercicio 1: Sistema de resortes acoplados - Procedimiento más detallado *)

    (* Definimos las constantes y las masas del problema con sus unidades explícitas *)

    k1 = 4; (* N/m *)

    k2 = 6; (* N/m *)

    m1 = 1; (* kg *)

    m2 = 1; (* kg *)

    f[t_] := 30 Cos[t]; (* N *)

    (* Definimos las variables como funciones del tiempo *)

    Clearx, y;

    (* Ecuaciones de movimiento basadas en la Segunda Ley de Newton *)

    (* Fuerza neta sobre m1 *)

    eq1 = m1 x ''[t] ⩵ - k1 x[t] + k2 y[t] - x[t];


    Print
    "Ecuación de movimiento para la masa m1 derivada de la Segunda Ley de Newton:";
    eq1 // TraditionalForm

    (* Fuerza neta sobre m2 *)

    eq2 = m2 y ''[t] ⩵ - k2 y[t] - x[t] + f[t];


    Print
    "Ecuación de movimiento para la masa m2 derivada de la Segunda Ley de Newton:";
    eq2 // TraditionalForm

    (* Sustituimos los valores conocidos para simplificar las ecuaciones *)

    eq1Simplified = eq1 /. k1 → 4, k2 → 6, m1 → 1;


    eq2Simplified = eq2 /. k2 → 6, m2 → 1;

    Print"Ecuaciones simplificadas:";
    eq1Simplified // TraditionalForm
    eq2Simplified // TraditionalForm

    (* Reorganizamos las ecuaciones para obtener


    un sistema de ecuaciones diferenciales acopladas *)

    Print"Reorganización del sistema de ecuaciones:";


    eq1Reorganized = x ''[t] ⩵ - 10 x[t] + 6 y[t];
    eq2Reorganized = y ''[t] ⩵ 6 x[t] - 6 y[t] + 30 Cos[t];
    eq1Reorganized // TraditionalForm
    eq2Reorganized // TraditionalForm

    (* Aplicación de la Transformada de Laplace a ambas ecuaciones *)


    2

    Print"Aplicando la Transformada de Laplace a las ecuaciones:";


    Clears, X, Y;

    laplaceEq1 = LaplaceTransformeq1Reorganized, t, s /.
    LaplaceTransformx[t], t, s → X[s], LaplaceTransformy[t], t, s → Y[s],
    x[0] → 0, y[0] → 0, Derivative[1][x][0] → 0, Derivative[1][y][0] → 0;
    laplaceEq2 = LaplaceTransformeq2Reorganized, t, s /.
    LaplaceTransformx[t], t, s → X[s], LaplaceTransformy[t], t, s → Y[s],
    x[0] → 0, y[0] → 0, Derivative[1][x][0] → 0, Derivative[1][y][0] → 0;

    laplaceEq1 // TraditionalForm
    laplaceEq2 // TraditionalForm

    (* Resolviendo el sistema de ecuaciones en el dominio de Laplace para X(s) y Y(s) *)

    Print"Resolviendo el sistema de ecuaciones en el dominio de Laplace:";


    solLaplace = SolvelaplaceEq1, laplaceEq2, X[s], Y[s];
    solLaplace // TraditionalForm

    (* Mostrar las expresiones de X(s) y Y(s) *)

    Xsol[s_] = X[s] /. solLaplace〚1〛;


    Ysol[s_] = Y[s] /. solLaplace〚1〛;

    Print"Soluciones en el dominio de Laplace para X(s) y Y(s):";


    Xsol[s] // TraditionalForm
    Ysol[s] // TraditionalForm

    (* Aplicando la Transformada Inversa de Laplace para encontrar x(t) e y(t) *)

    Print"Aplicando la Transformada Inversa de Laplace


    para obtener las soluciones en el dominio del tiempo:";
    xSol[t_] = InverseLaplaceTransformXsol[s], s, t;
    ySol[t_] = InverseLaplaceTransformYsol[s], s, t;

    xSol[t] // TraditionalForm
    ySol[t] // TraditionalForm

    (* Mostrar el resultado final de las posiciones en el tiempo *)

    Print"Las posiciones en el tiempo son:";


    Print"x(t) =";

    xSol[t] // TraditionalForm
    Print"y(t) =";

    ySol[t] // TraditionalForm
    3

    Ecuación de movimiento para la masa m1 derivada de la Segunda Ley de Newton:


    Out[124]//TraditionalForm=
    ′′
    x (t )  6 (y (t ) - x (t )) - 4 x (t )

    Ecuación de movimiento para la masa m2 derivada de la Segunda Ley de Newton:


    Out[127]//TraditionalForm=
    ′′
    y (t )  30 cos(t ) - 6 (y (t ) - x (t ))

    Ecuaciones simplificadas:
    Out[131]//TraditionalForm=
    ′′
    x (t )  6 (y (t ) - x (t )) - 4 x (t )
    Out[132]//TraditionalForm=
    ′′
    y (t )  30 cos(t ) - 6 (y (t ) - x (t ))

    Reorganización del sistema de ecuaciones:


    Out[136]//TraditionalForm=
    ′′
    x (t )  6 y (t ) - 10 x (t )
    Out[137]//TraditionalForm=
    ′′
    y (t )  6 x (t ) - 6 y (t ) + 30 cos(t )

    Aplicando la Transformada de Laplace a las ecuaciones:


    Out[142]//TraditionalForm=

    s 2 X (s )  6 Y (s ) - 10 X (s )
    Out[143]//TraditionalForm=

    30 s
    s 2 Y (s )  + 6 X (s ) - 6 Y (s )
    s2 + 1

    Resolviendo el sistema de ecuaciones en el dominio de Laplace:


    Out[146]//TraditionalForm=

    180 s 30 s s 2 + 10
    X (s ) → , Y (s ) → 
    s 6 + 17 s 4 + 40 s 2 + 24 s 6 + 17 s 4 + 40 s 2 + 24

    Soluciones en el dominio de Laplace para X(s) y Y(s):


    Out[150]//TraditionalForm=
    180 s

    s 6 + 17 s 4 + 40 s 2 + 24
    Out[151]//TraditionalForm=

    30 s s 2 + 10

    s 6 + 17 s 4 + 40 s 2 + 24

    Aplicando la Transformada Inversa de Laplace


    para obtener las soluciones en el dominio del tiempo:
    4

    Out[155]//TraditionalForm=

    cos(t )
    180 +
    9

    2 10 cos 2 4 - 10  t + 7 cos 2 4 - 10  t +2 10 cos 2 4 + 10  t - 7 cos 2 4 + 10  t

    18 4- 10 - 4+ 10 4- 10 + 4+ 10

    Out[156]//TraditionalForm=

    10 cos 2 4 - 10  t + 3 cos 2 4 - 10  t + 10 cos 2 4 + 10  t - 3 cos 2 4 + 10  t


    30 cos(t ) +
    4- 10 - 4+ 10 4- 10 + 4+ 10

    Las posiciones en el tiempo son:


    x(t) =
    Out[159]//TraditionalForm=

    cos(t )
    180 +
    9

    2 10 cos 2 4 - 10  t + 7 cos 2 4 - 10  t +2 10 cos 2 4 + 10  t - 7 cos 2 4 + 10  t

    18 4- 10 - 4+ 10 4- 10 + 4+ 10

    y(t) =
    Out[161]//TraditionalForm=

    10 cos 2 4 - 10  t + 3 cos 2 4 - 10  t + 10 cos 2 4 + 10  t - 3 cos 2 4 + 10  t


    30 cos(t ) +
    4- 10 - 4+ 10 4- 10 + 4+ 10

    También podría gustarte