激光激发黏弹表面波有限元数值模拟

研究黏弹性材料中激光激发的Rayleigh波的传播特征. 考虑到黏弹性材料的黏性特征,在频域内建立黏弹性材料中激光激发Rayleigh波的有限元数值模型. 在验证有限元频域数值模型正确性的基础上,模拟脉冲激光作用在黏弹性材料上激发出Rayleigh波,进而讨论激光激发的黏弹Rayleigh波的传播特征,并比较黏弹性材料与弹性材料中激光激发的Rayleigh波差异,同时分析了材料的黏性劲度参量变化对Rayleigh波特征的影响. 关键词： 表面波 激光超声 有限元方法 黏弹性     

黏弹材料动力学参数的宽频测试

利用改进的波速法,对黏弹性材料连续宽频范围的动力学参数进行测试。利用有限元方法,对波速法的测试过程进行了仿真,验证了利用长短棒波速法测量宽频动力学参数的有效性。采用可控脉冲生成技术,在激振器上产生了用于宽频测试的短脉冲信号。黏弹性长短棒在上述宽带短脉冲激励下作纵向强迫振动,利用激光测振仪测量长短两棒自由端的纵向振动速度,从而可得两振动信号在连续宽频范围的幅值比和相位差,进而可根据波速法原理计算得到材料在连续宽频范围的储能模量和损耗因子。测试结果表明,该系统通过较少次数测试,可直接计算得到1~5 kHz连续宽频范围的动力学参数,测试结果与黏弹仪数据吻合。通过对幅值比和相位差进行线性最小二乘拟合,可以进一步拓展测量的频率范围。该方法准确可靠、简便快速,具有实际应用价值。     

黏弹材料动力学参数的宽频测试

利用改进的波速法,对黏弹性材料连续宽频范围的动力学参数进行测试。利用有限元方法,对波速法的测试过程进行了仿真,验证了利用长短棒波速法测量宽频动力学参数的有效性。采用可控脉冲生成技术,在激振器上产生了用于宽频测试的短脉冲信号。黏弹性长短棒在上述宽带短脉冲激励下作纵向强迫振动,利用激光测振仪测量长短两棒自由端的纵向振动速度,从而可得两振动信号在连续宽频范围的幅值比和相位差,进而可根据波速法原理计算得到材料在连续宽频范围的储能模量和损耗因子。测试结果表明,该系统通过较少次数测试,可直接计算得到1~5 kHz连续宽频范围的动力学参数,测试结果与黏弹仪数据吻合。通过对幅值比和相位差进行线性最小二乘拟合,可以进一步拓展测量的频率范围。该方法准确可靠、简便快速,具有实际应用价值。     

黏弹细棒中直达波提取及其动态力学参数的宽频测试

黏弹性材料的动态力学参数可通过对材料样品中纵波测试并结合波速法计算得到。由于测试样品尺寸的限制,除了高损耗材料在高频范围外,不管是在样品端头或侧面都难以获得准确的直达波信号,导致参数测量不准确。针对这个问题,对一端固定一端自由的黏弹性细棒中的纵波传播过程进行了分析,提出了黏弹性细棒直达波提取的方法,并实现了黏弹性材料动态力学参数的宽频测试。采用宽带脉冲作为测试激励信号,利用两台激光测振仪对细棒的固定端和自由端的振动速度进行测试,计算得到细棒自由端处的直达波。然后,利用提取出的直达波信号进行波速法计算,得到了黏弹性细棒宽频范围内的储能模量和损耗因子,宽频脉冲测试结果与黏弹仪测试结果相吻合,验证了该方法的有效性。基于直达波提取的宽频测试方法不仅效率高,而且可以使波速法应用于更低的频段。     

一种新型的激光自混频测振的方法

1引言在激光器的许多应用中，我们都尽可能的避开出射光的自反馈现象．尽管如此，还是被来自各种各样的光学元器件所组成的共振腔产生的反射光所困扰．自反馈不仅增加了噪音强度，而且同时大大改变了激光器输出的特性曲线，这种作用的性质已被Lenstra等人‘’‘描述过．然而我们正是利用了这种“不利”的特性构成了自混频激光测振系统．2基本原理激光自混频干涉的简单装置如图1所示．首先由激光垂直照射在待测物体上，利用从待测物体上反馈的光对激光器输出功率调制，然后从尾光接收调制信号，这种可被提取的携带外界馈信息的光反馈调制后…     

大动态范围窄脉冲激光参数测量探头

设计了一种双探测器结构的窄脉冲激光参数测量探头,解决了大动态测量范围和高速放大电路高增益时信噪比低的矛盾。通过对探头前端的光学特性进行分析,给出了均匀照射时探头前端的照度传递系数,并且对非均匀照射情况下可能的最大模型误差进行了分析,指出增大漫透射陶瓷片与探测器的距离和减小出光孔径有利于降低模型测量误差,减小双探测器偏离光轴的安装距离有利于提高两个探测器的测量一致性。测试表明,对于10 ns脉宽的1.064μm激光,所研制的探头在50 nJ~8 mJ范围内的测量误差低于±9%,输出激光脉冲波形幅度在0.7~1.4 V之间。利用推导出的照度传递系数公式还测出了漫透射陶瓷片的透过率为0.60,测量误差为±1.7%。     

大动态范围窄脉冲激光参数测量探头

设计了一种双探测器结构的窄脉冲激光参数测量探头,解决了大动态测量范围和高速放大电路高增益时信噪比低的矛盾。通过对探头前端的光学特性进行分析,给出了均匀照射时探头前端的照度传递系数,并且对非均匀照射情况下可能的最大模型误差进行了分析,指出增大漫透射陶瓷片与探测器的距离和减小出光孔径有利于降低模型测量误差,减小双探测器偏离光轴的安装距离有利于提高两个探测器的测量一致性。测试表明,对于10 ns脉宽的1.064 μm激光,所研制的探头在50 nJ～8 mJ范围内的测量误差低于±9%,输出激光脉冲波形幅度在0.7～1.4 V之间。利用推导出的照度传递系数公式还测出了漫透射陶瓷片的透过率为0.60,测量误差为±1.7%。     

弹光调制器动态参数测量与高效驱动匹配研究

为了提高弹光调制器（photoelastic modulator,PEM）的驱动效率,提出了一种基于PEM动态参数测量的阻抗匹配参数计算方法。理论分析了PEM调制幅度与驱动电压大小的关系,建立了PEM及其谐振匹配网络的等效电路模型,在此基础上推导了PEM谐振驱动电压与各个参数的关系,并设计了针对PEM动态测试的测量系统与验证方法。对谐振频率为44.822 kHz的PEM进行各项特性曲线的理论仿真与试验比较,验证了动态测量系统的可靠性,并通过数值仿真得到PEM最优匹配参数。在最优匹配参数附近选取不同匹配参数,测量实际谐振驱动电压,测量结果与仿真曲线相关性达0.996 4,变化趋势在10μF前基本吻合,且在理论最优匹配参数下,实际谐振驱动电压峰值达510 V,优于其他驱动匹配方法,接近理论最大值,相对误差小于1.16%。通过PEM等效动态参数推算的最优匹配参数,可使PEM达到最大驱动效率。     

一种基于激光辐照热效应的薄膜参数反演方法

本文研究并建立了一种基于激光辐照热效应的薄膜参数反演方法.首先给出激光辐照薄膜产生温升问题的热传导理论模型,并利用拉普拉斯变换得到了膜层和基底温度场的解析解;然后以膜层和基底的导热系数为反演参数,基于非线性共轭梯度算法给出反演基本原理及流程,并推导得到了反演过程中灵敏度系数的解析表达式;以aluminum,silver,copper和gold四种金属薄膜为例,通过与有限元法的计算结果对比验证了温度场解析解的正确性;最后结合四种金属薄膜进行了参数反演,通过考察分析不同随机噪声等条件下的参数反演结果,验证了本文方法在薄膜参数反演精度与反演效率等方面的有效性.反演结果显示:本文方法具有较高的反演精度和效率,在迭代截止误差为10~(-7)时只需用少于20次迭代就能收敛;在测量数据中加入的随机噪声越小,反演的迭代收敛次数就越少,即使是在迭代初值与反演结果相差较大时,用包含5%随机噪声的测量数据反演也能快速收敛.本文提出的薄膜参数反演方法不仅适用于反演导热系数,也可扩展用于反演膜层反射系数或吸收率等参数,具有一定的适用性.本文方法对于激光加工或激光损伤过程中的参数反演及优化具有一定的指导意义.     

黏弹性材料的动态力学特性及其对覆盖层吸声性能的影响

本文利用标准化动态力学测量手段获得了某种高分子聚合物的动态杨氏模量,并根据时温等效原理对动态杨氏模量与声学测量在频段上的差异加以分析和转换,得到了500—7500 Hz频率范围内该黏弹性材料杨氏模量随频率变化的特性.基于所测得动态杨氏模量,采用有限元方法分析了均匀黏弹材料的吸声性能,并将仿真结果与样品声管实验数据进行对比,验证了测试所得参数的准确性.进一步仿真分析了含有局域共振结构的声学覆盖层吸声性能,并讨论了黏弹性材料的动态特性对其吸声性能的影响,提出了改进水声覆盖层低频宽带吸声特性的建议.     

PVDF水听器加速度响应特性的有限元分析

本文用有限元法分析了压电高聚物PVDF圆柱形水听器的材料参数与几何尺寸对其加速度响应的影响。针对理论分析,进行了水听器加速度响应实验,实验结果与理论分析基本相符。     

Analysis of brake squeal noise using the finite element method: A parametric study

Brake squeal noise has been under investigation by automotive manufacturers for decades due to consistent customer complaints and high warranty costs. J.D. Power surveys consistently show brake noise as being one of the most critical vehicle quality measurements. Furthermore, the development of methods to predict noise occurrence during the design of a brake system has been the target of many researchers in recent years.This paper summarizes the application of complex eigenvalue analysis in a finite element model of a commercial brake system. The effect of the operational parameters (friction coefficient, braking pressure and brake temperature) and wear on the dynamic stability of the brake system is examined. After identifying unstable frequencies and the behavior of the brake system under different conditions, the performance of some control methods are tested. Changes in material properties and the application of brake noise insulators are presented and their effects discussed.The results show that the effect of brake temperature changes the coupling mechanisms between rotor and pad, which in some cases can be useful in order to reduce the instabilities and generated noise. Wear is an operational condition that has an strong effect on the system instability, since stiffness properties of brake pads are influenced by the changes on geometry and on the friction material, leading to high-frequency noise generation when the system is in the end of its lifetime. Application of brake insulators requires a detailed investigation of the system since, for some cases, an increase on the system damping does not balance changes on stiffness, leading the system to instability and noise.     

A method for the low frequency qualification of reverberation test rooms using a validated finite element model

Low frequency behavior in small rooms is always a critical issue, but the recent extension of several standards to frequencies as low as 50 Hz opened an interesting debate in the scientific community as to which is the best (and most reliable) method to perform measurements. The present paper discusses the low frequency qualification of a typical reverberant test room in order to perform sound power measurements (carried out according to ISO 3741), by taking advantage of a finite element model of the room. Experimental measurements were first carried out in a standard reverberant chamber to demonstrate that the model provides accurate results in the range below 100 Hz. Statistical analysis of the results from measurements and simulations confirmed that, despite some small inaccuracies, the predicted results are in very good agreement with those measured both in terms of spectra and spatial distribution of the sound pressure level. Finally, the different steps of the low-frequency qualification of a reverberant test room are discussed. A selection procedure of the most suitable microphone positions is proposed, based on the results of the simulation, and, finally, on site measurements were carried out to validate the procedure.     

Birefringence study of photonic crystal fibers by using the full-vectorial finite element method

Modal solutions of photonic crystal fibers with equal and unequal circular air holes in a hexagonal matrix are presented, by using a rigorous full-vectorial finite element-based approach. The effective indices, mode field profiles, spot-sizes, modal hybridness, modal birefringence and group velocity dispersion values have been determined and presented. The effects of the pitch-distance, hole diameter, structural asymmetry, air hole arrangement and the operating wavelength on the modal birefringence are also reported. It is shown that a significant value of birefringence can be achieved by using only circular air holes, which would be easy to fabricate, and by operating it close to its modal cutoff. PACS 42.81.Qb; 42.81.Gs; 42.25.Bs; 31.15.Pf     

A contribution to simulating a three-dimensional larynx model using the finite element method

A three-dimensional model is presented to simulate the larynx during vocalization. The finite element method is used to calculate the airflow velocity and pressure along the larynx as well as tissue displacement. It is assumed that the larynx tissue is transversely isotropic and divided into three tissues: cover, ligament, and body. A contact-impact algorithm is incorporated to deal with the physics of the collision between both true vocal folds. The results show that the simulated larynx can reproduce the vertical and horizontal phase difference in the tissue movements and that the false vocal folds affect the pressure distribution over the larynx surfaces. The effects of exciting the larynx with different pressure drops are also investigated.     

Viscoelastic analysis of multi-body systems using the finite element method

A study of the effect of viscoelastic material damping on the dynamic response of multibody systems, consisting of interconnected rigid, elastic and viscoelastic components, is presented. The motion of each elastic or viscoelastic body is identified by using three sets of modes: rigid body, reference and normal modes. Rigid body modes describe translation and large angular rotation of a body reference. Reference modes are the result of imposing the body-axis conditions. Normal modes define the deformation of the body relative to the body reference. Constraints between different components are formulated by using a set of non-linear algebraic equations that can be introduced to the dynamic formulation by using a Lagrange multiplier technique or can be utilized to eliminate dependent co-ordinates by partitioning the constraint Jacobian matrix. In developing the system equations of motion of the viscoelastic component, an assumption of a linear viscoelastic model is made. A Kelvin-Voigt model is employed, wherein the stress is assumed to be proportional to the strain and its time derivative. The formulation yields a constant damping matrix and the damping forces depend only on the local deformation; thus, no additional coupling between the reference and elastic co-ordinates appears in the formulation when considering the viscoelastic effects. It is demonstrated, by a numerical example, that the viscoelastic material damping can have a significant effect on the dynamic response of multibody systems.     

Characterization of silicon nanophotonic devices using the finite element method

The Poynting vector and the full-vectorial H and E-field profiles are considered for use in nanoscale silicon waveguides in this article. This paper reveals that the mode profile of a circular silicon nanowire is not circular and also has a strong axial field component. From the analysis, the characteristics of single mode operation and the vector field profiles of both circular and planar silicon nanowires are presented. The modal birefringence of rectangular silicon nanowires and power density in low-index region of a slot-type waveguide and designs of a compact polarization rotator are also presented in this work.     

二维超声应变重建的有限元仿真与组织实验

近几年来,利用超声散射回波测量生物组织弹性系数的超声应变成象方法得到了很大的发展。本文建立了基于有限元分析的二维软组织超声散射模型,并通过这一模型获得了组织在压缩前后的仿真超声散射回波。对模型组织超声散射回波的应变重建实验结果表明,组织在压缩过程中垂直于压缩方向的形变会引入重建应变误差,采用中心部位的扫描线,减小压缩比能减小该误差,提高应变图象的质量。组织实验的结果表明,有限元散射模型能模拟实际组织应变成象的情况,为组织实验提供参考。     

A finite element model using multi-layered shell element in laser forming

Due to its enormously high flexibility laser forming has been gaining importance in recent years. This rapidness and flexibility demand very precise controlling strategies especially when simulating the process of large plates and challenging the limited computation power of the current workstation. A simple, robust and accurate modeling method of laser forming has been demonstrated to solve this problem. The simplified model is meshed by multi-layered shell element, simulated with a more real scanning method and fewer parameters. The intelligent meshing strategies have reduced the number of elements dramatically. Thus the simulation efficiency has been improved significantly. By comparing the simulation results under the simplified model with the results under the traditional model for laser forming, the applicability of proposed method has been proven. The method of these simplified models is also suitable to simulate complex finite element models, which take much time to simulate. It would throw some light on the thermal mechanically coupled-field simulation of large sheet.     

Numerical simulation of premixed combustion using an enriched finite element method

In this paper we present a novel discretization technique for the simulation of premixed combustion based on a locally enriched finite element method (FEM). Use is made of the G-function approach to premixed combustion in which the domain is divided into two parts, one part containing the burned and another containing the unburned gases. A level-set or G-function is used to define the flame interface separating burned from unburned gases. The eXtended finite element method (X-FEM) is employed, which allows for velocity and pressure fields that are discontinuous across the flame interface. Lagrange multipliers are used to enforce the correct essential interface conditions in the form of jump conditions across the embedded flame interface. A persisting problem with the use of Lagrange multipliers in X-FEM has been the discretization of the Lagrange multipliers. In this paper the distributed Lagrange multiplier technique is adopted. We will provide results from a spatial convergence analysis showing good convergence. However, a small modification of the interface is required to ensure a unique solution. Finally, results are presented from the application of the method to the problems of moving flame fronts, the Darrieus–Landau instability and a piloted Bunsen burner flame.