Acoustic analysis of a three-dimensional cylindrical shell model under electromagnetic vibration

This paper presents the acoustic analysis of a three-dimensional cylindrical shell model under electromagnetic vibration, a critical factor affecting the performance of electric motors in various applications such as automotive, aerospace, and industrial systems. The study provides a multidisciplinary approach that integrates electromagnetics, structural vibration, and acoustics, solved using the fast multipole boundary element method (FMBEM). The results summarize the validation of the analytical models and numerical simulations, offering insights into effective vibration reduction methods. The conclusions indicate that the 3-D numerical analysis using FMBEM aligns well with the analytical solution for the sound pressure in the exterior acoustic domain of the cylindrical shell model. The paper contributes valuable insights for the design of low-noise motors and the control of electromagnetic vibration and noise in electric motors.

1 Introduction

The performance of electric motors, especially in applications such as automotive, aerospace, and industrial systems, is heavily influenced by their electromagnetic vibration and noise. These aspects not only affect the operational efficiency but also the comfort and reliability of motor-driven systems. With the rapid development of electric vehicles and advanced industrial automation, there is an increasing demand for motors that are efficient, compact, and silent. Therefore, the accurate prediction and control of their electromagnetic vibration and noise have become paramount. A significant body of research has been dedicated to understanding and mitigating the sources of electromagnetic vibration and noise in permanent magnet synchronous motors (PMSMs) and other types of electric motors.

Studies by Ballo et al. [1] and Xing et al. [2, 3] have focused on developing simplified analytical models and numerical prediction models to forecast the noise and vibration in PMSMs at the design stage. The influence of electromagnetic forces on motor vibration [4, 5] has also been a central theme. Strategies to mitigate vibrations and noise have been explored [6, 7]. Experimental studies by Torregrossa et al. [8] and Zhao et al. [9] have validated theoretical models and numerical simulations, providing insights into the effectiveness of various vibration reduction methods. The detection and analysis of oscillations in rail vehicle systems has been significantly advanced [10–12], particularly focusing on pantograph control, which not only improved the understanding of signal processing in this context but also provided practical solutions for real-time applications, enhancing the safety and efficiency of rail transportation. A multidisciplinary approach by Chai et al. [13] and Wu et al. [14], combining electromagnetics, structural mechanics, and acoustics, has been employed to provide a comprehensive understanding of motor behavior. The dynamic behavior of the rotor and the acoustic performance of the entire motor system have been examined [15, 16]. Optimization studies by Mendizabal et al. [17] have provided guidelines for designing low-noise motors. Certain studies have focused on specific aspects such as the effects of laminations [18], axial forces [19], and the application of amorphous alloys in stators [20], offering specialized insights into motor design.

The finite element method (FEM) has been extensively employed by Mao et al. [21] and Wang et al. [22] to predict acoustics, fracture mechanics, electromagnetics, and vibrations. Additionally, [23, 24] investigated the natural frequencies of motor components. However, there are several problems with FEM when modeling infinite domains. The boundary element method (BEM) has been used to tackle potential problems because it offers good accuracy and easy mesh construction. Particularly for exterior acoustic problems, the Sommerfeld radiation condition at infinity is rapidly satisfied [25]. The boundary integral problem has been quantitatively solved using the Galerkin approach for BEM implementation [26, 27].

In order to directly resolve the equation system, the conventional boundary element method (CBEM) produces a dense and non-symmetric coefficient matrix that takes a long time to compute. The fast multipole method (FMM) [28–30], the fast direct solver [31, 32], and the adaptive cross approximation approach [33] are only a few of the methods that have been employed to expedite the resolution of the integral issue. Architects and designers are increasingly considering changing the structural geometry to reduce noise. There is much potential for radiated noise reduction with this structural-acoustic optimization [34–36]. FEM and BEM may be employed with some computer-aided engineering (CAE) software. However, contemporary CAE requires that the models produced by CAD software be converted into simulation-ready models as part of the preprocessing phase. The transfer of geometric model data by the CAE results in geometry errors. One proposed approach to this problem [12, 37] is to combine BEM with geometric modeling and numerical simulation using isogeometric analysis (IGA) [38–40]. IGABEM has been employed to tackle an extensive variety of problems, including elastic mechanics [41], potential problems [42–46], heat transfer problems [47], wave propagation [48–53], fracture mechanics [54], electromagnetics [55–60], and structural optimization [61–66].

In this study, the acoustic analysis under electromagnetic vibration is solved using the fast multipole boundary element method (FMBEM). With regard to the advantages of FMBEM over CBEM, please refer to Chen et al. [67].

2 Analytical solution of sound pressure in a cylindrical shell model for external acoustic analysis

A cylindrical shell with radius $a$ and infinite length is the subject of the investigation. Consider a region of length $l$ in the model where the cylindrical shell vibrates and outside of which there are no displacements of the cylindrical shell. The acoustic vibration coupling is not taken into account while analyzing the sound field. The region outside the cylindrical shell when $r>a$ is the acoustical analysis domain. The model is presented in Figure 1.

Figure 1

Figure 1. Cylindrical shell model for the acoustical analysis.

Assuming a known radial displacement $w$ in Equation 1:

$w\left(z,\theta ,t\right)=W_m{\gamma }_m\left(z\right)\cos \left(n\theta \right)e^{-\mathrm{j}\omega t},(1)$

where ${\gamma }_m(z)$ is the modal function in the axial direction, $\cos (n\theta )$ is the modal function in the circumferential direction, $W_m$ is the amplitude, and $\mathrm{j}$ is the imaginary unit. After that, the time term $e^{-\mathrm{j}\omega t}$ is omitted from the calculation in the frequency domain.

Equation 1 is derived for time $t$ in order to obtain the velocity expression, as shown in Equation 2.

$\dot{w}\left(z,\theta \right)=-\mathrm{j}\omega W_m{\gamma }_m\left(z\right)\cos \left(n\theta \right).(2)$

The sound field in the domain must follow Equation 3 after being excited by the displacement in Equation 1.

$p\left(r,z,\theta \right)=P\left(r\right){\gamma }_m\left(z\right)\cos \left(n\theta \right).(3)$

A Fourier transform of Equations 1, 3 in the $z$ direction yields Equations 4, 5.

$\tilde{w}\left(k_z,\theta ,t\right)=W_m{\tilde{\gamma }}_m\left(k_z\right)\cos \left(n\theta \right),(4)$

$\tilde{p}\left(r,k_z,\theta \right)=P\left(r\right){\tilde{\gamma }}_m\left(k_z\right)\cos \left(n\theta \right),(5)$

in which $k_z=\omega /c_z$ is the wave number, and we have Equation 6.

${\tilde{\gamma }}_m\left(k_z\right)={\int }_{-\infty }^{\infty }{\gamma }_m\left(z\right)e^{-\mathrm{j}{k}_{z}z}\mathrm{d}z.(6)$

In the cylindrical coordinate system, the acoustic Helmholtz equation has the form shown in Equation 7.

$\frac{{\partial }^{2}p}{\partial r^2}+\frac{1}{r}\frac{\partial p}{\partial r}+\frac{1}{r^2}\frac{{\partial }^{2}p}{\partial {\theta }^2}+\frac{{\partial }^{2}p}{\partial z^2}+k_f^{2}p=0,(7)$

where the wave number $k_f=\omega /c_f$.

A Fourier transform of Equation 7 in the $z$ direction yields Equation 8.

$\frac{{\partial }^2\tilde{p}}{\partial r^2}+\frac{1}{r}\frac{\partial \tilde{p}}{\partial r}+\left[k_f^2-k_z^2-{\left(\frac{n}{r}\right)}^2\right]\tilde{p}=0.(8)$

Substituting Equation 5 into Equation 8 yields Equation 9.

$\left\{\frac{{\partial }^{2}P\left(r\right)}{\partial r^2}+\frac{1}{r}\frac{\partial P\left(r\right)}{\partial r}+\left[k_f^2-k_z^2-{\left(\frac{n}{r}\right)}^2\right]P\left(r\right)\right\}{\tilde{\gamma }}_m\left(k_z\right)\cos \left(n\theta \right)=0.(9)$

We need Equation 10 to make Equation 9 constant.

$\frac{{\partial }^{2}P\left(r\right)}{\partial r^2}+\frac{1}{r}\frac{\partial P\left(r\right)}{\partial r}+\left[k_f^2-k_z^2-{\left(\frac{n}{r}\right)}^2\right]P\left(r\right)=0.(10)$

The homogeneous equation shown in Equation 10 is a Bessel equation with the solution in Equation 11.

$\begin{array}{cc}\hfill P\left(r\right)=& A{J}_n\left(k_{r}r\right)+B{Y}_n\left(k_{r}r\right),\hfill \\ \hfill k_r=& {\left(k_f^2-k_z^2\right)}^{1/2}.\hfill \end{array}(11)$

When $r\to \infty$, the radiation field tends to be in the form of a plane wave, as shown in Equation 12.

$P\left(r\right)\to A{e}^{\mathrm{j}{k}_{r}r-\mathrm{j}\omega t}=A\left(\cos \left(k_{r}r\right)+\mathrm{j}\sin \left(k_{r}r\right)\right)e^{-\mathrm{j}\omega t}.(12)$

Furthermore, based on the approximation of two classes of Bessel functions at infinity in Equation 13,

$\begin{array}{c}\hfill \underset{k,r\to \infty }{\lim }{J}_n\left(k_{r}r\right)\to \sqrt{\frac{2}{\pi k_{r}r}}\cos \left(k_{r}r-\frac{2n+1}{4}\pi \right),\\ \hfill \underset{k,r\to \infty }{\lim }{Y}_n\left(k_{r}r\right)\to \sqrt{\frac{2}{\pi k_{r}r}}\sin \left(k_{r}r-\frac{2n+1}{4}\pi \right).\end{array}(13)$

$B=\mathrm{j}A$ is needed to satisfy the infinity condition. Then, we have Equation 14.

$P\left(r\right)=A\left[J_n\left(k_{r}r\right)+\mathrm{j}{Y}_n\left(k_{r}r\right)\right]=A{H}_n^{\left(1\right)}\left(k_{r}r\right),(14)$

where $H_n^{(1)}$ denotes the Hankel functions of the $n$-th order first kind.

The continuity conditions shown in Equation 15 must be satisfied at the interface:

$-{\left.\frac{1}{\mathrm{j}\omega {\rho }_f}\frac{\partial p}{\partial r}\right|}_{r=a}=\dot{w}{|}_{r=a}.(15)$

Performing a Fourier inverse transform of Equation 15, we have Equation 16.

$\frac{\partial \tilde{p}}{\partial r}=-{\omega }^2{\rho }_f\tilde{w}.(16)$

Substituting Equations 4, 5 into Equation 16 yields Equation 17.

$A=-\frac{{\omega }^2{\rho }_f{W}_m}{k_r{H}_n^{\left(1\right)\prime }\left(k_{r}a\right)},(17)$

and then, we have Equation 18.

$\tilde{p}=-\frac{{\omega }^2{\rho }_f{H}_n^{\left(1\right)}\left(k_{r}r\right)}{k_r{H}_n^{\left(1\right)\prime }\left(k_{r}a\right)}{\tilde{\gamma }}_m\left(k_z\right)W_m\cos \left(n\theta \right).(18)$

Performing a Fourier inverse transform of Equation 18, we finally obtain an analytical solution for the sound pressure at any point in the exterior acoustic domain for the cylindrical shell model, as shown in Equation 19.

$p\left(r,z,\theta \right)=-\frac{{\omega }^2{\rho }_f}{2\pi }\left[{\int }_{-\infty }^{\infty }\frac{H_n^{\left(1\right)}\left(k_{r}r\right)}{k_r{H}_n^{\left(1\right)\prime }\left(k_{r}a\right)}{\tilde{\gamma }}_m\left(k_z\right)e^{\mathrm{j}{k}_{z}z}\mathrm{d}{k}_z\right]W_m\cos \left(n\theta \right).(19)$

Note that $r>a$.

3 Numerical analysis of sound pressure in three-dimensional external acoustic analysis

Consider the Helmholtz governing equation in time-harmonic acoustic analysis, as shown in Equation 20.

${\mathrm{\nabla }}^{2}p\left(\mathbf{x}\right)+k^{2}p\left(\mathbf{x}\right)=0,\forall \mathbf{x}\in \mathrm{\Omega },(20)$

where ${\mathrm{\nabla }}^2$ is the Laplace operator, $p(\mathbf{x})$ is the sound pressure at field point $\mathbf{x}$, $k=\omega /c$ is the wave number, $\omega$ is the angular frequency of the incoming wave, $c$ is the sound speed in the domain $\mathrm{\Omega }$, and $\mathrm{\Omega }$ is the domain for acoustic analysis. In the frequency domain, the time-dependent component $e^{-\mathrm{j}\omega t}$ can be excluded from the computation.

Applying Green’s second theorem to the Helmholtz equation yields the following integral equation, as shown in Equation 21.

$p\left(\mathbf{x}\right)+{\int }_{\mathrm{\Gamma }}\frac{\partial G\left(\mathbf{x},\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{y}\right)}p\left(\mathbf{y}\right)\mathrm{d}\mathrm{\Gamma }\left(\mathbf{y}\right)={\int }_{\mathrm{\Gamma }}G\left(\mathbf{x},\mathbf{y}\right)q\left(\mathbf{y}\right)\mathrm{d}\mathrm{\Gamma }\left(\mathbf{y}\right)+p^{\mathrm{i}\mathrm{n}\mathrm{c}}\left(\mathbf{x}\right),\mathbf{y}\in \mathrm{\Gamma },(21)$

where $\mathrm{\Gamma }=\partial \mathrm{\Omega }$ is the border of the domain $\mathrm{\Omega }$, $\mathbf{y}$ is the sound source point on the boundary $\mathrm{\Gamma }$, $\mathbf{x}$ is the field point affected by the sound source $\mathbf{y}$, $G(\mathbf{x},\mathbf{y})$ is the Green’s function, $\mathbf{n}(\mathbf{y})$ is the outward unit normal vector at point $\mathbf{y}$, $\partial ()/\partial \mathbf{n}=\nabla ()\cdot \mathbf{n}$ is the exterior derivative, $p^{\text{inc}}$ is the sound pressure of the incident wave, and $q(\mathbf{y})$ is the outward flux at point $\mathbf{y}$. $q(\mathbf{y})$ satisfies the Neumann boundary condition, as shown in Equation 22.

$q\left(\mathbf{y}\right)=\frac{\partial p\left(\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{y}\right)}=\mathrm{j}\rho \omega v_f\left(\mathbf{y}\right),(22)$

where $\rho$ is the density of the medium in which the sound wave propagates, and $v_f(\mathbf{y})$ is the normal velocity at point $\mathbf{y}$.

To obtain the unknown sound pressure on the boundary $\mathrm{\Gamma }$, the field point $\mathbf{x}$ may converge to the boundary, and thus, the Kirchhoff–Helmholtz conventional boundary integral equation (CBIE) is obtained from Equation 21, as shown in Equation 23.

$c\left(\mathbf{x}\right)p\left(\mathbf{x}\right)+{\int }_{\mathrm{\Gamma }}\frac{\partial G\left(\mathbf{x},\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{y}\right)}p\left(\mathbf{y}\right)\mathrm{d}\mathrm{\Gamma }\left(\mathbf{y}\right)={\int }_{\mathrm{\Gamma }}G\left(\mathbf{x},\mathbf{y}\right)q\left(\mathbf{y}\right)\mathrm{d}\mathrm{\Gamma }\left(\mathbf{y}\right)+p^{\mathrm{i}\mathrm{n}\mathrm{c}}\left(\mathbf{x}\right),\mathbf{x},\mathbf{y}\in \mathrm{\Gamma }.(23)$

The normal derivative boundary integral equation (NDBIE) of Equation 23 is given by Equation 24.

$c\left(\mathbf{x}\right)q\left(\mathbf{x}\right)+{\int }_{\mathrm{\Gamma }}\frac{{\partial }^{2}G\left(\mathbf{x},\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{y}\right)\partial \mathbf{n}\left(\mathbf{x}\right)}p\left(\mathbf{y}\right)\mathrm{d}\mathrm{\Gamma }\left(\mathbf{y}\right)={\int }_{\mathrm{\Gamma }}\frac{\partial G\left(\mathbf{x},\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{x}\right)}q\left(\mathbf{y}\right)\mathrm{d}\mathrm{\Gamma }\left(\mathbf{y}\right)+\frac{\partial p^{\text{inc}}\left(\mathbf{x}\right)}{\partial \mathbf{n}\left(\mathbf{x}\right)},\mathbf{x},\mathbf{y}\in \mathrm{\Gamma }.(24)$

In Equations 23, 24, $c(\mathbf{x})$ = 1/2 when the boundary around point $\mathbf{x}$ is smooth.

Applying only CBIE or NDBIE leads to non-uniqueness of the solution to the exterior sound field analysis, which can be solved by linearly combining CBIE and NDBIE. This is called the Burton–Miller method [68], as shown in Equation 25.

$\begin{array}{l}\hfill c\left(\mathbf{x}\right)p\left(\mathbf{x}\right)+\alpha c\left(\mathbf{x}\right)q\left(\mathbf{x}\right)+{\int }_{\mathrm{\Gamma }}\frac{\partial G\left(\mathbf{x},\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{y}\right)}p\left(\mathbf{y}\right)\mathrm{d}\mathrm{\Gamma }\left(\mathbf{y}\right)+\alpha {\int }_{\mathrm{\Gamma }}\frac{{\partial }^{2}G\left(\mathbf{x},\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{y}\right)\partial \mathbf{n}\left(\mathbf{x}\right)}p\left(\mathbf{y}\right)\mathrm{d}\mathrm{\Gamma }\left(\mathbf{y}\right)\\ \hfill ={\int }_{\mathrm{\Gamma }}G\left(\mathbf{x},\mathbf{y}\right)q\left(\mathbf{y}\right)\mathrm{d}\mathrm{\Gamma }\left(\mathbf{y}\right)+\alpha {\int }_{\mathrm{\Gamma }}\frac{\partial G\left(\mathbf{x},\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{x}\right)}q\left(\mathbf{y}\right)\mathrm{d}\mathrm{\Gamma }\left(\mathbf{y}\right)+p^{\mathrm{i}\mathrm{n}\mathrm{c}}\left(\mathbf{x}\right)+\alpha \frac{\partial p^{\text{inc}}\left(\mathbf{x}\right)}{\partial \mathbf{n}\left(\mathbf{x}\right)},\mathbf{x},\mathbf{y}\in \mathrm{\Gamma },\end{array}(25)$

in which $\alpha$ is the coupling coefficient, and we have Equation 26.

$\alpha =\left\{\begin{array}{cc}\mathrm{j},\hfill & k<1,\hfill \\ -\frac{\mathrm{j}}{k},\hfill & k\ge 1.\hfill \end{array}\right.(26)$

In Equation 25, for 2-D acoustic analysis, we have Equation 27.

$\begin{array}{cc}\hfill G\left(\mathbf{x},\mathbf{y}\right)=& \frac{\mathrm{j}}{4}{H}_0^{\left(1\right)}\left(kr\right),\hfill \\ \hfill \frac{\partial G\left(\mathbf{x},\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{x}\right)}=& -\frac{\mathrm{j}k}{4}{H}_1^{\left(1\right)}\left(kr\right)\frac{\partial r}{\partial \mathbf{n}\left(\mathbf{x}\right)},\hfill \\ \hfill \frac{{\partial }^{2}G\left(\mathbf{x},\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{y}\right)\partial \mathbf{n}\left(\mathbf{x}\right)}=& \frac{\mathrm{j}{k}^2}{4}{H}_2^{\left(1\right)}\left(kr\right)\frac{\partial r}{\partial \mathbf{n}\left(\mathbf{y}\right)}\frac{\partial r}{\partial \mathbf{n}\left(\mathbf{x}\right)}+\frac{\mathrm{j}k}{4r}{H}_1^{\left(1\right)}\left(kr\right)n_i\left(\mathbf{y}\right)n_i\left(\mathbf{x}\right),\hfill \end{array}(27)$

and for the 3-D acoustic problem, we have Equation 28.

$\begin{array}{cc}G\left(\mathbf{x},\mathbf{y}\right)=\hfill & \frac{{\mathrm{e}}^{\mathrm{j}kr}}{4\pi r},\hfill \\ \hfill \frac{\partial G\left(\mathbf{x},\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{x}\right)}=& -\frac{{\mathrm{e}}^{\mathrm{j}kr}}{4\pi r^2}\left(1-\mathrm{j}kr\right)\frac{\partial r}{\partial \mathbf{n}\left(\mathbf{x}\right)},\hfill \\ \hfill \frac{{\partial }^{2}G\left(\mathbf{x},\mathbf{y}\right)}{\partial \mathbf{n}\left(\mathbf{y}\right)\partial \mathbf{n}\left(\mathbf{x}\right)}=& \frac{{\mathrm{e}}^{\mathrm{j}kr}}{4\pi r^3}\left[\left(3-3\mathrm{j}kr-k^2{r}^2\right)\frac{\partial r}{\partial \mathbf{n}\left(\mathbf{y}\right)}\frac{\partial r}{\partial \mathbf{n}\left(\mathbf{x}\right)}+\left(1-\mathrm{j}kr\right)n_i\left(\mathbf{y}\right)n_i\left(\mathbf{x}\right)\right],\hfill \end{array}(28)$

where $H_n^{(1)}$ is the Hankel function of $n$-th order first kind, $r$ is the distance between point $\mathbf{x}$ and $\mathbf{y}$, $r=|\mathbf{x}-\mathbf{y}|$, and $n_i$ is the component of the coordinate.

In Equation 25, discretizing the boundary $\mathrm{\Gamma }$ yields the linear equations shown in Equation 29.

$\mathbf{H}\mathbf{p}=\mathbf{G}\mathbf{q}+{\mathbf{p}}_{\mathbf{i}},(29)$

where ${\mathbf{p}}_{\mathbf{i}}$ is the nodal pressure from $p^{\mathit{inc}}$.

The sound pressure at point $\mathbf{x}$ could be obtained by solving Equation 29.

4 Numerical example

In this part, we will calculate the sound pressure at some certain location using the analytical equation (Equation 19) and numerical analysis (using BEM), respectively. Both are implemented using our in-house Fortran code. The algorithm is crafted using the Fortran 90 programming language and compiled with the combination of Visual Studio 2022 and Intel^®oneAPI 2022 toolkit. It is executed on a PC with an Intel(R) Xeon(R) Bronze 3204 CPU $@$ 1.90 GHz and complemented by 128.0 GB of RAM. In this part, the coordinates of the points are described using the cylindrical coordinate system in Equation 19.

Consider the cylindrical shell model shown in Figure 1. The cylindrical shell of radius $a$ = 0.5 m and length $l$ = 0.4 m is closed at both ends for the three-dimensional (3-D) numerical analysis. Consider Equation 1. Let the vibration of the cylindrical shell be uniformly distributed in the axial direction; that is, let ${\gamma }_m$ = 1. Let the amplitude of the vibration $W_m$ = $2.24\times 10^{-8}$ m, the frequency $f=\omega /2\pi$ = [100, 1,000] Hz, and the order of the circumferential vibration order be $n$ = 0.

The cylindrical shell model is meshed with quadrilateral elements. Figures 2A, B gives the meshing scheme with element number of 2112 and 8320, separately. Figure 3 provides the numerical solutions of the sound pressure at ($\sqrt{2}/2$, 0.0, $\pi /4$) — (0.5, 0.5, 0.0) in the cartesian coordinates for the meshing schemes in Figure 2, where the center of the cylindrical shell is at (0.0, 0.0, 0.0). The comparison in Figure 3 is crucial, as it tests the sensitivity of the numerical solution to the mesh density. A fine mesh provides a more accurate representation of the continuous problem but at a higher computational cost. Conversely, a coarse mesh reduces computational demand but may lead to inaccuracies. The two mesh schemes have 2,112 and 8,320 elements, respectively. The agreement between these two schemes suggests that the coarser mesh (2,112 elements) is sufficiently refined for the problem at hand, capturing the essential features of the sound pressure distribution without the excessive computational cost associated with a very fine mesh. Good agreement between the two mesh schemes indicates that the numerical solution has converged to a value that accurately represents the sound pressure at the specified point. This convergence is a critical aspect of numerical methods, ensuring that the solution is reliable and independent of the discretization. The consistency of results from different mesh densities also serves as an indirect validation of the FMBEM approach. It demonstrates that the method is capable of providing accurate results even with relatively coarse meshes, which is particularly beneficial for complex geometries and larger models where fine meshing becomes impractical. Based on this, we will select the element size configuration shown in Figure 2A for the cylindrical shell model with radius $a$ = 0.5 m and length $l\ge$ 0.4 m in the subsequent analysis.

Figure 2

Figure 2. Meshing scheme for the cylindrical shell model with closed ends. Radius $a$ = 0.5 m and length $l$ = 0.4 m. (A) Model meshing. Number of elements = 2,112. (B) Model meshing. Number of elements = 8,320.

Figure 3

Figure 3. Numerical results of sound pressure at ($\sqrt{2}/2$, 0.0, $\pi /4$). Circumferential vibration order $n$ = 0, length $l$ = 0.4 m; number of elements = 2,112 vs. 8,320.

Consider the model presented in Figure 1. The cylindrical shell model for analytical analysis is of infinite length, in which only a segment of length $l$ vibrates. This infinite-length cylindrical shell is considered as the boundary of the exterior acoustic domain, meaning that the acoustic analysis does not take into account the interior region of the cylindrical shell. However, during the 3-D numerical analysis, the boundary of the exterior acoustic domain is not a cylindrical shell of infinite length but a shell of length $l$ with closed ends, which causes the space outside the $l$-length segment in the infinitely long cylindrical shell to be involved in the acoustic analysis. Thus, computing errors are introduced. For the cylindrical shell model with a radius of 0.5 m and a length of 0.4 m, the analytical and 3-D numerical solutions for the sound pressure at ($\sqrt{2}/2$, 0.0, $\pi /4$) are provided in Figure 4. The two curves do not fit well, as we have predicted. This discrepancy is expected and can be attributed to several factors. The analytical model assumes an infinite length for the cylindrical shell, which simplifies the boundary conditions but deviates from the actual finite length used in the numerical model. The numerical model includes the effects of the shell’s finite length, which introduces additional complexities not fully captured by the analytical model’s assumptions. The discrepancy highlights the limitations of analytical models when they are based on simplified assumptions that may not fully represent the physical system’s complexity. The comparison underscores the importance of validating analytical models with numerical simulations, especially when the models are used for design and optimization in engineering applications. The results suggest that the model configuration, particularly the length of the cylindrical shell, significantly affects the accuracy of the sound pressure predictions. This insight is crucial for selecting appropriate model parameters in both theoretical and numerical studies. The choice between analytical and numerical methods should be guided by the specific requirements of the problem, including the need for accuracy, computational resources, and the complexity of the model.

Figure 4

Figure 4. Numerical vs. analytical results of sound pressure at ($\sqrt{2}/2$, 0.0, $\pi /4$) for a cylindrical shell model of length $l$ = 0.4 m. Circumferential vibration order $n$ = 0.

To reduce the computational errors resulting from the axial external space of the cylindrical shell, we attempt to lengthen the shell model with closed ends in the 3-D numerical analysis. For models of various lengths, the numerical results are compared with the analytical solution to determine an appropriate model configuration. Figure 5 gives the sound pressure at ($\sqrt{2}/2$, 0.0, $\pi /4$) for cylindrical shell models of different lengths. The circumferential vibration order $n$ = 0. Figures 5A–F show that the 3-D numerical result agrees well with the analytical solution when the length l $\ge$ 4 m. The analytical model assumes an infinite length for the cylindrical shell, which simplifies the boundary conditions but deviates from the actual finite length used in the numerical model. The numerical model includes the effects of the shell’s finite length, which introduces additional complexities not fully captured by the analytical model’s assumptions. The discrepancy between the analytical and numerical solutions for shorter shells highlights the limitations of analytical models when they are based on simplified assumptions that may not fully represent the physical system’s complexity. The comparison underscores the importance of validating analytical models with numerical simulations, especially when the models are used for design and optimization in engineering applications.

Figure 5

Figure 5. Numerical vs. analytical results of sound pressure at ($\sqrt{2}/2$, 0.0, $\pi /4$) for cylindrical shell models of different lengths. Circumferential vibration order $n$ = 0. (A) Length $l$ = 1 m, number of elements = 2,880. (B) Length $l$ = 2 m, number of elements = 4,160. (C) Length $l$ = 4 m, number of elements = 6,784. (D) Length $l$ = 5 m, number of elements = 8,064. (E) Length $l$ = 9 m, number of elements = 13,312. (F) Length $l$ = 10 m, number of elements = 14,592.

The data in Figures 5C–F, 6 give the relative error of the 3-D numerical solution to the analytical result. Similar characteristics are displayed by the four curves, namely, a rising relative error between the analytical and numerical solutions with increasing frequency, indicating that higher frequencies are more challenging to model accurately. This could be due to the higher spatial and temporal resolution required for capturing the dynamics of higher-frequency waves. For the cylindrical shell model of length $l$ = 4 m, the relative error range is (−5%, 1%); for the model of $l$ = 5 m, it is (−4%, 1%); and for the models of $l$ = 9 m and $l$ = 10 m, it is about (−3%, −1%). The relative error decreases, and the computational accuracy generally improves as the model lengthens. The present shell model has a radius of 0.5 m. The model length is 10 times the radius when $l$ = 5 m and 20 times the radius when $l$ = 10 m. The results suggest that the length of the model plays a critical role in the accuracy of numerical simulations. The analytical model assumes an infinite length for simplicity, but the numerical model must account for the finite length, which introduces additional complexities. The findings indicate that a model length of at least 10 times the radius is necessary for accurate uncoupling analysis. Although increasing the model length improves accuracy, it also increases the computational cost in terms of memory usage and processing time. Therefore, a balance must be struck between accuracy and computational efficiency, which is achieved when the model length is approximately 10 times the radius. In our future work, acceleration algorithms may be employed to speed up the calculation.

Figure 6

Figure 6. Relative error, numerical results vs. analytical results. Circumferential vibration order $n$ = 0, length $l$ = 4 m, 5 m, 9 m, and 10 m.

Consider the cylindrical shell model shown in Figure 1. The cylindrical shell is closed at both ends for the 3-D numerical analysis. The radius $a$ = 0.5 m. The length $l$ = 5 m or $l$ = 10 m for the 3-D numerical analysis. Consider Equation 1. Let ${\gamma }_m$ = 1, the amplitude of the vibration $W_m$ = $2.24\times 10^{-8}$ m, and the frequency $f=\omega /2\pi$ = [100, 1,000] Hz. The analytical and 3-D numerical solutions for the sound pressure at ($\sqrt{2}/2$, 0.0, $\pi /4$) for different circumferential vibration orders $n$ are given in Figure 7. Figures 7A–F show that the 3-D numerical result agrees well with the analytical solution when the circumferential vibration order $n$ = 1, 3, or 4. The data in Figures 7A, 8 provide, as an example, the relative error between the 3-D numerical result and the analytical solution.

Figure 7

Figure 7. Sound pressure at ($\sqrt{2}/2$, 0.0, $\pi /4$) and the circumferential vibration order $n$ varies. Numerical results vs. analytical result, in which L5 is length $l$ = 5 m, and L10 is length $l$ = 10 m. (A) $n$ = 1. (B) $n$ = 2. (C) $n$ = 3. (D) $n$ = 4. (E) $n$ = 6. (F) $n$ = 10.

Figure 8

Figure 8. Numerical results. Sound pressure on the cylindrical shell with a length of 5 m. Frequency $f$ = 500 Hz. Circumferential vibration order $n$ varies. The black line represents the section for $\theta$ = $pi/4$. (A) $n$ = 0. (B) $n$ = 1. (C) $n$ = 2. (D) $n$ = 3. (E) $n$ = 4. (F) $n$ = 6. (G) $n$ = 8. (H) $n$ = 10.

Figure 7 shows that for n = 2, 6, and 10, the analytical solution goes to zero, and the numerical solution is less than 3.0E−12, which might be thought of as tending toward zero in numerical analysis. Note that the sound pressure at ($\sqrt{2}/2$, 0.0, $\pi /4$) is calculated in Figure 7. Consider Equation 19. When $\theta$ = $\pi /4$ and $n$ = 2, 6, or 10, $\cos (n\theta )$ = 0, meaning that the sound pressure $p(r,z,\theta )$ = 0 in Equation 19, which coincides with Figures 7B, E, F.

The agreement between analytical and numerical results for certain modes validates the theoretical model’s assumptions and the mathematical formulations that describe the sound pressure distribution. This is crucial for the reliability of predictions made using these models in practical applications. The results highlight the importance of considering different vibration modes in the analysis. The theoretical model predicts different behaviors for different modes, and the numerical simulations confirm these predictions, emphasizing the need for accurate modeling of vibration modes in the design and analysis of electric motors. The findings provide valuable insights for the design and optimization of electric motors, particularly in controlling the vibration modes that contribute to noise. Understanding which modes contribute to the sound pressure can guide the design of motor components to minimize noise emissions. The choice between analytical and numerical methods should be guided by the specific requirements of the problem, including the need for accuracy, computational resources, and the complexity of the model. The agreement between the two methods for certain modes suggests that, under these conditions, either method could be reliably used.

For the closed cylindrical shell model with a length of 5 m, Figures 8A–H provide the numerical sound pressure on the structure surface for various circumferential vibration orders, in which the frequency $f$ = 500 Hz. The black line in the figure represents the section for $\theta$ = $\pi /4$. In Equation 19, when $\theta$ = 0, $\cos (n\theta )$ = 1, and the sound pressure $p(r,z,\theta )$ reaches its maximum. This is in accordance with Figure 8I. Furthermore, in Equation 19, when $n$ = 0, $\cos (n\theta )$ = 1, indicating that the sound pressure $p(r,z,\theta )$ is constant along the cylindrical shell’s circumferential direction. This finding aligns with Figure 8A. Still in Equation 19, when $\theta =\pi /4$ and $n$ = 2, 6, or 10, $\cos (n\theta )$ = 0, indicating that the sound pressure $p(r,z,\theta )$ = 0. This conclusion is consistent with Figures 8C, F, H. Remaining in Equation 19, when $\theta =\pi /4$ and $n$ = 4 or 8, abs$(\cos (n\theta ))$ = 1, indicating that the sound pressure $p(r,z,\theta )$ is maximized. This is in accordance with Figures 8E, G. These results demonstrate good agreement between the results of our 3-D numerical analysis and the analytical solution. The agreement between the numerical results and the analytical predictions validates the analytical model’s accuracy in predicting the sound pressure distribution on the cylindrical shell under different vibration modes. This is essential for the model’s reliability in practical applications. The results highlight the importance of considering different vibration modes in the analysis. Each mode affects the sound pressure distribution, which is crucial for noise control and motor design.

5 Conclusion

This study concludes that the 3-D numerical analysis using FMBEM aligns well with the analytical solution for the sound pressure in the exterior acoustic domain of a cylindrical shell model, particularly when the model length is sufficient relative to its radius. The investigation reveals that computational accuracy improves as the model lengthens, suggesting a model length of at least 10 times the radius for accurate uncoupling analysis. The findings provide valuable insights for the design of low-noise motors and contribute to the understanding of electromagnetic vibration and noise control in electric motors.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material; further inquiries can be directed to the corresponding author.

Author contributions

YX: conceptualization, investigation, methodology, project administration, and writing–original draft. JW: formal analysis, resources, supervision, and writing–review and editing. SY: data curation, validation, and writing–original draft. GL: software, visualization, and writing–review and editing. KZ: funding acquisition, validation, and writing–review and editing

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The author(s) were sponsored by the Henan Provincial Key R&D and Promotion Project under Grant No. 232102220033, the National Natural Science Foundation of China under Grant No. 42207200, the Natural Science Foundation of Henan under Grant No. 222300420498, the Zhumadian 2023 Major Science and Technology Special Project under Grant No. ZMDSZDZX2023002, and the Postgraduate Education Reform and Quality Improvement Project of Henan Province under Grant No. YJS2023JD52.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Ballo F, Gobbi M, Mastinu G, Palazzetti R. Noise and vibration of permanent magnet synchronous electric motors: a simplified analytical model. IEEE Trans Transportation Electrification (2023) 9:2486–96. doi:10.1109/TTE.2022.3209917

CrossRef Full Text | Google Scholar

2. Xing Z, Wang X, Zhao W, Sun L, Niu N. Calculation method for natural frequencies of stator of permanent magnet synchronous motors based on three-dimensional elastic theory. IEEE Trans Energy Convers (2021) 36:755–66. doi:10.1109/TEC.2020.3030042

CrossRef Full Text | Google Scholar

3. Xing Z, Wang X, Zhao W. Fast calculation of electromagnetic vibration of surface-mounted pmsm considering teeth saturation and tangential electromagnetic force. IEEE Trans Ind Electronics (2024) 71:316–26. doi:10.1109/TIE.2023.3243298

CrossRef Full Text | Google Scholar

4. Zhu S, Zhao W, Ji J, Liu G, Mao Y, Liu T. Investigation of bread-loaf magnet on vibration performance in fscw pmsm considering force modulation effect. IEEE Trans Transportation Electrification (2021) 7:1379–89. doi:10.1109/TTE.2020.3035180

CrossRef Full Text | Google Scholar

5. Hong J, Wang S, Sun Y, Sun X, Cao H. Piecewise stagger poles with continuous skew edge for vibration reduction in surface-mounted pm synchronous machines. IEEE Trans Ind Electronics (2021) 68:8498–506. doi:10.1109/TIE.2020.3013535

CrossRef Full Text | Google Scholar

6. Lin F, Zuo SG, Deng WZ, Wu SL. Reduction of vibration and acoustic noise in permanent magnet synchronous motor by optimizing magnetic forces. J Sound Vibration (2018) 429:193–205. doi:10.1016/j.jsv.2018.05.018

CrossRef Full Text | Google Scholar

7. Hu S, Zuo S, Wu H, Liu M. An analytical method for calculating the natural frequencies of a motor considering orthotropic material parameters. IEEE Trans Ind Electronics (2019) 66:7520–8. doi:10.1109/TIE.2018.2883194

CrossRef Full Text | Google Scholar

8. Torregrossa D, Peyraut F, Fahimi B, M’Boua J, Miraoui A. Multiphysics finite-element modeling for vibration and acoustic analysis of permanent magnet synchronous machine. IEEE Trans Energy Convers (2011) 26:490–500. doi:10.1109/TEC.2010.2080681

CrossRef Full Text | Google Scholar

9. Zhao J, Wang Z, Liu H, Ning F, Hong X, Du J, et al. Modal analysis and structure optimization of permanent magnet synchronous motor. IEEE Access (2020) 8:151856–65. doi:10.1109/ACCESS.2020.3017679

CrossRef Full Text | Google Scholar

10. Mercorelli P. Denoising and harmonic detection using nonorthogonal wavelet packets in industrial applications. J Syst Sci Complexity (2007) 20:325–43. doi:10.1007/s11424-007-9028-z

CrossRef Full Text | Google Scholar

11. Mercorelli P. Detection of oscillations with application in the pantograph control. Recent Patents Electr Eng (2011) 4:71–7. doi:10.2174/1874476111104010071

CrossRef Full Text | Google Scholar

12. Mercorelli P. A denoising procedure using wavelet packets for instantaneous detection of pantograph oscillations. Mech Syst Signal Process (2013) 35:137–49. doi:10.1016/j.ymssp.2012.09.001

CrossRef Full Text | Google Scholar

13. Chai F, Li Y, Pei Y, Li Z. Accurate modelling and modal analysis of stator system in permanent magnet synchronous motor with concentrated winding for vibration prediction. IET Electric Power Appl (2018) 12:1225–32. doi:10.1049/iet-epa.2017.0813

CrossRef Full Text | Google Scholar

14. Wu Z, Fan Y, Chen H, Wang X, Lee CHT. Electromagnetic force and vibration study of dual-stator consequent-pole hybrid excitation motor for electric vehicles. IEEE Trans Vehicular Technology (2021) 70:4377–88. doi:10.1109/TVT.2021.3075461

CrossRef Full Text | Google Scholar

15. Park S, Kim W, Kim SI. A numerical prediction model for vibration and noise of axial flux motors. IEEE Trans Ind Electronics (2014) 61:5757–62. doi:10.1109/TIE.2014.2300034

CrossRef Full Text | Google Scholar

16. Čorović S, Miljavec D. Modal analysis and rotor-dynamics of an interior permanent magnet synchronous motor: an experimental and theoretical study. Appl Sci (2020) 10:5881. doi:10.3390/app10175881

CrossRef Full Text | Google Scholar

17. Mendizabal M, McCloskey A, Poza J, Zarate S, Iriondo J, Irazu L. Optimum slot and pole design for vibration reduction in permanent magnet synchronous motors. Appl Sci (2021) 11:4849. doi:10.3390/app11114849

CrossRef Full Text | Google Scholar

18. Wang H, Williams K. Effects of laminations on the vibrational behaviour of electrical machine stators. J Sound Vibration (1997) 202:703–15. doi:10.1006/jsvi.1996.0845

CrossRef Full Text | Google Scholar

19. Deng W, Zuo S. Axial force and vibroacoustic analysis of external-rotor axial-flux motors. IEEE Trans Ind Electronics (2018) 65:2018–30. doi:10.1109/TIE.2017.2739697

CrossRef Full Text | Google Scholar

20. Xu X, Han Q, Qin Z, Chu F. Analytical methods for the radial electromagnetic vibration of stator in permanent magnet motors with an amorphous alloy core. Mech Syst Signal Process (2020) 145:106909. doi:10.1016/j.ymssp.2020.106909

CrossRef Full Text | Google Scholar

21. Mao Y, Zhao W, Zhu S, Chen Q, Ji J. Vibration investigation of spoke-type pm machine with asymmetric rotor considering modulation effect of stator teeth. IEEE Trans Ind Electronics (2021) 68:9092–103. doi:10.1109/TIE.2020.3022530

CrossRef Full Text | Google Scholar

22. Wang S, Hong J, Sun Y, Cao H. Mechanical and magnetic pivot roles of tooth in vibration of electrical machines. IEEE Trans Energy Convers (2021) 36:139–48. doi:10.1109/TEC.2020.3007802

CrossRef Full Text | Google Scholar

23. Xing Z, Wang X, Zhao W, Wang F. Calculation of stator natural frequencies of electrical machines considering complex boundary conditions. IEEE Trans Industry Appl (2022) 58:7079–87. doi:10.1109/TIA.2022.3191301

CrossRef Full Text | Google Scholar

24. Zhao G, Sun L, Niu N, Wang X, Xing Z. Analysis of vibration characteristics of stators of electrical machines under actual boundary. Mech Syst Signal Process (2023) 185:109778. doi:10.1016/j.ymssp.2022.109778

CrossRef Full Text | Google Scholar

25. Sommerfeld A. Partial differential equations in physics. Academic Press (1949). v-vi. doi:10.1016/B978-0-12-654658-3.50003-3

CrossRef Full Text | Google Scholar

26. Engleder OSS. Stabilized boundary element methods for exterior Helmholtz problems. Numerische Mathematik (2008) 110:145–60. doi:10.1007/s00211-008-0161-y

CrossRef Full Text | Google Scholar

27. Chen L, Zhang Y, Lian H, Atroshchenko E, Ding C, Bordas S. Seamless integration of computer-aided geometric modeling and acoustic simulation: isogeometric boundary element methods based on catmull-clark subdivision surfaces. Adv Eng Softw (2020) 149:102879. doi:10.1016/j.advengsoft.2020.102879

CrossRef Full Text | Google Scholar

28. Greengard L, Rokhlin V. A fast algorithm for particle simulations. J Comput Phys (1987) 73:325–48. doi:10.1016/0021-9991(87)90140-9

CrossRef Full Text | Google Scholar

29. Coifman R, Rokhlin V, Wandzura S. The fast multipole method for the wave equation: a pedestrian prescription. IEEE Antennas Propagation Mag (1993) 35:7–12. doi:10.1109/74.250128

CrossRef Full Text | Google Scholar

30. Rokhlin V. Diagonal forms of translation operators for the Helmholtz equation in three dimensions. Appl Comput Harmonic Anal (1993) 1:82–93. doi:10.1006/acha.1993.1006

CrossRef Full Text | Google Scholar

31. Martinsson P, Rokhlin V. A fast direct solver for boundary integral equations in two dimensions. J Comput Phys (2005) 205:1–23. doi:10.1016/j.jcp.2004.10.033

CrossRef Full Text | Google Scholar

32. Martinsson P, Rokhlin V. A fast direct solver for scattering problems involving elongated structures. J Comput Phys (2007) 221:288–302. doi:10.1016/j.jcp.2006.06.037

CrossRef Full Text | Google Scholar

33. Bebendorf SMR. Adaptive low-rank approximation of collocation matrices. Computing (2003) 70:1–24. doi:10.1007/s00607-002-1469-6

CrossRef Full Text | Google Scholar

34. Kim NH, Dong J. Shape sensitivity analysis of sequential structural-acoustic problems using FEM and BEM. J Sound Vibration (2006) 290:192–208. doi:10.1016/j.jsv.2005.03.013

CrossRef Full Text | Google Scholar

35. Chen L, Zhao J, Lian H, Yu B, Atroshchenko E, Li P. A BEM broadband topology optimization strategy based on Taylor expansion and SOAR method-Application to 2D acoustic scattering problems. Int J Numer Methods Eng (2023) 124:5151–82. doi:10.1002/nme.7345

CrossRef Full Text | Google Scholar

36. Qu Y, Zhou Z, Chen L, Lian H, Li X, Hu Z, et al. Uncertainty quantification of vibro-acoustic coupling problems for robotic manta ray models based on deep learning. Ocean Eng (2024) 299:117388. doi:10.1016/j.oceaneng.2024.117388

CrossRef Full Text | Google Scholar

37. Simpson R, Bordas S, Trevelyan J, Rabczuk T. A two-dimensional isogeometric boundary element method for elastostatic analysis. Computer Methods Appl Mech Eng (2012) 209-212:87–100. doi:10.1016/j.cma.2011.08.008

CrossRef Full Text | Google Scholar

38. Hughes T, Cottrell J, Bazilevs Y. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Computer Methods Appl Mech Eng (2005) 194:4135–95. doi:10.1016/j.cma.2004.10.008

CrossRef Full Text | Google Scholar

39. Chen L, Lu C, Lian H, Liu Z, Zhao W, Li S, et al. Acoustic topology optimization of sound absorbing materials directly from subdivision surfaces with isogeometric boundary element methods. Computer Methods Appl Mech Eng (2020) 362:112806. doi:10.1016/j.cma.2019.112806

CrossRef Full Text | Google Scholar

40. Shen X, Du C, Jiang S, Sun L, Chen L. Enhancing deep neural networks for multivariate uncertainty analysis of cracked structures by POD-RBF. Theor Appl Fracture Mech (2023) 125:103925. doi:10.1016/j.tafmec.2023.103925

CrossRef Full Text | Google Scholar

41. Scott M, Simpson R, Evans J, Lipton S, Bordas S, Hughes T, et al. Isogeometric boundary element analysis using unstructured T-splines. Computer Methods Appl Mech Eng (2013) 254:197–221. doi:10.1016/j.cma.2012.11.001

CrossRef Full Text | Google Scholar

42. Takahashi T, Matsumoto T. An application of fast multipole method to isogeometric boundary element method for Laplace equation in two dimensions. Eng Anal Boundary Elem (2012) 36:1766–75. doi:10.1016/j.enganabound.2012.06.004

CrossRef Full Text | Google Scholar

43. Chen L, Cheng R, Li S, Lian H, Zheng C, Bordas S. A sample-efficient deep learning method for multivariate uncertainty qualification of acoustic-vibration interaction problems. Computer Methods Appl Mech Eng (2022) 393:114784. doi:10.1016/j.cma.2022.114784

CrossRef Full Text | Google Scholar

44. Zhang S, Yu B, Chen L. Non-iterative reconstruction of time-domain sound pressure and rapid prediction of large-scale sound field based on ig-drbem and pod-rbf. J Sound Vibration (2024) 573:118226. doi:10.1016/j.jsv.2023.118226

CrossRef Full Text | Google Scholar

45. Chen L, Lian H, Dong HW, Yu P, Jiang S, Bordas SP. Broadband topology optimization of three-dimensional structural-acoustic interaction with reduced order isogeometric FEM/BEM. J Comput Phys (2024) 509:113051. doi:10.1016/j.jcp.2024.113051

CrossRef Full Text | Google Scholar

46. Chen L, Lian H, Pei Q, Meng Z, Jiang S, Dong HW, et al. FEM-BEM analysis of acoustic interaction with submerged thin-shell structures under seabed reflection conditions. Ocean Eng (2024) 309:118554. doi:10.1016/j.oceaneng.2024.118554

CrossRef Full Text | Google Scholar

47. Cao G, Yu B, Chen L, Yao W. Isogeometric dual reciprocity bem for solving non-fourier transient heat transfer problems in fgms with uncertainty analysis. Int J Heat Mass Transfer (2023) 203:123783. doi:10.1016/j.ijheatmasstransfer.2022.123783

CrossRef Full Text | Google Scholar

48. Ginnis A, Kostas K, Politis C, Kaklis P, Belibassakis K, Gerostathis T, et al. Isogeometric boundary-element analysis for the wave-resistance problem using T-splines. Computer Methods Appl Mech Eng (2014) 279:425–39. doi:10.1016/j.cma.2014.07.001

CrossRef Full Text | Google Scholar

49. Zhang G, He Z, Qin J, Hong J. Magnetically tunable bandgaps in phononic crystal nanobeams incorporating microstructure and flexoelectric effects. Appl Math Model (2022) 111:554–66. doi:10.1016/j.apm.2022.07.005

CrossRef Full Text | Google Scholar

50. Zhang G, He Z, Gao XL, Zhou H. Band gaps in a periodic electro-elastic composite beam structure incorporating microstructure and flexoelectric effects. Archive Appl Mech (2023) 93:245–60. doi:10.1007/s00419-021-02088-9

CrossRef Full Text | Google Scholar

51. Chen L, Lian H, Xu Y, Li S, Liu Z, Atroshchenko E, et al. Generalized isogeometric boundary element method for uncertainty analysis of time-harmonic wave propagation in infinite domains. Appl Math Model (2023) 114:360–78. doi:10.1016/j.apm.2022.09.030

CrossRef Full Text | Google Scholar

52. Zhang G, Gao X, Wang S, Hong J. Bandgap and its defect band analysis of flexoelectric effect in phononic crystal plates. Eur J Mechanics-A/Solids (2024) 104:105192. doi:10.1016/j.euromechsol.2023.105192

CrossRef Full Text | Google Scholar

53. Zhang G, He Z, Wang S, Hong J, Cong Y, Gu S. Elastic foundation-introduced defective phononic crystals for tunable energy harvesting. Mech Mater (2024) 191:104909. doi:10.1016/j.mechmat.2024.104909

CrossRef Full Text | Google Scholar

54. Shen X, Du C, Jiang S, Zhang P, Chen L. Multivariate uncertainty analysis of fracture problems through model order reduction accelerated SBFEM. Appl Math Model (2024) 125:218–40. doi:10.1016/j.apm.2023.08.040

CrossRef Full Text | Google Scholar

55. Simpson R, Liu Z, Vázquez R, Evans J. An isogeometric boundary element method for electromagnetic scattering with compatible B-spline discretizations. J Comput Phys (2018) 362:264–89. doi:10.1016/j.jcp.2018.01.025

CrossRef Full Text | Google Scholar

56. Qu Y, Jin F, Yang J. Temperature effects on mobile charges in thermopiezoelectric semiconductor plates. Int J Appl Mech (2021) 13:2150037. doi:10.1142/s175882512150037x

CrossRef Full Text | Google Scholar

57. Xu Y, Li H, Chen L, Zhao J, Zhang X. Monte Carlo based isogeometric stochastic finite element method for uncertainty quantization in vibration analysis of piezoelectric materials. Mathematics (2022) 10:1840. doi:10.3390/math10111840

CrossRef Full Text | Google Scholar

58. Qu Y, Zhang G, Gao X, Jin F. A new model for thermally induced redistributions of free carriers in centrosymmetric flexoelectric semiconductor beams. Mech Mater (2022) 171:104328. doi:10.1016/j.mechmat.2022.104328

CrossRef Full Text | Google Scholar

59. Chen L, Wang Z, Lian H, Ma Y, Meng Z, Li P, et al. Reduced order isogeometric boundary element methods for CAD-integrated shape optimization in electromagnetic scattering. Computer Methods Appl Mech Eng (2024) 419:116654. doi:10.1016/j.cma.2023.116654

CrossRef Full Text | Google Scholar

60. Li H, Chen L, Zhi G, Meng L, Lian H, Liu Z, et al. A direct fe2 method for concurrent multilevel modeling of piezoelectric materials and structures. Computer Methods Appl Mech Eng (2024) 420:116696. doi:10.1016/j.cma.2023.116696

CrossRef Full Text | Google Scholar

61. Chen L, Liu C, Zhao W, Liu L. An isogeometric approach of two dimensional acoustic design sensitivity analysis and topology optimization analysis for absorbing material distribution. Computer Methods Appl Mech Eng (2018) 336:507–32. doi:10.1016/j.cma.2018.03.025

CrossRef Full Text | Google Scholar

62. Xu G, Li M, Mourrain B, Rabczuk T, Xu J, Bordas SP. Constructing IGA-suitable planar parameterization from complex CAD boundary by domain partition and global/local optimization. Computer Methods Appl Mech Eng (2018) 328:175–200. doi:10.1016/j.cma.2017.08.052

CrossRef Full Text | Google Scholar

63. Li S, Trevelyan J, Wu Z, Lian H, Wang D, Zhang W. An adaptive SVD-Krylov reduced order model for surrogate based structural shape optimization through isogeometric boundary element method. Computer Methods Appl Mech Eng (2019) 349:312–38. doi:10.1016/j.cma.2019.02.023

CrossRef Full Text | Google Scholar

64. Chen L, Lian H, Liu Z, Chen H, Atroshchenko E, Bordas S. Structural shape optimization of three dimensional acoustic problems with isogeometric boundary element methods. Computer Methods Appl Mech Eng (2019) 355:926–51. doi:10.1016/j.cma.2019.06.012

CrossRef Full Text | Google Scholar

65. Chen L, Lian H, Natarajan S, Zhao W, Chen X, Bordas S. Multi-frequency acoustic topology optimization of sound-absorption materials with isogeometric boundary element methods accelerated by frequency-decoupling and model order reduction techniques. Computer Methods Appl Mech Eng (2022) 395:114997. doi:10.1016/j.cma.2022.114997

CrossRef Full Text | Google Scholar

66. Lu C, Chen L, Luo J, Chen H. Acoustic shape optimization based on isogeometric boundary element method with subdivision surfaces. Eng Anal Boundary Elem (2023) 146:951–65. doi:10.1016/j.enganabound.2022.11.010

CrossRef Full Text | Google Scholar

67. Chen L, Lian H, Liu Z, Gong Y, Zheng C, Bordas S. Bi-material topology optimization for fully coupled structural-acoustic systems with isogeometric FEM-BEM. Eng Anal Boundary Elem (2022) 135:182–95. doi:10.1016/j.enganabound.2021.11.005

CrossRef Full Text | Google Scholar

68. Burton AJ, Miller GF. The application of integral equation methods to the numerical solution of some exterior boundary-value problems. Proc R Soc Lond (1971) 323:201–10. doi:10.1098/rspa.1971.0097

CrossRef Full Text | Google Scholar