热超弹性薄板非对称拉伸的分岔和稳定性

本文应用超弹性材料的有限变形理论分析了在面内等双向拉伸死载荷作用下不可压热超弹性方形薄板发生非对称变形的分岔及其稳定性问题.给出了方板变形的分岔曲线和临界载荷,发现对受面内等双向拉伸载荷作用的均匀方板,当拉伸载荷值较小时,方板双向等伸长变形,发生对称的拉伸变形;但当此载荷值大于某一临界值时,从方板的对称拉伸变形中分岔出非对称的变形,方板在两个方向的变形不再相等.通过变形发生分岔前后的能量比较发现,分岔后的对称变形是不稳定的,而非对称变形是稳定的.同时,给出了板中的应力分布曲线,并由不同温度下变形的分岔曲线和应力分布曲线讨论了温度对方板变形和板中的应力分布的影响.     

内压载荷作用下薄膜椭球的稳定性分析

超弹性橄榄状和南瓜状薄膜椭球在内压载荷作用下存在不同的分岔形式.对橄榄状薄膜椭球来说,细长比大于某一临界值时,在一定内压作用下会发生梨形分岔;小于该临界值时,薄膜椭球的分岔行为与圆管的局部起鼓现象相类似.对南瓜状薄膜椭球,无论圆扁,当内压达到某载荷值时都会发生梨形分岔.本文采用能量判据,分析了在压强控制和质量控制两种加载方式作用下,不同形状的薄膜椭球的均匀解及分岔解的稳定性.通过计算要考察的平衡状态及施加小扰动之后状态的能量差来判断当前状态是否稳定,结果表明,在压强控制下,P-V曲线下降段的均匀解和分岔解均为不稳定解.但在质量控制下,在P-V曲线下降段中只有均匀解出现时,均匀解为稳定解;而在均匀解和分岔解共存的区间内,均匀解为不稳定解,分岔解为稳定解.另外,P-V曲线两个上升段的均匀解则均为稳定解.     

Nonlinear Oscillations of a Nonresonant Cable under In-Plane Excitation with a Longitudinal Control

The nonlinear oscillations of a controlled suspended elastic cable under in-plane excitation are considered. Active control realized by longitudinal displacement of one support is introduced in order to reduce the transverse in-plane and out-of-plane vibrations. Linear and quadratic enhanced velocity feedback control laws are chosen and their effects on the cable motion are investigated using a two degree-of-freedom model. Perturbation analysis is performed to determine the in-plane steady-state solutions and their stability under an out-of-plane disturbance. The analysis is extended to the bifurcated two-mode steady-state oscillations in the region of parametric excitation. The dependence of the control effectiveness on the system parameters is investigated in the case of the first symmetric mode and the range of oscillation amplitudes in which the proposed control ensures a dissipation of energy is determined. Although control based only on in-plane response quantities is effective in reducing oscillations with a prevailing in-plane component, addition of out-of-plane measures has to be considered when the motion is characterized by two comparable components.     

On the Investigation of Some Parameter Identification and Experimental Modal Filtering Issues for Nonlinear Reduced Order Models

This paper discusses modal filtering of experimental data and the corresponding identification of linear and nonlinear parameters in reduced order space. Specifically, several experimental configurations will be discussed in order to provide insight into such identification issues as spatial discretization, observability, and the linear independence of the assumed filter or basis. The two experiments considered herein represent different measurement configurations of the same clamped–clamped beam. First, asymmetric inertial loading via asymmetric sensor location was considered, while the second scenario presents a symmetric sensor configuration. Several important conclusions can be drawn from the two experimental scenarios. First, by asymmetrically loading the beam, a corresponding asymmetric beam mode was excited yet not observable. In the second scenario, the symmetric distribution of sensors minimized the impact of the respective asymmetric mode. The resulting spatial information allowed for the proper filtering of the remnants of the asymmetric mode. Nonlinear parameters in modal space as well as the underlying linear parameters were successfully identified simultaneously in both experimental scenarios, although the usefulness of the asymmetrically loaded beam was limited. Finally, successful comparisons were made between the identified reduced order model and experimental response at the beam quarter point using the symmetric case and the beam midpoint using both experimental scenarios.     

Brittle fracture dynamics with arbitrary paths III. The branching instability under general loading

The dynamic propagation of a bifurcated crack under arbitrary loading is studied. Under plane loading configurations, it is shown that the model problem of the determination of the dynamic stress intensity factors after branching is similar to the anti-plane crack branching problem. By analogy with the exact results of the mode III case, the energy release rate immediately after branching under plane situations is expected to be maximized when the branches start to propagate quasi-statically. Therefore, the branching of a single propagating crack under mode I loading should be energetically possible when its speed exceeds a threshold value. The critical velocity for branching of the initial single crack depends only weakly on the criterion applied for selecting the paths followed by the branches. However, the principle of local symmetry imposes a branching angle which is larger than the one given by the maximum energy release rate criterion. Finally, it is shown that an increasing fracture energy with the velocity results in a decrease in the critical velocity at which branching is energetically possible.     

Stability of Cylindrical Composite Shells under Torsion

The stability of fiber-reinforced cylindrical shells under torsion is analyzed in the case where the principal directions of elasticity in the layers do not coincide with the coordinate directions. The solution to the linearized equations of the technical theory of anisotropic shells is obtained in the form of trigonometric series. It is shown that for some reinforcement configurations the critical loads may depend on the direction of the torsional moment. It is also established that the minimum (in absolute value) eigenvalue does not always correspond to the critical load. This fact should be taken into account not only in the case of torsion but also in more complicated cases of loading __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 10, pp. 100–107, October 2005.     

Stability,bifurcation and post-critical behavior of a homogeneously deformed incompressible isotropic elastic parallelepiped subject to dead-load surface tractions

We study the equilibrium homogeneous deformations of a homogeneous parallelepiped made of an arbitrary incompressible, isotropic elastic material and subject to a distribution of dead-load surface tractions corresponding to an equibiaxial tensile stress state accompanied by an orthogonal uniaxial compression of the same amount. We show that only two classes of homogeneous equilibrium solutions are possible, namely symmetric deformations, characterized by two equal principal stretches, and asymmetric deformations, with all different principal stretches. Following the classical energy-stability criterion, we then find necessary and sufficient conditions for both symmetric and asymmetric equilibrium deformations to be weak relative minimizers of the total potential energy. Finally, we analyze the mechanical response of a parallelepiped made of an incompressible Mooney–Rivlin material in a monotonic dead loading process starting from the unloaded state. As a major result, we model the actual occurrence of a bifurcation from a primary branch of locally stable symmetric deformations to a secondary, post-critical branch of locally stable asymmetric solutions.     

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.     

In-plane stability of arches

Classical buckling theory is mostly used to investigate the in-plane stability of arches, which assumes that the pre-buckling behaviour is linear and that the effects of pre-buckling deformations on buckling can be ignored. However, the behaviour of shallow arches becomes non-linear and the deformations are substantial prior to buckling, so that their effects on the buckling of shallow arches need to be considered. Classical buckling theory which does not consider these effects cannot correctly predict the in-plane buckling load of shallow arches. This paper investigates the in-plane buckling of circular arches with an arbitrary cross-section and subjected to a radial load uniformly distributed around the arch axis. An energy method is used to establish both non-linear equilibrium equations and buckling equilibrium equations for shallow arches. Analytical solutions for the in-plane buckling loads of shallow arches subjected to this loading regime are obtained. Approximations to the symmetric buckling of shallow arches and formulae for the in-plane anti-symmetric bifurcation buckling load of non-shallow arches are proposed, and criteria that define shallow and non-shallow arches are also stated. Comparisons with finite element results demonstrate that the solutions and indeed approximations are accurate, and that classical buckling theory can correctly predict the in-plane anti-symmetric bifurcation buckling load of non-shallow arches, but overestimates the in-plane anti-symmetric bifurcation buckling load of shallow arches significantly.     

Three-dimensional oscillations of suspended cables involving simultaneous internal resonances

The near resonant response of suspended, elastic cables driven by planar excitation is investigated using a three degree-of-freedom model. The model captures the interaction of a symmetric in-plane mode with two out-of-plane modes. The modes are coupled through quadratic and cubic nonlinearities arising from nonlinear cable stretching. For particular magnitudes of equilibrium curvature, the natural frequency of the in-plane mode is simultaneously commensurable with the natural frequencies of the two out-of-plane modes in 1:1 and 2:1 ratios. A second nonlinear order perturbation analysis is used to determine the existence and stability of four classes of periodic solutions. The perturbation solutions are compared with results obtained by numerically integrating the equations of motion. Furthermore, numerical simulations demonstrate the existence of quasiperiodic responses.A portion of this work was presented at the 1992 ASME Winter Annual Meeting, Anaheim, CA.     

Steady state and stability of a beam on a damped tensionless foundation under a moving load

In the first part of this paper we study the effect of damping on the multiple steady state deformations of an infinite beam resting on a tensionless foundation and under a point load moving with a sub-critical speed. Due to the non-linear characteristics of the problem, a guess on the deformed shape has to be made before a numerical search can be initiated. It is found that when the damping is present, all the steady state solutions are asymmetric. As the damping approaches zero, some of the steady state solutions become symmetric, while some others remain asymmetric. In the second part of the paper we propose to test the stability of these steady state deformations by a transient analysis on a long finite beam. Our numerical experiment indicates that among all these multiple steady state solutions only one of them is stable. This stable steady state deformation reduces to a symmetric solution when the damping approaches zero. Furthermore, it is found that this stable solution is also the one among all steady state solutions closest in shape to the linear solution based on a bilateral foundation model.     

Writhing instabilities of twisted rods: from infinite to finite length

We use three different approaches to describe the static spatial configurations of a twisted rod as well as its stability during rigid loading experiments. The first approach considers the rod as infinite in length and predicts an instability causing a jump to self-contact at a certain point of the experiment. Semi-finite corrections, taken into account as a second approach, reveal some possible experiments in which the configuration of a very long rod will be stable through out. Finally, in a third approach, we consider a rod of real finite length and we show that another type of instability may occur, leading to possible hysteresis behavior. As we go from infinite to finite length, we compare the different information given by the three approaches on the possible equilibrium configurations of the rod and their stability. These finite size effects studied here in a 1D elasticity problem could help us guess what are the stability features of other more complicated (2D elastic shells for example) problems for which only the infinite length approach is understood.     

Study of static and dynamic stability of flexible rods in a geometrically nonlinear statement

We study static and dynamic stability problems for a thin flexible rod subjected to axial compression with the geometric nonlinearity explicitly taken into account. In the case of static action of a force, the critical load and the bending shapes of the rod were determined by Euler. Lavrent’ev and Ishlinsky discovered that, in the case of rod dynamic loading significantly greater than the Euler static critical load, there arise buckling modes with a large number of waves in the longitudinal direction. Lavrent’ev and Ishlinsky referred to the first loading threshold discovered by Euler as the static threshold, and the subsequent ones were called dynamic thresholds; they can be attained under impact loading if the pulse growth time is less than the system relaxation time. Later, the buckling mechanism in this case and the arising parametric resonance were studied in detail by Academician Morozov and his colleagues.In this paper, we complete and develop the approach to studying dynamic rod systems suggested by Morozov; in particular, we construct exact and approximate analytic solutions by using a system of special functions generalizing the Jacobi elliptic functions. We obtain approximate analytic solutions of the nonlinear dynamic problem of flexible rod deformation under longitudinal loading with regard to the boundary conditions and show that the analytic solution of static rod system stability problems in a geometrically nonlinear statement permits exactly determining all possible shapes of the bent rod and the complete system of buckling thresholds. The study of approximate analytic solutions of dynamic problems of nonlinear vibrations of rod systems loaded by lumped forces after buckling in the deformed state allows one to determine the vibration frequencies and then the parametric resonance thresholds.     

Wrinkling of a Stretched Thin Sheet

When a thin rectangular sheet is clamped along two opposing edges and stretched, its inability to accommodate the Poisson contraction near the clamps may lead to the formation of wrinkles with crests and troughs parallel to the axis of stretch. A variational model for this phenomenon is proposed. The relevant energy functional includes bending and membranal contributions, the latter depending explicitly on the applied stretch. Motivated by work of Cerda, Ravi-Chandar, and Mahadevan, the functional is minimized subject to a global kinematical constraint on the area of the mid-surface of the sheet. Analysis of a boundary-value problem for the ensuing Euler?CLagrange equation shows that wrinkled solutions exist only above a threshold of the applied stretch. A sequence of critical values of the applied stretch, each element of which corresponds to a discrete number of wrinkles, is determined. Whenever the applied stretch is sufficiently large to induce more than three wrinkles, previously proposed scaling relations for the wrinkle wavelength and, modulo a multiplicative factor that depends on the Poisson ratio of the sheet and the applied stretch and arises from the more general and weaker nature of geometric constraint under consideration, root-mean-square amplitude are confirmed. In contrast to the scaling relations for the wrinkle wavelength and amplitude, the applied stretch required to induce any number of wrinkles depends on the in-plane aspect ratio of the sheet. When the sheet is significantly longer than it is wide, the critical stretch scales with the fourth power of the length-to-width ratio but, when the sheet is significantly wider than it is long, the critical stretch scales with the square of that same ratio.     

Elastic Critical Buckling of Welded Thin Steel Rectangular Plates∗

A theoretical approach for the study of the effect of residual stresses due to welding on the elastic critical buckling behavior of thin steel rectangular plates is described. A finite difference technique is utilized for the determination of in-plane residual stresses due to a weld, and the Rayleight-Ritz method is used for the critical buckling problem, with stresses due to external loads being superimposed on the residual stresses. A number of illustrative examples are included, showing the possible detrimental effect of residual stresses due to welding. An approximately linear relationship is shown to exist between the square of the natural frequency of lateral vibration and the drop in buckling strength, for certain loading conditions     

Constitutive branching analysis of cylindrical bodies under in-plane equibiaxial dead-load tractions

Finite homogeneous deformations of hyperelastic cylindrical bodies subjected to in-plane equibiaxial dead-load tractions are analyzed. Four basic equilibrium problems are formulated considering incompressible and compressible isotropic bodies under plane stress and plane deformation condition. Depending on the form of the stored energy function, these plane problems, in addition to the obvious symmetric solutions, may admit asymmetric solutions. In other words, the body may assume an equilibrium configuration characterized by two unequal in-plane principal stretches corresponding to equal external forces. In this paper, a mathematical condition, in terms of the principal invariants, governing the global development of the asymmetric deformation branches is obtained and examined in detail with regard to different choices of the stored energy function. Moreover, explicit expressions for evaluating critical loads and bifurcation points are derived. With reference to neo-Hookean, Mooney-Rivlin and Ogden-Ball materials, a broad numerical analysis is performed and the qualitatively more interesting asymmetric equilibrium branches are shown. Finally, using the energy criterion, a number of considerations are put forward about the stability of the computed solutions.     

Semi-analytical method for computing effective properties in elastic composite under imperfect contact

A parallel fiber-reinforced periodic elastic composite is considered with transversely iso-tropic constituents. Fibers with circular cross section are distributed with the same periodicity along the two perpendicular directions to the fiber orientation, i.e., the periodic cell of the composite is square. The composite exhibits imperfect contact, in particular, spring type at the interface between the fiber and matrix is modeled. Effective properties of this composite for in-plane and anti-plane local problems are calculated by means of a semi-analytic method, i.e. the differential equations that described the local problems obtained by asymptotic homogenization method are solved using the finite element method. Numerical computations are implemented and comparisons with exact solutions are presented.     

Computations of modes I and II stress intensity factors of sharp notched plates under in-plane shear and bending loading by the fractal-like finite element method

The fractal-like finite element method (FFEM) is used to compute the stress intensity factors (SIFs) for different configurations of cracked/notched plates subject to in-plane shear and bending loading conditions. In the FFEM, the large number of unknown variables in the singular region around a notch tip is reduced to a small set of generalised co-ordinates by performing a fractal transformation using global interpolation functions. The use of exact analytical solutions of the displacement field around a notch tip as the global interpolation functions reduces the computational cost significantly and neither post-processing technique to extract SIFs nor special singular elements to model the singular region are required. The results of numerical examples of various configurations of cracked/notched plates are presented and validated via published data. Also, new results for cracked/notched plate problems are presented. These results demonstrate the accuracy and efficiency of the FFEM to compute the SIFs for notch problems under in-plane shear and bending loading conditions.     

钢管混凝土圆弧桁架拱平面内徐变稳定分析

核心混凝土的徐变会增加钢管混凝土拱肋的屈曲前变形,降低结构的稳定承载力,因此只有计入屈曲前变形的影响,才能准确得到钢管混凝土拱的徐变稳定承载力。基于圆弧形浅拱的非线性屈曲理论,采用虚功原理,建立了考虑徐变和剪切变形双重效应的管混凝土圆弧桁架拱的平面内非线性平衡方程,求得两铰和无铰桁架拱发生反对称分岔屈曲和对称跳跃屈曲的徐变稳定临界荷载。探讨了钢管混凝土桁架拱核心混凝土徐变随修正长细比、圆心角和加载龄期对该类结构弹性稳定承载力的影响,为钢管混凝土桁架拱长期设计提供理论依据。