Deformation of an ellipsoidal ferrogel sample in a uniform magnetic field

The elongation of a ferroelastic material sample (whose initial shape is a sphere or an ellipsoid of revolution) under the action of an external magnetic field is studied in an in approximation of small strains. For a sphere, there is a classical estimate obtained under the assumption that elongating in the direction of the field, it becomes a spheroid and the stress and strain fields remain uniform. In the present calculation, it is assumed that the body is an ellipsoid (a sphere in a particular case) only in the absence of an external field; the shape of the sample in the presence of a field is not specified in advance but is found from the condition of balance of surface forces (elastic and magnetic). For the spherical case, the problem is solved exactly: it is shown, that the contour of the deformed body is described by a third-order algebraic equation. The case where the initial configuration is an ellipsoid of revolution is studied numerically. It is shown that in all versions, the refined solution leads to an appreciable increase in the elongation of the sample compared to the classical estimate. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 3, pp. 153–164, May–June, 2005.     

汶川地震触发陕西省境内地质灾害灾险情特征

在实地调查、汶川地震震后次生地质灾害隐患排查成果的基础上,结合县市地质灾害调查数据,从灾情和险情两个方面,对汶川地震触发地质灾害灾险情的致灾灾险情、灾种、类型等特征进行了较系统的分析对比。进一步比较了陕西在汶川地震触发地质灾害受灾省份中的情况,并对地质灾害趋势做了初步预测。研究表明,灾险情集中分布在龙门山构造带向陕西境内的延伸部分,险情比灾情分布范围广。灾险情的灾种都以崩塌为主,其次为滑坡,地震对地质体的震动触发作用明显。机关学校是汶川地震触发陕西境内地质灾害的主要受灾对象。陕西受汶川地震触发地质灾害的灾险情轻于四川、甘肃,三省地质灾害震前均较为严重,地震触发地质灾害的影响具有长期性,未来陕西境内地质灾害发生的频率可能增大。     

Steady sliding frictional contact problem for a 2d elastic half-space with a discontinuous friction coefficient and related stress singularities

The steady sliding frictional contact problem between a moving rigid indentor of arbitrary shape and an isotropic homogeneous elastic half-space in plane strain is extensively analysed. The case where the friction coefficient is a step function (with respect to the space variable), that is, where there are jumps in the friction coefficient, is considered. The problem is put under the form of a variational inequality which is proved to always have a solution which, in addition, is unique in some cases. The solutions exhibit different kinds of universal singularities that are explicitly given. In particular, it is shown that the nature of the universal stress singularity at a jump of the friction coefficient is different depending on the sign of the jump.     

Hydraulic analysis of free-surface flows with solidification

Summary The paper is concerned with a one-dimensional analysis of plane open-channel flow with continuous solidification. The process is of relevance for recent developments in the casting of steel and other metals. The bottom of the channel consists of a rotating casting roll and a horizontal cooling table, where the solidified material is withdrawn with given velocity. The study is restricted to the region downstream of the top of the casting roll. Surface tension is neglected. In the main part of the analysis inviscid fluid flow is considered since the Reynolds number is very large in the applications. It is found that the steady-state solutions are nonunique in a certain parameter range. In addition to a continuous solution, there are two solutions including hydraulic jumps, with one hydraulic jump being located on the casting roll, the other one on the cooling table. Regarding the stability of the non unique solutions, the evolution of disturbances is investigated numerically as an initial-value problem. It is concluded that the hydraulic jump on the cooling table is unstable, while the other discontinuous solution as well as the continuous solution are stable for sufficiently small disturbances. Which stable solution is attained in the steady state, depends on the history of the process. Friction at the liquid/solid interface is taken into account in the last part of the analysis. A constant friction coefficient is assumed. It is found that the history of the process determines the steady-state solution if, and only if, the friction coefficient is sufficiently small. For larger values of the friction coefficient, the steady-state solution is unique and independent of the history of the transient process. Furthermore, for sufficiently large friction coefficients, stable hydraulic jumps are found, in contrast to the inviscid case, also on the cooling table. Received 19 March 1999; accepted for publication 3 May 1999     

A kernel gradient free (KGF) SPH method

The finite particle method (FPM) is a modified SPH method with high order accuracy while retaining the advantages of SPH in modeling problems with free surfaces, moving interfaces, and large deformations. In both SPH and FPM, kernel gradient is necessary in kernel and particle approximation of a field function and its derivatives. In this paper, a new FPM is presented, which only involves kernel function itself in kernel and particle approximation. The kernel gradient is not necessary in the whole computation, and this approach is thus referred to as a kernel gradient free (KGF) SPH method. This is helpful when a kernel function is not differentiable or the resultant kernel gradient is not sufficiently smooth, and thus it is more general in selecting a kernel function. Moreover, different from the original FPM with an asymmetric corrective matrix, in the new FPM, the resultant corrective matrix is symmetric, and this is advantageous in particle approximations. A series of numerical examples have been conducted to show the efficiencies of KGF‐SPH including one‐dimensional mathematical tests of polynomial functions with equal or variable smoothing length and two‐dimensional incompressible fluid flow of shear cavity. It is found that KGF‐SPH is comparable with FPM in accuracy and is flexible as SPH. Copyright © 2015 John Wiley & Sons, Ltd.     

曲靖水城水库坝基膨胀土特性研究

在膨胀土地区构筑建筑物一直是岩土工程界探讨的热点, 尤其是在膨胀土地区修筑水库大坝更是水利工程中的难点, 为此, 本文结合水城水库实际工程, 对膨胀土的特性进行了系统试验研究。水城水库库区为一膨胀土地区, 拦河大坝为一心墙堆堆石坝, 最大坝高 5 9.2m。本项工作系统地进行了膨胀土的物理 -力学性质试验及分析, 尤其是膨胀率与干容重和含水量及外覆压力关系分析。得出了该地区膨胀土的膨胀性不仅与其自身矿物成分有着重要关系, 同时其是否表现为膨胀还取决于天然干容重的大小、含水量及外覆压力的大小, 并且得出了它们之间的相互关系规律。     

新疆下坂地水库坝基防渗墙试验研究

新疆下坂地水利枢纽工程地处高寒、高海拔的山区,坝基覆盖层厚度达147.95m,透水性强,易坍塌和漏浆,块石坚硬,坝基防渗处理难度罕见。针对这一工程问题进行了坝基防渗墙试验研究,结果表明,以大深度混凝土防渗墙为主墙,下接帷幕灌浆的墙幕结合方案。本文针对新疆下坂地深厚覆盖层坝基防渗墙进行了试验研究,找出了防渗墙有效的施工深度及合理的施工方法,从而为坝基防渗墙设计施工提供了可靠的依据。     

顾及高程差异的网络差分对流层误差内插模型

在GPS网络差分定位技术中,用户站与参考站间的高程差异会引起对流层误差改正数中存在系统性偏差,造成对流层改正精度降低,影响用户定位的效果。在常规两维线性组合模型的基础上,考虑高程因素影响,利用三维位置信息对内插模型参数进行约束,将对流程延迟改正数分为距离相关和高程相关两部分,通过卡尔曼滤波器分别估计,提出了一种能修正高程偏差的对流层三维线性组合模型。以重庆连续运行参考站系统中5个站作为研究对象,实验验证结果表明,在地形起伏较大的地区,三维线性组合模型可将对流层改正精度从20～40 cm提高到5～6 cm,相对于常规方法提高了4～6倍,消除了网络差分定位中高程差引起的对流层改正系统偏差。     

弹性半空间二维出平面自由波场的一维化时域算法

提出了一种计算出平面SH波斜入射时弹性半空间自由波场时域计算的一维化有限元方法。首先利用Snell定律确定平面波沿水平方向的传播规律,在用有限元法对弹性半空间进行离散化时,竖向单元尺寸根据波动有限元模拟精度要求确定,而水平向有限元网格尺寸根据水平向波的传播规律和采用的离散时间步长确定,使得有限元离散模型中任意节点的运动可以用水平向相邻节点的运动表示,从而将二维有限元节点运动方程组化为一维的形式。求解此一维方程组,可得到弹性半空间中一列节点的运动,再根据行波的传播规律,可确定全空间自由波场。理论分析和数值算例表明,该方法具有较高的精度和良好的稳定性。     

风机叶片涂层雨蚀研究

研究了风电叶片在雨水冲击下的失效过程,利用有限元方法建立雨水冲击的数值模型,分析了雨滴的冲击角度、直径、涂层性能等对涂层受力的影响,并对双雨滴及四雨滴耦合模型进行了研究。结果表明:风雨耦合流场中叶片最大载荷在叶尖处,涂层的冲击受力与雨滴冲击角度和雨滴直径的大小呈正相关。低模量涂层能大幅降低涂层受冲击时的拉伸应力,小幅增加其压缩应力。同一平面内的双雨滴耦合时,拉伸应力受影响的范围较大,压缩应力受影响的范围较小。不同高度差下的双雨滴耦合中,随着高度差的增大,材料拉伸和压缩应力都呈现先增后减的趋势,其中压缩应力随高度差的变化幅度明显大于拉伸应力。四雨滴耦合时压缩应力对间距的变化更加敏感。涂层失效的过程主要为:涂层表面拉应力疲劳产生微裂纹;液压渗透引起裂纹扩展,涂层质量减少,表面粗糙度增加;表面小范围破坏后形成的凹面结构促进涂层的失效。     

岩溶塌陷的类型、成因机制及防治途径──兼论武汉市岩溶塌陷勘查工作要点

按岩溶塌陷的类型, 对其成因机制作了分析。在此基础上, 提出了武汉市岩溶塌陷勘查与防治工作要点。     

Performance of a flapping foil flow energy harvester in shear flows

Through numerical simulations, we investigate the energy harvesting performance of a heaving/pitching foil in shear flow. With two-dimensional Navier–Stokes simulations, we examined the energy harvesting efficiencies of such a system in linear shear flows and compared the results with those in uniform flows. It is found that in low shear rates, the performance of the system in linear shear flow is slightly higher than that in uniform flow, whereas the energy harvesting efficiency is greatly diminished if the shear rate is sufficiently high (this effect is more pronounced in higher frequencies). This is attributed to the effects of linear shear on the vorticity generation and the synchronization between fluid forcing and foil motion – when a strong shear flow is introduced the lift force induced by the leading edge vortex that is in phase with the heaving motion of the foil is diminished. Furthermore, by studying the instability of the wake behind the foil, we confirm that the optimal performance of the foil in linear shear flows is associated with the same physical mechanism that controls the performance of the foil in uniform flows, i.e. the excitation of the most unstable modes in the wake when the oscillation frequency of the foil is close to the frequencies of these modes.     

多体系统接触碰撞问题的牛顿-欧拉线性互补方法

基于非光滑动力学方法的多体系统接触碰撞分析是目前多体系统动力学的研究热点.本文采用牛顿-欧拉方法建立多体系统接触、碰撞问题的动力学模型,给出一种牛顿-欧拉型线性互补公式.该建模方法与目前一般采用的拉格朗日建模方法的不同之处是约束条件中除了库仑摩擦、单边约束之外还含有光滑等式约束.在建立系统动力学模型时,首先解除摩擦约束和单边约束得到原系统对应的基本系统.牛顿-欧拉方法采用最大数目坐标建立基本系统的动力学方程,由于坐标不相互独立,因此基本系统中带有等式约束,其数学模型为一组微分代数方程.借助约束雅可比矩阵,在基本系统微分代数方程中添加摩擦接触和单边约束对应的拉氏乘子,就可以得到系统全局运动的具有变拓扑结构特征的动力学方程,再结合非光滑约束互补条件便可构成完备的系统动力学模型.完备的动力学模型由动力学微分方程以及等式约束和不等式约束组成.线性互补公式采用分块矩阵形式进行推导,简化了推导过程.数值计算采用基于线性互补的时间步进算法.时间步进算法是目前流行的非光滑数值算法,其突出特点是可以免去数值积分中繁琐的事件检测过程,而数值积分过程中通过对线性互补问题的求解可以确定系统的触-离状态.通过对典型的曲柄滑块间隙机构进行数值分析,验证本文方法的有效性.     

On Bending–Tension Coupling of Laminates

This paper deals with the evaluation of coupling between the in and out of plane behaviours of a laminate composed of identical anisotropic layers. The subject has been motivated by problems arising from coupling in modern practical industrial applications. Coupling is an effect which is often undesired by laminate designers, and some rules and solutions, recalled in the paper, are known today to counter it. Two subjects are analysed in the article: firstly, it is considered the problem of finding the stacking sequence angles which determine the maximum coupling of a laminate whose number of plies and elementary layer properties are known. Secondly, the sensitivity of uncoupling to an orientation defect is analysed, that is, how coupling changes when an error in the orientation of one layer is made: a simple and exact formula assessing the degree of coupling, which is the ratio of the real coupling with respect to the possible maximum, is given. The polar method, proposed by G. Verchery, has been used to develop the research, which is made in the theoretical framework of the Classical Laminated Plate Theory. This revised version was published online in July 2006 with corrections to the Cover Date.     

Marciniak–Kuczynski type modelling of the effect of Through-Thickness Shear on the forming limits of sheet metal

The Marciniak–Kuczynski (MK) forming limit model is extended in order to predict localized necking in sheet metal forming operations in which Through-Thickness Shear (TTS), also known as out-of-plane shear, occurs. An example of such a forming operation is Single Point Incremental Forming. The Forming Limit Diagram (FLD) of a purely plastic, isotropic hardening material with von Mises yield locus is discussed, for monotonic deformation paths that include TTS. If TTS is present in the plane containing the critical groove direction in the MK model, it is seen that formability is increased for all in-plane strain modes, except equibiaxial stretching. The increase in formability due to TTS is explained through a detailed study of some selected deformation modes. The underlying mechanism is a change of the stress mode in the groove that results in a delay of the onset of localized necking.     

Numerical calculation of multiple scattering with the layer model

The optical measurement technique based on Mie scattering has been applied to various areas, in which single scattering at low particle concentration is assumed. Nevertheless, since multiple scattering is usually unavoidable in online measurements, we present in this work a multiple scattering calculation method, in which a layer model is employed. The three-dimensional particle system is divided into a pile of layers the number of which is automatically determined, depending on the obscuration of the particle system. The calculation is found to be fast, reasonable and reliable.     

某软土深基坑工程时间效应有限元分析

软土地区的深基坑因土体的固结作用和流变性而具有了时间效应。本文以Biot固结有限元法为基础,用三元件模型中的第一个线性弹簧模拟固结作用,弹性模量考虑了开挖应力路径和应力历史的影响;另外一部分(KELVIN模型)来模拟土体的流变性,以实际变形的反演来得到两个参数的大致取值,再对基坑的变形情况以及进一步开挖进行分析。假定为正常固结饱和粘土,平面应变问题。通过对某饱和软粘土地基深基坑开挖工程实例的分析,得到的挡墙水平位移曲线与实测曲线很吻合,表明程序较好地反映出基坑的时间效应。     

近圆参考轨道卫星编队洛仑兹力控制

基于运动电荷在磁场中切割磁力线受到洛仑兹力作用的物理规律,分析了两种带电模式对经典轨道根数长期变化的影响:(1)卫星恒定带电模式;(2)前半个轨道周期卫星带电、后半周期不带电的非恒定带电模式.恒定带电模式可以有效地改变轨道升交点赤经、近地点幅角以及平近点角,对轨道半长轴、偏心率和倾角几乎不产生长期影响;而非恒定带电模式则可以有效地改变轨道偏心率.基于洛仑兹力作用下轨道根数长期变化规律以及轨道根数差描述的带电副星相对于不带电主星的运动,提出了利用洛仑兹力以及两种带电模式实现地球低轨近圆参考轨道卫星编队的控制策略,包括编队绕飞椭圆大小重构与编队中心漂移控制,解析求解了副星所需的带电量,并利用数值仿真验证了洛仑兹力控制的可行性.需要指出的是,洛仑兹力轨道控制无需消耗推进工质.      

张弦梁(桁架)结构的找形方法——分级卸载法

对张弦梁(桁架)结构的找形问题进行了深入研究,根据该类结构的受力特点,提出了分级卸载法。该方法无需假定结构零状态几何,而可以一次得到对应于初状态(即预应力态)几何的结构施工放样尺寸(即零状态几何),并通过分级加载对该方法的精度进行了检验。     

Comparison between three-dimensional linear and nonlinear tsunami generation models

The modeling of tsunami generation is an essential phase in understanding tsunamis. For tsunamis generated by underwater earthquakes, it involves the modeling of the sea bottom motion as well as the resulting motion of the water above. A comparison between various models for three-dimensional water motion, ranging from linear theory to fully nonlinear theory, is performed. It is found that for most events the linear theory is sufficient. However, in some cases, more-sophisticated theories are needed. Moreover, it is shown that the passive approach in which the seafloor deformation is simply translated to the ocean surface is not always equivalent to the active approach in which the bottom motion is taken into account, even if the deformation is supposed to be instantaneous.