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