A numerical scheme for calculating stresses and strain rates in glaciers

The stress and strain-rate fields in a glacier or other type of rock with known rheological properties can be calculated by using a numerical solution technique. Calculations are based on force-balance equations and the constitutive relation for polar ice, and do not involve any mathematical approximations. The geometry is prescribed. Basal velocities also are prescribed, either by specifying their values or by relying on a sliding relation. All other quantities are determined numerically.     

An application of Karmarkar's interior-point linear programming algorithm for multi-reservoir operations optimization

Optimization of multi-reservoir systems operations is typically a very large scale optimization problem. The following are the three types of optimization problems solved using linear programming (LP): (i) deterministic optimization for multiple periods involving fine stage intervals, for example, from an hour to a week (ii) implicit stochastic optimization using multiple years of inflow data, and (iii) explicit stochastic optimization using probability distributions of inflow data. Until recently, the revised simplex method has been the most efficient solution method available for solving large scale LP problems. In this paper, we show that an implementation of the Karmarkar's interior-point LP algorithm with a newly developed stopping criterion solves optimization problems of large multi-reservoir operations more efficiently than the simplex method. For example, using a Micro VAX II minicomputer, a 40 year, monthly stage, two-reservoir system optimization problem is solved 7.8 times faster than the advanced simplex code in MINOS 5.0. The advantage of this method is expected to be greater as the size of the problem grows from two reservoirs to multiples of reservoirs. This paper presents the details of the implementation and testing and in addition, some other features of the Karmarkar's algorithm which makes it a valuable optimization tool are illuminated.     

Evaluation of rainfall-runoff performance using the stochastic integral equation method

The stochastic integral equation method (S.I.E.M.) is used to evaluate the relative performance of a set of both calibrated and uncalibrated rainfall-runoff models with respect to prediction errors. The S.I.E.M. is also used to estimate confidence (prediction) interval values of a runoff criterion variable, given a prescribed rainfall-runoff model, and a similarity measure used to condition the storms that are utilized for model calibration purposes.Because of the increasing attention given to the issue of uncertainty in rainfall-runoff modeling estimates, the S.I.E.M. provides a promising tool for the hydrologist to consider in both research and design.     

二维不均匀介质中点源P-SV波响应的有限差分近似算法

本文从一阶方程组形式的波动方程出发,发展了一种计算二维不均匀介质中点源P-SV波响应的近似方法。该方法通过引入线分布的应力作为震源,利用二维有限差分方法计算出线源响应,然后再经过波形校正和几何扩散校正得出相应的近似点源响应。通过把波形和振幅与精确解比较表明,该方法具有较好的精度。由于有限差分方法对于介质中速度和密度的分布没有特殊要求,另一方面,本文所给出的震源可以适用于位错点源、爆炸源或集中力源,因此上述方法十分适合于研究横向不均匀介质中的近场强地运动、爆炸振动或地震勘探等问题。     

非线性弹性模型的改进

本文在邓肯-张模型进行深入探讨的基础上,提出了切线模量实质物理意义是剪切模量的观点,并分别采用应力法和应变法推导了剪切模量和体积模量(或变形模量)公式,还进行了考虑Biot固结非线性弹性油罐地基有限元实例分析。对不同的变形模量计算方法进行了比较,提出了改进的非线性弹性模型。     

求解非线性有限元方程的弧长法及在工程稳定分析中的应用

本文讨论了求解非线性有限元方程弧长法的基本原理及具体实施,并应用于压杆及滑块的稳定性分析。由于弧长法具有确定并通过结构临界应力状态的临界荷载值,因此可以应用于工程稳定性分析。算例表明,弧长法确定的临界荷载值是令人满意的。     

遥感筛选评价区化异常的初步探索

近年来,作者发展了七种以上重要找矿指示元素的野外快速分析法,此乃其一。文中提出一种新的金样冷浸技术,引进了新的光导光度检测方法,样品经溴化钠—硫酸—过氧化氢冷浸,用经丙酮—MIBK处理的泡沫塑料富集金,经混合解脱剂解脱并直接取解脱液用改进的微珠析出法或用光导光度法测定。取10g样品,可测定0．0005- 50g/t范围内的金。方法简便、快速,可在普查分队驻地进行,每一工作日可测定80多个样品,由于废弃了王水分解样品和泡沫塑料的灰化,具有广阔的应用前景。     

均匀离散分布的平面应变过程和Fry法应变测量

本文运用 Apple Macitonsh 计算机对无应变的均匀离散点分布进行系列平面应变模拟,应变叠加模式分别采用具普遍意义的成岩压实+顺层缩短(LPS)+压溶作用和成岩压实+简单剪切+压溶作用。对各应变阶段的变形点分布进行相应的 Fry 法应变测量,并配合低变形砂岩样品的 Fry 法应变分析实例.证实 Fry 法应变测量方法为一非常有用的应变测量方法,其结果不仅能较好地揭示全岩总应变特征,而且能揭示出许多应变叠加的信息,Fry 法揭示的全岩有限应变椭球主面的方位也较为真实可靠。而 Fry 法运用于应变分布不均匀的劈理化岩石中时,能揭示不同变形域的应变特征,从而达到应变分解的目的。     

Application of the weighted residual method to diffraction by 2-D canyons of arbitrary shape: I. Incident SH waves

A solution for the two-dimensional scattering and diffraction of plane SH waves by canyons of arbitrary shape in an elastic half space is presented. The wave field for arbitrary geometry in this paper is computed numerically by the method of weighted residues (moment method). The wave displacement field computed by the present residual method for the case of a semi-circular canyon was shown to agree analytically and numerically with that computed by the exact closed form series solution. The same observations about ground amplifications, their dependence on frequencies and orientations of the incident waves, can be stated here for canyons of arbitrary shape as previously made for circular canyons.