基于Hellinger-Reissner变分原理的应变梯度杂交元设计

从一般的偶应力理论出发，基于Hellinger-Reissner变分原理，通过对有限元 离散体系的位移试解引入非协调位移函数，得到了偶应力理论下有限元离散系统的能量相容 条件，并由此建立了应变梯度杂交元的应力函数优化条件. 根据该优化条件，构造了一 个C0类的平面4节点梯度杂交元，数值结果表明，该单元对可压缩和不可压缩状态的 梯度材料均可给出合理的数值结果，再现材料的尺度效应.     

全应变梯度挠曲电纳米梁有限单元法研究

挠曲电效应是一种存在于所有电介质材料中的特殊的力电耦合效应，本质上是应变梯度与电极化之间的线性耦合。然而，应变梯度会引入位移的高阶偏量，常给挠曲电问题的理论求解带来困难。且已有研究表明应变梯度弹性项会影响纳米结构中的力电耦合响应，但是现有的挠曲电研究大多忽略了应变梯度弹性的影响。因此，本文提出了一种既考虑应变梯度弹性，又考虑挠曲电效应的有效数值方法。基于全应变梯度弹性理论，建立了包含3个独立材料尺度参数的纳米欧拉梁的理论模型和有限元模型，提出了满足C₂弱连续的两节点六自由度单元。基于本文的有限单元法，以简支欧拉梁为例，通过分析讨论挠度、电势和能量效率，得到了挠曲电效应和应变梯度弹性项对梁的力电响应的影响。结果表明，挠曲电效应存在尺寸依赖性，且应变梯度弹性项在纳米电介质结构的挠曲电研究中的影响不可忽略。     

基于多尺度模型解释固体中的挠曲电效应

挠曲电效应是一种跨尺度的多场耦合现象。当前的宏观挠曲电理论均是基于应变梯度局部破坏晶体反演对称这一微观机理对该现象进行唯象描述。该宏观理论与基于晶格动力学及密度泛函理论的微观挠曲电理论模型之间存在较大差异。难以将两者结合用以跨尺度地研究材料中的挠曲电效应。针对该现状,本文基于前人提出的原子场理论,建立了一种新的多尺度挠曲电模型。并在该多尺度模型框架下解释了应变梯度诱发极化的微观机理。一方面,与基于连续介质力学的唯象理论不同,本文从材料微结构演化的角度推导了原子位移与极化的关系。另一方面,与通过晶格波假设原子位移的微观理论不同,本文得到的极化表达式更加真实和广义地解释了挠曲电效应。其能够适用于材料边界存在机械力作用,材料内部存在缺陷等复杂的情况。本文所建立的多尺度挠曲电模型能够为后续多尺度挠曲电效应的研究提供一些思路。     

微纳米尺度效应作用下充流微通道系统波动特性研究

结合非局部弹性应力/应变梯度耦合本构关系和流体非局部应力关系式,基于Euler梁理论,建立了充流微通道流固耦合波传导模型;根据耦合固体非局部应力/应变梯度弹性效应以及流体非局部效应,分别模拟了微通道和管腔内流体的尺度效应,推导得出了充流微通道在微纳米尺度的波动控制方程和边界条件。通过对控制方程的求解,分析了不同类型尺度效应对微通道的波动和振动特性的影响。结果显示,各类尺度效应对系统的动力学特性影响不同。微通道非局部弹性效应对波动产生阻尼,特别是对波长较短的波传导;而应变梯度弹性效应对波传导有促进作用,且该效应对波动的影响与波长无关;非局部效应和应变梯度效应对微通道刚度产生不同影响,非局部效应降低刚度,应变梯度效应增加刚度。     

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

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

非局部应变梯度理论下纳米梁的力学特性研究

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

基于精化锯齿理论的功能梯度夹心微板静弯曲模型

基于精化锯齿理论和新修正偶应力理论,建立了能够准确预测功能梯度夹心微板挠度、位移和应力的静弯曲模型。为了描述微板不同方向上的尺度效应,将两个正交材料尺度参数引入本文模型。以受双向正弦载荷作用的简支板为例,探究了夹心微板弯曲行为中尺度效应对结构刚度的影响。算例结果表明,当微板几何参数与材料尺度参数接近时,基于本文模型所测微板的最大弯曲挠度、局部位移和应力均小于传统精化锯齿理论给出的结果,捕捉到了尺度效应;尺度效应随着微板几何尺寸的增大而逐渐减弱,当微板几何尺寸远大于材料尺度参数时,尺度效应消失。此外,板的跨厚比和功能梯度变化指数也会对尺度效应产生一定影响。     

基于广义应变梯度理论的纳米梁挠曲电效应研究

挠曲电效应是应变梯度与电极化的耦合，它存在于所有的电介质材料中。在纳米电介质结构的挠曲电效应研究中，应变梯度弹性对挠曲电响应的影响一直以来被低估甚至被忽略了。根据广义应变梯度理论，应变梯度弹性中独立的尺度参数只有三个，而文献中所采用的一个或两个尺度参数的应变梯度理论只是它的简化形式。基于该理论，论文建立了考虑广义应变梯度弹性的三维电介质结构的理论模型，并以一维纳米梁为例研究了其弯曲问题的挠曲电响应及其能量俘获特性。结果表明，纳米梁的挠曲电响应存在尺寸效应，并且弹性应变梯度会影响结构挠曲电的尺寸效应，特别是当结构的特征尺寸低于尺度参数时。论文的工作为更进一步理解纳米尺度下的挠曲电机理和能量俘获特性提供理论基础和设计依据。     

纳米圆轴在弹性介质中的扭转振动分析

基于非局部应变梯度理论,考虑周围弹性介质的影响,研究纳米圆轴的扭转自由振动。首先通过Hamilton原理推导纳米圆轴扭转振动的控制方程及边界条件,接着采用微分求积法得到控制方程及三类边界条件的离散形式,最后由数值计算结果分析扭转振动特性。讨论了两个小尺度参数和弹性介质刚度的变化对振动频率的影响,并通过小尺度参数比对振动频率的影响分析两个尺度参数的耦合作用。研究结果表明,扭转自由振动频率随非局部参数增加而减小,随应变梯度尺度参数、弹性介质刚度增加而增大;当非局部参数大于应变梯度尺度参数时,小尺度效应体现为非局部效应,相反则体现为应变梯度效应。     

应变梯度理论进展

应变梯度理论是近10年来为解释材料在微米尺度下的尺寸效应现象而发展起来的一种新理论.首先综述了应变梯度理论近年的发展及其对材料力学行为研究方面的进展.其次主要介绍了不含高阶应力的一类应变梯度理论及其应用;最后对应变梯度理论的发展做了展望.      

Method of high-order lagrange multiplier and generalized variational principles of elasticity with more general forms of functionals

It is known[1]that the minimum principles of potential energy andcomplementary energy are the conditional variation principles underrespective conditions of constraints.By means of the method of La-grange multipliers,we are able to reduce the functionals of condi-tional variation principles to new functionals of non-conditionalvariation principles.This method can be described as follows:Mul-tiply undetermined Lagrange multipliers by various constraints,andadd these products to the original functionals.Considering these un-determined Lagrange multipliers and the original variables in thesenew functionals as independent variables of variation,we can see thatthe stationary conditions of these functionals give these undeter-mined Lagrange multipliers in terms of original variables.The sub-stitutions of these results for Lagrange multipliers into the abovefunctionals lead to the functionals of these non-conditional varia-tion principles.However,in certain cases,some of the undetermined Lagrangemultipliers ma     

A self-stabilized algorithm for enforcing constraints in multibody systems

A new algorithm is developed for the enforcement of constraints within the framework of nonlinear, flexible multibody system modeled with the finite element approach. The proposed algorithm exactly satisfies the constraints at the displacement and velocity levels, and furthermore, it achieves nonlinear unconditional stability by imposing the vanishing of the work done by the constraint forces when combined with specific discretizations of the inertial and elastic forces. Identical convergence rates are observed for the displacements, velocities, and Lagrange multipliers. The proposed algorithm is closely related to the stabilized index-2 or GGL method, although no additional multipliers are introduced in the present approach. These desirable characteristics are obtained without resorting to numerically dissipative algorithms. If high frequencies are present in the system, i.e. the system is physically stiff, dissipative schemes become necessary; the proposed algorithm is extended to deal with this situation.     

Divergence preserving reconstruction of the nodal components of a vector field from its normal components to edges

We have developed a new divergence preserving method for the reconstruction of the Cartesian components of a vector field from the orthogonal projection of a vector field to the normals to edges in two dimensional. In this method, discrete divergences computed from the nodal components and from the normal ones are exactly the same. Our new method consists of two stages. At the first stage, we use an extended version of the local procedure described in [J. Comput. Phys., 139 :406–409, 1998] to obtain a ‘reference’ nodal vector. This local procedure is exact for linear vector fields; however, the discrete divergence is not preserved. Then, we formulate a constrained optimization problem, in which this reference vector plays the role of a target, and the divergence constraints are enforced by using Lagrange multipliers. It leads to the solution of ‘elliptic’ like discrete equations for the cell‐centered Lagrange multipliers. The new global divergence preserving method is exact for linear vector fields. We describe all details of our new method and present numerical results, which confirm our theory. Copyright © 2016 John Wiley & Sons, Ltd.     

不规则区域平面弹性问题的正则区域重心插值配点法

将不规则区域嵌入到规则的矩形区域,在矩形区域上将弹性平面问题的控制方程采用重心Lagrange插值离散,得到控制方程矩阵形式的离散表达式。在边界节点上利用重心插值离散边界条件,规则区域采用置换法施加边界条件,不规则区域采用附加法施加边界条件,得到求解平面弹性问题的过约束线性代数方程组,采用最小二乘法进行求解,得到整个规则区域上的位移数值解。利用重心插值计算得到不规则区域内任意节点的位移值,计算精度可到10^-14以上。数值算例验证了所建立方法的有效性和计算精度。     

弹性理论中的临界变分及消除方法

临界变分现象是拉氏乘子法的固有特性，钱伟长应用高阶拉氏乘子消除了临界变分现象。本文将提出一种新的方法－凑合反推法，这种方法摒充了拉氏乘子法，把拉氏乘子所在的项目一个待定函数Ｆ代替。这样构成的泛函，作者称之为试泛函。而待定函数Ｆ的识别类似于拉氏乘子的识别。通过该法可以方便地构造出各种多变量广义变分原理，并且可以消除临界变分现象。     

基于对偶变量变分原理的完整约束系统保辛算法

基于对偶变量变分原理，选择积分区间两端位移为独立变量，构造了求解完整约束哈密顿动力系统的高阶保辛算法。首先，利用拉格朗日多项式对作用量中的位移、动量及拉格朗日乘子进行近似；然后，对作用量中不包含约束的积分项采用Gauss积分近似，对作用量中包含约束的积分项采用Lobatto积分近似，从而得到近似作用量；最后，在此近似作用量的基础上，利用对偶变量变分原理，将求解完整约束哈密顿动力系统问题转化为一组非线性方程组的求解。算法具有保辛性和高阶收敛性，能够在位移的插值点处高精度地满足完整约束。算法的收敛阶数及数值性质通过数值算例验证。     

Solution of the planar Newtonian stick–slip problem with the singular function boundary integral method

A singular function boundary integral method (SFBIM) is proposed for solving biharmonic problems with boundary singularities. The method is applied to the Newtonian stick–slip flow problem. The streamfunction is approximated by the leading terms of the local asymptotic solution expansion which are also used to weight the governing biharmonic equation in the Galerkin sense. By means of the divergence theorem the discretized equations are reduced to boundary integrals. The Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers, the values of which are calculated together with the singular coefficients. The method converges very fast with the number of singular functions and the number of Lagrange multipliers, and accurate estimates of the leading singular coefficients are obtained. Comparisons with the analytical solution and results obtained with other numerical methods are also made. Copyright © 2005 John Wiley & Sons, Ltd.     

Finite element analysis of shells of revolution using triangular discretization elements with corrective Lagrange multipliers

A thin shell algorithm based on a special triangular finite element is described. The column of nodal variable parameters of this element additionally contains corrective Lagrange multipliers, which allows one to improve the conditions of compatibility between the elements in the framework of the variational formulation of the problem.     

Involutory transformations and variational principles with multi-varoables in thin plate bending problems

In this paper, the generalizd variational principles of plate bending, froblems are established from their minimum potential energy principle and minimum complementary energy principle through the elimination of their constraints by means of the method of Lagrange multipliers. The involutory transformations are also introduced in order to reduce the order of differentiations for the variables in the variation. Funhermore, these involutory transformations become infacl the additional constraints in the varialion. and additional Lagrange multipliers may be used in order to remove these additional constraints. Thus, various multi-variable variational principles are obtained for the plate bending problems. However, it is observed that. nol all the constrainls ofva’iaticn can be removed simply by the ordinary method of linear Lagrange multipliers. In such cases, the method of high-order Lagrange multipliers are usedto remove iliose constrainls left over by ordinary linear multiplier method. And consequently. some funct ionals of more general forms are oblained for the generaleed variational principles of plate bending problems.