Reprise du message précédent :
comment avec CLOAD .??? tu peux m'expliquer ? stp  

j'ai mon fichier inp ya 7000noeuds et j'ai les valeurs de chaque contrainte de noueds c-t dire  
*Heading
** Job name: cubetorsion Model name: Model-1
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*Node
      1,  -1.33790886,           0.,           0.
      2,           0.,           0.,  0.100000001
      3,           0., -0.100000001,           0.
      4,           0.,           0.,          10.
      5,          -5.,           0.,           0.
      6,          -5.,           0.,          10.
      7,           0.,          -5.,           0.
      8,          -5.,          -5.,           0.
      9,          -5.,          -5.,          10.
     10,           0.,          -5.,          10.
     11,   1.33790886,           0.,           0.
     .
  jusque 8000noueds et j'ai une fonction qui me donne les valeure de les contraintes a chaque noeuds . donc le pb applique cette contrainte a chaques noeud  j'arrive pas ???. est ceque on peut faire ca directement ds le fichier .inp si oui comment?.................
 
 
 
 
sinon une autre qustion ya une subroutune sigini comment elle marche ????? et comment je peux appller une subroutune (ajouter la ds mon fichier .inp)      

Ben dans le fichier inp , pour applique une contrainte sur un noeud ou un ensemble de noeud, tu as la fonction CLOAD.
Mais dasn ton cas, c'est galère car ca change à chaque noeud
Je crois aue tu es bon pour utiliser la subroutine UTRACLOAD qui te permet d'appliquer un chargement variable à chaque noeud.
Regarde dasn l'aide comment elle se construit, c'est un peu délicat    
 
Pour appeller une subroutine, il faut utiliser le code ,USER comme option de la commande qui doit appeler la subroutine. Bon je connais pas la commade qui appelle cette subroutine, moi j'utilise une subroutine pour changer le maillage. Faut que tu regardes dans l'aide
Après, je rentre le nom du fichier fortran qui contient la subroutine en paramètre d'abaqus

jai regarde ds le aide mais j'ai rien compris  
  SUBROUTINE UTRACLOAD(ALPHA,T_USER,KSTEP,KINC,TIME,NOEL,NPT,
     1 COORDS,DIRCOS,JLTYP,SNAME)
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION T_USER(3), TIME(2), COORDS(3), DIRCOS(3,3)
      CHARACTER*80 SNAME
 
      user coding to define ALPHA and T_USER
 
 
      RETURN
      END

voila ma fonction  f(i)=1E-08*p(i)**5 - 1E-05*p(i)**4 + 0.0038*p(i)**3 - 0.5377*p(i)**2 + 32.733*p(i) - 1.
je dois appliquer cette fonction (suivant les profondeur z).je fais un programmme sur fortran qui me donne les resultat .mais le pb je sais pas comment je peux appliquer une contrain,te a neuds .

Ca te donne le début et la fin de la subroutine.  
Regarde la suite de l'aide. Les variables à indiquer sont alpha et T_user
 
Le alpha est l'amplitude du chargement. C'est ta fonction en fait.
T_user est la direction du chargement
 
Les autres variables sont données par abaqus. tu peux les utiliser.
Par conte, je viens de voir que tu peux pas travailler avec les noeuds avec cette subroutine.
 
Mais il y a un autre moyen en utilisant amplitude et Cload.
 
Est-ce que les noeuds qui recoivent la force sont sur la meme surface ?
Car tu pourrais avec fortran sortir la liste de valeurs de force correspondant à chaque noeud.
Imagine que tu as un ensemble de noeuds : 34 65 12 94 52 ....
avec fortran, tu lis cet ensemble. tu fais ton calcul et tu sors dasn un fichier les valeurs a b c d e f etc correspondant aux noeuds
Après avec la fonction amplitude, tu lis ces valeurs
Et enfin tu envoies cette amplitude dans la fonction Cload.
 
Enfin ca c'est le principe, je suis pas sur que ca marche.
De toute manière commence par lire en entier l'aide sur cload, amplitude et les subroutines. si tu sais pas comment ca marche, t'as aucune chance d'y arriver

j'ai vu ds le doc que je peux utilser la subroutune SIGINI. mais j'ai rien compris tu peux m'expliquer un peu . sinon le voila mais je sais pas ou je dois mettre ma fonction .f(i)=1E-08*p(i)**5 - 1E-05*p(i)**4 + 0.0038*p(i)**3 - 0.5377*p(i)**2 + 32.733*p(i) - 1.
sinon voila la subroutune SIGINI .    
 
 
 
SUBROUTINE SIGINI(SIGMA,COORDS,NTENS,NCRDS,NOEL,NPT,LAYER,
     1 KSPT,LREBAR,NAMES)
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION SIGMA(NTENS),COORDS(NCRDS)
      CHARACTER NAMES(2)*80
 
 
      user coding to define SIGMA(NTENS)
 
 
      RETURN
      END
Variables to be defined
 
SIGMA(1)
 
First stress component.
 
 
SIGMA(2)
 
Second stress component.
 
 
SIGMA(3)
 
Third stress component.
 
 
Etc.
 
Only NTENS stress values should be defined, where NTENS depends on the element type.
 
 
Variables passed in for information
 
COORDS
 
An array containing the initial coordinates of this point.
 
 
NTENS
 
Number of stresses to be defined, which depends on the element type.
 
 
NCRDS
 
Number of coordinates.
 
 
NOEL
 
Element number.
 
 
NPT
 
Integration point number in the element.
 
 
LAYER
 
Layer number (for composite shells and layered solids).
 
 
KSPT
 
Section point number within the current layer.
 
 
LREBAR
 
Rebar flag. If LREBAR=1, the current integration point is associated with element rebar. Otherwise, LREBAR=0.
 
 
NAMES(1)
 
Name of the rebar to which the current integration point belongs, which is the name given in the rebar or rebar layer definition (“Defining reinforcement,” Section 2.2.3 of the ABAQUS Analysis User's Manual, or “Defining rebar as an element property,” Section 2.2.4 of the ABAQUS Analysis User's Manual). If no name was given in the rebar or rebar layer definition, this variable will be blank. This variable is relevant only when LREBAR=1.
 
 
NAMES(2)
 
Element type name (see Section I.1, “ABAQUS/Standard Element Index,” of the ABAQUS Analysis User's Manual).
 
 
 
 

Ta fonction, elle donne la contrainte à un noeud donné non ?
c'est quoi le i de ta fonction ?
 
Cette subroutine va être appelée pour chaque élement.
Les variables qu'il faut que tu rentres sont les sigma.
Le problème, c'est qu'avec cette subroutine, tu définis les composantes d'efforts d'un élément et non d'un noeud !
 
Si tu veux l'utiliser, tu dois transformer ta formule pour qu'elle donne les contraintes dans un élément. Je sais pas quel type d'element tu as (2d, 3d, à interpolation linéaire ou pas etc) mais il est possible de connaitre le tenseur des contraintes dans un élément à partir des contraintes aux noeuds.
 
Cette subroutine demande en fait les valeurs de ton tenseur des contraintes pour l'élément. sigma(1), c'est sigma xx en fait, etc.  
C'est pas forcément évident

j'ai une fonction contra(z)=z+1.
cette fonction me donne les contraintes a chaque profondeur (z)
voila un exmpls simple pr bien comprendre  
voila comment il est un fichier .inp  
 
neouds             x                y              z      
1                    0                0              5
2                    5                3             0  
3
4
.                                                          
 
 
 
.
..
6000                                . ...................
je vais appliquer ma fonction sur chaque valeurs en z  
 
 
 
 
 
 
 
donc j'aurai 1 tableau  ou  y a que mes nouds et mes les valeurs de mes contraintes ...... le pb comment je peux traduire ca ds abaqus . mon tableau ce pas un fichier de type .inp
 
 
voila un exampl de mon inp
 
*Heading
** Job name: cubetorsion Model name: Model-1
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*Node
      1,  -1.33790886,           0.,           0.
      2,           0.,           0.,  0.100000001
      3,           0., -0.100000001,           0.
      4,           0.,           0.,          10.
      5,          -5.,           0.,           0.
................................................... jusque 8000neuds
*Element, type=C3D4
   1, 3510, 3511, 3512, 3513
   2, 3511, 3514, 3513, 3510
.............................................................
** Region: (Section-1:Picked)
*Elset, elset=_PickedSet13, internal, generate
     1,  28619,      1
** Section: Section-1
*Solid Section, elset=_PickedSet13, material=acier
1.,
*End Part
........................................................................................
 
 
 
 
 
 
 
 
 
 
 
 
 

Mais tes contraintes, elles sont suivant quelle direction ?
 
Pour appliquer une contrainte à un noeud en particulier, il faut utiliser Cload.
donc tu peux faire un programme fortran qui écrira dans un fichier texte :
*CLOAD
N,x,C
 
avec N : le numéro du noeud
x : la direction de la force
C : la valeur de la force
 
et tu fais ca pour chaque noeud
 
*CLOAD
1,1,0.313123
*CLOAD
2,1,0.34234
*CLOAD
3,1,-0.634636
 
Etc pour tous les noeuds. Après tu copies tous les Cload que ton programme fortran a sorti et tu les colles dans ton fichier inp

ma fonction ce une contrainte résuduelles

ma fonction ce une contrainte résudelles. sinon j'ai vu ds le doc la subroutune SIGINI mais tu peux m'expliquer cette SIGINI et ou je peux mettre ma fonction ?????

Mais une contrainte dans quelle direction ?
Et pourquoi tu veux absolument utiliser sigini ? Je t'ai dit que c'était une subtroutine qui s'appliquait sur des éléments et non des noeuds. Tu ne peux pas utiliser ta formule directement. Il faut que tu calcules les contraintes dans un éléments à partir de celles aux noeuds (que tu as avec ta formule) et intégrant aux points de Gauss !

mon ami j'ai un tableau ds lequel ya mes numero de noeuds et les contraintes resu suivant 2. (j'ai que s22)

Ah ben voilà, c'est simple
 
Tu fais autant de CLOAD que de noeuds.
 
*CLOAD
numéro du noeud,2,la valeur
 
Et tu répètes ça pour tous les noeuds.
Bien sur, tu vas pas le faire à la main, donc il faut que tu écrives un programme qui te sorte un fichier texte contenant tous ces CLOAD. Pas besoin de subroutine.

voila mon fichier inp .je voisq pas ou je dois mettre mes CLOAD
je fais
*CLOAD
1, 5
2, -4
3, 6  
....... jusque 6000noeuds
 
 
 
*Heading
** Job name: essai Model name: Model-1
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*Node
      1,          -5.,          -5.,          10.
      2,          -5.,         -2.5,          10.
      3,          -5.,           0.,          10.
      4,          -5.,          2.5,          10.
      5,          -5.,           5.,          10.
      6,          -5.,          -5.,          7.5
      7,          -5.,         -2.5,          7.5
      8,          -5.,           0.,          7.5
      9,          -5.,          2.5,          7.5
     10,          -5.,           5.,          7.5
     11,          -5.,          -5.,           5.
     12,          -5.,         -2.5,           5.
     13,          -5.,           0.,           5.
     14,          -5.,          2.5,           5.
     15,          -5.,           5.,           5.
     16,          -5.,          -5.,          2.5
     17,          -5.,         -2.5,          2.5
     18,          -5.,           0.,          2.5
     19,          -5.,          2.5,          2.5
     20,          -5.,           5.,          2.5
     21,          -5.,          -5.,           0.
     22,          -5.,         -2.5,           0.
     23,          -5.,           0.,           0.
     24,          -5.,          2.5,           0.
     25,          -5.,           5.,           0.
     26,         -2.5,          -5.,          10.
     27,         -2.5,         -2.5,          10.
     28,         -2.5,           0.,          10.
     29,         -2.5,          2.5,          10.
     30,         -2.5,           5.,          10.
     31,         -2.5,          -5.,          7.5
     32,         -2.5,         -2.5,          7.5
     33,         -2.5,           0.,          7.5
     34,         -2.5,          2.5,          7.5
     35,         -2.5,           5.,          7.5
     36,         -2.5,          -5.,           5.
     37,         -2.5,         -2.5,           5.
     38,         -2.5,           0.,           5.
     39,         -2.5,          2.5,           5.
     40,         -2.5,           5.,           5.
     41,         -2.5,          -5.,          2.5
     42,         -2.5,         -2.5,          2.5
     43,         -2.5,           0.,          2.5
     44,         -2.5,          2.5,          2.5
     45,         -2.5,           5.,          2.5
     46,         -2.5,          -5.,           0.
     47,         -2.5,         -2.5,           0.
     48,         -2.5,           0.,           0.
     49,         -2.5,          2.5,           0.
     50,         -2.5,           5.,           0.
     51,           0.,          -5.,          10.
     52,           0.,         -2.5,          10.
     53,           0.,           0.,          10.
     54,           0.,          2.5,          10.
     55,           0.,           5.,          10.
     56,           0.,          -5.,          7.5
     57,           0.,         -2.5,          7.5
     58,           0.,           0.,          7.5
     59,           0.,          2.5,          7.5
     60,           0.,           5.,          7.5
     61,           0.,          -5.,           5.
     62,           0.,         -2.5,           5.
     63,           0.,           0.,           5.
     64,           0.,          2.5,           5.
     65,           0.,           5.,           5.
     66,           0.,          -5.,          2.5
     67,           0.,         -2.5,          2.5
     68,           0.,           0.,          2.5
     69,           0.,          2.5,          2.5
     70,           0.,           5.,          2.5
     71,           0.,          -5.,           0.
     72,           0.,         -2.5,           0.
     73,           0.,           0.,           0.
     74,           0.,          2.5,           0.
     75,           0.,           5.,           0.
     76,          2.5,          -5.,          10.
     77,          2.5,         -2.5,          10.
     78,          2.5,           0.,          10.
     79,          2.5,          2.5,          10.
     80,          2.5,           5.,          10.
     81,          2.5,          -5.,          7.5
     82,          2.5,         -2.5,          7.5
     83,          2.5,           0.,          7.5
     84,          2.5,          2.5,          7.5
     85,          2.5,           5.,          7.5
     86,          2.5,          -5.,           5.
     87,          2.5,         -2.5,           5.
     88,          2.5,           0.,           5.
     89,          2.5,          2.5,           5.
     90,          2.5,           5.,           5.
     91,          2.5,          -5.,          2.5
     92,          2.5,         -2.5,          2.5
     93,          2.5,           0.,          2.5
     94,          2.5,          2.5,          2.5
     95,          2.5,           5.,          2.5
     96,          2.5,          -5.,           0.
     97,          2.5,         -2.5,           0.
     98,          2.5,           0.,           0.
     99,          2.5,          2.5,           0.
    100,          2.5,           5.,           0.
    101,           5.,          -5.,          10.
    102,           5.,         -2.5,          10.
    103,           5.,           0.,          10.
    104,           5.,          2.5,          10.
    105,           5.,           5.,          10.
    106,           5.,          -5.,          7.5
    107,           5.,         -2.5,          7.5
    108,           5.,           0.,          7.5
    109,           5.,          2.5,          7.5
    110,           5.,           5.,          7.5
    111,           5.,          -5.,           5.
    112,           5.,         -2.5,           5.
    113,           5.,           0.,           5.
    114,           5.,          2.5,           5.
    115,           5.,           5.,           5.
    116,           5.,          -5.,          2.5
    117,           5.,         -2.5,          2.5
    118,           5.,           0.,          2.5
    119,           5.,          2.5,          2.5
    120,           5.,           5.,          2.5
    121,           5.,          -5.,           0.
    122,           5.,         -2.5,           0.
    123,           5.,           0.,           0.
    124,           5.,          2.5,           0.
    125,           5.,           5.,           0.
*Element, type=C3D8
 1,  26,  27,  32,  31,   1,   2,   7,   6
 2,  27,  28,  33,  32,   2,   3,   8,   7
 3,  28,  29,  34,  33,   3,   4,   9,   8
 4,  29,  30,  35,  34,   4,   5,  10,   9
 5,  31,  32,  37,  36,   6,   7,  12,  11
 6,  32,  33,  38,  37,   7,   8,  13,  12
 7,  33,  34,  39,  38,   8,   9,  14,  13
 8,  34,  35,  40,  39,   9,  10,  15,  14
 9,  36,  37,  42,  41,  11,  12,  17,  16
10,  37,  38,  43,  42,  12,  13,  18,  17
11,  38,  39,  44,  43,  13,  14,  19,  18
12,  39,  40,  45,  44,  14,  15,  20,  19
13,  41,  42,  47,  46,  16,  17,  22,  21
14,  42,  43,  48,  47,  17,  18,  23,  22
15,  43,  44,  49,  48,  18,  19,  24,  23
16,  44,  45,  50,  49,  19,  20,  25,  24
17,  51,  52,  57,  56,  26,  27,  32,  31
18,  52,  53,  58,  57,  27,  28,  33,  32
19,  53,  54,  59,  58,  28,  29,  34,  33
20,  54,  55,  60,  59,  29,  30,  35,  34
21,  56,  57,  62,  61,  31,  32,  37,  36
22,  57,  58,  63,  62,  32,  33,  38,  37
23,  58,  59,  64,  63,  33,  34,  39,  38
24,  59,  60,  65,  64,  34,  35,  40,  39
25,  61,  62,  67,  66,  36,  37,  42,  41
26,  62,  63,  68,  67,  37,  38,  43,  42
27,  63,  64,  69,  68,  38,  39,  44,  43
28,  64,  65,  70,  69,  39,  40,  45,  44
29,  66,  67,  72,  71,  41,  42,  47,  46
30,  67,  68,  73,  72,  42,  43,  48,  47
31,  68,  69,  74,  73,  43,  44,  49,  48
32,  69,  70,  75,  74,  44,  45,  50,  49
33,  76,  77,  82,  81,  51,  52,  57,  56
34,  77,  78,  83,  82,  52,  53,  58,  57
35,  78,  79,  84,  83,  53,  54,  59,  58
36,  79,  80,  85,  84,  54,  55,  60,  59
37,  81,  82,  87,  86,  56,  57,  62,  61
38,  82,  83,  88,  87,  57,  58,  63,  62
39,  83,  84,  89,  88,  58,  59,  64,  63
40,  84,  85,  90,  89,  59,  60,  65,  64
41,  86,  87,  92,  91,  61,  62,  67,  66
42,  87,  88,  93,  92,  62,  63,  68,  67
43,  88,  89,  94,  93,  63,  64,  69,  68
44,  89,  90,  95,  94,  64,  65,  70,  69
45,  91,  92,  97,  96,  66,  67,  72,  71
46,  92,  93,  98,  97,  67,  68,  73,  72
47,  93,  94,  99,  98,  68,  69,  74,  73
48,  94,  95, 100,  99,  69,  70,  75,  74
49, 101, 102, 107, 106,  76,  77,  82,  81
50, 102, 103, 108, 107,  77,  78,  83,  82
51, 103, 104, 109, 108,  78,  79,  84,  83
52, 104, 105, 110, 109,  79,  80,  85,  84
53, 106, 107, 112, 111,  81,  82,  87,  86
54, 107, 108, 113, 112,  82,  83,  88,  87
55, 108, 109, 114, 113,  83,  84,  89,  88
56, 109, 110, 115, 114,  84,  85,  90,  89
57, 111, 112, 117, 116,  86,  87,  92,  91
58, 112, 113, 118, 117,  87,  88,  93,  92
59, 113, 114, 119, 118,  88,  89,  94,  93
60, 114, 115, 120, 119,  89,  90,  95,  94
61, 116, 117, 122, 121,  91,  92,  97,  96
62, 117, 118, 123, 122,  92,  93,  98,  97
63, 118, 119, 124, 123,  93,  94,  99,  98
64, 119, 120, 125, 124,  94,  95, 100,  99
*Nset, nset=_PickedSet2, internal, generate
   1,  125,    1
*Elset, elset=_PickedSet2, internal, generate
  1,  64,   1
** Region: (Section-1:Picked)
*Elset, elset=_PickedSet2, internal, generate
  1,  64,   1
** Section: Section-1
*Solid Section, elset=_PickedSet2, material=Material-1
1.,
*End Part
**  
**
** ASSEMBLY
**
*Assembly, name=Assembly
**  
*Instance, name=Part-1-1, part=Part-1
*End Instance
**  
*Nset, nset=_PickedSet5, internal, instance=Part-1-1, generate
  1,  25,   1
*Elset, elset=_PickedSet5, internal, instance=Part-1-1, generate
  1,  16,   1
*Nset, nset=_PickedSet6, internal, instance=Part-1-1
 1,
*Nset, nset=_PickedSet7, internal, instance=Part-1-1, generate
   1,  121,    5
*Elset, elset=_PickedSet7, internal, instance=Part-1-1, generate
  1,  61,   4
*Elset, elset=__PickedSurf4_S4, internal, instance=Part-1-1, generate
  4,  64,   4
*Surface, type=ELEMENT, name=_PickedSurf4, internal
__PickedSurf4_S4, S4
*End Assembly
**  
** MATERIALS
**  
*Material, name=Material-1
*Elastic
196000., 0.3
**  
** BOUNDARY CONDITIONS
**  
** Name: BC-1 Type: Displacement/Rotation
*Boundary
_PickedSet5, 1, 1
** Name: BC-2 Type: Displacement/Rotation
*Boundary
_PickedSet6, 3, 3
** Name: BC-3 Type: Displacement/Rotation
*Boundary
_PickedSet7, 2, 2
** ----------------------------------------------------------------
**  
** STEP: Step-1
**  
*Step, name=Step-1
*Static
1., 1., 1e-05, 1.
**  
** LOADS
**  
** Name: Load-1   Type: Pressure
*Dsload
_PickedSurf4, P, -500.
**  
** OUTPUT REQUESTS
**  
*Restart, write, frequency=0
**  
** FIELD OUTPUT: F-Output-1
**  
*Output, field, variable=PRESELECT
**  
** HISTORY OUTPUT: H-Output-1
**  
*Output, history, variable=PRESELECT
*End Step

Ben tu les mets n'importe où
Enfin non, il faut les mettre dans un step.
Tu les mets où tu veux, mais ils faut qu'il soient entre *STEP et *END STEP
 
Par exemple, avant ton OUTPUT REQUEST
 
Et dans CLOAD, il faut spécifier la direction en 2e position
*CLOAD
1,2,5 : Tu appliques un effort de 5 sur le noeud 1 dans la direction 2
2,2,-4
etc
 
Je crois que le nombre de lignes est limité pour un CLOAD. Il va sans doute falloir mettre plusieurs CLOAD

et voila les valeures de contrainte a chaque noeuds  
noeuds                                contraintes
 
  1 ,                                  -1.000000  
  2,                                    2.267927  
  3                                   -1.000000  
  4                                    276.2610  
  5                                  -1.000000  
  6                                   276.2610  
 

j'ai tjs cette eureur  
***ERROR: A CONCENTRATED LOAD HAS BEEN SPECIFIED ON NODE 0.  THIS NODE IS NOT  
           ACTIVE IN THE MODEL
  *output, field, variable=PRESELECT
  *output, history, variable=PRESELECT
  *endstep
 
 ***ERROR: Concentrated load is specified on 1 nodes that are not active. Check  
           loading definitions in the model. The nodes may have been deleted  
           because they were not connected to any elements.
 

Donc il faut que tu mettes :
*CLOAD
1,2,-1.000000  
2,2, 2.267927  
3,2,-1.000000  
 
Puisque tu m'as dit que ces contraintes étaient sur s22
 
Pour ton erreur, tu as du rentrer le nom d'un noeud qui n'existe pas

ce marche pas tjs . le noeud numero 1 il exite bien sur . mais ce marche pas ............ faut la subroutune sigini

mais je sais pas comment elle marche cette subroutune

Ben je te dis que tu as pas besoin de cette subroutine qui ne marche que pour des élements
Ou alors calcule tes contraintes par éléments.
 
Tu as mal défini Cload ou tes noeuds, c'est pas possible autrement

mon ami ce marche avec cette soubrut et j'ai trouve les bn resultat
merci

J'aimerais bien savoir comment tu as fait alors

voila mon fichier .dat ..... tu peux m'expliquer ce resultat stp  
LINEAR EQUATION SOLVER TYPE         DIRECT SPARSE
 SIGMA=  -3.02831216351296       -3.02831216351296       -3.02831216351296      
  -3.02831216351296       -3.02831216351296       -3.02831216351296      
 SIGMA=  -3.02831216351296       -3.02831216351296       -3.02831216351296      
  -3.02831216351296       -3.02831216351296       -3.02831216351296      
 SIGMA=  -3.02831216351296       -3.02831216351296       -3.02831216351296      
  -3.02831216351296       -3.02831216351296       -3.02831216351296  

Faut que tu regardes comment est défini sigma
Je crois que c'est un truc comme sigma11, sigma 22, sigma33, sigma12 etc
 
Pour être sur de tes résultats, demande une sortie de différents paramètres dans le fichier dat, en utilisant la fonction *NODE PRINT par exemple

tu peux m'expliquer la soubroutune SIGINI ,,... ses parametres stp ????

Ben c'est un peu vague, qu'est-ce que tu comprends pas ?
 
Et au risque de me répèter, je ne crois pas que cette subroutine soit bien adaptée à ton problème.
Il faut dans cette subroutine définir sigma, qui est le tenseur des contraintes au point materiel de calcul (material calculation point),
sachant que ces points n'ont pas l'air d'être les noeuds, sinon tu aurais le noeuds dans les variables que propose sigini.
 
En fait, tu peux savoir de quels points il parle, car tu as une fonction de sigini qui te sort les coordonnées ( COORDS ) sous forme d'un vecteur.
Mais j'ai pas l'impression que ce sont les noeuds, et comme ta fonction est définie sur les noeuds...

j'ai trouve bien les valeurs de ma contrai a chaque neuds mais le pb que  
en t=0 ou t=1 j'ai pas la meme valeurs de contraintes ......
0 -516.538
1 -73.0938
comment je peux faire pour trouver la meme valeure de contrainte ?

C'est quoi t ?

t ce le temps.

j'ai aucun chargement je plastifie pas .... j ai que ma géometrie et j'applelle ma subroutune sigini avec ma fonction et quand je regarde la valeur de s22 a un neud a t = 0 j'ai la contrainte ok mais a t=1 j'ai pas la meme valeur pourtant j'ai pas de chargemnt ... je comprend pas

Si les contraintes varient dans le temps, c'est que ton système n'est pas à l'équilibre. Il faut que tu regardes pourquoi elles varient, et à quel moment.
Pour ca, si tu l'as, tu peux sortir les valeurs de déplacement de contraintes dans le fichier odb pour le regarder avec abaqus viewer

en plus la valeur de s22 en t=0 ce doit etre -516MPa ok ce que j'ai trouve mais a t=1 .j 'ai -73.09 MPa je sais pas ce koi cette valeure

Ben ton calcul est faux alors

j'ai fais un example avec seulement 27 élements je peux t'envoyer ce fichier inp si tu veux

bah oui ce faux mais j'aimerai bien comprendre pourkoi ?

Désolé, mais j'ai pas trop le temos la, il me reste qu'une semaine pour boucler mon stage, alors tout mon temps de calcul alloué est occupé
 
Quand ca marche pas, la premiere chose à faire est de regarder ce qui marche pas. Pour celà, tu peux ecrire dans le fichier .dat les variables de déplacement, d'efforts, etc pour chaque noeud.
Ensuite, tu regardes en détail ce fichier et tu cherches ou ca coince, dans quel incrément, quelles sont les variables en cause etc.
Tu peux aussi sortir un fichier odb et regarder comment évolue ta structure avec abaqus viewer.
Il faut aussi vérifier chaque ligne du fichier inp.
 
Je sais c'est chiant, mais sans connaitre plus en détail ton problème, je peux difficilement faire plus. Je suis désolé de pas pouvoir mieux t'aider    
Mais sache que je suis resté bloqué 1 mois et demi pour UNE ligne de code fausse  

merci et bon courage

Toi aussi, et encore désolé de pas t'avoir aidé plus