用于C60纳米分子振动分析的碳-碳键合单元及坐标变换

C60是一种纳米分子,具有许多非常重要的物理化学性质,采用Raman光谱可以测得它的分子振动特征,正是这种振动特征可以明显地表征出它的分子结构和原子键合关系。本文基于分子力学的基本分析原理,给出一种用于描述碳-碳共价键在小位移情形下力能关系的计算单元,称为碳-碳键合单元,重点讨论该单元在实际建模过程中的坐标变换问题,将所建立的碳-碳键合单元用于C60纳米分子的振动分析,分别给出了C60纳米分子的Ag(1)、Ag(2)、Hg(1)三个特征振型和频率,并与群论的分析结果、Raman光谱实验结果进行比较,证明了所提出分析方法的正确性和实用性。     

碳纳米管的等价柔性节点空间框架模型

小变形情况下,碳纳米管C-C共价键间的相互作用可以用基于分子力学的宏观力学模型进行模拟.其中,基于分子结构力学的等价结构力学模型是最为有效的碳纳米管弹性参数的预测模型.现有的碳纳米管的等价结构力学模型是用具有刚性节点的空间框架结构模拟碳纳米管的原子晶格受力和变形的关系.根据碳纳米管的原子晶格的变形特点,本文首次提出了一个用柔性节点空间框架模拟碳纳米管原子晶格键角变化的分析模型,再通过应变能等价推导了柔性节点的等价抗弯刚度与分子力学中力常数的关系,从而给出了一个更精确的计算碳纳米管等价弹性参数的分子结构力学模型.文中用ANSYS计算了不同尺寸的锯齿型(zigzag)和扶手型(armchair)单壁碳纳米管的轴向杨氏模量、泊松比、剪切模量及径向杨氏模量,分析了碳纳米管的尺寸效应,并且与其它各种模型所得结果进行了比较.计算结果表明,本文所给碳纳米管的等价柔性节点空间框架模型不仅计算简单、高效,而且准确;并可以直接推广到多壁碳纳米管等价弹性模量的计算及碳纳米管的稳定和动力分析.     

纳米级石墨晶体层片力学性能的等效模型及分析

石墨晶体是由碳原子以sp2杂化轨道和邻近的三个碳原子形成共价键,在其平面内再构成网状层片结构。文中对基于原子键合力能关系的CCBE等效模型和经典的BEAM单元等效模型作了比较分析。其中CCBE包括两个节点,每个节点包含三个平移自由度。与经典BEAM单元相比,其两个节点均没有转动自由度。CCBE和BEAM单元的各个参数可以通过石墨晶体的相关实验获得。通过对石墨晶体层片的拉伸、剪切等效模型的柔度(刚度)系数计算分析,并与实验数据相对比,验证了CCBE等效建模方法的正确性和合理性。该等效建模方法为研究石墨的微观结构、性能等提供了一种有效的数值模拟方法。     

金沙江虎跳峡河段岸坡变形破坏的相关动力因子研究

石墨层和单臂碳纳米管都是以C---C共价键结合的. 在小变形条件下C---C键的势能可用谐和函 数来描述，这与梁单元的变形能具有相同的形式，因此可以用梁单元等效C---C键的作用. 提出了一种C---C键的等效梁单元有限元模型，该模型能够完备地替代谐和势描述C---C键的 伸长、面内键角变化、离面键角变化和扭转. 通过分析石墨层的典型受载情况得到了等效梁 单元的参数，以及等效梁单元参数与谐和势参数的关系，并用该模型计算了单臂碳纳米管的 杨氏模量和泊松比，计算结果为相关文献所验证.     

基于扩展有限元的应力强度因子的位移外推法

针对平面裂纹问题,阐述了扩展有限元法的单元位移模式、推导了扩展有限元法的控制方程、介绍了特殊单元的数值积分技术.基于最小二乘法,建立了应力强度因子位移外推法的计算公式.利用MATLAB编写计算程序,对平面裂纹问题用扩展有限元法进行了计算.基于扩展有限元法的计算结果,分别利用位移外推法和相互作用积分法,对平面裂纹的应力强度因子进行了计算.计算结果表明,位移外推法比相互作用积分法能更方便和准确地计算平面裂纹的应力强度因子.     

均匀热荷载作用下层合简支梁的弹性力学解

基于二维热弹性力学理论,研究均匀热荷载作用下层合简支梁的弹性力学解.首先导出均匀温度场中满足控制微分方程和两端简支边界条件的单层梁的弹性力学解,然后利用层间界面位移和应力必须连续的条件,递推得到底层梁与顶层梁间的位移和应力关系.最后根据层合梁上下表面的边界条件确定待定系数,带回递推公式得到整个层合梁的应力和位移分布.本文方法的计算结果有很好的收敛性.与有限元软件的结果对照说明了本文方法的精确性.最后,研究了不同的变温对层合简支梁的位移和应力的影响,结果显示每个层间界面在x方向的应力是不连续.随着温度的升高,梁的最大位移相应地增大.温度越高,位移沿厚度变化的速率越大.     

岩石试件端面摩擦效应数值模拟研究

试件端面摩擦效应直接影响试件内的塑性等效应变、侧向位移的分布和单元应力应变曲线。本文运用ANSYS中的接触单元模拟了平面应变状态下端面摩擦效应对塑性等效应变、侧向位移和单元应力应变曲线的影响,得到了不同摩擦系数时塑性等效应变及侧向位移的渐进变化形式。当接触面摩擦较小时,塑性等效应变图案为上下两个X形网络,侧向位移上下分布均匀;当接触面摩擦增大时,塑性等效应变网络向中部靠拢并且明显增大,侧向位移上下分布不均匀,中部较上下端面位移大;当试件端面侧向位移被限制,即摩擦力很大时,塑性等效应变网络变为一个X形局部化带,侧向位移分布更加不均匀,中部明显隆起。     

一种新的弹塑性杂交应力元法

本文基于一个改进的弹塑性的Hellinger/Reis■ner 混合变分原理构造了一种用于解弹塑性问题的四节点等参杂交应力元.新的模型中,在单元内增加了等效应力增量、塑性等效应变增量及不协调位移变量,从而使单元内的屈服准则及流动法则平均得到满足,不协调位移改进了单元应力精度.计算表明,新的模型可以提高弹塑性杂交法的精度和计算效率.     

带转动自由度的内参型非协调元研究

非协调元虽然破坏了单元间位移的连续性，却能很好地反映弯曲类变形，然而在不增加单元结点自由度的情况下，非协调元的计算精度总是滞留在某一水平，无法得到较大改变。基于修正后的位移型Reissner泛函中引入独立转动场的变分原理，采用连续介质力学中的转动自由度的定义，转动场采用结点真实转角来插值，结合平面四结点单元讨论了有效附加非协调位移的合理形式，引入了适用于任何四边形单元的非协调位移函数，从而建立了一种带转动自由度的平面四结点内参型非协调元模型。本文单元能通过分片检验，并易于与带转动自由度的梁单元相容．教值算例表明具有较高的计算精度。     

纳米硬度实验中的多位错生成分析

本文利用边界单元法及基于Peierls-Nabarro模型的位错理论，分析了理想纳米触头下多位错的生成，得到了滑移面上多位错的构形以及表面位移与载荷的关系曲线。根据得到的计算结果分析了薄膜厚度的影响，得到了与已有实验结果相同的结论。通过把表面看成一个包含在无穷大弹性体中的无穷大弹性体中无穷大裂纹来考虑表面的影响，从而可以直接采用已有的边界积分方程。该方法在连续介质力学的计算中引入了包含原子信息的层间势能函数，为分析多尺度的力学问题提供了有效的方法。     

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

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

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

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