Calcul des nouvelles vitesses et coodonnées des noeuds de l'élément enrichi (Eq 3.80 à 3.82)

In [29]:

    

Vd_cls = 15.68288727


    

In [30]:

    

Vd_enr = 16.59350136


    

In [31]:

    

Vg_cls = -15.68288727


    

In [32]:

    

Vg_enr = 16.59350136


    

In [33]:

    

dt = 4.0e-09


    

In [34]:

    

xd_t = 0.00510397


    

In [35]:

    

xg_t = 0.00489603


    

In [36]:

    

Vg =  Vg_cls - Vg_enr


    

In [37]:

    

Vd = Vd_cls + Vd_enr


    

In [38]:

    

xd_tpdt = xd_t + Vd * dt


    

In [39]:

    

xg_tpdt = xg_t + Vg * dt


    

In [40]:

    

print "Vitesse du noeud gauche : {:15.9g}".format(Vg)


    

    




Vitesse du noeud gauche :     -32.2763886

In [41]:

    

print "Vitesse du noeud droit : {:15.9g}".format(Vd)


    

    




Vitesse du noeud droit :      32.2763886

In [42]:

    

print "coord à t + dt du noeud gauche : {:15.9g}".format(xg_tpdt)


    

    




coord à t + dt du noeud gauche :   0.00489590089

In [43]:

    

print "coord à t + dt du noeud droit : {:15.9g}".format(xd_tpdt)


    

    




coord à t + dt du noeud droit :   0.00510409911

In [ ]: