MOTION OF A HANGING CHAIN AFTER AN INITIAL HORIZONTAL IMPACT AT THE FREE END

悬链在自由端受到冲击后的瞬态响应,不仅是一个有趣的数学问题,同时也具有一定的工程意义.在拉格朗日动力学微分方程的理论框架下,引入广义冲量并利用第二类拉格朗日方程对悬链在自由端受水平冲击力后的动力学响应进行了分析,得到了计算每节链段角速度的统一公式. 应用该公式能够方便地求解不同初始条件下具有较多链段数目的悬链在冲击作用下的瞬态响应问题.     

基于Slip-Link模型的凝胶大变形与扩散耦合瞬态问题分析

长链聚合物通过化学键,范德华力等相互作用交联构成的弹性体吸收溶剂后形成聚合物凝胶.聚合物凝胶的溶胀变形过程是溶剂在凝胶网络内的扩散、迁移与凝胶网络变形耦合作用的结果.通常,研究凝胶的扩散-变形耦合问题是基于经典的Flory-Rehner自由能模型进行的,然而该自由能函数由于忽略了凝胶网络链段相互缠绕引起的物理交联对凝胶弹性网络自由能的影响,在计算大变形与扩散耦合问题时误差较大.本文基于Edwards-Vilgis slip-link弹性模型和Flory-Huggins混合理论构建能反映聚合物网络链缠结拓扑限制作用的凝胶自由能函数,并结合凝胶的运动方程和扩散方程,得到能够反映链段缠结影响的凝胶大变形与扩散耦合瞬态方程,并基于此研究网络缠结对凝胶在受压状态下变形梯度和凝胶内溶剂的名义浓度的影响.     

爆炸荷载作用下弹性与阻尼支承梁的动力响应

为了提高梁的抗爆能力,目前通常是从材料、截面形状等方面考虑．通过在梁端部设置弹性和阻尼支承来提高梁的抗力．为此,应用拉格朗日方程建立了端部弹性与阻尼支承梁的动力方程,通过长作用时间和短作用时间两个爆炸动荷载的典型算例,分析了端部弹性与阻尼支承对梁动态响应的影响,并指出在短作用时间动荷载作用下在梁端部设置弹性与阻尼支承可有效提高梁的抗力．     

密度梯度柱壳链的弹性波传播特性研究

均匀圆柱壳链可以调控弹性波传播, 引入密度梯度有望进一步提高波形调控能力. 通过建立密度梯度柱壳链的细观有限元模型和连续介质模型, 研究了密度梯度柱壳链的弹性波传播特性. 通过将密度梯度柱壳链等效为变密度连续介质弹性杆, 建立了其在应力脉冲作用下的控制方程. 运用拉普拉斯积分变换方法, 考虑杆中密度遵循线性分布, 获得了方程的解析解. 以三角形应力脉冲作用为例, 通过与细观有限元模拟结果比较, 发现解析解可以较好地预测梯度柱壳链中载荷的演化趋势. 正梯度链中载荷峰值随着波传播逐渐增大, 负梯度链中载荷峰值随着波传播逐渐减小. 正梯度链支撑端峰值载荷高于均匀链, 负梯度链支撑端峰值载荷低于均匀链, 这表明相较于均匀柱壳链, 密度梯度柱壳链可以在更大范围内对波形进行调控. 线性密度梯度参数对梯度柱壳链的波形调控能力影响较大, 梯度参数越小, 传递到支撑端的峰值载荷越小; 相反, 梯度参数越大, 支撑端的峰值载荷越大. 建立的理论模型及其解析解为研究梯度柱壳链中应力波传播规律及揭示载荷调控机理提供了理论基础.      

双材料自由边界面端部问题的数值分析

对各向异性双材料自由边界面端部奇异性场问题进行了研究,利用有限元分析法所得到的各向异性双材料自由边界面端部的应力奇异性指数以及角分布函数,构造了一个自由边界面端部单元,据此建立了自由边界面端部奇异性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型.与四节点单元相结合,提出一种求解自由边界面端部广义应力强度因子的杂交元法.考核例结果表明:本文方法的数值解精度高,可应用于各向异性材料双材料自由边界面端部问题.     

DNA纳米管的扭转刚度模型

刚度是衡量材料弹性变形难易程度的一个定量表征参数,与DNA纳米管静动力学特性及其结构生物功能密切相关．本文致力于研究DNA纳米管的扭转刚度．首先,在六角形均匀封装条件下,考虑到单个DNA杆件弯扭组合问题的静不定特点,我们利用平衡方程、变形协调方程和弹性本构方程,合理预测了DNA纳米管扭转实验中单个DNA杆件的弯扭组合变形,由此给出了DNA纳米管扭转刚度预测的解析模型．最后的结果表明：随着DNA杆数的增加,DNA纳米管的弯曲刚度显著增加,而其扭转刚度却几乎不变,合理解释了扭转实验中发现的现象．有关结论为DNA折叠结构的设计和应用提供了参考．     

马氏体相变过程中形核、失稳及多界面微结构的动力学模拟

材料发生相变的过程中会出现失稳、滞后回线及多界面的微结构等复杂现象，而稳定性的丧失使其动力学方程的求解十分困难。对于形状记忆合金中的马氏体相变，相变过程中材料的等效杨氏模量变为负值，使得传统的动力学方程成为病态的，无法直接求解，必须要进行正则化。而相变的滞后回线与微结构的出现也说明经典的弹性理论不再适用，必须要引入新的能量项以能刻画这些现象。本文在非线性弹性理论的框架下，引入应变梯度界面能和位移非均匀能，利用变分原理建立了材料相变的一维动力学模型。高阶项的引入极大地改善了方程的性质，使数值求解成为可能。计算结果表明，该模型确能较有效地描述相变时的失稳与微结构。     

纤维与基体的轴对称界面端附近的奇异应力场

基于各向性弹性力学空间轴对称问题的基本方程,研究了纤维与基体的轴对称界面端的应力奇异性,并给出了界最佳 近的奇异应力场。研究结果表明,该轴对称界面端的应力奇异性与平面应变状态下相应模型的应力奇异性完全相同,材料对界面端附近奇异应力场的影响可用丰个双材料组合参数描述。     

微振激励下黏弹性阻尼器微观链结构力学模型

减小微振动对高精密仪器至关重要,利用黏弹性阻尼器进行微振动抑制是一个新兴而又具有挑战性的课题.本文采用分子链网络模型方法分析了黏弹性材料的微观分子链结构,综合考虑材料分子链结构中的网络链和自由链对黏弹性材料力学性能的影响,提出一种基于材料微观分子链结构的微振激励下黏弹性阻尼器力学模型.模型分别采用标准线性固体模型和Maxwell模型来描述网络链和自由链中单个链的力学性能,并分别采用8链网络模型和3链网络模型考虑两种类型分子链的综合效应,引入温频等效原理描述温度对微振激励下黏弹性阻尼器力学性能的影响.该模型能够描述温度和频率对黏弹性阻尼器动态力学性能的影响,并能够反映黏弹性材料的微观结构与材料力学性能的关系.为验证所提模型的有效性及考察黏弹性阻尼器在微振激励下的耗能能力和动态力学性能,在微振条件下对黏弹性阻尼器进行了动态力学性能试验.研究结果表明黏弹性阻尼器具有较好的微振耗能能力,其动态力学性能受温度和频率影响较大,所提的力学模型能够精确地描述微振激励下黏弹性阻尼器动态力学性能随温度和频率的变化关系.      

超大变形弹性细杆几何形态的拓扑描述

 以DNA和其他生物高分子链为背景的弹性杆模型以其极端细长性和超大变形而不同于传统 弹性力学的研究对象. 杆在外力作用下的几何形态取决于影响中心线形状的弯曲变形及截面 绕中心线的扭转变形. 讨论了超大变形弹性细杆几何形态的拓扑学描述，及其在弹性杆的 平衡和稳定性问题中的应用.     

Size-dependent effective elastic moduli of particulate composites with interfacial displacement and traction discontinuities

This work aims at estimating the size-dependent effective elastic moduli of particulate composites in which both the interfacial displacement and traction discontinuities occur. To this end, the interfacial discontinuity relations derived from the replacement of a thin uniform interphase layer between two dissimilar materials by an imperfect interface are reformulated so as to considerably simplify the characteristic expressions of a general elastic imperfect model which is adopted in the present work and include the widely used Gurtin–Murdoch and spring-layer interface models as particular cases. The elastic fields in an infinite body made of a matrix containing an imperfectly bonded spherical particle and subjected to arbitrary remote uniform strain boundary conditions are then provided in an exact, coordinate-free and compact way. With the aid of these results, the elastic properties of a perfectly bonded spherical particle energetically equivalent to an imperfectly bonded one in an infinite matrix are determined. The estimates for the effective bulk and shear moduli of isotropic particulate composites are finally obtained by using the generalized self-consistent scheme and discussed through numerical examples.     

平行于功能梯度材料夹层的币型裂纹起裂条件

分析了功能梯度材料中币型裂纹的扩展问题．裂纹平行于无限域中功能梯度材料夹层,受有与裂纹面成任意角度的拉应力．假定功能梯度材料夹层与两个半无限域均匀介质完全粘合,其弹性模量沿厚度方向变化．采用基于层状材料广义Kelvin基本解的边界元方法分析裂纹问题,给出了均布正应力和剪应力作用下裂纹的应力强度因子、将应力强度因子耦合于应变能密度断裂判据,讨论了裂纹体在拉伸应力作用下的起裂条件．     

Geometric effects of cross sections on equilibrium of helical and twisted ribbon

In the framework of elastic rod model, the Euler-Lagrange equations characterizing the equilibrium configuration of the polymer chain are derived from a free energy functional associated with the curvature, torsion, twisting angle, and its derivative with respect to the arc-length. The configurations of the helical ribbons with different crosssectional shapes are given. The effects of the elastic properties, the cross-sectional shapes,and the intrinsic twisting on the helical ribbons are discussed. The results show that the pitch angle of the helical ribbon decreases with the increase in the ratio of the twisting rigidity to the bending rigidity and approaches the intrinsic twisting. If the bending rigidity is much greater than the twisting rigidity, the bending and twisting of the helical ribbon always appear simultaneously.     

On flat thin elastic rods with rapidly varying periodic characteristics

A thin elastic rod is considered on a plane. The free shape of the rod is described by a periodic curve. It is shown that, under constant loads, its equilibrium shape tends to the equilibrium shape of a thin rectilinear rod when the frequency of the function describing the free shape increases infinitely. The problem under study is solved on the basis of modeling the three-dimensional shapes of circular DNA molecules by a thin rectilinear elastic rod.     

弹性杆基因模型的力学问题

概述弹性杆静力学与刚体动力学之间相似性的Kirchhoff理论．讨论其在分子生物学的弹性杆基因模型中的应用，以及与分析力学和运动稳定性理论有关的若干问题．     

Interfacial Energies for Incoherent Inclusions

We study a variational problem describing an incoherent interface between a rigid inclusion and a linearly elastic matrix. The elastic material is allowed to slip relative to the inclusion along the interface, and the resulting mismatch is penalized by an interfacial energy term that depends on the surface gradient of the relative displacement. The competition between the elastic and interfacial energies induces a threshold effect when the interfacial energy density is non-smooth: small inclusions are coherent (no mismatch); sufficiently large inclusions are incoherent. We also show that the relaxation of the energy functional can be written as the sum of the bulk elastic energy functional and the tangential quasiconvex envelope of the interfacial energy functional.     

Dynamic symmetry breaking and structure-preserving analysis on the longitudinal wave in an elastic rod with a variable cross-section

The longitudinal wave propagating in an elastic rod with a variable cross-section owns wide engineering background, in which the longitudinal wave dissipation determines some important performances of the slender structure. To reproduce the longitudinal wave dissipation effects on an elastic rod with a variable cross-section, a structure-preserving approach is developed based on the dynamic symmetry breaking theory. For the dynamic model controlling the longitudinal wave propagating in the elast...     

Separation phenomena at the interface of a finitely deformed composite sphere

The problem of the finite deformation of a composite sphere subjected to a spherically symmetric dead load traction is revisited focusing on the formation of a cavity at the interface between a hyperelastic, incompressible matrix shell and a rigid inhomogeneity. Separation phenomena are assumed to be governed by a vanishingly thin interfacial cohesive zone characterized by uniform normal and tangential interface force–separation constitutive relations. Spherically symmetric cavity shapes (spheres) are shown to be solutions of an interfacial integral equation depending on the strain energy density of the matrix, the interface force constitutive relation, the dead loading and the volume concentration of inhomogeneity. Spherically symmetric and non-symmetric bifurcations initiating from spherically symmetric equilibrium states are analyzed within the framework of infinitesimal strain superimposed on a given finite deformation. A simple formula for the dead load required to initiate the non-symmetrical rigid body mode is obtained and a detailed examination of a few special cases is provided. Explicit results are presented for the Mooney–Rivlin strain energy density and for an interface force–separation relation which allows for complete decohesion in normal separation.     

A Finite Strain Analysis of Cavity Formation at a Rigid Inhomogeneity

The formation of a cavity by inclusion-matrix interfacial separation is examined by analyzing the response of a plane rigid inclusion embedded in an unbounded incompressible matrix subject to remote equibiaxial dead load traction. A vanishingly thin interfacial cohesive zone, characterized by normal and tangential interface force-separation constitutive relations, is assumed to govern separation behavior. Rotationally symmetric cavity shapes (circles) are shown to be solutions of an interfacial integral equation depending on the strain energy density of the matrix, the interface force constitutive relation and the remote loading. Nonsymmetrical cavity formation, under rotationally symmetric conditions of geometry and loading, is treated within the theory of infinitesimal strain superimposed on a given finite strain state. Rotationally symmetric and nonsymmetric bifurcations are analyzed and detailed results, for the Mooney–Rivlin strain energy density and for an exponential interface force-separation law, are presented. For the nonsymmetric rigid body displacement mode, a simple formula for the critical load is presented. The effect on bifurcation behavior of interfacial shear stiffness and other interface parameters is treated as well. In particular we demonstrate that (i) for the smooth interface nonsymmetric bifurcation always precedes rotationally symmetric bifurcation, (ii) unlike rotationally symmetric bifurcation, there is no threshold value of interface parameter for which nonsymmetric bifurcation will not occur and (iii) interfacial shear may significantly delay the onset of nonsymmetric bifurcation. Also discussed is the range of validity of a nonlinear infinitesimal strain theory previously presented by the author (Levy [1]). This revised version was published online in July 2006 with corrections to the Cover Date.     

Modeling of properties and motions of an inhomogeneous one-dimensional continuum of a complicated Cosserat-type microstructure

We consider an approach to modeling the properties of the one-dimensional Cosserat continuum [1] by using the mechanical modeling method proposed by Il’yushin in [2] and applied in [3]. In this method, elements (blocks, cells) of special form are used to develop a discrete model of the structure so that the average properties of the model reproduced the properties of the continuum under study. The rigged rod model, which is an elastic structure in the form of a thin rod with massive inclusions (pulleys) fixed by elastic hinges on its elastic line and connected by elastic belt transmissions, is taken to be the original discrete model of the Cosserat continuum. The complete system of equations describing the mechanical properties and the dynamical equilibrium of the rigged rod in arbitrary plane motions is derived. These equations are averaged in the case of a sufficiently smooth variation in the parameters of motion along the rod (the long-wave approximation). It was found that the average equations exactly coincide with the equations for the one-dimensional Cosserat medium [1] and, in some specific cases, with the classical equations of motion of an elastic rod [4–6]. We study the plane motions of the one-dimensional continuum model thus constructed. The equations characterizing the continuum properties and motions are linearized by using several assumptions that the kinematic parameters are small. We solve the problem of natural vibrations with homogeneous boundary conditions and establish that each value of the parameter distinguishing the natural vibration modes is associated with exactly two distinct vibration mode shapes (in the same mode), each of which has its own frequency value.