Accelerated initial stiffness schemes for elastoplasticity

Iterative methods for the solution of non‐linear finite element equations are generally based on variants of the Newton–Raphson method. When they are stable, full Newton–Raphson schemes usually converge rapidly but may be expensive for some types of problems (for example, when the tangent stiffness matrix is unsymmetric). Initial stiffness schemes, on the other hand, are extremely robust but may require large numbers of iterations for cases where the plastic zone is extensive. In most geomechanics applications it is generally preferable to use a tangent stiffness scheme, but there are situations in which initial stiffness schemes are very useful. These situations include problems where a nonassociated flow rule is used or where the zone of plastic yielding is highly localized. This paper surveys the performance of several single‐parameter techniques for accelerating the convergence of the initial stiffness scheme. Some simple but effective modifications to these procedures are also proposed. In particular, a modified version of Thomas' acceleration scheme is developed which has a good rate of convergence. Previously published results on the performance of various acceleration algorithms for initial stiffness iteration are rare and have been restricted to relatively simple yield criteria and simple problems. In this study, detailed numerical results are presented for the expansion of a thick cylinder, the collapse of a rigid strip footing, and the failure of a vertical cut. These analyses use the Mohr–Coulomb and Tresca yield criteria which are popular in soil mechanics. Copyright © 2000 John Wiley & Sons, Ltd.     

煤矿采空区控制网起算点的变形分析与选择

大部分控制网建立时采用的起算点为50年代至80年代施测的三角点,起算点的可靠性和兼容性,直接关系到控制网建立的成败。通过大量生产实践数据详细介绍了起算点可靠性和兼容性分析、比较的过程和方法,可在今后生产中借鉴和应用。     

COMPARATIVE STUDY ON THE SEISMIC RETROFIT OF A MID-RISE STEEL BUILDING: STEEL BRACING VS ENERGY DISSIPATION

The paper presents a comparative study of an existing retrofit for a mid-rise steel building using additional stiff steel braced-frames against an alternate retrofit using ADAS (Added Damping and Stiffness) passive energy dissipation devices. The subject building, located near Alameda Park in downtown Mexico City, is a ten-storey office building that was built in the 1950s. The structure was damaged during the 1985 Michoacán Earthquake because of resonant response with the site. The building was later retrofitted using additional braced frames according to the seismic provisions of Mexico's 1987 Federal District Code. The retrofit scheme was planned to take the structure away from resonant responses and to inhibit structural damage. A proposed upgrade using ADAS energy dissipation devices was studied to compare energy dissipation against traditional stiffening using steel braces as retrofit options for mid-rise buildings in Mexico City's lake-bed zone. Different sets of analyses were carried out to compare both alternatives: (a) three-dimensional elastic analyses; (b) limit analyses and; (c) nonlinear dynamic analyses for postulated site ground motions for a M_s=8.1 earthquake. Initial costs of the retrofit schemes were also studied. The comparative studies suggest that a retrofit using ADAS devices would have a better dynamic performance than the one using steel braces. However, the steel bracing retrofit provides more strength and its initial cost of retrofit is less than that of the ADAS retrofit. © 1997 John Wiley & Sons, Ltd.     

带双腹板顶底角钢连接初始刚度研究

初始连接刚度在节点刚度分类和钢结构地震分析中占有很重要的地位,因此研究连接的初始刚度很有现实意义。主要研究带双腹板顶底角钢连接初始刚度的确定方法以及影响其性能的主要因素。     

流域水资源初始产权界定初探一以黑河流域中游为例

随着水资源短缺程度的加剧,可交易水权制度逐渐成为当今世界水资源管理制度发展的新方向.水资源初始产权是水资源管理的基础和前提.在干旱区内陆河流域,产权初始配置要在生态环境可持续的前提下,坚持公平、效益相结合的分配目标和配置原则.通过建立初始水权配置指标体系,对黑河流域进行实地问卷调查,在此基础上对用户调查结果进行专家判断修正,应用层次分析法研究了各分配指标的权重分布.以黑河流域张掖地区为例的界定结果表明,除上游的肃南县和山丹县外,其余各县市的实际用水量都超过了各自界定的初始水权,这种差别为水资源产权交易管理提供了现实依据,为开展流域上游的生态保护补偿提供了理论基础.     

Analysis on variation characteristics of geothermal response in Liaoning Province

Due to energy shortage and increasing environmental awareness, resources in shallow underground space have been rapidly exploited and utilized. So that studying variation characteristics of geothermal response in gneiss is necessary for effective and rational use of underground heat. Based on field test of thermal response in gneiss under hydrogeological survey project carried out in shallow geothermal energy development zone in Liaoning Province, this thesis analyzes mathematical statistics of geothermal response characteristics in main gneiss of Laoning Province. The initial formation temperature ranges from 10.80℃ to 15.80℃ according to field test. The statistical results show that in the condition of natural water content, the average thermal conductivity of Quaternary loose rocks comes as clay silty silty fine sand medium sand coarse sand gravelly sand. This order is consistent with thermal conductivity characteristics of gneiss obtained in the laboratory. Formation temperature recovery in different strata follows as granite medium sand clay. This order is opposite to the absolute value of temperature recovery curve slope of corresponding lithology, which shows that the stratum with higher temperature recovery rate has lower temperature recovery curve slope.     

微型抗滑桩单桩设计计算模型及算法研究

考虑微型桩与周围岩土体间摩擦力的作用,提出了一种新的微型抗滑桩单桩设计计算模型,并给出了具体算法。由微型桩加固滑坡体的变形特点,分析微型桩与岩土体之间的相互作用机制,将微型桩与周围岩土体的摩擦作用引入其受力分析中;根据微型桩上各部分受力特点的不同,将微型桩分成上部摩擦受拉段、中部滑坡推力作用段和下部锚固段3段进行分析,推导了微型桩总的变形控制方程及各分段的变形控制方程;采用初参数法对控制方程进行求解,得到了微型抗滑桩上的内力分布及变形规律。计算结果表明,在滑坡推力作用下,微型桩的变形主要发生在滑面附近及以上桩段,滑面附近桩段将产生较大的内力和弯曲变形,受拉段弯曲变形较小,近似水平移动;桩与岩土体间的摩擦力是微型桩与周围岩土体相互作用的重要组成部分,摩擦力的作用能显著减小微型桩的弯曲变形,有效控制滑坡体的位移。     

河流附近非水平初始面潜水渗流分析

考虑初始面非水平条件,建立河流边界控制下的半无限含水层渗流模型,提出非水平初始面及其相关参数（a、t_N）的简易计算方法,并探讨模型的适用条件。分析表明,该法适用于观测数据较少、地下水水位变幅相对于潜水含水层厚度较小、且观测孔距河岸较近条件下河流附近潜水含水层渗流计算。结合实例,采用查表法和配线法求解模型。计算结果表明:参数a和t_N的计算值不同,水位计算误差在0.02 m以内,较原公式精度有明显提高。新模型为初始面非水平条件下的潜水渗流计算提供一种简易方法。     

The relationship between nonconservative schemes and initial values of nonlinear evolution equations

For the nonconservative schemes of the nonlinear evolution equations, taking the one-dimensionalshallow water wave equation as an example, the necessary conditions of computational stability are given.Based on numerical tests, the relationship between the nonlinear computational stability and the construc-tion of difference schemes, as well as the form of initial values, is further discussed. It is proved throughboth theoretical analysis and numerical tests that if the construction of difference schemes is definite, thecomputational stability of nonconservative schemes is decided by the form of initial values.     

Field calibration of a formula for entrance mixing of river inflows to lakes: Lake Taupo,North Island,New Zealand

Field measurements were used to validate predictions for the initial dilution of negatively buoyant, cold‐water inflows to Lake Taupo, as part of a study to quantify mixing processes associated with the two largest inflows to the lake. The predictions were made using a formulation originally derived for positively buoyant, warm‐water inflows to cooling ponds. The formulation predicts the total dilution of an inflow during its inertia‐dominated phase between its entrance to the lake and the point where buoyancy forces are great enough to cause the inflow to plunge and form a submerged density current. In one of the measured inflows, the inflowing jet was free to entrain from both sides; in the other, entrainment was restricted on one side by attachment of the inflowing jet to the shoreline of a bay just upstream of the plunge point. In the former example, the unmodified coefficients from the cooling pond formulation provided an excellent prediction of initial dilution. In the latter example, entrainment was reduced and different coefficients were derived. In both examples the inflows remained attached to the lake bed throughout their course until their liftoff at depths of 45–55 m to form interflows. The difference between coefficients for the two inflows indicates that the coefficient values should be considered site‐specific. The formulation is not valid for inflows that separate from the bottom of the inflow channel before plunging. The entrance mixing formulation was incorporated in a more general model of lake stratification, DYRESM, which already includes a well‐documented routine for routing underflows down submerged channels on the bed of a lake. Application of the model to the inflows measured in Lake Taupo gave good results for two model outputs that were not involved in the calibration of the entrance mixing formulation, but that are affected by the result of the initial dilution calculation—the temperatures in the river plume after it has plunged, and the insertion depth.