结合材料界面端的三维应力奇异性

本文利用特殊有限元方法，开发了一个用来求解结合材料界面端三维应力奇异性问题的数值分析程序。该方法只需对界面端的角度方向进行离散即可求得应力奇异性。结合材料的应力奇异性取决于两种材料的材料常数和界面端形状。选用三个材料参数作为变量，用来研究结合材料三维应力奇异性随材料常数的变化规律。文中计算了几种重要而且常见的情况，并以此为基础建立了数据库。同时，还分析了应力奇异性随界面端形状的变化规律，并得到了应力函数的分布图。     

直角结合异材界面端应力强度系数的经验公式

由不同材料结合而成的材料(简称异材或双材料)的力学性能及其可靠性评价是工程中亟待解决的问题。表征界面端奇异应力场大小的应力强度系数是结合异材强度评价的依据，本文针对工程中最常见的直角结合异材，通过对大量不同材料组合的异材的边界元数值分析，提出了界面端应力强度系数的近似计算公式，无量纲化后的应力强度系数的值只与异材Dundurs参数a，卢有关，该公式具有较高的精度，可以作为一般工程上的应力强度系数的计算以及异材结构设计的依据。     

分析磁流变流体屈服应力微观力学模型

剪切屈服应力是磁流变流体固化强度的重要指标之一,具有重要的工程意义。基于磁性物理学理论,从磁流变流体在磁场作用下铁磁性固体颗粒极化成链的微观结构出发,探讨磁流变流体中固体颗粒间的相互作用机理,分析研究颗粒间的相互作用力,建立了一种微观力学模型,可用于分析磁流变流体在外加磁场作用下屈服应力及其影响因素的作用效果,揭示磁流变效应的微结构机理,为磁流变流体的性能优化、工程开发及应用磁流变流体提供理论依据。     

Smoothed Particle Hydrodynamics Model for Diffusion through Porous Media

A smoothed particle hydrodynamics (SPH) model is presented for the study of diffusion in spatially periodic porous media. The method of SPH is formulated to solve the convection–diffusion equation for tracer diffusion under steady state and transient conditions. Solutions obtained using SPH are compared with other available solutions and the model is used to calculate diffusion coefficients of spatially periodic porous media for the steady state diffusion problem. Diffusion coefficients are then used to calculate nondimensional diffusivities of the media. The effects of media properties on the values of nondimensional diffusivity are also presented.     

Analytic Model for DBF Under Multiple Particle Retention Mechanisms

Discrepancies between classical model (CM) predictions and experimental data for deep bed filtration (DBF) have been reported by various authors. In order to understand these discrepancies, an analytic continuum model for DBF is proposed. In this model, a filter coefficient is attributed to each distinct retention mechanism (straining, diffusion, gravity interception, etc.). It was shown that these coefficients generally cannot be merged into an effective filter coefficient, as considered in the CM. Furthermore, the derived analytic solutions for the proposed model (PM) were applied for fitting experimental data, and a very good agreement between experimental data and PM predictions were obtained. Comparison of the obtained results with empirical correlations allowed identifying the dominant retention mechanisms. In addition, it was shown that the larger the ratio of particle to pore sizes, the more intensive the straining mechanism and the larger the discrepancies between experimental data and CM predictions. Finally, the CM and PM were compared via statistical analysis. The obtained $p$ values allow concluding that the PM should be preferred especially when straining plays an important role.     

Modified Particle Detachment Model for Colloidal Transport in Porous Media

Particle detachment from the rock during suspension transport in porous media was widely observed in laboratory corefloods and for flows in natural reservoirs. A new mathematical model for detachment of particles is based on mechanical equilibrium of a particle positioned on the internal cake or matrix surface in the pore space. The torque balance of drag, electrostatic, lifting and gravity forces, acting on the particle from the matrix and the moving fluid, is considered. The torque balance determines maximum retention concentration during the particle capture. The particle torque equilibrium is determined by the dimensionless ratio between the drag and normal forces acting on the particle. The maximum retention function of the dimensionless ratio (dislodging number) closes system of governing equations for colloid transport with particle release. One-dimensional problem of coreflooding by suspension accounting for limited particle retention, controlled by the torque sum, allows for exact solution under the assumptions of constant filtration coefficient and porosity. The explicit formulae permit the calculation of the model parameters (maximum retention concentration, filtration and formation damage coefficients) from the history of the pressure drop across the core during suspension injection. The values for maximum retention concentration, as obtained from two coreflood tests, have been matched with those calculated by the torque balance on the micro scale.     

天津解放桥恢复开启工程应力监测

为使有80年历史的钢结构桥[CD2]天津市海河上的解放桥恢复开启功能,需对其进行全面整修,并对原桥的活动跨作较大的改造.为使解放桥恢复开启工作安全顺利实现,掌握活动跨承载下的力学性能,采用了实验应力分析电测法对活动跨进行称重、开启及预压试验进行监测.监测结果表明:解放桥在施工中是安全的,在恢复开启功能的运营中也是安全的.监测结果为类似工程的安全修建和科学施工提供了依据.     

一种测试动态应力集中的灰色建模方法

将灰色系统方法理论中的GM（1,1）模型,应用于动态应力集中测试的数据处理,根据应力集中处附近的应变测试值,能方便地求出应力集中处的最大应变值。结果表明,该方法简单,预测精度较高,应用方便。     

复杂应力状态下岩土材料非线性本构模型研究

基于统一强度理论,将洛德参数引入统一强度理论,推导出材料的统一强度参数,进而分析中间主应力以及主剪应力系数对材料统一强度参数的影响。在此基础上,对邓肯张双曲线模型进行改进,使得该模型能够反映复杂应力状态下的应力应变关系。采用粘土的平面应变试验进行验证,结果表明:在平面应变情况下,当b=0.4时,模型能够较好地反映粘土的应力应变关系,及其强度参数。分析了洛德参数以及中间主应力系数对模型的影响,进一步说明岩土材料存在主应力效应。该模型能够反映不同材料在复杂应力状态下的应力应变关系,有其更为广泛的适用性。     

模型复合材料弹塑性界面应力分析

由纤维增强弹塑性基体所产生的界面具有弹塑性力学行为。考虑到一般材料的塑性变形都遵循幂硬化规律,对模型复合材料的界面进行弹性和应变硬化状态下的变形规律及其应力分析。以纤维拔出试验为研究模型,将界面分成弹性区和塑性区。利用界面应力剪滞理论,分别建立弹性区和塑性区的界面力学基本方程。选择适当的位移函数满足基本方程及埋入纤维的边界条件,再按位移函数求出弹性区和塑性区的界面剪应力。推导出平均界面剪应力与纤维     

机匣模拟疲劳试验的相似条件研究

理论分析和试验结果表明，受内压的机匣与其相应的平板模型在开孔附近的等效应力分布非常相似；机匣模拟疲劳试验的应力相似条件下可表为，机匣与其模型在控制疲劳寿命的关键部位的等效应力水平和分布相同。     

Hypo-Elasticity Model Based upon the Logarithmic Stress Rate

Recently these authors have proved [46, 47] that a smooth spin tensor Ω^log can be found such that the stretching tensor D can be exactly written as an objective corotational rate of the Eulerian logarithmic strain measure ln V defined by this spin tensor, and furthermore that in all strain tensor measures only ln V enjoys this favourable property. This spin tensor is called the logarithmic spin and the objective corotational rate of an Eulerian tensor defined by it is called the logarithmic tensor-rate. In this paper, we propose and investigate a hypo-elasticity model based upon the objective corotational rate of the Kirchhoff stress defined by the spin Ω^log, i.e. the logarithmic stress rate. By virtue of the proposed model, we show that the simplest relationship between hypo-elasticity and elasticity can be established, and accordingly that Bernstein's integrability theorem relating hypo-elasticity to elasticity can be substantially simplified. In particular, we show that the simplest form of the proposed model, i.e. the hypo-elasticity model of grade zero, turns out to be integrable to deliver a linear isotropic relation between the Kirchhoff stress and the Eulerian logarithmic strain ln V, and moreover that this simplest model predicts the phenomenon of the known hypo-elastic yield at simple shear deformation. This revised version was published online in August 2006 with corrections to the Cover Date.     

气固两相流场的湍流颗粒浓度理论模型

本文进行了气固两相流动颗粒湍流扩散现象的理论分析，提出了颗粒湍流扩散系数和气流弥散效应二个颗粒湍流模化新概念，在此基础上建立了气固两相流场湍流颗粒浓度模型。理论模型包括离心力和其它外加力场作用下颗粒运动和浓度分布的计算方法。运用湍流颗粒浓度模型，对直管气固两相流动、受限射流气固两相流动和９０°弯管气固两相流动等三种流动做了数值模拟，计算获得颗粒速度、颗粒浓度等主要流动参数。讨论了湍流颗粒浓度模型的适用性。     

A Model for Spontaneous Imbibition as a Mechanism for Oil Recovery in Fractured Reservoirs

The flow of oil and water in naturally fractured reservoirs (NFR) can be highly complex and a simplified model is presented to illustrate some main features of this flow system. NFRs typically consist of low-permeable matrix rock containing a high-permeable fracture network. The effect of this network is that the advective flow bypasses the main portions of the reservoir where the oil is contained. Instead capillary forces and gravity forces are important for recovering the oil from these sections. We consider a linear fracture which is symmetrically surrounded by porous matrix. Advective flow occurs only along the fracture, while capillary driven flow occurs only along the axis of the matrix normal to the fracture. For a given set of relative permeability and capillary pressure curves, the behavior of the system is completely determined by the choice of two dimensionless parameters: (i) the ratio of time scales for advective flow in fracture to capillary flow in matrix $\alpha =\tau ^f/\tau ^m$ ; (ii) the ratio of pore volumes in matrix and fracture $\beta =V^m/V^f$ . A characteristic property of the flow in the coupled fracture–matrix medium is the linear recovery curve (before water breakthrough) which has been referred to as the “filling fracture” regime Rangel-German and Kovscek (J Pet Sci Eng 36:45–60, 2002), followed by a nonlinear period, referred to as the “instantly filled” regime, where the rate is approximately linear with the square root of time. We derive an analytical solution for the limiting case where the time scale $\tau ^{m}$ of the matrix imbibition becomes small relative to the time scale $\tau ^{f}$ of the fracture flow (i.e., $\alpha \rightarrow \infty $ ), and verify by numerical experiments that the model will converge to this limit as $\alpha $ becomes large. The model provides insight into the role played by parameters like saturation functions, injection rate, volume of fractures versus volume of matrix, different viscosity relations, and strength of capillary forces versus injection rate. Especially, a scaling number $\omega $ is suggested that seems to incorporate variations in these parameters. An interesting observation is that at $\omega =1$ there is little to gain in efficiency by reducing the injection rate. The model can be used as a tool for interpretation of laboratory experiments involving fracture–matrix flow as well as a tool for testing different transfer functions that have been suggested to use in reservoir simulators.     

Characterization and Control of Residual Stress and Curvature in Anodically Bonded Devices and Substrates with Etched Features

While anodic bonding is commonly used in a variety of microelectromechanical systems (MEMS) applications, devices and substrates that incorporate this processing technique are often subjected to significant residual stress and curvature that create post-processing and reliability issues. Here, using an anisothermal anodic bonding procedure, residual stresses and the resulting wafer curvature in these structures are controlled by varying the initial bond temperatures of the silicon and Pyrex wafers independently. Residual stresses are quantified by measuring bulk wafer curvature and, locally, stress concentrations are measured using infrared photoelasticity accompanied by 3-D thermomechanical finite element analysis. Based on the good agreement between numerical predictions and experimental results, this process can be used to determine the bulk post-bond wafer curvature and to reduce the likelihood of structural failure at these sites, by changing the residual stresses from tensile in nature, which may drive initiation and growth of cracks, to compressive, which can suppress such failures.     

不同应力分量下广义开尔文模型粘性系数探讨

对不同应力分量下的广义开尔文模型应力应变关系进行了研究,推导了在不同应力分量下的广义开尔文模型的粘性应变增量计算式;通过对这些粘性应变增量计算式的比较分析,得到结论:对于线性粘弹性模型,当应力张量引起粘性变形的规律与应力偏量和球应力分别引起粘性变形的规律相同时,它们的系数满足关系式Ek/ηk=Gsk/ηsk=Kmk/ηmk;否则,这个关系式不成立.现有文献采用应力张量表示的粘性变形有限元计算式隐含假定了球应力与应力偏量产生的粘性变形规律相同.对于复杂的工程材料而言,这种假定并不总是合适的.这在工程问题粘性分析时值得注意.     

一种非统一黏塑性本构的有限元实现

在与率相关黏塑性本构模型的背应力分量中加入静态回复项通常用来更精确地描述材料应力松弛效应．本文基于包含静态回复项的非统一黏塑性本构模型,提出了对应的隐式应力积分算法．该模型可简化为两个关于偏背应力和等效应力张量的方程组,通过牛顿迭代法进行求解．在应力积分算法中采用了Return-Mapping求解策略．此外还推导出了材料的一致切线刚度矩阵来判断有限元分析收敛性．在此基础上开发出了适用于ABAQUS工程软件的UMAT用户材料子程序．文中选取了高温下镍合金材料在伴随应力松弛的循环加载路径下的试验结果来验证算法的有效性．从比较结果中可以发现：用户子程序即使在大时间增量下也具有良好的收敛性,静态回复项的引入使得模型对应力松弛效应的模拟更加准确．     

Accounting for Buoyancy Effects in the Explicit Algebraic Stress Model: Homogeneous Turbulent Shear Flows

Most explicit algebraic stress models are formulated for turbulent shear flows without accounting for external body force effects, such as the buoyant force. These models yield fairly good predictions of the turbulence field generated by mean shear. As for thermal turbulence generated by the buoyant force, the models fail to give satisfactory results. The reason is that the models do not explicitly account for buoyancy effects, which interact with the mean shear to enhance or suppress turbulent mixing. Since applicable, coupled differential equations have been developed describing these thermal turbulent fields, it is possible to develop corresponding explicit algebraic stress models using tensor representation theory. While the procedure to be followed has been employed previously, unique challenges arise in extending the procedure for developing the algebraic representations to turbulent buoyant flows. In this paper the development of an explicit algebraic stress model (EASM) is confined to the homogeneous buoyant shear flow case to illustrate the methodology needed to develop the proper polynomial representations. The derivation is based on the implicit formulation of the Reynolds stress anisotropy at buoyant equilibrium. A five-term representation is found to be necessary to account properly for the effect of the thermal field. Thus derived, external buoyancy effects are represented in the scalar coefficients of the basis tensors, and structural buoyancy effects are accounted for in additional terms in the stress anisotropy tensor. These terms will not vanish even in the absence of mean shear. The performance of the new EASM, together with a two-equation (2-Eq) model, the non-buoyant EASM of Gatski and Speziale (1993) and a full second-order model, is assessed against direct numerical simulations of homogeneous, buoyant shear flows at two different Richardson numbers representing weak and strong buoyancy effects. The calculations show that this five-term representation yields better results than the 2-Eq model and the EASM of Gatski and Speziale where buoyancy effects are not explicitly accounted for. Received 5 March 2001 and accepted 15 January 2002     

直线振动筛相似模拟动应力研究

对27m2大型直线振动筛的相似试验模型进行了动应力测试,得到了试验模型筛筛体梁、板关键点处的最大动应力应变值。根据量纲分析法和方程分析法推导的工作状态参数相似关系,通过对原型筛与缩小的相似试验模型筛的试验结果对比分析,验证了相似模型与原型试验数据的一致性及可靠性,证实了利用缩小的相似试验模型进行动态测试及分析的有效性。本文结果为提高直线振动筛的使用寿命及结构的优化设计提供了理论依据,也为相似试验模型的发展研究提供了参考。