膜-基复合材料界面剪应力的三维半解析计算

 采用影响系数法对膜-基复合材料的界面剪应力三维半解析进行 了分析研究.利用三维有限元方法对薄膜的影响系数进行了计算. 将 基体作为半无限大体，利用其平面边界作用单位力时的位移场解析 解，得到基体的影响系数. 结果表明，对膜-基复合材料界面剪应力 进行三维半解析计算，克服了完全用三维有限元对其进行计算的限 制，为该类问题的分析提供了新的途径.     

膜-基复合材料界面剪应力分析的影响系数法

提出了一种计算膜 基复合材料界面剪应力的新方法。将薄膜、基体分别作为二维和三维问题进行有限元计算 ,利用影响系数 ,以剪应力为未知量 ,求出同一点处薄膜、基体的位移。通过薄膜、基体在界面上位移相等的条件建立线性方程组 ,联立求解 ,计算出界面剪应力。结果表明 ,本文提出的方法思路清晰 ,易于实现 ,精度令人满意 ,为膜 基复合材料结构强度分析提供了一种新的途径     

初应力对压电压磁板中水平剪切振动的影响

针对含有初始应力的压电压磁板中的水平剪切振动问题,从基本的线弹性理论出发,经过适当的变量代换,推导出了精确的解析解,得到了电磁学开路、电学开路磁学短路、电学短路磁学开路、电磁学短路四种情况下的频率控制方程;进一步应用二分法对方程进行求解,并数值模拟了初始应力以及不同的电磁学边界条件对振动频率、位移模态的影响．研究结果表明：初始应力和电磁学边界条件对压电压磁板的水平剪切振动频率都具有显著的影响,且当初始应力变化范围较小时（初始应力与有效刚度系数之比在5%以内）,振动频率与初始应力之间呈线性关系;电磁学边界条件能显著改变压电压磁板内振动位移的分布,但初始应力对其影响不大．     

磁电弹板双材料板中厚度-扭曲波的传播

研究了磁电弹双材料板(板的上下表面是机械自由、电学开路和磁学短路的)中厚度-扭曲波的传播. 基于磁电弹全耦合的本构方程, 首先推导了厚度扭曲波传播的一般解, 然后利用边界条件和界面条件得到了频率方程. 数值算例分析了压电(BaTiO_{3})-压磁(CoFe_{2}O_{4})双材料板中厚度扭曲波传播时, 位移、电势和磁势的变化特征以及模态数与几何尺寸的关系.      

无黏流体与固体界面的漏瑞利波

漏瑞利波存在于半无限无黏性流体和半无限固体媒质的界面处. 首先推导流固无限各向同性介质界面处漏瑞利波的特征方程和位移及应力的解析计算公式. 然后结合典型结构通过数值计算研究了漏瑞利波特性以及位移和应力在流体和固体中的分布规律. 数值计算结果表明漏瑞利波的相速度和衰减随流固密度比的增大而增大, 在流固界面上法向位移连续而切向位移不连续. 流固密度比对固体媒质中沿垂直于漏瑞利波的传播方向的位移、正应力和剪应力有比较大的影响,而对沿漏瑞利波的传播方向的正应力几乎没影响. 为利用漏瑞利波的无损检测与评价提供了理论基础.      

粘接–映射混合算法及其对颗粒材料中板的振动特性分析

结构与颗粒材料相互作用广泛存在于各工程领域,其研究过程中涉及的连续–离散耦合计算方法面对诸多挑战.本文提出了粘接–映射混合算法来研究连续体与离散介质耦合动力学问题.将连续体模型划分为内部区域及与颗粒接触的边界区域.边界区域采用粘接算法模拟连续体外部形状并使用高效的球形接触判断准则;提出一种包含Rayleigh阻尼映射的有限元映射质点弹簧算法来精确计算连续体内部区域内力和变形.二者相结合构成粘接–映射混合算法,并引入计算机集群和GPU(图形处理器)并行技术,对埋没于颗粒材料中受激振动固支方板的连续–离散耦合动力学问题进行了数值仿真研究.结果表明,粘接–映射混合算法有利于双层级并行算法的程序实现及优化,并在连续–离散耦合界面进行快速接触判断的同时实现对颗粒材料中方板位移、变形、振动形态等参数的研究.通过定幅扫频和定频变幅方式考察激振力频率和幅值对振动板非线性动力学行为的影响并观察到二倍周期现象,同时给出了该连续–离散耦合系统中颗粒体系的能量耗散特性.     

磁驱动中心对称挠曲电夹层板力电耦合性能分析

现代工业的发展对材料性能和结构尺寸提出更高的要求,机电器件的设计越来越偏向于小型化、高频化和智能化.最新研究成果表明,磁电耦合复合材料不仅能够以较强的磁电转换效率实现磁能、机械能和电能之间的相互转换,还可以避免结构与机械驱动源的直接接触,实现非接触调控,这对于制备多功能微纳米器件具有重要意义.文章基于Mindlin所发展的多物理场结构理论分析方法,结合宏观压磁理论和偶应力挠曲电理论,研究由单个挠曲电电介质层和两个对称压磁层构成的三明治型夹层板在外部横向磁场驱动下的动态力电耦合响应,其中通过引入曲率将经典力电耦合理论拓宽到中心对称材料.夹层板在正弦型全局磁场和均布局部磁场驱动下的动态数值算例表明:位移和电势具有一定的频率依赖性,当激振频率达到固有频率时,振幅达到最大值;此外,对称式驱动压磁层分布方式趋于提高多层复合板的力电耦合性能.文章理论模型和研究结果可为磁控机电器件的优化设计提供新的改进思路.     

考虑非局部效应的双层纳米板的磁弹性随机振动

研究了四边简支双层纳米板在外加磁场的作用下的磁弹性随机振动问题．基于非局部弹性理论和板壳 磁弹性理论建立了系统的磁弹性随机振动方程．通过模态分析法对其进行位移响应分析，得到了通入平稳随机 电流时双层纳米板位移响应均值、功率谱密度函数等数字特征．在此基础上，分析了非局部参数、磁场强度、 板厚比等对功率谱密度的影响．结果表明，非局部参数、磁场强度、板厚比等因素的变化会影响系统的振动能 量变化及振动响应带宽分布.     

粘接-映射混合算法及其对颗粒材料中板的振动特性分析

结构与颗粒材料相互作用广泛存在于各工程领域,其研究过程中涉及的连续-离散耦合计算方法面对诸多挑战.本文提出了粘接-映射混合算法来研究连续体与离散介质耦合动力学问题.将连续体模型划分为内部区域及与颗粒接触的边界区域.边界区域采用粘接算法模拟连续体外部形状并使用高效的球形接触判断准则;提出一种包含Rayleigh阻尼映射的有限元映射质点弹簧算法来精确计算连续体内部区域内力和变形.二者相结合构成粘接-映射混合算法,并引入计算机集群和GPU(图形处理器)并行技术,对埋没于颗粒材料中受激振动固支方板的连续-离散耦合动力学问题进行了数值仿真研究.结果表明,粘接-映射混合算法有利于双层级并行算法的程序实现及优化,并在连续-离散耦合界面进行快速接触判断的同时实现对颗粒材料中方板位移、变形、振动形态等参数的研究.通过定幅扫频和定频变幅方式考察激振力频率和幅值对振动板非线性动力学行为的影响并观察到二倍周期现象,同时给出了该连续-离散耦合系统中颗粒体系的能量耗散特性.      

杂交梁界面剪应力的理论分析和数值模拟

采用解析和数值方法研究FRP-混凝土杂交梁的界面应力问题。提出了杂交梁的新的力学模型和假设,克服了以往的分析模型中界面应力表达式非常复杂和界面应力的解析解与数值解相差较大的缺点,本文得到的FRP板加固梁的界面剪应力表达式与数值结果符合很好,并且具有简捷的表达式。利用有限元法研究了杂交梁各物理参数对界面剪应力的影响。研究表明,界面剪应力在FRP板的端部存在应力集中或应力奇性,这是造成杂交梁界面破坏的主要原因。这项研究对进行杂交结构的工程设计具有理论指导和参考价值。     

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

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

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


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.     

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).