三阶矢量有限元方法精确计算轴对称谐振腔高阶模

为了进一步提高数值求解谐振腔高阶模的精度，本文提出了三阶矢量有限元方法，并针对二阶矢量有限元轴对称谐振腔高阶模计算程序Cafe对曲线边界的计算能力较差和计算速度较慢的缺点做了改进. 在这些改进的基础上编制了三阶矢量有限元轴对称谐振腔高阶模计算程序meshmatrix3，得到了很好的结果.     

基于高阶叠层矢量基函数的E-H时域有限元方法分析谐振腔和波导结构

将高阶叠层矢量基函数用于E-H时域有限元方法,电场和磁场用相同的基函数展开并同时求解,时间离散采用Crank-Nicolson差分格式使得时间步长的选取摆脱稳定性条件的限制,同时采用完美匹配层来截断计算区域.对三维谐振腔及波导结构进行数值模拟与分析,结果表明,相较于低阶基函数,高阶叠层矢量基函数可以有效提高E-H时域有限元方法的计算精度.     

高阶间断伽辽金时域有限元方法分析三维谐振腔

从一阶麦克斯韦旋度方程出发,研究一种区域分解时域有限元目的——高阶间断伽辽金时域有限元目的.其中对时间的离散采用Crank-Nicolson差分格式,电场和磁场采用相同阶数的高阶矢量基函数展开.分析三维谐振腔问题,数值结果表明,目的 中时间步长的选取可以摆脱CFL稳定性条件的限制;此外,与基于常用Whitney矢量基函数的目的 相比,采用高阶矢量基函数可以明显地提高计算精度及计算效率.     

基于Fox-Li迭代法的偏振敏感谐振腔EI北大核心CSCD

矢量偏振光束的发展对偏振敏感谐振腔的模式计算提出了新的要求,基于Fox-Li迭代法和琼斯矢量理论,提出了用于偏振敏感谐振腔的矢量Fox-Li迭代法。该方法可以计算偏振敏感谐振腔内的矢量偏振模式,包括角向偏振TE01*模和径向偏振TM01*模。通过Matlab编程对球面偏振敏感谐振腔内的矢量偏振模式进行数值计算,计算结果与理论一致,证明了该方法的正确性。然后基于轴快流(FAF)CO2激光器平台,采用组合轴锥镜作为谐振腔的尾镜,得到了2.17kW的角向偏振模式输出,与矢量Fox-Li迭代法仿真结果吻合,实验证明了该方法的有效性和准确性。矢量Fox-Li迭代法对于偏振敏感谐振腔的分析和设计具有重要的指导意义。     

基于Fox-Li迭代法的偏振敏感谐振腔

矢量偏振光束的发展对偏振敏感谐振腔的模式计算提出了新的要求,基于Fox-Li迭代法和琼斯矢量理论,提出了用于偏振敏感谐振腔的矢量Fox-Li迭代法。该方法可以计算偏振敏感谐振腔内的矢量偏振模式,包括角向偏振TE01*模和径向偏振TM01*模。通过Matlab编程对球面偏振敏感谐振腔内的矢量偏振模式进行数值计算,计算结果与理论一致,证明了该方法的正确性。然后基于轴快流(FAF)CO2激光器平台,采用组合轴锥镜作为谐振腔的尾镜,得到了2.17kW的角向偏振模式输出,与矢量Fox-Li迭代法仿真结果吻合,实验证明了该方法的有效性和准确性。矢量Fox-Li迭代法对于偏振敏感谐振腔的分析和设计具有重要的指导意义。     

基于曲四面体单元剖分的CN E-H TDFEM方法

从采用Crank-Nicolson差分格式的基于麦克斯韦旋度方程的E-H时域有限元方法出发,将展开电场和磁场的叠层矢量基函数与曲四面体单元相结合.对金属球谐振腔及介质填充的圆柱形谐振腔的数值模拟表明：相较于规则四面体单元,在剖分单元数目相同的情况下,曲四面体单元离散表面弯曲结构可以获得更高的计算精度.同时,与0.5阶基函数结合曲四面体单元相比,1.0阶基函数与曲四面体单元结合可以用更少的单元数及未知量数目来获得更高的计算精度.     

基于有限元数值模拟永磁轨道上方高温超导块材悬浮力

采用有有限元计算方法和有限差分法对高温超导块材电磁场进行模拟,计算了永磁轨道上方超导块材所受到的磁悬浮力.文中给出了一种基于有限元的3D-模型数值计算方法；模型采用了电磁场磁矢量法(H-方法)；模型中采用Kim模型来描述高温超导块材的临界电流密度Jc；采用E-J幂指数来描述超导块材的本构关系.该计算模型采用Fortra...     

有源腔模式的有限元数值矩阵计算法

提出了一种计算有源腔模式分布的有限元数值矩阵方法。对腔镜进行基于菲涅耳–基尔霍夫衍射积分定理的有限元划分,利用MATLAB计算出光一次往返的传输矩阵。然后根据自再现理论,求出矩阵的特征值和特征向量。从而确定谐振腔的模式分布,给新型谐振腔的研究带来极大方便。     

有效抑制光子晶体加载矩形谐振腔中模式竞争的方法

研究了光子晶体加载矩形谐振腔,针对传统光子晶体加载无法完全抑制模式竞争的问题,提出了特性阻抗和介质加载方法进行改进,完全实现了谐振腔的高次模式的单模工作,为抑制高次模式工作的谐振腔中模式竞争的问题提供了一种新方法.在此基础之上,采用解析方法以及高频仿真软件HFSS,设计了一个采用TM₅₃₀高次模式工作的传统光子晶体加载的谐振腔,进一步采用了该新方法对其进行改进,并深入分析了采用这种方法抑制模式竞争的物理原因.计算结果表明,新方法能有效抑制模式竞争,实现谐振腔的高次工作模式的单模工作. 关键词： 光子晶体 谐振腔 模式竞争 特性阻抗     

同轴放电环形增益区谐振腔镜面的场分布

 对几类环形增益区谐振腔进行了比较，得出利用谐振腔两镜之一的全工作面输出的新型环形增益区谐振腔是适合的谐振腔。提出了一种计算全工作面输出的同轴放电环形增益区谐振腔镜面场分布的方法，利用几何方法对腔镜镜面进行划分,得出了基于菲涅尔-基尔霍夫衍射积分定理的积分区域，并对该区域进行了有限元划分，应用计算机编程计算出了同轴放电环形增益区谐振腔镜面处的光场分布和相位分布。近场分布在环形区中央部分近似地呈高斯分布，且沿径向出现高阶模场分布特点，沿角向无差别。     

Harmonic analysis of lossy, piezoelectric composite transducers using the plane wave expansion method

Periodic composite ultrasonic transducers offer many advantages but the periodic pillar architecture can give rise to unwanted modes of vibration which interfere with the piston like motion of the fundamental thickness mode. In this paper, viscoelastic loss is incorporated into a three-dimensional plane wave expansion model (PWE) of these transducers. A comparison with experimental and finite element data is conducted and a design to damp out these lateral modes is investigated. Scaling and regularisation techniques are introduced to the PWE method to reduce ill-conditioning in the large matrices which can arise. The identification of the modes of vibration is aided by examining profiles of the displacements, electrical potential and Poynting vector. The dispersive behaviour of a 2-2 composite transducer with high shear attenuation in the passive phase is examined. The model shows that the use of a high shear attenuation filler material improves the frequency band gap surrounding the fundamental thickness mode.     

A combined finite element and modal decomposition method to study the interaction of Lamb modes with micro-defects

This paper presents a combined finite element and modal decomposition method to study the interaction of Lamb waves with damaged area. The finite element mesh is used to describe the region around the defects. On the contrary to other hybrid models already developed, the interaction between Lamb waves and defects is computed in the temporal domain. Then, the modal decomposition method permits to determine the wave reflected and transmitted by the damaged area. Modal analysis allows also identifying the mode conversions induced by the defects. These numerical results agree with previous finite element results concerning the interaction of Lamb modes with a notch. Experiments, carried out with gauged defects on an aluminum plate, are also compared to numerical predictions to validate the simulation. Compared to classical techniques of simulation, this new method allows us to investigate the interaction of Lamb modes generated at high frequency-thickness product with micro-defects as corrosion pitting.     

分析各向异性光波导的一种新系统

报道一种灵活，交互和自动化的系统，它以矢量有限无法为基础，并在ＩＭＳＬ的ＰＤＥ／ｐｒｏｔｒａｎ题解环境中实现，用它分析各向异性光波异，得到了各向异性介矩形导中所用传播模式的数值结果，以及采用传播常数和模场分布二种形式评价这个方法的精度，并与早期的研究结果加以比较。     

低频宽带大功率“星型”柱面发射换能器

提出一种具有水平无指向性、低频、宽带、大功率特点的新结构换能器——“星型”柱面发射换能器,换能器由共享后质量块的6个复合棒换能器按“星型”方式组成柱面结构。采用有限元方法对换能器进行了分析设计并制作了换能器样机。在1—10 kHz内,换能器有3个主要工作模态,前两阶工作模态对换能器工作带宽的展宽有贡献,第三阶模态显示出较强的指向性,水池测试换能器具有1倍频程的工作带宽,最大声源级为199.1 dB。有限元模拟结果和实验结果符合较好。研究结果表明,新结构换能器在满足水平无指向性的前提下,利用复合棒的纵振及辐射面的弯曲振动实现了低频、宽带、大功率的要求,提供了一种设计该特性换能器的新思路。      

Vibration analysis of a new curved beam element

The common practice in developing a locking-free curved beam element is to ensure that its interpolation functions of displacement explicitly satisfy the inextensible bending mode condition for the membrane locking-free instead of the rigid body modes. In this paper, we study the impact of this practice on the dynamic characteristics of a finite element by conducting vibration analysis using our newly developed three-node locking-free curved beam element. In this case, the inextensible bending mode condition is satisfied explicitly, while the rigid body modes are satisfied implicitly to 4th-order accuracy. Numerical and experimental examples show that with the newly developed curved beam element, developed by using the implicit representation of a rigid body mode condition, it is possible to recover the rigid body modes of curved beams with low and medium slenderness ratios. This is even true for cases involving a half-circular element and the vibration of the curved beam is predicted with high accuracy.     

High order shear horizontal modes for minimum remnant thickness

Thickness mapping in aging structures suffering from corrosion is challenging especially when the structure is only partially accessible. In plates the high order shear horizontal guided wave modes all have a cutoff frequency thickness product below which they cannot propagate. This property is potentially attractive to estimate the minimum remnant thickness between two transducers. When using a source and a sensor array it is possible to control the number of modes being excited and the size of the region interrogated by the technique. Finite element simulations were used to show that by exciting multiple guided wave modes simultaneously and identifying which modes are received by a sensor array it is possible to estimate the minimum remaining thickness along the propagation path. Initial experimental results showed excellent agreement with the finite element simulations when the plate is uniform and with a thickness reduction between the source and the sensor arrays the minimum remnant thickness was underestimated by approximately 20%.     

A finite element method for determining the acoustic modes of irregular shaped cavities

A twenty node, isoparametric acoustic finite element model is developed for analysing the acoustic modes of irregular shaped cavities. The element is first used to analyse a rectangular enclosure. This indicates that good accuracy can be obtained if one element is used between each nodal plane. Experiments are described on a model van enclosure. A finite element analysis of the model shows very good agreement with the measured results.     

Coupled plasmon-LO phonon modes and Lindhard-Mermin dielectric function of n-GaAs

The wave vector dependence of coupled plasmon-LO phonon modes is studied with opaque Raman spectroscopy. The dispersion of the coupled modes is directly related to the wavevector dependence of the dielectric function of the electron gas and is influenced by finite temperatures and finite damping. The broadening of the wavevector due to absorption is also important. These effects are explained with the quantum mechanical Lindhard dielectric function for finite temperature including electron scattering in the relaxation time approximation of Mermin.     

Higher order representation of the beam cross section deformation in large displacement finite element analysis

Most existing beam formulations assume that the cross section of the beam remains rigid regardless of the amplitude of the displacement. The absolute nodal coordinate formulation (ANCF); however, allows for the deformation of the cross section and leads to a more general beam models that capture the coupling between different modes of displacement. This paper examines the effect of the order of interpolation on the modes of deformation of the beam cross section using ANCF finite elements. To this end, a new two-dimensional shear deformable ANCF beam element is developed. The new finite element employs a higher order of interpolation, and allows for new cross section deformation modes that cannot be captured using previously developed shear deformable ANCF beam elements. The element developed in this study relaxes the assumption of planar cross section; thereby allowing for including the effect of warping as well as for different stretch values at different points on the element cross section. The displacement field of the new element is assumed to be cubic in the axial direction and quadratic in the transverse direction. Using this displacement field, more expressions for the element extension, shear and the cross section stretch can be systematically defined. The change in the cross section area is measured using Nanson’s formula. Measures of the shear angle, extension, and cross section stretch can also be systematically defined using coordinate systems defined at the element material points. Using these local coordinate systems, expressions for a nominal shear angle are obtained. The differences between the cross section deformation modes obtained using the new higher order element and those obtained using the previously developed lower order elements are highlighted. Numerical examples are presented in order to compare the results obtained using the new finite element and the results obtained using previously developed ANCF finite elements.     

热可调液晶填充微结构聚合物光纤设计及特性分析

通过在微结构聚合物光纤(mPOF)的包层空气孔中填充液晶材料获得了高度热可调的带隙型mPOF. 带隙随温度的增加发生显著的蓝移, 带隙上边界的温度灵敏度可达-5.5 nm/℃. 采用全矢量有限元方法对模场特性以及基模有效模场面积的分析结果表明, 该填充液晶的mPOF在带隙的中心波长附近具有大的有效模场面积, 和相同结构未填充液晶的mPOF连接时具有较高的功率耦合效率. 研究结论为mPOF在温度传感领域的应用及各种可调光纤器件的制备提供了理论参考.