Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

Exposure dating and reinterpretation of coarse debris accumulations (‘rock glaciers’) in the Cairngorm Mountains,Scotland

Relict rock glaciers have considerable potential for contributing to palaeoclimatic reconstruction, but this potential is often undermined by lack of dating control and problems of interpretation. Here we reinvestigate and date four proposed ‘rock glaciers’ in the Cairngorm Mountains and show that the morphology of only one of these appears consistent with that of a true rock glacier produced by creep of underlying ice or ice‐rich sediment. All four features comprise rockslide or rock avalanche runout debris, and the possibility that all four represent unmodified runout accumulations cannot be discounted. Surface exposure dating of the four debris accumulations using cosmogenic ¹⁰Be produced uncertainty‐weighted mean ages of 15.4 ± 0.8 ka, 16.2 ± 1.0 ka, 12.1 ± 0.6 ka and 12.7 ± 0.8 ka. All four ages imply emplacement under cold stadial conditions, two prior to the Windermere Interstade of ca. 14.5–12.9 cal. ka BP and two during the Loch Lomond Stade of ca. 12.9–11.5 cal. ka BP. The above ages indicate that paraglacial rock‐slope failure on granite rockwalls occurred within a few millennia after deglaciation. The mean exposure ages obtained for runout debris at two sites – Strath Nethy (16.2 ± 1.0 ka) and Lairig Ghru (15.4 ± 0.8 ka) – are consistent with basal radiocarbon ages from Loch Etteridge, 22 km to the southwest (mean = 15.6 ± 0.3 cal. ka BP) and imply widespread deglaciation of the Cairngorms and adjacent valleys before 15 ka and possibly 16 ka. Copyright © 2008 John Wiley & Sons, Ltd.     

Uncoupling of coupled flows in soil—a finite element method

Coupled flow of water, chemicals, heat and electrical potential in soil are of significance in a variety of circumstances. The problem is characterized by the coupling between different flows, i.e. a flow of one type driven by gradients of other types, and by the dual nature of certain flows, i.e. combined convection and conduction. Effective numerical solutions to the problem are challenged due to the coupling and the dual nature. In this paper, we first present a general expression that can be used to represent various types of coupled flows in soil. A finite element method is then proposed to solve the generalized coupled flows of convection-conduction pattern. The unknown vector is first decomposed into two parts, a convective part forming a hyperbolic system and a conductive part forming a parabolic system. At each time step, the hyperbolic system is solved analytically to give an initial solution. To solve the multi-dimensional hyperbolic system, we assume that a common eigenspace exists for the coefficient matrices, so that the system can be uncoupled by transforming the unknown vector to the common eigenspace. The uncoupled system is solved by the method of characteristics. Using the solution of the hyperbolic system as the initial condition, we then solve the parabolic system by a Galerkin finite element method for space discretization and a finite difference scheme for time stepping. The proposed technique can be used for solving multi-dimensional, transient, coupled or simultaneous flows of convection-conduction type. Application to a flow example shows that the technique indeed exhibits optimality in convergence and in stability.     

海上矿石运输最优船队的组建

针对海上矿石运输，介绍了一种完成定量任务并能使投资费用最少为目标的最优船队组建方法，同时对易变参数进行了灵敏度分析，使决策更加科学化     

A new type of Shipboard Meteorological Satellite Receiving-processing System

INTRODUCTIONTheSMSRPSisstalledontheshipnavigatingovertheseaandreceivessatellitecloudmapsatanytimetoprovidereliablereal-timedataofmeteorologyandocean.Itisanimportantequipmentforsafeguardofshipnavigation.Chinahasthousandsofoceanicships.Butalmostallships,withonlyveryfewexception,arenotstalledbySMSRPS.Thecausesareasfollows:1.Thetechniqueisverycomplicated.2.Thedevelopmentcostisveryhigh.3.Ifweintroduceforeignequipment,thecostistooexpensive.4.Foreignequipmentanditssoftwarearecompletelyclo…     

About the specification of erosion flux for soft stratified cohesive sediments

The present paper deals with the specification of bed erosion flux that accounts for the effects of sediment-induced stratification in the water column. Owing to difficulties in measuring the bed shear stress _b and the erosive shear strength _s, we suggest a series of methods that combine laboratory and numerical experiments. A simplified turbulent transport model that includes these effects helps to quantify _b and _s. Focusing on soft stratified beds, the present study considers erosion rate formulas of the form =_f exp {[T_b-T_s]} where is a model constant (=1 for Gularte's (1978) formula and =1/2 for Parchure's (1984) formula). First, the bed erosive strength profile _s(Z) is adjusted by forcing the turbulent transport model with measured erosion rates. Second, three procedures are suggested to determine the erosion rate formula coefficients _f and : a global procedure and two different layer-by-layer procedures. Each procedure is applied to an erosion experiment conducted in a rotating annular flume by Villaret and Paulic (1986). The use of the layer-by-layer procedure based on a least squares fitting technique provides a closer fit than the global procedure. The present study points out the complementarity of experimental and numerical approaches and also suggests possible improvements in laboratory test procedures.     

第二松花江干流区地下水信息管理数据库系统开发与研制

介绍了第二松花江干流区地下水信息管理数据库系统的开发环境及总体设计。在分析与地下水有关的数据信息的基础上,提出了较为合理的地下水信息管理数据库系统的结构。将该系统按其功能分为5个子系统基础子系统、数据信息管理子系统、区域信息管理子系统、动态绘图子系统、水资源评价子系统,并对5个子系统的功能进行了详细阐述。     

中国兴凯湖北岸平原晚全新世花粉记录及泥炭沼泽形成

通过兴凯湖北岸平原泥炭剖面高分辩率花粉分析研究，对晚全新世花粉划分4个组合带。XKH-4组合时期(1857～1746aB．P．)为沼泽发育前期，这一时期花粉浓度小，陆生草本植物占优势，气候干冷。XKH-3组合时期(1746～1287aB．P．)为沼泽发育早期，这一时期花粉浓度较大，且水生植物花粉含量为剖面最高，喜温落叶阔叶植物大发展时期，气候温和湿润。XKH-2组合时期(1287～602aB．P．)为沼泽发育中期，这一时期花粉浓度最小，但以陆生草本植物为主，木本植物为辅，水生植物急剧减少，针叶植物出现两次高峰，气候向冷干方向发展。XKH-1组合时期(602aB．P．至今)为沼泽发育盛期，这一时期花粉浓度最大，陆生草本植物大发展时期，气候波动较大。     

右江盆地富宁地区石炮组化石及顶部凝灰岩年代学特征及意义

通过野外地质调查, 笔者在右江盆地富宁地区中三叠统板纳组与下伏下三叠统石炮组界面之下约100 m采获早三叠世双壳类和菊石等化石, 在T1/T2界面之间新发现一套厚约5 m凝灰岩层, 对其进行LA-ICP-MS U-Pb同位素定年, 自下而上分别获得(249.0±1.9) Ma、(247.9±3.4) Ma、(245.5±1.6) Ma同位素年龄, 在误差范围内一致, 首次获得石炮组顶界年龄为下三叠统奥伦尼克阶。本次新发现的凝灰岩层及化石可作为早—中三叠世等时地层标志进行区域对比, 丰富了右江盆地早三叠世火山岩分布范围, 认为右江盆地滇东南地区T1/T2的划分标志以该套凝灰岩作为标志较为合理。结合前人资料, 这套凝灰岩、生物化石为华南陆块T1/T2界线的划分提供了重要年代学及古生物证据, 认为安尼阶/奥伦尼克阶界线年龄应定为 245.5 Ma, 可能更为合理。