力学学报  2018 , 50 (4): 853-862 https://doi.org/10.6052/0459-1879-18-050

Orginal Article

基于S-R和分解定理的三维几何非线性无网格法

宋彦琦*, 周涛

中国矿业大学(北京)力学与建筑工程学院, 北京 100083

THREE-DIMENSIONAL GEOMETRIC NONLINEARITY ELEMENT-FREE METHOD BASED ON S-R DECOMPOSITION THEOREM

Song Yanqi*, Zhou Tao

School of Mechanics and Civil Engineering , China University of Mining & Technology , Beijing 100083, China

中图分类号:  O331

文献标识码:  A

通讯作者:  *宋彦琦, 教授, 主要研究方向: 固体力学理论及工程应用等方面的研究. E-mail:songyq@cumtb.edu.cn*宋彦琦, 教授, 主要研究方向: 固体力学理论及工程应用等方面的研究. E-mail:songyq@cumtb.edu.cn

收稿日期: 2018-03-5

接受日期:  2018-06-19

网络出版日期:  2018-07-18

版权声明:  2018 力学学报期刊社 力学学报期刊社 所有

基金资助:  国家自然科学基金重点项目(41430640)和中国矿业大学(北京)深部岩土力学与地下工程国家重点实验室开放基金项目(SKLGDUEK1728)资助

展开

摘要

S-R(strain-rotation)和分解定理克服了经典有限变形理论的一些缺点, 使其可以为几何非线性数值分析提供可靠的理论基础. 对于大变形问题, 由于无网格法(element-free method)避免了对单元网格的依赖, 从而从根本上避免了有限单元法(finite element method, FEM)的单元畸变问题, 保证了求解精度. 因此, 将无网格法和S-R和分解定理结合起来势必能建立一套更加合理可靠的几何非线性数值计算方法. 目前基于S-R 定理的无网格数值方法研究较少并且只能用于二维平面问题的求解, 但实际上绝大多数问题都必须以三维模型来进行处理, 因此建立适用于三维情况的S-R无网格法是非常有必要的. 本文给出了适用于三维情况的S-R 无网格法: 采用由更新拖带坐标法和势能率原理推导出来的增量变分方程, 利用基于全局弱式的无网格Galerkin 法(EFG)得到了用于求解三维空间问题的离散格式. 利用MATLAB编制三维S-R 无网格法程序, 对受均布载荷的三维悬臂梁和四边简支矩形板结构的非线性弯曲问题进行了计算. 最后将所得的数值结果与已有文献进行了比较, 验证了本文的三维S-R无网格数值算法的合理性、有效性和准确性. 本文的三维S-R无网格数值算法可以作为一种可靠的三维几何非线性数值分析方法.

关键词: S-R和分解 ; 三维无网格法 ; 几何非线性

Abstract

Due to its overcoming the deficiencies of classic finite deformation theories, Strain-Rotation (S-R) decomposition theorem can provide a reliable theoretical support for the geometrically nonlinear simulation. In addition, due to it’s independent of the elements and meshes, the element-free method has more advantages to solve large deformation problems compared to finite element method (FEM), so that the accuracy is guaranteed as a result of avoiding the element distortions. Therefore, a more reasonable and reliable geometric nonlinearity numerical method certainly will be established by combining the S-R decomposition theorem and element-free method. But the studies of element-free methods based on S-R decomposition theorem in current literature are limited to two-dimensional problems. In most cases, three-dimensional mathematical-physical models must be established for the practical problems. Therefore it is very necessary to establish a three-dimensional element-free method based on the S-R decomposition theorem. Present study extends the previously work by authors into three-dimensional case: The incremental variation equation is derived from updated co-moving coordinate formulation and principle of potential energy rate in this paper, and three-dimensional discretization equations are obtained by element-free Galerkin method (EFG). By using the MATLAB programs based on the proposed 3D S-R element-free method in present study, the nonlinear bending problems for three-dimensional cantilever beam and simply supported plates subjected to uniform load are numerical discussed. The reasonability, availability and accuracy of 3D S-R element-free method proposed by present paper are verified through comparison studies, and the numerical method in present work can provide a reliable way to analysis 3D geometric nonlinearity problems.

Keywords: S-R decomposition theorem ; three-dimensional element-free method ; geometric nonlinearity

0

PDF (3084KB) 元数据 多维度评价 相关文章 收藏文章

本文引用格式 导出 EndNote Ris Bibtex

宋彦琦, 周涛. 基于S-R和分解定理的三维几何非线性无网格法[J]. 力学学报, 2018, 50(4): 853-862 https://doi.org/10.6052/0459-1879-18-050

Song Yanqi, Zhou Tao. THREE-DIMENSIONAL GEOMETRIC NONLINEARITY ELEMENT-FREE METHOD BASED ON S-R DECOMPOSITION THEOREM[J]. Acta Mechanica Sinica, 2018, 50(4): 853-862 https://doi.org/10.6052/0459-1879-18-050

引言

固体几何非线性和物理非线性问题的建模和求解方法等一直都是非线性力学的重点课题之一[1]. 由于非线性问题在数学求解上的复杂性, 很难找到问题的解析解, 因此更多时候是从数值模拟方面寻找解答. 在传统固体几何非线性问题的数值模拟方法中, 一般都采用全拉格朗日格式(total Lagrange formulation, 简称 T.L. 格式)和更新拉格朗日格式(updated Lagrange formulation, 简称 U.L. 格式)的有限元方法(FEM)来进行求解[2]. 但由于FEM对单元具有依赖性, 对于大变形问题, 单元的畸变会严重影响FEM的求解精度[3]. 为了避免对单元网格的依赖, 无网格法(element-free methods)是一个很好的选择, 例如, Li和Belytschko[4]利用无网格Galerkin 法(EFG)[5]对金属成形过程中的接触冲击大变形问题进行了分析. 邵玉龙等[6]在一致性无单元Galerkin法(consistent EFG, CEFG)的基础之上, 发展了自适应分析方法, 进一步提高它求解局部高梯度问题的计算效率. 杨建军和郑健龙[7]基于无网格局部强弱法(meshless local strong-weak, 简称MLSW), 对不规则域问题进行了求解.

然而在传统几何非线性无网格数值方法中, 绝大多数都是采用Green应变张量作为有限变形理论的基础. 例如Chen 等[8]、Li 等[9]、Liew等[10], 采用了无网格RKPM (reproducing kernel particle method)法对一些结构的大变形问题进行了数值计算. 又例如最近几年, Zhang 等[11]、Do等[12]分别利用无网格IMLS-Ritz法(improved moving least-square ritz method)和无网格RPIM 法(radial point interpolation method)对功能梯度碳纳米管增强复合材料板的大变形弯曲问题进行了数值讨论. 虽然Green应变张量克服了Cauchy应变张量在变形体作纯转动时产生虚假应变这一缺点, 但是由于没有导出相应的转动张量和应变张量相匹配, 因此会限制其在大位移和大转动问题中的应用[13,14]. 陈至达等[13,15-16] 提出的S-R(strain-rotation)和分解定理, 克服了Green应变张量的缺点, 应变和局部转动是同时发生的, 并且应变是唯一确定的[17]. 因此S-R和分解定理可以为几何非线性数值方法提供精确完备的数学力学基础.

陈至达等[18]基于S-R和分解定理, 采用拖带坐标描述法和势能率变分原理, 给出了用以数值求解的位移增量方程, 为S-R 和分解定理在数值方法中的应用奠定了基础. 到目前为止, 基于S-R和分解定理的数值模拟方法的研究大多数以FEM为主: 秦忠[19]提出了基于S-R定理的节理单元和曲壳单元, 编制了非线性有限元计算程序SNAP, 对岩体蠕变和壳体大变形问题进行了求解; 尚勇等[20,21,22,23] 进一步拓展了S-R定理和拖带坐标描述法在有限元中的应用, 建立了连续体大变形变边界变分公式和大变形接触问题的计算模型; 李平等[24,25]提出了更新拖带坐标法(updated co-moving system, 简称U.C. 法)的大变形有限元计算方法, 编制了程序NASRAC, 并与采用U.L. 格式的程序ADINA 进行了比较; 高立堂等[26,27]采用李平的更新拖带坐标有限元法对火灾作用下的钢筋混凝土板进行了数值计算; 何满潮等[28]开发了有限元程序“深部软岩工程大变形力学设计分析系统”, 对巷道的不同埋深下进行开挖的巷道围岩大变形进行数值模拟. 除了FEM之外, 一些学者将S-R定理也应用到了其他一些数值方法中. 谢和平[29,30]导出了基于S-R定理的非线性边界元计算格式, 编写了NBEM程序; 高亚楠等[31]采用S-R 和分解理论对原有非连续变形分析(DDA)程序进行了修正; 罗丹[32]、 陈芳祖和罗丹[33]、宋彦琦等[34] 都采用了更新拖带坐标描述法, 导出了基于S-R 定理的二维无网格Galerkin法离散方程, 对二维平面大变形梁结构算例进行了求解; 最近, Fan等[2]同样基于S-R 定理和更新拖带坐标描述法, 利用动量守恒定律给出了动态增量变分方程, 并结合数值流形法(numerical manifold method, NMM), 建立了二维情况下的S-R 动态数值流形法(S-R-D-based NMM); 随后Fan等[35,36]又根据自己推导的动态增量变分方程, 给出了三维情况下的离散格式[36], 并结合DDA, 分别建立了基于S-R 定理的二维DDA方法(SR-2D-DDA)[35] 和三维[36]DDA方法(SR-3D-DDA).

从前面的文献综述可以看出, 基于S-R定理的无网格法的研究依然较少, 由于在大变形问题分析中, 无网格法相对于FEM更有优势, 因此目前结合S-R定理和无网格数值方法的研究依然值得关注. 另外文献[32,33,34]只是对二维问题进行了讨论, 并未给出基于S-R和分解定理的三维无网格法离散格式. 基于以上认识, 本文将S-R无网格法进一步延伸到三维情况: 采用由更新拖带坐标法和势能率原理推导出来的增量变分方程, 利用基于全局弱势的无网格Galerkin法得到了用于求解三维空间问题的离散格式, 其本质边界条件的施加是通过直接修改系统离散方程中的刚度矩阵和总体力向量来实现的. 利用MATLAB编制三维S-R 无网格法程序, 对受均布载荷的三维悬臂梁和四边简支矩形板结构的非线性弯曲问题进行了计算. 最后将所得的数值模拟结果与现有文献进行比较, 验证了本文的三维S-R无网格数值模拟算法的合理性和准确性.

1 S-R和分解定理

S-R和分解定理[13]证明: “给定一个物理可能的位移函数, 此函数在形变体点集域内是单值连续的, 处处具有一阶导数, 则此运动变换总可以分解为正交与对称两个子变换的直和. 正交变化体现点集的转动, 而对称变换体现点集的形变.”

为了描述欧式空间 E3中任意变形体的运动, 可以选取两个参考系, 如图1所示. (1) 全局固定参考系 Xi(i=1,2,3); (2) 嵌含在变形体中的自然拖带坐标系 xi(i=1,2,3). 其坐标线会随着物体形变而变形.

图1   三维欧式空间中变形体示意图.

Fig.1   A deformable body in three-dimension Euclidean space

图1中, rR分别为变形前和变形后的位置矢量; u为位移矢量; g0igi(i=1,2,3)分别为局部拖带坐标系在一点变形前后的基矢量. 位置矢量(r, R)和位移矢量(u)满足如下关系

R=r+u(1)

变形前后拖带坐标系中的基矢量分别定义如下

g0i=rxi,i=1,2,3(2)

gi=Rxi,i=1,2,3(3)

将式(2)和式(3)代入式(1)得到

Rxi=rxi+uxi(4)

在曲线坐标系下

uxi=xiujg0j=uj|ig0j(5)

其中

uj|i=ujxi+Γikj0uk(6)

Γikj0=12gjl0g0lixk+g0lkxi-g0ikxl(7)

其中, uj|i称为一阶逆变分量 ujxi的协变导数或者绝对导数; Γikj0称为第二类Christoffel符号. 进而可以得到变形前后局部拖带坐标系基矢量的转换关系

gi=Fijg0j(8)

Fij=δij+uj|i(9)

其中, F称为局部坐标系基矢的变换函数(或称为位移梯度函数); δij为Kronecker符号.

换句话说, S-R和分解定理可以表述为: 任何一个可逆的线性微分变换函数 F都存在一个唯一的直和分解

F=S+R(10)

其中, S为对称的子变换, 表示应变张量; R为正交的子变换, 表示转动张量. 文献[13]证明了和分解定理的存在性和唯一性, 这里只是给出其中相关的数学结果.

应变张量和转动张量分别为

Sji=12ui|j+ui|jT-LkiLjk(1-cosθ)(11)

Rji=δji+Ljisinθ+LkiLjk(1-cosθ)(12)

其中, Lji为转轴单位矢量

Lji=12sinθui|j-ui|jT(13)

其中, θ为平均整旋角

2 基于更新拖带坐标法的增量变分方程

2.1 增量变分方程

假设从0到 t时刻的所有时间点的物理量已经得到, 随后需要求解的是 t+Δt时刻的物理量, 这一求解过程是具有代表性的, 这样可以重复应用该过程直到整个求解路径完成. 由势能率原理, t+Δt时刻的变分为[25]

δt+ΔtJ=Ω(t+Δt)t+Δtσjiδt+ΔtVj||idΩ-t+ΔtR=0(15)

式中,δ为变分符号;J为泛函;σij为应力张量;表示在 t+Δt时刻速度 t+ΔtVi关于 t+Δt时刻拖带系协变基矢 t+Δtgi的协变导数; t+ΔtR为外力虚功

Ȃt+ΔtR=Sp(t+Δt)t+ΔtPiδt+ΔtVidS+Ω(t+Δt)t+Δtρfiδt+ΔtVidΩ(16)

其中, Pf分别为面力和体力.

假设所考虑的问题为准静力学问题, 根据速度梯度和应力的线性近似逼近有

t+ΔtVk||i=δlk-Δuk||jVj||i(17)

t+Δtσji=tσji+Δtσ̇ji(18)

其中, Δuk||j是位移增量 Δuktgi系下的协变导数; σ̇ji是以 t时刻拖带系 tgi作为基准的拖带应力速率. 将式(17)和式(18)代入式(15)中并利用关系 Δuj||k=ΔtVj||k得到

Δtσ̇jiδVj||idΩ-ΔtσjiVj||kδVk||idΩ=t+ΔtR-σjiδVj||idΩ(19)

2.2 更新拖带坐标法

在更新拖带坐标法[25]中, 是以 t时刻初始拖带系 tg0i内的速度分量 tV0i作为基本插值量. 如图2所示, 设质点在初始时刻 t0时位于空间位置 P点, P处的初始拖带基矢为 g0i.t时刻, 质点移至空间位置 P'处, P'处的初始拖带系基矢为 tg0i, 瞬时拖带系基矢为 tgi. 在由 t时刻至 t+Δt时刻的增量步上, 选用 t时刻的初始拖带坐标系 tg0i(一般取为标准直角坐标系)作为参考基准. 这样, t时刻瞬时拖带系 tgi上的 tVi||j表示为

tVi||j=tVk0|lxitxktxlxj(20)

tVk0|l同样有广义欧拉运动学公式

tVk0|l=S̅̇lk+R̅̇lk(21)

式中

S̅̇lk=12iVk0|l+iVk0|lT(22)

R̅̇lk=Ltkθ̅̇=12iVk0|l+iVk0|lT(23)

图2   更新拖带坐标系.

Fig.2   Update of the three-dimension co-moving coordinate

其中, 分别表示应变速率和转动速率; 顶标“‒”表示相应物理量是在 t时刻初始拖带系 tg0i下的物理量.

将式(20)运用到式(19)中

Δtσ̅̇jiδtVj0|idΩ-Δtσ̅jitVj0|kδtVk0|idΩ=t+ΔtR-σ̅jiδtVj0|idΩ(24)

上式中

σ̅ji=σlkxltxjtxixk(25)

σ̅̇ji=σ̇lkxltxjtxixk(26)

tg0i系为基准的速率型物性方程

σ̅̇ji=C̅jlikS̅̇kl(27)

其中

C̅jlik=Cntmsxntxjtxixmxttxltxkxs(28)

对于各向同性体有 C̅jlik=Cjlik.

将式(27)代入式(24)中

ΔtC̅jlikS̅̇klδS̅̇ijdΩ-Δtσ̅jiS̅̇kjδS̅̇ikdΩ-Δtσ̅jiR̅̇kjδR̅̇ikdΩ-Δtσ̅jiS̅̇kjδR̅̇ik+σ̅jiR̅̇kjδS̅̇ikdΩ+Δtσ̅̇jiδR̅̇ijdΩ=t+ΔtR-σ̅jiδS̅̇ijdΩ(29)

如果不考虑体力矩的影响, 假设初始系为正交系, 文献[25]利用如下结果

σ̅̇ji=σ̅kiS̅̇jk-S̅̇kiσ̅jk(30)

对式(29)进行简化

ΔtC̅jlikS̅̇klδS̅̇ijdΩ-Δtσ̅jiS̅̇kjδS̅̇ikdΩ-Δtσ̅jiR̅̇kjδR̅̇ikdΩ-Δtσ̅jiR̅̇kjδS̅̇ik+σ̅kjS̅̇jiδR̅̇ikdΩ=t+ΔtR-σ̅jiδS̅̇ijdΩ(31)

此式可作为全局弱形式无网格法的增量变分方程.

更新拖带坐标法最大的特点是, 在每一次迭代步结束后, 重新选择初始拖带系, 不仅保证了应力速率的客观性, 自动消除了转动的影响, 使应力分量不随时间改变, 而且大大简化了增量变分方程求解中的系数矩阵, 降低了增量变分方程的求解难度,利于算法实现和推广应用[34].

3 基于无网格Galerkin法的系统离散方程

在无网格方法中, 问题域被离散成一组节点 xk(k=1,2,,n)来进行求解. 根据无网格Galerkin 法, 对于 t时刻节点位移增量 tΔu和速度 tV, 采用移动最小二乘(MLS) 插值形函数进行近似

tΔu=k=1ntϕk(x)tΔuk(32)

tV=k=1ntϕk(x)tV(33)

其中, tϕkt时刻节点 k的MLS的形函数; tΔuktV分别是 t时刻节点 k的位移增量参数和速度参数.

类似于式(6), 对于 Vj||i同样有

Vj||i=Vjxi+ΓiljVl(34)

将式(32)和式(33)分别代入式(6)和式(34)中

tΔuj|i=k=1ntϕkxitΔukj+Γiljk=1ntϕktΔukl(35)

tVj|i=k=1ntϕkxitVkj+Γiljk=1ntϕktVkl(36)

另外

ΔS̅=ΔtS̅̇,ΔR̅=ΔtR̅̇(37)

从而可以将式(22), (23)和式(37)改写成矩阵形式

ΔS̅=BstΔu,ΔR̅=BrtΔu(38)

S̅̇=BstV,R̅̇=BrtV(39)

以上两式中

tΔu=tΔu11,tΔu12,tΔu13,,tΔun1,tΔun2,tΔun3(3n×1)T(40)

tV=tV11,tV12,tV13,,tVn1,tVn2,tVn3(3n×1)T(41)

Bs=B1s,B2s,,Bns(6×3n)(42)

Br=B1r,B2r,,Bnr(6×3n)(43)

ΔS̅=ΔS̅11,ΔS̅22,ΔS̅33,2ΔS̅21,2ΔS̅32,2ΔS̅316×1T(44)

ΔR̅=2ΔR̅21,2ΔR̅32,2ΔR̅313×1T(45)

S̅̇=S̅̇11,S̅̇22,S̅̇33,2S̅̇21,2S̅̇32,2S̅̇316×1T(46)

R̅̇=2R̅̇21,2R̅̇32,2R̅̇313×1T(47)

Bis=Nkx1+Γ111NkΓ121NkΓ131NkΓ212NkNkx2+Γ222NkΓ232NkΓ133NkΓ323NkNkx3+Γ333NkNkx2+(Γ211+Γ112)NkNkx1+(Γ221+Γ122)Nk(Γ231+Γ132)Nk(Γ312+Γ213)NkNkx3+(Γ322+Γ223)NkNkx2+(Γ332+Γ233)NkNkx3+(Γ311+Γ113)Nk(Γ321+Γ123)NkNkx1+(Γ331+Γ133)Nk(48)

Bir=Nkx2+(Γ211+Γ112)Nk-Nkx1+(Γ221-Γ122)Nk(Γ231-Γ132)Nk(Γ312-Γ213)NkNkx3+(Γ322+Γ223)Nk-Nkx2+(Γ332-Γ233))NkNkx3+(Γ311-Γ113)Nk(Γ321+Γ123)Nk-Nkx1+(Γ331+Γ133)Nk(49)

在更新拖带坐标法中, 如果初始拖带系和直线直角系胚

Γjki=Γjki0=0(50)

将式(38)和式(39)代入式(31)中, 并根据 δtVT的任意性, 得到增量变分方程的矩阵形式

tKL-tKNtΔu=t+ΔtR-tF(51)

其中

tKL=Ω(t)BsTDBsdΩ(52)

tKN=Ω(t)BsTσ1BsdΩ-Ω(t)BrTσ2BrdΩ-Ω(t)BsTσ3BrdΩ-Ω(t)BrTσ3TBsdΩ(53)

t+ΔtR=sp(t)ϕTPdS+Ω(t)ϕTρfdS(54)

tF=Ω(t)BsTσdΩ(55)

上式中, D是物性系数矩阵

σ={σ11,σ22,σ33,σ21,σ32,σ31}T(56)

σ1=σ1100σ2120σ3120σ220σ122-1mmσ322000σ330σ232σ132σ122σ2120σ11+σ224σ314σ3240σ232σ322σ134σ22+σ334σ124σ1320σ312σ234σ214σ11+σ334(57)

σ2=-14σ11+σ22-σ31σ32-σ13σ22+σ33σ12σ23σ21σ11+σ33(58)

σ3=-12σ210-12σ3112σ12-12σ320012σ2312σ1314(σ11-σ22)-14σ31-14σ32-14σ1314(σ22-σ33)14σ12-14σ2314σ2114(σ11-σ33)(59)

在本文中, 采用直接法来施加本质边界条件, 即通过直接修改系统增量方程(51)中的刚度矩阵和总体力向量. 在本文的无网格法数值计算过程中, 其它一些必要的参数设置见表1.

表1   本文三维SR-EFG的必要参数设定

Table 1   The essential parameters for 3D-SR-EFG

NameType
nodes distributionsregular nodes
influence domainsphere: scaling factor d=2.5
basis functionslinear basis function:
pT(x)={1,x1,x2,x3},m=4
ui(xi)=a1(xi)+a2(xi)xi1+
a3(xi)xi2+a4(xi)xi3
weight functionCSRBF2:
W=0,r>1(1-r)6(6+36r+82r2+72r3+30r4+5r5),r#x2264;1
numerical integrationGauss integration: 3×3×3

新窗口打开

4 数值算例

4.1 悬臂梁结构的非线性弯曲计算

图3所示为一上表面受均布载荷 q的三维悬臂梁结构. 已知梁的长宽高分别为 L=254mm, D=H=25.4mm, 杨氏模量 E=82.74MPa, 泊松比 . 图4给出了该算例中梁右端中点的无量纲挠度值与无量纲均布载荷值之间的关系, 其中包括本文数值方法收敛性的验证, 以及本文结果与罗丹[32] 和Holden[37]的结果、小变形线性理论的结果之间的比较. 由不同节点数之间的对比可以反映出本文三维S-R无网格法在该算例中的收敛性特点, 可以发现当采用 21×3×3个节点时, 计算结果已经趋于稳定. 通过与罗丹[32]的结果进行比较, 可以发现本文的结果与罗丹[32] 的结果之间存在微小的差异. 在罗丹的所有算例中, 结构都是有厚度的, 因此实际上都应该是三维空间结构. 当他们使用自己的二维S-R无网格法来求解这些三维问题时, 必定需要满足如下假设(以这里的悬臂梁为例): 所有物理量与厚度方向的坐标( X2)是没有关系的, 即所有物理量只是另外两个方向坐标( X1,X3)的函数; 沿厚度方向的物理分量( σi3,εi3,u3,i=1,2,3)都为零. 本文的三维S-R 无网格数值法未对该算例中的问题作以上的假设近似简化处理, 因此本文方法对于该算例的求解应该具有更高的精度. 而由于Holden[37] 采用不同的有限变形理论, 因此与之比较起来差异性会更加明显一点. 通过该算例的求解, 可以验证本文三维S-R 无网格法的收敛性和合理性.

图3   受均布载荷的空间悬臂梁.

Fig.3   3D cantilever beam subjected to a uniform load

图4   梁右端中点挠度与载荷的关系

Fig.4   The right end central deflection of beam versus load

4.2 矩形板结构的非线性弯曲计算

图5所示为一上表面受均布载荷q的三维矩形板结构示意图, 其长度为b, 宽度为a, 厚度为h. 在该算例中, 假定板为各向同性材料(µ=0.316), 并且a=b, 其四边受简支约束. 如图6, 给出了平板中面中点的无量纲挠度与无量纲均布载荷强度之间的关系曲线, 与之作为对比的结果是来自于Shen[38]、Levy[39]和Gorgi[40]. 首先可以看到在该算例中采用节点数11×11×3时结果就已经收敛了. 另外从曲线图中可以看出, 对于该算例本文计算得到的结果小于Shen[38]、Levy[39]、Gorgi[40]给出的结果. Shen[38] 、Levy[39]、Gorgi[40]都是在位移几何关系中加入冯·卡门(Von Kármán)应变项来实现非线性的考虑, 只是三者采用了不同板壳理论和不同的偏微分方程求解方法. Levy[39]采用的是经典板壳理论(classic plate theory, CPT), 通过三角级数展开法对方程进行了求解. Gorgi[40]则利用一阶剪切变形理论(first order shear deformation theory, FSDT), 通过有限差分法和有限单元法进行的求解. Shen[38]是基于Reddy的高阶剪切变形理论(higher order shear deformation plate theory, HSDT), 最后通过摄动法对非线性方程进行求解. 但是经典板壳理论和绝大多数剪切变形板壳理论包括一阶剪切变形理论和Reddy的高阶剪切变形理论, 会采用一些假设, 如假设沿厚度方向的正应变εz=0. 而本文方法所采用的理论基础未做任何近似假设, 因此同本文的三维方法相比较而言, Shen[38] 、Levy[39] 和Gorgi[40]三者的准确性应该相对较低. 在传统板壳几何非线性分析中, 都是采用Von Kármán应变项来实现非线性的考虑[41]. 而作者认为Von Kármán应变在本质上为Green 应变张量, 而Green 应变张量理论最致命的缺点是没有导出相应的转动张量和应变张量相匹配[13,14]. 总结以上的讨论, 该算例同样验证了本文三维S-R 无网格法的收敛性, 并且同Shen[38]、Levy[39] 和Gorgi[40] 之间的差异性也是可以解释的.

图5   受均布载荷的矩形板结构示意图.

Fig.5   3D rectangle plate subjected to a uniform load

图6   平板中面中点挠度与载荷关系(各向同性板).

Fig.6   The central deflection of rectangle plate’s mid-plane versus load (isotropic plate)

同样以图5为例, 若板为正交各向异性材料, 杨氏模量有 E1=2.069×104MPa, E2=8.826×103MPa, 剪切模量有 G12=G13=G23=2.551×103MPa, 泊松比 µ=0.32, 其长宽高满足 a=b=304.8mm, h=3.5mm. 图7给出该算例的平板中面中点的无量纲挠度与无量纲均布载荷强度之间的关系. 同样可以发现当节点数为 11×11×3时, 计算结果就已经收敛. 这里还可以发现, 对于正交各向异性板, 本文方法的计算结果是大于Shen[38] 给出的结果, 这是与前文的各向同性板算例不一样的. 相对于Zaghoul 等[42] 给出的经典解(CPT) 和Shen[38] 的结果, 本文的三维S-R 无网格法计算结果与Zaghoul[42] 给出的实验值最接近. 因此充分说明了本文三维S-R无网格法的准确性, 也说明了S-R 和分解理论相对于其它有限变形理论的优越性.

图7   平板中点挠度与载荷的关系(正交各向异性板).

Fig.7   The central deflection of rectangle plate’s mid-plane versus load (orthotropic plate)

5 结 论

本文将基于S-R和分解定理的二维无网格数值方法扩展到三维情况, 给出了3D-SR-EFG数值计算方法, 并通过对梁、板结构的非线性弯曲问题典型算例进行了对比分析, 验证了该方法在求解三维几何非线性问题时的合理性、有效性和准确性. 另外, 对于本文提出的方法在与其他方法进行对比时所表现出来的一些差异性, 也给出了合理的解释. 在受均布载荷的四边简支正交各向异性方形板的算例中, 本文提出的数值计算方法是在所有对比方法中最接近实验结果的一种方法, 说明了3D-SR-EFG方法可以作为一种更加可靠的几何非线性数值计算方法.

The authors have declared that no competing interests exist.


参考文献

[1] 谢和平. 非线性力学理论与实践. 徐州: 中国矿业大学出版社, 1997

[本文引用: 1]     

(Xie Heping.Theory of Nonlinear Mechanics and Its Application. Xuzhou: China University of Mining and Technology Press, 1997 (in Chinese))

[本文引用: 1]     

[2] Fan H, Zheng H, He S.

S-R decomposition based numerical manifold method

.Computer Methods in Applied Mechanics and Engineering, 2016, 304: 452-478

[本文引用: 2]     

[3] Liu GR, Gu YT.An Introduction to Meshfree Methods and Their Programming. Netherlands: Springer, 2005

[本文引用: 1]     

[4] Li G, Belytschko T.

Element-free Galerkin method for contact problems in metal forming analysis

.Engineering Computations, 2001, 18: 62-78

[本文引用: 1]     

[5] Belytschko T, Lu YY, Gu L.

Element-free Galerkin methods

.International Journal for Numerical Methods in Engineering, 1994, 37(2): 229-256

[本文引用: 1]     

[6] 邵玉龙, 段庆林, 高欣.

自适应一致性高阶无单元伽辽金法

. 力学学报, 2017, 49(1): 105-116

[本文引用: 1]     

(Shao Yulong, Duan Qinglin, Gao Xin, et al.

Adaptive consistent high order element-free Galerkin method

.Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(1): 105-116 (in Chinese))

[本文引用: 1]     

[7] 杨建军, 郑健龙.

无网格局部强弱法求解不规则域问题

. 力学学报, 2017, 49(3): 659-666

[本文引用: 1]     

(Yang Jianjun, Zheng Jianlong.

Meshless local strong-weak (MLSW) method for irregular domain problems

.Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(3): 659-666 (in Chinese))

[本文引用: 1]     

[8] Chen JS, Pan C, Wu CT, et al.

Reproducing Kernel particle methods for large deformation analysis of non-linear structures

.Computer Methods in Applied Mechanics & Engineering, 1996, 139(1-4): 195-227

[本文引用: 1]     

[9] Li S, Hao W, Liu WK.

Numerical simulations of large deformation of thin shell structures using meshfree methods

.Computational Mechanics, 2000, 25(2-3): 102-116

[本文引用: 1]     

[10] Liew KM, Ng TY, Wu YC.

Meshfree method for large deformation analysis-a reproducing kernel particle approach

.Engineering Structures, 2002, 24(5): 543-551

[本文引用: 1]     

[11] Zhang LW, Song ZG, Liew KM.

Nonlinear bending analysis of FG-CNT reinforced composite thick plates resting on Pasternak foundations using the element-free IMLS-Ritz method

.Composite Structures, 2015, 128: 165-175

[本文引用: 1]     

[12] Do VNV, Lee CH.

Bending analyses of FG-CNTRC plates using the modified mesh-free radial point interpolation method based on the higher-order shear deformation theory

.Composite Structures, 2017, 168: 485-497

[本文引用: 1]     

[13] 陈至达. 有理力学:非线性连续体力学. 徐州: 中国矿业大学出版社, 1988

[本文引用: 5]     

(Chen Zhida. Rational Mechanics.Xuzhou: China University of Mining and Technology Press, 1988 (in Chinese))

[本文引用: 5]     

[14] 陈至达. 杆、板、壳大变形理论. 北京: 科学出版社, 1994

[本文引用: 2]     

(Chen Zhida.Finite Deformation Theory of Link, Plate and Shell. Beijing: Science Press, 1994 (in Chinese))

[本文引用: 2]     

[15] 陈至达.

连续介质有限变形力学几何场论

. 力学学报, 1979, 15(2): 107-117

[本文引用: 1]     

(Chen Zhida.

Geometric field theory of finite deformation mechanics for continuum

.Chinese Journal of Theoretical and Applied Mechanics, 1979, 15(2): 107-117 (in Chinese))

[本文引用: 1]     

[16] 陈至达.

连续介质有限变形力学几何场论(续)(弹性有限变形能量原理)

. 清华大学学报(自然科学版), 1979, 19: 45-57

[本文引用: 1]     

(Chen Zhida.

Geometric field theory of finite deformation mechanics for continuum (continued)(elastic finite deformation energy principle)

.Journal of Tsinghua University (Science and Technology), 1979, 19: 45-57 (in Chinese))

[本文引用: 1]     

[17] Chen M, Liang J, Chen X, et al.

On uniqueness, existence and objectivity of S-R decomposition theorem

.Applied Mathematics and Mechanics (English Edition), 1997, 18(9): 817-823

[本文引用: 1]     

[18] 陈至达.

连续体力学中有限变形与转动的计算增量法

. 应用数学和力学, 1985, 6(9): 769-774

[本文引用: 1]     

(Chen Zhida.

Increment method of Finite deformation and rotational computing in continuum mechanics

.Applied Mathematics and Mechanics, 1997, 18(9): 817-823 (in Chinese))

[本文引用: 1]     

[19] 秦忠. 基于新的大变形理论的非线性有限元及其应用. [博士论文]. 北京: 中国矿业大学(北京), 1986

[本文引用: 1]     

(Qin Zhong. Nonlinear finite element and its applications based on new large deformation thoery. [PhD Thesis]. Beijing: China University of Mining and Technology (Beijing), 1986 (in Chinese))

[本文引用: 1]     

[20] 尚勇. 带有接触摩擦边界弹塑性大变形问题的理论和有限元分析. [博士论文]. 北京: 中国矿业大学(北京), 1987

[本文引用: 1]     

(Shang Yong. The theory and finite element analysis of the elastoplasticity large deformation problem with contact friction boundary. [PhD Thesis]. Beijing: China University of Mining and Technology (Beijing), 1987 (in Chinese))

[本文引用: 1]     

[21] 尚勇, 陈至达.

橡胶环直径方向受压大变形问题

. 计算结构力学及其应用, 1989, 6(4): 13-20

[本文引用: 1]     

(Shang Yong, Chen Zhida.

Large deformation problem of rubber ring compressed in diameter direction

.Journal of Computational Structural Mechanics and Applications, 1989, 6(4): 13-20 (in Chinese))

[本文引用: 1]     

[22] 尚勇, 陈至达.

带有摩擦的单边接触大变形问题的研究(I)——增量变分方程

. 应用数学和力学, 1989, 10(12): 1049-1058

[本文引用: 1]     

(Shang Yong, Chen Zhida.

On large deformation unilateral contact problem with friction(I)-Incremental variational equation

.Applied Mathematics and Mechanics, 1989, 10(12): 1049-1058 (in Chinese))

[本文引用: 1]     

[23] 尚勇, 陈至达.

带有摩擦的单边接触大变形问题的研究(II)——非线性有限元解及应用

. 应用数学和力学, 1990, 11(1): 1-11

[本文引用: 1]     

(Shang Yong, Chen Zhida.

On large deformation unilateral contact problem with friction(II)-nonlinear finite element technique and its application

.Applied Mathematics and Mechanics, 1990, 11(1): 1-11 (in Chinese))

[本文引用: 1]     

[24] 李平. 非线性大变形有限元分析的更新拖带坐标法及其应用. [博士论文]. 北京: 中国矿业大学(北京), 1991

[本文引用: 1]     

(Li Ping. The updated co-moving coordinate formulation for the nonlinear lager deformation finite element analysis and application. [PhD Thesis]. Beijing: China University of Mining and Technology (Beijing), 1991 (in Chinese))

[本文引用: 1]     

[25] Li P, Chen ZD.

The updated co-moving coordinate formulation of continuum mechanics based on the SR decomposition theorem

.Computer Methods in Applied Mechanics and Engineering, 1994, 114(1-2): 21-34

[本文引用: 4]     

[26] 高立堂, 李晓东, 陈礼刚.

火灾下钢筋混凝土板的热弹塑性有限元分析——基于S-R分解原理(II: 算例分析)

. 固体力学学报, 2006, 27: 138-142

[本文引用: 1]     

(Gao Litang, Li Xiaodong, Chen Ligang, et al.

Thermal-elastic-plastic finite element analysis of reinforced slabs under fire-based on S-R decomposition theorem(II: Analysis of examples)

.Acta Mechanica Solida Sinica, 2006, 27: 138-142 (in Chinese))

[本文引用: 1]     

[27] 高立堂, 宋玉普, 董毓利.

火灾下钢筋混凝土板的热弹塑性有限元分析——基于S-R分解原理(I: 理论)

. 计算力学学报, 2007, 24(1): 86-90

[本文引用: 1]     

(Gao Litang, Song Yupu, Dong Yuli.

Thermal-elastic-plastic finite element analysis of reinforced slabs under fire-based on S-R decomposition theorem(I: Theories)

.Chinese Journal of Computational Mechanics, 2007, 24(1): 86-90 (in Chinese))

[本文引用: 1]     

[28] 何满潮, 郭宏云, 陈新.

基于和分解有限变形力学理论的深部软岩巷道开挖大变形数值模拟分析

. 岩石力学与工程学报, 2010, 29(s2): 4050-4055

[本文引用: 1]     

(He Manchao, Guo Hongyun, Chen Xin, et al.

Numerical simulation analysis of large deformation of deep soft rock engineering based on solar decomposition theorem

.Chinese Journal of Rock Mechanics and Engineering, 2010, 29(s2): 4050-4055 (in Chinese))

[本文引用: 1]     

[29] 谢和平, 陈至达.

非线性大变形边界元法探讨

. 中国矿业大学学报, 1986(4): 101-109

[本文引用: 1]     

(Xie Heping, Chen Zhida.

The boundary element method of nolinear large deformation

.Journal of China University of Mining & Technology, 1986(4): 101-109 (in Chinese))

[本文引用: 1]     

[30] 谢和平.

非线性大变形问题的边界元分析法

. 应用数学和力学, 1988, 9(12): 1087-1096

[本文引用: 1]     

(Xie Heping.

Analysis of nonlinear large deformation problems by boundary element method

.Applied Mathematics and Mechanics, 1988, 9(12): 1087-1096 (in Chinese))

[本文引用: 1]     

[31] Gao Y, Feng G, Yeung MCR.

Modification of discontinuous deformation analysis method based on finite deformation theory

.Chinese Journal of Rock Mechanics and Engineering, 2011, 30(11): 2360-2365

[本文引用: 1]     

[32] 罗丹. 基于S-R和分解定理的几何非线性问题的无网格Galerkin法分析. [硕士论文]. 长沙: 湖南大学, 2011

[本文引用: 5]     

(Luo Dan.Based on S-R decomposition theorem analysis of element free Galerkin method on geometric nonlinear problems. [Master Thesis]. Changsha: Hunan University, 2011 (in Chinese))

[本文引用: 5]     

[33] 陈芳祖, 罗丹.

基于S-R和分解定理的无网格Galerkin法求解几何非线性问题

. 湖南大学学报(自科版), 2012, 39(1): 42-46

[本文引用: 2]     

(Chen Fangzu, Luo Dan.

Element free Galerkin method for geometrically nonlinear problems based on the S-R decomposition theorem

.Journal of Hunan Unviersity (Natural Sciences), 2012, 39(1): 42-46 (in Chinese))

[本文引用: 2]     

[34] 宋彦琦, 郝亮钧, 李向上.

基于S-R和分解定理的几何非线性问题的数值计算分析

. 应用数学和力学, 2017, 38(9):1029-1040

[本文引用: 3]     

(Song Yanqi, Hao Liangjun, Li Xiangshang.

Numerical analysis geometrically nonlinear problems based on The S-R decompostion theorem

.Applied Mathematics and Mechanics, 2017, 38(9): 1029-1040 (in Chinese))

[本文引用: 3]     

[35] Fan H, Zheng H, Zhao J.

Discontinuous deformation analysis based on strain-rotation decomposition

.International Journal of Rock Mechanics & Mining Sciences, 2017, 92: 19-29

[本文引用: 2]     

[36] Fan H, Zheng H, Zhao J.

Three-dimensional discontinuous deformation analysis based on strain-rotation decomposition

.Computers & Geotechnics, 2018, 95: 191-210

[本文引用: 3]     

[37] Bathe KJ, Ramm E, Wilson EL.

Finite element formulations for large deformation dynamic analysis

.International Journal for Numerical Methods in Engineering, 2010, 9(2): 353-386

[本文引用: 2]     

[38] Shen HS.

Large deflection of composite laminated plates under transverse and in-plane loads and resting on elastic foundations

.Composite Structures, 1999, 45(2): 115-123

[本文引用: 8]     

[39] Levy S.Bending of rectangular plates with large deflections. NACA Tech, Note No. 846, 1942

[本文引用: 6]     

[40] Gorji M.

On large deflection of symmetric composite plates under static loading

.Journal of Mechanical Engineering Science, 1986, 200(1): 13-19

[本文引用: 6]     

[41] 刘人怀, 薛江红.

复合材料层合板壳非线性力学的研究进展

. 力学学报, 2017, 49(3): 487-506

[本文引用: 1]     

(Liu Renhuai, Xue Jianghong.

Development of nonlinear mechanics for laminated composite

.Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(3): 487-506 (in Chinese))

[本文引用: 1]     

[42] Zaghloul SA, Kennedy JB.

Nonlinear behavior of symmetrically laminated plates

.Journal of Applied Mechanics, 1975, 42(1): 234-236

[本文引用: 2]     

/