Stability of rod-shaped nanoparticles embedded in an elastic matrix

In many biological tissues as well as in some technical materials we find nano-sized rod-shaped particles embedded in a relatively soft matrix. Loss of stability of equilibrium, i.e. buckling, is one of the possible failure modes of such materials. In the present paper different kinds of load transfer between matrix and reinforcing particles, which are typical for rod-shaped nanostructures in biological tissues, are considered with respect to stability of equilibrium. Two regimes of matrix stiffnesses leading to different modes of buckling, and a transition regime in between, have been found: soft matrix materials leading to the so-called ‘flip mode’ (also called ‘tilt mode’) and hard matrix materials resulting in ‘bending mode’ buckling. The transition regime is of particular interest for biological tissues. Numerical and semi-analytical as well as asymptotic concepts are employed leading to results for estimating the critical load intensities both in the form of closed form solutions and diagrams. The analytical solutions are compared with results of finite element analyses. From these comparisons indications are gained for deciding which of the different analytical approaches should be chosen for a particular nanostructure configuration in terms of the associated buckling modes.     

再入过载的匹配渐进展开分析与卸载方法研究

探月飞船返回地球时将以第二宇宙速度再入地球大气层,面临极其严苛的气动环境,因此对于再入气动过载的分析具有重要意义.再入运动方程是一组非线性很强的常微分方程,数值方法计算量大,不适用于在线任务.因此,本文采用一种近似解法对气动过载进行分析.首先,基于匹配渐进展开方法将再入纵向运动解在大气外层区域与内层区域分别展开,得到统一形式的闭型近似解,在此基础上分段求解气动过载,并与精确解进行对比分析.其次,利用闭型近似解,通过当前状态反解虚拟初始条件,在此基础上提出初次过载峰值的解析预测方法,并分析了不同条件下预测的相对误差变化规律.最后,基于过载峰值的解析预测对飞船的初次再入过程进行卸载,将飞船在再入过程中耗散的总能量进行重新分配,并通过蒙特卡罗飞行仿真试验验证了卸载方法的有效性.     

Nonlinear dynamic buckling of shallow circular arches under a sudden uniform radial load

The structural behavior of a shallow arch is highly nonlinear, and so when the amplitude of the oscillation of the arch produced by a suddenly applied load is sufficiently large, the oscillation of the arch may reach a position at its primary unstable equilibrium path or secondary bifurcation unstable equilibrium path, leading the arch to buckle dynamically. This paper presents an analytical study of the nonlinear dynamic in-plane buckling of a shallow circular arch under a uniform radial load that is applied suddenly and with an infinite duration. The principle of conservation of energy is used to establish the criterion for dynamic buckling of the arch, and the analytical solution for the dynamic buckling load is derived. It is shown that under a suddenly applied uniform radial load, a shallow pinned–fixed arch has a unique possible dynamic buckling load, while shallow pinned–pinned and fixed–fixed arches may have two possible dynamic buckling loads: a lower dynamic buckling load and an upper dynamic buckling load. The dynamic buckling loads of a shallow arch under a suddenly applied uniform radial load with infinite duration are found to be lower than their static counterparts, and to increase with an increase of the arch included angle and slenderness. The effect of static preloading on the dynamic buckling of an arch is also investigated. It is found that the pre-applied static load decreases the dynamic buckling load of the arch, but increases the sum of the pre-applied load and the dynamic buckling load.     

Buckling analysis of nanobeams with exponentially varying stiffness by differential quadrature method

We present the application of differential quadrature(DQ) method for the buckling analysis of nanobeams with exponentially varying stiffness based on four different beam theories of Euler-Bernoulli, Timoshenko, Reddy, and Levison.The formulation is based on the nonlocal elasticity theory of Eringen. New results are presented for the guided and simply supported guided boundary conditions. Numerical results are obtained to investigate the effects of the nonlocal parameter,length-to-height ratio, boundary condition, and nonuniform parameter on the critical buckling load parameter. It is observed that the critical buckling load decreases with increase in the nonlocal parameter while the critical buckling load parameter increases with increase in the length-to-height ratio.     

轴压和扭转复合载荷作用下氮化硼纳米管屈曲行为的分子动力学模拟

采用分子动力学方法模拟了氮化硼纳米管在轴压和扭转复合荷载作用下的屈曲和后屈曲行为.在各加载比例下,给出了初始线性变形阶段和后屈曲阶段原子间相互作用力的变化,确定了屈曲临界荷载关系.通过对屈曲模态的细致研究,从微观变形机理上分析了纳米管对不同外荷载力学响应的差异.研究结果表明,扶手型和锯齿型纳米管均呈现出非线性的屈曲临界荷载关系,复合加载下的屈曲行为具有强烈的尺寸依赖性.温度升高将导致屈曲临界荷载的下降,且温度的影响随加载比例的变化而变化.无论在简单加载或复合加载中,同尺寸的碳纳米管均比氮化硼纳米管具有更强地抵抗屈曲荷载的能力.     

Dynamic column buckling of multi-walled carbon nanotubes under axial impact load

This paper studies the dynamic column buckling of multi-walled carbon nanotubes (MWNTs) under axial impact load. The analysis is based on the continuum mechanics model and a simplified model for the van der Waals forces between adjacent layers. By introducing initial imperfections for MWNTs and applying the method of preferred mode, a buckling condition is derived for the buckling load and associated buckling mode. In particular, explicit expressions are obtained for double-walled carbon nanotubes (DWNTs). Finally, numerical calculations are worked out for a DWNT and a five-layer MWNT with different length-to-radius ratios.     

Load-response determination for imperfect column using vibratory data

In addition to the initial crook to a column, the structural integrity is often not as reliable as expected due to the immaturity of current engineering techniques. Thus, the actual load response is sometimes not consistent with the design predictions for the column. For design considerations, it is necessary to establish experimentally an analytical method for determining the load-deflection curves and buckling load. In this paper, a dynamic method is described for building the load-deflection curves up to its buckling load for an imperfect column. The proposed method does not require the application of axial load and is feasible for arbitrary types of boundary conditions. The load-deflection curves and buckling load are determined from the measured natural frequencies and the vibration modes.     

Multi-objective optimization of ring stiffened cylindrical shells using a genetic algorithm

In this paper, the genetic algorithm (GA) method is used for the multi-objective optimization of ring stiffened cylindrical shells. The objective functions seek the maximum fundamental frequency and minimum structural weight of the shell subjected to four constraints including the fundamental frequency, the structural weight, the axial buckling load, and the radial buckling load. The optimization process contains six design variables including the shell thickness, the number of stiffeners, the width and height of stiffeners, the stiffeners eccentricity distribution order, and the stiffeners spacing distribution order. The real coding scheme is used for representing the solution string, while the generation number-based adaptive penalty function is applied for penalizing infeasible solutions. In analytical solution, the Ritz method is applied and the stiffeners are treated as discrete elements. Some examples of simply supported cylindrical shells with nonuniform eccentricity distribution and nonuniform rings spacing distribution are provided to demonstrate the optimality of the solution obtained by the GA technique. The effects of objective weighting coefficients and bounding values of the design variables on the optimum solution are studied for various cases. The results show that the optimal solution can vary with the weighting coefficients significantly. It is also found that extreme reduction and augmentation in turn in the structural weight and fundamental frequency can be simultaneously achieved by selecting suitable stiffeners’ geometrical parameters and distributions. Furthermore, the bounding values of the design variables have great effects on the optimum results.     

弧长法在斜直井内管柱非线性屈曲分析中的应用

对于斜直井内底部-段管柱的后屈曲问题,基于受径向约束管柱的微分求积(DQ,Differential Quadrature)单元,构建了弧长迭代法.给出详细的迭代步骤和迭代初值的确定方法,对不同端部侧向约束条件下的管柱非线性屈曲进行迭代计算.并与现有文献中的近似解析解、实验结果和纯载荷增量迭代法的数值计算结果进行比较.结果显示,本文方法克服了有限单元法在处理管柱自重时的困难,同时能自动调节增量步长,跟踪管柱非线性后屈曲平衡路径的全过程.计算效率高、收敛性好、易于实施,可以用来分析斜直井内管柱的非线性屈曲问题.     

Vibration and buckling of a double-beam system under compressive axial loading

On the basis of the Bernoulli–Euler beam theory, the properties of free transverse vibration and buckling of a double-beam system under compressive axial loading are investigated in this paper. It is assumed that the two beams of the system are simply supported and continuously joined by a Winkler elastic layer. Explicit expressions are derived for the natural frequencies and the associated amplitude ratios of the two beams, and the analytical solution of the critical buckling load is obtained. The influences of the compressive axial loading on the responses of the double-beam system are discussed. It is shown that the critical buckling load of the system is related to the axial compression ratio of the two beams and the Winkler elastic layer, and the properties of free transverse vibration of the system greatly depend on the axial compressions.     

轴向冲击载荷作用下双壁碳纳米管的动力屈曲

应用改进的有限元方法,建立考虑层间范德华力作用的壳-弹簧非线性有限元模型,基于B-R运动准则,系统地研究了双壁碳纳米管的动力屈曲问题,得到了轴向冲击载荷作用下双壁碳纳米管的临界动力屈曲载荷和临界动力失效载荷. 研究结果表明,在动力屈曲过程中,双壁碳纳米管层间距的变化非常小,各管的变形相互协调;碳纳米管中应力波的传播导致碳纳米管出现非对称屈曲模态,可明显观测到四个环向波瓣,沿着碳纳米管的轴线方向,四个波瓣的波峰和波谷交替变化. 对碳纳米管动力屈曲问题的研究表明,冲击载荷的大小和持续时间对碳纳米管的动力屈曲有 关键词： 碳纳米管 动力屈曲 冲击载荷     

Buckling of 2D-FG Cylindrical Shells under Combined External Pressure and Axial Compression

This paper presents the stability of two-dimensional functionally graded (2D-FG) cylindrical shells subjected to combined external pressure and axial compression loads, based on classical shell theory. The material properties of functionally graded cylindrical shell are graded in two directional (radial and axial) and determined by the rule of mixture. The Euler's equation is employed to derive the stability equations, which are solved by GDQ method to obtain the critical mechanical buckling loads of the 2D-FG cylindrical shells. The effects of shell geometry, the mechanical properties distribution in radial and axial direction on the critical buckling load are studied and compared with a cylindrical shell made of 1D-FGM. The numerical results reveal that the 2D-FGM has a significant effect on the critical buckling load.     

带缺口加强圈的圆柱壳屈曲特性分析

以低温液体罐车外筒体为原型,采用有限元方法研究了带缺口加强圈的圆柱壳的屈曲特性。着重考察了加强圈缺口间的夹角、加强圈之间的间距、边界条件对此类结构屈曲载荷的影响。计算结果表明:圆柱壳的屈曲载荷与加强圈缺口间夹角不是简单的线性关系,当加强圈缺口间夹角增大到一定程度后,圆柱壳的屈曲载荷几乎不再变化;加强圈布置的均匀度不但会影响到圆柱壳的屈曲载荷,同时也对其屈曲模态产生影响;在不同的边界条件下,圆柱壳也表现出不同的屈曲特性。     

Vibration and stability of rectangular strip-plates

The problem considered is that of free vibration and stability of a simply supported rectangular strip-plate subjected to constant in-plane forces. The relevant continuity conditions at the interface between the adjacent regions, which play a significant role in this type of problem, are derived. The eigenfrequencies and the buckling load are estimated by the method of the new quotient which is based on a variational statement proposed by Nemat-Nasser. The results are compared with those obtained by means of the usual Rayleigh quotient and the exact solution. The good accuracy obtained by the application of the method of the new quotient is demonstrated by means of numerical examples.     

On the nonlinear axisymmetric dynamic buckling behavior of clamped functionally graded spherical caps

Here, the dynamic thermal buckling behavior of functionally graded spherical caps is studied considering geometric nonlinearity based on von Karman's assumptions. The formulation is based on first-order shear deformation theory and it includes the in-plane and rotary inertia effects. The material properties are graded in the thickness direction according to the power-law distribution in terms of volume fractions of the material constituents. The effective material properties are evaluated using homogenization method. The governing equations obtained using finite element approach are solved employing the Newmark's integration technique coupled with a modified Newton–Raphson iteration scheme. The pressure load corresponding to a sudden jump in the maximum average displacement in the time history of the shell structure is taken as the dynamic buckling load. The present model is validated against the available isotropic case. A detailed numerical study is carried out to highlight the influences of shell geometries, power law index of functional graded material and boundary conditions on the dynamic buckling load of shallow spherical shells.     

Exact closed-form solutions for the natural frequencies and stability of elastically connected multiple beam system using Timoshenko and high-order shear deformation theory

A general procedure for the determination of the natural frequencies and buckling load for a set of beam system under compressive axial loading is investigated using Timoshenko and high-order shear deformation theory. It is assumed that the set beams of the system are simply supported and continuously joined by a Winkler elastic layer. The model of beam includes the effects of axial loading, shear deformation and rotary inertia. In the special case of identical beams, explicit expressions for the natural frequencies and the critical buckling load are determined using a trigonometric method. The influences of the compressive axial loading and the number of beams in the system on the natural frequencies and the critical buckling load are discussed. These results are of considerable practical interest and have wide application in engineering practice of frameworks.     

Prediction of compressive post-buckling behavior of single-walled carbon nanotubes in thermal environments

In the present investigation, the axial buckling and post-buckling configurations of single-walled carbon nanotubes (SWCNTs) are studied including the thermal environment effect. For this purpose, Eringen’s nonlocal elasticity continuum theory is implemented into the classical Euler–Bernoulli beam theory to represent the SWCNTs as a nonlocal elastic beam model. A closed-form analytical solution is carried out to analyze the static response of SWCNTs in their post-buckling state in which the axial buckling load is assumed to be beyond the critical axial buckling load. Common sets of boundary conditions, named simply supported–simply supported (SS–SS), clamped–clamped (C–C), and clamped–simply supported (C–SS), are considered in the investigation. Selected numerical results are given to represent the variation of the carbon nanotube’s mid-span deflection with the applied axial load corresponding to various nonlocal parameters, length-to-diameter aspect ratios, temperature changes, and end supports. Moreover, a comparison between the post-buckling behaviors of SWCNTs at low- and high-temperature environments is presented. It is found that the size effect leads to a decrease of the axial buckling load especially for SWCNTs with C–C boundary conditions. Also, it is revealed that the value of the temperature change plays different roles in the post-buckling response of SWCNTs at low- and high-temperature environments.     

An asymptotic outer solution applied to the Keller box method

The Keller box method (“Numerical Solutions of Partial Differential Equations, Vol. 2” (B. Hubbard Ed.), pp. 327–350, Academic Press, New York, 1970) was applied to incompressible flow past a flat plate to demonstrate that the basic computation region must extend outward from the wall until the outer boundary conditions are effectively obtained. The Keller box method was modified to include an asymptotic outer solution for the case of the self-similar solution for compressible flow in a boundary layer. Initial application of the basic and modified Keller box methods to incompressible flow past a flat plate showed similar rates of convergence but smaller RMS error for the same basic range of the independent variable when the asymptotic outer solution is applied. Furthermore, extension of the solution beyond the range of the independent variable for the numerical solution using the resulting asymptotic solution produced RMS error at least as small as the RMS error on the range of the numerical solution. Also, when the asymptotic solution was applied, a smaller range of independent variables could be used in the numerical solution to obtain the same RMS error. Numerical results for compressible flow were qualitatively the same as for the case with the incompressible velocity profile except the rate of iterative convergence was slightly slower. Application of asymptotic outer solution for incompressible flow at a two dimensional stagnation point produced similar results with smaller relative improvements. For compressible flow with smaller favorable pressure gradients than the stagnation point and with adverse pressure gradients, significant improvements were again obtained. Examination of the errors associated with the asymptotic solution reveals that greatest success is obtained for flows with thicker boundary layers and shows that the boundary layer at a two dimensional stagnation point is too thin for small error in the asymptotic solution. Despite relatively large errors in the asymptotic solutions for boundary layer in strong favorable pressure gradients where the boundary layer is thin, the boundary layer solutions generally showed improvement in error and reduction in computation times.     

Column buckling of doubly parallel slender nanowires carrying electric current acted upon by a magnetic field

Axial buckling of current-carrying double-nanowire-systems immersed in a longitudinal magnetic field is aimed to be explored. Each nanowire is affected by the magnetic forces resulted from the externally exerted magnetic field plus the magnetic field resulted from the passage of electric current through the adjacent nanowire. To study the problem, these forces are appropriately evaluated in terms of transverse displacements. Subsequently, the governing equations of the nanosystem are constructed using Euler–Bernoulli beam theory in conjunction with the surface elasticity theory of Gurtin and Murdoch. Using a meshless technique and assumed mode method, the critical compressive buckling load of the nanosystem is determined. In a special case, the obtained results by these two numerical methods are successfully checked. The roles of the slenderness ratio, electric current, magnetic field strength, and interwire distance on the axial buckling load and stability behavior of the nanosystem are displayed and discussed in some detail.     

On the durable critic load in creep buckling of viscoelastic laminated plates and circular cylindrical shells

Based on the first order shear deformation theory and classic buckling theory, the paper investigates the creep buckling behavior of viscoelastic laminated plates and laminated circular cylindrical shells. The analysis and elaboration of both instantaneous elastic critic load and durable critic load are emphasized. The buckling load in phase domain is obtained from governing equations by applying Laplace transform, and the instantaneous elastic critic load and durable critic load are determined according to the extreme value theorem for inverse Laplace transform. It is shown that viscoelastic approach and quasi-elastic approach yield identical solutions for these two types of critic load respectively. A transverse disturbance model is developed to give the same mechanics significance of durable critic load as that of elastic critic load. Two types of critic loads of boron/epoxy composite laminated plates and circular cylindrical shells are discussed in detail individually, and the influencing factors to induce creep buckling are revealed by examining the viscoelasticity incorporated in transverse shear deformation and in-plane flexibility.