The dynamic analysis of a cracked Timoshenko beam by the spectral element method

The aim of this paper is to introduce a new finite spectral element of a cracked Timoshenko beam for modal and elastic wave propagation analysis. The proposed approach deals with the spectral element method. This method is suitable for analyzing wave propagation problems as well as for calculating modal parameters of the structure. In the paper, the results of the change in modal parameters due to crack appearance are presented. The influence of the crack parameters, especially of the changing location of the crack, on the wave propagation is examined. Responses obtained at different points of the beam are presented. Proper analysis of these responses allows one to indicate the crack location in a very precise way. This fact is very promising for the future work in the damage detection field.     

渐变折射率波导传播特性微扰计算的代数递推公式

应用量子力学的超位力定理（HVT）和赫尔曼－费恩曼（Hellmann-Feynman）定理（HFT），导出了用微扰法计算渐变折射率平面波导传播常数各级近似值的代数递推公式，并给出了由传播常数近似值计算模场近似表达式的代数式。计算了几个典型实例，此法的特点是便于直接计算任意级近似的传播常数及模场分布，且可推广到二维波导的计算。     

水平变化双层波导声传播问题的耦合简正波解

为了验证现有模型的精度,导出了全反射下边界双层波导中简正波耦合矩阵的解析表达式,并将其应用到全局矩阵耦合简正波模型(Direct Global Matrix Coupled-Mode)中,使得该模型可以提供水平变化双层波导问题的标准解。文中首先利用COUPLE的简正波及耦合矩阵数值解验证了该简正波及耦合矩阵解析表达式的正确性;其次,采用改进的DGMCM模型求解了双层波导海山声传播损失,结果表明,改进后的DGMCM模型可以非常精确地求解水平变化双层波导问题,可作为求解此类问题的标准模型使用。     

Sound propagation in inhomogeneous waveguides with sound-speed profiles using the multimodal admittance method

The multimodal admittance method and its improvement are presented to deal with various aspects in underwater acoustics, mostly for the sound propagation in inhomogeneous waveguides with sound-speed profiles, arbitrary-shaped liquid-like scatterers, and range-dependent environments. In all cases, the propagation problem governed by the Helmholtz equation is transformed into initial value problems of two coupled first-order evolution equations with respect to the modal components of field quantities(sound pressure and its derivative), by projecting the Helmholtz equation on a constructed orthogonal and complete local basis. The admittance matrix, which is the modal representation of Direchlet-to-Neumann operator, is introduced to compute the first-order evolution equations with no numerical instability caused by evanescent modes. The fourth-order Magnus scheme is used for the numerical integration of differential equations in the numerical implementation. The numerical experiments of sound field in underwater inhomogeneous waveguides generated by point sources are performed. Besides, the numerical results computed by simulation software COMSOL Multiphysics are given to validate the correction of the multimodal admittance method. It is shown that the multimodal admittance method is an efficient and stable numerical method to solve the wave propagation problem in inhomogeneous underwater waveguides with sound-speed profiles, liquid-like scatterers, and range-dependent environments. The extension of the method to more complicated waveguides such as horizontally stratified waveguides is available.     

Full-vector analysis of a realistic photonic crystal fiber

We analyze the guiding problem in a realistic photonic crystal fiber, using a novel full-vector modal technique. This is a biorthogonal modal method based on the non-self-adjoint character of the electromagnetic propagation in a fiber. Dispersion curves of guided modes for different fiber structural paremeters are calculated, along with the two-dimensional transverse intensity distribution of the fundamental mode. Our results match those achieved in recent experiments in which the feasibility of this type of fiber was shown.     

An accurate variational analysis of single-mode diffused channel waveguides

We present a procedure based on the variational principle to obtain the modal propagation characteristics of single-mode diffused channel waveguides with arbitrary index profiles. The only assumption in our method, in addition to the scalar approximation, is that the trial field is separable along the depth and the width directions. The method is iterative, with each cycle of iteration requiring the solution of two planar waveguide propagation problems. The convergence is usually achieved within 2–3 cycles of iteration. Comparisons with other methods show that our method gives much better accuracy.     

ON THE RELATION BETWEEN COMPLEX MODES AND WAVE PROPAGATION PHENOMENA

This paper discusses the well-known, but often misunderstood, concept of complex modes of dynamic structures. It shows how complex modes can be interpreted in terms of wave propagation phenomena caused by either localized damping or propagation to the surrounding media. Numerical simulation results are presented for different kinds of structures exhibiting modal and wave propagation characteristics: straight beams, an L-shaped beam, and a three-dimensional frame structure. The input/output transfer relations of these structures are obtained using a spectral formulation known as the spectral element method (SEM). With this method, it is straightforward to use infinite elements, usually known as throw-off elements, to represent the propagation to infinity, which is a possible cause of modal complexity. With the SEM model, the exact dynamic behavior of structures can be investigated. The mode complexity of these structures is investigated. It is shown that mode complexity characterizes a behavior that is half-way between purely modal and purely propagative. A coefficient for quantifying mode complexity is introduced. The mode complexity coefficient consists of the correlation coefficient between the real and imaginary parts of the eigenvector, or of the operational deflection shape (ODS). It is shown that, far from discontinuities, this coefficient is zero in the case of pure wave propagation in which case the plot of the ODS in the complex plane is a perfect circle. In the other extreme situation, a finite structure without damping (or with proportional damping), where the mode shape (or the ODS) is a straight line on the complex plane, has a unitary complexity coefficient. For simple beam structures, it is shown that the mode complexity factor can also be calculated by curve-fitting the mode to an ellipse and computing the ratio of its radii.     

共振域矩形介质光栅的逆向设计

采用模态法研究了共振域光栅分束器及闪耀光栅等光学器件的设计原理及逆向设计方法,给出了-1级非偏振闪耀光栅的典型设计示例。研究表明,自准直角入射时,光束在光栅内传输将激发产生离散模,通过调节离散模透过光栅传输时的累积相位差,可对衍射光进行调控,实现偏振(非偏振)分束、偏振(非偏振)闪耀等功能。采用严格耦合波法对基于模态法逆向设计的-1级非偏振石英闪耀光栅的几何参数进行验证,计算结果显示两者非常吻合,证明了该方法的正确性。     

共振域矩形介质光栅的逆向设计

采用模态法研究了共振域光栅分束器及闪耀光栅等光学器件的设计原理及逆向设计方法,给出了-1级非偏振闪耀光栅的典型设计示例。研究表明,自准直角入射时,光束在光栅内传输将激发产生离散模,通过调节离散模透过光栅传输时的累积相位差,可对衍射光进行调控,实现偏振(非偏振)分束、偏振(非偏振)闪耀等功能。采用严格耦合波法对基于模态法逆向设计的-1级非偏振石英闪耀光栅的几何参数进行验证,计算结果显示两者非常吻合,证明了该方法的正确性。     

All propagation modes of large-core multimode optical fibers with an arbitrary core profile

Plastic optical fibers are typical large-core multimode optical fibers. High-bandwidth graded-index plastic optical fibers, which support a great number of propagation modes compared with conventional silica-glass multimode optical fibers, were developed in the 1990s. However, because they support a great number of propagation modes, their modal analyses have been limited to the WKB analysis. We obtain all the propagation modes of a convex-index large-core multimode optical fiber by use of the finite-element method, which has a strong advantage for arbitrary core profiles.     

Modal Doppler theory of an arbitrarily accelerating continuous-wave source applied to mode extraction in the oceanic waveguide

A Doppler-based method for using a moving narrow-band source to extract the modes of acoustic propagation in a range-independent shallow ocean waveguide over a partial-water-column spanning vertical line array (VLA) is introduced. Because the modal components propagate at distinct frequencies in the case of uniform radial source motion, the modal depth functions may be isolated and extracted from a frequency decomposition of the field. Because Doppler broadening due to radial source accelerations degrades the effectiveness of the extraction method, the method incorporates a technique to compensate for Doppler broadening. As the basis for the compensation technique, a theory is introduced for describing the VLA field from an accelerating cw source. By connecting the range of the source at the time a signal feature is emitted (the retarded time) to the range of the source at the time the signal feature arrives at the receiver (the contemporary time), the theory incorporates the Doppler effects associated with the finite group velocities of the modal components. The mode extraction method and compensation technique are applied to simulation and ocean data.     

分层半空间表面非线性瑞利波的激发*

针对非线性瑞利波在均匀分层半空间结构中的激发和传播规律进行研究。根据摄动理论和模态分解将分层半空间结构中瑞利波的二次谐波声场表示为二倍频瑞利波模式的线性组合,经由互易关系得到各模式的展开系数表达式。对不同分层半空间结构中瑞利波二次谐波的激发和传播特性进行讨论,结果表明相速度匹配的瑞利波模式其二次谐波分量随传播距离线性增长,非匹配模式的二次谐波分量则沿传播方向周期震荡传播。此外,文中定义非线性参数表征瑞利波模式产生的非线性程度,这有利于选择出具有明显非线性效应的匹配点,为后续检测工作提供理论依据,具有指导意义。     

Eigenmodes of index-modulated layers with lateral PMLs

Maxwell equations are solved in a layer comprising a finite number of homogeneous isotropic dielectric regions ended by anisotropic perfectly matched layers (PMLs). The boundary-value problem is solved and the dispersion relation inside the PML is derived. The general expression of the eigenvalues equation for an arbitrary number of regions in each layer is obtained, and both polarization modes are considered. The modal functions of a single layer ended by PMLs are found, and their orthogonality relation is derived. The present method is useful to simulate scattering problems from dielectric objects as well as propagation in planar slab waveguides. Its potential to deal with more complex problems such as the scattering from an object with arbitrary cross section in open space using the multilayer modal method is briefly discussed.     

光子晶体光纤的正交函数模型

提出了一种用于分析光子晶体光纤的正交函数模型。采用一种新型超格子的构造方法,将光子晶体光纤的横向介电常量表示为两种周期性结构叠加。将横向电场以厄密一高斯函数展开,利用正交函数的性质,将全矢量波动方程转化为矩阵本征值问题,求解本征值问题可得到模式的传输常量及模场分布。利用此模型举例讨论了椭圆孔三角格子光子晶体光纤的模场分布和偏振特性以及三角格子光子晶体光纤的色散特性和有效面积等传输特性。作为一种普适的模型,此方法还可适用于四方结构、蜂窝结构及椭圆孔等多种结构光子晶体光纤。     

基于变量变换伽辽金法光波导半矢量分析

基于变量变换伽辽金法(VTGM),建立了半矢量分析光波导本征值问题的数理模型;获得了矩形光波导、脊形光波导的本征模的模场分布及归一化传播常数,计算结果和已发表的结果非常接近.半矢量VTGM采用平面映射边界条件,避免了边界截断,又能反映模场的偏振特性,是一种较为理想的光子器件分析方法.     

Solutions of the full- and quasi-vector wave equations based on H-field for optical waveguides by using mapped Galerkin method

The mapped Galerkin method in solving the full-vector and quasi-vector wave equations in terms of transverse magnetic fields (H-formulation) for optical waveguides with step-index profiles is described. By transforming the whole x-y space onto a unit square and using two-dimensional Fourier series expansion, the modal distributions and propagation constants for optical waveguides are obtained in the absence of boundary truncation. Results for step-index circular fiber, buried rectangular waveguide, and optical rib waveguide are presented. Solutions are good agreed with exact solutions and numerical results by using vector nonlinear iterative method, Fourier operator transform method, and vector beam propagation method.     

发动机圆形管道内噪声衰减谱模态分析

本文应用模态分析方法建立了剪切流存在条件下,发动机多段声衬圆形管道声传播工程计算模型,对管内各模态频谱和总噪声衰减频谱进行了算例计算,并与有关文献试验数据进行了对比。结果表明,多段声衬圆形管道中声传播工程计算方法是可行的,从而为发动机前短舱管内声传播研究提供了一种模态分析工程预测方法。     

Photonic crystal fibers with material anisotropy

In this paper we are modeling the interplay of material and form birefringence in photonic crystal fibers. We introduce an efficient numerical method for the calculation of the modal structure. Our approach relies solving the fully vectorial wave equation for the transverse magnetic field and the respective propagation constants using a plane wave expansion. The method accounts for a simple form of material anisotropy. Our analysis is relevant to certain application areas, and in particular to fiber sensing, where material birefringence arises for instance due to transversally applied mechanical stress. We analyze the influence of material birefringence on the modal birefringence and the state of polarization of the fundamental mode.     

Lamb waves propagation in elastic plane layers with a joint strip

The Lamb waves are used for the ultrasonic characterization of welds because of their relative long-range propagation. In this paper, a simplified model of a weld-strip between two identical semi-infinite elastic layers is investigated. The reflected and transmitted ultrasonic fields are expressed by modal series whose coefficients are obtained by application of orthogonality relation. Comparisons with solutions obtained by finite elements wave propagation simulations are made. The energy balance between the incident and the scattered waves is also used to verify the accuracy of the obtained modal amplitudes.