黄河口裸置管线在位稳定性数值模拟研究

基于Biot固结理论,并结合接触算法,考虑了海水流经管线时的尾流效应,建立了波浪荷载作用下裸置管线与周围土体相互作用的有限元模型,使计算结果更加符合实际情况。将此模型应用于黄河口埕岛油田裸置管线在位稳定性计算,分析了尾流效应对裸置管线稳定性的影响,在此基础上提出了适用于黄河口不同地质分区的裸置管线稳定性校核曲线。     

Modeling wave-induced pore pressure and effective stress in a granular seabed

The response of a sandy seabed under wave loading is investigated on the basis of numerical modeling using a multi-scale approach. To that aim, the discrete element method is coupled to a finite volume method specially enhanced to describe compressible fluid flow. Both solid and fluid phase mechanics are upscaled from considerations established at the pore level. Model’s predictions are validated against poroelasticity theory and discussed in comparison with experiments where a sediment analog is subjected to wave action in a flume. Special emphasis is put on the mechanisms leading the seabed to liquefy under wave-induced pressure variation on its surface. Liquefaction is observed in both dilative and compactive regimes. It is shown that the instability can be triggered for a well-identified range of hydraulic conditions. Particularly, the results confirm that the gas content, together with the permeability of the medium are key parameters affecting the transmission of pressure inside the soil.     

Instability of frictional contact states in infinite layers

The research reported in this paper is focused on the instability of equilibrium and steady sliding states of elastic orthotropic layers in the presence of unilateral obstacles with Coulomb friction with emphasis on a divergence type instability called directional instability of frictional contact states that cannot occur in isotropic layers. Analytic expressions and numerical solutions are provided for the instability mode and for the coefficient of friction at the onset of instability. A parametric study is done to investigate how this coefficient of friction and the instability mode vary with changes of the system parameters. For certain combinations of material data, significantly low coefficients of friction were required for the onset of instability. A finite element model that approximates the continuum and a lumped model that captures some of the features of the continuum are presented.     

地下爆炸冲击波引起的基底滑移隔震建筑的动力响应

对埋置在土体中的圆形衬砌内发生意外爆炸引起的附近地表的滑移隔震建筑 (考虑了土 结构相互作用 )的动力响应进行了研究。在土 结构相互作用分析中 ,分别采用有限元、无穷元和集中质点体系建立了衬砌周围的有限土区、四周的无限土区以及基底滑移隔震建筑的计算模型 ;在基底滑移隔震分析中 ,提出了用连续摩擦力模型取代传统的库仑摩擦力模型 ,避免跟踪啮合滑移的过渡边界并可减小计算误差。通过数值模拟 6层基底滑移隔震建筑在遭受地下爆炸冲击波作用且考虑土 结构相互作用所产生的动力响应 ,分析了基底滑移隔震效果以及不同隔震参数对隔震效果的影响 ,为相关的理论研究和工程实践提供了依据。     

Wave-current-induced soil response in marine sediments

This letter presents a new analytical approximation for wave (current)-induced dynamic soil response in marine sediments. In the model, the third-order approximation for wave-current interactions is employed for the flow model, while Biot's dynamic poro-elastic model is used to simulate the porous flow in a seabed. The newly analytical solution is validated with the field observations. Based on the solution, effects of currents and wave-nonlinearity on soil response are examined and a parametric study will be carried out to examine the influence of currents on the liquefaction potential.     

The frequency response of dynamic friction: Model comparisons

Any linearised theory of the initiation of friction-excited vibration via instability of the state of steady sliding requires information about the dynamic friction force in the form of a frequency response function for sliding friction. Recent measurements of this function for an interface consisting of a nylon pin against a glass disc are used to probe the underlying constitutive law. Results are compared to linearised predictions from the simplest rate-state model of friction, and a rate-temperature model. In both cases the observed variation with frequency is not compatible with the model predictions, although there are some significant points of similarity. The most striking result relates to variation of the normal load: any theory embodying the Coulomb relation F∝N would predict behaviour entirely at variance with the measurements, even though the steady friction force obtained during the same measurements does follow the Coulomb law.     

Analytical solutions of tsunamis generated by underwater earthquakes

Tsunamis induced by underwater earthquakes are theoretically analyzed by applying the linear potential theory. Special attention is placed on the initial state of tsunami. For instantaneous seabed deformations, analytical wave solutions induced by three fundamental seabed deformations at initial stage are derived rather than integral expressions in past studies. These analytical solutions constitute a fundamental base for analyzing waves generated by arbitrary seabed displacement with the help of Fourier analysis. Tsunamis induced by non-instantaneous seabed deformation are analyzed as well. For the sake of examining the contributions of all wave components involved in the tsunami waveform, the amplitude density is proposed to examine the effects of deformation width, water depth, harmonic mode and rising time on waveforms. Results show that a larger ratio of water depth to deformation width results in a greater difference between initial waveform and seabed deformation, and the effect of the rising time is significant in deeper-water configuration. For cosine and sine seabed liftings, the effects of higher harmonic modes might be ignored.     

On the interaction of water waves and seabed by the porous medium model

The interaction of water waves and seabed is studied by using Yamamoto's model, which takes into account the deformation of soil skeletal frame, compressibility of pore fluid flow as well as the Coulumb friction. When analyzing the propagation of three kinds of stress waves in seabed, a simplified dispersion relation and a specific damping formula are derived. The problem of seabed stability is further treated analytically based on the Mohr-Coulomb theory. The theory is finally applied to the coastal problems in the Lian-Yun Harbour and compared with observations and measurements in soil-wave tank with satisfactory results. The project supported by the National Science Foundation of China     

平面波输入下海水-海床-结构动力相互作用分析

海洋工程结构的地震反应分析是保证海洋工程结构地震安全的重要环节.由于其所处的复杂环境, 该问题涉及到流固耦合和土-结相互作用.本文基于海水、饱和海床、基岩流固耦合统一计算框架,采用Davidenkov模型和修正的Masing准则考虑饱和海床的非线性,在脉冲SV波垂直入射下, 进行了海域场地和海洋工程结构的动力响应分析. 首先,对比分析了线性自由场和非线性自由场输入情形的海域场地非线性反应,结果表明线性自由场输入时反应不合理,自由场分析和场地分析应该采用相一致的本构模型. 然后,对比分析了海床分别为线性和非线性情形时,海域场地以及海水-海床-结构体系的反应特征. 与线性海床情形相比,非线性对海床反应的影响主要由如下两方面因素控制: 一方面,非线性导致饱和海床模量减小, 饱和海床与基岩间的波阻抗比减小,由基岩到饱和海床间的反射系数和透射系数增加, 导致反应增大; 另一方面,非线性导致阻尼加大, 使海床反应减小. 对于本文算例而言,阻尼对非线性海床结果的影响占主导作用.      

Comparison and analysis of two Coulomb friction models on the dynamic behavior of slider-crank mechanism with a revolute clearance joint

The objective of this study is to investigate the effects of the Coulomb dry friction model versus the modified Coulomb friction model on the dynamic behavior of the slider-crank mechanism with a revolute clearance joint. The normal and tangential forces acting on the contact points between the journal and the bearing are described by using a Hertzian-based contact force model and the Coulomb friction models, respectively. The dynamic equations of the mechanism are derived based on the Lagrange equations of the first kind and the Baumgarte stabilization method. The frictional force is solved via the linear complementarity problem (LCP) algorithm and the trial-and-error algorithm. Finally, three numerical examples are given to show the influence of the two Coulomb friction models on the dynamic behavior of the mechanism. Numerical results show that due to the stick friction, the slider-crank mechanism may exhibit stick-slip motion and can balance at some special positions, while the mechanism with ideal joints cannot.     

Effects of bottom shear stresses on the wave-induced dynamic response in a porous seabed：PORO-WSSI （shear） model

When ocean waves propagate over the sea floor,dynamic wave pressures and bottom shear stresses exert on the surface of seabed.The bottom shear stresses provide a horizontal loading in the wave-seabed interaction system,while dynamic wave pressures provide a vertical loading in the system.However,the bottom shear stresses have been ignored in most previous studies in the past.In this study,the effects of the bottom shear stresses on the dynamic response in a seabed of finite thickness under wave loading will be examined,based on Biot’s dynamic poro-elastic theory.In the model,an "u-p" approximation will be adopted instead of quasi-static model that have been used in most previous studies.Numerical results indicate that the bottom shear stresses has certain influences on the wave-induced seabed dynamic response.Furthermore,wave and soil characteristics have considerable influences on the relative difference of seabed response between the previous model(without shear stresses) and the present model(with shear stresses).As shown in the parametric study,the relative differences between two models could up to 10% of p0,depending on the amplitude of bottom shear stresses.     

带有摩擦耗能元件的框架结构动力分析

本文提出粘性屈服模型来模拟摩擦耗能元件的力-速度关系,该模型是连续变化的,克服了库仑摩擦模型不连续导致数值计算复杂的缺点,在进行摩擦耗能体系的动力分析中,采用缩减自由度技术,并作适当的变换,则带有摩擦耗能元件体系的动力分析归结为求解微分代数方程,本文采用增量型Rosenbrock二级三阶半隐式Runge-Kutta法求解该方程,以考虑框架和支撑的材料和几何非线性。对带有摩擦耗能元件的钢框架进行了弹塑性动力分析,研究了支撑刚度与结构层刚度的比值、摩擦力的大小以及地震波类型等参数对体系的影响。     

Thermo-elastic dynamic instability (TEDI) in frictional sliding of two elastic half-spaces

In some simplified 1D models, we recently studied the coupling of TEI (thermoelastic instability) and DI (dynamic instability), finding that thermal effects can render unstable the otherwise neutrally stable natural elastodynamic modes of the system, giving rise to a new family of instability which we called TEDI.Here, we study the general case of two sliding elastic half-planes, finding again a relatively weak coupling between thermal and dynamic effects, and the general family of instability TEDI class is found to modify both the otherwise separated TEI and DI classes. The growth factor, the phase velocity and the migrating speeds of the perturbations are wavelength-dependent, and it is difficult to give a complete picture given the high number of materials’ parameters, and the dependence on speed, friction coefficient, and the underlying uniform pressure. However, a set of results are given for “large” and “small” mismatch of shear wave speeds in the materials, and as a function of (i) friction coefficient; (ii) sliding speed V₀; (iii) wavenumber parameter γ. In the case of small mismatch, generalized Rayleigh waves exists already under frictionless conditions, the critical f for instability is zero. DI dominates over TEI typically for large wavenumbers, where the growth factors increase without limit and hence become eventually meaningless, requiring regularizations for example with rate-state dependent friction laws. TEI growth factors vice versa have a maximum at a certain wavenumber and therefore are always well posed. Larger coupling effects are noticed for two materials with large mismatch, but significantly only for sliding speeds comparable with the wave speed. In general, TEI growth factors increase with speed, whereas DI growth factors increase with speed for similar materials and decrease when the mismatch between materials is large.     

The effects of Coulomb friction on the performance of centrifugal pendulum vibration absorbers

We investigate analytically and experimentally the effects of Coulomb friction on the performance of centrifugal pendulum vibration absorbers (CPVAs), which are used to reduce torsional vibrations in rotating machinery. The analysis is based on perturbation methods applied to the nonlinear equations of motion for a rotor subjected to an engine order applied torque and equipped with a circular path CPVA with viscous and Coulomb damping. The experimental work is based on quantifying parameters for the damping model using free vibration measurements with a viscous and Coulomb damping identification scheme that is enhanced to better handle measurement noise, and running tests for steady-state operation under a range of loading conditions. The level of Coulomb damping is varied by adjusting the friction of the absorber connection bearing. Good agreement is found between the analytical predictions and the experimental data. It is shown that the absorber sticks up to a level of excitation that allows it to release, after which the Coulomb damping acts in the expected manner, resulting in lowered response amplitudes. The results obtained are of general use in assessing absorber performance when dry friction is present in absorber suspensions.     

浅水区海底管道周围海床渗透压力计算研究

我国海上油田开采起步较晚,大部分油田处于浅水区,因此,在设计管道时,应充分考虑由浅水区波浪引起的管道周围海床渗流力。根据浅水波相关假设,考虑自由水面非线性影响,推导出椭圆余弦波的波面方程,在此基础上进一步得到一个关于速度势的表达式,并根据该表达式得出作用于海床表面的波压公式。考虑海床土的压缩性,推导出一阶近似椭圆余弦波作用下浅水区埋置管道周围海床的渗流压力解析解,最后将计算结果与大型水槽试验及以往研究成果作对比。结果表明,在椭圆余弦波的作用下,由一阶椭圆余弦波理论得到的计算结果与试验结果规律基本一致,与相似工况下的现有理论成果数值基本相同,具有一定的可行性和工程价值。     

A set-valued force law for spatial Coulomb–Contensou friction

The aim of this paper is to develop a set-valued contact law for combined spatial Coulomb–Contensou friction, taking into account a normal friction torque (drilling friction) and spin. The set-valued Coulomb–Contensou friction law is derived from a non-smooth velocity pseudo potential. A higher-order Runge–Kutta time-stepping method is presented for the numerical simulation of rigid bodies with Coulomb–Contensou friction. The algebraic inclusion describing the contact problem is solved with an Augmented Lagrangian approach. The theory and numerical methods are applied to the Tippe-Top. The analysis and numerical results on the Tippe-Top illustrate the importance of Coulomb–Contensou friction for the dynamics of systems with friction.     

Experimental evidence of flutter and divergence instabilities induced by dry friction

Flutter and divergence instabilities have been advocated to be possible in elastic structures with Coulomb friction, but no direct experimental evidence has ever been provided. Moreover, the same types of instability can be induced by tangential follower forces, but these are commonly thought to be of extremely difficult, if not impossible, practical realization. Therefore, a clear experimental basis for flutter and divergence induced by friction or follower-loading is still lacking. This is provided for the first time in the present article, showing how a follower force of tangential type can be realized via Coulomb friction and how this, in full agreement with the theory, can induce a blowing-up vibrational motion of increasing amplitude (flutter) or an exponentially growing motion (divergence). In addition, our results show the limits of a treatment based on the linearized equations, so that nonlinearities yield the initial blowing-up vibration of flutter to reach eventually a steady state. The presented results give full evidence to potential problems in the design of mechanical systems subject to friction, open a new perspective in the realization of follower-loading systems and of innovative structures exhibiting ‘unusual’ dynamical behaviors.     

Wave induced sliding at the interface between a layered elastic medium and a half-space

The sliding interface between an unrestrained elastic half-space and a grounded layered half-space excited by an incident harmonic wave is investigated. The problem is formulated considering various possible boundary conditions and boundary inequalities at the sliding interface. The Coulomb friction model without distinction between the static and kinetic coefficients of friction is considered to govern the sliding condition. Three possible bands at the interface, namely slip, stick, and separation, are considered. The interface is assumed to be preloaded under normal and shear stresses. The solution is developed by modifying the problem of welded interface, which then is reduced to a set of algebraic equations. The effects of the incident angle, layer thickness, friction coefficient and externally applied stresses on the drift velocity of the unrestrained half-space are studied numerically for a pair of materials. It is shown that the sliding interface, and hence the drift velocity of unrestrained half-space is noticeably influenced by the layered medium. These results are expected to be useful for the development of a new kind of ultrasonic drive in future.     

Non-linear dynamic response of a rotating radial Timoshenko beam with periodic pulse loading at the free-end

Consideration is given to the dynamic response of a Timoshenko beam under repeated pulse loading. Starting with the basic dynamical equations for a rotating radial cantilever Timoshenko beam clamped at the hub in a centrifugal force field, a system of equations are derived for coupled axial and lateral motions which includes the transverse shear and rotary inertia effects, as well. The hyperbolic wave equation governing the axial motion is coupled with the flexural wave equation governing the lateral motion of the beam through the velocity-dependent skew-symmetric Coriolis force terms. In the analytical formulation, Rayleigh-Ritz method with a set of sinusoidal displacement shape functions is used to determine stiffness, mass and gyroscopic matrices of the system. The tip of the rotating beam is subjected to a periodic pulse load due to local rubbing against the outer case introducing Coulomb friction in the system. Transient response of the beam with the tip deforming due to rub is discussed in terms of the frequency shift and non-linear dynamic response of the rotating beam. Numerical results are presented for this vibro-impact problem of hard rub with varying coefficients of friction and the contact-load time. The effects of beam tip rub forces transmitted through the system are considered to analyze the conditions for dynamic stability of a rotating blade with intermittent rub.     

Limit analysis and homogenization of porous materials with Mohr–Coulomb matrix. Part I: Theoretical formulation

The present two-part study aims at investigating the specific effects of Mohr–Coulomb matrix on the strength of ductile porous materials by using a kinematic limit analysis approach. While in the Part II, static and kinematic bounds are numerically derived and used for validation purpose, the present Part I focuses on the theoretical formulation of a macroscopic strength criterion for porous Mohr–Coulomb materials. To this end, we consider a hollow sphere model with a rigid perfectly plastic Mohr–Coulomb matrix, subjected to axisymmetric uniform strain rate boundary conditions. Taking advantage of an appropriate family of three-parameter trial velocity fields accounting for the specific plastic deformation mechanisms of the Mohr–Coulomb matrix, we then provide a solution of the constrained minimization problem required for the determination of the macroscopic dissipation function. The macroscopic strength criterion is then obtained by means of the Lagrangian method combined with Karush–Kuhn–Tucker conditions. After a careful analysis and discussion of the plastic admissibility condition associated to the Mohr–Coulomb criterion, the above procedure leads to a parametric closed-form expression of the macroscopic strength criterion. The latter explicitly shows a dependence on the three stress invariants. In the special case of a friction angle equal to zero, the established criterion reduced to recently available results for porous Tresca materials. Finally, both effects of matrix friction angle and porosity are briefly illustrated and, for completeness, the macroscopic plastic flow rule and the voids evolution law are fully furnished.