Calcul des nouvelles tailles des parties sectionnées des éléments enrichis (Eq 3.88 et 3.89)

In [1]:

    

Lg_n = 0.000103908305743


    

In [2]:

    

Ld_n = 0.000103908305743


    

In [3]:

    

Vd_cls = 15.68288727


    

In [4]:

    

Vd_enr = 16.59350136


    

In [5]:

    

Vg_cls = -15.68288727


    

In [6]:

    

Vg_enr = 16.59350136


    

In [7]:

    

RHOg_n = 7664.92796304


    

In [8]:

    

RHOd_n = 7664.92796304


    

In [9]:

    

dt = 4.0e-09


    

In [10]:

    

Lg_npun = Lg_n + ((Vd_cls - Vd_enr) * 0.5 - (Vg_cls - Vg_enr) * 0.5) * dt


    

In [11]:

    

print "Taille gauche à t+dt = {:15.9g}".format(Lg_npun)


    

    




Taille gauche à t+dt =  0.000103971037

In [12]:

    

Ld_npun = Ld_n + ((Vd_cls + Vd_enr) * 0.5 - (Vg_cls + Vg_enr) * 0.5) * dt


    

In [13]:

    

print "Taille droite à t+dt = {:15.9g}".format(Ld_npun)


    

    




Taille droite à t+dt =  0.000103971037

In [14]:

    

RHOg_npun = RHOg_n * Lg_n / Lg_npun


    

In [15]:

    

RHOd_npun = RHOd_n * Ld_n / Ld_npun


    

In [16]:

    

print "Densite gauche à t+dt = {:15.9g}".format(RHOg_npun)


    

    




Densite gauche à t+dt =      7660.30328

In [17]:

    

print "Densite droite à t+dt = {:15.9g}".format(RHOd_npun)


    

    




Densite droite à t+dt =      7660.30328

In [ ]: