Fmincon issue: Failure in initial objective function evaluation. FMINCON cannot continue.

3 次查看（过去 30 天）
Cristopher Angola Ararat
移动：Walter Roberson 2023-6-29
I am trying to minimize the absolute differences between frequencies in ETABSV21 through API using a self created objective function. But with the code i have this error is showed: Failure in initial objective function evaluation. FMINCON cannot continue.
Bellow is the function code
function fo=objetivo(x)
global SapModel PropFrame PropMaterial Analyze AnalysisResultsSetup Results
PropMaterial.SetMPIsotropic('conc_losa',x(1)*24855.58,0.2,0.0000099);
PropMaterial.SetMPIsotropic('conc_muros',x(2)*24855.58,0.2,0.0000099);
Analyze.RunAnalysis();
AnalysisResultsSetup.DeselectAllCasesAndCombosForOutput;
AnalysisResultsSetup.SetCaseSelectedForOutput('MODAL');
NumberResults = 0;
LoadCase = {''};
StepType = {''};
StepNum = 0;
Period = 0;
Frequency = 0;
CircFreq = 0;
EigenValue = 0;
[ret,NumberResults,~,~,~,Period,Frequency,~,~]=Results.ModalPeriod(NumberResults, LoadCase, StepType, StepNum, Period, Frequency, CircFreq, EigenValue)
frecuencias=double(Frequency)
fo=mean(abs((frecuencias(1)-2.5781))/2.5781 + abs((frecuencias(3)-2.7930))/2.7930+abs((frecuencias(2)-4.8438))/4.8438);
xo = x;
save('xo', 'xo')
end
NOW this is how i call the function
LI=[20000.58/24855.58];
CI=[24855.58/24855.58];
LS=[33855.58/24855.58];
opciones=optimoptions('fmincon','Algorithm','sqp','Display','iter','PlotFcn',@optimplotfval)
[x,fval,~] = fmincon(@objetivo,CI,[],[],[],[],LI,LS,[],opciones);
  4 个评论
显示 2更早的评论
Torsten
Torsten 2023-6-28
You only have one unknown CI = x(1). Thus refering to x(2) in the line
PropMaterial.SetMPIsotropic('conc_muros',x(2)*24855.58,0.2,0.0000099);
makes no sense.
Cristopher Angola Ararat
Cristopher Angola Ararat 2023-6-28
Thank you so much, however,
I do have this:
LI=[20000.58/24855.58, 20000.58/24855.58];
CI=[24855.58/24855.58, 24855.58/24855.58];
LS=[33855.58/24855.58, 33855.58/24855.58];
Its still the same error

请先登录，再进行评论。

请先登录，再回答此问题。

回答（1 个）

Torsten
Torsten 2023-6-28
移动：Walter Roberson 2023-6-29
It seems that
PropMaterial.SetMPIsotropic
cannot be accessed in "objetivo".
Are you sure it has been written into the global variable "PropMaterial" somewhere before in your code ?
What do you get if you type
PropMaterial
right after
function fo=objetivo(x)
global SapModel PropFrame PropMaterial Analyze AnalysisResultsSetup Results
?
Maybe empty ?

类别

在 Help Center 和 File Exchange 中查找有关 Choose a Solver 的更多信息

标签

产品


版本

R2023a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by