Discussion :

[onsmateriau]

Bonjour,
je travaille sur la simulation d'un essai de traction sous déformation imposée c.à.d j'impose un cycle triangulaire (charge-décharge) de déformation en fonction du temps (figure 1). Pour ce cycle j'ai pas trouvé un problème avec fortran de le programmer; je le programme de la façon suivante:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
PARAMETER (pente=1.34d-04,periode=1044.35d0,temp1=522.175d0)
...
      DO WHILE (temp.LE.periode)  ! DEBUT BOUCLE TEMPS
 
        IF (temp.LE.temp1) THEN
            pendefo=pente
        ELSE IF (temp.GT.temp1) THEN        
            pendefo=-pente
        END IF
 
c Déformation totale IMPOSEE 'depstot'
      depstot=pendefo*deltat
.....

Toutefois, si je veux imposer le chargement réel comme en expérimental (figure 2) je ne sais pas comment le faire.
Merci d'avance
Pouvez vous m'aider svp?Pièce jointe 203472Pièce jointe 203473Pièce jointe 203472Pièce jointe 203473

[dva2tlse]

Bonjour,
j'aurais bien aimé t'aider un peu puisque c'est mon boulot de calculateur de structures, donc je connais un peu la RDM et les essais, et le fortran; malheureusement la pièce jointe est illisible.
David

[onsmateriau]

Merci pour votre réponse David,
je vous joints de nouveau les 2 figures: la 1ère est celle du cycle triangulaire (figure1) dont je l'ai appliqué auparavant (celui imposé dans le programme fortran que j'ai présenté) et le 2ème (figure2) est celui que je souhaite imposer (celui expérimental).

[dva2tlse]

Ok super, j'étais encore devant mon PC quand la notification de votre message est arrivée.
Ça ne va pas être bien compliqué, il suffit d'échantillonner le triangle que vous voulez appliquer, je vais le faire sommairement, libre à vous d'affiner par la suite : (j'écris comme du code pour avoir une police fixe, et puis il y a vraiment du code par la suite)

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
 
 temps (s) epstot (x 100% je suppose)
   0         0.
 100         0.01
 200         0.01
 300         0.016
 400         0.04
 500         0.07
 600         0.06
 700         0.035
 800         0.02
 900         0.005
1000         0.

Ensuite il suffit de rentrer ces valeurs dans une table, par une instruction DATA, juste après les déclarations du début du programme, et avant la première instruction exécutable

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
c----------------------------------------------------------------------| -----|
      real(eps(11))
      DATA eps/O., 0.01, 0.01,  0.016, 0.04, 0.07, 
     +             0.06, 0.035, 0.02,  0.02, 0.005, 0./
      T=0
      deltaT=100
      defoinit=0.
      DO WHILE (T.LE.1000)  ! DEBUT BOUCLE TEMPS
c dans chaque intervalle :
          N=int(T/100.)
          pente=(eps(N+1)-eps(N))/deltaT
          defofinale=defoinit+pente*deltaT
          print*, defofinale ! intermédiaire
          T=T+deltaT
          defoinit=defofinale
      ENDDO
c Déformation totale IMPOSEE 'depstot'
      print*, defofinale ! finale
 
c----------------------------------------------------------------------| -----|

Si vous avez réussi à faire le premier petit bout de code que j'ai vu, vous devriez savoir adapter ça à votre besoin.
À+,
David

EDIT: Il y a douze réels dans le tableau, donc il faut déclarer :
real eps(12)
EDIT2: Avant de boucler sur le prochain pas de temps, juste avant le "ENDDO" donc, il faut mettre à jour :
defoinit=defofinale
Bye.

[onsmateriau]

Merci pour votre réponse David,
j'obtiens ce message d'erreur
la table 'eps' accepte une seule valeur sinon j'obtiens cette erreur.
Pouvez vous m'aider?
Merci encore une fois

[dva2tlse]

Bonjour,
vous avez de la chance que je rentre tôt du boulot.
J'ai été obligé de faire tourner le programme chez moi pour trouver l'erreur, donc je vais le coller en bas de ce message et l'envoyer en pièce jointe. (mais en ".txt" sinon ça ne part pas)
Le premier terme de eps atit un 0 (lettre "o" majuscule) et non pas 0 (chiffre zéro).
Maintenant ça marche pour les chiffres mais je ne sais pas comment le dessiner. La façon de le faire dépend de ce que peut faire votre machine, et je ne suis pas du tout un as dans ce domaine.
Bon courage,
David

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
c echo; echo; echo; gfortran FDEV.f -o FDEV
c ./FDEV
      program FDEV
      implicit none
      integer deltaT
      integer N
      integer T
      real defofinale
      real defoinit
      real depstot
      real eps(12)
      real pente
      DATA eps/0.0, 0.01, 0.01,  0.016, 0.04, 0.07, 
     +              0.06, 0.035, 0.02,  0.02, 0.005, 0./
c
      T=0
      deltaT=100
      defoinit=0.
      DO WHILE (T.LE.1000)  ! DEBUT BOUCLE TEMPS
c dans chaque intervalle :
          N=int(T/100.)
          pente=(eps(N+1)-eps(N))/deltaT
          defofinale=defoinit+pente*deltaT
          print'(A,F8.3)', 'FDEV: defo=',  defofinale
          defoinit=defofinale
          T=T+deltaT
      ENDDO
      END
cEOF