冲压板材拉伸筋阻力的一种有效数值计算方法

将弹塑性有限变形的拟流动角点本构理论和厚向各向异性屈服函数引入弹塑性动力显式有限元列式,对板材通过拉伸筋的变形过程及拉伸筋阻力进行了数值模拟,并与有关实验结果进行了比较很好的一致性,表明了该计算方法的有效性。进而,数值研究了拉伸筋形状、界面摩擦状况以及材料的厚向各向异性对拉伸筋阻力的影响,为实际覆盖件冲压成形中拉伸筋的设置提供了重要的定量依据。     

超高强度钢AF1410塑性流动特性及其本构关系

在本文中,为揭示超高强度钢AF1410的塑性流动性,并研究其塑性流动本构关系,利用CSS4410电子万能试验机和改进的Hopkinson拉压杆技术,对AF1410钢在温度从100K到600K,应变率从0.001/s到2000/s,塑性应变超过20%的塑性流动特性进行了试验研究。结果表明,拉伸加载下AF1410钢屈服强度低于压缩屈服强度,且随应变率增加,拉压屈服强度差值越来越大;该材料塑性流动应力对应变率敏感性低,而对温度较为敏感;随应变率的提高,该材料拉伸失效应变减小,但温度对失效应变无明显影响。最后基于位错的运动学关系,借助试验数据,获得了AF1410钢的塑性流动物理概念本构模型,并通过与经典J-C模型的结果对比对该物理概念本构模型进行了评估分析,表明该物理概念本构模型在较宽温度和应变率范围能较好的预测AF1410钢的塑性流动应力。     

6061铝合金板材作为汽车顶盖在压力成形过程中的有限元分析

针对6061铝合金车身板材成形性能较差的难点,本文通过试验获得了6061铝合金板材的材料参数,并采用冲压仿真软件 DYNAFORM 以典型覆盖件(某型号轿车的顶盖)为例分析了铝合金顶盖在压力成形过程中的厚度变化、应力和应变分布,说明了6061铝合金板对于汽车顶盖的冲压成形具有良好的适应性。用 Belytschko-Tsay 薄壳单元的 Barlat 三参数模拟结果进一步验证了6061铝合金轿车复杂覆盖件数值模拟的有效性和可行性,其结果可为6061铝合金板材的冲压成形工艺的设计提供依据。     

Q235B钢动态本构及在LS-DYNA中的应用

采用万能材料试验机和分离式霍普金森拉杆（SHTB）装置,对我国钢结构建筑中最常用的Q235B钢进行准静态拉伸实验、高温拉伸实验和动态拉伸实验。基于实验数据对LS-DYNA常用的3种动态材料模型Cowper-Symonds本构模型、Johnson-Cook本构模型、Zerilli-Armstrong本构模型进行了拟合,通过Taylor杆实验对3种本构模型进行验证和对比分析。结果表明:Q235B钢具有较为明显的高温软化和应变率强化效应;Cowper-Symonds本构模型可以较好地适用于工程领域低速碰撞的模拟;Johnson-Cook本构模型可适用于较大应变率范围内的模拟;不推荐Zerilli-Armstrong本构模型在工程低速碰撞领域中使用。     

金属铍的压缩变形行为

利用材料试验机及Hopkinson杆装置系统开展热等静压金属铍在不同温度下的静动态压缩力学行为研究,获得了温度、应变率对金属铍屈服强度和加工硬化行为的影响规律。结果表明:金属铍在压缩应力状态下呈现出良好的塑性,同时其力学性能具有显著的应变率敏感性与热软化效应,屈服强度和流动应力随应变率提高呈明显增大趋势,随着温度升高逐渐降低。同时,室温下其加工硬化行为随着应变增大表现为分段硬化特征,随温度升高则趋于理想塑性。最后,采用修正的Johnson-Cook本构模型对实验结果进行了拟合,模型计算结果与实验结果吻合较好。     

????????????????????????????

基于群体与适应度的概念，应用改进的PSO算法，从随机解出发，提出了基于PSO 算法的本构模型参数识别方法. 该方法解决了橡胶类材料大应变时硬化现象的本构模型参数 的确定这一难题. 首先通过单轴拉伸本构模型实验，在针对硬化实验曲线存在拐点的情 况下，应用PSO算法进行拟合，最后利用简单剪切实验进行验证. 结果表明该方法科学可行， 且具有速度快、精度高、易于收敛等优点. 并且有效地解决了本构模型参数识别的困难，可 广泛应用于各种复杂材料.     

PMMA低、中应变率单向拉伸力学性能的实验研究

利用MTS810试验机和自行研制的中应变率材料试验机对有机玻璃(PMMA)准静态(应变率为2.00×10-2s-1,2.01×10-3s-1和2.38×10-4s-1)和中应变率(应变率为18.6s-1,2.81s-1,6.54×10-1s-1和2.92×10-1s-1)单向拉伸力学性能进行了实验研究。应变率为2.38×10-4s-1时,应力-应变曲线中存在软化段。结果显示,在上述应变率范围内,PMMA的力学性能具有明显的应变率相关性,表现为应变率强化、应变率硬化和高速脆性。PMMA的拉伸力学性能在应变率达到18.6s-1时出现了更大的应变率敏感性。采用包含一个非线性弹簧和六个松弛模式的粘弹性模型分析PMMA力学性能的应变率相关性,得到PMMA在低、中应变率下单向拉伸的本构方程,理论结果与实验结果符合较好。     

高应变率下TiC_P/Ti复合材料的动态拉伸性能

采用旋转盘式杆-杆型动态拉伸试验机对TiCP颗粒增强钛基复合材料及其基体钛合金的动态拉伸性能进行了研究.同时为了比较,在MTS810试验机上做了两种材料的准静态试验.试验结果表明,复合材料及基体材料屈服后至材料的迅速失效,几乎没有应变硬化效应;复合材料的抗拉强度和屈服强度较基体明显提高,但延性明显下降;钛合金基体和复合材料均有明显的应变率强化效应,但复合材料的应变率强化效应明显高于基体;建立了复合材料率相关的本构关系.最后从位错等微观角度分析了复合材料的强化机理、复合材料的应变率敏感性以及复合材料应变率敏感性高于基体的原因.     

用损伤理论方法预测铝合金薄板成型极限

应用各向异性损伤理论研究2024-T3铝合金薄板的成形极限,通过构造有限元单胞模型预测薄板结构的极限应变.单胞模型由两相材料组成:铝合金基体和金属强化物.基体采用全耦合弹塑性-损伤本构方程描述,而金属强化物则视为弹脆性材料.采用所提出的缩颈准则,得到了双轴拉伸状态下铝合金薄板的极限应变,和实验结果比较两者吻合较好.研究结果揭示有限元单胞模型可以提供铝合金的细观损伤机理信息,当忽略材料的损伤影响,采用金属薄板成型理论的研究结果将过高估计薄板的极限应变.     

6061-T6铝合金动态拉伸本构关系及失效行为

采用HMH-206高速材料试验机开展了6061-T6铝合金在0.001～100 s?1应变率范围内的静、动态拉伸力学性能实验,分析了其应力-应变响应特征和应变率敏感性,讨论了应变率对6061-T6铝合金流动应力和应变率敏感性指数的影响,并基于实验结果对Johnson-Cook本构模型进行了修正。结合缺口试件的实验结果和模拟数据,得到了材料的Johnson-Cook失效模型参数,并对模型的准确性和适用性进行了验证。结果表明,在拉伸载荷作用下,6061-T6铝合金表现出明显的应变硬化特征和应变率敏感性,其流动应力随应变率的升高而提高,修正的Johnson-Cook本构模型可以描述材料的动态塑性流动行为,建立的Johnson-Cook失效模型能够表征材料的断裂失效行为。     

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

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

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

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