粘弹性固体精细积分有限元算法.pdf

第21卷第1期 　 计 算 力 学 学 报 　 V o l. 2 1 , N o . 1 2004 年 2 月 Ch inese Journa l of Com puta tiona l M echan ics February 2004 文章编号:(2004) 粘弹性固体的精细积分有限元算法 1 2 3 刘靖华 , 　范益群 , 　钟万勰 ( 1. 上海交通大学 建筑工程与力学学院, 上海 200240; 2. 上海隧道工程与轨道交通设计研究院, 上海 200070; 3. 大连理工大学 工程力学系, 辽宁 大连, 116024) 摘　要: 粘弹性固体本构方程的数学表达式分为微分型和积分型两种, 其数值求解主要是时域上离散计算。文中 从微分型表达式出发导出其状态空间方程的数学表达式, 通过严格推导论证了它与微、积分型表达式的等价性; 引入状态空间方程, 从而利用精细积分格式来求解粘弹性固体本构方程; 给出了粘弹性固体本构方程的精细积 分有限元算法, 为求解粘弹性固体本构方程的数值解提供了一个新的途径, 具有计算简便, 求解精度高等优点。 关键词: 粘弹性固体本构方程; 状态空间方程; 精细积分有限元 中图分类号: O 241. 8　　　文献标识码: A N 步划分为 2 个精细步长, 其精度可以达到计算机 1　引　言 的精度, 具有理论简单, 编程计算容易, 计算精度高 一般而言, 粘弹性固体本构方程数学表达式分 等优点。最后, 本文由有限元初应变法出发, 形成一 为积分型和微分型, 其数值求解, 主要是时域上的 套求解粘弹性问题的精细积分有限元体系。 离散计算[ 13 ] 。对积分型的粘弹性体本构方程而言, 2　粘弹性固体本构方程的状态空间 若直接采用数值积分, 则为求算由整个应力历史所 表示及其与微、积分型本构 描述的积分项, 一般必须保存每个离散时刻的应力 等大量信息, 如果在空间也做大规模的数值离散, 方程间的转换 将给计算带来很大困难。对于微分型的粘弹性体本 　　不失一般性, 普遍线性粘弹性体的应力、应变 构方程, 可以通过有限差分而离散为多步递推公 速率关系为[ 12 ] 式, 但如微分方程的阶数较高, 则也会出现信息贮 (n) (n - 1) 　　　 + a + …+ a = 1 n 存量过大的问题。对于某些积分型的本构方程, 可 (n) (n- 1) ( )