FEM Analysis of the Three-Dimensional Buckling Problem for a Clamped Thick Rectangular Plate Made of a Viscoelastic Composite

The buckling instability of a thick rectangular plate made of a viscoelastic composite material is studied. The investigation is carried out within the framework of the three-dimensional linearized theory of stability. The plate edges are clamped and the plate is compressed through the clamps. Moreover, it is assumed that the plate has an initial infinitesimal imperfection, and, as a buckling criterion, the state is taken where this imperfection starts to increase indefinitely at fixed finite values of external compressive forces. From this criterion, the critical time is determined. The corresponding boundary-value problems are solved by employing the three-dimensional FEM and the Laplace transform. The material of the plate is assumed orthotropic, viscoelastic, and homogeneous. Numerical results related to the critical time are presented.     

Free vibration analysis of cracked postbuckled plate

In this paper, a methodology is introduced to address the free vibration analysis of cracked plate subjected to a uniaxial inplane compressive load for the first time. The crack, assumed to be open and at the edge is modeled by a massless linear rotational spring. The governing differential equations are derived using the Mindlin theory, taking into account the effect of initial imperfection. The response is assumed to be consisting of static and dynamic parts. For the static part, differential equations are discretized using the differential quadrature element method and resulting nonlinear algebraic equations are solved by an arc-length strategy. Assuming small amplitude vibrations of the plate about its buckled state and exploiting the static solution in the linearized vibration equations, the dynamic equations are converted into a non-standard eigenvalue problem. Finally, natural frequencies and modal shapes of the cracked buckled plate are obtained by solving this eigenvalue problem. To ensure the validity of the suggested approach an experimental setup and a numerical finite element model have been made to analyze the vibration of a cracked square plate with simply supported boundary conditions. Also, several case-studies of cracked buckled plate problem have been solved utilizing the proposed method, and effects of selected parameters have been studied. The results show that the applied load and geometric imperfection as well as the position, size and depth of the crack have different impact on natural frequencies of the plate.     

Fiber Buckling in a Viscoelastic Matrix

Within the framework of the three-dimensional linearized theory of stability, an approach for investigating fiber buckling in the structure of unidirectional fibrous viscoelastic composites is developed. For simplicity, a small fiber concentration is considered, and the buckling problem for a single elastic fiber in an infinite viscoelastic matrix is investigated. In this case, it is assumed that the fiber has an insignificant initial periodical imperfection, and the growth of this imperfection with time is studied. The state where this imperfection starts to grow indefinitely is taken as a fiber-buckling criterion, and the critical time is determined from this criterion.     

Subcritical and critical states of a crack with failure zones

The problem on the stress–strain state near a mode I crack in an infinite plate is solved in the frame of a cohesive zone model. The complex variable method of Muskhelishvili is used to obtain the crack opening displacements caused by the cohesive traction, which models the failure zone at the crack tip, as well as by the external load. The finite stress condition and logarithmic singularity of the derivative of the separation with respect to the coordinate at the tip of a physical crack are taken into account.The cohesive traction distribution is sought in a piecewise linear form, nodal values of which are being numerically chosen to satisfy the traction-separation law. According to this law, the cohesive traction is coupled with the corresponding separation and fracture toughness. The tips of the physical crack and cohesive zone (geometric variables) along with the discrete cohesive traction are used as the problem parameters determining the stress-strain state. If the crack length is included in the set, then the critical crack size can be found for the given loading intensity.The obtained determining system of equations is solved numerically. To find the initial point for a standard numerical algorithm, the asymptotic determining system is derived. In this system, the geometric variables can be easily eliminated, which make it possible to linearize the system.In the numerical examples, the one-parameter traction-separation laws are used. Influence of the shape parameters of the law on the critical crack size and the corresponding cohesive length is studied. The possibility of using asymptotic solutions for determining the critical parameters is analysed. It is established that the critical crack length slightly depends on the shape parameter, while the cohesive length shows a strong dependence on the shape of cohesive laws.     

Local Buckling around an Interfacial Crack in a Viscoelastic Sandwich Plate

Buckling around an interfacial crack in a clamped sandwich plate is studied. The layers of the plate are assumed to be linearly viscoelastic, isotropic, and homogeneous. The investigations are carried out within the framework of a piecewise homogeneous body model with the use of a three-dimensional linearized theory of stability. The corresponding boundary value problems are solved numerically by employing the FEM and the Laplace transform. The calculated critical times are presented for various problem parameters. In this case, the upper and lower layers are assumed to be viscoelastic and are described by Rabotnov operators, whereas the midlayer is regarded as purely elastic. The influence of rheological parameters on the critical time is also analyzed.     

Near-surface buckling in stability of a system consisting of a moderately rigid substrate, a viscoelastic bond layer, and an elastic covering layer

Within the frame work of the three-dimensional linearized theory of stability of deformable bodies (TLTSDB), the near-surface buckling instability of a system consisting of a half-plane (substrate), a viscoelastic bond layer, and an elastic covering layer is suggested. The equations of the TLTSDB are obtained from the three-dimensional geometrically non linear equations of viscoelasticity theory by using the boundary-form perturbation technique. By employing the Laplace transform, a method for solving the problem is developed. It is supposed that the covering layer has an insignificant initial imperfection. The stability of the system is considered lost if the imperfection starts to increase and grows indefinitely. Numerical results for the critical compressive force and the critical time are presented. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 42, No. 4, pp. 517–530, July–August, 2006.     

Local near-surface buckling of a system consisting of an elastic (viscoelastic) substrate,a viscoelastic (elastic) bond layer,and an elastic (viscoelastic) covering layer

Within the framework of a piecewise homogenous body model and with the use of a three-dimensional linearized theory of stability (TLTS), the local near-surface buckling of a material system consisting of a viscoelastic (elastic) half-plane, an elastic (viscoelastic) bond layer, and a viscoelastic (elastic) covering layer is investigated. A plane-strain state is considered, and it is assumed that the near-surface buckling instability is caused by the evolution of a local initial curving (imperfection) of the elastic layer with time or with an external compressive force at fixed instants of time. The equations of TLTS are obtained from the three-dimensional geometrically nonlinear equations of the theory of viscoelasticity by using the boundary-form perturbation technique. A method for solving the problems considered by employing the Laplace and Fourier transformations is developed. It is supposed that the aforementioned elastic layer has an insignificant initial local imperfection, and the stability is lost if this imperfection starts to grow infinitely. Numerical results on the critical compressive force and the critical time are presented. The influence of rheological parameters of the viscoelastic materials on the critical time is investigated. The viscoelasticity of the materials is described by the Rabotnov fractional-exponential operator. Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 43, No. 6, pp. 771–788, November–December, 2007.     

初始缺陷和拉-弯耦合对于叠层板的振动、屈曲和非线性动力稳定性的影响

研究了复合材料叠层板的初始缺陷和拉伸＿弯曲耦合对于其振动、屈曲和非线性动力稳定性的影响·推导了控制方程·这是一个修正的非线性Ｍａｔｈｉｅｕ方程·进行了５种典型复合材料的数值计算,它们是玻璃环氧Ｓｃｏｔｃｈ＿１００２,芳纶环氧Ｋｅｖｌａｒ＿４９,硼环氧Ｂ４／５５０５,石墨环氧Ｔ３００／５２０８,和ＡＳ／３５０１·结果表明：由于初始缺陷耦合效应的存在,使叠层板对于进入参数共振更加敏感,并且其振幅大于无初始缺陷或者无耦合效应的叠层板·对于不同复合材料的叠层板,尤其是层数较少的板,耦合效应是不相同的·在板结构的屈曲和动力稳定性设计中,如果忽略了耦合效应的影响,其不安全性将超过１０％以上     

Compressive and thermal postbuckling behaviors of laminated plates with piezoelectric fiber reinforced composite actuators

This paper studied compressive postbuckling under thermal environments and thermal postbuckling due to a uniform temperature rise for a shear deformable laminated plate with piezoelectric fiber reinforced composite (PFRC) actuators based on a higher order shear deformation plate theory that includes thermo-piezoelectric effects. The material properties are assumed to be temperature-dependent and the initial geometric imperfection of the plate is considered. The compressive and thermal postbuckling behaviors of perfect, imperfect, symmetric cross-ply and antisymmetric angle-ply laminated plates with fully covered or embedded PFRC actuators are conducted under different sets of thermal and electric loading conditions. The results reveal that, the applied voltage usually has a small effect on the postbuckling load–deflection relationship of the plate with PFRC actuators in the compressive buckling case, whereas the effect of applied voltage is more pronounced for the plate with PFRC actuators, compared to the results of the same plate with monolithic piezoelectric actuators.     

具有初始缺陷弹性板的稳定性分析

本文以Marguerre方程为基础,用奇异性理论研究了初始挠度缺陷以及横向载荷对弹性板屈曲后分叉解的影响。借助于普适开折的原理,在单特征值局部邻域内将该问题的失稳分析转化为三次代数方程的讨论,从而确定出分叉解的性态。同时绘出了在不同参数下的分叉解文,讨论了几何缺陷和横向载荷对特征值的影响。     

The problem of an interface crack with a rigid patch plate on part of its edge

A mixed boundary-value problem is solved for a piecewise-homogeneous elastic body with a rectilinear semi-infinite crack on the line where the materials are joined. A rigid patch plate (a reinforcing plate) of specified shape is attached to the upper edge of the crack on a finite interval adjacent to the crack tip. The edges of the crack are loaded with specified stresses. The body is stretched at infinity by a specified longitudinal stress. External forces with a given principal vector and moment act on the patch plate. The problem reduces to a Riemann-Hilbert boundary-value matrix problem with a piecewise-constant coefficient, the solution of which is explicitly constructed using a Gaussian hypergeometric function. The angle of rotation of the patch plate and the complex potentials describing the stress state of the body are found and the stress state of the body close to the ends of the patch plate, one of which is also simultaneously the crack tip, is investigated. Numerical examples are presented that illustrate the effect of the initial force parameters, the length of the patch plate and other parameters of the body on the angle of rotation of the patch plate and the stress state of the body.     

The inverse bending problem for a thin plate with an elastic imperfection

We consider the inverse problem consisting of determining the unknown shape of an elastic imperfection contained in a thin plate from the condition of equal strength in the stressed state along the phase interface surface. It is shown that such a state is attained in the case of an elliptic imperfection whose shape depends on the values of the applied moments and the mechanical properties of the component phases. It is established that for the geometry found for the imperfection the sum of the moments is constant and the second invariant of the deviator of the stress tensor is superharmonic over the entire plate. Numerical computations are carried out. In special cases the results obtained coincide with known data. One figure. Bibliography: 5 titles. Translated fromTeoreticheskaya i Prikladnaya Mekhanika, No. 22, pp. 34–40, 1991.     

Nonlinear vibration and instability of functionally graded nanopipes with initial imperfection conveying fluid

In this paper, the nonlinear vibration and instability of a fluid-conveying nanopipe made of functionally graded (FG) materials with consideration of the initial geometric imperfection are investigated. The material properties are assumed to vary smoothly along the radial direction according to a power-law exponent form. The fluid-conveying FG nanopipe is modeled as a Euler-Bernoulli beam, and the governing equation is derived based on the nonlocal strain gradient theory incorporating the effects of Von-Karman geometrical nonlinearity and initial imperfection. The nonlinear frequency and critical fluid velocity are achieved via He's Hamiltonian approach. After verifying the present model with comparison of several previous studies, the effect of several different system parameters including the amplitude of the nonlinear oscillator, the initial geometric imperfection, size-dependent parameters, and the power-law index on the frequency response of the fluid-conveying FG nanopipe are explored. Moreover, the critical velocity of the conveying fluid under different system parameters is also investigated and discussed in detail. The developed size-dependent nonlinear model is expected to provide a possible theoretical way to guide the application of FG nanopipe as micro/nanofluidic devices.     

Imperfection Sensitivity or Insensitivity of Zero-stiffness Postbuckling … that is the Question

Zero-stiffness postbuckling of a structure is characterized by a secondary load-displacement path along which the load remains constant. In sensitivity analysis of the (initial) postbuckling path it is usually considered as a borderline case between imperfection sensitivity and imperfection insensitivity. However, it is unclear whether zero-stiffness postbuckling as such is imperfection sensitive or insensitive. In this paper, Koiter's initial postbuckling analysis is used as a tool for sensitivity analysis. Distinction between two kinds of imperfections is made on the basis of the behavior of the equilibrium path of the imperfect structure. New definitions of imperfection insensitivity of the postbuckling behavior are provided according to the classification of imperfections. A structure with two degrees of freedom with a zero-stiffness postbuckling path is studied, considering four different imperfections. The results from this example show that zero-stiffness postbuckling is a case of transition from imperfection sensitivity to imperfection insensitivity for imperfections of the first kind and that it is imperfection insensitive for imperfections of the second kind. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stress intensity factors for a finite plate with an inclined crack by boundary collocation

In this paper, we combine the Muskhelishvili＇s complex variable method and boundary collocation method, and choose a set of new stress function based on the stress boundary condition of crack surface, the higher precision and less computation are reached. This method is applied to calculating the stress intensity factor for a finite plate with an inclined crack. The influence of θ （the obliquity of crack） on the stress intensity factors, as well as the number of summation terms on the stress intensity factor are studied and graphically represented.     

Invariant Integrals for the Equilibrium Problem for a Plate with a Crack

We consider the equilibrium problem for a plate with a crack. The equilibrium of a plate is described by the biharmonic equation. Stress free boundary conditions are given on the crack faces. We introduce a perturbation of the domain in order to obtain an invariant Cherepanov–Rice-type integral which gives the energy release rate upon the quasistatic growth of a crack. We obtain a formula for the derivative of the energy functional with respect to the perturbation parameter which is useful in forecasting the development of a crack (for example, in study of local stability of a crack). The derivative of the energy functional is representable as an invariant integral along a sufficiently smooth closed contour. We construct some invariant integrals for the particular perturbations of a domain: translation of the whole cut and local translation along the cut.     

平面弹性裂纹分析的一种有效边界元方法

提出了一种简单而有效的平面弹性裂纹应力强度因子的边界元计算方法．该方法由Crouch与Starfield建立的常位移不连续单元和闫相桥最近提出的裂尖位移不连续单元构成A·D2在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界．算例(如单向拉伸无限大板中心裂纹、单向拉伸无限大板中圆孔与裂纹的作用)说明平面弹性裂纹应力强度因子的边界元计算方法是非常有效的．此外,还对双轴载荷作用下有限大板中方孔分支裂纹进行了分析．这一数值结果说明平面弹性裂纹应力强度因子的边界元计算方法对有限体中复杂裂纹的有效性,可以揭示双轴载荷及裂纹体几何对应力强度因子的影响．     

Influence of Initial Stresses on Stress Intensity Factors at Crack Tips in a Composite Strip

The influence of initial tension or compression along cracks on the stress intensity factor (SIF) at crack tips under the action of additional normal forces on crack edges is studied for infinite bodies. A strip made of a composite material is considered. The strip ends are simply supported, and the strip contains a crack whose edges are parallel to its face planes. The strip is first stretched or compressed along crack edges, and then additional uniformly distributed normal forces are applied to the crack edges. The influence of the initial tension (compression) on the SIF caused by the additional normal forces is studied. The corresponding boundary-value problems are modelled with the use of the three-dimensional linearized theory of elasticity. All the investigations are carried out numerically by employing the finite-element method. The values of SIF are calculated by the energy release method.     

Extension of a piezoceramic bimorph with a crack crossing the interface

We consider the boundary-value problem of electroelasticity for a composite plate weakened by a crack crossing the line joining the media. The initial boundary-value problem, is reduced to a mixed system of singular integral and algebraic equations. We present the calculation results characterizing the variation in the stress intensity factors as a function of the opening angle of a segmented crack for different types of loading.Translated from Mekhanika Kompozitnykh Materialov, Vol. 33, No. 4, pp. 482–488, July–August 1997.