Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (247)

Search Parameters:
Keywords = Newton–Raphson method

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
22 pages, 1846 KiB  
Article
Estimation of Uncertain Parameters in Single and Double Diode Models of Photovoltaic Panels Using Frilled Lizard Optimization
by Süleyman Dal and Necmettin Sezgin
Electronics 2025, 14(4), 796; https://doi.org/10.3390/electronics14040796 - 18 Feb 2025
Abstract
Renewable energy sources are increasingly crucial for sustainable development. Photovoltaic (PV) systems, which convert solar energy into electricity, offer an environmentally friendly solution. Enhancing energy efficiency and minimizing environmental impacts in these systems heavily rely on parameter optimization. In this study, the Frilled [...] Read more.
Renewable energy sources are increasingly crucial for sustainable development. Photovoltaic (PV) systems, which convert solar energy into electricity, offer an environmentally friendly solution. Enhancing energy efficiency and minimizing environmental impacts in these systems heavily rely on parameter optimization. In this study, the Frilled Lizard Optimization (FLO) algorithm is proposed as a novel approach, integrating the newton-raphson method into the root mean square error (RMSE) objective function process to address nonlinear equations. Extensive analyses conducted on RTC France, STM6-40/36, and Photowatt PWP201 modules demonstrate the superior performance of the FLO algorithm using MATLAB R2022a software. The RMSE values were calculated as 0.0030375 and 0.011538 for SDM and DDM in the RTC France dataset, 0.012036 for the STM6-40/36 dataset and 0.0097545 for the Photowatt-PWP201 dataset, respectively, indicating significantly lower error margins compared to other optimisation methods. Additionally, comprehensive evaluations were carried out using error metrics such as individual absolute error (IAE), relative error (RE) and mean absolute error (MAE), supported by detailed graphical representations of measured and predicted parameters. Current-voltage (I-V) and power-voltage (P-V) characteristic curves, as well as convergence behaviors, were systematically analyzed. This study introduces an innovative and robust solution for parameter optimization in PV systems, contributing to both theoretical and industrial applications. Full article
(This article belongs to the Special Issue Machine Learning Techniques for Image Processing)
29 pages, 11632 KiB  
Article
An Improved Unscented Kalman Filter Applied to Positioning and Navigation of Autonomous Underwater Vehicles
by Jinchao Zhao, Ya Zhang, Shizhong Li, Jiaxuan Wang, Lingling Fang, Luoyin Ning, Jinghao Feng and Jianwu Zhang
Sensors 2025, 25(2), 551; https://doi.org/10.3390/s25020551 - 18 Jan 2025
Viewed by 779
Abstract
To enhance the positioning accuracy of autonomous underwater vehicles (AUVs), a new adaptive filtering algorithm (RHAUKF) is proposed. The most widely used filtering algorithm is the traditional Unscented Kalman Filter or the Adaptive Robust UKF (ARUKF). Excessive noise interference may cause a decrease [...] Read more.
To enhance the positioning accuracy of autonomous underwater vehicles (AUVs), a new adaptive filtering algorithm (RHAUKF) is proposed. The most widely used filtering algorithm is the traditional Unscented Kalman Filter or the Adaptive Robust UKF (ARUKF). Excessive noise interference may cause a decrease in filtering accuracy and is highly likely to result in divergence by means of the traditional Unscented Kalman Filter, resulting in an increase in uncertainty factors during submersible mission execution. An estimation model for system noise, the adaptive Unscented Kalman Filter (UKF) algorithm was derived in light of the maximum likelihood criterion and optimized by applying the rolling-horizon estimation method, using the Newton–Raphson algorithm for the maximum likelihood estimation of noise statistics, and it was verified by simulation experiments using the Lie group inertial navigation error model. The results indicate that, compared with the UKF algorithm and the ARUKF, the improved algorithm reduces attitude angle errors by 45%, speed errors by 44%, and three-dimensional position errors by 47%. It can better cope with complex underwater environments, effectively address the problems of low filtering accuracy and even divergence, and improve the stability of submersibles. Full article
(This article belongs to the Section Navigation and Positioning)
Show Figures

Figure 1

Figure 1
<p>Illustration of the UT Transform: (<b>a</b>) actual; (<b>b</b>) UT Transform.</p>
Full article ">Figure 2
<p>Flowchart of the RHAUKF algorithm.</p>
Full article ">Figure 3
<p>Underwater vehicle trajectories.</p>
Full article ">Figure 4
<p>In the first case, the attitude angle error curves under the three algorithms: (<b>a</b>) heading angle errors; (<b>b</b>) pitch angle errors; (<b>c</b>) roll angle errors.</p>
Full article ">Figure 5
<p>The root mean square errors of the attitude angle parameters: (<b>a</b>) the errors in the time period <span class="html-italic">T</span><sub>1</sub>; (<b>b</b>) the errors in the time period <span class="html-italic">T</span><sub>2</sub>.</p>
Full article ">Figure 6
<p>Deviations of attitude angle parameters: (<b>a</b>) BIAS of heading angle; (<b>b</b>) BIAS of elevation angle; (<b>c</b>) BIAS of roll angle.</p>
Full article ">Figure 7
<p>The attitude angle error curves under the three algorithms in the second case: (<b>a</b>) heading angle errors; (<b>b</b>) pitch angle errors; (<b>c</b>) roll angle errors.</p>
Full article ">Figure 8
<p>The root mean square errors of the attitude angle parameters obtained in <span class="html-italic">T</span><sub>1</sub> and <span class="html-italic">T</span><sub>2</sub>: (<b>a</b>) errors in the time period <span class="html-italic">T</span><sub>1</sub>; (<b>b</b>) errors in the time period <span class="html-italic">T</span><sub>2</sub>.</p>
Full article ">Figure 9
<p>The deviations of the attitude angle parameters obtained by the three algorithms in the second case under <math display="inline"><semantics> <mrow> <msub> <mstyle mathsize="normal"> <mi>T</mi> </mstyle> <mstyle mathsize="normal"> <mn>1</mn> </mstyle> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <msub> <mstyle mathsize="normal"> <mi>T</mi> </mstyle> <mstyle mathsize="normal"> <mn>2</mn> </mstyle> </msub> </mrow> </semantics></math>: (<b>a</b>) BIAS of heading angle; (<b>b</b>) BIAS of elevation angle; (<b>c</b>) BIAS of roll angle.</p>
Full article ">Figure 10
<p>The attitude angle error curves of the three algorithms in the third case: (<b>a</b>) heading angle errors; (<b>b</b>) pitch angle errors; (<b>c</b>) roll angle errors.</p>
Full article ">Figure 11
<p>The root mean square errors of the attitude angle parameters: (<b>a</b>) errors in the <math display="inline"><semantics> <mrow> <msub> <mstyle mathsize="normal"> <mi>T</mi> </mstyle> <mstyle mathsize="normal"> <mn>1</mn> </mstyle> </msub> </mrow> </semantics></math> time period; (<b>b</b>) errors in the <math display="inline"><semantics> <mrow> <msub> <mstyle mathsize="normal"> <mi>T</mi> </mstyle> <mstyle mathsize="normal"> <mn>2</mn> </mstyle> </msub> </mrow> </semantics></math> time period.</p>
Full article ">Figure 12
<p>The deviations of the attitude angle parameters obtained by the three algorithms in the third case: (<b>a</b>) BIAS of heading angle, (<b>b</b>) BIAS of elevation angle, (<b>c</b>) BIAS of roll angle.</p>
Full article ">Figure 13
<p>Three-axis velocity error curves for the six methods: (<b>a</b>) eastward velocity errors; (<b>b</b>) northward velocity errors; (<b>c</b>) vertical velocity errors.</p>
Full article ">Figure 14
<p>The root mean square errors of the three-way velocities: (<b>a</b>) RMSEs of eastward velocities; (<b>b</b>) RMSEs of northward velocities; (<b>c</b>) RMSEs of vertical velocities.</p>
Full article ">Figure 15
<p>The deviations of three-way velocities: (<b>a</b>) the errors in the <math display="inline"><semantics> <mrow> <msub> <mstyle mathsize="normal"> <mi>T</mi> </mstyle> <mstyle mathsize="normal"> <mn>1</mn> </mstyle> </msub> </mrow> </semantics></math> time period; (<b>b</b>) the errors in the <math display="inline"><semantics> <mrow> <msub> <mstyle mathsize="normal"> <mi>T</mi> </mstyle> <mstyle mathsize="normal"> <mn>2</mn> </mstyle> </msub> </mrow> </semantics></math> time period.</p>
Full article ">Figure 16
<p>Three directional position error curves under the three algorithms: (<b>a</b>) eastward position errors; (<b>b</b>) northward position errors; (<b>c</b>) celestial position errors.</p>
Full article ">Figure 17
<p>The root mean square errors of the three-way positions: (<b>a</b>) RMSEs of eastward position; (<b>b</b>) RMSEs of northward position; (<b>c</b>) RMSEs of vertical position.</p>
Full article ">Figure 18
<p>Deviations of three-way positions: (<b>a</b>) deviation curves during the period <math display="inline"><semantics> <mrow> <msub> <mstyle mathsize="normal"> <mi>T</mi> </mstyle> <mstyle mathsize="normal"> <mn>1</mn> </mstyle> </msub> </mrow> </semantics></math>; (<b>b</b>) deviation curves during the period <math display="inline"><semantics> <mrow> <msub> <mstyle mathsize="normal"> <mi>T</mi> </mstyle> <mstyle mathsize="normal"> <mn>2</mn> </mstyle> </msub> </mrow> </semantics></math>.</p>
Full article ">
20 pages, 18408 KiB  
Article
Vibration-Based Damage Prediction in Composite Concrete–Steel Structures Using Finite Elements
by Mario D. Cedeño-Rodríguez, Sergio J. Yanez, Erick I. Saavedra-Flores, Carlos Felipe Guzmán and Juan Carlos Pina
Buildings 2025, 15(2), 200; https://doi.org/10.3390/buildings15020200 - 10 Jan 2025
Viewed by 782
Abstract
The prediction of structural damage through vibrational analysis is a critical task in the field of composite structures. Structural defects and damage can negatively influence the load-carrying capacity of the beam. Therefore, detecting structural damage early is essential to preventing catastrophic failures. This [...] Read more.
The prediction of structural damage through vibrational analysis is a critical task in the field of composite structures. Structural defects and damage can negatively influence the load-carrying capacity of the beam. Therefore, detecting structural damage early is essential to preventing catastrophic failures. This study addresses the challenge of predicting damage in composite concrete–steel beams using a vibration-based finite element approach. To tackle this complex task, a finite element model to a quasi-static analysis emulating a four-point pure bending experimental test was performed. Notably, the numerical model equations were carefully modified using the Newton–Raphson method to account for the stiffness degradation resulting from material strains. These modified equations were subsequently employed in a modal analysis to compute modal shapes and natural frequencies corresponding to the stressed state. The difference between initial and damaged modal shape curvatures served as the foundation for predicting a damage index. The approach effectively captured stiffness degradation in the model, leading to observable changes in modal responses, including a reduction in natural frequencies and variations in modal shapes. This enabled the accurate prediction of damage instances during construction, service, or accidental load scenarios, thereby enhancing the structural and operational safety of composite system designs. This research contributes to the advancement of vibration-based methods for damage detection, emphasizing the complexities in characterizing damage in composite structural geometries. Further exploration and refinement of this approach are essential for the precise classification of damage types. Full article
Show Figures

Figure 1

Figure 1
<p>Scheme of the experimental four-point bending test setup [<a href="#B29-buildings-15-00200" class="html-bibr">29</a>].</p>
Full article ">Figure 2
<p>Composite concrete–steel beam FE model.</p>
Full article ">Figure 3
<p>Concrete slab finite element model.</p>
Full article ">Figure 4
<p>Modeling of the composite beam steel components. (<b>a</b>) Structural steel beam FE mesh. (<b>b</b>) Steel deck FE mesh. (<b>c</b>) Shear stud FE mesh.</p>
Full article ">Figure 5
<p>Menetrey–Willam concrete model stress-strain response.</p>
Full article ">Figure 6
<p>Boundary conditions for the FE composite beam model.</p>
Full article ">Figure 7
<p>Loading protocol.</p>
Full article ">Figure 8
<p>Load-midspan deflection curve obtained from FE simulations compared with experimental test results from Meruane et al. [<a href="#B29-buildings-15-00200" class="html-bibr">29</a>].</p>
Full article ">Figure 9
<p>Contour plots for the modal displacements.</p>
Full article ">Figure 10
<p>Natural frequency reduction percentages.</p>
Full article ">Figure 11
<p>Numerical model MSC and damage index.</p>
Full article ">Figure 12
<p>Numerical damage index.</p>
Full article ">
23 pages, 8821 KiB  
Article
Basins of Convergence in a Multi-Perturbed CR3BP
by Alicia Herrero, Santiago Moll-Lopez, José-A. Moraño, Erika Vega-Fleitas and Daniel Villalibre
Appl. Sci. 2025, 15(1), 106; https://doi.org/10.3390/app15010106 - 26 Dec 2024
Viewed by 562
Abstract
The circular restricted three-body problem (CR3BP) is analyzed to introduce additional factors into the dynamic model, such as radiation forces, flattening of the primary bodies, relativity effects, and the presence of natural satellites. The introduction of these factors increases the accuracy when obtaining [...] Read more.
The circular restricted three-body problem (CR3BP) is analyzed to introduce additional factors into the dynamic model, such as radiation forces, flattening of the primary bodies, relativity effects, and the presence of natural satellites. The introduction of these factors increases the accuracy when obtaining the position of the Lagrange points and the basins of convergence of the system. The Newton–Raphson methodis used to implement a searching algorithm. Finally, an application to the Sun–Mars system including the presence of Phobos and Deimos is developed. Full article
(This article belongs to the Section Aerospace Science and Engineering)
Show Figures

Figure 1

Figure 1
<p>Basins of convergence without radiation forces for the Sun–Earth system (<math display="inline"><semantics> <mrow> <msub> <mi>q</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>q</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>1</mn> <mspace width="0.166667em"/> <mrow> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>μ</mi> <mo>=</mo> <mn>0.05</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mi>n</mi> <mo>=</mo> <mn>0.5</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>).</p>
Full article ">Figure 2
<p>Basins of convergence for the Sun–Earth system obtained for <math display="inline"><semantics> <mrow> <msub> <mi>q</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>0.15</mn> <mspace width="0.166667em"/> <mrow> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>q</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>0.25</mn> <mspace width="0.166667em"/> <mrow> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>μ</mi> <mo>=</mo> <mn>0.5</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mi>n</mi> <mo>=</mo> <mn>0.25</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 3
<p>Basins of convergence for the Sun–Earth system obtained for <math display="inline"><semantics> <mrow> <msub> <mi>q</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>0.15</mn> <mspace width="0.166667em"/> <mrow> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>q</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>0.25</mn> <mspace width="0.166667em"/> <mrow> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>μ</mi> <mo>=</mo> <mn>0.5</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mi>n</mi> <mo>=</mo> <mn>0.95</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>Basins of convergence for the Sun–Earth system considering the effects of oblateness with <math display="inline"><semantics> <mrow> <mi>μ</mi> <mo>=</mo> <mn>0.1</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <msub> <mn>1</mn> <mn>1</mn> </msub> </msub> <mo>=</mo> <mn>0.7</mn> <mspace width="0.166667em"/> <mrow> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <msub> <mn>2</mn> <mn>1</mn> </msub> </msub> <mo>=</mo> <mn>0.5</mn> <mspace width="0.166667em"/> <mrow> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <msub> <mn>1</mn> <mn>2</mn> </msub> </msub> <mo>=</mo> <msub> <mi>σ</mi> <msub> <mn>2</mn> <mn>2</mn> </msub> </msub> <mo>=</mo> <mn>0</mn> <mspace width="0.166667em"/> <mrow> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </mrow> </semantics></math>.</p>
Full article ">Figure 5
<p>Basins of convergence for the Sun–Earth system considering the effects of oblateness with <math display="inline"><semantics> <mrow> <mi>μ</mi> <mo>=</mo> <mn>0.1</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <msub> <mn>1</mn> <mn>1</mn> </msub> </msub> <mo>=</mo> <mn>0.5</mn> <mspace width="0.166667em"/> <mrow> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <msub> <mn>2</mn> <mn>1</mn> </msub> </msub> <mo>=</mo> <mn>0.7</mn> <mspace width="0.166667em"/> <mrow> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <msub> <mi>σ</mi> <msub> <mn>1</mn> <mn>2</mn> </msub> </msub> <mo>=</mo> <msub> <mi>σ</mi> <msub> <mn>2</mn> <mn>2</mn> </msub> </msub> <mo>=</mo> <mn>0</mn> <mspace width="0.166667em"/> <mrow> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </mrow> </semantics></math>.</p>
Full article ">Figure 6
<p>Basins of convergence for the Sun–Earth system considering relativistic effects with <math display="inline"><semantics> <mrow> <mi>μ</mi> <mo>=</mo> <mn>0.5</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mi>ϵ</mi> <mo>=</mo> <mn>0</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 7
<p>Basins of convergence for the Sun–Earth system considering relativistic effects with <math display="inline"><semantics> <mrow> <mi>μ</mi> <mo>=</mo> <mn>0.5</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mi>ϵ</mi> <mo>=</mo> <mn>1</mn> <mspace width="0.166667em"/> <mo>[</mo> <mo>−</mo> <mo>]</mo> </mrow> </semantics></math>.</p>
Full article ">Figure 8
<p>Basins of convergence for the Sun–Mars system in the classical formulation and for a broad size domain.</p>
Full article ">Figure 9
<p>Basins of convergence for the Sun–Mars system in the classical formulation and for a reduced size domain.</p>
Full article ">Figure 10
<p>Basins of convergence for the Sun–Mars system in the classical formulation, zooming in on Mars.</p>
Full article ">Figure 11
<p>Basins of convergence for the Sun–Mars system considering triaxiality effects and a broad size domain.</p>
Full article ">Figure 12
<p>Basins of convergence for the Sun–Mars system considering triaxiality effects and a reduced size domain.</p>
Full article ">Figure 13
<p>Basins of convergence for the Sun–Mars system considering triaxiality effects, zooming in on Mars.</p>
Full article ">Figure 14
<p>Basins of convergence for the Sun–Mars system considering triaxiality and radiation effects for a broad size domain.</p>
Full article ">Figure 15
<p>Basins of convergence for the Sun–Mars system considering triaxiality and radiation effects for a reduced size domain.</p>
Full article ">Figure 16
<p>Basins of convergence for the Sun–Mars system considering triaxiality and radiation effects, zooming in on Mars.</p>
Full article ">Figure 17
<p>Iterations needed until convergence for the Sun–Mars system when triaxiality effects are included [broad domain].</p>
Full article ">Figure 18
<p>Iterations needed until convergence for the Sun–Mars system when triaxiality effects are included [reduced domain].</p>
Full article ">Figure 19
<p>Iterations needed until convergence for the Sun–Mars system when triaxiality effects are included [zooming in on Mars].</p>
Full article ">
39 pages, 3488 KiB  
Article
Parameter Extraction for Photovoltaic Models with Flood-Algorithm-Based Optimization
by Yacine Bouali and Basem Alamri
Mathematics 2025, 13(1), 19; https://doi.org/10.3390/math13010019 - 25 Dec 2024
Viewed by 795
Abstract
Accurately modeling photovoltaic (PV) cells is crucial for optimizing PV systems. Researchers have proposed numerous mathematical models of PV cells to facilitate the design and simulation of PV systems. Usually, a PV cell is modeled by equivalent electrical circuit models with specific parameters, [...] Read more.
Accurately modeling photovoltaic (PV) cells is crucial for optimizing PV systems. Researchers have proposed numerous mathematical models of PV cells to facilitate the design and simulation of PV systems. Usually, a PV cell is modeled by equivalent electrical circuit models with specific parameters, which are often unknown; this leads to formulating an optimization problem that is addressed through metaheuristic algorithms to identify the PV cell/module parameters accurately. This paper introduces the flood algorithm (FLA), a novel and efficient optimization approach, to extract parameters for various PV models, including single-diode, double-diode, and three-diode models and PV module configurations. The FLA’s performance is systematically evaluated against nine recently developed optimization algorithms through comprehensive comparative and statistical analyses. The results highlight the FLA’s superior convergence speed, global search capability, and robustness. This study explores two distinct objective functions to enhance accuracy: one based on experimental current–voltage data and another integrating the Newton–Raphson method. Applying metaheuristic algorithms with the Newton–Raphson-based objective function reduced the root-mean-square error (RMSE) more effectively than traditional methods. These findings establish the FLA as a computationally efficient and reliable approach to PV parameter extraction, with promising implications for advancing PV system design and simulation. Full article
(This article belongs to the Special Issue Applied Mathematics and Intelligent Control in Electrical Engineering)
Show Figures

Figure 1

Figure 1
<p>Equivalent circuit of single-diode model.</p>
Full article ">Figure 2
<p>Equivalent circuit of double-diode model.</p>
Full article ">Figure 3
<p>Equivalent circuit of three-diode model.</p>
Full article ">Figure 4
<p>Equivalent circuit of PV module model.</p>
Full article ">Figure 5
<p>Concept of objective function <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math> calculation.</p>
Full article ">Figure 6
<p>Concept of objective function <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math> calculation.</p>
Full article ">Figure 7
<p>The flowchart of FLA.</p>
Full article ">Figure 8
<p>RMSE evolution of different algorithms for the SDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 9
<p>RMSE evolution of different algorithms for the SDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 10
<p>Average CPU times of different algorithms for the SDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 11
<p>RMSE evolution of the FLA for the SDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 12
<p>RMSE evolution of different algorithms for the DDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 13
<p>RMSE evolution of different algorithms for the DDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 14
<p>Average CPU times of different algorithms for the DDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 15
<p>RMSE evolution of the FLA for the DDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 16
<p>RMSE evolution of different algorithms for the TDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 17
<p>RMSE evolution of different algorithms for the TDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 18
<p>Average CPU times of different algorithms for the TDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 19
<p>RMSE evolution of the FLA for the TDM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 20
<p>The I-V curves of the FLA for the (<b>a</b>) SDM, (<b>b</b>) DDM, and (<b>c</b>) TDM of RTC France, using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 21
<p>The P-V curves of the FLA for the (<b>a</b>) SDM, (<b>b</b>) DDM, and (<b>c</b>) TDM of RTC France, using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 22
<p>RMSE evolution of different algorithms for the PVM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 23
<p>RMSE evolution of different algorithms for the PVM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 24
<p>RMSE evolution of the FLA for the PVM using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 25
<p>The relative error values of the simulated current data and the experimental current data using the FLA for Photowatt-PWP201.</p>
Full article ">Figure 26
<p>The (<b>a</b>) I-V and (<b>b</b>) P-V curves of the FLA for Photowatt-PWP201, using <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>1</mn> </msub> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>O</mi> <mi>b</mi> <msub> <mi>j</mi> <mn>2</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">
23 pages, 5375 KiB  
Article
Power Flow Analysis of Ring AC/DC Hybrid Network with Multiple Power Electronic Transformers Based on Hybrid Alternating Iteration Power Flow Algorithm
by Zhen Zheng, Chenhong Huang, Xiaoli Ma, Wenwen Chen, Yinan Huang, Min Wang and Dongqian Pan
Processes 2025, 13(1), 7; https://doi.org/10.3390/pr13010007 - 24 Dec 2024
Viewed by 410
Abstract
AC/DC hybrid distribution networks with power electronic transformers (PETs) as distribution hubs are in line with the future development direction of the AC/DC hybrid distribution network. Unlike traditional transformers, power electronic transformers introduce new node types and may transform the network topology from [...] Read more.
AC/DC hybrid distribution networks with power electronic transformers (PETs) as distribution hubs are in line with the future development direction of the AC/DC hybrid distribution network. Unlike traditional transformers, power electronic transformers introduce new node types and may transform the network topology from radial to ring structures. These changes render traditional power flow calculation methods inadequate for achieving satisfactory results in AC/DC hybrid networks. In addition, existing commercial power flow calculation software packages are mainly based on the traditional AC power flow calculation method, which have limited support for the DC network. Especially when the DC network is coupled with the AC network, it is difficult to achieve a unified calculation of its power flow. To address these challenges, this paper proposes a novel power flow calculation method for ring AC/DC hybrid distribution networks with power electronic transformers. The proposed method is based on the alternating iterative method to ensure compatibility with mature AC power flow calculation programs in commercial software, thereby improving the feasibility of engineering applications. Firstly, the steady-state power flow calculation model of PET is constructed by analyzing that the working principle and control modes of power electronic transformer are proposed based on the source-load attributes of its connected subnetworks. According to the characteristics of the power electronic transformer, AC distribution network, and DC distribution network, a hybrid alternating iteration method combining the high computational accuracy of the Newton–Raphson (NR) method with the high efficiency of the Zbus Gaussian method in dealing with ring networks is proposed. On this basis, the power flow calculation model of the AC/DC hybrid distribution network with power electronic transformers is established. Finally, the simulation of the constructed 44-node ring AC/DC hybrid distribution network example is carried out. The simulation results show that the proposed method can not only converge reliably when the convergence accuracy is 1 × 10−6 p.u., but also ensure that the voltage magnitudes of all nodes are above 0.96 p.u. whose maximum offset value is 0.789% when the outputs of the connected distributed generations fluctuate, which verifies the effectiveness and accuracy of the proposed method. Full article
(This article belongs to the Section Energy Systems)
Show Figures

Figure 1

Figure 1
<p>Traditional AC distribution networks without DC bus.</p>
Full article ">Figure 2
<p>AC/DC hybrid distribution networks with non-isolated bidirectional converters.</p>
Full article ">Figure 3
<p>AC/DC hybrid distribution networks with PETs.</p>
Full article ">Figure 4
<p>Physical structure of PET.</p>
Full article ">Figure 5
<p>Model of the AC port based on VSC.</p>
Full article ">Figure 6
<p>Model of the DC port based on the DC/DC converter.</p>
Full article ">Figure 7
<p>A flowchart of the power flow calculation process for AC/DC hybrid distribution networks with PETs.</p>
Full article ">Figure 8
<p>Improved IEEE 33-bus system connected to the DC system.</p>
Full article ">Figure 9
<p>Partition results.</p>
Full article ">Figure 10
<p>Active power of all ports of two PETs.</p>
Full article ">Figure 11
<p>Voltage magnitude of all ports of two PETs.</p>
Full article ">Figure 12
<p>Modulation ratio of AC ports of two PETs.</p>
Full article ">Figure 13
<p>Outputs of PVs and WTs in three cases.</p>
Full article ">Figure 14
<p>Comparisons of the bus voltage magnitudes between three cases.</p>
Full article ">Figure 15
<p>Bus voltage offset distribution of case 2 and case 3.</p>
Full article ">Figure 16
<p>Comparisons of the active power of all ports of two PETs between three cases.</p>
Full article ">
20 pages, 4161 KiB  
Article
Transformer–BiLSTM Fusion Neural Network for Short-Term PV Output Prediction Based on NRBO Algorithm and VMD
by Xiaowei Fan, Ruimiao Wang, Yi Yang and Jingang Wang
Appl. Sci. 2024, 14(24), 11991; https://doi.org/10.3390/app142411991 - 21 Dec 2024
Viewed by 699
Abstract
In order to solve the difficulties that the uncertain characteristics of PV output, such as volatility and intermittency, will bring to the development of microgrid scheduling plans, this paper proposes a Transformer–Bidirectional Long Short-Term Memory (BiLSTM) neural network PV power generation forecasting fusion [...] Read more.
In order to solve the difficulties that the uncertain characteristics of PV output, such as volatility and intermittency, will bring to the development of microgrid scheduling plans, this paper proposes a Transformer–Bidirectional Long Short-Term Memory (BiLSTM) neural network PV power generation forecasting fusion model based on the Newton–Raphson optimization algorithm (NRBO) and Variational Modal Decomposition (VMD). Firstly, the principle of the VMD technique and the Gray Wolf Optimization (GWO) algorithm’s key parameter optimization method for VMD are introduced. Then, the Transformer decoder partially fuses the BiLSTM network and retains the encoder to obtain the body of the prediction model, followed by explaining the principle of the NRBO algorithm. And finally, the VMD-NRBO-Transformer-BiLSTM prediction model and hyperparameter selection are evaluated by the NRBO algorithm. The algorithm sets up a multi-model comparison experiment, and the results show that the prediction model proposed in this paper has the best prediction accuracy and the optimal evaluation index. Full article
(This article belongs to the Section Energy Science and Technology)
Show Figures

Figure 1

Figure 1
<p>Flowchart of GWO-optimized VMD parameters.</p>
Full article ">Figure 2
<p>The structure of the Transformer model.</p>
Full article ">Figure 3
<p>The structure of BiLSTM.</p>
Full article ">Figure 4
<p>The structure of the Transformer–BiLSTM model.</p>
Full article ">Figure 5
<p>VMD-NRBO-Transformer-BiLSTM overall framework.</p>
Full article ">Figure 6
<p>Partial unprocessed original input data.</p>
Full article ">Figure 7
<p>Partial normalization of the original input data.</p>
Full article ">Figure 8
<p>Waveforms of each sequence after VMD.</p>
Full article ">Figure 9
<p>The iteration curve of the NRBO algorithm.</p>
Full article ">Figure 10
<p>Comparison of training performance of four prediction models.</p>
Full article ">Figure 11
<p>Comparison of testing performance of four prediction models.</p>
Full article ">Figure 12
<p>Prediction errors of each model on the test set. (<b>a</b>) Prediction errors of the Transformer model on the test set; (<b>b</b>) prediction errors of the Transformer–BiLSTM model on the test set; (<b>c</b>) prediction errors of the VMD-Transformer-BiLSTM model on the test set; and (<b>d</b>) prediction errors of the VMD-NRBO-Transformer-BiLSTM model on the test set.</p>
Full article ">Figure 13
<p>Compass diagrams of error evaluation metrics for each model. (<b>a</b>) Compass plot of MAE for each model; (<b>b</b>) compass plot of MAPE for each model; (<b>c</b>) compass plot of MSE for each model; (<b>d</b>) compass plot of RMSE for each model; and (<b>e</b>) compass plot of 1-R<sup>2</sup> for each model.</p>
Full article ">
16 pages, 1360 KiB  
Article
Modeling and Analysis of Thermoelastic Damping in a Piezoelectro-Magneto-Thermoelastic Imperfect Flexible Beam
by Ayman M. Alneamy, Sayantan Guha and Mohammed Y. Tharwan
Mathematics 2024, 12(24), 4011; https://doi.org/10.3390/math12244011 - 20 Dec 2024
Viewed by 1069
Abstract
This research addresses the phenomena of thermoelastic damping (TED) and frequency shift (FS) of a thin flexible piezoelectro-magneto-thermoelastic (PEMT) composite beam. Its motion is constrained by two linear flexible springs attached to both ends. The novelty behind the proposed study is to mimic [...] Read more.
This research addresses the phenomena of thermoelastic damping (TED) and frequency shift (FS) of a thin flexible piezoelectro-magneto-thermoelastic (PEMT) composite beam. Its motion is constrained by two linear flexible springs attached to both ends. The novelty behind the proposed study is to mimic the uncertainties during the fabrication of the beam. Therefore, the equation of motion was derived utilizing the linear Euler–Bernoulli theory accounting for the flexible boundary conditions. The beam’s eigenvalues, mode shapes, and the effects of the thermal relaxation time (t1), the dimensions of the beam, the linear spring coefficients (KL0 and KLL), and the critical thickness (CT) on both TED and FS of the PEMT beam were investigated numerically employing the Newton–Raphson method. The results show that the peak value of thermoelastic damping (Qpeak1) and the frequency shift (Ω) of the beam increase as t1 escalates. Another observation was made for the primary fundamental mode, where an increase in the spring coefficient KLL leads to a further increase in Ω. On the other hand, the opposite trend is noted for the higher modes. Indeed, the results show the possibility of using the proposed design in a variety of applications that involve damping dissipation. Full article
(This article belongs to the Special Issue Mathematical and Computational Methods for Mechanics and Engineering)
Show Figures

Figure 1

Figure 1
<p>Schematics showing (<b>a</b>) the proposed PTFRC composite beam and (<b>b</b>) the beam connected to a flexible supports at both ends and its parameters.</p>
Full article ">Figure 2
<p>Impact of <math display="inline"><semantics> <mrow> <mi>K</mi> <mi>L</mi> <mi>L</mi> </mrow> </semantics></math> on the scaled <math display="inline"><semantics> <msup> <mi>Q</mi> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> </semantics></math> varying against <span class="html-italic">H</span> for fixed <math display="inline"><semantics> <mrow> <mi>L</mi> <mo>=</mo> <mn>1000</mn> </mrow> </semantics></math> µm, <math display="inline"><semantics> <mrow> <mi>K</mi> <mi>L</mi> <mn mathvariant="italic">0</mn> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>8</mn> <msub> <mi>t</mi> <mn>0</mn> </msub> </mrow> </semantics></math> for (<b>a</b>) mode 1, (<b>b</b>) mode 2, (<b>c</b>) mode 3, (<b>d</b>) mode 4, and (<b>e</b>) mode 5.</p>
Full article ">Figure 3
<p>Scaled <math display="inline"><semantics> <msup> <mi>Q</mi> <mrow> <mo>−</mo> <mn>1</mn> </mrow> </msup> </semantics></math> varying against <span class="html-italic">H</span> for fixed beam length <math display="inline"><semantics> <mrow> <mi>L</mi> <mo>=</mo> <mn>1000</mn> </mrow> </semantics></math> µm, left linear spring coefficient <math display="inline"><semantics> <mrow> <mi>K</mi> <mi>L</mi> <mn mathvariant="italic">0</mn> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math>, and right linear spring coefficient <math display="inline"><semantics> <mrow> <mi>K</mi> <mi>L</mi> <mi>L</mi> <mo>=</mo> <mn>2</mn> </mrow> </semantics></math> for different values of the time relaxation (<b>a</b>) <math display="inline"><semantics> <mrow> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>4</mn> <msub> <mi>t</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, (<b>b</b>) <math display="inline"><semantics> <mrow> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>8</mn> <msub> <mi>t</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, and (<b>c</b>) <math display="inline"><semantics> <mrow> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>12</mn> <msub> <mi>t</mi> <mn>0</mn> </msub> </mrow> </semantics></math>.</p>
Full article ">Figure 4
<p>The impact of <math display="inline"><semantics> <mrow> <mi>K</mi> <mi>L</mi> <mi>L</mi> </mrow> </semantics></math> on the frequency <math display="inline"><semantics> <mo>Ω</mo> </semantics></math> while varying the beam thickness <span class="html-italic">H</span> and taking into account a fixed <math display="inline"><semantics> <mrow> <mi>A</mi> <mi>R</mi> <mo>=</mo> <mn>100</mn> </mrow> </semantics></math>, <math display="inline"><semantics> <mrow> <msub> <mi>t</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>8</mn> <msub> <mi>t</mi> <mn>0</mn> </msub> </mrow> </semantics></math>, and <math display="inline"><semantics> <mrow> <mi>K</mi> <mi>L</mi> <mn mathvariant="italic">0</mn> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math> for a variational mode of (<b>a</b>) mode 1, (<b>b</b>) mode 2, (<b>c</b>) mode 3, (<b>d</b>) mode 4, and (<b>e</b>) mode 5.</p>
Full article ">Figure 5
<p>The influence of the time relaxiation <math display="inline"><semantics> <msub> <mi>t</mi> <mn>1</mn> </msub> </semantics></math> on the frequency <math display="inline"><semantics> <mo>Ω</mo> </semantics></math> while varying the beam thickness <span class="html-italic">H</span> and taking into account a fixed <math display="inline"><semantics> <mrow> <mi>A</mi> <mi>R</mi> <mo>=</mo> <mn>100</mn> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mi>K</mi> <mi>L</mi> <mn mathvariant="italic">0</mn> <mo>=</mo> <mi>K</mi> <mi>L</mi> <mi>L</mi> <mo>=</mo> <mn>0.5</mn> </mrow> </semantics></math> for a variational mode of (<b>a</b>) mode 1, (<b>b</b>) mode 2, (<b>c</b>) mode 3, (<b>d</b>) mode 4, and (<b>e</b>) mode 5.</p>
Full article ">
18 pages, 5056 KiB  
Article
Accurate Dynamic Analysis Method of Cable-Damper System Based on Dynamic Stiffness Method
by Hui Jiao, Bin Xu, Zhengkai Jiang, Can Cui and Haoxiang Yang
Buildings 2024, 14(12), 4007; https://doi.org/10.3390/buildings14124007 - 17 Dec 2024
Viewed by 493
Abstract
To suppress large vibrations of the cable in cable-stayed bridges, it is common to install transverse dampers near the end of the cable. This paper focuses on the cable-damper system; based on the dynamic stiffness method, an accurate dynamic analysis method considering cable [...] Read more.
To suppress large vibrations of the cable in cable-stayed bridges, it is common to install transverse dampers near the end of the cable. This paper focuses on the cable-damper system; based on the dynamic stiffness method, an accurate dynamic analysis method considering cable parameters, damper parameters, and cable forces is proposed. First, a mechanical analysis model is established which is closer to the cable with a transverse damper installed in the bridge. The model considers the cable bending stiffness, sag, inclination angle, cable force, damper stiffness, damping coefficient, and damper installation height. Then, the characteristic frequency equation of the cable-damper system is established, and a solution method that combines the initial value method and Newton–Raphson method is proposed. This method is confirmed to provide more accurate frequency analysis for the cable-damper system. Finally, using this method, the effect of the damper parameters on the dynamic characteristics of the system is investigated. Full article
(This article belongs to the Special Issue Advances and Applications in Structural Vibration Control)
Show Figures

Figure 1

Figure 1
<p>Cable-damper system mechanical model.</p>
Full article ">Figure 2
<p>Force boundary condition and dynamic displacement boundary conditions of the cable.</p>
Full article ">Figure 3
<p>Variation curve of characteristic frequency function of cable-damper system with frequency.</p>
Full article ">Figure 4
<p>Schematic diagram of experimental setting and photograph of experimental setting [<a href="#B21-buildings-14-04007" class="html-bibr">21</a>].</p>
Full article ">Figure 5
<p>Shaozhou Bridge site photos.</p>
Full article ">Figure 6
<p>The overall layout of Shaozhou Bridge and external damper installation diagram.</p>
Full article ">Figure 7
<p>Acceleration history and spectrogram of Shaozhou Bridge cables.</p>
Full article ">Figure 8
<p>Influence of damper stiffness and damping coefficient on cable-damper system frequency.</p>
Full article ">Figure 9
<p>Influence of damping coefficient on cable frequency.</p>
Full article ">Figure 10
<p>Influence of damper stiffness on cable frequency.</p>
Full article ">Figure 11
<p>Effect of damper installation height on cable frequency.</p>
Full article ">
17 pages, 1673 KiB  
Article
Nonlinear Thermomechanical Low-Velocity Impact Behaviors of Geometrically Imperfect GRC Beams
by Tao Zhang, Qiang Li, Jia-Jia Mao and Chunqing Zha
Materials 2024, 17(24), 6062; https://doi.org/10.3390/ma17246062 - 11 Dec 2024
Viewed by 520
Abstract
This paper studies the thermomechanical low-velocity impact behaviors of geometrically imperfect nanoplatelet-reinforced composite (GRC) beams considering the von Kármán nonlinear geometric relationship. The graphene nanoplatelets (GPLs) are assumed to have a functionally graded (FG) distribution in the matrix beam along its thickness, following [...] Read more.
This paper studies the thermomechanical low-velocity impact behaviors of geometrically imperfect nanoplatelet-reinforced composite (GRC) beams considering the von Kármán nonlinear geometric relationship. The graphene nanoplatelets (GPLs) are assumed to have a functionally graded (FG) distribution in the matrix beam along its thickness, following the X-pattern. The Halpin–Tsai model and the rule of mixture are employed to predict the effective Young modulus and other material properties. Dividing the impact process into two stages, the corresponding impact forces are calculated using the modified nonlinear Hertz contact law. The nonlinear governing equations are obtained by introducing the von Kármán nonlinear displacement–strain relationship into the first-order shear deformation theory and dispersed via the differential quadrature (DQ) method. Combining the governing equation of the impactor’s motion, they are further parametrically solved by the Newmark-β method associated with the Newton–Raphson iterative process. The influence of different types of geometrical imperfections on the nonlinear thermomechanical low-velocity impact behaviors of GRC beams with varying weight fractions of GPLs, subjected to different initial impact velocities, are studied in detail. Full article
(This article belongs to the Special Issue Functionally Graded Graphene Nanocomposite Materials and Structures)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) <span class="html-italic">N</span>-layered geometrically imperfect GRC beam subjected to low-velocity impact and (<b>b</b>) the schematic diagram of the functionally graded distributed GPL along the thickness of the beam.</p>
Full article ">Figure 2
<p>Types of imperfection: (<b>a</b>) sine, (<b>b</b>) global, and (<b>c</b>) local.</p>
Full article ">Figure 3
<p>Comparison results: low-velocity impact force history of the intact CNTRC beam with FG-X pattern [<a href="#B52-materials-17-06062" class="html-bibr">52</a>].</p>
Full article ">Figure 4
<p>Influences of the imperfect mode and amplitude on the (<b>a</b>) time-dependent impact force and (<b>b</b>) nonlinear thermomechanical central deflection of a geometrically imperfect GRC beam.</p>
Full article ">Figure 5
<p>Effects of temperature variation on the (<b>a</b>) time-dependent impact force and (<b>b</b>) nonlinear thermomechanical central deflection of a GRC beam with different geometrical imperfections.</p>
Full article ">Figure 6
<p>Effects of the initial impact velocity of the impactor on the (<b>a</b>) time-dependent impact force and (<b>b</b>) nonlinear thermomechanical central deflection of a GRC beam with different geometrical imperfections.</p>
Full article ">Figure 7
<p>Effects of the weight fraction of GPLs and imperfect amplitude on the (<b>a</b>) time-dependent impact force and (<b>b</b>) nonlinear thermomechanical central deflection of a GRC beam with global imperfections.</p>
Full article ">Figure 8
<p>Effects of the weight fraction of GPLs and imperfect amplitude on the (<b>a</b>) time-dependent impact force and (<b>b</b>) nonlinear thermomechanical central deflection of a GRC beam with sine imperfections.</p>
Full article ">Figure 9
<p>Effects of the weight fraction of GPLs and imperfect amplitude on the (<b>a</b>) time-dependent impact force and (<b>b</b>) nonlinear thermomechanical central deflection of a GRC beam with local imperfections.</p>
Full article ">
23 pages, 5711 KiB  
Article
Modeling and Dynamic Analysis of Double-Row Angular Contact Ball Bearing–Rotor–Disk System
by Haibiao Zhang, Zhen Li, Haijian Liu, Tao Liu and Qingshan Wang
Lubricants 2024, 12(12), 441; https://doi.org/10.3390/lubricants12120441 - 10 Dec 2024
Viewed by 778
Abstract
This article presents a general numerical method to establish a mathematical model of a bearing–rotor–disk system. This mathematical model consists of two double-row angular contact ball bearings (DRACBBs), a rotor and a rigid disk. The mathematical model of the DRACBB is built on [...] Read more.
This article presents a general numerical method to establish a mathematical model of a bearing–rotor–disk system. This mathematical model consists of two double-row angular contact ball bearings (DRACBBs), a rotor and a rigid disk. The mathematical model of the DRACBB is built on the basis of elastic Hertz contact by adopting the Newton Raphson iteration method, and three different structure forms are taken into account. The rotor is modeled by employing a finite element method in conjunction with Timoshenko beam theory, and the rigid disk is modeled by applying the lumped parameter method. The mathematical model of the bearing–rotor–disk system is constructed by the coupling of the bearing, rotor and disk, and the dynamic response of the bearing–rotor–disk system can be solved by employing the Newmark-β method. The validation of the above mathematical model is demonstrated by comparing the proposed results with the results from the existing literature and finite element software. The dynamic characteristics of the DRACBBs and the dynamic response of the bearing–rotor–disk system are investigated by parametric study. A dynamic characteristic analysis of the DRACBB is conducted to ensure the optimal structure form of the DRACBB under complex external loads, and it can provide a reference for the selection of the structural forms of DRACBBs. Full article
Show Figures

Figure 1

Figure 1
<p>The structural forms of DRACBBs.</p>
Full article ">Figure 2
<p>The coordinate systems of a DRACBB.</p>
Full article ">Figure 3
<p>The geometric location relationship of a DRACBB.</p>
Full article ">Figure 4
<p>The geometric description of a two-node beam element.</p>
Full article ">Figure 5
<p>The geometric description of the disk.</p>
Full article ">Figure 6
<p>The dynamic model of the bearing–rotor–disk system.</p>
Full article ">Figure 7
<p>The comparison results of the contact force of DRACBBs.</p>
Full article ">Figure 8
<p>The contact characteristics of DRACBBs under purely axial force <span class="html-italic">F<sub>z</sub></span>.</p>
Full article ">Figure 9
<p>The dynamic stiffness of DRACBBs under purely axial force <span class="html-italic">F<sub>z</sub></span>.</p>
Full article ">Figure 10
<p>The contact characteristics of DRACBBs under purely radial force <span class="html-italic">F<sub>x</sub></span> and <span class="html-italic">F<sub>y</sub></span>.</p>
Full article ">Figure 11
<p>The dynamic stiffness of DRACBBs under radial force <span class="html-italic">F<sub>x</sub></span> and <span class="html-italic">F<sub>y</sub></span>.</p>
Full article ">Figure 12
<p>The contact characteristics of DRACBBs under torque <span class="html-italic">M<sub>x</sub></span> and <span class="html-italic">M<sub>y</sub></span>.</p>
Full article ">Figure 13
<p>The dynamic stiffness of DRACBBs under torque <span class="html-italic">M<sub>x</sub></span> and <span class="html-italic">M<sub>y</sub></span>.</p>
Full article ">Figure 14
<p>The dynamic responses of the bearing–rotor–disk system under unbalanced excitation.</p>
Full article ">Figure 14 Cont.
<p>The dynamic responses of the bearing–rotor–disk system under unbalanced excitation.</p>
Full article ">Figure 15
<p>FFT of the dynamic response of the bearing–rotor–disk system.</p>
Full article ">Figure 15 Cont.
<p>FFT of the dynamic response of the bearing–rotor–disk system.</p>
Full article ">
14 pages, 730 KiB  
Article
Fired Heaters Optimization by Estimating Real-Time Combustion Products Using Numerical Methods
by Ricardo Sánchez, Argemiro Palencia-Díaz, Jonathan Fábregas-Villegas and Wilmer Velilla-Díaz
Energies 2024, 17(23), 6190; https://doi.org/10.3390/en17236190 - 9 Dec 2024
Viewed by 690
Abstract
Fired heaters upstream of distillation towers, despite their optimal thermal efficiency, often suffer from performance decline due to fluctuations in fuel composition and unpredictable operational parameters. These heaters have high energy consumption, as fuel properties vary depending on the source of the crude [...] Read more.
Fired heaters upstream of distillation towers, despite their optimal thermal efficiency, often suffer from performance decline due to fluctuations in fuel composition and unpredictable operational parameters. These heaters have high energy consumption, as fuel properties vary depending on the source of the crude oil. This study aims to optimize the combustion process of a three-gas mixture, mainly refinery gas, by incorporating more stable fuels such as natural gas and liquefied petroleum gas (LPG) to improve energy efficiency and reduce LPG consumption. Using real-time gas chromatography-mass spectrometry (GC-MS) data, we accurately calculate the mass fractions of individual compounds, allowing for more precise burner flow rate determinations. Thermochemical data are used to calculate equilibrium constants as a function of temperature, with the least squares method, while the Newton–Raphson method solves the resulting nonlinear equations. Four key variables (X4,X6,X8, and X11), representing H2,CO,O2, and N2, respectively, are defined, and a Jacobian matrix is constructed to ensure convergence within a tolerance of 1 ×106 over a maximum of 200 iterations, implemented via Python 3.10.4 and the scipy.optimize library. The optimization resulted in a reduction in LPG consumption by over 50%. By tailoring the fuel supply to the specific thermal needs of each processing unit, we achieved substantial energy savings. For instance, furnaces in the hydrocracking unit, which handle cleaner subproducts and benefit from hydrogen’s adiabatic reactions, require much less energy than those in the primary distillation unit, where high-impurity crude oil is processed. Full article
(This article belongs to the Section I: Energy Fundamentals and Conversion)
Show Figures

Figure 1

Figure 1
<p>Daily Gas Sample Data Recorded and Analyses.</p>
Full article ">Figure 2
<p>Equilibrium Constants <math display="inline"><semantics> <msub> <mi>K</mi> <mi>p</mi> </msub> </semantics></math> at 3600 °R.</p>
Full article ">Figure 3
<p>Levels fraction of harmful emissions.</p>
Full article ">Figure 4
<p>RG Density chromatography and model results.</p>
Full article ">Figure 5
<p>RG LHV chromatography and model results.</p>
Full article ">Figure 6
<p>Newton Rhapson applied to <math display="inline"><semantics> <msub> <mi>M</mi> <mrow> <mi>R</mi> <msub> <mi>G</mi> <mn>29</mn> </msub> </mrow> </msub> </semantics></math>.</p>
Full article ">Figure 7
<p>Newton Rhapson applied to <math display="inline"><semantics> <msub> <mi>M</mi> <mrow> <mi>R</mi> <msub> <mi>G</mi> <mn>87</mn> </msub> </mrow> </msub> </semantics></math>.</p>
Full article ">
21 pages, 20484 KiB  
Article
Structure and Strength Optimization of the Bogdan ERCV27 Electric Garbage Truck Spatial Frame Under Static Loading
by Kostyantyn Holenko, Oleksandr Dykha, Eugeniusz Koda, Ivan Kernytskyy, Orest Horbay, Yuriy Royko, Yevhen Fornalchyk, Oksana Berezovetska, Vasyl Rys, Ruslan Humenuyk, Serhii Berezovetskyi, Mariusz Żółtowski, Adam Baryłka, Anna Markiewicz, Tomasz Wierzbicki and Hydayatullah Bayat
Appl. Sci. 2024, 14(23), 11012; https://doi.org/10.3390/app142311012 - 27 Nov 2024
Viewed by 736
Abstract
Taking into account the requirements to reduce the release of harmful emissions into the environment, the EU’s environmental standards when transitioning to the Euro 7 standard in 2025 will actually lead vehicles having to operate without producing emissions in all driving situations. Carmakers [...] Read more.
Taking into account the requirements to reduce the release of harmful emissions into the environment, the EU’s environmental standards when transitioning to the Euro 7 standard in 2025 will actually lead vehicles having to operate without producing emissions in all driving situations. Carmakers believe that the new, much stricter regulations will mark the end of the internal combustion engine era. For example, in 2030, the manufacturer SEAT will cease its activities, leaving behind the Cupra brand, which will be exclusively electric in the future. This trend will apply not only to private vehicles (passenger cars), but also to utility vehicles, which is the subject of our research, namely the spatial tubular frame in the Bogdan ERCV27 garbage truck, presented in the form of a solid model. The peculiarity of the studied model is the installation of a battery block behind the driver’s cabin, causing an additional load to be placed on the spatial frame of the garbage truck, which in terms of its architecture is more like the body of a bus. During the conditions involving various modes of operation of a full-scale Bogdan ERCV27 garbage truck sample, questions about the strength and uniformity of its load-bearing spatial frame inevitably arise, which are decisive, even at the stage of designing and preparing the technical documentation. The main static load mode, which, despite its name, also covers dynamic conditions, was modeled using the appropriate coefficient kd = 2.0. The maximum stresses on the model during the “bending” mode were 381.13 MPa before structure optimization and 270.5 MPa as a result of the improvement measures. The spatial frame mass was reduced by 4.13%. During the “torsion” mode, the maximum deformation values were 12.1–14.5 mm, which guarantees the normal operation of the aggregates and units of the truck. Full article
Show Figures

Figure 1

Figure 1
<p>Bogdan ERCV27 garbage truck.</p>
Full article ">Figure 2
<p>The overall dimensions of the Bogdan ERCV27 garbage truck.</p>
Full article ">Figure 3
<p>Ansys module showing the Bogdan ERCV27 garbage truck.</p>
Full article ">Figure 4
<p>Meshed beam model of the Bogdan ERCV27 garbage truck: (<b>a</b>) rear part; (<b>b</b>) front part.</p>
Full article ">Figure 5
<p>Flat one-mass model, with suspension and wheels.</p>
Full article ">Figure 6
<p>Scheme of support for the “bending” mode.</p>
Full article ">Figure 7
<p>Scheme of support for “torsion” mode, with the hanging left rear wheel.</p>
Full article ">Figure 8
<p>Scheme of support for “bending” mode in Ansys Static Structural environment.</p>
Full article ">Figure 9
<p>Acceleration in the case of the “bending” mode.</p>
Full article ">Figure 10
<p>Mass distribution along the frame.</p>
Full article ">Figure 11
<p>Maximum stress before optimization.</p>
Full article ">Figure 12
<p>Maximum stress in the center frame part.</p>
Full article ">Figure 13
<p>Maximum stress in the cabin frame part.</p>
Full article ">Figure 14
<p>Safety factor map: (<b>a</b>) front part; (<b>b</b>) rear part.</p>
Full article ">Figure 15
<p>Garbage truck frame structural modifications (aka optimization steps): (<b>a</b>) roof dome; (<b>b</b>) channels; (<b>c</b>) front suspension mounting brackets; (<b>d</b>) sheets on the reverse side; (<b>e</b>) suspension brackets; and (<b>f</b>) rear axle suspension mounting bracket.</p>
Full article ">Figure 15 Cont.
<p>Garbage truck frame structural modifications (aka optimization steps): (<b>a</b>) roof dome; (<b>b</b>) channels; (<b>c</b>) front suspension mounting brackets; (<b>d</b>) sheets on the reverse side; (<b>e</b>) suspension brackets; and (<b>f</b>) rear axle suspension mounting bracket.</p>
Full article ">Figure 16
<p>Maximum stress after optimization.</p>
Full article ">Figure 17
<p>Safety factor after optimization.</p>
Full article ">Figure 18
<p>Frame deformation under “torsion” mode conditions: (<b>a</b>) front left wheel hanging; (<b>b</b>) rear right wheel hanging.</p>
Full article ">Figure 19
<p>Frame deformation under “torsion” mode conditions: (<b>a</b>) truck stands on the front right and rear left wheels; (<b>b</b>) front right and middle left wheels are supported.</p>
Full article ">
27 pages, 699 KiB  
Article
Estimating the Lifetime Parameters of the Odd-Generalized-Exponential–Inverse-Weibull Distribution Using Progressive First-Failure Censoring: A Methodology with an Application
by Mahmoud M. Ramadan, Rashad M. EL-Sagheer and Amel Abd-El-Monem
Axioms 2024, 13(12), 822; https://doi.org/10.3390/axioms13120822 - 25 Nov 2024
Viewed by 632
Abstract
This paper investigates statistical methods for estimating unknown lifetime parameters using a progressive first-failure censoring dataset. The failure mode’s lifetime distribution is modeled by the odd-generalized-exponential–inverse-Weibull distribution. Maximum-likelihood estimators for the model parameters, including the survival, hazard, and inverse hazard rate functions, are [...] Read more.
This paper investigates statistical methods for estimating unknown lifetime parameters using a progressive first-failure censoring dataset. The failure mode’s lifetime distribution is modeled by the odd-generalized-exponential–inverse-Weibull distribution. Maximum-likelihood estimators for the model parameters, including the survival, hazard, and inverse hazard rate functions, are obtained, though they lack closed-form expressions. The Newton–Raphson method is used to compute these estimations. Confidence intervals for the parameters are approximated via the normal distribution of the maximum-likelihood estimation. The Fisher information matrix is derived using the missing information principle, and the delta method is applied to approximate the confidence intervals for the survival, hazard rate, and inverse hazard rate functions. Bayes estimators were calculated with the squared error, linear exponential, and general entropy loss functions, utilizing independent gamma distributions for informative priors. Markov-chain Monte Carlo sampling provides the highest-posterior-density credible intervals and Bayesian point estimates for the parameters and reliability characteristics. This study evaluates these methods through Monte Carlo simulations, comparing Bayes and maximum-likelihood estimates based on mean squared errors for point estimates, average interval widths, and coverage probabilities for interval estimators. A real dataset is also analyzed to illustrate the proposed methods. Full article
Show Figures

Figure 1

Figure 1
<p>PDF for OGE-IWD.</p>
Full article ">Figure 2
<p>HRF for OGE-IWD.</p>
Full article ">Figure 3
<p>Description of the PFFC scheme.</p>
Full article ">Figure 4
<p>The KD, box, TTT, Q-Q, P-P, SF, PDF, and violin plots for the data set.</p>
Full article ">
17 pages, 4864 KiB  
Article
Lubrication Analysis of a Mechanical Seal Considering the Mixed Lubricant State of Gas and Liquid During External Power Shutdown of a Reactor Cooling Pump
by Youngjun Park, Gwanghee Hong, Sanghyun Jun, Jeongmook Choi, Taegyu Kim, Minsoo Kang and Gunhee Jang
Lubricants 2024, 12(12), 406; https://doi.org/10.3390/lubricants12120406 - 21 Nov 2024
Viewed by 660
Abstract
We proposed a method to calculate the pressure, opening force, and leakage rate in a mechanical seal under the mixed lubricant state of a gas and liquid for the mechanical seal in a reactor cooling pump (RCP) during external power loss. We calculated [...] Read more.
We proposed a method to calculate the pressure, opening force, and leakage rate in a mechanical seal under the mixed lubricant state of a gas and liquid for the mechanical seal in a reactor cooling pump (RCP) during external power loss. We calculated the pressure by solving the nonlinear finite element equation composed of the linear Reynolds equation of an incompressible liquid lubricant and the nonlinear Reynolds equation of a compressible gas lubricant using the Newton–Raphson method. In addition, we calculated the temperature distribution by solving the two-dimensional energy equation utilizing the finite element method. Additionally, we included the turbulence effect in the incompressible liquid lubricant and the turbulence and slip effects in the compressible gas lubricant. The accuracy of the developed program was validated by comparing the simulated opening force and leakage rate of both the mechanical seal with the liquid lubricant and the mechanical seal with the gas lubricant with prior research. Our analysis shows that in high-temperature environments, the increase in the gas region at the lubrication surface leads to a decrease in pressure and opening force and an increase in the leakage rate. Conversely, as the outer pressure increases, the gas region decreases, resulting in an increase in pressure, opening force, and leakage rate. Full article
Show Figures

Figure 1

Figure 1
<p>Mechanical seal configuration in a reactor coolant pump (RCP).</p>
Full article ">Figure 2
<p>Governing equation of the element matrix according to gas–liquid flow conditions.</p>
Full article ">Figure 3
<p>Numerical method for calculating pressure, opening force, and leakage of a mechanical seal under gas–liquid flow conditions.</p>
Full article ">Figure 4
<p>Comparison of leakage rate according to rotation speed [<a href="#B23-lubricants-12-00406" class="html-bibr">23</a>].</p>
Full article ">Figure 5
<p>Comparison of opening force (<b>a</b>) and leakage rate (<b>b</b>) with pressure changes [<a href="#B24-lubricants-12-00406" class="html-bibr">24</a>].</p>
Full article ">Figure 6
<p>Geometry of a stationary seal for RCP.</p>
Full article ">Figure 7
<p>Finite element model of a mechanical seal for RCP.</p>
Full article ">Figure 8
<p>Pressure distribution (<b>a</b>), temperature distribution (<b>b</b>), and fluid flow distribution (<b>c</b>) under normal operating conditions.</p>
Full article ">Figure 9
<p>Pressure distribution (<b>a</b>), temperature distribution (<b>b</b>), and fluid flow distribution (<b>c</b>) under external power shutdown conditions.</p>
Full article ">Figure 10
<p>Distribution of fluid flow according to outer pressure at (<b>a</b>) 313K and (<b>b</b>) 583K.</p>
Full article ">Figure 11
<p>Opening force according to outer pressure.</p>
Full article ">Figure 12
<p>Leakage rate according to outer pressure.</p>
Full article ">
Back to TopTop