一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Preface: Current research progress on mechanics of high speed rail

The construction and operation of high-speed rail（HSR）grid within the past two decades in China,in terms of the scale,may be possibly comparable to any national-wide construction in the history of China,even the Great Wall.By counting railways with commercial train service at the speed of200 km/h,China has the world’s longest HSR network with over 19 369.8 km of track in service today and this number     

ACTA MECHANICA SINICA 2014 ANNUAL SUBJECT INDEX

正~~     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).     

Physiology,mechanics, and rheology

Biological systems possess rather specialized mechanical properties, acquired as part and parcel of the evolutionary development of the system as a whole. Their optimization permits the system to function physiologically in the context of a biologically essential, but mechanically often widely varying environment with adequate efficiency. The system's environment is its source of food and shelter; it represents the space in which it forages or preys on other creatures and in which it has to defend itself against still others. Thus, the system has to develop an adequately pliant, rheologically matched, energy-use efficient, mechanical interface between it and its surroundings. This must be an interface that both effectively excludes, but also effectively admits, the external. Internally, as well, it has to adapt the mechanical properties of cell and connective tissue to physiological function and the efficient performance of useful work. This will be illustrated by way of examples. Blood rheology is briefly discussed and put into the context of clinical hemorheology and epithelial protection; and function, by way of a mucus coating or a mucociliary clearance system, is reviewed in some detail. The importance of all aspects of rheological matching is demonstrated.Delivered as a Keynote Lecture at the Golden Jubilee Conference of the British Society of Rheology and Third European Rheology Conference, Edinburgh, 3–7 September, 1990.     

人造地震动合成中的位移误差分析

三角函数级数法是合成人工地震动常用算法之一,但是通过对加速度积分求取位移时,却存在与零线漂移相类似的位移漂移现象。     

Rainwater rivulets on a cable subject to windFilets d'eau de pluie sur un câble soumis au vent

Rainwater rivulets appear on inclined cables of cable-stayed bridges when wind and rain occur simultaneously. In a restricted range of parameters this is known to cause vibrations of high amplitudes on the cable. The mechanism underlying this effect is still under debate but the role of rainwater rivulets is certain. We use a standard lubrication model to analyse the dynamics of a water film on a cylinder under the effect of gravity and wind load. A simple criterion is then proposed for the appearance and position of rivulets, where the Froude number is the control parameter. Experiments with several geometries of cylinder covered with water in a wind tunnel show the evolution of the rivulets with the Froude number. Comparison of the prediction by the model with these experimental data shows that the main mechanism of rivulet formation and positioning is captured. To cite this article: C. Lemaitre et al., C. R. Mecanique 334 (2006).     

Hermitian space-time finite elements for estuarine mass transport

Space-time finite element solutions of the convection–dispersion equation using higher-order nodal continuity and Hermitian polynomial shape functions are described. Five separate elements ranging from a complete linear element with C^0,0 nodal continuity to a complete first-order Hermitian element with C^1,1 nodal continuity are subjected to detailed analysis. Wave deformation analyses identify the source of leading or trailing edge oscillations, trailing edge oscillations being the major source of difficulty. These observations are confirmed by numerical experiments which further demonstrate the potential of higher-order nodal continuity. The performance of the complete first-order Hermitian element is quite satisfactory and measurably superior to the linear element.     

Cristallisation par onde acoustique : le cas de l'héliumCrystallization by acoustic waves: the case of helium

We establish a theoretical model to explain the nucleation of a crystal of helium by an acoustic over-pressure. We explain the interfacial laws for this ultra-fast cristallization, close to the sound speed. Assuming spherical symmetry and taking into account the experimental data, we recover the dynamics of the growth and melting during an over-pressure impulse. To cite this article: M. Ben Amar et al., C. R. Mecanique 331 (2003).     

Scale shifting laws from pico to macro in consecutive segments by use of transitional functions

Systems with parts that vary in size from pico to macro inclusive are vulnerable of being incapacitated when a single part fails owing to deterioration of material properties. The majority of system failure can be attributed to incompatibility of integrated parts that were designed individually for general purpose. Total reliability calls for all parts, small and large, to be compatible in life spans. Mass, when regarded as energized matter, can vary as a function of time. This, in retrospect, explains why non-equilibrium and non-homogeneity cannot be avoided for multiscale shifting laws. A consistent and scale invariant definition of energy dissipation gives rise to mass pulsation, a common mechanism that seems to be applicable to living and non-living organisms.Scale shifting laws are developed from the use of transitional functions that stand for the mass ratios related to absorption energies and dissipation energies . The notations j and j + 1 stand for two successive scales: pi-na, na-mi, and mi-ma. Hence, the mass ratios , and can be referred to as the transitional inhomogeneity coefficients. They make up the scale shifting laws . Connection of the accelerated test data at the different scales, say from pico to nano to micro to macro, can be made by application of the definition of a scale invariant energy density dissipation function.On physical grounds, the segmented non-equilibrium and non-homogeneous test data can be connected through a velocity dependent mass and energy relation. Energy and power efficiency are defined to explore the macroscopic experiences to those at the lower scales. The time evolution properties of the material can also be derived as a package to include the accelerated test data, a procedure normally referred to as validation. The separation of derive-first and test-later, can never be abridged without ambiguities. Hence, total reliability of a system with many parts is advocated by judiciously matching the nine primary variables consisting of the initial disorder sizes, the time rates, and increments of the absorbed and dissipated energy density. The nine controllable variables consisting of life span distribution, energy, and power efficiencies for the three scale ranges are of secondary consideration.     

Direct measurement of static yield properties of cohesive suspensions

A technique of yield stress investigation based upon the combined use of two devices (an applied stress rheometer and an instrument for measuring the propagation velocity of small amplitude, torsional shear waves) is described. Investigations into the low shear rate rheological properties of illitic suspensions are reported for shear rates, typically, in the range 10^–4— 10^–1 s^–1 under applied stresses in the range 0.01 — 10 Nm^–2 and involving shear strains between 10^–1 and 10^–4. Results are presented which demonstrate that the technique does not invoke the excessive structural disruption of material associated with applied shear rate based methods (direct and otherwise) and the widely encountered problem of wall slip at the surface of rotational measuring devices is avoided using miniature vane geometries. Results are compared with those obtained using smooth-walled cyclindrical measuring devices in both applied stress and applied shear rate instruments.Yield measurements are considered in relation to the structural properties of the undisturbed material state and shear moduli obtained by studying the propagation of small amplitude (10^–5 rad), high frequency (~ 300 Hz) torsional shear waves through the test materials are reported. Experimental techniques and instrument modifications to permit these measurements are described.     

An experimental study of exit pressures for polymer melts

A slit die apparatus is used to measure exit pressures for five different polymer melts. Viscosity data obtained from the same apparatus agree well with values obtained from a cone-and-plate rheometer or a capillary rheometer. Except for a PVC sample where thermal degradation was found to occur, the exit pressures obtained by linear extrapolation of the measured pressure profiles are all positive, and increase with increasing shear stress. The values of the first normal stress difference calculated according to the exit pressure theory are of the right order of magnitude and in some cases correlate satisfactorily with values measured in a cone-and-plate rheometer. However, the high sensitivity of the exit pressure values to the method of extrapolation and the wild scatter of exit pressure data for some materials make it difficult to use the exit pressure method as a routine procedure for accurate determination of the first normal stress difference.