非连续变形计算力学模型的粘弹性分析方法_算例分析.docVIP

收稿日期:2002-02-28;修改稿收到日期:2002-09-23. 基金项目:国家自然科学基金大连理工大学青 年教师培养基金;国家重点基础研究发展规划(2002CB412708;教育部优秀青年教师奖励计划 资助项目. 作者简介:黎　勇(1970-,男,博士,博士后,讲师; 栾茂田(1962-,男,博士,教授,博士生导师; 冯夏庭(1964-,男,博士,教授,博士生导师;王泳嘉(1935-,男,教授,博士生导师. 第21卷第4期2004年8月　计算力学学报　 Chinese Journal of Computational Mechanics V ol.21,N o .4A ug ust 2004 文章编号:1007-4708(200303-0448-07 非连续变形计算力学模型的粘弹性分析方法 Ⅱ:算例分析 黎　勇 1,2,3 ,　栾茂田1,2,　冯夏庭3　,王泳嘉1,　叶祥记 1 (1.大连理工大学土木水利学院,辽宁大连116024; 2.大连理工大学海岸和近海工程国家重点实验室,辽宁大连116024; 3.东北大学资源与土木工程学院,辽宁沈阳110006 摘　要:应用作者们所提出的非连续变形计算力学模型粘弹性分析方法对具体的多体系统在动力外荷载作用下的响应进行了数值模拟,探讨了系统中物体的粘性阻尼对整个系统的变形、应力和接触应力等的影响。与通用软件A BAQ U S 的计算结果进行了比较研究,两种方法所得到的结果基本一致,但所发展的方法克服了A BA QU S 中接触应力分段均匀分布的缺陷,因而所得到的接触应力较为合理,同时也说明了所提出的方法为由刚体、弹性体和粘弹性体所构成的复杂多体系统的数值分析提供了有力的工具。 关键词:非连续变形;计算力学;粘弹性;接触力元;广义有限单元;相互作用;多体系统;界面;应力集中中图分类号:O 342;T U 4　　　文献标识码:A 1　引　言 随着科学技术的发展和工程建设的需要,数值分析与模拟由不同特性的物体所构成的复杂多体系统的静动力特性变得日趋重要。目前所采用的数值分析方法有:(1不考虑物体间的相互作用而把多体系统当作连续介质系统用有限元法和边界元法 [1-3] ;(2基于连续介质理论,通过设置界面单元 或弹簧来模拟物体间的相互作用而模拟和分析多体系统行为的方法如界面单元-有限元法和刚体弹簧模型等 [4-6] ;(3基于离散介质系统,在物体(或单 元间通过设置弹簧模拟多体系统的离散单元法和非连续变形分析方法等 [7-9] 。显然后者比较合理,但 是在物体间设置附加弹簧将影响系统的变形和应力,并不能直接给出物体间的相互作用力。作者们前文所发展的非连续变形计算力学模型粘弹性分析方法 [10] 能够克服这些缺陷,并能够分析由刚体、 弹性体和粘弹性体所构成的复杂多体系统的静动力响应。本文应用所提出的粘弹性分析方法对具体 的多体系统进行了数值分析,探讨了物体材料特性对系统的变形、应力和接触应力的影响,并与ABAQUS [11] 的计算结果进行了比较研究。 2　多体系统模型与算例条件 研究对象为如图1(a所示的由3个物体所构成的多体系统,其中物体A 为刚体且固定;物体B 为粘弹性体,其质量密度、弹性模量和泊松比分别为2.5t/m 3,10M Pa 和0.35;物体C 为弹性体,其质量密度、弹性模量和泊松比分别为3.0t/m 3、1000M Pa 和0.1。物体A 与B 间的接触界面CS1和物体B 与C 间的接触界面CS2分别为完全粗糙和光滑(即相当于连续和非连续接触界面。该系统所受到的外荷载参见图1(a,并按图2所示的正弦波变化,计算中的广义有限单元剖分参见图1(b,共有291个广义结点,258个广义有限单元,34个接触对。为了探讨物体B 的粘性阻尼对该多体系统动力响应的影响,采用了两种工况,在工况1中物体B 的变形阻尼系数和运动阻尼系数分别为0.01和1;在工况2中物体B 的变形阻尼系数和运动阻尼系数分别为0.1和10。分别来用作者们所开发的非连续变形计算力学模型的粘弹性分析程序和通用有限元分析软件ABAQU S 进行了计算,对计算结果进行比较了对比分析研究如下。 图1　多体系统的构形及单元剖分 Fig.1　Configuration and element s ubdivis ion of th e multi-body sys tem 图2　水平荷载时程曲线 Fig.2　T ime curve of horiz on tal load 图3　不同工况下物体B 内a 点处的变形时程曲线Fig.3　T ime h istories of displacements at Point a in Body B 3　数值计算与比较研究 这里通过数值计算分析比较了系