动力方程的一种直接积分方法在负刚度条件下的收敛性和稳定性

研究了动力方程的一种直接积分法在负刚度条件下的收敛性和稳定性。在负刚度条件下,该方法是收敛的和无条件稳定的,它不仅适于分析有阻尼系统的动力反应,而且还适于分析无阻尼系统的动力反应。     

动力方程求解的显式积分格式及其稳定性与适用性

文献（１）给出了一种求解有阻尼体系动力方程的显式积分格式,文中以数值表格的形式给出了格式的稳定性条件,本文对该格式的稳定性问题作了进一步的分析,并给出了其计算稳定性条件的表示式。本文还着重讨论了基于这一显式积分格式的推导过程而派生出的另一形式的积分格式的稳定性,并指出了该派生格式的适用性问题。     

逐步积分法求解复阻尼结构运动方程的稳定性问题

平均常加速度法是无条件稳定的，但用它来解复阻尼振动方程时，现了不稳定现象，我们通过理论分析和数值计算结果的比较，探明了这种现象不是由于逐步积分算法本身有不稳定的问题，而是由于复阻尼振动方程解集本身含有不稳定子集的缘故。     

Implicit time integration for pseudodynamic tests: Convergence and energy dissipation

The convergence and energy-dissipation characteristics of an unconditionally stable implicit time integration scheme that has been adopted for pseudodynamic testing are examined in this paper. A convergence criterion is derived for general multiple-degree-of-freedom softening systems. Furthermore, it is shown that undesired loading and unloading cycles can be avoided in numerical iterations by scaling down the incremental corrections. Finally, it is proved that the total energy dissipation introduced by the residual convergence errors and proposed numerical correction is always positive for any softening structures.     

An error analysis of truncated starting conditions in step-by-step time integration: Consequences for structural dynamics

Most of the step-by-step time integration algorithms for structural dynamics require an initial acceleration vector to be specified, in addition to displacement and velocity vectors. A consistent initial acceleration vector may be calculated by solving the equations of motion at the initial time, while a truncated initial acceleration vector is obtained by setting the acceleration values to zero. Although the truncated starting procedure decreases computational effort, it is shown to affect accuracy adversely. For the structural dynamics algorithms considered herein, the rate at which the numerical solution converges to the exact solution is O(Δt) when the truncated starting procedure is used, compared to O(Δt²) when consistent initial acceleration values are used.     

求解振动方程的一种显式积分格式及其精度与稳定性

介绍了一种求解有阻尼体系振动方程的显式积分格式及其逐步求解的过程，并对其计算精度和稳定性进行了分析。该方法不但能同时求得体系的位移、速度和加速度反应，而且所得到的加速度反应的精度能满足工程需要。     

Numerical integration of a quasi-geostrophic atmospheric model with an asymmetric zonal current

Summary A two-level quasi-geostrophic model on a -plane, with an initial state characterized by an asymmetric zonal current was integrated numerically for a period of 14 days. Various field variables were computed and compared with the corresponding distributions obtained in two experiments with the same model and with initial states characterized by a constant zonal flow and by a flow varying sinusoidally with latitude.The research in this report has been sponsored wholly by the Air Force Cambridge Research Laboratories under Grant AF EOAR 63-108 through the European Office of Aerospace Research (OAR), US Air Force.     

Nodal domain integration model of two-dimensional advection-diffusion processes

The nodal domain integration method is applied to a two-dimensional advection-diffusion process in an anisotropic inhomogeneous medium. The domain is discretised into the union of irregular triangle finite elements with vertex-located nodal points and a linear trial function is used to approximate the governing flow equation's state variable in each element. Non-linear parameters are assumed quasi-constant for small durations in time in each element. The resulting numerical model represents the Galerkin and subdomain integration weighted residual methods and the integrated finite difference method as special cases. Both Dirichlet and Neumann boundary conditions are accommodated in a manner similar to the Galerkin finite element approach.     

海洋水体指数模型及其稳定性研究

基于遥感光学影像数据,针对海洋水体的特性,本文提出一种指数模型:在海洋水体平均反射率基础上,构建新的水体指数模型Ocean Water Index(OWI).主要以华南沿海地区、英国泽西岛为实验区,以Landsat 8 Surface Reflectance Tier 1影像、高分四号影像作为主要的数据源,与常用的8种...     

Developing land use scenario dynamics model by the integration of system dynamics model and cellular automata model

Models of land use change are useful tools for un-derstanding the analysis of the cause and conse-quences of land use changes, assessing the impacts of land use system on ecological system and supporting land use planning and policy[1,2]. Modeling land use scenario changes and its potential impacts on the structure and function of the ecosystem in the typical regions are regarded as one of the good ways to un-derstand the interactive mechanism between land use system and ecological system[3―10…     

非差非组合电离层加权PPP-RTK/INS紧组合模型

PPP-RTK是一种可实现快速模糊度固定的精密单点定位（PPP）技术,它的组成包括由多个接收机构成的网端用于生成各类精密改正产品,以及单站构成的用户端以实现快速高精度定位.当网端是中等尺度参考网,即参考站之间的站间距在百余公里,有必要引入电离层加权模型,其优势在于能够增强PPP-RTK函数模型.此外,从非差非组合观测值层面出发,其能够兼容单/双频用户.然而,单凭卫星导航定位技术很难在复杂环境下提供连续导航定位服务.本文从理论到应用,基于非差非组合电离层加权PPP-RTK模型,构建了单/双频PPP-RTK/INS紧组合模型.在距离最近参考站约70 km的区域进行了两次车载实验,并分别利用单/双频PPP-RTK以及PPP-RTK/INS紧组合模型的实时处理结果进行分析.结果表明,双频PPP-RTK/INS紧组合在半城市环境能从分米级定位提升至厘米级定位,重收敛时间从11 s提升至3.6 s.以水平偏差小于0.2 m评估定位可用率,双频PPP-RTK/INS紧组合模型在半城市环境和复杂城市环境中分别可实现100%和96.97%.考虑车道级定位,单频PPP-RTK/INS紧组合模型在半城市和复杂城市环境内分别可实现分米以及亚米级定位,其在水平偏差在0.5 m内的可用率分别为98.97%和76.1%.

     

Nodal domain integration model of one-dimensional advection—diffusion

The nodal domain integration method is applied to a one-dimensional advection—diffusion mathematical model without a source term. Comparison of the resulting numerical model to the well known Galerkin finite element, subdomain, and finite difference domain models indicates that a single numerical statement can be developed which includes the Galerkin finite element, subdomain, and finite difference models as special cases.     

A stability analysis of a conduit flow model for lava dome eruptions

Periodic variations in magma discharge rate and ground deformation have been commonly observed during lava dome eruptions. We performed a stability analysis of a conduit flow model by Barmin et al. [Barmin, A., Melnik, O., Sparks, R.S.J., 2002. Periodic behavior in lava dome eruptions. Earth and Planetary Science Letters 199 (1-2), 173–184], in which the periodic variations in magma flow rate and chamber pressure are reproduced as a result of the temporal and spatial changes of the magma viscosity controlled by the kinetics of crystallization. The model is reduced to a dynamical system where the time derivatives of the magma flow rate (dQ/dt) and the chamber pressure (dP/dt) are functions of Q and P evaluated at a shifted time t − t_?. Here, the time delay t_? represents the time for the viscosity of fluid particle to increase in a conduit. The dynamical system with time delay is approximated by a simple two-dimensional dynamical system of Q and P where t_? is given as a parameter. The results of our linear stability analyses for these dynamical systems indicate that the transition from steady to periodic flow depends on nonlinearities in the steady state relation between Q and P. The steady state relation shows a sigmoidal curve in Q − P phase plane; its slope has negative values at intermediate flow rates. The steady state solutions become unstable, and hence P and Q oscillate periodically, when the negative slope of the steady state relation ([dP/dQ]_S) exceeds a critical value; that is [dP/dQ]_S < − t_?γ/(2V_ch), where V_ch is the chamber volume and γ is an elastic constant which is related to the rigidity of chamber wall. We also found that the period and the pattern of oscillation of the conduit flow primarily depend on a quantity defined by LV_ch/r⁴, where L is the conduit length and r is the conduit radius.     

Improvements in turbulence model realizability for enhanced stability of ocean forecast and its importance for downstream components

Ocean Dynamics - The numerical stability of ocean circulation models is of high significance in operational forecasting. A substantial improvement in numerical stability of the 3D-ocean model HBM...     

Uncertainty and sensitivity in a bank stability model: implications for estimating phosphorus loading

Eutrophication of aquatic ecosystems is one of the most pressing water quality concerns in the United States and around the world. Bank erosion has been largely overlooked as a source of nutrient loading, despite field studies demonstrating that this source can account for the majority of the total phosphorus load in a watershed. Substantial effort has been made to develop mechanistic models to predict bank erosion and instability in stream systems; however, these models do not account for inherent natural variability in input values. To quantify the impacts of this omission, uncertainty and sensitivity analyses were performed on the Bank Stability and Toe Erosion Model (BSTEM), a mechanistic model developed by the US Department of Agriculture – Agricultural Research Service (USDA‐ARS) that simulates both mass wasting and fluvial erosion of streambanks. Generally, bank height, soil cohesion, and plant species were found to be most influential in determining stability of clay (cohesive) banks. In addition to these three inputs, groundwater elevation, stream stage, and bank angle were also identified as important in sand (non‐cohesive) banks. Slope and bank height are the dominant variables in fluvial erosion modeling, while erodibility and critical shear stress had low sensitivity indices; however, these indices do not reflect the importance of critical shear stress in determining the timing of erosion events. These results identify important variables that should be the focus of data collection efforts while also indicating which less influential variables may be set to assumed values. In addition, a probabilistic Monte‐Carlo modeling approach was applied to data from a watershed‐scale sediment and phosphorus loading study on the Missisquoi River, Vermont to quantify uncertainty associated with these published results. While our estimates aligned well with previous deterministic modeling results, the uncertainty associated with these predictions suggests that they should be considered order of magnitude estimates only. Copyright © 2016 John Wiley & Sons, Ltd.     

基于人工边坡稳定性的降雨型滑坡预警预报模型设计与应用

由削坡建房遗留的人工边坡存在大量滑坡隐患问题,在降雨引发土质边坡自身动力变化分析条件下,以稳定性评价建模为基础,提出降雨型滑坡动力学预警预报模型.文中以广东省梅州市花岗岩地区为例,使用GIS技术构建了1727个预警分析单元,并进行关键地质环境因子赋值及与气象站点数据关联;按坡高、坡度等参数,分别构建16个边坡失稳动力学...     

区域地震前兆数据集成共享模型

根据前兆监测台网的体系架构、运行机制、数据集成的概念与共享模式,提出了区域地震前兆数据集成共享模型,并在山东省地震局和中国地震台网中心进行了原型测试.该模型对各区域台网中心前兆集成数据库的建设、数据共享服务系统的开发,具有指导意义.     

Global stability of the ballooning mode in a cylindrical model

Ballooning disturbances in a finite-pressure plasma in a curvilinear magnetic field are described by the system of coupled equations for the Alfvén and slow magnetosonic modes. In contrast to most previous works that locally analyzed the stability of small-scale disturbances using the dispersion relationship, a global analysis outside a WKB approximation but within a simple cylindrical geometry, when magnetic field lines are circles with constant curvature, is performed in the present work. This model is relatively simple; nevertheless, it has the singularities necessary for the formation of the ballooning mode: field curvature and non-uniform thermal plasma pressure. If the disturbance finite radial extent is taken into account, the instability threshold increases as compared to a WKB approximation. The simplified model used in this work made it possible to consider the pattern of unstable disturbances at arbitrary values of the azimuthal wavenumber (k _y ). Azimuthally large-scale disturbances can also be unstable, although the increment increases with decreasing azimuthal scale and reaches saturation when the scales are of the order of the pressure nonuniformity dimension.     

基于稀疏反演的OBS数据多次波压制方法

自由表面多次波压制是海底地震仪（Ocean Bottom Seismometer,OBS）数据处理和成像中的难点,OBS数据多次波能量强,周期长,严重影响深层一次反射波的处理和成像.不同于常规拖缆观测系统,OBS数据站点一般相隔较远,仅仅利用检波点稀疏的波场信息难以压制OBS数据中的自由表面多次波.本文采用拖缆数据与OBS数据联合,利用稀疏反演估计（Estimation of Primaries and Multiples by Sparse Inversion,EPSI）方法,研究了OBS数据自由表面多次波压制理论,分析了OBS多次波产生的机理,详细推导了拖缆数据与OBS数据联合预测OBS多次波的EPSI方法基本原理.通过利用拖缆数据的信息,实现了OBS检波点稀疏数据多次波的压制问题.EPSI方法通过稀疏反演直接估计一次反射波,避免了SRME（Surface Related Multiple Elimination）方法中自适应相减对有效信号的损害,保真了一次反射有效信号,理论模拟OBS数据验证了方法的有效性.