Global dynamics of the cable under combined parametrical and external excitations

The chaotic dynamics and global bifurcations of the suspended elastic cable under combined parametric and external excitations are investigated. The non-linear equations of motion of the elastic cable to small vibration of one support are derived. The averaged equations are obtained by using the method of multiple scales. Based on the averaged equations, the theory of normal form and Maple program are used to obtain the explicit expressions of normal form associated with a double zero and a pair of pure imaginary eigenvalues. On the basis of the normal form, global bifurcation analysis of the parametrically and externally excited suspended elastic cable is given by a global perturbation method developed by Kovacic and Wiggins. The chaotic motion of the elastic cable is also found by numerical simulation.     

Non-local finite element analysis of damped beams

Non-local viscoelastic beam models are used to analyse the dynamics of beams with different boundary conditions using the finite element method. Unlike local damping models the internal force of the non-local model is obtained as weighted average of state variables over a spatial domain via convolution integrals with spatial kernel functions that depend on a distance measure. In the finite element analysis, the interpolating shape functions of the element displacement field are identical to those of standard two-node beam elements. However, for non-local damping, nodes remote from the element do have an effect on the energy expressions, and hence on the damping matrix. The expressions of these direct and cross damping matrices may be obtained explicitly for some common spatial kernel functions and Euler–Bernoulli beam theory. Alternatively numerical integration may be applied to obtain solutions. Examples are given where the eigenvalues are compared to the exact solution for a pinned–pinned beam to demonstrate the convergence of the finite element method. The results for beams with other boundary conditions are used to demonstrate the versatility of the finite element technique.     

Shape sensitivity analysis and the energy momentum tensor for the kinematic and static models of torsion

The aim of this paper is to bridge shape sensitivity analysis and configurational mechanics by means of a widespread use of the shape derivative concept. This technique will be applied as a systematic procedure to obtain the Eshelby’s energy momentum tensor associated to the problem under consideration. In order to highlight special features of this procedure and without loss of generality, we focus our attention in the application of shape sensitivity analysis to the problem of twisted straight bars within the framework of linear elasticity.Kinematic and static variational formulations as well as the direct method of sensitivity analysis are used to perform shape derivatives of both models. Integral expressions of first and second order shape derivatives of the total potential energy and the complementary potential energy with respect to an arbitrary transverse cross-section shape change, are achieved. These integral expressions put in evidence the relationship between shape sensitivity analysis and the first and second order Eshelby’s energy momentum tensors. Also, the null divergence property of these tensors is easily proved by comparing, in each case, the domain and boundary integral shape derivative arrived at. Finally, an example with a known exact solution, corresponding to an elastic bar with elliptical transverse cross-section submitted to twist, is presented in order to illustrate the usefulness of these tensors to compute the corresponding shape derivatives.     

A three-dimensional inverse problem for inhomogeneities in elastic solids

The Newtonian potential is used to solve an inverse problem in which we seek the shape of an inhomogeneity in an infinite elastic matrix under uniform applied stresses at infinity such that certain stress components are uniform on the boundary of the inhomogeneity. It is shown that ellipsoids furnish the solution of this inverse problem. Exact and general expressions for the stress and displacement are given explicitly for points in the elastic matrix outside the inhomogeneity. The solution of the corresponding plane deformation problem is found as a limiting case. Several applications are presented, and results from the literature are confirmed as special cases.     

Optimal shape of a twisted and compressed rod

We formulate and solve the problem of determining the shape of an elastic rod stable against buckling and having minimal volume. The rod is loaded by a concentrated force and a couple at its ends. The equilibrium equations are reduced to a single nonlinear second-order equation. The eigenvalues of the linearized version of this equation determine the stability boundary. By using Pontryagin's maximum principle we determine the optimal shape of the rod.     

Elastic constants for granular materials modeled as first-order strain-gradient continua

Formulation of a stress–strain relationship is presented for a granular medium, which is modeled as a first-order strain-gradient continuum. The elastic constants used in the stress–strain relationship are derived as an explicit function of inter-particle stiffness, particle size, and packing density. It can be demonstrated that couple-stress continuum is a subclass of strain-gradient continua. The derived stress–strain relationship is simplified to obtain the expressions of elastic constants for a couple-stress continuum. The derived stress–strain relationship is compared with that of existing theories on strain- gradient models. The effects of inter-particle stiffness and particle size on material constants are discussed.     

On the stability of an elastic column positioned on an elastic half space

Summary Stability of a heavy elastic column loaded by a concentrated force at the top is analysed. It is assumed that the column is fixed to a rigid circular plate that is positioned on a homogeneous, isotropic, linearly elastic half-space. The constitutive equations for the column are assumed in the form that allows axial compressibility and takes into account the influence of shear stresses. It is shown that eigenvalues of the linearized equations determine the bifurcation points of the full non-linear system of equilibrium equations. Also, the type of bifurcation at the lowest eigenvalue is examined and shown that it could be both super-and sub-critical. The post-critical shape of the column is determined by numerical integration of the equilibrium equations. Received 13 June 1998; accepted for publication 12 November 1998     

Effective properties of elastic periodic composite media with fibers

A series solution to obtain the effective properties of some elastic composites media having periodically located heterogeneities is described. The method uses the classical expansion along Neuman series of the solution of the periodic elasticity problem in Fourier space, based on the Green's tensor, and exact expressions of factors depending on the shape of the inclusions. Some properties of convergence of the solution are presented, more specifically concerning the elasticity tensor of the reference medium, showing that the convergence occurs even for empty fibers. The solution is extended for rigid inclusions. A comparison is made with previous exact solutions for a fiber composite made of cylindrical fibers with circular cross-sections and with previous estimates. Different examples are presented for new situations concerning the study of fiber composites: composites with elliptic cross-sections and multi-phase fibrous composites.     

Fourth-order quasi-harmonic equations of state for solids of cubic and tetragonal symmetry

General expressions are given for the dependence of the pressure and the effective elastic moduli on deformation and temperature in the form of a Taylor series expansion with respect to elastic and thermal strains. The temperature dependence of these expressions is derived within the quasi-harmonic approximation of lattice dynamics. The expressions are developed in terms of the Lagrangian strain and an alternative strain measure identical with the Eulerian strain for a pure deformation. They are then used to obtain the third- and fourth-order equations of state for crystals of cubic and tetragonal symmetry and to relate the parameters entering these equations to quantities which are commonly (or may be potentially) measured experimentally. It is shown that available ultrasonic data are not completely sufficient to evaluate the parameters of fourth-order equations of state. For tetragonal symmetry, this problem is still in abeyance; while in the cubic case, it is possible to estimate the fourth-order parameters from shock-wave data and so to give illustrative numerical applications of our equations. Finally, the third- and fourth-order Hugoniots and isotherms of Cu and Ag are calculated in terms of both the Lagrangian and Eulerian strain measures.     

Applications of the method of complex functions to dynamic stress concentrations

The complex function method used in the solution of static stress concentration around an irregularly shaped cavity in an infinite elastic plane is generalized to the case of dynamic loading. This paper presents the solutions of two dimensional elastic wave equations in terms of complex wave functions, and general expressions for boundary conditions for steady state incident waves. Dynamic stresses around a cavity of arbitrary shape are then expressed in series of complex ‘domain functions’, the coefficient of the series can be determined by truncating a set of infinite algebraic equations. Results of dynamic stress concentration factors for circular and elliptical cavities are given in this paper.     

Vibrations of a continuous web on elastic supports

We consider an infinite, homogenous linearly elastic beam resting on a system of linearly elastic supports, as an idealized model for a paper web in the middle of a cylinder-based dryer section. We obtain closed-form analytical expressions for the eigenfrequencies and the eigenmodes. The frequencies increase as the support rigidity is increased. Each frequency is bounded from above by the solution with absolutely rigid supports, and from below by the solution in the limit of vanishing support rigidity. Thus in a real system, the natural frequencies will be lower than predicted by commonly used models with rigid supports.     

CONSTITUTIVE MODELING OF ROLLED SHAPE MEMORY ALLOY SHEETS TAKING INTO ACCOUNT PRE-TEXTURE AND ANISOTROPIC HARDENING

The drawing or rolling process endows polycrystal shape memory alloy with a crys- tallographic texture, which can result in macroscopic anisotropy. The main purpose of this work is to develop a constitutive model to predict the thermomechanical behavior of shape memory alloy sheets, which accounts for the crystallographic texture. The total macroscopic strain is decom- posed into elastic strain and macro-transformation strain under isothermal condition. Considering the transformation strain in local grains and the orientation distribution function of crystallo- graphic texture, the macro-transformation strain and the effective elastic modulus of textured polycrystal shape memory alloy are developed by using tensor expressions. The kinetic equation is established to calculate the volume fraction of the martensite transformation under given stress. Furthermore, the Hill＇s quadratic model is developed for anisotropic transformation hardening of textured SMA sheets. All the calculation results are in good agreement with experimental data, which show that the present model can accurately describe the macro-anisotropic behaviors of textured shape memory alloy sheets.     

On reflected interactions in elastic solids containing inhomogeneities

Interactions in linear elastic solids containing inhomogeneities are examined using integral equations. Direct and reflected interactions are identified. Direct interactions occur simply because elastic fields emitted by inhomogeneities affect each other. Reflected interactions occur because elastic fields emitted by inhomogeneities are reflected by the specimen boundary back to the individual inhomogeneities. It is shown that the reflected interactions are of critical importance to analysis of representative volume elements. Further, the reflected interactions are expressed in simple terms, so that one can obtain explicit approximate expressions for the effective stiffness tensor for linear elastic solids containing ellipsoidal and non-ellipsoidal inhomogeneities. For ellipsoidal inhomogeneities, the new approximation is closely related to that of Mori and Tanaka. In general, the new approximation can be used to recover Ponte Castañeda–Willis׳ and Kanaun–Levin׳s approximations. Connections with Maxwell׳s approximation are established.     

基于扭曲Kagome点阵结构的一模材料设计

零能模式超材料指弹性矩阵的特征值中有若干为零的弹性材料,根据零特征值的个数可将其分类为一模至五模材料。当前,针对五模材料已有较深入研究,并在水声和弹性波调控方面获得重要应用,而对其他类型零能模式材料的研究尚未展开。本文对扭曲Kagome周期桁架这样一类欠约束点阵材料的有效弹性性质进行了研究,结果表明通过调节点阵材料的微观几何构型和杆件刚度,该类结构能够涵盖一系列一模材料谱系。针对给定一模弹性张量,发展了软-硬模式分离的微结构逆向优化设计策略。通过特定一模材料中的波传播现象对有效性质预测和微结构设计进行了数值验证。     

含正交排列夹杂和缺陷材料的等效弹性模量和损伤

研究含正交排列夹杂和缺陷材料的等效弹性模量和损伤,推导了以Eshelby－Mori－Tanaka方法求解多相各向异性复合材料等效弹性模量的简便计算公式,针对含三相正交椭球状夹杂的正交各向异性材料,得到了由细观参量（夹杂的形状、方位和体积分数）表示的等效弹性模量的解析表达式．在此基础上,提出了一个宏细观结合的正交各向异性损伤模型,从而建立了以细观量为参量的含损伤材料的应力应变关系．最后,对影响材料损伤的细观结构参数进行了分析．     

二类复特征值表达式的适用性研究

对于线性系统，根据其质量、刚度或阻尼等矩阵的性质了解特征值，是十分有意义的问题．对于有阻尼陀螺系统，本文首先通过模态正交性推导了有关矩阵的若干关系式．然后分析了复特征值的两类表达式的适用性，着重讨论了通过一元二次方程研究特征值性质的现行处理方法存在的困难．提出一种确定特征值的补充方法，指出可能存在虚数特征值的系统的性质．算例表明结论正确     

Crack edge displacement and elastic constant determination for an orthotropic material

A method for determining the elastic constants of an isotropic material, based on crack edge displacement data, is extended to an orthotropic material. Complex potentials are used to obtain the stresses and displacements for plane strain. Mode I crack problems in three mutually orthogonal planes are considered and solved. In particular, the expressions of crack edge displacements are obtained in an explicit form. An iterative statistical identification method, based on a Bayesan approach, is used to identify the elastic constants of an orthotropic medium from the Mode I crack displacements measured from the mid-point of the crack. Some graphics are displayed to illustrate the convergence of the pertinent parameters and the approach of the analytical displacements to their experimental values.     

Determination of optimal shape of a moving punch with friction taken into account

The optimization problem for the contact interaction between a rigid punch and an elasticmediumis considered. It is assumed that that the punch is under the action of some prescribed forces and momenta and moves along a surface bounding a half-space filled with an elastic medium. It is also assumed that themotion is quasistatic and the friction forces arising in the region of contact are taken into account. The punch shape is considered as the desired design variable, and the integral functional characterizing the discrepancy between the pressure distribution in the region of contact that corresponds to the optimized shape of the punch and a given goal distribution of pressure is taken as the minimizing criterion. The optimal shape can be determined efficiently by solving the following two problems: first, to obtain the optimal pressure distribution and then to solve a boundary value problemfor the elastic half-space under the action of the obtained normal pressure and friction forces. By way of example, the optimal shape is analytically determined for a punch of rectangular shape in horizontal projection.     

Electromechanical response of (3–0, 3–1) particulate,fibrous, and porous piezoelectric composites with anisotropic constituents: A model based on the homogenization method

An analytical framework based on the homogenization method has been developed to predict the effective electromechanical properties of periodic, particulate and porous, piezoelectric composites with anisotropic constituents. Expressions are provided for the effective moduli tensors of n-phase composites based on the respective strain and electric field concentration tensors. By taking into account the shape and distribution of the inclusion and by invoking a simple numerical procedure, solutions for the electromechanical properties of a general anisotropic inclusion in an anisotropic matrix are obtained. While analytical forms are provided for predicting the electroelastic moduli of composites with spherical and cylindrical inclusions, numerical evaluation of integrals over the composite microstructure is required in order to obtain the corresponding expressions for a general ellipsoidal particle in a piezoelectric matrix. The electroelastic moduli of piezoelectric composites predicted by the analytical model developed in the present study demonstrate excellent agreement with results obtained from three-dimensional finite-element models for several piezoelectric systems that exhibit varying degrees of elastic anisotropy.