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

In [22]:

    

Lg_n = 0.000103908305743


    

In [36]:

    

Ld_n = 0.000103908305743


    

In [37]:

    

Vd_cls = 15.68288727


    

In [38]:

    

Vd_enr = 16.59350136


    

In [39]:

    

Vg_cls = -15.68288727


    

In [40]:

    

Vg_enr = 16.59350136


    

In [41]:

    

RHOg_n = 7664.92796304


    

In [42]:

    

RHOd_n = 7664.92796304


    

In [43]:

    

dt = 4.0e-09


    

In [44]:

    

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


    

In [45]:

    

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


    

    




Taille gauche à t+dt =  0.000103971037

In [50]:

    

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


    

In [51]:

    

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


    

    




Taille droite à t+dt =  0.000103971037

In [48]:

    

RHOg_npun = RHOg_n * Lg_n / Lg_npun


    

In [49]:

    

RHOd_npun = RHOd_n * Ld_n / Ld_npun


    

In [52]:

    

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


    

    




Densite gauche à t+dt =      7660.30328

In [53]:

    

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


    

    




Densite droite à t+dt =      7660.30328

In [ ]: