Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 9 vistas

    Ejemplo Sim Ula Do

    Cargado por

    Akiles castro
    Este documento describe el movimiento de una partícula cargada sometida a campos eléctrico y magnético perpendiculares. Se utilizan las ecuaciones de Lagrange para derivar las ecuaciones de movimiento y se resuelve numéricamente para graficar la trayectoria en forma de espiral confinada a un plano.

    Copyright:

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

    Ejemplo Sim Ula Do

    Cargado por

    Akiles castro
    9 vistas4 páginas
    Este documento describe el movimiento de una partícula cargada sometida a campos eléctrico y magnético perpendiculares. Se utilizan las ecuaciones de Lagrange para derivar las ecuaciones de movimiento y se resuelve numéricamente para graficar la trayectoria en forma de espiral confinada a un plano.

    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?

    Este documento describe el movimiento de una partícula cargada sometida a campos eléctrico y magnético perpendiculares. Se utilizan las ecuaciones de Lagrange para derivar las ecuaciones de movimiento y se resuelve numéricamente para graficar la trayectoria en forma de espiral confinada a un plano.

    Copyright:

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

    Ejemplo Sim Ula Do

    Cargado por

    Akiles castro
    Este documento describe el movimiento de una partícula cargada sometida a campos eléctrico y magnético perpendiculares. Se utilizan las ecuaciones de Lagrange para derivar las ecuaciones de movimiento y se resuelve numéricamente para graficar la trayectoria en forma de espiral confinada a un plano.

    Copyright:

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

    Movimiento en campos eléctrico y magnético cruzados

    En esta página, vamos a estudiar el movimiento de una partícula de masa m, y carga q, sometida a la
    acción simultánea de un campo eléctrico E y de un campo magnético B, ambos uniformes y perpendicu
    lares entre sí. Queremos determinar a rayectoria de la partícula y hacer una animación del movimiento

    Solución
    Se trata de un problema de 3 grados de libertad. Usaremos como coordenadas generalizadas las
    cartesianas, y supongamos que el campo magnético B tiene la dirección del eje Z, el campo eléctrico
    E la dirección del eje Y y el vector velocidad v está en el plano XY. La partícula cargada parte de la
    posición inicial (x0 , y0 ) con velocidad inicial (v0 x , v0 y ).
    Se trata de un problema con 3 grados de libertad y vamos usar como coordenadas generalizadas las
    Cartesias. Es sabido que en general el Lagrangiano de una partícula en un CEM está dado por
    m 2 2 2
    L= x + y + z  - q ϕ + q A · v (1)
    2
    Los campos E(r, t) y B(r, t) son funciones contínuas en sus argumentos. Estos campos se relacionan
    con el potencial escalar ϕ (r, t) y el potencial vectorial A(r, t) a través de (Pag 39. Ec 2.30)
    ∂A
    E = - ∇ϕ - ; B = ∇ ×A (2)
    ∂t
    Para nuestro caso, como E = E j y B = B k son uniformes, estos potenciales se reducen a
    1
    ϕ = - E y; A = B×r (3)
    2
    1
    BB = {0, 0, B}; r[t_] := {x[t], y[t], z[t]}; A = Cross[BB, r[t]]
    2 producto vectorial

    1 1
    - B y[t], B x[t], 0
    2 2
    m
    v[t_] := r '[t]; L = v[t].v[t] - q (- EE y[t]) + q A.v[t]
    2
    1 1 1
    EE q y[t] + q - B y[t] x′ [t] + B x[t] y′ [t] + m x′ [t]2 + y′ [t]2 + z′ [t]2 
    2 2 2
    2 EjemploSimulado.nb

    Ecuaciones de Movimiento
    ⅆ ∂L ∂L ⅆ ∂L ∂L ⅆ ∂L ∂L
     - = 0;  - = 0;  - ∂z = 0 (4)
    ⅆt ∂x ∂x ⅆt ∂y ∂y ⅆt ∂z
    ecuX = D[D[L, x '[t]], t] - D[L, x[t]] ⩵ 0
    ⋯ deriva deriva
    ecuY = D[D[L, y '[t]], t] - D[L, y[t]] ⩵ 0
    ⋯ deriva deriva
    ecuZ = D[D[L, z '[t]], t] - D[L, z[t]] ⩵ 0
    ⋯ deriva deriva

    - B q y′ [t] + m x′′ [t] ⩵ 0

    - EE q + B q x′ [t] + m y′′ [t] ⩵ 0

    m z′′ [t] ⩵ 0

    Solución Numérica
    Introducimos las siguientes definiciones ω = q B/m, vd = E /B, ω vd = q E/m

    In[1]:= sistema = {- w y′ [t] + x′′ [t] ⩵ 0, - w vd + w x′ [t] + y′′ [t] ⩵ 0, z′′ [t] ⩵ 0};

    In[2]:= w = 1; vd = 0.05;

    In[3]:= cond = {x[0] == - 0.8, x '[0] ⩵ 0, y[0] ⩵ 0, y '[0] ⩵ 0.3, z[0] ⩵ 0, z '[0] ⩵ 0};

    In[4]:= sol = NDSolve[sistema ⋃ cond, {x, y, z}, {t, 0, 10 π / w}]


    resolvedor diferencial numérico

    Out[4]= x → InterpolatingFunction Domain: {{0., 31.4}} ,


    Output: scalar

    y → InterpolatingFunction Domain: {{0., 31.4}} ,


    Output: scalar

    z → InterpolatingFunction Domain: {{0., 31.4}} 


    Output: scalar

    Las soluciones son:

    In[5]:= X[t_] := x[t] /. sol[[1]];


    Y[t_] := y[t] /. sol[[1]];
    Z[t_] := z[t] /. sol[[1]];
    EjemploSimulado.nb 3

    In[8]:= Fig1 = Plot[{X[t], Y[t], Z[t]}, {t, 0, 6 π / w}, PlotLegends → "Expressions"]


    representación gráfica leyendas de representación

    0.6

    0.4

    0.2

    X(t)
    5 10 15
    Out[8]= Y(t)
    -0.2
    Z(t)
    -0.4

    -0.6

    -0.8

    Se observa que la coordenada Z se mantiene siempre igual a cero (consecuencia de las condiciones
    iniciales y no de una ligadura) y por lo tanto el movimiento está confinado al plano XY. Costruyamos la
    trayectoria.

    In[51]:= Clear[trayectoria]
    borra

    In[56]:= trayectoria[t_] := ParametricPlot[{X[tp], Y[tp]},


    gráfico paramétrico
    {tp, 0, t}, PlotRange → {{- 0.8, 1}, {- 0.3, 0.4}}, PlotStyle → Dashed]
    rango de representación estilo de repre⋯ rayado

    Hagamos la animación del movimiento

    In[57]:= particula[t_] := Graphics[{Gray, Disk[{X[t], Y[t]}, 0.02]}];


    gráfico gris disco
    TrayPart[t_] := Show[trayectoria[t], particula[t]]
    muestra
    4 EjemploSimulado.nb

    In[59]:= Manipulate[TrayPart[t], {t, 0.001, 10 π / w}]


    manipula

    20.1685

    0.4

    0.3
    Out[59]=
    0.2

    0.1

    -0.5 0.5 1.0


    -0.1

    -0.2

    -0.3

    También podría gustarte