GNSS船姿测量方法及对多波束测深精度的影响分析

GNSS船姿测量以其观测误差不随时间累积的特点得到了广泛研究和应用,本文基于三天线GNSS船姿测量方式,构建了波束脚印误差与姿态误差间的关系模型,设计仿真实验分析了基线长度对姿态误差的影响,以及不同水深环境下姿态误差与GNSS定位误差的关系。为突破传统RTK在测量距离上的限制,本文采用PPP、PPK、MBD (动态参考站差分)三种方法进行GNSS船姿计算,并通过海上实验与高精度惯性导航系统进行对比分析,结果表明使用MBD测姿结果要优于PPK和PPP模式,得到的航偏角、横摇角、纵摇角标准差均在0.1°左右,可满足通常情况下多波束测深对姿态精度的要求。     

不同排列方式和间隙比下双圆柱流致振动特性研究

采用基于ALE的流固耦合分析方法,将圆柱体振动模型简化为平面双自由度质量-弹簧-阻尼系统,进行了计算参数无关性分析。对不同θ角和间隙比T下的双圆柱流致振动进行数值模拟,并与单圆柱振动计算结果进行比较,得到了不同θ角和间隙比T下中心圆柱的升阻力系数、流场分布和运动轨迹变化情况。通过分析发现小间隙比下由于两圆柱间的相互干扰,改变了中心圆柱的尾涡结构,导致了升力和阻力系数随不同θ角而发生变化;小间隙比下由于两圆柱间的相互干扰影响,均使中心圆柱的轨迹在不同θ角下发生了较大的变化,顺流向的振动频率与单圆柱振动时相比出现了多倍频,导致其运动轨迹发生了变化。     

非球形轴对称粒子光散射的散射矩阵元素展开系数的数值计算

扼要地介绍了非球形粒子光散射的散射矩阵元素按广义球函数展开的基本理论,以及球形粒子散射矩阵元素间的关系。在此基础上较详细地阐述了计算展开系数的方法以及在程序研制中所遇到的一些问题。最后,对非球形的Ｃｈｅｂｙｓｈｅｖ粒子光散射的散射矩阵元素的展开系数进行了计算,结果表明数值技术和计算的正确性。     

Effect of connecting rain barrels on the storage size reduction

In this article, the possibility of sharing rain barrels and the potential benefit of reducing storage size through physical and non‐physical connections of rain barrels in a community are investigated. Using the concepts of homogeneous/heterogeneous users in rainwater harvesting systems (RWHS), two simple cases of a community composed of four prospective users are examined. The first is performed with the users who have the same mean and variance in water demands (homogeneous users), and the second is with the users with different means and variances (heterogeneous users). To take account for the rainfall characteristics in different places, historical records from six cities in the USA are used for storage–reliability–yield analysis. The result indicates that required total storage can be reduced by connecting multiple rain barrels. In addition, a significant difference is found between homogeneous and heterogeneous user groups. Homogeneous users do not achieve a substantial benefit from connecting their rain barrels; these users may even be disadvantaged by sharing. In contrast, heterogeneous users receive benefit by reducing the total required storage. Most benefit is expected between users with maximum difference in mean water demands. The reduction in storage size was as considerable as 37% in this study. The quantity of storage reduction depends on locations and target reliabilities. Knowledge of the benefits and limitations of rain barrel connections can improve RWHS performance through ability to customize a network plan for individual users. Copyright © 2011 John Wiley & Sons, Ltd.     

Preconditioned IDR(s) iterative solver for non‐symmetric linear system associated with FEM analysis of shallow foundation

Non‐associated flow rule is essential when the popular Mohr–Coulomb model is used to model nonlinear behavior of soil. The global tangent stiffness matrix in nonlinear finite element analysis becomes non‐symmetric when this non‐associated flow rule is applied. Efficient solution of this large‐scale non‐symmetric linear system is of practical importance. The standard Krylov solver for a non‐symmetric solver is Bi‐CGSTAB. The Induced Dimension Reduction [IDR(s)] solver was proposed in the scientific computing literature relatively recently. Numerical studies of a drained strip footing problem on homogenous soil layer show that IDR(s = 6) is more efficient than Bi‐CGSTAB when the preconditioner is the incomplete factorization with zero fill‐in of global stiffness matrix K_ep (ILU(0)‐K_ep). Iteration time is reduced by 40% by using IDR(s = 6) with ILU(0)‐K_ep. To further reduce computational cost, the global stiffness matrix K_ep is divided into two parts. The first part is the linear elastic stiffness matrix K_e, which is formed only once at the beginning of solution step. The second part is a low‐rank matrix Δ, which is re‐formed at each Newton–Raphson iteration. Numerical studies show that IDR(s = 6) with this ILU(0)‐K_e preconditioner is more time effective than IDR(s = 6) with ILU(0)‐K_ep when the percentage of yielded Gauss points in the mesh is less than 15%. The total computation time is reduced by 60% when all the recommended optimizing methods are used. Copyright © 2013 John Wiley & Sons, Ltd.     

Analytical layer‐element method for non‐axisymmetric consolidation of multilayered soils

A numerically efficient and stable method is developed to analyze Biot's consolidation of multilayered soils subjected to non‐axisymmetric loading in arbitrary depth. By the application of a Laplace–Hankel transform and a Fourier expansion, the governing equations are solved analytically. Then, the analytical layer‐element (i.e. a symmetric stiffness matrix) describing the relationship between generalized displacements and stresses of a layer is exactly derived in the transformed domain. Considering the continuity conditions between adjacent layers, the global stiffness matrix of multilayered soils is obtained by assembling the inter‐related layer‐elements. Once the solution in the Laplace–Hankel transformed domain that satisfies the boundary conditions has been obtained, the actual solution can be derived by the inversion of the Laplace–Hankel transform. Finally, numerical examples are presented to verify the theory and to study the influence of the layered soil properties and time history on the consolidation behavior. Copyright © 2011 John Wiley & Sons, Ltd.     

Rayleigh 阻尼作用下斜拉桥扭转振动

在侧向风力或侧向水平地源力作用下,斜拉桥容易发生扭转振动。本文采用非线性Rayleigh阴尼,它能反映在低速振动时有激发而加速振动,至一定高速后又衰减的完整振动过程。对非线性的统辖方程及定解条件,先使之无量纲化,再按非线性项的系数这个小参数展开未知函数,得到线性化的各级近似方程。把待求函数展成梁轴坐标的Fourier级数,其系数为待定的时间函数,应用Fourier级数的正交性,得到求解未知时间函数的二阶常微分方程,完全类似单自由度质点在线性粘滞性阻尼作用下的统辖方程,其解为已知。     

Global spectra of energy and enstrophy and their fluxes during July 1979

Transient and stationary spectra of kinetic energy (KE), available potential energy (APE) and enstrophy (EN), and their spectral fluxes as a function of the two-dimensional wavenumbern were computed for July 1979. Triangular truncation at zonal wavenumber 42 was used for computation. The slopes of various spectra in the wavenumber range 14≤n≤25 were obtained by fitting a straight line in log-log scale by the least square method. The transientKE, APE andEN spectra in the lower (upper) troposphere had slopes −2·21 (−2·30), −2·65 (−2·64) and −0·36 (−0·46), respectively. The effect of stationary and divergent motion on the slope values was investigated. The possible correlation between the slope and percentage of transient component in the combined energy and enstrophy was examined to identify the transient motion of the atmosphere with the two-dimensional homogeneous isotropic turbulence. The vertically averaged slope of kinetic energy and enstrophy in the lower (upper) troposphere was close to the value at 700 (200) hPa level. The spectral fluxes of kinetic energy and enstrophy in the wavenumber range 14≤n≤25 satisfied, to a very rough approximation, the criteria of inertial subrange. The stationary fluxes were small. The estimated stationary-transient component of flux was larger, comparable and less than the corresponding transient flux of APE, KE and EN. Representative levels for computation of energy and enstrophy spectra and their fluxes in the lower and upper troposphere were identified.     

地震反应分析中的主时段

在计算反应谱及土层反应时,发现反应谱的形态和大小只受时程的主要部分控制,由此给出主时段的定义和确定主时段的方法。主时段的应用不仅可以大大缩短计算时间,而且能够成为持时的一种有物理意义的度量。     

区域化变量的空间聚类分析

传统的聚类分析方法利用数据对平均值的偏差，研究变量之间的相互关系。当它应用于区域化变量的分类时，因缺乏对区域化变量的空间结构的考虑，对合理解释区域化变量之间的相互关系不利。空间聚类分析方法比较好地解决了这个问题。本文介绍了该方法的基本原理及其算法，并说明了它在非条件模拟及在某矿区中的应用情况