First-order gradient damage theory

Taking the strain tensor, the scalar damage variable, and the damage gradient as the state variables of the Helmholtz free energy, the general expressions of the firstorder gradient damage constitutive equations are derived directly from the basic law of irreversible thermodynamics with the constitutive functional expansion method at the natural state. When the damage variable is equal to zero, the expressions can be simplified to the linear elastic constitutive equations. When the damage gradient vanishes, the expressions can be simplified to the classical damage constitutive equations based on the strain equivalence hypothesis. A one-dimensional problem is presented to indicate that the damage field changes from the non-periodic solutions to the spatial periodic-like solutions with stress increment. The peak value region develops a localization band. The onset mechanism of strain localization is proposed. Damage localization emerges after damage occurs for a short time. The width of the localization band is proportional to the internal characteristic length.     

塑性损伤本构模型的自由能势函数研究

自由能势函数是不可逆热力学量及其共轭量的标量泛函,是广义标准材料模型的重要组成.论文根据损伤状态、过程、作用和伪势的定义和定理,分析了损伤耗散势的数学和物理含义,探讨了损伤耗散势的理论基础和描述耗散势的一般理论方法.研究了耗散势、损伤作用和自由能势函数的相互关系,给出了满足上述基本原理的自由能表达式,基于增量最小原理的变分形式得到了损伤理论基本方程的数值离散方程,实验模拟验证了该文理论和数值方法的正确性和有效性.     

Enhancing the Bolzon–Schrefler–Zienkiewicz Constitutive Model for Partially Saturated Soil

This paper presents an enhanced version of the elasto-plastic model for partially saturated soil first proposed by Bolzon, Schrefler and Zienkiewicz in 1996, “BSZ” model, which uses the effective stress tensor and suction as independent stress variables. It is recalled that the effective stress tensor proposed by Lewis and Schrefler in 1982 is thermodynamically consistent and, compared with other choices of stress tensors, results particularly suitable for partially saturated soil mechanics. A hydraulic constitutive relationship and a hydraulic hysteresis are introduced in the model, to take into account the irreversible deformation during cyclic drying and wetting until structural collapse. For this reason the plastic rate of strain is split into the sum of two components: one depending on the effective stress tensor and the other one on suction. This is the new feature of the BSZ model. This enhanced model is then cast into a thermodynamical framework at macroscopic level and it is shown that it is possible to derive the constitutive law from the Helmholtz free energy and a dissipation function, both for associative and non- associative plasticity. Finally the model predictions have been compared with experimental data for Sion slime, with particular emphasis on the deviatoric part, and model predictions of hysteretic behaviour have been investigated in case of a wetting and drying cycle on compacted betonite–kaolin.     

AN ANALYTICAL SOLUTION OF SECOND－ORDER WAVE FORCE ON A VERTICAL CIRCULAR CYLINDER

ＡＮＡＮＡＬＹＴＩＣＡＬＳＯＬＵＴＩＯＮＯＦＳＥＣＯＮＤ－ＯＲＤＥＲＷＡＶＥＦＯＲＣＥＯＮＡＶＥＲＴＩＣＡＬＣＩＲＣＵＬＡＲＣＹＬＩＮＤＥＲＺｈｏｕＺｈｉ－ｌｉ（邹志利）（ＤａｌｉａｎＵｎｉｖｅｒｓｉｔｙｏｆＴｅｃｈｎｏｌｏｇｙ,Ｄａｌｉａｎ）Ｄａｉ...     

Width effects on hydrodynamics of pendulum wave energy converter

Based on two- and three-dimensional potential flow theories, the width effects on the hydrodynamics of a bottom-hinged trapezoidal pendulum wave energy converter are discussed. The two-dimensional eigenfunction expansion method is used to obtain the diffraction and radiation solutions when the converter width tends to be infinity. The trapezoidal section of the converter is approximated by a rectangular section for simplification. The nonlinear viscous damping effects are accounted for by including a drag term in the two- and three-dimensional methods. It is found that the three- dimensional results are in good agreement with the two-dimensional results when the converter width becomes larger, especially when the converter width is infinity, which shows that both of the methods are reasonable. Meantime, it is also found that the peak value of the conversion efficiency decreases as the converter width increases in short wave periods while increases when the converter width increases in long wave periods.     

An approximate analytical solution of the cylindrical detonation wave generated by the linear explosion

The solution of the cylindrical detonation wave generated by the linear explosion was obtained by numerical method in ref.[1].In this paper,when the ratio of specific heatγ》1 by using the enlargement coordinate method,the first-order analytical solutions are obtained.The perturbation parameter isε=1/γ~2 .The correction of these solutions is checked at the end of this paper.     

一种冲击波作用下结构毁伤算法研究

针对爆炸冲击波与建筑物结构相互作用过程,分析了冲击波与结构碎块作用机理,发展了一种能够模拟建筑物结构破坏及冲击波传播过程的计算模型和方法。采用建筑物结构工程毁伤载荷作为判据,处理结构在冲击波作用下的破坏问题;利用流固耦合界面算法处理结构运动引起的泄压效应,利用“虚拟网格通气技术”处理结构碎块对冲击波的阻碍作用,模拟了冲击波作用下典型建筑物的毁伤过程及冲击波传播过程。结果表明,该模型在模拟冲击波与结构的作用过程中,压力计算结果与非结构动网格模拟结果符合较好;在典型建筑物毁伤过程的数值模拟中,计算得到的建筑物毁伤效果和冲击波超压分布与建筑物物理毁伤特点符合。

     

Influences of microscope size effect of matrix on plastic potential and void growth of porous materials

Based on approximate theoretical analyses on a typical spherical cell containing a spherical microvoid, the influences of matrix materials' microscopic scale on the macroscopic constitutives potential theory of porous material and microvoid growth have been investigated in detail. By assuming that the plastic deformation behavior of matrix materials follows the strain gradient (SG) plastic theory involving the stretch and rotation gradients, the ratio (λ=l/a) of the matrix materials' intrinsic characteristic lengthl to the microvoid radiusa is introduced into the plastic constitutives potential and the void growth law. The present results indicate that, when the radiusa of microvoids is comparable with the intrinsic characteristic lengthl of the matrix materials, the influence of microscopic size effect on neither the constitutive potential nor the microvoid evolution predicted can be ignored. And when the void radiusa is much lager than the intrinsic characteristic lengthl of the matrix materials, the present model can retrogress automatically to the improved Gurson model that takes into account the strain hardening effect of matrix materials. Project supported by the National Natural Science Foundation of China (No. A10102006).     

Perturbation solution to the nonlinear problem of oblique water exit of an axisymmetric body with a large exit-angle

In this paper,a nonlinear,unsteady3-D free surface problem of the oblique water exitof an axisymmetric body with a large water exit-angle was investigated by means of theperturbation method in which the complementary angle a of the water exit angle waschosen as a small parameter.The original3-D problem was solved by expanding it into apower series of a and reduced to a number of2-D problems.The integral expressions forthe first three order solutions were given in terms of the complete elliptic functions of thefirst and second kinds.The zeroth-order solution didn‘t turn out to be a linear problem asusual but a nonlinear one corresponding to the vertical water exit for the same body.Computational results were presented for the free surface shapes and the forces exerted upto the second order during the oblique water exit of a series of ellipsoids with various ratiosof length to diameter at different Froude numbers.     

混凝土动态冲击问题的一种欧拉数值方法

针对混凝土的动态冲击问题给出了混凝土的Holmquist Johnson Cook本构关系 ,本构关系中的等效强度作为压力、应变率和损伤的函数 ;压力表示为体积应变的函数 ,并考虑了永久破碎的效应 ;积累损伤作为等效塑性应变、塑性体积应变和压力的函数。结合三维弹塑性流体力学欧拉数值方法 ,提出了混凝土的Holmquist Johnson Cook本构关系与欧拉程序相结合的计算方法。介绍了带混凝土Holmquist Johnson Cook本构关系的三维弹塑性流体力学欧拉程序。     

A combined first principles and analytical determination of the modulus of cohesion,surface energy,and the additional constants in the second strain gradient elasticity

Mindlin’s (1965) second strain gradient theory due to its competency in capturing the effects of edges, corners, and surfaces is of particular interest. Formulation in this framework, in addition to the usual Lamé constants, requires the knowledge of sixteen additional materials constants. To date, there are no successful experimental techniques for measuring these material parameters which reflect the discrete nature of matter. The present work gives an accurate remedy for the atomistic calculations of these parameters by utilizing the first principles density functional theory (DFT) for the calculations of the atomic force constants combined with an analytical formulation. It will be shown that writing the consistency conditions obtained from the equivalency between the atomistic crystal lattice dynamics of the bulk material and its counterpart in the second strain gradient elasticity is insufficient for the calculations of all the additional constants. As it will be discussed, there are two missing conditions which are then provided by consideration of the free standing film problem that bring the surface effect into account. As a consequence of surface effect consideration, the modulus of cohesion which is one of the important additional constants is calculated. Moreover, an analytical expression for the surface energy in terms of the modulus of cohesion, Lamé constants, materials characteristic lengths, and the film thickness is presented. If the film thickness is much bigger than the magnitude of the characteristic lengths of the material, then the surface energy would no longer depend on the film thickness.     

Parallel finite element simulation of mooring forces on floating objects

The coupling between the equations governing the free‐surface flows, the six degrees of freedom non‐linear rigid body dynamics, the linear elasticity equations for mesh‐moving and the cables has resulted in a fluid‐structure interaction technology capable of simulating mooring forces on floating objects. The finite element solution strategy is based on a combination approach derived from fixed‐mesh and moving‐mesh techniques. Here, the free‐surface flow simulations are based on the Navier–Stokes equations written for two incompressible fluids where the impact of one fluid on the other one is extremely small. An interface function with two distinct values is used to locate the position of the free‐surface. The stabilized finite element formulations are written and integrated in an arbitrary Lagrangian–Eulerian domain. This allows us to handle the motion of the time dependent geometries. Forces and momentums exerted on the floating object by both water and hawsers are calculated and used to update the position of the floating object in time. In the mesh moving scheme, we assume that the computational domain is made of elastic materials. The linear elasticity equations are solved to obtain the displacements for each computational node. The non‐linear rigid body dynamics equations are coupled with the governing equations of fluid flow and are solved simultaneously to update the position of the floating object. The numerical examples includes a 3D simulation of water waves impacting on a moored floating box and a model boat and simulation of floating object under water constrained with a cable. Copyright © 2003 John Wiley & Sons, Ltd.     

爆轰波碰撞聚能无网格MPM法数值模拟

爆轰聚能经常被应用于工程爆破切割,为了探索新的聚能爆破形式,本文应用无网格MPM法对圆柱筒装药和导爆索多点起爆形式的爆轰波碰撞聚能问题进行三维数值模拟,采用显式积分算法完成爆轰波形成和汇聚过程以及爆轰压力场分布的数值计算,并将数值模拟结果与实验结果做了对比分析。由此验证了无网格MPM法的准确性,也为炸药聚能爆破提供了新的思路和高效的操作方法。     

大跨度钢桁桥模型的精细化损伤定位模拟和试验研究

在役大跨钢桁桥数日众多,其杆件往往较长,所以识别出损伤杆并找出局部损伤所处位置对安全评定和加固尤为重要。针对这一问题,参考我国在役钢桁梁桥,基于健康监测的试验目的和相似理论设计制作了贝雷梁式筒支钢桁桥Benchmark模型。根据其结构形式、具体尺寸,并经试验结果修正,建立了基准Matlab有限元模型。环境激励下,首先隔...     

Characterization of long time fluctuations of forces exerted on an oscillating circular cylinder at KC=10

Flow dynamics, in-line and transverse forces exerted on an oscillating circular cylinder in a fluid initially at rest are studied by numerical resolution of the two-dimensional Navier-Stokes equations. The Keulegan-Carpenter number is held constant at KC=10 and Re is increased from 40 to 500. For the different flow regimes, links between flow spatio-temporal symmetries and force histories are established. Besides simulations of long duration show that in two ranges of Re, forces exhibit low frequency fluctuations compared to the cylinder oscillation frequency. Such observations have been only mentioned in the literature and are more deeply examined here. In both ranges, force fluctuations correspond to oscillations of the front and rear stagnation points on the cylinder surface. However, they occur in flow regimes whose basic patterns (V-shaped mode or diagonal mode) have different symmetry features, inducing two distinct behaviors. For 80≤Re≤100, fluctuations are related to a spectral broadening of the harmonics and to a permutation between three vortex patterns (V-shaped, transverse and oblique modes). In the second range 150≤Re≤280, amplitude fluctuations are correlated to the appearance of low frequency peaks interacting with harmonics of the cylinder frequency. Fluctuations are then a combination of a wavy fluctuation and an amplitude modulation. The carrier frequency corresponding to the wavy fluctuation depends on Re and is related to a fluid characteristic time; the modulation frequency is independent of Re and equal to 1/4 of the cylinder oscillation frequency.     

Study on wave dispersion characteristics of piezoelectric sandwich nanoplates considering surface effects

In this study, the wave propagation properties of piezoelectric sandwich nanoplates deposited on an orthotropic viscoelastic foundation are analyzed by considering the surface effects (SEs). The nanoplates are composed of a composite layer reinforced by graphene and two piezoelectric surface layers. Utilizing the modified Halpin-Tsai model, the material parameters of composite layers are obtained. The displacement field is determined by the sinusoidal shear deformation theory (SSDT). The Euler-Lagrange equation is derived by employing Hamilton’s principle and the constitutive equations of piezoelectric layers considering the SEs. Subsequently, the nonlocal strain gradient theory (NSGT) is used to obtain the equations of motion. Next, the effects of scale parameters, graphene distribution, orthotropic viscoelastic foundation, and SEs on the propagation behavior are numerically examined. The results reveal that the wave frequency is a periodic function of the orthotropic angle. Furthermore, the wave frequency increases with the increase in the SEs.

     

热载荷作用下梯度多孔材料梁弯曲及过屈曲力学行为的研究

基于Bernoulli-Euler梁理论,引入物理中面解耦了复合材料结构的面内变形与横向弯曲特性,研究了梯度多孔材料矩形截面梁在热载荷作用下的弯曲及过屈曲力学行为．假设沿梁厚度方向材料的性质是连续变化的,利用能量法推导了矩形截面梁的控制微分方程和边界条件,并用打靶法对无量纲化的控制方程进行数值求解．利用计算得到的结果分析了材料的性质、热载荷、边界条件对矩形截面梁非线性力学行为的影响．结果表明,对称材料模型下,固支梁与简支梁均显示出了典型的分支屈曲行为特征,而其临界屈曲热载荷值均会随着孔隙率系数的增加而单调增加．非对称材料模型下,固支梁仍显示出分支屈曲行为特征,但其临界屈曲热载荷不再随着孔隙率系数的变化而单调变化;而对于两端简支梁,发生了弯曲变形,弯曲挠度随载荷的增大而增大．     

An experimental study on the effects of uniform injection through one perforated surface of a square cylinder on some aerodynamic parameters

It is well known that injection/suction (transpiration) through a perforated surface is an efficient way of influencing the characteristics of a turbulent boundary layer. Injection application creates a thicker boundary layer on a flat plate and it thus decreases drag. In aeronautical applications, suction is frequently used to delay boundary layer separation. This paper presents an experimental study on the effects of uniform injection through one perforated surface of a square cylinder on the pressure distribution and drag coefficient in a two-dimensional turbulent flow. For this purpose, surface pressure measurements around a square cylinder have been performed at three different Reynolds numbers in a wind tunnel. The parameters taken into account were injection rate, position of perforated surface (i.e., front, top, and rear), and pressure coefficient and drag coefficient. The results show that variation in pressure coefficient around the square cylinder and drag coefficient were influenced by the position of perforated surface and by injection rate.     

On the influence of numerical schemes and subgrid–stress models on large eddy simulation of turbulent flow past a square cylinder

Influence of finite difference schemes and subgrid‐stress models on the large eddy simulation calculation of turbulent flow around a bluff body of square cylinder at a laboratory Reynolds number, has been examined. It is found that the type and the order of accuracy of finite‐difference schemes and the subgrid‐stress model for satisfactory results are dependent on each other, and the grid resolution and the Reynolds number. Using computational grids manageable by workstation‐level computers, with which the near‐wall region of the separating boundary layer cannot be resolved, central‐difference schemes of realistic orders of accuracy, either fully conservative or non‐conservative, suffer stability problems. The upwind‐biased schemes of third order and the Smagorinsky eddy‐viscosity subgrid model can give reasonable results resolving much of the energy‐containing turbulent eddies in the boundary layers and in the wake and representing the subgrid stresses in most parts of the flow. Noticeable improvements can be obtained by either using higher order difference schemes, increasing the grid resolution and/or by implementing a dynamic subgrid stress model, but each at a cost of increased computational time. For further improvements, the very small‐scale eddies near the upstream corners and in the laminar sublayers need to be resolved but would require a substantially larger number of grid points that are out of the range of easily accessible computers. Copyright © 2002 John Wiley & Sons, Ltd.     

含时滞反馈的楔式制动系统动力学分析

考虑含时滞反馈的影响,建立楔式制动系统动力学模型,运用多尺度方法对黏滑界面附近区域进行受迫主共振求解,分析时滞量、楔角与系统刚度对系统幅频响应的影响,应用Routh-Hurwitz判据分析系统稳定性的影响因素。基于解析解的分析表明:稳态幅值和稳定性边界都随时滞量发生周期性变化,周期内较大的时滞量引起鞍结分岔,并发展至不稳定多解;楔角和系统刚度增加引起主共振振幅增大,并扩大了不稳定区域。