Design of shaped piezoelectric modal sensors for cantilever beams with intermediate support by using differential transformation method

In this paper a solution to the problem of finding the shape of piezoelectric modal sensors for a cantilever beam with intermediate support is proposed by using the differential transformation method (DTM). A general expression for designing the shape of a piezoelectric modal sensor is presented, in which the output signal of the designed sensor is proportional to the response of the target mode. Other modes are filtered out. The modal sensor shape is expressed as a linear function of the second spatial derivative of the structural mode shape function. By using boundary condition and continuity condition equations at intermediate support, the closed-form series solution of the second spatial derivative of the mode shapes can be determined based on DTM. Then the shapes of the designed modal sensors are obtained. Finally, numerical examples are given to demonstrate the feasibility of the proposed modal sensors for the cantilever beam with intermediate support.     

任意边界板结构的声辐射模态误差传感器的设计

获取声辐射模态伴随系数是基于声辐射模态理论进行主动结构声控制(ASAC)的重要环节。以往PVDF分布式传感器的设计难点是振速展开受边界条件的限制,其设计过程往往是针对特定边界条件展开的。本文在声辐射模态理论和两维分布式传感器的压电方程的基础上,将板表面振速分布用Legendre多项式展开,给出了两维板结构的PVDF传感器形状与边界条件无关的设计方法。这样设计得到的传感器能应用于任意边界条件和任意振速分布的两维板结构,且实时性好,拓宽了其应用范围。本文还分别以固定边界条件板及在该板中任取一小区域两种情况为例,证明了该设计方法的可行性。     

Free vibration analysis of elastically connected multiple-beams by using the Adomian modified decomposition method

The Adomian modified decomposition method (AMDM) is employed in this paper to investigate the free vibrations of N elastically connected parallel Euler–Bernoulli beams, which are continuously joined by a Winkler-type elastic layer. The proposed AMDM method can be used to analyze the vibration of beam system consisting of an arbitrary number of beams. By using boundary conditions the natural frequencies and corresponding mode shapes can be easily obtained simultaneously. The numerical results for different boundary conditions, beam numbers and the stiffness of the Winkler-type elastic layer are presented. It is shown that the AMDM offers an accurate and effective method of free vibration analysis of multiple-connected beams with arbitrary boundary conditions.     

Boundary conditions for the bending of a piezoelectric beam

For beam bending in transversely isotropic piezoelectric media, the reciprocal theorem and the general solution of piezoelasticity are applied in a novel way to obtain the appropriate stress and mixed boundary conditions accurate to all orders for the beam of general edge geometry and loadings. By generalizing the method developed by Gregory and Wan, a set of necessary conditions on the edge-data for the existence of a rapidly decaying solution is established. The prescribed edge-data of the beam must satisfy these conditions in order that they could generate a decaying state within the beam. When stress and mixed conditions are imposed on the beam edge, these decaying state conditions for the case of bending deformation of piezoelectric beam are derived explicitly. They are then used for the correct formulation of boundary conditions for the beam theory solution (or the interior solution). Besides, an analytical solution of elastic beam is formulated to verify validity of our boundary conditions. For the stress data, our boundary conditions coincide with those obtained in conventional forms of beam theories. More importantly, the appropriate boundary conditions with two sets of mixed edge-data are obtained for the first time.     

Implicit spectrally-accurate method for moving boundary problems using immersed boundary conditions concept

A fully implicit, spectral algorithm for the analysis of moving boundary problem is described. The algorithm is based on the concept of immersed boundary conditions (IBC), i.e., the computational domain is fixed while the time dependent physical domain is submerged inside the computational domain, and is described in the context of the diffusion-type problems. The physical conditions along the edges of the physical domain are treated as internal constraints. The method eliminates the need for adaptive grid generation that follows evolution of the physical domain and provides sharp resolution of the location of the boundary. Various tests confirm the spectral accuracy in space and the first- and second-order accuracy in time. The computational cost advantage of the IBC method as compared with the more traditional algorithm based on the mapping concept is demonstrated.     

Complex FEM modal solver of optical waveguides with PML boundary conditions

A full-wave modal analysis of two-dimensional, lossy and anisotropic optical waveguides using the finite element method (FEM) is presented. In order to describe the behavior of radiating fields, anisotropic perfectly matched layer boundary conditions are applied for the first time in modal solvers. The approach has been implemented using high order edge elements. The resulting sparse eigenvalue algebraic problem is solved through the Arnoldi method. Application to an antiresonant reflecting optical waveguide is reported.     

消声器传递损失预测的边界元模态展开混合方法

将边界元法和解析方法结合形成一种混合方法用于计算消声器的传递损失,消声器被划分成若干个子结构,解析方法和边界元方法被分别用于计算规则结构和不规则结构的阻抗矩阵,不同子结构之间通过阻抗矩阵连接起来。为减少计算时间,采用一种基于模态配点法的简化方法。对单级膨胀腔、双级膨胀腔和穿孔管阻性消声器的传递损失进行了计算,混合方法计算结果与解析方法和三维数值方法计算结果吻合良好。分析了混合方法的计算效率并与传统子结构方法进行了比较,混合方法能明显节省计算时间。     

消声器传递损失预测的边界元模态展开混合方法

将边界元法和解析方法结合形成一种混合方法用于计算消声器的传递损失,消声器被划分成若干个子结构,解析方法和边界元方法被分别用于计算规则结构和不规则结构的阻抗矩阵,不同子结构之间通过阻抗矩阵连接起来。为减少计算时间,采用一种基于模态配点法的简化方法。对单级膨胀腔、双级膨胀腔和穿孔管阻性消声器的传递损失进行了计算,混合方法计算结果与解析方法和三维数值方法计算结果吻合良好。分析了混合方法的计算效率并与传统子结构方法进行了比较,混合方法能明显节省计算时间。     

有限元+边界元——修正的模态分解法预报水下旋转薄壳的振动和声辐射

本文提出了一种有限元+边界元——修正的模态分解法用于预报水下结构的振动和声辐射问题.通过引入结构的剩余模态计入被忽略的高阶模态的准静态响应,使模态分解法的收敛性大大改善.文中以轴对称力作用下水下弹性薄球壳振动和声辐射问题为例进行了数值计算,给出了ka=1,2,3和4时球壳表面位移、声压和远场辐射声的计算结果.与原有的方法比较可见,该方法提高了模态分解法的收敛性,特别对远场辐射声指向性尤为明显。     

A method for determining foil thicknesses in TEM by using convergent beam electron diffraction under weak beam conditions

In the present work, convergent beam electron diffraction was studied in zirconium (a material of intermediate atomic number) at 300 keV, under weak beam diffraction conditions. For a particular thickness, the details in an observed low order disc were matched to those calculated using the multibeam dynamical theory. This presents the possibility of determining foil thickness over a wide range, with an estimated experimental accuracy of ≈7% or less. In contrast to other convergent beam techniques, the present method, which uses weak beam conditions, can employ commonly-occurring low order reflections to extract thicknesses.

A simple equation based on the two beam approximation, is derived to determine foil thickness (to within ≈ 10%) without resorting to detailed image matching. This equation can be used for a rough estimate of foil thickness while carrying out TEM observations.     

不同边界条件下的封闭矩形声腔的结构-声耦合分析

分析了由2块四边弹性支承的弹性板及4块刚性板构成的封闭矩形声腔的结构-声耦合特性。通过在弹性板边界施加假想的连续分布弹簧系统模拟板的不同边界和连接条件,利用汉密尔顿函数和瑞利-李兹方法,充分考虑了两弹性板之间以及弹性板与腔体内声场之间的耦合。腔体内的声压响应和弹性板的振动速度与Kim的实验结果能较好的吻合,验证了本文理论推导的正确性。最后,通过对耦合声场的分析计算,表明两弹性板之间为弱耦合,耦合声场主要由受到外激励的结构所决定;支承板的线弹簧的刚度变化对耦合声场的影响较旋转弹簧大。     

Elliptic PDE formulation and boundary conditions of the spherical harmonics method of arbitrary order for general three-dimensional geometries

The inherent complexity of the radiative transfer equation makes the exact treatment of radiative heat transfer impossible even for idealized situations and simple boundary conditions. Therefore, a wide variety of efficient solution methods have been developed for the RTE. Among these solution methods the spherical harmonics method, the moment method, and the discrete ordinates method provide means to obtain higher-order approximate solutions to the equation of radiative transfer. Although the assembly of the governing equations for the spherical harmonics method requires tedious algebra, their final form promises great accuracy for any given order, since it is a spectral method (rather than finite difference/finite volume in the case of discrete ordinates). In this study, a new methodology outlined in a previous paper on the spherical harmonics method (P_N) is further developed. The new methodology employs successive elimination of spherical harmonic tensors, thus reducing the number of first-order partial differential equations needed to be solved simultaneously by previous P_N approximations (=(N+1)²). The result is a relatively small set (=N(N+1)/2) of second-order, elliptic partial differential equations, which can be solved with standard PDE solution packages. General boundary conditions and supplementary conditions using rotation of spherical harmonics in terms of local coordinates are formulated for the general P_N approximation for arbitrary three-dimensional geometries. Accuracy of the P_N approximation can be further improved by applying the “modified differential approximation” approach first developed for the P₁-approximation. Numerical computations are carried out with the P₃ approximation for several new two-dimensional problems with emitting, absorbing, and scattering media. Results are compared to Monte Carlo solutions and discrete ordinates simulations and a discussion of ray effects and false scattering is provided.     

Numerical method of designing three-dimensional open cloaks with arbitrary boundary shapes

In this paper, three-dimensional (3-D) open cloaks were designed based on the coordinate transformation method. When the transformation center point is close to the boundary of the cloaks, the material parameter tensors in the near-boundary areas of the cloaks approximate to those in the background, so the near-boundary areas can be designed to be open. Full wave simulations based on finite element method verified the open cloaks we designed. Open cloaks can exchange information and materials with background media. Using the design method proposed in this paper, the boundary of open cloaks can be arbitrarily shaped, which greatly enhances the applicability of open cloaks.     

Free vibration of an embedded single-walled carbon nanotube with various boundary conditions using the RMVT-based nonlocal Timoshenko beam theory and DQ method

The nonlocal Timoshenko beam theories (TBTs), based on the Reissner mixed variation theory (RMVT) and principle of virtual displacement (PVD), are derived for the free vibration analysis of a single-walled carbon nanotube (SWCNT) embedded in an elastic medium and with various boundary conditions. The strong formulations of the nonlocal TBTs are derived using Hamilton's principle, in which Eringen's nonlocal constitutive relations are used to account for the small-scale effect. The interaction between the SWCNT and its surrounding elastic medium is simulated using the Winkler and Pasternak foundation models. The frequency parameters of the embedded SWCNT are obtained using the differential quadrature (DQ) method. In the cases of the SWCNT without foundations, the results of RMVT- and PVD-based nonlocal TBTs converge rapidly, and their convergent solutions closely agree with the exact ones available in the literature. Because the highest order with regard to the derivatives of the field variables used in the RMVT-based nonlocal TBT is lower than that used in its PVD-based counterpart, the former is more efficient than the latter with regard to the execution time. The former is thus both faster and obtains more accurate solutions than the latter for the numerical analysis of the embedded SWCNT.     

Off-diagonal Bethe ansatz solution of the XXX spin chain with arbitrary boundary conditions

Employing the off-diagonal Bethe ansatz method proposed recently by the present authors, we exactly diagonalize the XXX spin chain with arbitrary boundary fields. By constructing a functional relation between the eigenvalues of the transfer matrix and the quantum determinant, the associated T–Q$T\text{–}Q$ relation and the Bethe ansatz equations are derived.     

A unified approach for predicting sound radiation from baffled rectangular plates with arbitrary boundary conditions

This paper discusses sound radiation from a baffled rectangular plate with each of its edges arbitrarily supported in the form of elastic restraints. The plate displacement function is universally expressed as a 2-D Fourier cosine series supplemented by several 1-D series. The unknown Fourier expansion coefficients are then determined by using the Rayleigh-Ritz procedure. Once the vibration field is solved, the displacement function is further simplified to a single standard 2-D Fourier cosine series in the subsequent acoustic analysis. Thus, the sound radiation from a rectangular plate can always be obtained from the radiation resistance matrix for an invariant set of cosine functions, regardless of its actual dimensions and boundary conditions. Further, this radiation resistance matrix, unlike the traditional ones for modal functions, only needs to be calculated once for all plates with the same aspect ratio. In order to determine the radiation resistance matrix effectively, an analytical formula is derived in the form of a power series of the non-dimensional acoustic wavenumber; the formula is mathematically valid and accurate for any wavenumber. Several numerical examples are presented to validate the formulations and show the effect of the boundary conditions on the radiation behavior of planar sources.     

Non-conformal domain decomposition method with second-order transmission conditions for time-harmonic electromagnetics

Non-overlapping domain decomposition (DD) methods provide efficient algorithms for solving time-harmonic Maxwell equations. It has been shown that the convergence of DD algorithms can be improved significantly by using high order transmission conditions. In this paper, we extend a newly developed second-order transmission condition (SOTC), which involves two second-order transverse derivatives, to facilitate fast convergence in the non-conformal DD algorithms. However, the non-conformal nature of the DD methods introduces an additional technical difficulty, which results in poor convergence in many real-life applications. To mitigate the difficulty, a corner-edge penalty method is proposed and implemented in conjunction with the SOTC to obtain truly robust solver performance. Numerical results verify the analysis and demonstrate the effectiveness of the proposed methods on a few model problems. Finally, drastically improved convergence, compared to the conventional Robin transmission condition, was observed for an electrically large problem of practical interest.