数值分析中粘弹材料剪切模量松弛函数的拟合方法.docVIP

数值分析中粘弹材料剪切模量松弛函数的拟合方法 数值分析中粘弹材料剪切模量松弛函数的拟合方法 FITTING METHOD OF SHEAR MODULUS RELAXATION FUNCTION FOR VISCOELASTIC MATERIALS IN NUMERICAL ANALYSIS 帅词俊 　段吉安　　王　炯 (中南大学机电工程学院,长沙410083) SHUAI CiJun　DUAN JiAn　WANG Jiong (College ofMechanical and Electrical Engineering,Central South University,Changsha410083,China) 摘要　粘弹材料在进行有限元分析时,应先将经验或实验获得的剪切模量松弛函数拟合为相应的粘弹物理模型。 文中针对粘弹材料的经验剪切模量松弛函数Kohlrausch-William-Watts (KWW)方程,基于Matlab软件,利用线性方程组插 值求解法,实现对粘弹材料的广义Maxwell模型拟合,得出用于粘弹有限元分析的Maxwell模型参数的计算方法。并以1 050℃的粘弹光纤玻璃为例,给出一组计算实例。 关键词　粘弹材料　剪切模量松弛函数　Maxwell模型　参数拟合 中图分类号　O345　TB115 Abstract　While adopting finite elementmethod to analyze the stress and strain of viscoelastic material, the data of relaxationmod- ulus that is derived from experiments or experiences should be fitted for its corresponding viscoelastic physical model·Based on the Mat- lab software, experimental shear modulus relaxation function of viscoelastic materials is fitted to the generalizedMaxwell model by means of interpolation method of linear equation group·And have drawn the Maxwell model parameters computing method used for viscoelastic finite element analysis·Taking the optic fibre glass of 1 050℃as an example, one group of data is provided· Key words　Viscoelastic material;Shear modulus relaxation function;Maxwell model;Parameter fitting Corresponding author:SHUAI CiJun,E-mail:shuaicijun@163·com,Tel:+86-731-8836858 The project supported by the National Natural Science Foundation of China(No. Manuscript received in revised form 1　引言 粘弹材料已广泛应用于通信、建筑、医学、航天等 各个领域。对粘弹材料力学属性的描述有不同的物理 模型和方法。因广义Maxwell粘弹物理模型能够准确 地解释和描述粘弹材料的应力应变本构关系,所以在 现有数值分析软件中(如Ansys、MARC)都采用广义 Maxwell模型。 在粘弹材料的本构方程中剪切松弛模量是材料的 重要粘弹属性参数,Maxwell剪切松弛模量表达式是一 组松弛时间离散分布的指数加权的和[1~3]。要得到某 种材料的剪切松弛模量表达式,必须得到该材料的松 弛时间分布以及相应的加权系数。这些参数被称为 Maxwell模型参数。这此模型参数很难从实验或经验 中直接获得。 从实验或经验中获得的粘弹材料松弛模量表达式 常用扩展指数形式的Kohlrausch-William-Watts function (KWW函数)表达,可以通过查找某种材料的KWW指 数扩展因子和粘度及初始剪切模量大小确定材料的 KWW剪切模量松弛函数表达式。在进行数值分析前, 必需把KWW函数拟合为相应的广义Maxwel