nbhkdz.com冰点文库

对于线性粘弹性沥青混凝土的卷积积分的数值计算2

时间:2013-03-04


对于线性粘弹性沥青混凝土的卷积积分的数值计算

By Sungho Mun*
摘要 此研究关注表示电流的卷积积分的数值制剂的发展,它代表以前所有的历史沥青混凝土的响应时间。为了解 决涉及线性粘弹性的初始/边界值问题, 本构方程, , 表现在沥青混凝土从频率扫描测试中的卷积积分形式, 从而得到松弛模量和数值配制成的增量。松弛模量是线性弹簧和阻尼器组成的力学模型。机械维谢尔模型 是由 Prony 用于积分基础公式的上的线性粘弹性松弛模量的数值而制定的。对于评价目的,是涉及的十字 头张力(或恒定的应变速率单调加载) ,锯齿形状加载,和一个光束前端负荷下的数值例子纳入数值积分的 验证。 关键词:卷积积分,粘弹性,维谢尔模型,数值配方,Prony 系列 1.介绍 应用程序的初始/边界值问题的粘弹性行为已经被文献报道。粘弹性本构模型的数值积分公式材料已经在包 括泰勒等人的作品中讨论过。 (1970 年) 措赫尔等人 , (1997 年) Kaliske 和 Rothert 1997 年) 潘· , ( , 艾 哈迈德(1998 年),SIMO 和休斯(1998 年),Hiterhoelzl(2000 年)。基本上,该方法已被纳入有限 元模拟。总之,SIMO 和 Hughes(1998)提出的方法的递推公式为一体,采用了类似的公式到目前这个 数值积分公式。对于粘弹性的基本理论,一些研究人员的作品可以咨询, ,如 Schapery(1974 年) ,轮渡 (1980) ,克里斯坦森(1982) ,和 Tschoegl(1989) 。粘弹性材料区别于广义的材料,这些材料的弹性描 述了弹性和抗老化材料的线性粘弹行为。这个理论是一个洞察材料普遍的响应行为随时间变化的。在本文 中,一个机械的本构方程开发并接受机械载荷的固体的卷积积分的松弛模量由 Prony 系列取代。这里给出 这一提法的的时间 historydependent 变形和负载被施加到。基于 Prony 系列,依赖于时间的线性正交各向 异性粘弹性本构关系,制定厮磨等人所提出的递推公式。此外,计算的切线模量从增量的应变和应力之间 的关系。 2.力学模型 粘弹性材料的松弛模量是用来定义的应激反应,从各自的张量分量的应变历史输入如下:

其中, IJ=应激反应, Eijk=松弛模量 =应变输入。 KL 一种粘弹性材料配方下面的开发相关的松弛模量的力学模型组成的弹簧和阻尼。首先,被给定为一个一维 的分析基本方法。其次,相应的三维分析描述被跟随。此外,被认为在制定 thermorheologically 简单的材 料。这个假设是随后所有松弛曲线可以叠加水平平移的方式沿时间轴以形成一个单一的主曲线。 因此,热 机械耦合的过程被简化,并且不考虑温度的变化(例如,等温过程) 。在下面的广义模型,维谢尔模型,进 行了讨论。如图 1 中所示,左弹簧中的应力,&= Maxwell 元件的弹簧与一个的缓冲器相结合的应力&,来自微分方程:

图1

维谢尔型号

其中, =压力 M ¥m=粘度系数 EM = 在第 m 个内的松弛模量。 维歇特模型(或广义的 Maxwell 模型)上的总应力是通过以下方式获得的求和形式:

在解决上述微分方程,拉普拉斯变换的方法是有用的。其方法是转换成差分方程的代数方程,这使得获得 拉普拉斯变换的解决方案(例如,反演问题) 。拉普拉斯转换式。 (2)?(4) ,然后消除应力 给出 m

其中,E=卡森变换的松弛弹性模量(例如 E,SE,,和第 n 个麦克斯韦元素的弛豫时间定义:P=n/E ) 这可以被转换成时域模量从拉普拉斯反演松弛模量在 s 域方面。

上面的方程被称为 Prony 系列相应的维谢尔模型(或广义 Maxwell 模型) 。一维 Prony 级数可以配制的三 维表达。

3。本构方程的数值积分 在这项工作中所考虑的粘弹性材料本构方程推导出式的松弛模量的 Prony 系列。 (7) 。下面的卷积积分 在一般的各向异性应力和应变之间的关系:

为了轻松处理多个指数,上述项建议的是潘和艾哈迈德(1998) 。

这个定义可以用来改写的 EQ。 (8)作为

这里

&的定义,以制定一个简单的数值实现。此外,被构造的应力, ij其中由两个组件如瞬时弹性应力,Eijkl0 千升,和应力响应的基础上的粘性应变,&。 状态变量是正确的选择,目前的粘弹性时间制定整合过程中的一个重要方面。存储和更新的状态变量从一 个时间增量到下一个不知道他们的整个时间的历史。上述状态变量的选择:

上述状态变量可以被重新复发的方法“ (SIMO 和休斯,1998 年) 。下列标准的指数函数和对任何常量,DT 和 a。

从这样的说法,

状态变量的方程。 (12)可以是半分组,通过使用上述的可添加的属性。采取的最后形式,如下所示:

如式中所示。 (15) ,本发明的状态变量, :荷航 (tn +1 处) ,是从以前的变量, 荷航(TN) ,更新的时 间之前不知道时间的历史数据中,tn。对于进一步的数值算法,它被假定的菌株, 千升,与增量的线性变 化 。 此 外 , 该 公 式 可 以 重 写 相 对 于 共 同 的 指 数 项 。 在 这 种 情 况 下 ,

因此,本构法归纳如下:

4。切线模量 在准备用于导出的切线模量,下面的制剂是来自于两个不同的方法。然后,结果是被相互比较的一次微分 从方程。 (8) ,如下所示:

小增量的假设,

应力张量被认为是通过构的关系,式作为一个功能的应变。 (19) ,在所选择的状态变量,如 ?形式 (i)和 (ⅰ)n +1,和切线模量可以推导如下:

二 次 的 切 线 模 量 计 算 方 法 , 从 EQ 。 17 ) 的 基 础 上 它 相 对 于 ( tn +1 处 ) 区 别 。 因 此 , ( KL

如果使用的定义式被重写方程公式(9) ,上述切线模量相当于与第一种方法的结果。 5。验证数值计算 为了验证算法,选择实施例中,是已知的解析解。为了简化,onedimension 例讨论在实施例(1)和(2) , 和一个叠加技术采用解析来表示的应变历史。这种方法是基于假设的线性度。最后,实施例 3 有关的一个 简单的剪切试验。 实施例 1 是有关一个恒定应变速率单调加载(或恒定十字头应变速率拉伸)的输入。应变历史在当前时间 是: 其中,C 是一个常数的应变率,和 H(t)的亥维赛阶梯函数,H = 0 时,当 t <0 时,和 H =1 当 t> 0。 使用一阶导数的应变函数式(21) ,和由 Prony 级数,方程表示的松弛模量(6) ,压力可以从方程获得公 式(8) 。

其中 C 是 2。 对于普罗尼的弛豫时间相对于系数,精确的材料性能的沥青混凝土从扫频测试施加正弦加载沥青试样,并 嵌合在参考温度,25℃,得到示于表 1。图。图 2 示出数值模拟和分析的解决方案之间的比较的基础上的 恒定 strainrate 的单调加载,并且它指出增量的 0.01 和 0.1 倍。 锯齿形状的应变输入被选择作为实施例 2 使用相同的材料性能的 Prony 算法在图系数。 3。的实施例 2 的 解析解,可以得到从的线性叠加技术取决于对三个不同的区域,如 0<t <t1 时。对于为 0 <t<t1 时,方程。 (22)可用于与不同的常数 C 的值

图 2。恒应变速率单调加载

(25) 其中 C 是 0.5 的情况下的正斜率和-C 是在负,t1= 2.5 和 t2= 5 秒。

图 3 锯片齿形应变加载

图 3 显示了优良的协议之间的理论分析和数值的解决方案,甚至想到不同的增量。 要说明如何可以验证本算法的正交各向异性材料中,选择是基于实施例 3 的材料特性,这在表 2 中所示, 使用由 Poon 和 Ahmad(1998 年) ,和他们定义 Eijkl 作为平衡松弛模量和 eijkl0 使用一个单一的弛豫时 间,为玻璃状松弛模量。实施例 3 是一个简单的材料样品的剪切的情况下进行描述为施加的剪切应变历史 输入:

可以推导出的解析解如下:

术语的增量的时间,例如 0.01 和 0.1,计算出的数值解是在良好的协议与解析解,如图 4。 6。反演问题 不同于应变输入问题,应变输入和应激反应,式构关系,因为在所施加的应力的情况下,在试样中的反转 过程是必需的。 (1) ,是无效的。从式(18)和(19) , (tn +1 处)可以从以下的假设计算:

“-0”是指在纵坐标上的原点轴从非常小的变化因此,初始应变值是零。利用上述假设,对于小的增量方 程, (19)是可以改变变分形式。

最后 对于电流反转算法的验证,被认为是与前端负荷悬臂梁问题,对应于规定的松弛模量,它需要蠕变柔量, 。 基本上,输入从应力历史计算应变响应,使用上面描述的技术反演该光束具有长度 L,12m4 转动惯量的时 刻,I,1/12 的 m4 束进行的前端负荷: 其中 P0=1 兆帕和 t1= 1 秒。 据的弹性的解决方案中,位移,WL,由于前端负荷可以以下方式获得:

为了得到粘弹性体的溶液,施加弹性粘弹性对应原理。此对应原理表明静态弹性的解决方案可以被转换成 准静态粘弹性解决方案。 该过程涉及更换弹性模量的粘弹性能, 拉普拉斯变换和反相转换后的粘弹性变量, 得到时域解决方案。使用这个属性公式(32)被改变的拉普拉斯变换如下: 其中,D=1/E 和 D 是卡森变换的蠕变,符合(D`=SD) 。

图 5。提示荷载工况的梁 当方程的反转时,最终的时域解得到的粘弹性位移:

与解析解的数值解相比,它是必要蠕变柔量,D(t)的计算,从所描述的 Prony 系列,在表 1 中所示,在 松弛模量的方面。因此,使用下面的关系 D = 1 / E, ,蠕变柔可以得到。涉及的拉普拉斯变换和反演的蠕 变柔量的计算有关的问题的解决通过软件包 Mathematica 的。 如图 5 中所示, 数值解是增量时间, 0.1, 如 0.01,和 0.0001 增量的不同而不同,在压力 applied.The 为零表明,当前的解决方案收敛到解析解的压力 的持续时间。此外,该 最小递增的时间提供正确的斜率趋势,对分析解决方案,在所施加的压力历史上输入. 7。结论 基于普罗尼指数级数表示数值积分方法被提出并证明是非常有效和准确的。四个实施例中的每一个都被选 择,以验证本发明方法。结果发现是相当不敏感的增量在实施例 1 至 3 倍。但是,在实施例 4 中,自应力 时间历史输入需要的反演技术来计算应变响应是依赖于时间增量。反演问题最小增量值的使用产生非常精 确的结果。这一框架可以容易地并入通用的有限元法分析的数值积分公式。 参考文献 克里斯滕森,R.M. (1982)。粘弹性理论:介绍。学术出版公司,伦敦。 Ferry, J.D. (1980).(1980 年)。粘弹特性的聚合物。第 3 版,Wiley 出版社,纽约州。 。固体的推进剂中的有限元程序 ABAQUS 的 UMAT 根据“组织法” ,研究报告 Hinterhoelzl,R.(2000) 。 ”计算力学方法。 kaliske,M.和 Rothert 的,H.(1997)“的制定和执行情况的三维粘弹性小和有限应变。 19,第 228-239 页。 潘,H.和艾哈迈德,M.F. (1998 年)“A 材料各向异性点时间整合过程,热流变简单,粘弹性的固体。 。 ” 计算力学,第。 21,第 236-242 页。 。复合材料的粘弹行为与分析。第 4 章复合材料 2。 (的主编 Sendeckjy,G. P.) , Schapery,R.A. (1974) 科学出版社。

厮磨,J.C.和休斯,T.J.R. (1998 年) 。计算非弹性施普林格出版社,纽约州。 泰勒, 匹斯特, 。 “热力学分析的粘弹性固体” J.编号。 。 冰毒。 Engng, RL, KS, Goudrequ, (1970 年) GL 卷。2 期,页 45-59。 。唯象理论的线性粘弹行为。斯普林格,柏林。 Tschoegl,N.W. (1989) 措赫尔, :MA,树丛,SE,和艾伦,DH(1997 年)“A 三维有限元热黏弹性各向异性介质配方。。 J. 。 ” 冰毒。工程,卷。 40,第 2267 至 2288 年。 (9 月 8 日,2005/Accepted2006 年 4 月 27 日)

标题 沥青混凝土的粘弹性数值计算卷积积分 日志 KSCE 土木工程学报 第 10 卷,第 3 期,195 页-200 页 盖子日期 2006-05-01


赞助商链接