峰后岩石非Darcy渗流的分岔行为研究

煤矿采动围岩大多处于峰后应力状态或破碎状态，其渗流一般不符合Darcy定律，为非Darcy渗流系统．峰后岩石非Darcy渗流系统的失稳和分岔是煤矿突水和煤与瓦斯突出动力灾害发生的根源．文中用谱截断方法建立了Ahmed-Sunada型非Darcy渗流系统的降阶动力学方程，再由变量代换得到以无量纲变量表示的平衡态附近的演化方程，分析了系统的分岔条件，给出了系统的各种吸引子图案，并结合采矿工程实际，用非线性数学的观点揭示了煤矿突水和煤与瓦斯突出的机理．研究表明：当非Darcy渗流系统渗流特性和边界压力的初始值满足一定条件时，系统由平衡转向不稳定，即存在跨临界Hopf分岔和切分岔，并且，系统的动力学响应不随渗透特性连续变化，即该系统存在突变性．     

一种计算变质量破碎岩体渗透性参量的方法

建立了一种基于压力梯度和渗流速度时间序列提取变质量破碎岩体渗透性参量(渗透率、非Darcy流β因子、加速度系数)的方法。首先,根据破碎岩体渗透性参量之间的幂指数关系和Forchheimer关系,建立了采样时刻渗透率的代数方程,并利用Newton切线法求得此代数方程的根;其次,对渗透性参量的参考值及幂指数进行优化。通过算例分析了该方法用于计算变质量破碎岩体渗透性参量的可行性及准确性。研究结果表明:随着细小颗粒的迁移流失,破碎岩体在恒定的压力梯度下,渗流速度和渗透率增加,非Darcy流β因子和加速度系数降低。     

考虑损伤的内变量黏弹-黏塑性本构方程

基于Rice 不可逆内变量热力学框架,在约束构型空间中讨论材料的蠕变损伤问题. 通过给定具体的余能密度函数和内变量演化方程推导出考虑损伤的内变量黏弹-黏塑性本构方程. 通过模型相似材料单轴蠕变加卸载试验对一维情况下的本构方程进行参数辨识和模型验证,本构方程能很好地描述黏弹性变形和各蠕变阶段.不同的蠕变阶段具有不同的能量耗散特点. 受应力扰动后,不考虑损伤的材料系统能自发趋于热力学平衡态或稳定态. 在考虑损伤的整个蠕变过程中,材料系统先趋于平衡态再背离平衡态发展. 能量耗散率可作为材料系统热力学状态偏离平衡态的测度;能量耗散率的时间导数可用于表征系统的演化趋势;两者的域内积分值可作为结构长期稳定性的评价指标.      

考虑损伤的内变量黏弹--黏塑性本构方程

基于Rice不可逆内变量热力学框架,在约束构型空间中讨论材料的蠕变损伤问题.通过给定具体的余能密度函数和内变量演化方程推导出考虑损伤的内变量黏弹--黏塑性本构方程.通过模型相似材料单轴蠕变加卸载试验对一维情况下的本构方程进行参数辨识和模型验证,本构方程能很好地描述黏弹性变形和各蠕变阶段.不同的蠕变阶段具有不同的能量耗散特点.受应力扰动后,不考虑损伤的材料系统能自发趋于热力学平衡态或稳定态.在考虑损伤的整个蠕变过程中,材料系统先趋于平衡态再背离平衡态发展.能量耗散率可作为材料系统热力学状态偏离平衡态的测度;能量耗散率的时间导数可用于表征系统的演化趋势;两者的域内积分值可作为结构长期稳定性的评价指标.     

水泥注浆渗透机理初探

本文在岩体等效渗流模型的基础上，结合注浆浆液对渗流过程的影响，推导出岩体注浆过程的控制方程，并给出确定注浆参数的一般方法．     

裂隙岩体非稳态渗流数值模型及其应用

为解决裂隙岩体非稳态渗流问题, 发展了一种新的数值模型. 对于单裂隙渗 流求解, 其控制方程是基于一定假设的简化Navier-Stokes方程, 数值方法采用有限差分法 和流体体积法. 在裂隙网络中, 交界处渗流可以由专门的控制方程求解. 计算结果表明, 该 数值模型既可以大幅提高非稳态渗流的计算效率, 还可以避免孤立裂隙所带来的影响. 最后, 通过两个工程算例验证该数值模型的适用性.     

爆破开挖扰动下锚固节理岩质边坡位移突变特征与能量机理

针对白鹤滩水电站左岸坝基河谷底部边坡岩体爆破开挖,采用现场岩体位移监测、锚索轴力监测及数值模拟的手段,研究了爆破开挖扰动下锚固节理岩质边坡的位移突变特征及其能量机理。研究结果表明:对于深切河谷底部高地应力边坡岩体爆破开挖,爆炸荷载挤压及地应力作用下,岩体所积聚的应变能快速释放,导致了节理岩质边坡的位移突变,突变位移包括节理张开位移和岩体回弹位移两部分;地应力水平越高、岩体弹性模量越低,总的突变位移量越大;预应力锚索主要通过抑制节理张开位移来控制边坡岩体的位移突变,锚索预应力等级越高,其吸能和释能速率越高,对节理岩体位移突变的控制效果越好,当锚索的预应力等级高到一定程度后,节理岩体的突变位移不再随锚索预应力等级的升高而显著减小。     

岩体结构非平衡演化的有效应力原理及长期稳定性分析

开挖卸荷后的天然岩体往往处于非平衡演化状态, 将直接影响岩体工程结构的正常运行、长期稳定和安全. 时效变形和损伤演化是岩体结构非平衡演化的核心. 在赖斯(Rice) 内变量热力学理论框架下, 提出了岩体结构非平衡演化的有效应力原理, 指出有效应力是总应力中能有效驱动结构演化的部分. 将内变量率形式的非弹性应变率方程和能量耗散率函数表示为有效应力形式, 并提出非弹性余能概念. 给定具体的余能密度函数和内变量演化方程, 得到了考虑损伤的内变量黏塑性应变率方程. 通过相似材料加卸载蠕变试验结果进行参数辨识, 并分别计算了内变量率形式和有效应力形式的黏塑性应变率、能量耗散率和非弹性余能, 并对其进行比较分析. 结果表明:在过渡蠕变和稳态蠕变阶段两种形式的方程计算的黏塑性应变率几乎相等, 但在加速蠕变阶段两者相差较大;非弹性余能和能量耗散率全域积分分别从驱动结构非平衡演化的内在潜力和实际效果的角度表征了结构的非平衡演化状态和演化趋势, 能量耗散率积分更合适用于评价岩体工程结构的长期稳定性. 最后以深埋地下洞室作为工程算例, 并对其长期稳定性进行分析.      

??????????????????????о?

根据岩心实验提出了两参数的非线性渗流方程,并在此基础上结合Buckly-Leverett 方程建立了非线性非混相驱替方程. 通过对该方程求解分析可知,与线性达西渗流相比,非 线性渗流的驱油相率较差;同一含水饱和度下的含水率更高;而储层中各点的压力梯度更大, 导致所需注采压差更大. 说明有效的开发低渗透油藏需要比中高渗油藏更小的井距和更高的 注采压力体系.     

可储能原则在一类岩爆问题中的研究和应用

从能量的角度对一类岩爆问题--应变型岩爆进行一些探讨.对此类岩爆能量的来源、能量的去处及能量释放的原因和大小进行了完整的论述.提出了岩体的可储能原则--岩体可存储的能量并非定值,而是随应力状态的变化而改变,当应力发生突变,岩体内实际存储的能量大于可储能时,多余的能量必将以各种形式瞬间释放,其中包括动能.根据此准则对简化了的一维、三维模型进行了详细的分析和计算,并结合离散元软件3DEC再现了三维卸载岩爆的破坏过程.     

Bifurcations of a cantilevered pipe conveying steady fluid with a terminal nozzle

This paper studies interactions of pipe and fluid and deals with bifurcations of a cantilevered pipe conveying a steady fluid, clamped at one end and having a nozzle subjected to nonlinear constraints at the free end. Either the nozzle parameter or the flow velocity is taken as a variable parameter. The discrete equations of the system are obtained by the Ritz-Galerkin method. The static stability is studied by the Routh criteria. The method of averaging is employed to examine the analytical results and the chaotic motions. Three critical values are given. The first one makes the system lose the static stability by pitchfork bifurcation. The second one makes the system lose the dynamical stability by Hopf bifurcation. The third one makes the periodic motions of the system lose the stability by doubling-period bifurcation. The project supported by the Science Foundation of Tongji University and Tongji University and National Key Projects of China under Grant No. PD9521907.     

输液管的非线性振动、分叉与混沌——现状与展望

围绕输液管的非线性振动、分叉与混沌问题，对近几年来的主要研究工作加以综述并提出预测，其中包括运动方程中非线性项的归纳与讨论、非线性动力分析的一些现代计算方法、定常流速下输液管的分叉与混沌行为、振荡流速下输液管的参数共振以及今后值得进一步研究的某些问题．     

Multiple scales analysis for double Hopf bifurcation with?1:3?resonance

The dynamical behavior of a general n-dimensional delay differential equation (DDE) around a 1:3 resonant double Hopf bifurcation point is analyzed. The method of multiple scales is used to obtain complex bifurcation equations. By expressing complex amplitudes in a mixed polar-Cartesian representation, the complex bifurcation equations are again obtained in real form. As an illustration, a system of two coupled van der Pol oscillators is considered and a set of parameter values for which a 1:3 resonant double Hopf bifurcation occurs is established. The dynamical behavior around the resonant double Hopf bifurcation point is analyzed in terms of three control parameters. The validity of analytical results is shown by their consistency with numerical simulations.     

Nonlinear dynamical stability analysis of the circular three-dimensional frame

The three-dimensional frame is simplified into flat plate by the method of quasiplate. The nonlinear relationships between the surface strain and the midst plane displacement are established. According to the thin plate nonlinear dynamical theory, the nonlinear dynamical equations of three-dimensional frame in the orthogonal coordinates system are obtained. Then the equations are translated into the axial symmetry nonlinear dynamical equations in the polar coordinates system. Some dimensionless quantities different from the plate of uniform thickness are introduced under the boundary conditions of fixed edges, then these fundamental equations are simplified with these dimensionless quantities. A cubic nonlinear vibration equation is obtained with the method of Galerkin. The stability and bifurcation of the circular three-dimensional frame are studied under the condition of without outer motivation. The contingent chaotic vibration of the three-dimensional frame is studied with the method of Melnikov. Some phase figures of contingent chaotic vibration are plotted with digital artificial method.     

A Perturbation Method for Nonlinear Two-Dimensional Maps

A new perturbation method for a weakly nonlinear two-dimensional discrete-time dynamical system is presented. The proposed technique generalizes the asymptotic perturbation method that is valid for continuous-time systems and detects periodic or almost-periodic orbits and their stability. Two equations for the amplitude and the phase of solutions are derived and their fixed points correspond to limit cycles for the starting nonlinear map. The method is applied to various nonlinear (autonomous or not) two-dimensional maps. For the autonomous maps we derive the conditions for the appearance of a supercritical Hopf bifurcation and predict the characteristics of the corresponding limit cycle. For the nonautonomous maps we show amplitude-response and frequency-response curves. Under appropriate conditions, we demonstrate the occurrence of saddle-node bifurcations of cycles and of jumps and hysteresis effects in the system response (cusp catastrophe). Modulated motion can be observed for very low values of the external excitation and an infinite-period bifurcation occurs if the external excitation increases. Analytic approximate solutions are in good agreement with numerically obtained solutions.     

Nonlinear vibration of slightly curved pipe with conveying pulsating fluid

The nonlinear governing motion equation of slightly curved pipe with conveying pulsating fluid is set up by Hamilton’s principle. The motion equation is discretized into a set of low dimensional system of nonlinear ordinary differential equations by the Galerkin method. Linear analysis of system is performed upon this set of equations. The effect of amplitude of initial deflection and flow velocity on linear dynamic of system is analyzed. Curves of the resonance responses about \(\varOmega \approx {\omega _\mathrm{{1}}}\) and \(\varOmega \approx \mathrm{{2}}{\omega _\mathrm{{1}}}\) are performed by means of the pseudo-arclength continuation technique. The global nonlinear dynamic of system is analyzed by establishing the bifurcation diagrams. The dynamical behaviors are identified by the phase diagram and Poincare maps. The periodic motion, chaotic motion and quasi-periodic motion are found in this system.     

一类生物流体力学连续系统的分岔研究

连续动力系统的非线性动力学研究,由于其应用的广泛性与问题的复杂性,近年来越来越受到重视。本文对一类生物流体力学中的连续系统-动脉局部狭窄时血液流动的分岔特性进行了研究,采用有限差分方法,将由偏微分方程组描述的边境动力系统约化为由常微分方程组描述的高维离散动力系统。求得了离散动力系统的平衡解并分析其稳定性,同时讨论了流场中变量空间分布的变化情况。求得了离散动力系统的前三个Lyapunov指数,以此作为系统是否发生混沌的判别条件。     

Bifurcation in most probable phase portraits for a bistable kinetic model with coupling Gaussian and non-Gaussian noises

The phenomenon of stochastic bifurcation driven by the correlated non-Gaussian colored noise and the Gaussian white noise is investigated by the qualitative changes of steady states with the most probable phase portraits. To arrive at the Markovian approximation of the original non-Markovian stochastic process and derive the general approximate Fokker-Planck equation (FPE), we deal with the non-Gaussian colored noise and then adopt the uni¯ed colored noise approximation (UCNA). Subsequently, the theoretical equation concerning the most probable steady states is obtained by the maximum of the stationary probability density function (SPDF). The parameter of the uncorrelated additive noise intensity does enter the governing equation as a non-Markovian effect, which is in contrast to that of the uncorrelated Gaussian white noise case, where the parameter is absent from the governing equation, i.e., the most probable steady states are mainly controlled by the uncorrelated multiplicative noise. Additionally, in comparison with the deterministic counterpart, some peculiar bifurcation behaviors with regard to the most probable steady states induced by the correlation time of non-Gaussian colored noise, the noise intensity, and the non-Gaussian noise deviation parameter are discussed. Moreover, the symmetry of the stochastic bifurcation diagrams is destroyed when the correlation between noises is concerned. Furthermore, the feasibility and accuracy of the analytical predictions are verified compared with those of the Monte Carlo (MC) simulations of the original system.

     

On the Resonant Generation of Weakly Nonlinear Stokes Waves in Regions with Fast Varying Topography and Free Surface Current

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｉｎｔｅｒａｃｔｉｏｎｏｆｓｕｒｆａｃｅｗａｔｅｒｗａｖｅｓｗｉｔｈａｍｂｉｅｎｔｃｕｒｒｅｎｔｓａｎｄｕｎｄｕｌａｔｉｎｇｓｅａｂｅｄｔｏｐｏｇｒａｐｈｙｉｓｏｆｆｕｎｄａｍｅｎｔａｌｉｍｐｏｒｔａｎｃｅｔｏｃｏａｓｔａｌｅｎｇｉｎｅｅｒｓａｎｄｓｅｄｉｍｅｎｔｏｌｏｇｉｓｔｓ.Ｆｏｒｅｘａｍｐｌｅ,ｔｈｅｒｅｓｏｎａｎｔｇｅｎｅｒａｔｉｏｎｏｆｓｕｒｆａｃｅｗａｖｅｓｉｎａｃｕｒｒｅｎｔｏｖｅｒｏｔｈｅｒｔｉｄａｌｌｙｏｒｗａｖｅｉｎｄｕｃｅｄｂｅｄｆｏｒｍｓ,ｓ…     

NONLINEAR DYNAMICAL BIFURCATION AND CHAOTIC MOTION OF SHALLOW CONICAL LATTICE SHELL

The nonlinear dynamical equations of axle symmetry are established by the method of quasi-shells for three-dimensional shallow conical single-layer lattice shells. The compatible equations are given in geometrical nonlinear range. A nonlinear differential equation containing the second and the third order nonlinear items is derived under the boundary conditions of fixed and clamped edges by the method of Galerkin. The problem of bifurcation is discussed by solving the Floquet exponent. In order to study chaotic motion, the equations of free oscillation of a kind of nonlinear dynamics system are solved. Then an exact solution to nonlinear free oscillation of the shallow conical single-layer lattice shell is found as well. The critical conditions of chaotic motion are obtained by solving Melnikov functions, some phase planes are drawn by using digital simulation proving the existence of chaotic motion.