弹性波传播模拟的Chebyshev谱元法

通过在每一个单元中采用谱展开近似,Chebyshev谱元法兼具了有限元处理边界及复杂结构的灵活性和谱方法的快速收敛特性,为弹性波传播的数值模拟提供了一种有效工具。从加权余量原理出发,详细阐述了Chebyshev谱元法用于弹性波传播模拟的基本理论及相应数学公式.给出了使用Chebyshev正交多项式展开得到的,存在等参变换时有关单元质量矩阵和单元刚度矩阵的精确积分公式。同时应用逐元技术极大地减少了内存和计算需求.最后,两个数值算例被用于验证这种谱元方法的高精度和强适应性。     

3维Mortar谱元法模拟Kelvin-Helmhtz不稳定性混合层

采用Mortar谱元法和多处理器并行计算技术模拟了Kelvin-Helmhtz界面不稳定性湍流的混合发展过程,通过对混合层动量厚度、能谱和总动能的计算,评估了Kelvin-Helmhtz混合层的演化机理。计算结果表明：3维Mortar谱元法具有高计算精度和光滑区域的指数收敛特性,可以有效模拟混合层流动的湍流混合和演化,能够捕捉到涡的合并现象和大涡到小涡的级联过程;初期的混合层层流运动发展成具有连续谱结构的湍流运动过程,实现了Kelvin-Helmhtz界面不稳定性混合层流动从2维发展到3维的转捩特征,总湍流统计动能的变化反映了粘性耗散过程的作用。通过对Kelvin-Helmhtz 3维界面不稳定性混合层流动和3维层流向湍流转捩过程的数值模拟,程序的有效性得到了验算,表明谱元法应用于湍流混合模拟是可行的。     

3维Mortar谱元法模拟Kelvin-Helmhtz不稳定性混合层

采用Mortar谱元法和多处理器并行计算技术模拟了Kelvin-Helmhtz界面不稳定性湍流的混合发展过程,通过对混合层动量厚度、能谱和总动能的计算,评估了Kelvin-Helmhtz混合层的演化机理。计算结果表明:3维Mortar谱元法具有高计算精度和光滑区域的指数收敛特性,可以有效模拟混合层流动的湍流混合和演化,能够捕捉到涡的合并现象和大涡到小涡的级联过程;初期的混合层层流运动发展成具有连续谱结构的湍流运动过程,实现了Kelvin-Helmhtz界面不稳定性混合层流动从2维发展到3维的转捩特征,总湍流统计动能的变化反映了粘性耗散过程的作用。通过对Kelvin-Helmhtz 3维界面不稳定性混合层流动和3维层流向湍流转捩过程的数值模拟,程序的有效性得到了验算,表明谱元法应用于湍流混合模拟是可行的。     

Influence of spectral dispersion parameters on the beam smoothing performance of a diffractive optical element

In this paper, a theoretical model to analyze the beam smoothing performance, realized by the smoothing by spectral dispersion (SSD) and the diffractive optical element (DOE), is set up with the use of the spatial frequency spectra of the far-field intensity distribution. The influences of the parameters of the SSD, for example, the frequency of the electrical field, the integration time and the angular dispersion of the diffraction grating, on the beam smoothing performance of the DOE are numerically analyzed. Furthermore, the phase distortion is considered, and the simulated results show the phase distortion should be changed during the integration time in order to obtain a good performance of beam smoothing.     

光谱色散匀滑技术对衍射光学器件束匀滑性能的影响

模拟分析了光谱色散匀滑技术对衍射光学器件束匀滑性能的影响。光谱色散匀滑技术的使用提高了衍射光学器件的抗波前畸变能力,降低了对加工精度的高要求。光谱色散匀滑技术与衍射光学器件的联用有望满足惯性约束聚变对束匀滑性能的高要求。     

光谱色散匀滑技术对衍射光学器件束匀滑性能的影响

 模拟分析了光谱色散匀滑技术对衍射光学器件束匀滑性能的影响。光谱色散匀滑技术的使用提高了衍射光学器件的抗波前畸变能力，降低了对加工精度的高要求。光谱色散匀滑技术与衍射光学器件的联用有望满足惯性约束聚变对束匀滑性能的高要求。     

Waves and energy in random elastic guided media through the stochastic wave finite element method

Energy propagation in random viscoelastic media is considered in this Letter. The forced response of uncertain waveguide subject to time harmonic loading is treated. This energy model is based on a spectral approach called the “Stochastic Wave Finite Element” (SWFE) method which is detailed in this Letter. Assuming that the random properties are spatially homogeneous in the media, the SWFE is a hybridization of the deterministic wave finite element and a parametric probabilistic approach. The proposed model is applicable in a wide frequency band with reduced time consumption. Numerical examples show the effectiveness of the proposed approach to predict the statistics of kinematic and quadratic variables of guided wave propagation. The results are compared to Monte Carlo simulations.     

多网格谱元法及其在高性能计算中的应用*

针对传统谱元法在每个单元内只能存在单一均匀介质,应用在复杂非均匀介质的波传播模拟中可能造成极大计算规模的问题,发展了多网格谱元法。该方法在谱元法单元内引入独立的辅助网格,用于精细描述单元内的介质和外力分布变化,在较稀疏的主网格上进行波场的求解。基于声波和弹性波方程推导了多网格谱元法公式,并对几种典型模型进行了波场的数值模拟。与传统谱元法的对比结果表明,此算法在复杂非均匀介质的弹性波传播模拟中可以利用较少的网格点数达到不低于传统算法的精度。此外,实现了并行化的多网格谱元法,获得了较好的并行效率。     

Simulation on laser-induced surface acoustic wave on isotropic cylinders by finite element method

Based on the thermoelastic theory, a finite element model is developed to simulate the process of laser inducing ultrasonic field in isotropic cylinders, which can take the temperature dependence of thermal parameters into account. Using the finite element model, we have simulated the ultrasonic fields induced by a pulse laser line source impacting on the generatrix of aluminum cylinders with different diameters. And the intact waveforms of surface acoustic wave (SAW including cylindrical Rayleigh and Whispering gallery (WG) modes) are presented, which are in very good agreement with the calculated and experimental waveforms in other literatures. Furthermore, the dispersion properties of cylindrical Rayleigh waves are analyzed by the method of phase spectral analysis, and the results show that with the increasing frequency, the phase velocity of cylindrical Rayleigh wave rapidly increases to the maximum value, and then gradually decreases to that of plane Rayleigh wave. With the diameter of cylinder decreasing, the maximum value of phase velocity and the corresponding frequency increase.     

石墨烯力常数和弹性波波速随温度变化规律研究

在哈里森键联轨道法框架下,考虑到原子的短程相互作用和原子的非简谐振动,建立石墨烯弹性模型并求出原子振动的简谐系数和非简谐系数。在此基础上,确定了石墨烯的力常数和弹性波波速随温度和形变势参量的变化关系式,探讨了原子非简谐振动对它们的影响。结果表明：（1）石墨烯的力常数和弹性波波速均随温度升高而非线性减小,但有不同的变化特征：在C11、C12、C44这三个力常数中,以C12随温度的变化最大,以C44的变化最小;纵波波速随温度的变化要大于横波;（2）石墨烯的三个力常数随形变势参量的变化有不同的特征：C11随η和γ的增大而增大、而C12则随η的增大而减小,随γ的增大而增大;C44随η的增大而增大,而几乎不随γ而变;（3）若不考虑原子振动的非简谐项,则石墨烯的力常数和弹性波波速均为常数。考虑到非简谐项后,不仅它们均随温度的升高而减小,而且,非简谐情况的力常数Cij和弹性波波速vi与简谐近似的值Cij0、vi0的差|Cij－Cij0|、|vi－vi0|均随温度的升高而增大,即温度愈高,非简谐效应愈显著。     

石墨烯力常数和弹性波波速随温度变化规律研究

在哈里森键联轨道法框架下,考虑到原子的短程相互作用和原子的非简谐振动,建立石墨烯弹性模型并求出原子振动的简谐系数和非简谐系数。在此基础上,确定了石墨烯的力常数和弹性波波速随温度和形变势参量的变化关系式,探讨了原子非简谐振动对它们的影响。结果表明：（1）石墨烯的力常数和弹性波波速均随温度升高而非线性减小,但有不同的变化特征：在C11、C12、C44这三个力常数中,以C12随温度的变化最大,以C44的变化最小;纵波波速随温度的变化要大于横波;（2）石墨烯的三个力常数随形变势参量的变化有不同的特征：C11随η和γ的增大而增大、而C12则随η的增大而减小,随γ的增大而增大;C44随η的增大而增大,而几乎不随γ而变;（3）若不考虑原子振动的非简谐项,则石墨烯的力常数和弹性波波速均为常数。考虑到非简谐项后,不仅它们均随温度的升高而减小,而且,非简谐情况的力常数Cij和弹性波波速vi与简谐近似的值Cij0、vi0的差|Cij－Cij0|、|vi－vi0|均随温度的升高而增大,即温度愈高,非简谐效应愈显著。     

Modified null field method for elastic wave scattering from a partially debonded fiber-SH-waves

Debonded fibers influence the macro-mechanical behavior of fiber-reinforced composite materials. Debonded fibers contribute to the initiation and growth of cracks at the fiber/matrix interface. To examine such problems, the scattering of elastic SH-waves (problem of anti-plane strain) from debonded fibers is studied with a numerical method which can handle the mixed boundary conditions on the partially bonded fiber. A modification of the null field of T-matrix method is developed for this purpose. The modification is achieved by the introduction of a mathematical surface. The simultaneous solution of the integral representations of the field scattered by the mathematical surface and the actual fiber surface give rise to sufficient equations that permit solution of the debonded fiber problem. The scattering cross-section and the far field amplitude are calculated as a function of frequency, fiber properties and debonding area. This will be of interest in structural applications where such cross-sections can be used to compute the dynamical effective properties of damaged composites.     

石墨烯力常数和弹性波波速随温度变化规律研究

在哈里森键联轨道法框架下,考虑到原子的短程相互作用和原子的非简谐振动,建立石墨烯弹性模型并求出原子振动的简谐系数和非简谐系数.在此基础上,确定了石墨烯的力常数随温度和形变势参量的变化关系式,以及弹性波波速随温度和形变势参量的变化关系式,并探讨了原子非简谐振动对它们的影响.结果表明:(1)石墨烯的力常数和弹性波波速均随温度升高而非线性减小,但有不同的变化特征:在C_(11)、C_(12)、C_(44)这三个力常数中,以C_(12)随温度的变化最大,以C_(44)的变化最小;纵波波速随温度的变化要大于横波;(2)石墨烯的三个力常数随形变势参量的变化有不同的特征:C_(11)随η和γ的增大而增大、而C_(12)则随η的增大而减小,随γ的增大而增大;C_(44)随η的增大而增大,而几乎不随γ而变;(3)若不考虑原子振动的非简谐项,则石墨烯的力常数和弹性波波速均为常数.考虑到非简谐项后,不仅它们均随温度的升高而减小,而且,非简谐情况的力常数C_(ij)和弹性波波速v_i与简谐近似的值C_(ij)~0、v_i~0的差|C_(ij)-C_(ij)~0|、|v_i-v_i~0|均随温度的升高而增大,即温度愈高,非简谐效应愈显著.     

Pure elastic resonance scattering of an obliquely incident plane acoustic wave by a submerged infinite cylinder

I.IntroductionTheResonanceScattcringThcory(RST)asanimportantmethodtoana1yzeacousticscat-teringofsubmergede1asticbodieswascxtcnsivelyuscdforcy1indrica1andspherica1geometries,inc1udingsolidstructuresandho11owoncs.RSTismain1yappliedtostudytheresonancescat-teringspcctra(oracousticspcctra)ofanc1astictargctcxcitedby'theinciderttacousticwavesoastOidentifythetargct.Tocxtracttheresonancescattcringspcctra,anappropriatebackgroundsuchasrigidbackgroundorsoftonemust-besubtractcd.Expcrimenta1ly,theReson…     

井间地震波场数值模拟和弹性波逆时偏移

该文利用时域有限差分法数值模拟了井间地震波场,探讨了井间地震弹性波逆时偏移成像的若干问题.基于多层地层模型,根据数值模拟的人工合成阵列信号,对比了3种不同成像条件下井间地震逆时偏移的结果,并对不同倾角的狭缝模型进行逆时偏移成像.通过对不同成像条件的对比,发现震源归一化成像条件能够对检波器一端的地质结构成像进行增强.由于...     

对于具有大范围运动和非线性变形的柔性梁的有限元动力学建模

研究具有大范围运动和非线性变形的柔性梁的有限元动力学建模.采用有限元方法对梁结构进行离散,利用Lagrange方程建立系统的精确动力学方程.该方程不仅增加了新的表征纵向、横向、侧向弯曲变形,以及扭转变形的耦合项,同时包含了变形运动与大范围运动之间的相互耦合项.     

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

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

Picosecond supercontinuum generation seeded by a weak continuous wave

We experimentally investigate the spectral details of a picosecond supercontinuum pumped at 1064 nm and seeded by a weak continuous wave (～20000 times weaker than the pulse peak power) at several power levels in photonic crystal fibers. Seeding at different wavelengths leads to different spectral details and the effects to the bandwidth of supercontinuum are also distinct. Spectra can be widened when seeded by a continuous wave at 1070 nm and narrowed by ～100 nm when seeded at 1080 nm. The influence is enhanced by increasing the average seeded power.     

弹性介质的Lagrange动力学与地震波方程

本文将地球介质看作是弹性介质, 从弹性体的Navier方程出发, 建立均匀弹性介质和非均匀弹性介质的分析动力学方程——Lagrange方程, 利用弹性介质的Lagrange方程导出匀弹性介质和非均匀弹性介质的地震波方程, 为用Lagrange分析动力学研究地球介质中地震波传播规律和解决地震勘探中的有关问题提供基础. 关键词： 地震勘探 弹性介质 Lagrange 方程 地震波方程     

Simulating regular wave propagation over a submerged bar by boundary element method model and Boussinesq equation model

Numerical models based on the boundary element method and Boussinesq equation are used to simulate the wave transform over a submerged bar for regular waves.In the boundary-element-method model the linear element is used,and the integrals are computed by analytical formulas.The Boussinesq-equation model is the well-known FUNWAVE from the University of Delaware.We compare the numerical free surface displacements with the laboratory data on both gentle slope and steep slope,and find that both the models simulate the wave transform well.We further compute the agreement indexes between the numerical result and laboratory data,and the results support that the boundary-element-method model has a stable good performance,which is due to the fact that its governing equation has no restriction on nonlinearity and dispersion as compared with Boussinesq equation.