A variational approach for large deflection of ends supported nanorod under a uniformly distributed load,using intrinsic coordinate finite elements

This paper presents a variational formulation that can be used for large deflection analysis of ends supported nanorod including the coupled effects of nonlocal elasticity and surface stress under a uniformly distributed load. The variational formulation involving the strain energy due to bending of nonlocal elasticity including the surface stress effect and virtual work done by a uniformly distributed load, is expressed in terms of the intrinsic coordinates. The Lagrange multiplier technique is applied to impose the boundary conditions which accomplished in the formulation. The validity of the variational approach is ensured by Euler's equation, which identical to the one derived by the force equilibrium consideration of an infinitesimal nanorod segment. The finite element method and Newton–Raphson iterative procedure based on the variational formulation are used to solve a system of nonlinear equations. Moreover, the very large deflection configurations of ends supported nanorod are highlighted in this study.     

Exact analytical and numerical solutions of stability problems for a straight composite bar subjected to axial compression and torsion

Based on linearized equations of the theory of elastic stability of straight composite bars with a low shear rigidity, which are constructed using the consistent geometrically nonlinear equations of elasticity theory for small deformations and arbitrary displacements and a kinematic model of Timoshenko type, exact analytical solutions of nonclassical stability problems are obtained for a bar subjected to axial compression and torsion for various modes of end fixation. It is shown that the problem of direct determination of the critical parameter of the compressive load at a given torque parameter leads to transcendental characteristic equations that are solvable only if bar ends have cylindrical hinges. At the same time, we succeeded in obtaining solutions to these equations in terms of wave formation parameters of the bar; these parameters, in turn, enabled us to find the parameter of the critical load at any boundary conditions. Also, an algorithm for numerical solution of the problems stated is proposed, which is based on reducing the problems to systems of integroalgebraic equations with Volterra-type operators and on solving these equations by the method of mechanical quadratures (finite sums). It is demonstrated that such numerical solutions exist only for certain ranges of parameters of the bar and of the parameter of torque. In the general case, they can not be obtained by the numerical method used. It is also shown that the well-known solutions of the stability problem for a bar subjected to torsion or to compression with torsion are in correct. Translated from Mekhanika Kompozitnykh Materialov, Vol. 45, No. 2, pp. 167–200, March–April, 2009.     

解任意四边形板弯曲问题的样条有限元法

关于用样条函数解板的弯曲问题,[1]在1979年讨论了矩形板和菱形板的弯曲;[2]在1981年对简支边界条件的矩形板,用振动梁函数和B样条函数组合作为插值函数,得到了效率更高的算式;[3]在1984年对[2]作了补充,采用拉格朗日乘子法,得到了在各种边界条件下平板弯曲的近似解,但所讨论的仍然是矩形板.     

Extension of the scaled boundary finite element method to plate bending problems

The scaled boundary finite element method (SBFEM) is extended to the static analysis of thin plates in the framework of Kirchhoff's plate theory. The governing equations are transformed into scaled boundary coordinates. Applying a discrete form of the Kantorovich reduction method results in a set of ordinary differential equations, which can be solved in a closed-form analytical manner. The element stiffness matrices for bounded and unbounded media can be computed, using appropriate subsets of the analytical solution. Examples show the efficiency of the method, applied to plate bending problems. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

非均匀Reissner板弯曲的精确元法

本文在阶梯折算法和精确解析法的基础上,提出构造有限元的新方法——精确元法.该方法不用变分原理,可适用于任意变系数正定和非正定偏微分方程.利用该方法,得到Reissuer板弯曲的一个非协调单元,它具有十五个自由度.由于节点位移参数仅含有挠度和转角,因此处理任意边界条件非常容易.文中给出证明,位移和内力均收敛于精确解.由精确元法所得到的单元不仅能用于厚板,也可用于薄板.文末给出四个算例.算例表明,利用本文的方法,可获得满意的结果,并有较高的数值精度.     

基于非局部弹性应力场理论的纳米尺度效应研究:纳米梁的平衡条件、控制方程以及静态挠度

该文成功地解答了3个关于非局部应力理论用于纳米梁的问题:(ⅰ) 在绝大多数研究中,非局部效应增加导致纳米结构体刚度下降,其现象表现为弯曲挠度增加,固有频率减少,屈曲载荷下降,但为什么Eringen 的非局部弹性理论给出了完全相反的结论;(ⅱ) 为什么在某些研究结果中,非局部效应消失或是对研究结果无影响,比如纳米悬臂梁在集中载荷作用下的弯曲挠度; (ⅲ) 在高阶控制方程中,为什么高阶边界条件不存在．通过应用非局部弹性理论和精确变分原理分析纳米梁的弯曲问题,推导出全新的平衡条件、控制方程、边界条件和静态响应．这些方程和条件包含了与之前的相关研究结果符号相反的高阶微分项,这一差别导致了纳米效应对结构体的影响结果完全相反． 还证明之前为大家所公认的纳米梁静态或动态平衡条件实际上没有达到平衡,只有用等效弯矩代替非局部弯矩时,才可达到平衡．这些结论通常是可以被其它方法,比如应变梯度理论、耦合应力模型以及相关实验所证明．     

Analysis of notches and cracks in circular Kirchhoff plates using the scaled boundary finite element method

A new formulation of the scaled boundary finite element method (SBFEM) is presented for the analysis of circular plates in the framework of Kirchhoff's plate theory. Essential for the SBFEM is, that a domain is described by the mapping of its boundary with respect to a scaling centre. The governing partial differential equations are transformed into scaled boundary coordinates and are reduced to a set of ordinary differential equations, which can be solved in a closed-form analytical manner. If the scaling centre is selected at the root of an existent crack or notch, the SBFEM enables the effective and precise calculation of singularity orders of cracked and notched structures. The element stiffness matrices for bounded and unbounded media are derived. Numerical examples show the performance and efficiency of the method, applied to plate bending problems. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Torsional rigidity of reinforced concrete bars with arbitrary sectional shape

In this paper, the torsional rigidity of arbitrarily shape bar made of different materials is studied on the basis of theory of elasticity and finite element approach. With additional boundary conditions for the common boundaries of different materials from the continuous conditions of deformation and traction across the interior boundary, the torsion function can be solved numerically from the second boundary-value problem of potential theory. The traction jump boundary conditions across the interior surfaces are enforced in the alternate finite element approach. Several examples are shown to check the computational approach proposed, and the approach, at last, is applied to calculate the torsional rigidity of reinforced concrete bar and some multiply connected cross sections such as tower leg section of the Tsing Ma Bridge and other engineering structures.     

非均匀薄板弯曲的精确元法

本文在阶梯折算法的基础上,提出构造有限元的新方法——精确元法.它不用一般变分原理,可适用于任意变系数正定和非正定偏微分方程.利用该方法,得到薄板弯曲一个非协调三角形单元,它具有6个自由度.文中给出证明,位移和内力均收敛于精确解,并有很好的精度.文末给出算例.算例表明利用本文的方法,内力和位移均可获得满意的结果.     

Modelling groundwater changes due to fluctuating dam discharge

In this contribution, two numerical methods are used to predict the free surface changes in a sand bar due to fluctuations in river stage. One is a fixed-mesh, finite-element seepage formulation including Biot's consolidation theory, and the other is a boundary element method solution of the Laplace equation. Both models give overall predictions that are in good agreement with field data recorded at an instrumented sand bar in the Colorado River subjected to stage fluctuations from operation of the Glen Canyon Dam. The boundary element method appears to offer significant advantage in data preparation and computational times over the finite-element method for the problem studied in this paper.     

3D elasticity analytical solution for bending of FG micro/nanoplates resting on elastic foundation using modified couple stress theory

This paper addresses a 3D elasticity analytical solution for static deformation of a simply-supported rectangular micro/nanoplate made of both homogeneous and functionally graded (FG) material within the framework of modified couple stress theory. The plate is assumed to be resting on a Winkler–Pasternak elastic foundation, and its modulus of elasticity is assumed to vary exponentially along thickness. By expanding displacement components in double Fourier series along in-plane coordinates and imposing relevant boundary conditions, the boundary value problem (BVP) of plate system, including its governing partial differential equations (PDEs) of equilibrium are reduced to BVP consisting only ordinary ones (ODEs). Parametric studies are conducted among displacement and stress components developed in the plate and FG material gradient index, length scale parameter, and foundation stiffnesses. From the numerical results, it is concluded that the out-of-plane shear stresses are not necessarily zero at the top and bottom surfaces of plate. The results of this investigation may serve as a benchmark to verify further bending analyses of either homogeneous or FG micro/nanoplates on elastic foundation.     

具有弹性边拱及拉杆支承的双曲扁壳(一)

本文利用变分原理建立了具有弹性边拱及拉杆支承的双曲扁壳的平衡方程式及相应的边界条件和角点条件,这里假定边拱只在其本身平面内有刚度,边拱的扭转刚度和垂直于其平面的弯曲刚度都略去不计,本文研究了不许自由外伸的角点铰支条件,以及能够自由外伸的角点简支条件,前者相当于周边有拉杆限制角点外伸位移的情况,后者相当于周边无拉杆的情况.对于前者而言,本文近似地假定边拱沿弧方向的抗拉伸刚度为无穷大,亦即假定扁壳的边界切向位移为零,边拱只通过其垂直于扁壳平面的弯曲来产生弹性支承的作用.这些支承条件是近似地符合当前双曲扁壳屋盖的设计条件的.本文利用双三角级数解法求得具有弹性边拱及拉杆支承的方形底球面扁壳在自重载荷下的正确解.其特点在于先将边界条件积分处理使先满足角点条件,然后求解平面应力微分方程使满足积分后的边界条件.本文的结果直接给出拉杆中的拉力,对于具体设计问题是有用的.本文提出的积分形式的边界条件方法,对于弹性支承的边界问题在板壳方面的应用中是有它的普遍实用意义的.本文还给出了具有弹性边拱支承的方形底扁球壳的数值结果,角点为铰支或简支的,选取的参数值为λ=11.5936.计算结果表明级数收敛很快,并得出了边拱的弹性变形对壳体内力、内力矩及挠度分布规律的影响.     

An asymptotic approach to analysis of 3D beam stress-strain elastic states

A regular asymptotic approach to analysis of 3D beam stress–strain states is proposed based on the linear theory of elasticity and the method of integrodifferential relations. Using the integral formulation of Hooke's law and polynomial expansions of unknown stress and displacement functions with respect to transversal Cartesian coordinates the initial system of partial differential equations is reduced to a countable system of ordinary differential equations with constant coefficients. For rectilinear beams with rectangular cross–sections the consistent boundary value problems describing independently the compression and stretch, bends, and torsion states are derived. To find equilibrium stress and admissible displacement fields satisfying boundary conditions an effective numerical algorithm is worked out. Integral and local criteria for explicit bilateral estimates of resulted solution quality are proposed. The numerical results are presented and discussed. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

弹性薄板弯曲问题的等价的边界积分方程

用非解析开拓数学方法建立平面弹性薄板弯曲问题理论中具有间接变量的等价边界积分方程,并采用变分法进行了严格的说明.以往出现的三种间接变量边界积分方程,它们都不是等价的,对此我们进行了深入的讨论.     

Elasticity solutions for functionally graded rectangular plates with two opposite edges simply supported

England (2006) [13] proposed a novel method to study the bending of isotropic functionally graded plates subject to transverse biharmonic loads. His method is extended here to functionally graded plates with materials characterizing transverse isotropy. Using the complex variable method, the governing equations of three plate displacements appearing in the expansions of displacement field are formulated based on the three-dimensional theory of elasticity for a transverse load satisfying the biharmonic equation. The solution may be expressed in terms of four analytic functions of the complex variable, in which the unknown constants can be determined from the boundary conditions similar to that in the classical plate theory. The elasticity solutions of an FGM rectangular plate with opposite edges simply supported under 12 types of biharmonic polynomial loads are derived as appropriate sums of the general and particular solutions of the governing equations. A comparison of the present results for a uniform load with existing solutions is made and good agreement is observed. The influence of boundary conditions, material inhomogeneity, and thickness to length ratio on the plate deflection and stresses for the load x²yq are studied numerically.     

Variational formulation of a modified couple stress theory and its application to a simple shear problem

A variational formulation is provided for the modified couple stress theory of Yang et al. by using the principle of minimum total potential energy. This leads to the simultaneous determination of the equilibrium equations and the boundary conditions, thereby complementing the original work of Yang et al. where the boundary conditions were not derived. Also, the displacement form of the modified couple stress theory, which is desired for solving many problems, is obtained to supplement the existing stress-based formulation. All equations/expressions are presented in tensorial forms that are coordinate-invariant. As a direct application of the newly obtained displacement form of the theory, a simple shear problem is analytically solved. The solution contains a material length scale parameter and can capture the boundary layer effect, which differs from that based on classical elasticity. The numerical results reveal that the length scale parameter (related to material microstructures) can have a significant effect on material responses.      

Axisymmetrical vibrations of shells of revolution under abrupt loading

A frequency method is proposed for solving the problem of the vibrations of shells of revolution taking into account the energy dissipation under arbitrary force loading and on collision with a rigid obstacle. The Laplace transform is taken of the equation of the vibrations of a shell of revolution with non-zero initial conditions. For the inhomogeneous differential equation obtained, a variational method is used to solve the boundary-value problem, which consists of finding the Laplace-transformed boundary transverse and longitudinal forces and bending moments as functions of the boundary displacements. The equations of equilibrium of nodes, i.e. the corresponding equations of the finite-element method, are then compared, using results obtained earlier [1–4]. Amplitude-phase-frequency characteristics (APFCs) for the shell cross-sections selected are plotted. An inverse Laplace transformation is carried out using the clear relationship between the extreme points of the APFCs and the coefficients of the corresponding terms of the series in an expansion vibration modes [3]. In view of the fact that the proposed approach is approximate, numerical testing is used.     

Variational formulation of a modified couple stress theory and its application to a simple shear problem

A variational formulation is provided for the modified couple stress theory of Yang et al. by using the principle of minimum total potential energy. This leads to the simultaneous determination of the equilibrium equations and the boundary conditions, thereby complementing the original work of Yang et al. where the boundary conditions were not derived. Also, the displacement form of the modified couple stress theory, which is desired for solving many problems, is obtained to supplement the existing stress-based formulation. All equations/expressions are presented in tensorial forms that are coordinate-invariant. As a direct application of the newly obtained displacement form of the theory, a simple shear problem is analytically solved. The solution contains a material length scale parameter and can capture the boundary layer effect, which differs from that based on classical elasticity. The numerical results reveal that the length scale parameter (related to material microstructures) can have a significant effect on material responses.     

Imposing boundary conditions in Sinc method using highest derivative approximation

Sinc approximate methods are often used to solve complex boundary value problems such as problems on unbounded domains or problems with endpoint singularities. A recent implementation of the Sinc method [Li, C. and Wu, X., Numerical solution of differential equations using Sinc method based on the interpolation of the highest derivatives, Applied Mathematical Modeling 31 (1) 2007 1–9] in which Sinc basis functions are used to approximate the highest derivative in the governing equation of the boundary value problem is evaluated for structural mechanics applications in which interlaminar stresses are desired. We suggest an alternative approach for specifying the boundary conditions, and we compare the numerical results for analysis of a laminated composite Timoshenko beam, implementing both Li and Wu’s approach and our alternative approach for applying the boundary conditions. For the Timoshenko beam problem, we obtain accurate results using both approaches, including transverse shear stress by integration of the 3D equilibrium equations of elasticity. The beam results indicate our approach is less dependent on the selection of the Sinc mesh size than Li and Wu’s SIHD. We also apply SIHD to analyze a classical laminated composite plate. For the plate example, we experience difficulty in obtaining a complete system of equations using Li and Wu’s approach. For our approach, we suggest that additional necessary information may be obtained by applying the derivatives of the boundary conditions on each edge. Using this technique, we obtain accurate results for deflection and stresses, including interlaminar stresses by integration of the 3D equilibrium equations of elasticity. Our results for both the beam and the plate problems indicate that this approach is easily implemented, has a high level of accuracy, and good convergence properties.     

弹性理论中各种变分原理的分类

本文按弹性理论中各种变分原理的约束条件的不同,对所有变分原理进行分类.我们在前文中业已指出,应力应变关系这样的约束条件是不能用拉氏乘子法解除的.剩下的可能约束条件共有四种:(1)平衡方程,(2)应变位移关系,(3)边界外力已知的边界条件,和(4)边界位移已知的边界条件.弹性理论的各种变分原理中,有的只有一种约束条件,有的有两种或三种,最多只能有四种约束条件.这样一共可能有15种变分原理,但是每种变分原理既可以用应变能A表示,又可以用余能B表示.这样,我们一共应有30种形式完全不同的变分原理,我们全部列出了这三十种形式的变分原理.