Penalty finite element method for nonlinear dynamic response of viscous fluid-saturated biphase porous media

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｐｏｒｏｕｓｍｅｄｉａｍｏｄｅｌｓｄｅｄｕｃｅｄｆｒｏｍｍｉｘｔｕｒｅｔｈｅｏｒｙａｒｅａｔｔｒａｃｔｉｎｇａｔｔｅｎｔｉｏｎｏｆｍｏｒｅａｎｄｍｏｒｅｒｅｓｅａｒｃｈｅｒｓｂｅｃａｕｓｅｉｔｉｓｂａｓｅｄｏｎｔｈｅｆｒａｍｅｏｆｃｏｎｔｉｎｕｕｍｍｅｃｈａｎｉｃｓ[1]ａｎｄｗｅｒｅｓｕｃｃｅｓｓｆｕｌｌｙｕｓｅｄｔｏｄｅｐｉｃｔｔｈｅｍｅｃｈａｎｉｃａｌｂｅｈａｖｉｏｒｓｏｆｓｏｉｌｓａｎｄｂｉｏｌｏｇｉｃａｌｓｏｆｔｔｉｓｓｕｅｓ[2 ,3].Ｗｉｔｈｔｈｉｓｔｙｐｅｏ…     

STUDIES OF MELNIKOV METHOD AND TRANSVERSAL HOMOCLINIC ORBITS IN THE CIRCULAR PLANAR RESTRICTED THREE－BODY PROBLEM

ＳＴＵＤＩＥＳ　ＯＦ　ＭＥＬＮＩＫＯＶ　ＭＥＴＨＯＤ　ＡＮＤ　ＴＲＡＮＳＶＥＲＳＡＬ　ＨＯＭＯＣＬＩＮＩＣ　ＯＲＢＩＴＳ　ＩＮ　ＴＨＥ　ＣＩＲＣＵＬＡＲ　ＰＬＡＮＡＲ　ＲＥＳＴＲＩＣＴＥＤ　ＴＨＲＥＥ－ＢＯＤＹ　ＰＲＯＢＬＥＭＺｈｕＲｕｚｅｎｇ（朱如...     

考虑耦合质量影响的饱和多孔介质动力响应分析的显式有限元法

根据Biot饱和多孔介质动力方程,采用解耦技术,提出了考虑耦合质量Pd影响的饱和多孔介质中动力响应分析的显式有限元法。文中建立并推导了显式有限元的公式,编制了相应的计算程序并进行了实例计算。计算结果与解析解进行了对比,两者符合很好,表明本文方法是处理饱和多孔介质动力问题的一种有效方法。文中还分析了耦合质量ρa对固相和液相动位移的影响。     

集中载荷作用下圆拱梁的静分叉问题

采用有限元法中的伪弧长算法研究了集中载荷作用下圆拱的大范围非线性问题,给出了临界载荷与圆心角间的关系曲缄以及极值分叉与简单分叉的分界点,并对屈曲后的变形进行了追踪,文中首先简述伪弧长算法,然后给出了计算结果。     

A finite amplitude analysis of tunnel deformation by the finite element method with highly viscous fluid

A finite element method for highly viscous fluid is used to calculate the velocity and stress fields in the surrounding soft rock of a tunnel. In order to fit the calculated values with the measured displacement of tunnel wall, we inverted the boundary forces and the mechanical parameters of the surrounding rocks.     

Studying the stability of equilibrium solutions in a planar circular restricted four-body problem

The Newtonian circular restricted four-body problem is considered. We obtain nonlinear algebraic equations determining equilibrium solutions in the rotating frame and find six possible equilibrium configurations of the system. Studying the stability of equilibrium solutions, we prove that the radial equilibrium solutions are unstable, while the bisector equilibrium solutions are stable in Lyapunov’s sense if the mass parameter satisfies the conditions μ ∈ (0, μ₀, where μ₀ is a sufficiently small number, and μ ≠ μ_j, j = 1, 2, 3. We also prove that, for μ = μ₁ and μ = μ₃, the resonance conditions of the third order and the fourth order, respectively, are satisfied and, for these values of μ, the bisector equilibrium solutions are unstable and stable in Lyapunov’s sense, respectively. All symbolic and numerical calculations are done with the Mathematica computer algebra system. Published in Neliniini Kolyvannya, Vol. 10, No. 1, pp. 66–82, January–March, 2007.     

The bending of elastic circular ring of non-homogeneous and variable cross section under the actions of arbitrary loads

On the basis of the stepped reduction method suggested in [1], we investigate the problem of the bending of elastic circular ring of non-homogeneous and variable cross section under the actions of arbitrary loads. The general solution of this problem is obtained so that it can be used for the calculations of strength and rigidity of practical problems such as arch, tunnel etc. In order to examine results of this paper and explain the application of this new method, an example is brought out at the end of this paper. Circular ring and arch are commonly used structures in engineering. Timoshenko, S.^[2], Barber, J. R.^[3], Tsumura Rimitsu^[4] et al. have studied these problems of bending, but, so far as we know, it has been solely restricted to the general solution of homogeneous uniform cross section ring. The only known solution for the problems with variable cross section ones has been solely restricted to the solution of special case of flexural rigidity in linear function of coordinates. On account of fundamental equations of the non-homogeneous variable cross section problem being variable coefficients, it is very difficult to solve them. In this paper, we use the stepped reduction method suggested in [1] to transform the variable coefficient differential equation into equivalent constant coefficient one. After introducing virtual internal forces, we obtain general solution of an elastic circular ring with non-homogeneity and variable cross section under the actions of arbitrary loads.     

The annular twist mapping of the restricted three-body problem

This paper presents a new derivation of the twist mapping in the planar restricted problem. It differs from other treatments in the use of a novel canonical transformation which allows for the utilization of symplectic reduction techniques.     

辛解析奇异元在动荷载作用下V型切口问题中的应用

利用辛解析奇异单元,结合时域精细算法,研究了动荷载作用下的含平面V型切口问题。时域上,采用时域精细算法,并结合自适应算法控制展开项数,保证了计算精度。空间域上,切口尖端附近采用辛解析奇异单元,其余区域采用常规有限单元,避免了局部网格加密。本文使用的辛解析奇异单元不需要过渡单元和局部网格加密,且能够通过奇异单元内部的参数关系直接给出切口尖端的应力强度因子,不需要复杂的后处理过程。数值结果表明,本文方法具有良好的精度和稳定性,可以准确地计算动态应力强度因子。     

Self-similar solutions of the problem of formation of a hydraulic fracture in a porous medium

The problem of hydraulic fracture formation in a porous medium is investigated in the approximation of small fracture opening and inertialess incompressible Newtonian fluid fracture flow when the seepage through the fracture walls into the surrounding reservoir is asymptotically small or large. It is shown that the system of equations describing the propagation of the fracture has self-similar solutions of power-law or exponential form only. A family of self-similar solutions is constructed in order to determine the evolution of the fracture width and length, the fluid velocity in the fracture, and the length of fluid penetration into the porous medium when either the fluid flow rate or the pressure as a power-law or exponential function of time is specified at the fracture entrance. In the case of finite fluid penetration into the soil the system of equations has only a power-law self-similar solution, for example, when the fluid flow rate is specified at the fracture entrance as a quadratic function of time. The solutions of the self-similar equations are found numerically for one of the seepage regimes.     

Finite element analysis of wave propagation in fluid-saturated porous media

Ｔｈｅｄｙｎａｍｉｃｔｒａｎｓｉｅｎｔｒｅｓｐｏｎｓｅａｎａｌｙｓｉｓｏｆｐｏｒｏｕｓｍｅｄｉａｐｌａｙｓａｖｅｒｙｉｍｐｏｒｔａｎｔｒｏｌｅｉｎａｌｏｔｏｆｅｎｇｉｎｅｅｒｉｎｇｐｒａｃｔｉｃｅｓｓｕｃｈａｓｔｒａｎｓｉｅｎｔｃｏｎｓｏｌｉｄａｔｉｏｎ,ｎｏｉｓｅｃｏｎｔｒｏｌ,ｅａｒｔｈｑｕａｋｅｅｎｇｉｎｅｅｒｉｎｇａｎｄｂｉｏｅｎｇｉｎｅｅｒｉｎｇ．Ｂｉｏｔ［１］ｏｒｉｇｉｎａｌｌｙｄｉｓｃｕｓｓｅｄｔｈｅｗａｖｅｐｒｏｐａｇａｔｉｏｎｐｒｏｂｌｅｍｉｎｆｌｕｉｄ＿ｓａｔｕｒａｔｅｄｐｏ…     

The equations of motion of a system under the action of the impulsive constraints

In order to solve the problem of motion for the system with n degrees of freedom under the action of p impulsive constraints, we must solve the simultaneous equations consisting of n+p equations. In this paper, it has been shown that the undetermined multipliers in the equations of impact can be cancelled for the cases of both the generalized coordinates and the quasi-coordinates. Thus there are only n-p equations of impact. Combining these equations with p impulsive constraint equations, we have simultaneous equations consisting of n equations. Therefore, only n equations are necessary to solve the problem of impact for the system subjected to impulsive constraints. The method proposed in this paper is simpler than ordinary methods.     

基于显式有限元方法的两相介质弹塑性动力反应计算分析

针对增量形式的流体饱和两相多孔介质弹塑性波动方程组,运用基于显式逐步积分格式的时域显式有限元方法对该波动方程组进行求解,并应用基于SMP破坏准则的弹塑性动力本构模型描述两相介质的动力反应性质,对两相介质在输入地震波作用下的弹塑性动力反应进行计算和分析,将计算结果与相应的弹性动力反应的计算结果进行对比;对本文应用的弹塑性...     

Analytical investigation of the fluid flow in the interface region between a porous medium and a clear fluid in channels partially filled with a porous medium

In this paper analytical solutions for the steady fully developed laminar fluid flow in the parallel-plate and cylindrical channels partially filled with a porous medium and partially with a clear fluid are presented. The Brinkman-extended Darcy equation is utilized to model the flow in a porous region. The solutions account for the boundary effects and for the stress jump boundary condition at the interface recently suggested by Ochoa-Tapia and Whitaker. The dependence of the velocity on the Darcy number and on the adjustable coefficient in the stress jump boundary condition is investigated. It is shown that accounting for a jump in the shear stress at the interface essentially influences velocity profiles.     

非饱和多孔介质有限元分析的基本控制方程与变分原理

本文在对问题研究现状进行阐述的基础上较系统地给出了骨架可变形非饱和多孔介质的全耦合分析模型，模型中考虑了孔隙气体，水（油）流动对介质力学性能的影响，多孔介质的饱和度，渗透系数与毛吸压力的关系，由实验给出，所导出的控制方程以固体骨架的位移与孔隙流体压力为基本未知量，由于问题的非自共轭特征，文中构造了非饱和介质动力问题的参数变分形式，并在此基础上给出有限元离散方程。     

基于显式有限元方法的两相多孔介质地震反应研究

基于流体饱和两相多孔介质的弹性波动方程组,运用显式逐步积分格式与局部透射人工边界相结合的时域显式有限元方法对该波动方程组进行求解,对两相多孔介质在输入地震波作用下的弹性动力反应进行了计算和分析,以揭示两相多孔介质弹性地震反应的规律和性质.计算结果表明:两相介质弹性地震反应时程的波形与入射地震波的波形相同,且弹性地震反应的峰值出现的时刻对应于入射地震波的峰值出现的时刻.本文的数值计算同时表明了时域显式有限元方法在进行两相多孔介质地震反应计算分析时的有效性.     

On the stability of the rotational motion of a rigid body having a liquid filled cavity under finite initial disturbance

In this paper the problem of the stability of rotational motion of a rigid body which has a liquid filled cavity and a fixed point is investigated without any approximation. Criteria of stability and instability under finite disturbance are obtained. The region of stability is found out explicitly.     

流体饱和两相多孔介质拟静态问题的有限元解法

给出基于混合物理论的流体饱和两相多孔介质模型,该模型由一可变形固体 一流体相组成。采用Ｇａｌｅｒｋｉｎ加权残值法导出求解拟静态问题的有限元公式,并编制了二维有限元程序。用程序分析了一维和二维问题,得到合理的结果。