Discussion :

[foufi86]

Bonjour

j'ai un programme pour l'identification inverse des paramètres mécaniques que je l'ai utilisé avant et ce marche très bien maintenant je veut l'utiliser mais il marche pas j'ai pas compris ce qu'il arrive. Tous ce que j'ai changé c'est qu'avant j'ai utilisé un essai d'indentation mais maintenant j'ai utilisé un essai de traction charge décharge. Dés que je lance la compilation j'ai cette erreur :

Abaqus JOB Job-9 
Abaqus 6.12-1 
Abaqus License Manager checked out the following licenses: 
Abaqus/Standard checked out 5 tokens. 
<1019 out of 1024 licenses remain available>. 
Begin Analysis Input File Processor 
12/3/2016 9:24:12 AM 
Run pre.exe 
12/3/2016 9:24:20 AM 
Abaqus Error: Analysis Input File Processor exited with an error. 
Abaqus/Analysis exited with errors 
??? Attempted to access c(:,3); index out of bounds because size(c)=[0,0].

Error in ==> fonctapproche1 at 14
[Cmax,J]=max(-c(:,3));

Error in ==> fminsearchbnd>intrafun at 259
fval = feval(params.fun,reshape(xtrans,params.xsize),params.args{:});

Error in ==> fminsearch at 205
fv(:,1) = funfcn(x,varargin{:});

Error in ==> fminsearchbnd at 229
[xu,fval,exitflag,output] = fminsearch(@intrafun,x0u,options,params);

Error in ==> minas at 8
[lam,fval,exitflag,output] = fminsearchbnd(@fonctapproche1,lam0,LB,UB,optio);

y'a t'il quelque'un qui peut m'aider et merci beaucoup

le programme principale est le suivant :

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
30
31
32
33
34
function f = fonctapproche1(lam)
 
global Uexp RFexp
%lecture des données expérimentaux 
fid =fopen('cyclique300%.m','r');
M = fscanf(fid,'%e',[1,1]);
xy=fscanf(fid,'%e%e',[2,M]);
xy1=xy';
Uexp=xy1(:,1) ; %vesteur colonne U2- déplacement 
RFexp=xy1(:,2) ;  %vecteur colonne RF2- réaction
mef(lam);
rf4;
[Umax,I]=max(Uexp);
[Cmax,J]=max(-c(:,3));
ta=size(Uexp);
tb=size(-c(:,3));
Fint1=interp1(Uexp(1:I),RFexp(1:I),-c(1:J,3),'linear');
Fint2=interp1(Uexp(I+1:ta),RFexp(I+1:ta),-c(J:tb,3),'linear');
b=0;
e=0;
Fsim1=-c(1:J,5);
Fsim2=-c(J:tb,5);
tc=size(Fint1,1);
td=size(Fint2,1);
for i=1:tc
    X=(Fint1(i)-Fsim1(i)).^2;
    b=b+X;
end
for j=1:td;
    Y=(Fint2(j)-Fsim2(j)).^2;
    e=e+Y;
end
f=e+b;
plot(-c(J:tb,3),Fint2,-c(1:J,3),Fint1,-c(:,3),-c(:,5));

Merci d'avance

[foufi86]

Rebonjour
voici la fonction qui fait extraire les données de la simulation pour les générer dans un fichier .dat j'ai pas compris ou est le problème exactement une aide svp je me suis bloqué

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
%ouvrir fichier .dat
fid=fopen('Job-15.dat');
i=0;
 c=[];    
 while 1
           i=i+1;
            tline = fgetl(fid);
 
            xy1=fscanf(fid,'%c',[1,16]);
            a=xy1;
            C = strcmp(xy1,'            NOTE');
            y=i;
            if strcmp(xy1,'            NOTE')==1,
               y=i;    
               xy=fscanf(fid,'%12e  ',[1,5]);
               a=xy';
%                react1(1,:)=rx(xy(2:3,:)'); 
  react1=xy(1,:);
 
               c= [c;react1];
            end
 
         if ~ischar(tline), break, end
%             disp(tline);
 
        end
        fclose(fid);

[Jerome Briot]

La variable c n'existe pas car, a priori, il n'y a pas de ligne commençant par ' NOTE' dans le fichier Job-15.dat

Réfère toi à la documentation d'Abaqus. Ce n'est pas un problème concernant MATLAB.

[foufi86]

Merci pour votre réponse

Bon j'ai résolu le problème. Mais ce n'est pas un problème pour ABAQUS car j'utilise Matlab pour lancer ABAQUS. En plus c'est Matlab qui va générer le fichier. dat.
Juste j'ai changé le dossier et ça marche très bien