I have developed a program that utilizes parfor to run parallel workers. The structure of the program is as follows:
parfor ix = 1:numel(param_files)
% Read the parameter files
% Run fmincon optimization
end

The program functioned normally with all param_files when executed using local workers. However, when I attempted to run it on a parallel cluster, it became unresponsive for many hours while processing one of the parameter files during the optimization phase. As a result, I had to manually stop the parallel cluster. (Note: The hardware capabilities of the parallel cluster server are equivalent to those of my local PC.)I have developed a program that utilizes parfor to run parallel workers. The structure of the program is as follows:
parfor ix = 1:numel(param_files)
% Read the parameter files
% Run fmincon optimization
end

The program functioned normally with all param_files when executed using local workers. However, when I attempted to run it on a parallel cluster, it became unresponsive for many hours while processing one of the parameter files during the optimization phase. As a result, I had to manually stop the parallel cluster. (Note: The hardware capabilities of the parallel cluster server are equivalent to those of my local PC.) I have developed a program that utilizes parfor to run parallel workers. The structure of the program is as follows:
parfor ix = 1:numel(param_files)
% Read the parameter files
% Run fmincon optimization
end

The program functioned normally with all param_files when executed using local workers. However, when I attempted to run it on a parallel cluster, it became unresponsive for many hours while processing one of the parameter files during the optimization phase. As a result, I had to manually stop the parallel cluster. (Note: The hardware capabilities of the parallel cluster server are equivalent to those of my local PC.) paralell-cluster, parallel computing toolbox, parfor MATLAB Answers — New Questions

​

Hello all,
I hope you are all doing well. I have established a 2DOF differetial equation. Basically, the equations are the 2 spring-mass-damper system. But I am wondering how to define the variable X to be a dimentionless variable, such as X/L. In my equations, the variable X should be dimensionless but when I exclude the force (F_o the dimensionless force), the results of the equation are not dimensionless which I have already checked. My codes are attached. I think you can easily run it.
I appreciate your support.
Best wishes,
Yu
clc;
clear all;

tspan = 0:0.25:200;

X0 = [0 0 0 0];
%parameters——-
mu = 0.02; % mass ratio
f_opt = 1/(1+mu); %frequency ratio
xi_2 = sqrt(3*mu/8*(1+mu));
Omega_1 = 0.188; % 0.03 Hz
Omega_2 = Omega_1*f_opt; % frequency of the second DOF
xi_1 = 0.01;

%matrix——–
M = [1+mu mu;
1 1];
C = [2*xi_1*Omega_1 0;
0 2*xi_2*Omega_2];
K = [Omega_1^2 0;
0 Omega_2^2];

% state space model——————-
O = zeros(2,2);
I = eye(2);
A = [O I; -inv(M)*K -inv(M)*C];
B = [O; inv(M)];
E = [I O];
D = zeros(2,2);

% solve the equations——————-

options = odeset(‘RelTol’,1e-10,’AbsTol’,1e-10);

[t,X] = ode45(@(t,X) QZSdamper(t,A,B,X),tspan,X0,options);

x = X(:,1:4);

% F_o = 2*Omega_2.^2.*alpha.*(sqrt(1-gamma.^2)+X(:,2)).*(-1+sqrt(1/(X(:,2).^2+2.*sqrt(1-gamma^2).*X(:,2)+1)));
%
% f_s = Omega_2.^2*(X(:,2)-2*alpha*(sqrt(1-gamma^2)+X(:,2))*(-1+sqrt(1/(X(:,2).^2+2*sqrt(1-gamma^2)*X(:,2)+1))));
%
% f_1 = f_s(:,1);
%
% figure,
% plot(t,f_1);

figure,
plot(t,x(:,1)); xlabel(‘Time/s’),ylabel(‘Dimensionless displacement of primary structure’)
figure,
plot(t,(x(:,2)))

sys = ss(A,B,E,D);
figure,
bodeplot(sys(1,1)) % from input 1 to output 3
bp.FrequencyScale = "linear";

[mag,phase,wout] = bode(sys(1,1));

mag = squeeze(mag);
phase = squeeze(phase);
fout = wout/(2*pi);
BodeTable = table(fout,mag,phase);

function dXdt = QZSdamper(t,A,B,X)

mu = 0.025;
alpha = 1;
gamma = 2*alpha/(2*alpha+1);

f_opt = 1/(1+mu);
Omega_1 = 0.188; % 0.03 Hz
Omega_2 = Omega_1*f_opt;

w = 0.180;

F = 0.001*sin(w*t);

F_o = 2*Omega_2^2*alpha*(sqrt(1-gamma^2)+X(2))*(sqrt(1/(X(2)^2+2*sqrt(1-gamma^2)*X(2)+1))-1);

% F_o = Omega_2^2*(sqrt(1-gamma^2)+2*alpha*(X(2))*(1/((X(2)-sqrt(1-gamma^2))^2+2*sqrt(1-gamma^2)*(X(2)-sqrt(1-gamma^2))+1)^(1/2)-1))

F = [F;F_o];

dXdt = A*X+B*(F);

endHello all,
I hope you are all doing well. I have established a 2DOF differetial equation. Basically, the equations are the 2 spring-mass-damper system. But I am wondering how to define the variable X to be a dimentionless variable, such as X/L. In my equations, the variable X should be dimensionless but when I exclude the force (F_o the dimensionless force), the results of the equation are not dimensionless which I have already checked. My codes are attached. I think you can easily run it.
I appreciate your support.
Best wishes,
Yu
clc;
clear all;

tspan = 0:0.25:200;

X0 = [0 0 0 0];
%parameters——-
mu = 0.02; % mass ratio
f_opt = 1/(1+mu); %frequency ratio
xi_2 = sqrt(3*mu/8*(1+mu));
Omega_1 = 0.188; % 0.03 Hz
Omega_2 = Omega_1*f_opt; % frequency of the second DOF
xi_1 = 0.01;

%matrix——–
M = [1+mu mu;
1 1];
C = [2*xi_1*Omega_1 0;
0 2*xi_2*Omega_2];
K = [Omega_1^2 0;
0 Omega_2^2];

% state space model——————-
O = zeros(2,2);
I = eye(2);
A = [O I; -inv(M)*K -inv(M)*C];
B = [O; inv(M)];
E = [I O];
D = zeros(2,2);

% solve the equations——————-

options = odeset(‘RelTol’,1e-10,’AbsTol’,1e-10);

[t,X] = ode45(@(t,X) QZSdamper(t,A,B,X),tspan,X0,options);

x = X(:,1:4);

% F_o = 2*Omega_2.^2.*alpha.*(sqrt(1-gamma.^2)+X(:,2)).*(-1+sqrt(1/(X(:,2).^2+2.*sqrt(1-gamma^2).*X(:,2)+1)));
%
% f_s = Omega_2.^2*(X(:,2)-2*alpha*(sqrt(1-gamma^2)+X(:,2))*(-1+sqrt(1/(X(:,2).^2+2*sqrt(1-gamma^2)*X(:,2)+1))));
%
% f_1 = f_s(:,1);
%
% figure,
% plot(t,f_1);

figure,
plot(t,x(:,1)); xlabel(‘Time/s’),ylabel(‘Dimensionless displacement of primary structure’)
figure,
plot(t,(x(:,2)))

sys = ss(A,B,E,D);
figure,
bodeplot(sys(1,1)) % from input 1 to output 3
bp.FrequencyScale = "linear";

[mag,phase,wout] = bode(sys(1,1));

mag = squeeze(mag);
phase = squeeze(phase);
fout = wout/(2*pi);
BodeTable = table(fout,mag,phase);

function dXdt = QZSdamper(t,A,B,X)

mu = 0.025;
alpha = 1;
gamma = 2*alpha/(2*alpha+1);

f_opt = 1/(1+mu);
Omega_1 = 0.188; % 0.03 Hz
Omega_2 = Omega_1*f_opt;

w = 0.180;

F = 0.001*sin(w*t);

F_o = 2*Omega_2^2*alpha*(sqrt(1-gamma^2)+X(2))*(sqrt(1/(X(2)^2+2*sqrt(1-gamma^2)*X(2)+1))-1);

% F_o = Omega_2^2*(sqrt(1-gamma^2)+2*alpha*(X(2))*(1/((X(2)-sqrt(1-gamma^2))^2+2*sqrt(1-gamma^2)*(X(2)-sqrt(1-gamma^2))+1)^(1/2)-1))

F = [F;F_o];

dXdt = A*X+B*(F);

end Hello all,
I hope you are all doing well. I have established a 2DOF differetial equation. Basically, the equations are the 2 spring-mass-damper system. But I am wondering how to define the variable X to be a dimentionless variable, such as X/L. In my equations, the variable X should be dimensionless but when I exclude the force (F_o the dimensionless force), the results of the equation are not dimensionless which I have already checked. My codes are attached. I think you can easily run it.
I appreciate your support.
Best wishes,
Yu
clc;
clear all;

tspan = 0:0.25:200;

X0 = [0 0 0 0];
%parameters——-
mu = 0.02; % mass ratio
f_opt = 1/(1+mu); %frequency ratio
xi_2 = sqrt(3*mu/8*(1+mu));
Omega_1 = 0.188; % 0.03 Hz
Omega_2 = Omega_1*f_opt; % frequency of the second DOF
xi_1 = 0.01;

%matrix——–
M = [1+mu mu;
1 1];
C = [2*xi_1*Omega_1 0;
0 2*xi_2*Omega_2];
K = [Omega_1^2 0;
0 Omega_2^2];

% state space model——————-
O = zeros(2,2);
I = eye(2);
A = [O I; -inv(M)*K -inv(M)*C];
B = [O; inv(M)];
E = [I O];
D = zeros(2,2);

% solve the equations——————-

options = odeset(‘RelTol’,1e-10,’AbsTol’,1e-10);

[t,X] = ode45(@(t,X) QZSdamper(t,A,B,X),tspan,X0,options);

x = X(:,1:4);

% F_o = 2*Omega_2.^2.*alpha.*(sqrt(1-gamma.^2)+X(:,2)).*(-1+sqrt(1/(X(:,2).^2+2.*sqrt(1-gamma^2).*X(:,2)+1)));
%
% f_s = Omega_2.^2*(X(:,2)-2*alpha*(sqrt(1-gamma^2)+X(:,2))*(-1+sqrt(1/(X(:,2).^2+2*sqrt(1-gamma^2)*X(:,2)+1))));
%
% f_1 = f_s(:,1);
%
% figure,
% plot(t,f_1);

figure,
plot(t,x(:,1)); xlabel(‘Time/s’),ylabel(‘Dimensionless displacement of primary structure’)
figure,
plot(t,(x(:,2)))

sys = ss(A,B,E,D);
figure,
bodeplot(sys(1,1)) % from input 1 to output 3
bp.FrequencyScale = "linear";

[mag,phase,wout] = bode(sys(1,1));

mag = squeeze(mag);
phase = squeeze(phase);
fout = wout/(2*pi);
BodeTable = table(fout,mag,phase);

function dXdt = QZSdamper(t,A,B,X)

mu = 0.025;
alpha = 1;
gamma = 2*alpha/(2*alpha+1);

f_opt = 1/(1+mu);
Omega_1 = 0.188; % 0.03 Hz
Omega_2 = Omega_1*f_opt;

w = 0.180;

F = 0.001*sin(w*t);

F_o = 2*Omega_2^2*alpha*(sqrt(1-gamma^2)+X(2))*(sqrt(1/(X(2)^2+2*sqrt(1-gamma^2)*X(2)+1))-1);

% F_o = Omega_2^2*(sqrt(1-gamma^2)+2*alpha*(X(2))*(1/((X(2)-sqrt(1-gamma^2))^2+2*sqrt(1-gamma^2)*(X(2)-sqrt(1-gamma^2))+1)^(1/2)-1))

F = [F;F_o];

dXdt = A*X+B*(F);

end ode45, dimensionless expressions, differential equations MATLAB Answers — New Questions

​

Hi there,
My task is generating a list of all edges in mesh. (ps. I mean not edges of geometry, i saw some function there, which shows edges of a geometry, they are performed bases on interacting between mesh and geometry.)
For a mesh i can easily get access to elements (3d : tetraeder). I am considering break down each element into four faces. It’s easy for me to get edges from faces (2nd step).
(elements -> faces -> edges)
now the problem is in 1st step. A 3d quadratic case, that means an element consists of ten nodes. I have no idea which node is vertices which is in middle of edges.
The way to create mesh :
cubic = [100, 50, 20];
model = createpde(thermal="transient");
geo = multicuboid(cubic(1), cubic(2),…
cubic(3), ‘Zoffset’, 0);
model.Geometry = geo;
generateMesh(model, GeometricOrder=’quadratic’);
pdemesh(model)
Is there any tips for me?
Thank you all.

Best regards
KunHi there,
My task is generating a list of all edges in mesh. (ps. I mean not edges of geometry, i saw some function there, which shows edges of a geometry, they are performed bases on interacting between mesh and geometry.)
For a mesh i can easily get access to elements (3d : tetraeder). I am considering break down each element into four faces. It’s easy for me to get edges from faces (2nd step).
(elements -> faces -> edges)
now the problem is in 1st step. A 3d quadratic case, that means an element consists of ten nodes. I have no idea which node is vertices which is in middle of edges.
The way to create mesh :
cubic = [100, 50, 20];
model = createpde(thermal="transient");
geo = multicuboid(cubic(1), cubic(2),…
cubic(3), ‘Zoffset’, 0);
model.Geometry = geo;
generateMesh(model, GeometricOrder=’quadratic’);
pdemesh(model)
Is there any tips for me?
Thank you all.

Best regards
Kun Hi there,
My task is generating a list of all edges in mesh. (ps. I mean not edges of geometry, i saw some function there, which shows edges of a geometry, they are performed bases on interacting between mesh and geometry.)
For a mesh i can easily get access to elements (3d : tetraeder). I am considering break down each element into four faces. It’s easy for me to get edges from faces (2nd step).
(elements -> faces -> edges)
now the problem is in 1st step. A 3d quadratic case, that means an element consists of ten nodes. I have no idea which node is vertices which is in middle of edges.
The way to create mesh :
cubic = [100, 50, 20];
model = createpde(thermal="transient");
geo = multicuboid(cubic(1), cubic(2),…
cubic(3), ‘Zoffset’, 0);
model.Geometry = geo;
generateMesh(model, GeometricOrder=’quadratic’);
pdemesh(model)
Is there any tips for me?
Thank you all.

Best regards
Kun pde_mesh, edge_of_mesh MATLAB Answers — New Questions

​

Hello everyone,
I’m having problems with the heavy computational load of a Simulink model while running it in real-time. The profiler indicates that the model’s switches are consuming a significant amount of computation time. Does anyone have any recommendations for me, or have they encountered similar issues? If so, please assist me in resolving the issue.
I was wondering if somebody has any alternative model for a MOSFET (or a switch acts as a circuit breaker) that makes the simulations a bit less non-linear. That would be highly appreciated.
Thank youHello everyone,
I’m having problems with the heavy computational load of a Simulink model while running it in real-time. The profiler indicates that the model’s switches are consuming a significant amount of computation time. Does anyone have any recommendations for me, or have they encountered similar issues? If so, please assist me in resolving the issue.
I was wondering if somebody has any alternative model for a MOSFET (or a switch acts as a circuit breaker) that makes the simulations a bit less non-linear. That would be highly appreciated.
Thank you Hello everyone,
I’m having problems with the heavy computational load of a Simulink model while running it in real-time. The profiler indicates that the model’s switches are consuming a significant amount of computation time. Does anyone have any recommendations for me, or have they encountered similar issues? If so, please assist me in resolving the issue.
I was wondering if somebody has any alternative model for a MOSFET (or a switch acts as a circuit breaker) that makes the simulations a bit less non-linear. That would be highly appreciated.
Thank you real-time simulation, switch, mosfet, speedgoat MATLAB Answers — New Questions

​