基于凹坑深度的复合材料低速冲击损伤分析

对复合材料层合板进行了低速冲击实验,建立了冲击凹坑深度与冲击能量的关系。依据凹坑深度反推冲击能量,并用能量确定冲锤质量和冲击速度,从而可对层合板进行动态数值模拟。冲击凹坑深度与冲击能量的关系表明,凹坑深度的变化是与冲击能量的变化过程相适应的,在此基础上对损伤的分布形式及大小做了详细的分析。采用ANSYS有限元程序对复合材料层合板横向低速冲击进行模拟,采用瞬态分析方法来研究层合板的冲击与损伤过程。对冲击后的试验件进行了C扫描损失检测。计算和试验结果表明,此方法是可行的,特别适合于层合板冲击后的损伤评估。     

低速冲击激励下嵌入黏弹性阻尼芯层的纤维金属混杂层合板动态响应预测模型

本文首次从解析角度建立了低速冲击激励下嵌入黏弹性阻尼芯层的纤维金属混杂层合板动态响应预测模型. 首先,结合经典层合板理论和冯$\cdot$卡门假设,建立了嵌入黏弹性芯层的纤维金属混杂层合板弹性损伤本构关系. 然后,将层合板受冲击时的变形分成接触和拉伸两个区域,在接触区域内,对金属层采用 Von Mises 失效准则,纤维层采用 Tsai-Hill 失效准则和对黏弹性层采用指数 Drucker-Prager 失效准则判断层合板损伤情况. 考虑不同材料层对冲击动态响应的贡献来修正两个变形区域的位移公式,进而计算结构因弹性变形产生的应变能,以及接触区域因塑性变形消耗的能量,实现每次失效事件发生后各层材料的能量、位移和冲击接触力的理论求解,并给出了结构动态响应分析的具体流程图. 最后,以嵌入 Zn33 黏弹性芯层的 TA2 钛合金混杂 T300 碳纤维/树脂层合板为研究对象,开展落锤冲击实验. 验证结果表明,理论预测与测试获得的冲击接触力、位移响应以及冲击载荷-位移曲线吻合较好,且关注的峰值点计算误差最大不超过 9%,进而验证了所提出的理论模型的有效性.      

A DYNAMIC RESPONSE PREDICTION MODEL OF FIBER-METAL HYBRID LAMINATED PLATES EMBEDDED WITH VISCOELASTIC DAMPING CORE UNDER LOW-VELOCITY IMPACT EXCITATION 1)

本文首次从解析角度建立了低速冲击激励下嵌入黏弹性阻尼芯层的纤维金属混杂层合板动态响应预测模型. 首先,结合经典层合板理论和冯$\cdot$卡门假设,建立了嵌入黏弹性芯层的纤维金属混杂层合板弹性损伤本构关系. 然后,将层合板受冲击时的变形分成接触和拉伸两个区域,在接触区域内,对金属层采用 Von Mises 失效准则,纤维层采用 Tsai-Hill 失效准则和对黏弹性层采用指数 Drucker-Prager 失效准则判断层合板损伤情况. 考虑不同材料层对冲击动态响应的贡献来修正两个变形区域的位移公式,进而计算结构因弹性变形产生的应变能,以及接触区域因塑性变形消耗的能量,实现每次失效事件发生后各层材料的能量、位移和冲击接触力的理论求解,并给出了结构动态响应分析的具体流程图. 最后,以嵌入 Zn33 黏弹性芯层的 TA2 钛合金混杂 T300 碳纤维/树脂层合板为研究对象,开展落锤冲击实验. 验证结果表明,理论预测与测试获得的冲击接触力、位移响应以及冲击载荷-位移曲线吻合较好,且关注的峰值点计算误差最大不超过 9%,进而验证了所提出的理论模型的有效性.     

基于连续丝束剪切技术的变角度复合材料层合板的热屈曲分析

连续丝束剪切(Continuous Tow Shearing, CTS)铺放技术是一种新的变角度层合板制作技术，这种新技术能显著减少丝束铺放过程中产生的丝束重叠和间隙等缺陷，然而，采用CTS技术铺设时，层合板的厚度将随着纤维角度的变化而变化．本文基于一阶剪切变形理论并应用Chebyshev-Ritz法对这种变厚度的变角度复合材料层合板的热屈曲问题进行了研究．假设纤维方向角沿板的长度方向按照线性变化，获得了均匀温度变化时变厚度层合板的临界热屈曲荷载．通过与现有文献的比较验证了本文方法的正确性，并进一步讨论了纤维铺设技术、纤维方向角的变化以及边界条件的不同对变角度复合材料层合板的临界屈曲温度的影响．研究结果表明，在体积相同的情况下，采用CTS铺设较传统的自动丝束铺放(AFP)可以进一步提升变角度层合板的临界屈曲温度．本文的研究结果可为变角度复合材料的设计提供一定的参考．     

灵芝孢子孢壁的弹性模量和硬度研究

采用纳米压痕技术和数值模拟研究灵芝孢子孢壁的弹性模量和硬度.利用原位纳米力学测试与分析系统,测试灵芝孢子孢壁的弹性模量和硬度.得到了载荷--位移曲线图和硬度、弹性模量随压痕深度变化的值.并用有限元方法模拟压痕过程,利用ANSYS软件,按照灵芝孢子孢壁和Berkovich压头的结构,建立了二维计算模型,得到纳米压痕的等效应力分布以及压痕过程中加载和卸载时的载荷--位移曲线.考察了摩擦、压头尖端半径对模拟结果的影响.结果显示:灵芝孢子孢壁的平均弹性模量为2.0GPa,硬度为0.13GPa.模拟结果在趋势上与实验结果有较好的吻合,与理论分析的载荷--位移关系基本一致.摩擦、压头尖端半径小于100nm时对模拟结果不会造成明显影响.研究结果为分析孢子的破壁机理提供必要参数.     

双模量面板泡沫铝芯圆形层合板的非线性弯曲

双模量材料是典型的拉压弹性模量不同的材料,在均匀外载荷作用下,双模量面板泡沫铝芯圆形层合板相当于三种不同材料组成的层合板。采用弹性理论建立了双模量面板泡沫铝芯圆形层合板在均布载荷作用下的静力平衡方程,利用该静力平衡方程确定了层合板的中性面位置。在此基础上建立了双模量面板泡沫铝芯圆形层合板的大挠度弯曲微分方程组,求得了层合板中心挠度与均布载荷的关系式。该方法计算结果与有限元计算结果的最大误差仅为3.8%,这说明该方法是可靠的。算例分析表明不考虑面板拉压弹性模量相异时其计算结果与实际情况相差较大,超过了工程上所允许的计算误差5%。所以,在计算双模量面板泡沫铝芯圆形层合板的非线性弯曲时,不宜采用相同弹性模量弹性理论,而应该采用拉压弹性模量不同的弹性理论。     

具损伤压电智能层合板的动力分析

基于应变能等效原理、高阶剪切变形理论和Hamilton变分原理,考虑复合材料铺设层内的损伤效应,建立了具损伤压电智能层合板的运动控制方程,并运用Galerkin方法进行求解.数值算例中,讨论了损伤效应、厚跨比及压电层厚度与层合板总厚度之比对四边简支压电智能层合板自由振动频率的影响和外部控制电压对其动力响应的影响.     

用能量法分析层状材料弹性接触问题

利用能量法分析了层状材料（薄膜／基体）弹性接触问题，得到了具有一阶精度的闭合解，给出了求解薄膜弹性模量和泊松比的表达式，并与有限元的数值解进行了比较。二者比较结果表明：在工程材料范围内，理论解与数值解相差在6％以内；同时表明单相材料中剪切模量与弹性模量之间的关系也适用层状材料中的薄膜材料。在数值解的基础上，讨论了薄膜厚度与压头半径的比值对求解精度的影响，发现此比值对精度影响不大。通过对层状材料等效泊松比与等效弹性模量的定义，给出了用压痕实验测定薄膜泊松比与弹性模量的方法。     

基于Legendre多项式的双压电层合板三维压电谱元法模拟

谱单元作为一种高阶单元具有计算效率高和精度高的特点。本文在基于Legendre正交多项式的三维谱单元基础上提出了三维压电谱单元模型,用于压电层合板静力和动力性能模拟研究。压电层合板结构中的位移和电势自由度均离散到三维压电谱单元Gauss-Lobatto-Legendre(GLL)配置节点上,并且未对电势沿压电层厚度方向上的变化做任何假定。通过提高谱单元中沿压电层合板厚度方向上的形函数阶数的方法,来削弱三维谱单元在模拟薄板结构中出现的剪切闭锁现象。为验证单元的计算精度,取双压电层合板结构进行静力和动力行为模拟验证,并将计算结果与现有文献中的其它单元模型及有限元结果进行对比。结果表明,压电谱单元可有效模拟压电层合板的静力和动力行为,且提高谱单元形函数阶数可提高数值模拟精度。     

考虑剪切效应层合板的Hamilton体系及辛几何方法

基于ＹＮＳ层合板理论,通过对混合能变分原理的修正,建立了层合板问题的Ｈａｍｉｌｔｏｎ正则方程。在辛几何数学框架下,采用共轭辛正交归一关系给出精确解。并与经典层板理论进行了比较。     

A New Formula for Evaluating Indentation Toughness in Ceramics

The Lawn-Evans-Marshall model for Vickers-induced crack driving force is widely applied for determining fracture toughness of ceramic materials, yet its simplification that the volumes of Vickers impression and equivalent expanding cavity merely depend on the half-diagonal length of residual impression will bring significant error. In this work, the model is modified by considering the influence of elastic recovery in depth after completely unloading on the volumes. As a result, a new expression for crack driving force is obtained and correspondingly a new fracture toughness formula is presented based on the modified model, which is expressed as the well-known Anstis formula with a different constant multiplied by the ratio of residual impression depth to the maximum indentation depth. Both finite element analysis and experimental results demonstrate the effectiveness of the proposed formula and, especially, that it works much better than Anstis formula on materials with large ratio of hardness to elastic modulus.     

Analytical solution to continuous contact problem for a functionally graded layer loaded through two dissimilar rigid punches

A continuous contact problem of functionally graded layer resting on an elastic semi-infinite plane, which is loaded with through two different blocks is addressed in this study. The elasticity theory and integral transformation techniques are used in solution of the problem. The problem is reduced to a system of singular integral equations, and solved numerically by the aid of appropriate Gauss–Chebyshev integration formula. It is assumed that the elastic semi-infinite homogeneous plane is isotropic and all surfaces are frictionless and continuous. The shear modulus and the mass density of the FG layer vary exponentially along the thickness direction.     

An analytical solution for a linear viscoelastic layer loaded with a cylindrical punch: Evaluation of the rebound indentation test with application for assessing viability of articular cartilage

A general closed-form solution for the so-called rebound indentation test’ is obtained for a cylindrical flat-ended punch indenting a linear viscoelastic layer lying on a rigid substrate. Under the assumption of time-independent Poisson's ratio, we derive closed-form analytical expressions for the contact force (in a displacement controlled regime) and for the indentation displacement (in a load-controlled regime) and we consider in detail the case of standard viscoelastic solid. Our results indicate that the rebound displacement (in other words the indentation displacement in the load-controlled stage) is independent of the relaxed elastic modulus and Poisson's ratio, and also of the layer's thickness. Our analytical solution can be used for layered samples of arbitrary materials exhibiting viscoelastic properties; however, since the rebound indentation test has been recently suggested for assessing the viability of biomedical materials, we have applied our theoretical framework to the identification of materials parameters from experiments on articular cartilage. In this context, we have found a pretty good agreement for the rebound deformation, even until the strain becomes relatively large.     

Asymptotic modeling of the impact of a spherical indenter on an elastic half-space

The impact of a rigid sphere on a homogeneous, isotropic elastic half-space in the absence of friction and adhesion is considered. The influence of the superseismic stage immediately following the moment of first contact upon the impact process is investigated in the frame of the Hertzian impact theory. The first order asymptotic approximation for the contact force in a three-dimensional dynamic contact problem with the slowly moving contact zone boundary is obtained and the corresponding asymptotic model of impact is developed. The motion of the indenter as it indents and rebounds from the elastic medium is analytically described. Explicit formulas are derived for the peak indentation depth, contact time, and rebound velocity as functions of the initial impact velocity, indenter mass, and characteristics of the elastic half-space.     

New Double Indentation Technique for Measurement of the Elasticity Modulus of Thin Objects

In this paper we introduce a new method to determine the Young’s modulus of thin (biological) samples. The method is especially suitable for small objects with a thickness of a few hundred micrometers. Such specimens cannot be examined with existing tests: compression and tensile tests need well-known geometry and boundary conditions while classic indentation tests need relatively thick pieces of material. In order to determine the elastic modulus we use the indentation theory as proposed by Sneddon and correct it with a finite element calculated κ factor to compensate for the small thickness. In order to avoid material deformations at the contact zone between the sample bottom and the sample stage, we replace the sample stage by a second indentation needle. In this way the sample can be clamped between two identical needles and a virtual mirror plane is introduced. The new method was used on four test-materials and results agreed well with the outcome of a standard compression method applied on large samples of the same materials. As an application example the technique was applied on thin biological samples, namely middle ear ossicles of rabbits.     

Cylindrical nano-indentation on metal film/elastic substrate system with discrete dislocation plasticity analysis: A simple model for nano-indentation size effect

The cylindrical nano-indentation on metal film/elastic substrate is computationally studied using two-dimensional discrete dislocation plasticity combined with the commercial software ANSYS®, with a focus on the storage volume for geometrically necessary dislocations (GNDs) inside the films and the nano-indentation size effect (NISE). Our calculations show that almost all GNDs are stored in a rectangular area determined by the film thickness and the actual contact width. The variations of indentation contact width with indentation depth for various film thicknesses and indenter radii are fitted by an exponential relation, and then the GND density underneath the indenter is estimated. Based on the Taylor dislocation model and Tabor formula, a simple model for the dependence of the nano-indentation hardness of the film/substrate system on the indentation depth, the indenter radius and the film thickness is established, showing a good agreement with the present numerical results.     

具有弹性层的两圆管间夹入弹性薄膜时的回转接触

应用弹性接触理论及Fourier变换,考虑接触域的粘着及微小滑状态,对具有弹性层的两圆筒间夹入弹性薄膜的非对称回转接触问题进行了0研究,讨论了接触宽度,弹性层厚度及摩擦系数等对弹性薄膜传送速度和圆筒公称圆周速度之比的影响。     

Smooth static and dynamic indentation of a cantilever beam

The static and dynamic indentation of structural elements such as beams and plates continue to be intriguing problems, especially for scenarios where large area contacts are expected to occur. Standard methods of indentation analyses use a beam theory solution to obtain an overall load–displacement relationship and then a Hertzian contact solution to calculate local stresses under the indenter. However, these techniques are only applicable in a fairly limited class of problems: the stress distribution in the contact region will differ significantly from a Hertzian one when the contact length exceeds the thickness of the beam. The indentation models developed herein are improvements over existing GLOBAL/LOCAL models for static and dynamic indentation of cantilever beams. Maximum contact stresses, beam displacements, and contact force time histories are obtained and compared with the predictions of current static and dynamic indentation models. The validity of the solutions presented herein is further assessed by comparing the results obtained to the predictions of modified beam theory solutions.     

Contact vibration analysis of the functionally graded material coated half-space under a rigid spherical punch

This paper studies the contact vibration problem of an elastic half-space coated with functionally graded materials (FGMs) subject to a rigid spherical punch. A static force superimposing a dynamic time-harmonic force acts on the rigid spherical punch. Firstly, we give the static contact problem of FGMs by a least-square fitting approach. Next, the dynamic contact pressure is solved by employing the perturbation method. Lastly, the dynamic contact stiffness with different dynamic contact displacement conditions is derived for the FGM coated half-space. The effects of the gradient index, coating thickness, internal friction, and punch radius on the dynamic contact stiffness factor are discussed in detail.     

考虑等效曲率的超二次曲面单元非线性接触模型

基于连续函数包络的超二次曲面单元可有效地描述自然界和工业生产中的非球体颗粒形态, 并通过非线性迭代方法精确计算单元间的接触力. 对于具有复杂几何形态的超二次曲面单元, 线性接触模型不能准确地计算不同接触模式下的作用力. 考虑超二次曲面单元相互作用时不同颗粒形状及表面曲率的影响, 本文发展了相应的非线性黏弹性接触模型. 该模型将不同接触模式下的法向刚度和黏滞力统一表述为单元间局部接触点处等效曲率半径的函数; 切向接触作用则借鉴基于Mohr-Coulomb摩擦定律的球体单元非线性接触模型的计算方法. 为检验超二次曲面单元接触模型的可靠性, 对球形颗粒间的法向碰撞、椭球体颗粒间的斜冲击过程、圆柱体的静态堆积和椭球体的动态卸料过程进行离散元模拟, 并与有限元数值结果及试验结果进行对比验证. 计算表明, 考虑接触点处等效曲率半径的超二次曲面非线性接触模型可准确地计算单元间的接触碰撞作用, 并合理地反映非球形颗粒体系的运动规律. 在此基础上进一步分析了不同长宽比和表面尖锐度对卸料过程中颗粒流动特性的影响, 为非球形颗粒材料的流动特性分析提供了一种有效的离散元方法.