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

In [1]:

    

Vd_cls = 15.68288727


    

In [2]:

    

Vd_enr = 16.59350136


    

In [3]:

    

Vg_cls = -15.68288727


    

In [4]:

    

Vg_enr = 16.59350136


    

In [5]:

    

dt = 4.0e-09


    

In [6]:

    

xd_t = 0.00510397


    

In [7]:

    

xg_t = 0.00489603


    

In [8]:

    

Vg =  Vg_cls - Vg_enr


    

In [9]:

    

Vd = Vd_cls + Vd_enr


    

In [10]:

    

xd_tpdt = xd_t + Vd * dt


    

In [11]:

    

xg_tpdt = xg_t + Vg * dt


    

In [12]:

    

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


    

    




Vitesse du noeud gauche :     -32.2763886

In [13]:

    

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


    

    




Vitesse du noeud droit :      32.2763886

In [14]:

    

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


    

    




coord à t + dt du noeud gauche :   0.00489590089

In [15]:

    

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


    

    




coord à t + dt du noeud droit :   0.00510409911

In [ ]: