计算机辅助工程分析中模拟材料老化效应的方法和系统的制作方法

专利名称：计算机辅助工程分析中模拟材料老化效应的方法和系统的制作方法
技术领域：
本发明涉及时间流变材料(例如，橡胶、泡沫等)，更具体地说，涉及一种用于在计 算机辅助工程分析(例如，有限元分析(FEA)、无网格分析、有限差分析等)中模拟时间流变 材料的老化效应的方法和系统，用于辅助用户(例如，工程师、科学家等)在改进和设计至 少部分包含时间流变材料的工程结构时作出决策。
背景技术：
橡胶之类的材料(例如，人造橡胶、泡沫)多年来已经被用于各种行业(例如，汽 车、航空、建筑等)的工程结构的许多部分中。橡胶之类的材料通常具有时间流变(流变和 老化)的特性。时间流变特性的幅度随时间改变，经常很大使得它们不能被忽视。随着现代计算机技术的进步，通常使用计算机辅助工程分析来设计和改进工程结 构(例如，汽车、飞机等)。为了捕捉包括橡胶之类的材料的老化效应在内的结构响应，需要 确定材料特性的数值表示，然后将它结合到计算机辅助工程分析的应用模块或者软件中。 一种数值表示形式是材料本构方程，它与材料应力-应变或者力-位移关系有关。通常，这 通过对感兴趣的材料进行样本测试、接下来将测试结果与特定的本构方程相关联来完成。 但是，现有技术方法中所使用的橡胶之类的材料的本构方程通常是不合适的。有一些是基 于对不能压缩的弹性特性或者某些特定情况(例如neo-Hookean，Mooney-Rivlin)的假设 的。其它(例如，Ogden等式)可包括不能压缩的和可压缩的材料特性，但是没有老化效应。现有技术方法的缺点使模拟泡沫胶(可压缩和黏弹性材料)的特性产生了问题， 这种泡沫胶广泛地用于汽车行业。特别是，设计汽车的一个行业标准是在碰撞或者汽车撞 击事件时确保汽车所有者的安全。为了确保满足这样的要求，汽车制造商需要对每个汽车 模型的原形进行物理碰撞测试。在某些情况下，需要进行不同场景的多次测试。物理撞击 测试的执行不仅很昂贵，并且还具有它自己的一些困难和挑战(例如，可测量性、精确度、 可靠性、可重复性等)。因此，撞击事件的计算机模拟(也就是，计算机辅助工程分析)已经 被广泛用于替代物理原形撞击测试或者至少将物理原形撞击测试的数量最小化。在物理撞击测试中，一个或多个撞击假人被放置在汽车中，以代表真人所有者 (也就是，司机或者司机加乘客)，以研究安全要求。撞击假人由通常使用泡沫材料的多个 部件(例如，头、躯干、四肢等)构成。由于撞击假人通常在物理撞击测试中被毁坏，之后被 修复并重复用于另一次测试，撞击假人的不同部分可能包含不同老化程度的泡沫材料。为 了更准确地模拟这样的情况(也就是，在一个撞击假人中具有多种老化程度的泡沫材料)， 计算机辅助工程应用软件需要能够采用包含老化效应的材料本构方程来计算结构响应。因 此，期望有一种用于在有限元分析中模拟时间流变材料的老化效应的方法和系统。

发明内容
本发明公开了一种用于在计算机辅助工程分析中模拟时间流变材料的材料老化 效应的方法和系统。根据一方面，执行一组材料特性测试，以获得感兴趣的时间流变材料的材料特性。这一组材料特性测试包括对多个不同老化程度的材料样本重复基本类似的测 试，该老化程度是例如天、周、月等。可以通过单轴向、双轴向、或者弯曲测试程序来执行每 次测试。针对被测试的样本，可以得到一系列松弛测试数据。在单轴向测试程序中，可以执 行拉伸或者压缩松弛测试，在不同的时间测量轴向应力。在双轴向测试中，可使用圆形(例 如，圆形、卵圆形、或者椭圆形)隔膜的涨大。保持隔膜的变形恒定，并测量涨大压力(也就 是，松弛测试数据)。在松弛测试数据表或者图中以对数-对数标度来组织或者绘制每次材 料特性测试的一系列松弛测试数据。通常，数据图表包括多个系列的松弛测试数据，每个系 列表示感兴趣的时间流变材料的特定老化程度。在本发明的另一方面，创建一种新的材料本构方程来表示包含材料老化效应的时 间流变材料。所述新的材料本构方程被配置为包括了现有技术方法(例如neo-Hookean， Mooney-Rivlin,Ogden等)中使用的现有本构方程的通用和统一的公式。在新的材料本构 方程中包括一对与时间相关的材料老化效应参数。第一参数与材料的松弛测试数据相关， 表示为材料老化效应参数a，而第二参数与时间相关，表示为材料老化效应参数3。这两 个参数可以由松弛测试数据在材料特性测试中获得的对数-对数标度表中的垂直和水平 移位以及匹配来确定。表中的每个系列的松弛测试数据表示所获得的特定老化程度样本的 结果。每对数据的移位和匹配得到一组特定的材料老化效应参数。使用其中配置有这样的材料本构方程的计算机辅助工程分析(CAE)应用模块，时 间推进模拟中的时间流变材料特性的影响或者行为便被包含在内。通常，获得的一组第一 和第二与时间相关的材料老化效应参数被输入到其上安装有CAE应用模块的计算机内。通过以下结合附图对具体实施方式
的详细描述，本发明的其他目的、特征和优点 将会变得显而易见。


参照以下的描述、后附的权利要求和附图，将会更好地理解本发明的这些和其它 特征、方面和优点，其中图1A是单轴向拉伸状况下的人造橡胶棒样本的示意图；图1B是双轴向拉伸状态下的人造橡胶棒样本的示意图；图2A是根据本发明的一个实施例的一系列示范性模拟松弛测试数据的示意图；图2B是图2A的一系列模拟松弛测试数据的总曲线；图3是根据本发明的实施例在单轴向测试中获得的一组示范性松弛测试数据的 示意图；图4A是根据本发明的实施例使用双轴向测试程序从一组材料特性测试中得到的 示范性松弛测试数据的示意图；图4B是图4A的松弛测试数据移位和匹配后的结果示意图，该移位和匹配用于确 定一组第一和第二与时间相关的材料老化效应参数；图5A和5B共同示出了根据本发明的一个实施例的包含第一和第二与时间相关的 材料老化效应参数的时间流变材料本构方程；图6是根据本发明的实施例用于在有限元分析中模拟材料老化效应的示范性过 程的流程图；以及
图7是计算机设备的主要组件的功能框图，本发明的实施例可在该计算机设备中 实施。
具体实施例方式在此将参照图1A-7讨论本发明的实施例。但是本技术领域的人员将会理解，此处 参照附图给出的详细描述用作解释的目的，而本发明延伸到这些有限的实施例之外。首先参照图6，示出了根据本发明的实施例用于在有限元分析中模拟材料老化效 应的示范性方法600的流程图。方法600可以在计算机中结合材料特性测试实验来实施。 方法600最好结合图1A-5B来理解。刚开始，在步骤602，对不同老化程度的感兴趣的时间流变材料执行一组材料特性 测试。可以通过应用单轴向拉伸/压缩(图1A)或者双轴向拉伸(图1B)来执行材料特性 测试。如图1A所示，样本110受到一对力P 112的作用，力P 112形成单轴向拉伸状况，而 图1B中的样本140在两个方向上都在拉伸力F 142的作用下处于双轴向拉伸状况。可以 通过在双轴向拉伸测试设备中使圆形隔膜样本膨胀来得到双轴向拉伸，这是现有技术中的 已知程序。材料特性测试经配置，使得可以对不同老化程度的样本重复进行基本相同的测 试，来得到不同老化程度(例如，天、周、月等)下的时间流变材料特性。一种方法是长期执 行材料特性测试，使得一组样本中每一者都变成不同的老化程度。在步骤604中，每一材料特性测试被配置成通过保持各样本的预定的恒定应力或 者位移(例如，伸长或者延长)，来获得一系列松弛测试数据。例如，使单轴向拉伸测试中 的样本110保持为预定的应力要求最初具有某幅度的一对力P。随着时间的推进，样本110 的材料特性改变，因此要求有不同幅度的力来保持相同的应力。类似地，在双轴向测试中， 初始压力使隔膜样本140膨胀，使得在样本中心两个方向上都达到预定应力。然后压力被 调节，从而在双轴向测试过程中随时间保持相同的应力。松弛测试数据被记录。在一个实施例中，图3所示的力与时间的图绘制了从单轴向测试程序得到的多个 示范性松弛测试数据310a-e。每个松弛测试数据310a-e表示不同老化程度(也就是，0天、 18天、31天、60天和109天)的样本的结果。在另一个实施例中，在图4A的图表400中绘 制了双轴向测试程序的几个示范性松弛测试数据410a_e。图表400是对数-对数标度的 压力_时间图，其中绘制了五个材料特性测试的结果410a_e。每个结果410a_e表示一系 列松弛测试数据。长期执行这一组材料特性测试，以包括不同老化程度的样本的影响。如 图400所示，获得不同老化程度(也就是，1天、1周、1个月、3个月和6个月)下的五个测 试的结果410a_e。回到图6，在步骤606，方法600通过将所获得的测试结果进行移位和匹配，来确定 一组第一和第二与时间相关的材料老化效应参数。例如，在图400中，通过在压力和时间两 个方向上的移位，来匹配每对的五个测试结果410a_e。压力方向由箭头412指示，时间方 向由箭头414指示。在移位和匹配后，在图4B所示的图430中绘制出五个移位后的结果 440a_eo为了进一步解释本发明如何工作，图2A-2B示出了数字实例。针对压缩负载下四 个不同老化程度(也就是，、=0，、(1)，、(2)和ta(3))的单轴向样本计算出所模拟的松 弛测试数据210a_d。绘制出公称应力对比材料特性测试中的松弛时间的曲线。材料的伸展
6率为0.8，对应于样本的20%的压缩应力。材料是可压缩的，且压缩可以数字化表示(n = 10)。长期的材料常数M，Cj和bj设置如下M = 2，Q = 30，C2 = 3，h = 2且b2 = -2。松 弛函数由Prony系列表示如下gj=Cj∑Rke-r该关系式也通过常数&将黏弹特性与弹性特性相联系起来。材料常数以应 力为单位，是衰变常数。黏弹性材料常数是k Rffl Yffl1 0. 50 0. 012 0. 50 0. 103 0. 50 1. 004 0. 50 10. 00模拟的松弛测试数据210a_d在图2A中绘制出。使用模拟的松弛测试数据210a_d 的垂直和水平移位，可以得到如图2B所示的总曲线220。不同老化时间下的一组第一和第 二与时间相关的材料老化效应参数a (ta)和0 (ta)是a (ta)3 (ta)
ta = 01. 001. 00
ta(l)1. 250. 50
t“2)1. 500. 10
ta(3)2. 000. 05位于这些老化时间之间的a (ta)和3 (ta)的值可以从两个相邻值之间推断得到。接下来，在方法600的步骤608中，创建有限元分析应用模块(例如，经编程或者 汇编)，以包含时间流变材料本构方程，该方程使用了在材料特性测试中得到的这一组第一 和第二与时间相关的材料老化效应参数。接下来，包含有时间流变材料本构方程(例如，经 编程)的有限元分析应用模块被用于模拟至少部分包含时间流变材料的结构的材料老化 效应。根据本发明的一个实施例，图5A-B示出了示范性的时间流变材料本构方程。公式 502定义了黏弹高度可压缩材料依据伸长率、的柯西应力(CauchystresshiU)。材料的 应变被定义为(入i_l)，而下标“i”表示空间的三个方向中的每一个。材料常数(”…和！! 由长期测试数据(也就是t =⑴)和未老化样本的松弛测试数据的松弛函数g^t)确定。松弛函数gj(t)可以被看作是一系列指数函数，可以被假定具有公式504的形式。 Rm是无因次的材料常数，是衰减常数。当松弛函数被另一个松弛函数代替 时，公式502变成公式506，因为(^是可以因子分解得到的公共项。对于时间流变材料的材料老化效应，公式508包括第一和第二与时间相关的材 料老化效应参数a (ta)和3 (ta)，ta是老化时间。导出的公式512、514和516示出了 公式508与从材料特性测试得到两个参数的过程中的垂直移位log(a (ta))和水平移位 log(0 (ta))的相互关系。黏弹性材料的本构方程经常可以由第一类的Volterra积分方程式描述。通常此 处给出的方法可以扩展到线性、准线性、Hill-foam黏弹性本构方程，并可扩展到第一类的
7Vol terra积分方程式。根据一方面，本发明涉及一个或多个能够执行在此描述的功能的计算机系统。计 算机系统700的例子在图7中示出。计算机系统700包括一个或多个处理器，例如处理器 704。处理器704连接到计算机系统内部通信总线702。关于该示范性的计算机系统，有各 种软件实现的描述。在读完这一描述后，相关技术领域的人员将会明白如何使用其它计算 机系统和/或计算机架构来实施本发明。计算机系统700还包括主存储器708，优选随机存取存储器(RAM)，还可包括辅助 存储器710。辅助存储器710包括例如一个或多个硬盘驱动器712和/或一个或多个可移 除存储驱动器714，它们代表软磁盘机、磁带驱动器、光盘驱动器等。可移除的存储驱动器 714用已知的方式从可移除存储单元718中读取和/或向可移除存储单元718中写入。可 移除存储单元718代表可以由可移除存储驱动器714读取和写入的软盘、磁带、光盘等。可 以理解，可移除存储单元718包括其上存储有计算机软件和/或数据的计算机可读媒介。在可选实施例中，辅助存储器710可包括其它类似的机制，允许计算机程序或者 其它指令被装载到计算机系统700。这样的机制包括例如可移动存储单元722和接口 720。 这样的例子可包括程序盒式存储器和盒式存储器接口(例如，视频游戏设备中的那些)、可 移动存储芯片(例如可擦除的可编程只读存储器(EPR0M))、通用串行总线(USB)闪存、或 者PR0M)以及相关的插槽、以及其它可移动存储单元722和允许软件和数据从可移动存储 单元722传递到计算机系统700的接口 720。通常，计算机系统700由操作系统(OS)软件 控制和管理，操作系统执行例如进程调度、存储器管理、网络连接和I/O服务。可能还设有连接到总线702的通信接口 724。通信接口 724允许软件和数据在计 算机系统700和外部设备之间传递。通信接口 724的例子包括调制解调器、网络接口(例 如以太网卡)、通信端口、个人计算机存储卡国际协会(PCMCIA)插槽和卡等等。通过通信接 口 724传输的软件和数据是信号728的形式，该信号可以是电子、电磁、光学或者其他可以 被通信接口 724接收的信号。计算机700基于一组特定的规则(也就是，协议)通过数据 网络与其它计算设备通信。通用协议的其中一种是在互联网中通用的TCP/IP (传输控制协 议/互联网协议)。通常，通信接口 724将数据文件组合处理成较小的数据包以通过数据网 络传输，或将接收到的数据包重新组合成原始的数据文件。此外，通信接口 724处理每个数 据包的地址部分以使其到达正确的目的地，或者中途截取发往计算机700的数据包。在这 份文件中，术语“计算机程序媒介”和“计算机可用媒介”都用来指代媒介，例如可移动存储 驱动器714和/或设置在硬盘驱动器712中的硬盘。这些计算机程序产品是用于将软件提 供给计算机系统700的手段。本发明涉及这样的计算机程序产品。计算机系统700还包括输入/输出(I/O)接口 730，它使得计算机系统700能够接 入显示器、键盘、鼠标、打印机、扫描器、绘图机、以及类似设备。计算机程序(也被称为计算机控制逻辑)作为应用模块706存储在主存储器708 和/或辅助存储器710中。也可通过通信接口 724接收计算机程序。这样的计算机程序被 执行时，使得计算机系统700执行如在此所讨论的本发明的特征。特别地，当执行该计算机 程序时，使得处理器704执行本发明的特征。因此，这样的计算机程序代表计算机系统700 的控制器。在本发明采用软件实现的实施例中，该软件可存储在计算机程序产品中，并可使用可移动存储驱动器714、硬盘驱动器712、或者通信接口 724加载到计算机系统700中。应 用模块706被处理器704执行时，使得处理器704执行如在此所述的本发明的功能。主存储器708可被加载一个或多个应用模块706 (例如，其内编程或配置有时间流 变材料本构方程的有限元分析应用模块)，所述应用模块706可被一个或多个处理器704执 行以实现期望的任务，所述处理器可具有或不具有通过I/O接口 730输入的用户输入。在 运行中，当至少一个处理器704执行一个应用模块706时，结果被计算并存储在辅助存储器 710 (也就是，硬盘驱动器712)中。有限元分析(例如，在FEA模型的撞击假人中包括材料 老化效应的汽车抗撞击模拟)的状态以文字或者图形表示的方式通过I/O接口报告给用 户。虽然参照特定的实施例对本发明进行了描述，但是这些实施例仅仅是解释性的， 并不用于限制本发明。本技术领域的人员可得到暗示，对具体公开的示范性实施例做出各 种修改和改变。例如，虽然已经将时间流变材料描述为人造橡胶，但是也可以包括其它材 料，例如橡胶、泡沫等。此外，已经将材料特性测试描述为样本的单轴向拉伸测试，但是也可 以使用压缩测试来代替。最后，发明人想要表明的是，时间流变材料这一术语已由位于美 国加利福利亚州马丁内斯的凯泽医疗公司(Kaiser Permanente Medical Group)的Lydia T.Lee博士和James P. Berry博士提出。总之，本发明的范围不限于在此公开的特定示范性 实施例，对本技术领域人员来说暗含的所有修改都将被包括在本申请的精神和范围以及后 附权利要求的范围内。
9
权利要求
一种用于在计算机辅助工程分析中模拟材料老化效应的方法，其特征在于，所述方法包括为时间流变材料执行多次材料特性测试，每次测试被配置以获得不同老化程度下时间流变材料的一系列松弛测试数据，其中每一系列的松弛测试数据通过将时间流变材料的样本保持在预定的位移或者应变来测量；通过将所述多次材料特性测试中每一对的一系列松弛测试数据相关联，确定一组第一和第二与时间相关的材料老化效应参数；以及配置计算机辅助工程分析模块，使得可以采用时间流变材料本构方程来数学表示材料老化效应，其中所述时间流变材料本构方程中包含所述一组第一和第二与时间相关的材料老化效应参数，所述计算机辅助工程分析应用模块用于执行至少部分包含所述时间流变材料的工程结构的计算机辅助工程分析，其中所述计算机辅助工程分析用于辅助用户在改进所述工程结构的设计时作出决策。
2.根据权利要求1所述的方法，其特征在于，所述每次材料特性测试都包括单轴向测 试程序，松弛测试数据是样本的轴向应力。
3.根据权利要求1所述的方法，其特征在于，所述每次材料特性测试都包括双轴向测 试程序，松弛测试数据是样本的膨胀压力。
4.根据权利要求1所述的方法，其特征在于，所述时间流变材料包括人造橡胶。
5.根据权利要求1所述的方法，其特征在于，将所述多次材料特性测试的每一对的一 系列松弛测试数据相关联还包括将所述多次材料特性测试的每一对的一系列松弛测试数据绘制成具有纵轴和横轴的 二维图表中的第一和第二曲线，其中纵轴表示松弛测试数据，横轴表示对数-对数标度的 时间；以及通过在图表中垂直和水平移位，将第一和第二曲线彼此匹配。
6.根据权利要求5所述的方法，其特征在于，所述第一与时间相关的材料老化效应参 数是匹配所述每一对的垂直移位量。
7.根据权利要求5所述的方法，其特征在于，所述第二与时间相关的材料老化效应参 数是匹配所述每一对的水平移位量。
8.根据权利要求1所述的方法，其特征在于，所述时间流变材料本构方程定义了所述 时间流变材料的应变_应力关系式。
9.一种使用在计算机上运行的计算机辅助工程分析应用模块模拟至少部分包含时间 流变材料的工程结构的材料老化效应的系统，其特征在于，所述系统包括一计算机，被配置成通过从用户输入接收一组第一和第二与时间相关的材料老化效应 参数来执行工程结构的计算机辅助工程分析，其中所述一组第一和第二与时间相关的材料 老化效应参数与计算机辅助工程分析应用模块中配置的时间流变材料本构方程相关；以及 所述计算机根据所述时间流变材料本构方程执行所述包含材料老化效应的工程结构的计 算机辅助工程分析。
10.根据权利要求9所述的系统，其特征在于，所述一组第一和第二与时间相关的材料 老化效应参数通过如下方式确定为时间流变材料执行多次材料特性测试，每次测试被配 置以获得不同老化程度下的时间流变材料的一系列松弛测试数据，其中每一系列的松弛测试数据是通过将时间流变材料的样本保持在预定的位移或者应变来测得的；以及将所述多 次材料特性测试的每一对的一系列松弛测试数据相关联。
11.根据权利要求10所述的系统，其特征在于，将所述多次材料特性测试的每一对的 一系列松弛测试数据相关联还包括将所述多次材料特性测试的每一对的一系列松弛测试数据绘制成具有纵轴和横轴的 二维图表中的第一和第二曲线，其中纵轴表示松弛测试数据，横轴表示对数-对数标度的 时间；以及通过在图表中垂直和水平移位，将第一和第二曲线彼此匹配。
12.一种用于模拟至少部分包含时间流变材料的工程结构内的材料老化效应的计算机 程序产品，所述计算机程序产品嵌入在计算机可读存储媒介中且包括计算机指令，所述计 算机指令用于从用户输入接收一组第一和第二与时间相关的材料老化效应参数，其中所 述第一和第二与时间相关的材料老化效应参数与计算机辅助工程分析应用模块中配置的 时间流变材料本构方程相关；以及根据所述时间流变材料本构方程执行所述包含材料老化 效应的工程结构的计算机辅助工程分析。
13.根据权利要求12所述的计算机程序产品，其特征在于，所述一组第一和第二与时 间相关的材料老化效应参数通过如下方式确定为时间流变材料执行多次材料特性测试， 每次测试被配置以获得不同老化程度下的时间流变材料的一系列松弛测试数据，其中每一 系列的松弛测试数据是通过将时间流变材料的样本保持在预定的位移或者应变来测得的； 以及将所述多次材料特性测试的每一对的一系列松弛测试数据相关联。
14.根据权利要求13所述的计算机程序产品，其特征在于，将所述多次材料特性测试 的每一对的一系列松弛测试数据相关联还包括将所述多次材料特性测试的每一对的一系列松弛测试数据绘制成具有纵轴和横轴的 二维图表中的第一和第二曲线，其中纵轴表示松弛测试数据，横轴表示对数-对数标度的 时间；以及通过在图表中垂直和水平移位，将第一和第二曲线彼此匹配。
全文摘要
本发明公开了一种用于在计算机辅助工程(CAE)分析中模拟时间流变材料的材料老化效应的方法和系统。根据一方面，为感兴趣的时间流变材料执行一组材料特性测试。每次测试获得一系列材料特性例如不同老化程度的松弛测试数据。通过将时间流变材料的样本保持在预定的应变，来测量松弛测试数据。通过将每一对测试之间的一系列松弛测试数据进行移位和匹配，来确定一组第一和第二与时间相关的材料老化效应参数。接下来，通过对至少部分包含时间流变材料的工程结构执行CAE分析，该组第一和第二与时间相关的材料老化效应参数、以及其中包含有时间流变材料本构方程的有限元分析应用模块被用于模拟材料老化效应。
文档编号G06F17/50GK101866377SQ20101015506
公开日2010年10月20日 申请日期2010年3月25日 优先权日2009年4月20日
发明者威廉姆·W·冯, 约翰·O·霍尔奎斯特 申请人:利弗莫尔软件技术公司