Quote from zcdc0429 on 4. December 2025, 16:53

Hi All,

I have encountered something really strange here.  And I am writing this post to seek some assistance.

I was running the MATLAB script that I wrote successfully on my old Laptop, which has MATLAB R2022b installed. I downloaded the file on my new laptop, which has MATLAB R2025b installed, and the script does not work anymore. I got to successfully entered the simulation loop; however, all values are -1. I check the .sim file I have, and everything seems fine (can be opened by the QBlade.exe and run). Thus, I am really clueless here.

I am using MinGW64 Compiler (version 8.1.0), and I don’t have Visual Studio C++ Compiler installed. Below is my code. Could someone look at it and tell me why? Thanks in advance!

 

%% Helix of NREL5MW in script

clear

close all

addpath(‘.\Functions’);

%clc

%% Define paths

UserPath = ‘.\wakeMixing’;

QBladePath = ‘.\QBladeEE_2.0.9\’; % change

SourcePath = [UserPath ‘Source\’];

DllPath = [QBladePath ‘QBladeEE_2.0.9.dll’]; % change

simFile = [SourcePath ‘NREL5MW_1turbine.sim’];

addpath(‘.\Functions’);

loadlibrary(DllPath,’QBladeLibInclude.h’,’alias’,’QBladeDLL’)

m = libfunctions(‘QBladeDLL’);

if isempty(m)

fprintf(‘Error’)

end

%% Data file (Chage this accordingly)

simTime = 3000; % in timestep, actual time is simTime*timestep(Q-blade define)

timeStep = 0.1; % same with the Q-blade setting

simLen = simTime * timeStep; % seconds

%% Load project and Initialize simulation

%this is setup using relative path and depends on the location of this file

calllib(‘QBladeDLL’,’createInstance’,2,64);

calllib(‘QBladeDLL’,’setLibraryPath’,DllPath);

calllib(‘QBladeDLL’,’loadSimDefinition’,simFile);

calllib(‘QBladeDLL’,’initializeSimulation’)

% Variables

valuestr = ‘Rotational Speed [rpm]’;

valuestr2 = ‘Gen. HSS Torque [Nm]’;

valuestr3 = ‘Tip Speed Ratio [-]’;

Azimu1 = ‘Azimuthal Position Blade 1 [deg]’;

Azimu2 = ‘Azimuthal Position Blade 2 [deg]’;

Azimu3 = ‘Azimuthal Position Blade 3 [deg]’;

Pit1 = ‘Pitch Angle Blade 1 [deg]’;

Pit2 = ‘Pitch Angle Blade 2 [deg]’;

Pit3 = ‘Pitch Angle Blade 3 [deg]’;

PowerVar = ‘Aerodynamic Power [kW]’;

CpVar = ‘Power Coefficient [-]’;

Moop1Var = ‘Aero. OOP RootBend. Mom. Blade 1 [Nm]’;

Mip1Var = ‘Aero. IP RootBend. Mom. Blade 1 [Nm]’;

Moop2Var = ‘Aero. OOP RootBend. Mom. Blade 2 [Nm]’;

Mip2Var = ‘Aero. IP RootBend. Mom. Blade 2 [Nm]’;

Moop3Var = ‘Aero. OOP RootBend. Mom. Blade 3 [Nm]’;

Mip3Var = ‘Aero. IP RootBend. Mom. Blade 3 [Nm]’;

%% Set Turbulent Wind

U_inflow = 10; % Inflow wind speed, same with the Q-blade setting

D_NREL5MW = 126; % Rotor diameter

Hub_NREL5MW = 90; % Hub height

%% Defining Torque Control Setting

% This need to be changed when inflow windspeed is varied

K = 2.24;

N = 97; % Gearbox ratio

%% Simulation

% pre-define array to speed up code

TSR_store = zeros(simTime, 1);

Power_store = zeros(simTime, 1);

% blade 1

Moop1_store = zeros(simTime, 1);

Mip1_store = zeros(simTime, 1);

Mflap1_store = zeros(simTime, 1);

Medge1_store = zeros(simTime, 1);

% blade 2

Moop2_store = zeros(simTime, 1);

Mip2_store = zeros(simTime, 1);

Mflap2_store = zeros(simTime, 1);

Medge2_store = zeros(simTime, 1);

% blade 3

Moop3_store = zeros(simTime, 1);

Mip3_store = zeros(simTime, 1);

Mflap3_store = zeros(simTime, 1);

Medge3_store = zeros(simTime, 1);

Cp_store = zeros(simTime, 1);

%% Simulation

% start simulation

tic

f = waitbar(0,’Initializing Simulation’);

disp(“Initialized”)

for i = 1:1:simTime

calllib(‘QBladeDLL’,’advanceTurbineSimulation’);

% Get current value

omega = calllib(‘QBladeDLL’,’getCustomData_at_num’,valuestr, 0, 0);

genTorqueQB = calllib(‘QBladeDLL’,’getCustomData_at_num’,valuestr2, 0, 0);

TSR = calllib(‘QBladeDLL’,’getCustomData_at_num’,valuestr3, 0, 0);

Azimuth1 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Azimu1, 0, 0);

Azimuth2 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Azimu2, 0, 0);

Azimuth3 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Azimu3, 0, 0);

Pitch1 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Pit1, 0, 0);

Pitch2 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Pit2, 0, 0);

Pitch3 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Pit3, 0, 0);

Power = calllib(‘QBladeDLL’,’getCustomData_at_num’, PowerVar, 0, 0);

Cp = calllib(‘QBladeDLL’,’getCustomData_at_num’, CpVar, 0, 0);

Moop1 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Moop1Var, 0, 0);

Mip1 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Mip1Var, 0, 0);

Moop2 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Moop2Var, 0, 0);

Mip2 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Mip2Var, 0, 0);

Moop3 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Moop3Var, 0, 0);

Mip3 = calllib(‘QBladeDLL’,’getCustomData_at_num’, Mip3Var, 0, 0);

% ==================== Control

% I. Torque control to maintain optimal TSR of 9

omega_g = omega*N; % rotor to generator

genTorque = K.*(omega_g*(2*pi/60))^2;

% II. Send control signal to qblade

calllib(‘QBladeDLL’,’setControlVars_at_num’,[genTorque 0 …

0 0 0],0);

% ==================== Store values

% omega_store(i,:) = omega;

% genTorqueQB_store(i,:) = genTorqueQB;

% genTorque_store(i,:) = genTorque;

TSR_store(i) = TSR;

Power_store(i) = Power;

Cp_store(i) = Cp;

% blade 1

Moop1_store(i) = Moop1;

Mip1_store(i) = Mip1;

Mflap1_store(i) = Moop1*cosd(Pitch1) + Mip1*sind(Pitch1);

Medge1_store(i) = -Moop1*sind(Pitch1) + Mip1*cosd(Pitch1);

% blade 2

Moop2_store(i) = Moop2;

Mip2_store(i) = Mip2;

Mflap2_store(i) = Moop2*cosd(Pitch2) + Mip2*sind(Pitch2);

Medge2_store(i) = -Moop2*sind(Pitch2) + Mip2*cosd(Pitch2);

% blade 3

Moop3_store(i) = Moop3;

Mip3_store(i) = Mip3;

Mflap3_store(i) = Moop3*cosd(Pitch3) + Mip3*sind(Pitch3);

Medge3_store(i) = -Moop3*sind(Pitch3) + Mip3*cosd(Pitch3);

AzimuthAngles(i,:) = [Azimuth1 Azimuth2 Azimuth3];

PitchAngles(i,:) = [Pitch1 Pitch2 Pitch3];

waitbar(i/simTime, f, sprintf(‘Simulation Running: %.1f%%’, (i/simTime)*100));

end

close(f)

calllib(‘QBladeDLL’,’closeInstance’);

toc