How to use Aspen Plus - Matlab Link?

71 次查看（过去 30 天）

Zain Javed 2020-2-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/504450-how-to-use-aspen-plus-matlab-link

评论： Walter Roberson 2024-8-26，18:56

Please help me I am new to Matlab and Aspen Plus but have basic knowledge so please guide in detail. Please tell me what steps do I follow to link Aspen file and Matlab.

The Matlab cod file is:

%% Created by Ing. Andrés Felipe Abril. Universidad Nacional de Colombia. Departamento de Ingeniería Química.
%% Linking
Aspen = actxserver('Apwn.Document.36.0'); %34.0 ---> V8.8; 35.0 ---> V9.0; and 36.0 ---> V10.0
[stat,mess]=fileattrib; % get attributes of folder (Necessary to establish the location of the simulation)
Simulation_Name = 'Reactive_Distillation';% Aspeen Plus Simulation Name
Aspen.invoke('InitFromArchive2',[mess.Name '\' Simulation_Name '.bkp']);
Aspen.Visible = 1; % 1 ---> Aspen is Visible; 0 ---> Aspen is open but not visible
Aspen.SuppressDialogs = 1; % Suppress windows dialogs.
Aspen.Engine.Run2(1); % Run the simulation
while Aspen.Engine.IsRunning == 1 % 1 --> If Aspen is running; 0 ---> If Aspen stop.
    pause(0.5);
end
%% Example of Application
Reflux_Ratio = [1, 3, 5, 7];
for i = 1:length(Reflux_Ratio)
    Aspen.Tree.FindNode("\Data\Blocks\RC-101\Input\BASIS_RR").Value = Reflux_Ratio(i); % Column Reflux
    Aspen.Reinit; % Reinit simulation
    Aspen.Engine.Run2(1); %Run the simulation. (1) ---> Matlab isnt busy; (0) Matlab is Busy;
    time = 1;
    while Aspen.Engine.IsRunning == 1 % 1 --> If Aspen is running; 0 ---> If Aspen stop.
        pause(0.5);
        time = time+1;
        if time==15 % Control of simulation time.
            Aspen.Engine.Stop;
        end
    end
    Simulation_Convergency = Aspen.Tree.FindNode("\Data\Results Summary\Run-Status\Output\PCESSTAT").Value; % 1 Doesn't Convergence; 0 Converge
    if Simulation_Convergency == 0 && time < 10
        Duty(i) = Aspen.Tree.FindNode("\Data\Blocks\RC-101\Output\REB_UTL_DUTY").Value; %Duty value of reactive column
    else
        Duty(i) = inf; % Its Penalized if simulation doesn't converge;
        
    end
end
%% Plotting
figure()
plot(Reflux_Ratio, Duty, 'ok');
Aspen.Close;
Aspen.Quit;

8 个评论
显示 6更早的评论隐藏 6更早的评论

MINH 2024-8-26，10:29

function [Qc, QR, Tc, TR, purity, Dt,diameter ] = retrieve_results(Aspen)

% Initialize output variables

Qc = NaN;

QR = NaN;

Tc = NaN;

TR = NaN;

purity = NaN;

Dt = NaN;

try

% Retrieve nodes

Qc_node = Aspen.Tree.FindNode("\Data\Blocks\C1\Output\COND_DUTY");

QR_node = Aspen.Tree.FindNode("\Data\Blocks\C1\Output\REB_DUTY");

Tc_node = Aspen.Tree.FindNode("\Data\Blocks\C1\Output\TOP_TEMP");

TR_node = Aspen.Tree.FindNode("\Data\Blocks\C1\Output\REB_TOUT");

purity_node = Aspen.Tree.FindNode("\Data\Blocks\C1\Output\COMPONENT\BENZENE\2");

BOTTOM_TEMP_node = Aspen.Tree.FindNode("\Data\Blocks\C1\Output\BOTTOM_TEMP");

FEED_TEMP_node = Aspen.Tree.FindNode("\Data\Blocks\C1\Output\TEMP\1");

diameter_node = Aspen.Tree.FindNode("\Data\Blocks\C1\Input\CA_DIAM\INT-1\CS-1");

% Debugging: Print node paths

fprintf('Qc_node path: %s\n', Qc_node.Name);

fprintf('QR_node path: %s\n', QR_node.Name);

fprintf('Tc_node path: %s\n', Tc_node.Name);

fprintf('TR_node path: %s\n', TR_node.Name);

fprintf('Purity_node path: %s\n', purity_node.Name);

fprintf('BOTTOM_TEMP_node path: %s\n', BOTTOM_TEMP_node.Name);

fprintf('FEED_TEMP_node path: %s\n', FEED_TEMP_node.Name);

fprintf('Diameter_node path: %s\n', diameter_node.Name);

% Check if nodes exist

if isempty(Qc_node)

error('Node for condenser duty not found');

end

if isempty(QR_node)

error('Node for reboiler duty not found');

end

if isempty(Tc_node)

error('Node for condenser temperature not found');

end

if isempty(TR_node)

error('Node for reboiler temperature not found');

end

if isempty(purity_node)

error('Node for top product purity not found');

end

if isempty(BOTTOM_TEMP_node)

error('Node for bottom temperature not found');

end

if isempty(FEED_TEMP_node)

error('Node for feed temperature not found');

end

if isempty(diameter_node)

error('diameter_node not found');

end

% Retrieve and convert node values

Qc = Qc_node.Value * 0.00419; % Convert to kW

QR = QR_node.Value * 0.00419; % Convert to kW

Tc = Tc_node.Value + 273.15; % Convert to Kelvin

TR = TR_node.Value + 273.15; % Convert to Kelvin

purity = purity_node.Value;

BOTTOM_TEMP = BOTTOM_TEMP_node.Value + 273.15; % Convert to Kelvin

FEED_TEMP = FEED_TEMP_node.Value + 273.15; % Convert to Kelvin

diameter= diameter_node.Value;

MINH 2024-8-26，10:30

I send you my diameter retrival code, if you want more please like the picture of process flowsheet in aspen just let me know

Walter Roberson 2024-8-26，18:56

All of the places that you use the Name property are within a try block. You do not show the catch block code.

You should probably put a breakpoint at the assignment to Qc_node and single-step looking at the results of the FindNode calls -- looking to see what size() and class() is being returned, and testing whether the result has an accessible Name property.

请先登录，再进行评论。

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