高边坡结构可靠度的二次处理有限元分析

本文依据强度折减理论,利用MIDAS/GTS有限元软件,分析计算了高边坡结构的安全系数K,找到边坡滑裂带的位置。在此基础上,对有限元输入数据和输出结果进行二次处理,建立基本随机变量c,f与滑裂带中单元的最大(和最小)主应力σ1(和σ3)的拟合关系f1(和f3),将其代入高边坡结构的功能函数Z中,使Z由隐式形式变为显式。基于该显式表示的Z,利用Monte Carlo法计算滑裂带中所有失效单元的可靠指标β1,并将其单元面积A1作为权重系数,经过加权平均得到边坡结构的整体可靠指标β。上述方法使得结合有限元软件计算边坡结构的整体可靠度得以简化。经实例分析可知,本文提出的方法是合理可行的,可使边坡结构整体可靠性分析得以简化,也可为高边坡结构整体可靠性分析提供理论参考。     

金华山软岩铁路隧道施工过程围岩屈服接近度分析

隧道施工过程中围岩处于复杂应力状态下,隧道围岩屈服区演化特征的确定对于围岩稳定性分析和开挖支护方案优化具有重要的意义。采用屈服接近度指标衡量围岩破坏接近程度可以合理地描述复杂应力状态下围岩的应力危险性,对Mohr-Coulomb类岩体材料的屈服接近度函数进行了相应的推导,并在非线性有限元用户子程序上编程予以实现。介绍了赣州-龙岩铁路DKl33+095～DKl38+237段软弱围岩单线隧道正台阶步施工方案以及湿喷纤维混凝土支护方案。为了对该隧道施工过程中隧道围岩屈服区的演化特征进行合理评价,采用非线性有限元法对软弱围岩条件下的铁路隧道施工过程进行了数值模拟,分析了施工过程中隧道围岩屈服接近度分布特征,判定了隧道台阶步施工过程中隧道围岩的稳定性。分析结果表明:该隧道施工过程中围岩破坏区主要发生在下台阶步施工过程中;屈服接近度指标比传统的塑性区分布提供的信息更加丰富,有利于工程技术人员定量地评价隧道开挖支护方案。     

香港花岗岩风化分级化学指标体系与风化壳分带

火成岩的化学风化过程实质上就是碱、碱土金属组分水解淋失与脱硅富铁铝化的地球化学过程。不同风化程度的岩石处于整个化学风化过程的不同阶段,因而具有不同的化学成分和风化度指标。笔者在对香港九龙两个地方5个钻孔共计150件不同深度、不同风化程度的样品化学全分析结果与岩土风化程度宏观判别结果的相互对比、选择和确认的基础上,建立了香港地区花岗岩风化岩土风化程度工程分级的化学指标体系。这一研究成果的提出不仅有利于提高花岗岩风化岩土工程分类的科。学性、可靠性和实用性,而且可以为风化岩土工程特性与风化作用关系研究提供理论依据。     

煤层底板原位综合测试及突水危险性评估

在实施井下钻孔水压致裂地应力测量的同时,评价岩层阻水性能,确定导升带高度,做到一孔多用。测试区实际观测表明,泥岩、粉砂岩、中砂岩、灰岩阻水性能由大到小,导升高度则由小变大,平均导升高度约32m; 地应力测量表明,该测试区原始地应力状态为SH＞Sv＞Sh,最大水平主应力方向NE至NEE。突水临界指数评价表明,-550m水平采煤工作面掘进时突水的可能性较大,应当采取防范措施。     

Transient MHD stagnation flow of a non-Newtonian fluid due to impulsive motion from rest

The transient boundary layer flow and heat transfer of a viscous incompressible electrically conducting non-Newtonian power-law fluid in a stagnation region of a two-dimensional body in the presence of an applied magnetic field have been studied when the motion is induced impulsively from rest. The non-linear partial differential equations governing the flow and heat transfer have been solved by the homotopy analysis method and by an implicit finite-difference scheme. For some cases, analytical or approximate solutions have also been obtained. The special interest are the effects of the power-law index, magnetic parameter and the generalized Prandtl number on the surface shear stress and heat transfer rate. In all cases, there is a smooth transition from the transient state to steady state. The shear stress and heat transfer rate at the surface are found to be significantly influenced by the power-law index N except for large time and they show opposite behaviour for steady and unsteady flows. The magnetic field strongly affects the surface shear stress, but its effect on the surface heat transfer rate is comparatively weak except for large time. On the other hand, the generalized Prandtl number exerts strong influence on the surface heat transfer. The skin friction coefficient and the Nusselt number decrease rapidly in a small interval 0<t^*<1 and reach the steady-state values for t^*≥4.     

正交异性体某些断裂动力学问题

首先对具有任意自相似指数的正交异性弹性动力学问题推导了解的一般表示,给出了一般解法,然后用这一方法对若干具体问题进行求解,利用本文方法可以迅速将所论问题化为半平面上的Ｒｉｍａｎｎ－Ｈｉｌｂｅｒｔ问题,并可以相当简单地得到问题的闭合解。     

Pore Network Analysis of Resistivity Index for Water-Wet Porous Media

Pore network analysis is used to investigate the effects of microscopic parameters of the pore structure such as pore geometry, pore-size distribution, pore space topology and fractal roughness porosity on resistivity index curves of strongly water-wet porous media. The pore structure is represented by a three-dimensional network of lamellar capillary tubes with fractal roughness features along their pore-walls. Oil-water drainage (conventional porous plate method) is simulated with a bond percolation-and-fractal roughness model without trapping of wetting fluid. The resistivity index, saturation exponent and capillary pressure are expressed as approximate functions of the pore network parameters by adopting some simplifying assumptions and using effective medium approximation, universal scaling laws of percolation theory and fractal geometry. Some new phenomenological models of resistivity index curves of porous media are derived. Finally, the eventual changes of resistivity index caused by the permanent entrapment of wetting fluid in the pore network are also studied.Resistivity index and saturation exponent are decreasing functions of the degree of correlation between pore volume and pore size as well as the width of the pore size distribution, whereas they are independent on the mean pore size. At low water saturations, the saturation exponent decreases or increases for pore systems of low or high fractal roughness porosity respectively, and obtains finite values only when the wetting fluid is not trapped in the pore network. The dependence of saturation exponent on water saturation weakens for strong correlation between pore volume and pore size, high network connectivity, medium pore-wall roughness porosity and medium width of the pore size distribution. The resistivity index can be described succesfully by generalized 3-parameter power functions of water saturation where the parameter values are related closely with the geometrical, topological and fractal properties of the pore structure.     

多体系统指标2运动方程HHT方法违约校正1）

采用Cartesian绝对坐标建模方法,完整约束多体系统运动方程是指标3的微分--代数方程（differentialalgebraic equations,DAEs）,数值求解指标3的DAEs属于高指标问题,通过对位置约束方程求导,可使运动方程的指标降为2.位置约束方程求导得到的是速度约束方程.直接求解指标3的运动方程,速度约束方程得不到满足,而且高指标DAEs的数值求解存在一些问题.论文首先采用HHT（Hilber--Hughes--Taylor）直接积分方法求解降指标得到的指标2运动方程,此时速度约束方程参与离散计算,从机器精度上讲速度约束自然得到满足,而位置约束方程没有参与计算,存在“违约”.针对违约问题,采用基于Moore--Penrose广义逆理论的违约校正方法,消除位置约束方程的违约.指标2运动方程HHT方法违约校正,将HHT方法和违约校正方法很好地结合,在数值求解指标2运动方程的过程中,位置约束方程和速度约束方程都不存在违约问题,而且新方法没有引入新的未知数向量,离散得到的非线性方程组的方程数量与原指标2运动方程的方程数量相同,求解规模没有扩大.新方法的实用和有效性通过算例的数值实验得到验证,数值实验也说明新方法保持了HHT方法本身具有的数值阻尼可以控制和二阶精度的特性.最后从非线性方程组的求解规模和计算速度上与其他方法进行了比较分析,说明新方法的优势所在.     

Multiplicity Results for Periodic Solutions to Second-Order Difference Equations

By using the Z _p geometrical index theory, some sufficient conditions on the multiplicity results of periodic solutions to the second-order difference equations

The Acceleration Severity Index in the impact of a vehicle against permanent road equipment support structures

The acceleration Severity Index (ASI), described in European Standard EN12767 (The passive safety of support structures for road equipment. Requirements, classification and the test method) is regarded as the most important indicator of impact on the occupants. The requirements for experiments are described, having in mind that the results depend on many factors. One of them is the selection of a vehicle to be used in the crash test. To perform numerical vehicle crash simulation, the finite-element models of permanent road equipment support structures were developed using the LS Dyna software available. To examine the response of the vehicle upon the impact, the acceleration severity index curves were calculated and visualised in Matlab.