Solution of the hyperbolic mild‐slope equation using the finite volume method

A finite volume solver for the 2D depth‐integrated harmonic hyperbolic formulation of the mild‐slope equation for wave propagation is presented and discussed. The solver is implemented on unstructured triangular meshes and the solution methodology is based upon a Godunov‐type second‐order finite volume scheme, whereby the numerical fluxes are computed using Roe's flux function. The eigensystem of the mild‐slope equations is derived and used for the construction of Roe's matrix. A formulation that updates the unknown variables in time implicitly is presented, which produces a more accurate and reliable scheme than hitherto available. Boundary conditions for different types of boundaries are also derived. The agreement of the computed results with analytical results for a range of wave propagation/transformation problems is very good, and the model is found to be virtually paraxiality‐free. Copyright © 2003 John Wiley & Sons, Ltd.     

Objective acoustic voice-quality parameters in the computer phonetogram

A phonetogram is a plot of the dynamic range of the voice as a function of fundamental frequency. Traditionally, the phonetogram only records the sound-pressure level (SPL) of the threshold of phonation and the upper limit in SPL the voice can reach with sustained phonation as a function of the fundamental frequency (F₀). In this study, the phonetogram is extended by including acoustic voice-quality parameters. Three additional parameters are tested: jitter, shimmer, and crest factor. For each individual voice, the variation in the three parameters is evaluated over the entire phonetogram area. By averaging individual phonetograms, standard or norm phonetograms are obtained revealing general patterns in voice-quality variation. These patterns reflect the complex relations between F₀, SPL, and the acoustic voice-quality parameters just mentioned. Jitter and shimmer distributions over the phonetogram show that large variations in perturbation values can be expected when production conditions are altered. Highest perturbation values are found for the low F₀ and low SPL phonations. For all voices, a gradual decrease of the crest factor is found with increasing F₀, reflecting the degrading of spectral complexity with F₀. The crest-factor parameter can also be used to mark those areas in the phonetogram where the fundamental dominates the waveform and where flow control is the main SPL regulating mechanism in voice production. The strong quality variations within the phonetogram stress the need for accurate control of F₀ and SPL in objective voice-quality measurement.     

用Talbot效应及云纹干涉法测试镜面反射体曲率和斜率的灵敏度分析

对用Talbot效应及云纹干涉法测试镜面反射体曲率和斜率的灵敏度进行了分析,同时给出了由于光栅稀密对测试的影响和相应的解释。这对从事这项工作或工程上利用此技术测试物体曲率和斜率提供了一个可行的方案。     

旋转锥镜相机研测动态曲率和斜率

采用旋转锥镜相机，进行激光散斑剪切照相，研测缓慢变形物体的动态曲率和斜率。拍摄一张散斑图即可获得曲率和斜率动态变化的全过程。     

Bahadur efficiency and robustness of studentized score tests

We derive the exact Bahadur slopes of studentized score tests for a simple null hypothesis in a one-parameter family of distributions. The Student's t-test is included as a special case for which a recent result of Rukhin (1993, Sankhy Ser. A, 55, 159–163) was improved upon. It is shown that locally optimal Bahadur efficiency for one-sample location models with a known or estimated scale parameter is attained within the class of studentized score tests. The studentized test has an asymptotic null distribution free of the scale parameter, and the optimality of likelihood scores does not depend on the existence of a moment generating function. We also consider the influence function and breakdown point of such tests as part of our robustness investigation. The influence of any studentized score test is bounded from above, indicating certain degree of robustness of validity, but a bounded score function is needed to cap the influence from below and to ensure a high power breakdown point. We find that the standard Huber-type score tests are not only locally minimax in Bahadur efficiency, but also very competitive in global efficiency at a variety of location models.     

MODELLING OF WAVE PROPAGATION IN THE NEARSHORE REGION USING THE MILD SLOPE EQUATION WITH GMRES-BASED ITERATIVE SOLVERS

The mild slope equation in its linear and non-linear forms is used for the modelling of nearshore wave propagation. The finite difference method is used to descretize the governing elliptic equations and the resulting system of equations is solved using GMRES-based iterative method. The original GMRES solution technique of Saad and Schultz is not directly applicable to the present case owing to the complex coefficient matrix. The simpler GMRES algorithm of Walker and Zhou is used as the core solver, making the upper Hessenberg factorization unneccessary when solving the least squares problem. Several preconditioning-based acceleration strategies are tested and the results show that the GMRES-based iteration scheme performs very well and leads to monotonic convergence for all the test-cases considered.     

三维非连续变形分析方法碰撞分析及其应用

将具有完备动力学理论的非连续变形分析(DDA)方法应用于块体碰撞研究。基于三维DDA(3D DDA)方法,按时步输出块体碰撞过程速度变化和接触嵌入量,进而得到块体碰撞恢复系数、冲量、冲击力。以此为参考指标,采用斜抛、面-面对心等碰撞模型,验证3D DDA方法模拟块体碰撞的有效性,并将3D DDA方法应用于多米诺骨牌倾倒、滚石边坡成灾及防护等算例分析,探讨了多米诺骨牌倾倒机制、滚石启动及运动行为、滚石灾害防护方案。结果表明:多米诺骨牌间距越大,同一块体被碰撞时间越迟,其最终稳定时间也越迟,与下一块体碰撞的动能越大;滚石运动呈侧向平动及转动三维运动特征,每一次碰撞,均引起动能、轨迹或状态的显著变化;滚石拦挡设施弹簧刚度越大,越先达到最大冲击力,最大冲击力随弹簧刚度的增加而减小;可结合树木阻挡效应,耗散滚石动能,降低滚石飞跃高度,使滚石灾害减轻或控制在防护范围以内。     

基于持续同调的边坡土颗粒应力链失稳分析

采用颗粒离散元方法和持续同调理论,研究了内排土场堆叠至不同高度时的边坡稳定性。为便于研究,现采用一水平金属板向下施加压力,代替不同厚度土层的重力荷载,对边坡在竖向荷载作用下的失稳破坏过程进行了颗粒离散元模拟。研究了二维边坡土颗粒速度总矢量、边坡失稳破坏时滑移开裂面的角度以及边坡坡顶y方向的平均速度等宏观响应过程,并构建了自然堆积下边坡堆积体颗粒的法向力链无向网络模型。最后,用持续同调方法对边坡坡顶颗粒接触力链网络的拓扑特征进行分析,获得条码图,建立了岩体结构持续同调特征与失稳演化的关系。本文为研究边坡失稳拓扑识别提供了一种新方法,从而可以有效预测边坡的失稳破坏。     

Double Feature Extraction Method of Ship-Radiated Noise Signal Based on Slope Entropy and Permutation Entropy

In order to accurately identify various types of ships and develop coastal defenses, a single feature extraction method based on slope entropy (SlEn) and a double feature extraction method based on SlEn combined with permutation entropy (SlEn&PE) are proposed. Firstly, SlEn is used for the feature extraction of ship-radiated noise signal (SNS) compared with permutation entropy (PE), dispersion entropy (DE), fluctuation dispersion entropy (FDE), and reverse dispersion entropy (RDE), so that the effectiveness of SlEn is verified, and SlEn has the highest recognition rate calculated by the k-Nearest Neighbor (KNN) algorithm. Secondly, SlEn is combined with PE, DE, FDE, and RDE, respectively, to extract the feature of SNS for a higher recognition rate, and SlEn&PE has the highest recognition rate after the calculation of the KNN algorithm. Lastly, the recognition rates of SlEn and SlEn&PE are compared, and the recognition rates of SlEn&PE are higher than SlEn by 4.22%. Therefore, the double feature extraction method proposed in this paper is more effective in the application of ship type recognition.     

时程稳定性系数确定的边坡逐孔起爆孔间微差降振时间

为了获得边坡逐孔爆破最佳降振微差时间,以某个实际边坡逐孔微差爆破施工现场为原型,先利用ANSYS建立二维静态模型,借助有限元折减法确定自然状态下的潜在滑动面和静态安全系数;基于已确定的二维潜在滑动面重新建立同尺寸同性质的三维逐孔微差爆破动态模型,利用LS-DYAN进行动力分析,整个过程分别设置同排3个炮孔0、17、25、42和65 ms等5种不同孔间微差起爆方式;同时,对该施工现场进行排、孔间(25 ms,17 ms)、(25 ms,25 ms)、(25 ms,42 ms)、(25 ms,65 ms)等4种微差时间控制的等比例相似小炮测振实验。提取模拟结果中3个炮孔同时起爆时滑面单元的应力数值代入极限平衡法计算公式,绘制了冲击载荷作用下边坡稳定性系数曲线,通过对曲线的理论分析发现,最佳降振微差时间约为48 ms;而三维数值模拟和测振实验结果均显示,孔间微差时间取42 ms时降振效果较佳。这说明,边坡稳定性系数曲线给出的微差时间与模拟和实验结果较为接近,可为今后边坡逐孔微差爆破降振研究提供参考。