基于Lenosky原子作用势单层石墨烯片的力学模型

基于对Lenosky 碳-碳共价键作用势连续化得到的单层石墨烯的势能和Hamilton 原理,导出了单层石墨烯的动力学方程. 使用该数学模型及Galerkin 方法,研究了矩形单层石墨烯片的静力挠曲问题. 结果显示,石墨烯片的几何尺寸较小时,弯曲刚度对结构的受力影响较大,可用板理论来描述;随着结构尺寸的增大,弯曲刚度的影响迅速降低;当矩形石墨烯片的短边尺寸大于10 nm 时,可以忽略弯曲刚度对结构的影响,使用薄膜理论来描述单层石墨烯的力学性质.      

Nonlinear softening and hardening nonlocal bending stiffness of an initially curved monolayer graphene

In the present article, the governing nonlinear nonlocal elastic equations are obtained for a monolayer graphene with an initial curvature and the related softening and hardening bending stiffness is analytically calculated. The effects of large deformation, initial curvature, discreteness and direction of chiral vector on the bending stiffness of the monolayer graphene are discussed in detail. A behavior more complex than previously reported in the literature emerges. It is found that the bending stiffness of graphene strongly depends on the initial configuration, showing not obvious maxima and minima, and suggesting the possibility of a smart tuning.     

INVESTIGATIONS ON THE MECHANICAL CHARACTERISTICS OF NANOBEAMS BASED ON THE NONLOCAL STRAIN GRADIENT THEORY1)

基于非局部应变梯度理论,建立了一种具有尺度效应的高阶剪切变形纳米梁的力学模型. 其中,考虑了应变场和一阶应变梯度场下的非局部效应. 采用哈密顿原理推导了纳米梁的控制方程和边界条件,并给出了简支边界条件下静弯曲、自由振动和线性屈曲问题的纳维级数解. 数值结果表明,非局部效应对梁的刚度产生软化作用,应变梯度效应对纳米梁的刚度产生硬化作用,梁的刚度整体呈现软化还是硬化效应依赖于非局部参数与材料特征尺度的比值. 梁的厚度与材料特征尺度越接近,非局部应变梯度理论与经典弹性理论所预测结果之间的差异越显著.     

THIRD-ORDER GENERALISED BEAM THEORY AND CALCULATION METHOD

一阶广义梁理论描述通过运用加入弯曲、扭转和畸变函数的普通非耦合微分方程组解决棱柱状结构行为.二阶广义梁理论,是添加上偏离力效果的微分方程. 通过引入纵向膜弯矩和膜剪应变虚功到广义梁理论系统当中,完全展开的三阶广义梁方程组将以一串大型离散迭代函数且能转化为可用于数值分析的若干切线刚度矩阵形式出现. 通过膜应力派生出三阶分项^ijrkv_σ和^ijrkv_τ并结合先进数值技术寻求全解,三阶广义梁理论提供了一种严谨和高效的数值工具用于调查薄壁结构后屈曲大变形行为.     

A size-dependent Kirchhoff micro-plate model based on strain gradient elasticity theory

A size-dependent Kirchhoff micro-plate model is developed based on the strain gradient elasticity theory. The model contains three material length scale parameters, which may effectively capture the size effect. The model can also degenerate into the modified couple stress plate model or the classical plate model, if two or all of the material length scale parameters are taken to be zero. The static bending, instability and free vibration problems of a rectangular micro-plate with all edges simple supported are carried out to illustrate the applicability of the present size-dependent model. The results are compared with the reduced models. The present model can predict prominent size-dependent normalized stiffness, buckling load, and natural frequency with the reduction of structural size, especially when the plate thickness is on the same order of the material length scale parameter.     

弹性地基上四边自由矩形厚板非线性大挠度弯曲

基于Ｒｅｉｓｓｎｅｒ－Ｍｉｎｄｌｉｎ一阶剪切变形板理论,采用摄动－Ｇａｌｅｒｋｉｎ混合法,给出双参数弹性地基上四边自由矩形中厚板在对称分布局部荷载作用下的大挠度弯曲渐近解,满足全部自由边界条件和控制方程,同时讨论弹性地基刚度系数对自由矩形厚板大挠度弯曲的影响。     

Compression analysis of rectangular elastic layers bonded between rigid plates

An elastic layer bonded between two rigid plates has higher compression stiffness than the elastic layer without bonding. While the finite element method can be applied to calculate the stiffness, the compression stiffness of bonded rectangular layers derived through a theoretical approach in this paper provides a convenient way for parametric study. Based on two kinematics assumptions, the governing equation for the mean pressure is derived from the equilibrium equations. Using the approximate shear boundary condition, the mean pressure is solved and the compression stiffness of the bonded rectangular layer is then established in an explicit single-series form. Through the solved pressure, the horizontal displacements are derived from the corresponding equilibrium equations, from which the shear stress on the bonding surface can be found. It is found that the effect of the rectangular aspect on the compression stiffness is significant only when Poisson’s ratio is near 0.5. For the smaller Poisson’s ratio, the compression stiffness of the rectangular layer can be approximated by the formula for the infinite-strip layer of the same shape factor.     

THE TRANSVERSE NONLINEAR FREE VIBRATION AND THE BUCKLING OF A COMPRESSIVE BAR ON AN ELASTIC FOUNDATION

基于Hamilton 原理,运用假设时间模态法,得到了弹性基础上压杆的横向非线性自由振动与屈曲的位移型常微分控制方程. 考虑一端固定另一端可移简支边界条件,采用打靶法得到了结构第一至第三阶结构频率与一阶屈曲载荷的数值结果. 结果表明：随轴心压力增加,结构频率减小;随弹性基础刚度增加,结构频率与屈曲载荷均增加;弹性基础刚度对结构频率的影响随振型阶数增加在减小;在小振幅的情形下,不同振型对一阶屈曲载荷的影响很小.     

Analytical solution of rectangular plate with in-plane variable stiffness

The bending problem of a thin rectangular plate with in-plane variable stiffness is studied. The basic equation is formulated for the two-opposite-edge simply supported rectangular plate under the distributed loads. The formulation is based on the assumption that the flexural rigidity of the plate varies in the plane following a power form, and Poisson’s ratio is constant. A fourth-order partial differential equation with variable coefficients is derived by assuming a Levy-type form for the transverse displacement. The governing equation can be transformed into a Whittaker equation, and an analytical solution is obtained for a thin rectangular plate subjected to the distributed loads. The validity of the present solution is shown by comparing the present results with those of the classical solution. The influence of in-plane variable stiffness on the deflection and bending moment is studied by numerical examples. The analytical solution presented here is useful in the design of rectangular plates with in-plane variable stiffness.     

LOCATING BEAM DAMAGE BASED ON INITIAL FLEXURAL STIFFNESS DETECTION

应变模态变化率指标在服役梁结构的损伤定位方面已有应用,但现有研究大多忽视了梁上的初始局部抗弯刚度具有离散性的情况,因此难以区分真实损伤和初始离散性造成的局部刚度变化,对于实际梁结构的应用效果不佳. 先提出了一种通过求解线性方程组来得到梁上各区间真实初始抗弯刚度的方法,然后采用应变模态变化率指标来进行损伤定位. 研究结果表明,该方法可以处理梁上初始局部抗弯刚度具有离散性的情况,实现准确的损伤定位.     

Vibration of quadrilateral embedded multilayered graphene sheets based on nonlocal continuum models using the Galerkin method

Free vibration analysis of quadrilateral multilayered graphene sheets(MLGS) embedded in polymer matrix is carried out employing nonlocal continuum mechanics.The principle of virtual work is employed to derive the equations of motion.The Galerkin method in conjunction with the natural coordinates of the nanoplate is used as a basis for the analysis.The dependence of small scale effect on thickness,elastic modulus,polymer matrix stiffness and interaction coefficient between two adjacent sheets is illustrated.The non-dimensional natural frequencies of skew,rhombic,trapezoidal and rectangular MLGS are obtained with various geometrical parameters and mode numbers taken into account,and for each case the effects of the small length scale are investigated.     

具有非轴对称刚度转轴的分岔

研究具有非轴对称刚度转轴的１／２亚谐共振和分岔,首先用Ｈａｍｉｌｔｏｎ原理导出运动微分方程,这是刚度系数周期性变化的参数激励方程,然后用多尺度法求得平均方程分岔响应方程和定常解,最后用奇异性理论分析分岔响应方程和定常解的稳定性,得到了局部分岔集和不同区域的不同分岔响应曲线。     

Cyclic density functional theory: A route to the first principles simulation of bending in nanostructures

We formulate and implement Cyclic Density Functional Theory (Cyclic DFT) — a self-consistent first principles simulation method for nanostructures with cyclic symmetries. Using arguments based on Group Representation Theory, we rigorously demonstrate that the Kohn-Sham eigenvalue problem for such systems can be reduced to a fundamental domain (or cyclic unit cell) augmented with cyclic-Bloch boundary conditions. Analogously, the equations of electrostatics appearing in Kohn-Sham theory can be reduced to the fundamental domain augmented with cyclic boundary conditions. By making use of this symmetry cell reduction, we show that the electronic ground-state energy and the Hellmann-Feynman forces on the atoms can be calculated using quantities defined over the fundamental domain. We develop a symmetry-adapted finite-difference discretization scheme to obtain a fully functional numerical realization of the proposed approach. We verify that our formulation and implementation of Cyclic DFT is both accurate and efficient through selected examples.The connection of cyclic symmetries with uniform bending deformations provides an elegant route to the ab-initio study of bending in nanostructures using Cyclic DFT. As a demonstration of this capability, we simulate the uniform bending of a silicene nanoribbon and obtain its energy-curvature relationship from first principles. A self-consistent ab-initio simulation of this nature is unprecedented and well outside the scope of any other systematic first principles method in existence. Our simulations reveal that the bending stiffness of the silicene nanoribbon is intermediate between that of graphene and molybdenum disulphide — a trend which can be ascribed to the variation in effective thickness of these materials. We describe several future avenues and applications of Cyclic DFT, including its extension to the study of non-uniform bending deformations and its possible use in the study of the nanoscale flexoelectric effect.     

A general tetrakaidecahedron model for open-celled foams

A micro-mechanics model for non-isotropic, open-celled foams is developed using an elongated tetrakaidecahedron (Kelvin model) as the repeating unit cell. Assuming the cell edges possess axial and bending rigidity, the mechanics of deformation of the elongated tetrakaidecahedron lead to a set of equations for the Young’s modulus, Poisson’s ratio and tensile strength of the foam in the principal material directions. These equations are written as a function of the cell edge lengths and cross-section properties, the inclination angle and the strength and stiffness of the solid material. This micro-mechanics model employs an elongated Kelvin model geometry which is more general than that employed by previous authors, as the size and shape of the repeating unit cell are defined by specifying three independent dimensions. As a result, the model accounts for an additional variation in the unit cell shape which is not accounted for in the previous models. The effect of this additional shape parameter on the non-isotropic stiffness and strength behavior is demonstrated and the advantages of this more general micro-mechanics model are illustrated.     

Atomistic simulation of free transverse vibration of graphene,hexagonal SiC,and BN nanosheets

Free transverse vibration of monolayer graphene, boron nitride (BN), and silicon carbide (SiC) sheets is inves-tigated by using molecular dynamics finite element method. Eigenfrequencies and eigenmodes of these three sheets in rectangular shape are studied with different aspect ratios with respect to various boundary conditions. It is found that aspect ratios and boundary conditions affect in a similar way on nat-ural frequencies of graphene, BN, and SiC sheets. Natural frequencies in all modes decrease with an increase of the sheet's size. Graphene exhibits the highest natural frequen-cies, and SiC sheet possesses the lowest ones. Missing atoms have minor effects on natural frequencies in this study.     

Karman型正交异性矩形薄板非线性弯曲的统一求解方法

本文对Karman型四边支承正交异性薄板在5种不同边界条件下的几何非线性弯曲进行了统一分析。所设的位移函数均为梁振动函数。它们精确地满足边界条件，利用Galerkin方法和位移函数的正交属性，转换控制方程为非线性代数方程。用“稳定化双共轭梯度法”求解稀疏矩阵线性方程组以及“可调节参数的修正迭代法”求解非线性代数方程组，最后给出了相应的数值结果。     

Post-buckling behaviour of slender structures with a bi-linear bending moment-curvature relationship

Certain classes of slender structures of complex cross-section or fabricated from specialised materials can exhibit a bi-linear bending moment-curvature relationship that has a strong influence on their global structural behaviour. This condition may be encountered, for instance, in (a) non-linear elastic or inelastic post-buckling problems if the cross-section stiffness may be well approximated by a bi-linear model; (b) multi-layered structures such as stranded cables, power transmission lines, umbilical cables and flexible pipes where the drop in the bending stiffness is associated with an internal friction mechanism. This paper presents a mathematical formulation and an analytical solution for such slender structures with a bi-linear bending moment versus curvature constitutive behaviour and subjected to axial terminal forces. A set of five first-order non-linear ordinary differential equations are derived from considering geometrical compatibility, equilibrium of forces and moments and constitutive equations, with hinged boundary conditions prescribed at both ends, resulting a complex two-point boundary value problem. The variables are non-dimensionalised and solutions are developed for monotonic and unloading conditions. The results are presented in non-dimensional graphs for a range of critical curvatures and reductions in bending stiffness, and it is shown how these parameters affect the structure's post-buckling behaviour.     

基于非局部弹性理论的纳米板横向振动

本文利用非局部弹性理论研究了单层石墨烯的纳米板的横向自由振动响应．通过迭代法推导了非局部应力表达,进一步通过哈密顿原理推导了纳米板的控制方程,应用纳维解法得到四边简支纳米板振动固有频率的数值解,并将本文研究结果与已有文献结果进行对比,进一步讨论了小尺寸效应,以及纳米板的三维尺寸和半波数对振动频率的影响．结果表明：非局部效应的存在使得纳米板的等效刚度和固有频率降低;半波数的增加则使得纳米板的固有频率提高．相关分析结果对基于二维纳米材料的新设备的设计和优化具有重要意义．     

类石墨烯二维原子晶体的微态理论模型

基于微态(Micromorphic) 连续介质理论,提出了针对类石墨烯二维原子晶体的新力学模型. 该模型以有限大小的布拉维单胞为基元体,考虑基元粒子的宏观位移和微观变形,依据微态理论基本方程,推导了全局坐标系下模型的主导方程. 然后针对布拉维单胞中含有两个原子的类石墨烯晶体,通过分析单胞中声子振动模式与基元体自由度的关系,获得了微态形式下声子色散关系的久期方程,并根据二维晶体声子色散特性对久期方程进行了简化,进而确定了类石墨烯晶体模型的本构方程. 最后,以石墨烯和单层六方氮化硼为例,利用简化的表达式拟合了它们面内声子色散关系数据,计算了模型材料的常数,石墨烯模型的等效杨氏模量、泊松比分别为1.05 TPa 和0.197,氮化硼分别为0.766 TPa 和0.225,均与已有的实验值相符合.