粘结单元在模拟FRP层合板低速冲击响应中的应用

纤维增强树脂基复合材料层合板（fibre reinforced plastic composites，FRP）在航空、航天、交通、造船等诸多工程中得到了日益广泛的应用，而其在冲击载荷下的响应和破坏特别是分层一直为学术界所关注。本文中对FRP层合板在冲击载荷下的响应和破坏进行数值模拟，并通过引入粘结层重点研究其分层破坏。首先，介绍一种基于改进的粘结区域方法的粘结层损伤模型；其次，详细介绍了有限元模型建模过程和建模细节；最后，对有限元模型进行验证，并分析分层损伤发生的原因。模拟结果表明，该模型不仅能准确预测FRP层合板在低速冲击载荷下的载荷-时间曲线和载荷-位移曲线，还能成功地预测其分层破坏。     

弱围压与垂直冲击强扰动下块系岩体超低摩擦效应试验研究

随开采深度增加,深部岩体结构趋于破碎,临近工作面待开采岩体受弱围压作用,遇冲击载荷强扰动时,易诱发超低摩擦型冲击地压。采用自主研制的弱围压块系岩体超低摩擦模拟试验装置,以5个上下叠加花岗岩块体模拟深部临近工作面破碎岩体,通过螺栓加载方式对块系岩体施加弱围压作用、垂直方向施加冲击载荷作用模拟强扰动。以垂直加速度差值变化作为超低摩擦效应强度特征参数,分析得到了冲击载荷强度、围压对块系岩体超低摩擦效应影响规律。研究结果表明:垂直冲击载荷作用下,块系岩体振动分为受迫振动和自由振动两个阶段,随冲击载荷增大,块体垂直加速度差值峰值增大,垂直加速度差值衰减时间延长,此时超低摩擦效应强度和持续时间均显著增加;围压与冲击载荷共同作用下,存在围压临界值使得超低摩擦效应强度达到最大;随围压增加,超低摩擦效应最大强度发生时刻提前,当试验围压达到最大值时,超低摩擦效应最大强度发生时刻趋于一致。     

延性材料的损伤断裂

以细观统计损伤理论为基础 ,对材料受冲击载荷作用下初期损伤的成核成长效应进行了定量计算。得到了不同动载荷作用下从损伤成核到成长阶段转换的特征时间 ,并把它应用于Johnson提出的延性材料的损伤断裂模型 ,减小了原模型中的模型参数初始损伤度的随意性 ,增加了Johnson模型的确定性与可预测性。用改进后的模型计算不同冲击速度下OFHC铜的损伤断裂 ,计算结果与实测结果吻合较好。     

隧道围岩强度不均地段塌方成因及其处理方法

采用微面模型理论和损伤力学方法,建立了节理岩体的弹塑性损伤耦合微面模型. 在节理岩体的微 面上,将岩体视为由节理面与岩石组成的二元介质,以节理连通率作为岩体沿该方向的面积 损伤变量,考虑微面法向拉应力和压应力下的不同塑性变形和损伤耦合作用机制,基于塑性 理论建立了节理岩体的微面塑性损伤增量本构关系. 采用微面物理量与宏观物理量的几何约 束模型,根据微面方向积分导出了节理岩体的宏观弹塑性增量本构关系. 编制了节理岩体微 面模型的MARC有限元子程序,对节理岩体的单轴拉伸、压缩试验和泥浆压力作用下的井壁稳 定问题进行了数值模拟研究. 数值计算结果表明,该模型能很好地揭示载荷作用下节理岩体 的各向异性非弹性变形和次生节理演化过程.     

深部节理岩体塑性损伤耦合微面模型

采用微面模型理论和损伤力学方法,建立了节理岩体的弹塑性损伤耦合微面模型. 在节理岩体的微面上,将岩体视为由节理面与岩石组成的二元介质,以节理连通率作为岩体沿该方向的面积损伤变量,考虑微面法向拉应力和压应力下的不同塑性变形和损伤耦合作用机制,基于塑性理论建立了节理岩体的微面塑性损伤增量本构关系. 采用微面物理量与宏观物理量的几何约束模型,根据微面方向积分导出了节理岩体的宏观弹塑性增量本构关系. 编制了节理岩体微面模型的MARC有限元子程序,对节理岩体的单轴拉伸、压缩试验和泥浆压力作用下的井壁稳定问题进行了数值模拟研究. 数值计算结果表明,该模型能很好地揭示载荷作用下节理岩体的各向异性非弹性变形和次生节理演化过程.      

无限冰介质中爆炸裂纹扩展机理与数值模拟

根据冰体力学性能,提出了无限淡水冰介质在爆炸冲击载荷下的压碎区及裂隙区的理论计算方法.利用大型有限元软件ANSYS/LS-DYNA模拟了爆炸载荷下无限区域冰体内裂纹扩展过程,模拟结果与实验结果基本一致,表明采用的损伤断裂模型能够较准确的反映无限冰介质动态力学性能及裂纹成形过程.数值模型可近似计算冰体中裂隙区的半径.数值计算值与理论计算值误差不超过5%,数值模拟结果与实验结果吻合,从而确定了无限区域冰体在-15℃时的基本破坏特征及范围.     

冲击载荷作用下岩体破坏规律的数值流形方法模拟研究

为了更好地模拟岩体这种典型的连续与非连续介质共存的复杂结构体在冲击载荷作用下的破坏规律,在近年来新出现的数值流形方法基本理论的基础上,利用其算法程序对岩体中存在的不连续面如节理、裂隙等对岩体破裂效果的影响进行了数值模拟分析,计算结果充分显示了数值流形方法在模拟这类问题中的准确性和有效性。     

多种冲击载荷条件下的人体肋骨骨折有限元分析

旨在研究不同碰撞载荷条件下基于不同失效模型的人体肋骨骨折机理. 为此采用已验证的人体有限元胸部模型来分析人体肋骨骨折现象. 该模型基于人体解剖学结构,包含了人体胸椎、腰椎、肋骨、胸骨、肋间软骨、胸腹部器官和其他的软组织,定义的生物材料参数都基于已有的文献记载. 基于人体在损伤生物力学领域内一些较为典型的肋骨骨折失效模型,根据已发表文献中的人体标本实验载荷条件模拟了人体肋骨结构在不同冲击载荷下的骨折现象,并与这些实验结果进行了对比分析. 所引用的实验结果包括单根肋骨强度结构实验和人体胸部正面碰撞块冲击实验. 从文中有限元仿真分析的结果来看,针对不同的载荷条件,不同肋骨骨折失效模型的适用性各不相同. 该人体胸部有限元模型可用于车辆交通事故中冲击载荷条件下的人体肋骨损伤生物力学研究.      

碳纳米管/碳纤维增强复合材料层合板低速冲击响应和破坏的数值模拟

碳纳米管/碳纤维增强复合材料（carbon nanotube/carbon fibre reinforced plastic,CNT/CFRP）是一种多尺度复合材料,比传统CFRP有更好的综合性能和更广阔的应用前景。对CNT/CFRP在低速冲击下的响应和破坏进行了数值模拟研究。首先,基于先前的研究通过引入基体增韧因子、残余强度因子并改进损伤耦合方程,建立了新的FRP动态渐进损伤模型;然后,利用新建立的本构模型并结合黏结层损伤模型,对4种碳纳米管含量的增韧碳纤维增强树脂基复合材料层合板在5个能量下的冲击实验进行了数值模拟;最后,将模拟结果与文献中的相关实验结果进行了比较,并讨论了冲击速度的影响。结果表明:新建立的FRP本构模型能够预测CNT/CFRP层合板在低速冲击载荷作用下的响应、破坏过程和分层形貌,模拟得到的载荷-位移曲线和破坏形貌与实验吻合较好;冲击速度会影响CNT/CFRP层合板拉伸和压缩破坏的比例,相同的冲击能量下,更大的冲击速度会造成更多的拉伸破坏。     

水中爆炸冲击波载荷作用下舰船结构动态响应的数值模拟

采用大型通用有限元软件MSC.DYTRAN,对某型水面舰船全船结构在水下爆炸冲击波载荷作用下的动态响应进行了数值模拟。分别从结构的变形损伤形式、能量吸收和冲击环境等几个方面研究了舰船结构在水下爆炸载荷作用下的破坏机理和响应特征。     

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

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

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

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

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

正~~     

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