以电网络为模型的粘弹性波动有限元分析法.pdfVIP

以电网络为模型的粘弹牲波动 有限元分析法 蔡抽宝0 俞集辉0 杨瑞芳③ ①重庆大学电气工程系，②重庆大学生物工程研究院，重庆：400044 摘要：阐迷本文作者已证明了的两个重要蛄诧：任意复杂的广业曲线性粘弹性模 型均可等效为广义的yoigt链；线性粘弹性介质三雄本构关系的力学模型孵络与六端口 Rc电网络互为比拟．藉助于电网搏理论和滤波方法将弹性动力擘育浪觅分析的理论和 方法推广到了粘弹性波动的有限元分析中．导出了适用于一般的巍挂粘弹性三堆波动 有限元分析的基本公式及其计算方法．多方法，多角度地翊明彳奉史曲研克为基于线 性粘弹性模型的生物组织和工程材料的动力学分析导出?宽适应的、系砘韵和高精度 的有限元快速分析法． 关键词：电路理论，粘弹挫理论，有限元 KEYWORbS：ElectricalCircuitTheory，Viscoelastici“1h—j碍，Pi湖teElement l引 言 粘弹性理论是由电磁场理论的创始者J．C．Maxweij于)857车奠定的“．后来，在 I；atanson、 几乎停滞了五十年之后的两个世纪之交又由Voigt、KeI,Ⅶ、Duhem、 Zaremba和其他一些科学家进行了研究。值得指出的是，我国的著名科学家钱学森旧1在 粘弹性对应原理的研究方面作出了杰出的贡献。在生物软组织辐弹性力卓的理论研究和 科学实验方面，著名的华裔学者冯元桢【j1先生作出了杰出的贡献。 近年来，随着航空和航天技术的发展，新型的赢分子复合材斟、钻弹’{：兰咀尼材料和 智能型阻尼材料得到了进一步的发展。由于在固体燃料工业中存在一些急待解决的问 题，世界上有关国家在这个领域中进行了广泛韵实验工作，惶是在这方面的大量工作部 是保密的，并末出版．出版了的那些工作通常局限于某些范围F?，例如，_}lI静力，一维 或二维的，温度和时间区间是育限的等等14J。对于各向异性牯弹滗物质还发育台理的实 验方法。或许就是这样的一些实验和模型的问题，以及数学处理上的困难，牯弹性波动 数值分析事实上还没有得到充分发展，落后子工程应用和科学研究的迫切需要。 国内外已发表的粘弹l生动力学响应有限元分析法均采用增量法．并采用总体刚度矩 阵形成线性或非线性方程组的求解模式．介质的本构方程尚局限于线睦冬向同性粘弹 性，时间步长的选择较困难．这就限制了计算效率的提高和算法的适用苞嗣。线性粘弹 性工程材料的动力学分析，在边界元分忻法的研究方面己见开创性的文献报导p】，存在 的问题是，在复频域里用边界元分析法求出粘弹性动力学响应的数值解，再用拉普拉斯 变换数值反演法求得时域数值解，这样的算法在计算的复杂性上存在无法避免的问题， 即存在大量的．在数值计算中最违忌的复数运算，从而导致在计算效率上抵消了边界元 法的优越性，时间步长的选择较困难，故不能在计算开始时保l正边界元解的收敛性。 本文旨在鼯助于电网络理论和滤波方法，导出不必建立刚度矩阵的总体精点位移盲 接递推有限元新算法，流变介质的复杂程度可以从线l生备向同炷逐步放宽到赎端各向异 |生，描述生物组织或工程材料力学性能的广义线性桔弹性模型可以是两稃，基本元件与= 种基本模型的任意组合。时间步硅的选择简洁明了，可确保计算结果收敛一稳态解，即 算法是十分稳健的。 2电网络与广义的线性粘弹性模型 依据相似性原理，推广应用电网络的对偶原理，将力学模型网络转换成电网络，并 应用电网络的综合理论，本文作者已导出下述两个结论”j． 2·l i 广义的VOgt链 结论1：由线性弹簧元件，阻尼器与三种基本模型组成的任意复杂的广义的线性粘 弹性模型，均可等效为图1所示的与一端口RC电网络互为比拟的广义的Voigt链， En — 罔1广义的Voigt链 J面系统函数的一般形式如式(1)所示。 1 1 n H。