›› 2012, Vol. 33 ›› Issue (11): 3497-3502.

• 数值分析 • 上一篇    下一篇

四节点矩形弹簧元及其特性研究

张青波,李世海,冯 春   

  1. 中国科学院力学研究所,北京 100190
  • 收稿日期:2011-08-24 出版日期:2012-11-12 发布日期:2012-11-14
  • 通讯作者: 李世海,男,1958年生,博士,研究员,主要从事非连续介质力学数值计算方法及滑坡灾害防治中的关键力学问题等方面的研究工作。E-mail: shli@imech.ac.cn E-mail: zhangqingbo@imech.ac.cn
  • 作者简介:张青波,男,1986年生,博士研究生,研究方向为地震作用下边坡的破坏机理及数值模拟
  • 基金资助:

    国家重点基础研究发展计划(973)项目(No. 2010CB731506);国家自然科学基金项目(No. 11002146)

Study of four-node rectangular spring element and its properties

ZHANG Qing-bo,LI Shi-hai,FENG Chun   

  1. Institute of Mechanics, Chinese Academy of Sciences, Beijing 100190, China
  • Received:2011-08-24 Online:2012-11-12 Published:2012-11-14

摘要: 弹簧元法是一种将单元离散为一系列弹簧的数值计算方法。不同的单元具有不同的离散方式,确定相应的离散弹簧的刚度系数表达式是弹簧元法的关键。将四节点矩形单元离散为6个基本弹簧,每个基本弹簧包括法向弹簧和切向弹簧两个派生弹簧,并用泊松弹簧和纯剪弹簧描述单元的泊松效应和剪切效应,用有限元的单元刚度矩阵标定各弹簧的刚度系数,实现了一种四节点矩形弹簧元的构造形式。该单元的同类弹簧具有相同的表达形式。法向与切向弹簧的刚度表达式中分别含有法向和切向弹簧刚度待定系数。通过改变待定系数的值可使该单元分别对应于有限元的常应变、双线性及Wilson非协调单元。将上述弹簧元方法进行理论推导,并应用于基于连续介质的离散单元法(CDEM)的核心计算进行简单算例验证,证明了提出方法的正确性。通过以上研究发现,四节点矩形弹簧元有以下特点:对于相同问题,不同单元有不同的计算精度;对于梁弯曲问题,应用该单元可显著提高离散单元法的求解精度;改变待定系数的值,可得到更高或者更低精度的单元。

关键词: 四节点矩形单元, 弹簧元法, 有限元法, 弹簧刚度

Abstract: The spring element method (SEM) is a numerical method that uses a spring system to describe an element. Different elements can be described as different spring systems; and the definition of the spring stiffness expressions in the systems is the key point of the spring element method. The four-node rectangular element is described by 6 basic springs, each of which contains two derived springs: normal spring and tangential spring. Poisson spring and pure shear spring are used to describe Poisson and shear effects of the element. Thus a four-node rectangular spring element is presented. Compared with the element stiffness matrix of finite element method, the stiffness expression of each spring is obtained. Springs of the same kind have the same expressions. The stiffness expressions of the normal and tangential springs have corresponding coefficients to be decided. By varying the coefficients in the stiffness expressions of springs, expressions of constant strain, bilinear or Wilson incompatible finite element are achieved by this element. The accuracy of the SEM is verified by theoretical derivation; and this method is applied to the continuum-based discrete element method (CDEM) for case verification. The features of the four-node rectangular spring element are as follows. Different accuracies can be found in different elements. This element can significantly improve the accuracy of the bending problem of beam. Elements with different accuracies can be achieved by using different coefficients.

Key words: four-node rectangular element, spring element method (SEM), finite element method (FEM), spring stiffness

中图分类号: 

  • O 241
[1] 曹洪, 胡瑶, 骆冠勇. 滤管两端均不在含水层层面的承压不 完整井近似计算方法研究[J]. 岩土力学, 2019, 40(7): 2774-2780.
[2] 王翔南, 李全明, 于玉贞, 喻葭临, 吕禾, . 基于扩展有限元法对土体滑坡破坏过程的模拟[J]. 岩土力学, 2019, 40(6): 2435-2442.
[3] 郑安兴, 罗先启, 陈振华, . 基于扩展有限元法的岩体水力劈裂耦合模型[J]. 岩土力学, 2019, 40(2): 799-808.
[4] 王冬勇, 陈曦, 于玉贞, 吕彦楠, . 基于二阶锥规划有限元增量加载法的条形浅基础极限承载力分析[J]. 岩土力学, 2019, 40(12): 4890-4896.
[5] 郑安兴,罗先启,. 危岩水力劈裂分析的扩展有限元法[J]. , 2018, 39(9): 3461-3468.
[6] 宋 佳,古 泉,许成顺,杜修力,. 饱和土动力方程全显式有限元法在 OpenSees中的实现与应用[J]. , 2018, 39(9): 3477-3485.
[7] 宋 佳,杜修力,许成顺,孙宝印,. 饱和土场地-桩基-地上结构体系的地震响应研究[J]. , 2018, 39(8): 3061-3070.
[8] 李 宁,郭双枫,姚显春,. 再论岩质高边坡稳定性分析方法[J]. , 2018, 39(2): 397-406.
[9] 罗先启,郑安兴,. 岩体裂隙模拟的扩展有限元法应用研究[J]. , 2018, 39(2): 728-734.
[10] 刘忠玉, 张家超, 郑占垒, 关 聪. 考虑Hansbo渗流的二维Biot固结有限元分析[J]. 岩土力学, 2018, 39(12): 4617-4626.
[11] 刘振平,杜根明,蔡 洁,周 凡,刘 建,卞 康,. 基于MeshPy的3DGIS与三维有限元数值计算无缝耦合方法[J]. , 2018, 39(10): 3841-3852.
[12] 邹德高,刘 锁,陈 楷,孔宪京,余 翔,. 基于四叉树网格和多边形比例边界有限元方法的岩土工程地震响应非线性静动力分析[J]. , 2017, 38(S2): 33-40.
[13] 林良庆,陈福全. 大直径贝诺特桩套管上拔机制分析[J]. , 2017, 38(8): 2385-2394.
[14] 闫富有,常 键,刘忠玉. 黏弹-双曲线Drucker-Prager塑性模型应力更新隐式算法[J]. , 2017, 38(6): 1797-1804.
[15] 刘振平,刘 建,何雨微,贺怀建,卞 康,. 3D GIS与有限元模拟无缝耦合方法及其在隧道工程中的应用研究[J]. , 2017, 38(3): 866-874.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
[1] 刘小文,常立君,胡小荣. 非饱和红土基质吸力与含水率及密度关系试验研究[J]. , 2009, 30(11): 3302 -3306 .
[2] 黄建华,宋二祥. 大型锚碇基础围护工程冻结帷幕力学性态研究[J]. , 2009, 30(11): 3372 -3378 .
[3] 姚仰平,冯 兴,黄 祥,李春亮. UH模型在有限元分析中的应用[J]. , 2010, 31(1): 237 -245 .
[4] 王朝阳,许 强,倪万魁. 原状黄土CT试验中应力-应变关系的研究[J]. , 2010, 31(2): 387 -391 .
[5] 高 阳,张庆松,徐帮树,李 伟. 海底采煤顶板支承压力分布规律与影响因素研究[J]. , 2010, 31(4): 1309 -1313 .
[6] 刘泉声,胡云华,刘 滨. 基于试验的花岗岩渐进破坏本构模型研究[J]. , 2009, 30(2): 289 -296 .
[7] 赵尚毅,郑颖人,李安洪,邱文平,唐晓松,徐 俊. 多排埋入式抗滑桩在武隆县政府滑坡中的应用[J]. , 2009, 30(S1): 160 -164 .
[8] 周 扬,周国庆. 土体一维冻结问题温度场半解析解[J]. , 2011, 32(S1): 309 -0313 .
[9] 陈保国 ,孙金山 ,张 磊. 上埋式钢筋混凝土拱涵受力特性及地基处理研究[J]. , 2011, 32(5): 1500 -1506 .
[10] 邢万波 ,周 钟 ,唐忠敏 ,孙 钢. 基于ν-SVR和改进PSO算法的反分析方法及应用[J]. , 2009, 30(S2): 540 -546 .