Singularities of homogeneous deformations of constrained hyper-elastic materials

Singularity analysis is performed for homogeneous deformations of any hyper-elastic, constrained anisotropic material, under any type of conservative quasi-static loading. Critical conditions for branching of the equilibrium paths are defined and their post-critical behavior is discussed. Classification of the simple (cuspoids) and compound (umbilics) singularities of the total potential energy function is effected. The theory is implemented into an umbilic elliptical singularity of an isotropic and totally inextensible unit cube under normal loading.     

Constitutive branching in compressible elastic materials

Summary Branching analysis for the homogeneous deformations of a compressible elastic unit cube under dead loading is performed. Critical conditions for branching of the equilibrium paths are derived and the post-critical equilibrium paths are described. Special attention is given to the compound branching.     

DIVERGENCE INSTABILITY OF VARIABLE-ARC-LENGTH ELASTICA PIPES TRANSPORTING FLUID

A flexible elastic pipe transporting fluid is held by an elastic rotational spring at one end, while at the other end, a portion of the pipe may slide on a frictional support. Regardless of the gravity loads, when the internal flow velocity is higher than the critical velocity, large displacements of static equilibrium and divergence instability can be induced. This problem is highly nonlinear. Based on the inextensible elastica theory, it is solved herein via the use of elliptic integrals and the shooting method. Unlike buckling with stable branching of a simply supported elastica pipe with constant length, the variable arc-length elastica pipe buckles with unstable branching. The friction at the support has an influence in shifting the critical locus over the branching point. Alteration of the flow history causes jumping between equilibrium paths due to abrupt changes of direction of the support friction. The elastic rotational restraint brings about unsymmetrical bending configurations; consequently, snap-throughs and snap-backs can occur on odd and even buckling modes, respectively. From the theoretical point of view, the equilibrium configurations could be formed like soliton loops due to snapping instability.     

超细长弹性杆的分析力学问题

超细长弹性杆作为DNA等生物大分子链的力学模型，其平衡和稳定性问题已成为力学与分子生物学交叉的研究热点．虽然在Kirchhoff动力学比拟的基础上，用分析力学方法讨论弹性杆的文章已见诸文献，但尚未形成弹性杆分析力学的严格理论．本文研究了超细长弹性杆分析力学的若干基础性问题．对杆截面的自由度、虚位移、约束方程及约束力等基本概念给出严格的定义和表达式．建立弹性杆平衡的D’Alembert-Lagrange原理、Jourdain原理和Gauss原理；从D’Alembert-Lagrange原理导出Hamilton原理．从变分原理出发导出Lagrange方程、Nielsen方程、Appell方程和Hamilton正则方程；对于受约束的弹性杆，导出了带乘子的Lagrange方程．讨论了Lagrange方程的首次积分．对于杆中心线存在尖点的情形，导出了微段杆平衡的近似方程。     

非圆截面弹性细杆的平衡稳定性与分岔

本文研究存在初始曲率或挠率的非圆截面弹性细杆的平衡及稳定性问题,在两端受力矩单儿作用的条件下,杆的平衡微分方程可转换为用欧拉角表述的一阶自治系统,并有可能利用相平面的奇点理论分析弹性细杆平衡状态的稳定性,文中对杆截面的对称性,以及杆的初始曲率和挠率对平衡状态性的影响进行了定性分析,导出了解析形式的稳定性判据,揭示了杆平衡状态的列态分岔现象。     

A multiscale model for the ductile fracture of crystalline materials

In this paper, a multiscale model that combines both macroscopic and microscopic analyses is presented for describing the ductile fracture process of crystalline materials. In the macroscopic fracture analysis, the recently developed strain gradient plasticity theory is used to describe the fracture toughness, the shielding effects of plastic deformation on the crack growth, and the crack tip field through the use of an elastic core model. The crack tip field resulting from the macroscopic analysis using the strain gradient plasticity theory displayes the 1/2 singularity of stress within the strain gradient dominated region. In the microscopic fracture analysis, the discrete dislocation theory is used to describe the shielding effects of discrete dislocations on the crack growth. The result of the macroscopic analysis near the crack tip, i.e. a new K-field, is taken as the boundary condition for the microscopic fracture analysis. The equilibrium locations of the discrete dislocations around the crack and the shielding effects of the discrete dislocations on the crack growth at the microscale are calculated. The macroscopic fracture analysis and the microscopic fracture analysis are connected based on the elastic core model. Through a comparison of the shielding effects from plastic deformation and the discrete dislocations, the elastic core size is determined.     

Theory of images in plane-layered axisymmetric elastostatics

A fundamental theorem which automatically transforms the solution due to the presence of an arbitrary axisymmetric singularity in an unbounded homogeneous isotropic elastic solid into the corresponding solution for two perfectly bonded isotropic semi-infinite elastic solids is systematically applied in a stepwise fashion to obtain the complete image system when an arbitrary axisymmetric singularity is operative in or near a thick elastic layer which separates two other dissimilar isotropic semi-infinite elastic solids. The solution is closed and well structured. As an illustration, the distant effect in the interface layer produced by an influencing normal point force is luminously revealed to be two-dimensional, consisting of a combination at the origin of a bending hot spot and an infinite line of centres of dilatation. We conclude the paper with a complete theory of images for a free elastic layer under the influence of both axisymmetric and asymmetric singularities. We find that if the influencing displacement field is the gradient of a harmonic function, then the calculation of the induced elastic field reduces to the operation of differentiation only.     

The theory of Kirchhoff rods as an exact consequence of three-dimensional elasticity

A one-dimensional model of a linearly elastic thin rod is deduced from three-dimensional elasticity by regarding the Kirchhoff hypotheses as internal constraints prevailing in a three-dimensional tubular region. It follows from such an assumption that the displacement and the strain fields are linear in the cross-sectional coordinates. A constitutive relation that exhibits the maximal symmetry compatible with the assumed constraints is chosen and the equilibrium equations in terms of displacements are obtained.     

pH敏感凝胶平衡溶胀的有限元研究

1943年Flory提出第一个凝胶溶胀理论以来,凝胶溶胀理论的发展极大地促进了凝胶科学的发展．在诸多凝胶溶胀理论中,Flory-Rehner凝胶溶胀理论是最为经典的理论之一,成功地预测了聚合物胶体的诸多溶胀变形特性．但该模型是基于一种最为简单的自由连接链模型,存在精度较差的局限．因此,本文采用能够表征网络缠结拓扑等微观结构影响的弹性自由能模型,即由Edwards-Vilgis提出的Slip-Link模型,构造自由能模型并形成pH敏感凝胶的平衡溶胀理论,并基于Abaqus有限元分析了微观结构及几何限制作用对pH敏感凝胶力学行为的影响．分析结果可为微流体控制阀的设计提供参考．     

Local Existence for the FENE-Dumbbell Model of Polymeric Fluids

We study the well-posedness of a multi-scale model of polymeric fluids. The microscopic model is the kinetic theory of the finitely extensible nonlinear elastic (FENE) dumbbell model. The macroscopic model is the incompressible non-Newton fluids with polymer stress computed via the Kramers expression. The boundary condition of the FENE-type Fokker-Planck equation is proved to be unnecessary by the singularity on the boundary. Other main results are the local existence, uniqueness and regularity theorems for the FENE model in certain parameter range.     

Nonlinear in-plane thermal buckling of rotationally restrained functionally graded carbon nanotube reinforced composite shallow arches under uniform radial loading

The nonlinear in-plane instability of functionally graded carbon nanotube reinforced composite (FG-CNTRC) shallow circular arches with rotational constraints subject to a uniform radial load in a thermal environment is investigated. Assuming arches with thickness-graded material properties, four different distribution patterns of carbon nanotubes (CNTs) are considered. The classical arch theory and Donnell’s shallow shell theory assumptions are used to evaluate the arch displacement field, and the analytical solutions of buckling equilibrium equations and buckling loads are obtained by using the principle of virtual work. The critical geometric parameters are introduced to determine the criteria for buckling mode switching. Parametric studies are carried out to demonstrate the effects of temperature variations, material parameters, geometric parameters, and elastic constraints on the stability of the arch. It is found that increasing the volume fraction of CNTs and distributing CNTs away from the neutral axis significantly enhance the bending stiffness of the arch. In addition, the pretension and initial displacement caused by the temperature field have significant effects on the buckling behavior.     

A simple and effective nonlinear elastic one-dimensional model for the structural analysis of masonry arches

In this paper the static response of a masonry arch is studied by way of a one-dimensional nonlinear elastic model in which masonry is regarded as a material with bounded tensile and compressive strengths. By following an approach analogous to that followed in the theory of bending of elastic beams, the equilibrium problem for the arch leads to a free-boundary, nonlinear differential problem. An approximate solution to such problem can be pursued by means of an ad hoc iterative procedure, illustrated in detail herein. The results obtained in three case studies are compared with some numerical and experimental results available in the literature. In addition, the case of an actual arch undergoing spreading of the springings is considered, and the distribution and possible evolution of the cracking pattern discussed.     

三维切口应力奇性指数计算

将三维切口根部的位移渐近展开式引入线弹性力学平衡方程,导得关于切口应力奇性指数的特征微分方程组。再采用插值矩阵法,一次性地计算出三维切口的各阶应力奇性指数,它们具有同阶精度,并可同时获取相应的特征角函数。算例显示该法是分析三维切口应力奇异指数的一个有效的路径。计算结果表明,若直接用平面应变理论预测三维切口应力奇性指数将导致部分重要的奇性指数丢失。     

Dynamics of Holonomic Rigid-Body Systems

The theory behind analytical algorithms for computer-algebra systems is discussed. A computer system is described which is capable of optimizing input data, constructing equations of equilibrium and motion, formulating and solving the basic mechanics problems for a broad class of holonomic systems with elastic and dissipative constraints on the basis of the Lagrange-D'Alembert principle __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 7, pp. 130–138, July 2005.     

Singularity analysis of a two-dimensional elastic cable with 1：1 internal resonance

A two-degree-of-freedom bifurcation system for an elastic cable with 1：1 internal resonance is investigated in this paper. The transition set of the system is obtained with the singularity theory for three cases. The whole parametric plane is divided into several different persistent regions by the transition set. The bifurcation diagrams in different persistent regions are obtained.     

The singularity in weakly nonlinear fracture mechanics

Material failure by crack propagation essentially involves a concentration of large displacement-gradients near a crack's tip, even at scales where no irreversible deformation and energy dissipation occurs. This physical situation provides the motivation for a systematic gradient expansion of general nonlinear elastic constitutive laws that goes beyond the first order displacement-gradient expansion that is the basis for linear elastic fracture mechanics (LEFM). A weakly nonlinear fracture mechanics theory was recently developed by considering displacement-gradients up to second order. The theory predicts that, at scales within a dynamic lengthscale ℓ from a crack's tip, significant logr displacements and 1/r displacement-gradient contributions arise. Whereas in LEFM the 1/r singularity generates an unbalanced force and must be discarded, we show that this singularity not only exists but is also necessary in the weakly nonlinear theory. The theory generates no spurious forces and is consistent with the notion of the autonomy of the near-tip nonlinear region. The J-integral in the weakly nonlinear theory is also shown to be path-independent, taking the same value as the linear elastic J-integral. Thus, the weakly nonlinear theory retains the key tenets of fracture mechanics, while providing excellent quantitative agreement with measurements near the tip of single propagating cracks. As ℓ is consistent with lengthscales that appear in crack tip instabilities, we suggest that this theory may serve as a promising starting point for resolving open questions in fracture dynamics.     

Qualitative study of cavitated bifurcation for a class of incompressible generalized neo-Hookean spheres

IntroductionIn 1 958,GentandLindleyobservedthephenomenonofsuddenvoidnucleationinsolidsexperimentallyintensioningahomogenousclose_grainedvulcanizedrubbercylinderforthefirsttime.ButthemathematicalmodelonvoidnucleationandgrowthhasnotbeendescribedasabifurcationproblembasedonthetheoryofnonlinearelasticmechanicsbyBall[1]until1 982 .Inrecentyears,manyinvestigationshavebeenmadeonthisaspect.Theproblemofcavitatedbifurcationforincompressibleisotropichyperelasticmaterialswithpower_lawtypehasbeeninvestig…     

Modeling viscoelastic dielectrics

Dielectric elastomers, as an important category of electroactive polymers, are known to have viscoelastic properties that strongly affect their dynamic performance and limit their applications. Very few models accounting for the effects of both electrostatics and viscoelasticity exist in the literature, and even fewer are capable of making reliable predictions under general loads and constraints. Based on the principles of non-equilibrium thermodynamics, this paper develops a field theory that fully couples the large inelastic deformations and electric fields in deformable dielectrics. Our theory recovers existing models of elastic dielectrics in the equilibrium limit. The mechanism of instantaneous instability, which corresponds to the pull-in instability often observed on dielectric elastomers, is studied in a general non-equilibrium state. The current theoretical framework is able to adopt most finite-deformation constitutive relations and evolution laws of viscoelastic solids. As an example, a specific material model is selected and applied to the uniform deformation of a dielectric elastomer. This model predicts the stability criteria of viscoelastic dielectrics and its dependence on loading rate, pre-stress, and relaxation. The dynamic response, as well as the hysteresis behavior of a viscoelastic dielectric elastomer under cyclic electric fields, is also studied.     

Probability density functions of contact forces for cohesionless frictional granular materials

A theory is developed for the probability density functions of contact forces for cohesionless, frictional granular materials in quasi-static equilibrium. This theory is based on a maximum information entropy principle, with an expression for information entropy that is appropriate for granular materials. Entropy is maximized under the constraints of a prescribed stress and that the normal component of the contact force is compressive and that the tangential component of the contact force is limited by Coulomb friction. The theory results in a dependence of the probability density function for the tangential contact forces on the friction coefficient. The theoretical predictions are compared with results from discrete element simulations on isotropic, two-dimensional assemblies under hydrostatic stress. Good qualitative agreement is found for means and standard deviations of contact forces and the shape of the probability density functions, while the quantitative agreement is fairly good. Discrepancies between theory and simulations, such as the difference in shape of the probability density function for the normal force and the observed dependence on elastic properties of the exponential decay rate of tangential forces, are attributed to the fact that the method does not take into account any kinematics, which are essential in relation to elastic effects.     

Stress Singularity in an Elastic Cylinder of Cylindrically Anisotropic Materials

The issue of stress singularity in an elastic cylinder of cylindrically anisotropic materials is examined in the context of generalized plane strain and generalized torsion. With a viewpoint that the singularity may be attributed to a conflicting definition of anisotropy at r=0, we study the problem through a compound cylinder in which the outer cylinder is cylindrically anisotropic and the core is transversely isotropic. By letting the radius of the core go to zero, the cylinder becomes one with the central axis showing no conflict in the radial and tangential directions. Closed-form solutions are derived for the cylinder under pressure, extension, torsion, rotation and a uniform temperature change. It is found that the stress is bounded everywhere, and singularity does not occur if the anisotropy at r=0 is defined appropriately. This revised version was published online in July 2006 with corrections to the Cover Date.