CN102308195B - 橡胶材料变形行为预测装置和橡胶材料变形行为预测方法 - Google Patents

橡胶材料变形行为预测装置和橡胶材料变形行为预测方法 Download PDF

Info

Publication number
CN102308195B
CN102308195B CN2010800065557A CN201080006555A CN102308195B CN 102308195 B CN102308195 B CN 102308195B CN 2010800065557 A CN2010800065557 A CN 2010800065557A CN 201080006555 A CN201080006555 A CN 201080006555A CN 102308195 B CN102308195 B CN 102308195B
Authority
CN
China
Prior art keywords
elastomeric material
dimensional model
deformational behavior
rubber
filler
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN2010800065557A
Other languages
English (en)
Other versions
CN102308195A (zh
Inventor
浜谷悟司
芥川惠造
岛广志
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Bridgestone Corp
Original Assignee
Bridgestone Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Bridgestone Corp filed Critical Bridgestone Corp
Publication of CN102308195A publication Critical patent/CN102308195A/zh
Application granted granted Critical
Publication of CN102308195B publication Critical patent/CN102308195B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/32—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/0202—Control of the test
    • G01N2203/0212—Theories, calculations
    • G01N2203/0216—Finite elements
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Immunology (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Pathology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Biochemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

所公开的橡胶材料变形行为预测方法即使从微观层面也可以以良好的精度对橡胶材料变形行为进行分析。更具体地，该橡胶材料变形行为预测方法创建橡胶内已混入有填料的橡胶材料的三维模型，对构成所述三维模型的橡胶层部分提供已基于根据分子动力学得到的厚度信息和温度信息确定了应变和应力之间关系的结构条件，并且分析橡胶材料变形行为。在所述橡胶材料变形行为预测方法中，优选地，对于应用了所述结构条件的三维模型使用有限元方法来分析橡胶材料变形行为。

Description

橡胶材料变形行为预测装置和橡胶材料变形行为预测方法
技术领域
本发明涉及用于预测通过向橡胶添加诸如炭黑和二氧化硅等的填料所形成的橡胶材料的变形行为的装置和方法，尤其涉及即使从微观层面也能够精确地分析橡胶材料的变形行为的变形行为预测方法。
背景技术
已知可以通过向橡胶添加诸如炭黑和二氧化硅等的填料来获得补强效果，并且通过向橡胶添加填料所形成的橡胶材料适用于诸如车辆用轮胎等的橡胶制品。通过实验来测量这种橡胶材料在被施力时的变形行为，并且对测量结果进行评价以设计该橡胶材料中所包含的填料量。
近来，随着诸如有限元方法(FEM)等的数值分析方法和与计算器有关的环境的发展，已提出了通过创建橡胶材料的填料部分和橡胶部分的三维模型来分析变形行为的各种方法。此外，作为能够在缩短分析所需的时间的情况下精确地分析橡胶材料的变形行为的方法，有包括以下步骤的一种方法：利用透射电子显微镜(TEM)对实际的橡胶材料中填料的配置摄像；通过计算机断层成像法(CT法)将由此获得的数据重建为三维基本模型；并且通过有限元方法(FEM)来预测橡胶材料的变形行为(参见日本特开2006-200937(专利文献1))。
顺便提及，作为在橡胶材料包含橡胶和填料的情况下通过FEM进行计算的最简单模型，已知有将通过对橡胶和填料单独进行测量所获得的机械特性应用于三维模型的方法。
然而，认为存在于填料周围的聚合物吸附有填料并且具有与从聚合物本身所获得的机械特性不同的机械特性，并且随着近来测量技术的进步，这种现象已被确认。例如，通过利用原子力显微镜(AFM)测量弹性模量，可以得知存在于填料周围的聚合物与该聚合物本身相比较具有较大的弹性模量。
现有技术文献
专利文献
专利文献1：日本特开2006-200937
发明内容
发明要解决的问题
因此，为了提高对橡胶材料的变形行为的分析准确度，需要对存在于填料周围的橡胶部分、即吸附有填料的橡胶层部分应用与从单体橡胶获得的材料常数不同的常数。
然而，在诸如橡胶等的展现非线性的材料中，难以对吸附有填料的橡胶层部分应用适当的材料常数。此外，尽管可以通过使用AFM来测量吸附有填料的橡胶层部分的弹性模量和粘弹性，但随着橡胶材料的变形量越大越难以测量这两者。
此外，尽管可以使用通过分子动力学方法来对填料和橡胶建模并且使其变形以获得机械特性的方法，但利用分子动力学方法的可测量的时间量程为几纳秒。根据与高分子材料的粘弹性有关的温度-时间转换法则，这与极低温度范围内的粘弹性相对应，因而难以准确地获得机械特性。
如上所述，在准确地分析包括吸附有填料的橡胶层部分的橡胶材料的变形行为方面，仍然存在问题。
考虑到上述事实，本发明的目的在于解决上述传统技术的问题，并且提供即使从微观层面也能够准确地分析橡胶材料的变形行为的橡胶材料变形行为预测方法。此外，本发明的另一目的在于提供可以应用上述方法的橡胶材料变形行为预测装置。
用于解决问题的方案
本发明人为了实现上述目的而进行了深刻的研究，结果发现了通过使用以下的橡胶材料变形行为预测方法，即使从微观层面也可以准确地分析橡胶材料的变形行为，在该橡胶材料变形行为预测方法中：使用各自用于构成橡胶材料的橡胶和填料的二值转换图像来生成三维模型；基于分子动力学方法将该三维模型转换成包括橡胶部分、填料部分和吸附有填料的橡胶层部分的三层形式；并且对该橡胶层部分应用用于根据基于分子动力学方法所获得的厚度信息和温度信息来指定应力和应变之间的关系的结构条件。结果，本发明人做出了本发明。
更具体地，根据本发明的一种橡胶材料变形行为预测方法包括：生成通过向橡胶添加填料所形成的橡胶材料的三维模型；对构成所述三维模型的橡胶层部分应用用于根据基于分子动力学方法所获得的厚度信息和温度信息来指定应力和应变之间的关系的结构条件；以及分析所述橡胶材料的变形行为。
根据基于本发明的橡胶材料变形行为预测方法的优选示例，通过对应用了所述结构条件的三维模型应用有限元方法来分析所述橡胶材料的变形行为。
此外，根据本发明的一种橡胶材料变形行为预测装置包括：用于获得各自示出通过向橡胶添加填料所形成的橡胶材料的截面形状的多个切片图像的部件；用于将各个所述切片图像转换成二值转换图像以区别所述橡胶材料中所包含的橡胶部分和填料部分的部件；用于在所述二值转换图像中设置吸附有填料的橡胶层部分并且将所述二值转换图像转换成三层图像的部件；用于层叠所述三层图像以生成所述橡胶材料的三维模型的部件；用于对所述三维模型应用基于分子动力学方法所指定的结构条件的部件；用于使用应用了所述结构条件的三维模型来分析所述橡胶材料的变形行为的部件；以及用于呈现对所述橡胶材料的变形行为的分析结果的部件。
根据基于本发明的橡胶材料变形行为预测装置的优选示例，所述基于分子动力学方法所指定的结构条件是用于根据三层值来指定所述橡胶部分、所述填料部分以及所述橡胶层部分中的应力和应变之间的关系的结构条件。
根据基于本发明的橡胶材料变形行为预测装置的另一优选示例，所述基于分子动力学方法所指定的结构条件是用于根据基于分子动力学方法所获得的所述橡胶层部分的厚度信息和温度信息来指定应力和应变之间的关系的结构条件。
根据基于本发明的橡胶材料变形行为预测装置的另一优选示例，用于分析所述橡胶材料的变形行为的部件通过对应用了所述结构条件的三维模型应用有限元方法来进行分析。
根据基于本发明的橡胶材料变形行为预测装置的另一优选示例，用于呈现对所述橡胶材料的变形行为的分析结果的部件基于所述分析结果计算应变分布或应力分布，区别应变分布区域或应力分布区域，并且指定各个区域的位置。
发明的效果
根据本发明，可以通过使用以下的橡胶材料变形行为预测方法来从微观层面准确地分析橡胶材料的变形行为，在该橡胶材料变形行为预测方法中：使用各自用于构成橡胶材料的橡胶和填料的二值转换图像来生成三维模型；基于分子动力学方法将该三维模型转换成包括橡胶部分、填料部分和吸附有填料的橡胶层部分的三层形式；并且对该橡胶层部分应用用于根据基于分子动力学方法所获得的厚度信息和温度信息来指定应力和应变之间的关系的结构条件。利用该结构，可以预测由于橡胶材料的应变所引起的材料特性的变化，并且控制诸如橡胶材料的断裂强度和损耗等的材料特性。
附图说明
图1是示出表示根据本发明的变形行为预测装置的一个示例的橡胶材料变形行为预测系统的结构的图。
图2是用于解释构成图1所示的橡胶材料变形行为预测系统的计算机的电气系统的结构的主要部分的图。
图3是示出通过将切片图像转换成二值形式所获得的二值转换图像的图。
图4是示出通过将二值转换图像转换成三值形式所获得的三层图像的图。
图5是示出由基于分子动力学方法所获得的由填料和聚合物所形成的体系的示意图。
图6是示出相对于填料表面的距离和扩散系数之间的关系的图。
图7是示出橡胶的扩散系数和温度之间的关系的图。
图8是示出对橡胶材料的变形行为的分析结果的截面图像图。
图9是示出通过FEM计算所获得的橡胶材料的应力-应变曲线的图。
图10是通过将图9所示的关系转换成杨氏(Young’s)模量和应变之间的关系所获得的图。
具体实施方式
以下将参考附图来详细说明根据本发明的变形行为预测方法和变形行为预测装置。图1是示出用作根据本发明的变形行为预测装置的一个示例的橡胶材料变形行为预测系统的结构的图。图2是用于解释构成图1所示的橡胶材料变形行为预测系统的计算机的电气系统的结构的主要部分的图。
图1所示的变形行为预测装置1包括CT扫描器(计算机断层成像扫描器)2和计算机3。CT扫描器2和计算机3经由线缆4相互连接。
CT扫描器2具有透射电子显微镜(TEM)和试样台。CT扫描器2通过透射电子显微镜对放置在试样台上的要分析的橡胶材料摄像，并且通过计算机断层成像法(CT法)将通过摄像所获得的数据重建为三维基本模型。此外，CT扫描器2相对于预定平面按预定间隔对重建得到的三维基本模型进行切片，以生成多个切片图像数据。注意，在本发明中，尽管没有示出，但可以通过利用聚焦离子束(FIB)对橡胶材料的表面进行处理(例如，蚀刻)、并且利用扫描电子显微镜(SEM)观察处理后的表面来获得多个切片图像。在这种情况下，代替CT扫描器2，本发明具有能够发出聚焦离子束的扫描电子显微镜系统。
计算机3具有：键盘5，用于输入与分析有关的各种条件；计算机本体6，用于根据预先存储的各处理程序来分析橡胶材料的变形行为；和显示器7，用于显示计算机本体6的计算结果等。计算机3使用CT扫描器2所生成的切片图像数据来分析橡胶材料的变形行为等。计算机本体6包括软盘驱动单元(以下称为FDU)9，其中用作记录介质的软盘(以下称为FD)8可以插入FDU 9中或者从FDU 9移除。
如图2所示，计算机3包括：中央处理单元(CPU)10，用于管理整个装置的操作；ROM 11，用于预先存储包括用于控制计算机3的控制程序的各种程序和各种参数等；RAM 12，用于暂时存储各种数据；外部I/O控制单元14，其连接至与线缆4相连接的连接器13，用于经由连接器13从CT扫描器2获得切片图像数据；硬盘驱动器(HDD)15，用于存储所获得的切片图像数据；软盘I/F单元16，用于与FDU 9所配置的FD 8进行数据输入/输出；显示器驱动器17，用于控制在显示器7上显示各种信息；和操作输入检测单元18，用于检测对键盘5的键操作。
CPU 10、RAM 12、ROM 11、HDD 15、外部I/O控制单元14、软盘I/F单元16、显示器驱动器17和操作输入检测单元18经由系统总线BUS相互连接。利用该结构，CPU 10可以对RAM 12、ROM 11和HDD 15进行存取；经由软盘I/F单元16对FDU 9所设置的FD 8进行存取；控制经由外部I/O控制单元14进行的数据发送/接收；并且经由显示器驱动器17将各种信息显示在显示器7上。此外，CPU 10可以随时掌握对键盘5的键操作。
应当注意，可以使用FDU 9来从FD 8读取处理程序和数据等或者向FD 8写入处理程序和数据等。因此，可以预先将处理程序和数据等记录在FD 8中，并且经由FDU 9执行记录在FD 8中的处理程序。此外，可以将记录在FD 8中的处理程序安装到HDD 15并且执行所安装的处理程序。此外，作为记录介质，可以采用磁带、诸如CD-ROM和DVD等的光盘以及诸如MD和MO等的磁光盘。在采用这些记录介质的情况下，代替FDU 9而使用读/写单元，或者使用与FDU 9兼容的读/写单元。
此外，图1所示的橡胶材料变形行为预测系统的结构仅是例子，并且根据应用，可以按照需要修改已知的结构。
以下将详细说明根据本发明的橡胶材料变形行为预测方法。根据本发明的橡胶材料变形行为预测方法提供了以下的橡胶材料变形行为预测方法：获得各自示出通过向橡胶添加填料所形成的橡胶材料的截面形状的多个切片图像；将各个切片图像转换成二值转换图像以区别该橡胶材料中所包含的橡胶部分和填料部分，由此生成三维模型；并且，在该变形行为预测方法中，基于分子动力学方法将该三维模型形成为具有橡胶部分、填料部分和吸附有填料的橡胶层部分的三层形式，并且对该橡胶层部分应用用于根据基于分子动力学方法所获得的厚度信息和温度信息来指定应力和应变之间的关系的结构条件。通过对橡胶层部分应用根据分子动力学方法所获得的特定结构条件，可以在预测橡胶材料的变形行为时考虑到该橡胶层部分，由此可以准确地分析该橡胶材料的变形行为。
在根据本发明的变形行为预测方法中，首先获得各自示出通过向橡胶添加填料所形成的橡胶材料的截面形状的多个切片图像，以了解该橡胶材料的内部构造。在图1所示的橡胶材料变形行为预测系统中，用户对要分析的橡胶材料应用胶体金标记，并且将该橡胶材料放置在CT扫描器2所配置的试样台上。然后，对CT扫描器2进行用于开始处理的操作以拍摄切片图像。注意，CT扫描器2由具有使用透射电子显微镜法(TEMT)的计算机结构的测量装置构成。CT扫描器2在使透射电子显微镜和放置有橡胶材料的试样台相对于彼此以在预定角度范围内(例如，-60度到+60度的范围内)的各预定角度(例如，每2度)转动移动的情况下进行扫描，以拍摄橡胶材料的连续倾斜的图像。CT扫描器2通过使用拍摄到的多个角度图像(例如，61个图像)的图像数据来获得这些图像之间的转动轴，并且通过计算机断层成像法将这些图像重建为三维基本模型。然后，CT扫描器2生成切片图像，其中这些切片图像是通过在与重建得到的三维基本模型的各个平面平行的平面处按预定间隔(例如，4nm间隔)对该三维基本模型进行切片所获得的。
应当注意，橡胶材料是通过向橡胶添加填料所获得的，并且根据应用，可以从橡胶工业领域通常使用的橡胶和填料中选择上述的橡胶和填料。此外，没有特别限制橡胶材料的形状，但优选地，使橡胶材料形成为利用CT扫描器2可以容易地获得切片图像的形状，并且更具体地，使橡胶材料形成为诸如立方体和平行六面体等的六面体形状。
接着，在根据本发明的变形行为预测方法中，将通过以上处理所获得的各自示出橡胶材料的截面形状的多个切片图像转换成二值图像，以判断该橡胶材料中所包含的橡胶部分和填料部分。在上述切片图像中，由于在各自用于构成橡胶材料的橡胶部分和填料部分之间透射率实质上不同，通常为填料部分深(浓度值大)并且橡胶部分浅(浓度值小)。因此，可以基于切片图像中各像素的密集度来判断切片图像中的橡胶部分和填料部分。可以通过实验等预先指定用于判断切片图像中的橡胶部分和填料部分的浓度值。
在图1所示的橡胶材料变形行为预测系统中，将通过实验等预先指定的用于判断切片图像中的橡胶部分和填料部分的浓度值设置为阈值h；将切片图像中的各像素的浓度值分别与阈值h进行比较；将各个像素转换成二值数字以获得二值转换图像，由此生成二值转换图像数据。图3是示出通过将切片图像转换成二值形式所获得的二值转换图像的图。
应当注意，为了进一步精确地提取橡胶材料内的填料部分，用于转换成二值图像的处理包括：将切片图像中的各个像素的浓度值与阈值h进行比较；并且使位于预定数量(例如，5个以上)的浓度值大于或等于阈值h的像素沿着垂直方向和水平方向连续配置的部分处的像素变黑并且使其它的像素变白以获得二值转换图像，由此生成二值转换图像数据。此外，可以将有色的二值转换图像数据的格式转换成变黑部分中的各个像素的值为“1”并且变白部分中的各个像素的值为“0”的二值转换图像数据。
在根据本发明的变形行为预测方法中，通过指定吸附有填料的橡胶层部分来将二值转换图像转换成三层图像。如后面所述，通过转换成三层图像，可以基于分子动力学方法来将橡胶材料的三维模型转换成三层模型。在图1所示的橡胶材料变形行为预测系统中，通过进行将值为“0”且相邻像素的值为“1”的像素的值变为“2”的这种转换，将格式已转换为如上所述的二值转换图像数据进一步格式转换成三层图像。更具体地，吸附有填料的橡胶层部分的像素的值为“2”。图4是示出通过将二值转换图像转换成三值转换图像所获得的三层图像，其中填料部分、橡胶部分和吸附有填料的橡胶层部分分别为黑色、白色和灰色。应当注意，如上所述的用于确定橡胶层部分的方法仅是一个例子，并且可以根据应用基于像素的大小来确定橡胶层部分。例如，通过代替将橡胶层部分限制为与填料部分的像素邻接的橡胶部分的像素而将其它的橡胶部分设置为橡胶层部分以使得橡胶层部分具有特定厚度，或者通过将橡胶层部分转换成更多个部分，可以进一步精确地分析变形行为。
接着，在根据本发明的变形行为预测方法中，层叠三层图像以生成橡胶材料的三维模型。在图1所示的橡胶材料变形行为预测系统中，将格式转换得到的三层图像以符合诸如获得了相应切片图像的位置等的条件的方式层叠，并且形成了橡胶材料的三维结构，由此可以生成将三层图像中的各像素设置为格子单元的三维模型。更具体地，在如上所述的三维模型的格式转换得到的三层图像中，像素值为“0”的部分表示橡胶部分；像素值为“1”的部分表示填料部分；并且像素值为“2”的部分表示吸附有填料的橡胶层部分。注意，对于所生成的橡胶材料的三维模型，通过在三层图像之间进行能够通过将具有同一值的像素放置在同一格子区域内来合并这些三层图像的图像处理，可以生成橡胶材料的计算用立体图像以产生显示器7上的立体图像。
在根据本发明的变形行为预测方法中，由于基于分子动力学方法将上述的三维模型形成为包括橡胶部分、填料部分和吸附有填料的橡胶层部分的三层，因此需要对通过层叠三层图像所生成的三维模型应用基于分子动力学方法所指定的结构条件。更具体地，基于分子动力学方法来指定结构条件，以使得对吸附有填料的橡胶层部分应用与应用于橡胶部分的结构条件不同的结构条件。
以下将说明用于确定应用于橡胶材料的吸附有填料的橡胶层部分的结构条件的方法的一个例子。首先，基于分子动力学方法来构造由填料表面和橡胶(聚合物)所形成的体系。图5是示出基于分子动力学方法所获得的由填料和聚合物所形成的体系的示意图。然后，计算该体系内橡胶的均方位移，以基于该计算结果获得扩散系数。图6是示出相对于填料表面的距离和扩散系数之间的关系的图，其中，y轴表示扩散系数并且x轴表示相对于填料表面的距离。注意：x轴与图5相对应；x轴的两端表示各个填料表面的位置；并且x轴的标度分别表示相对于一侧的填料表面的距离。根据图6可以理解以下：吸附有填料的橡胶层部分由于填料的键合而与橡胶部分相比具有较小的扩散系数。此外，与该计算平行地，计算聚合物本身的均方位移相对于温度的变化率，以获得该聚合物的扩散系数的温度依赖性。图7是示出橡胶的扩散系数和温度之间的关系的图，其中，y轴表示扩散系数并且x轴表示温度。根据图7可以理解以下：橡胶的扩散系数随着温度上升而增大，并且特别地，从玻璃化转变点附近(-52℃附近)开始急剧增大。更具体地，这表示，在扩散系数较小的橡胶层部分中，橡胶处于低温状态。该事实证实了存在于填料附近的聚合物相为玻璃态的一般观念。如通过这些结果可以明确理解，通过将结构条件从传统上应用于橡胶部分的温度不变的结构条件变为温度变化的结构条件，使得可以输入与相对于填料的距离相对应的材料特性变化，由此可以从微观层面精确地分析实际的橡胶材料的变形行为。
应当注意，优选对三维模型的橡胶层部分应用用于指定应力和应变之间的关系的结构条件，并且如上所述，还优选应用用于根据基于分子动力学方法所获得的厚度信息和温度信息来指定应力和应变之间的关系的结构条件。更具体地，对与三维模型的橡胶层部分相对应的格子区域应用在与橡胶层部分的扩散系数相对应的温度下实际测量出的应力和应变之间的关系，作为结构条件。注意，用于指定应力和应变之间的关系的结构条件包括如日本特开2006-200937所述的广义门尼-里夫林(Mooney-Rivlin)方程式和广义奥格登(Ogden)方程式以及日本特开2005-345413所公开的以下方程式(1)：
G = P · exp ( β · k ) · cosh ( 2 β · I 1 ) - 1 exp ( β · k ) · cosh ( 2 β · I 1 ) + 1 + Q · T ΔT · ∂ S ∂ I 1 · · · ( 1 )
[其中，G表示杨氏模量；S表示橡胶变形时的熵变化；P和Q各自表示与弹性模量有关的系数；I1表示应变的不变量；T表示绝对温度；β等于1/(kΔT)；k表示玻尔兹曼(Boltzmann)常数；并且ΔT表示与橡胶的玻璃化转变温度之差]。
在预先获得了用于指定填料的应力和应变之间的关系的结构条件作为三维模型的填料部分的结构条件的情况下，优选对三维模型的填料部分的格子区域应用这些结构条件。可选地，可以使用通过利用实验预先测量填料的硬度所获得的实际值作为结构条件，或者使用根据填料的结晶部和非晶部之间的比所获得的估计值作为结构条件。此外，通常，橡胶材料中所包含的填料比橡胶硬并且具有与橡胶的杨氏模量(弹性模量)相比较大的杨氏模量，由此作为针对三维模型的填料部分的结构条件，可以应用比根据应用于橡胶部分的结构条件所获得的杨氏模量大了预定倍数(例如，1000倍)的杨氏模量。
此外，优选对三维模型的橡胶部分应用用于指定应力和应变之间的关系的结构条件，并且可以优选使用如日本特开2006-200937所述的广义Mooney-Rivlin方程式和广义Ogden方程式以及上述方程式(1)等。
接着，在根据本发明的变形行为预测方法中，使用应用了上述结构条件的三维模型来分析橡胶材料的变形行为。在图1所示的橡胶材料变形行为预测系统中，用户经由键盘5输入要分析的三维模型和分析用的条件，并且将这两者提供至计算机3，然后进行用于分析橡胶材料的变形行为的处理。注意，作为分析用的条件，可以指定三维模型变化的方向以及三维模型在该指定方向上进行了拉伸、压缩或剪切变形时的变化率。此外，通过分析处理，可以分析当三维模型在被指定为分析用的条件的方向上进行了拉伸、压缩或剪切变形时的三维模型的应变和内部应力分布以及整体三维模型的应力值。由于应用了上述结构条件的三维模型具有与实际橡胶材料的构造相同的构造，因此可以通过使用有限元方法(FEM)分析橡胶材料的变形行为来精确地分析橡胶材料内的弹性模量和应力分布。
以下将详细说明根据本发明的橡胶材料变形行为预测装置。根据本发明的橡胶材料变形行为预测装置涉及用于进行上述的变形行为预测方法的装置，并且包括图1所示的橡胶材料变形行为预测系统。更具体地，根据本发明的变形行为预测装置包括：用于获得各自示出通过向橡胶添加填料所形成的橡胶材料的截面形状的多个切片图像的部件；用于将各个切片图像转换成二值转换图像以区别橡胶材料中所包含的橡胶部分和填料部分的部件；用于在二值转换图像中设置吸附有填料的橡胶层部分以将二值转换图像转换成三层图像的部件；用于层叠三层图像以生成橡胶材料的三维模型的部件；用于对三维模型应用基于分子动力学方法所指定的结构条件的部件；用于使用应用了结构条件的三维模型来分析橡胶材料的变形行为的部件；以及用于呈现对橡胶材料的变形行为的分析结果的部件。
在根据本发明的变形行为预测装置中用于对三维模型应用基于分子动力学方法所指定的结构条件的部件中，优选地，基于分子动力学方法所指定的结构条件是用于根据三层值来指定橡胶部分、填料部分以及橡胶层部分中的应力和应变之间的关系的结构条件，也可以是用于根据基于分子动力学方法所获得的橡胶层部分的厚度信息和温度信息来指定应力和应变之间的关系的结构条件。
在根据本发明的变形行为预测装置中，优选地，用于分析橡胶材料的变形行为的部件通过对应用了上述结构条件的三维模型应用有限元方法来进行分析。此外，优选地，用于呈现对橡胶材料的变形行为的分析结果的部件基于该分析结果计算应变分布或应力分布，区别应变分布区域或应力分布区域，并且指定各个区域的位置。
示例
以下将通过给出实施例来更加详细地说明本发明。然而，本发明不限于以下示例。
在图1所示的橡胶材料变形行为预测系统中，利用三维透射电子显微镜对通过向100质量部分的橡胶添加30质量部分的炭黑所形成的橡胶材料进行摄像，以生成该橡胶材料的三维模型。由于将由此获得的三维模型成形为三层，因此对橡胶部分和填料部分的格子区域应用由上述方程式(1)所表示的表现弹性模量的温度和应变依赖性的本构方程式以及基于实际测量值或上述的估计测量值所获得的杨氏模量(弹性模量)作为结构条件，并且对橡胶层部分的格子区域应用用于根据基于分子动力学方法所获得的厚度信息和温度信息来指定应力和应变之间的关系的结构条件。此外，对应用了结构条件的三维模型进行FEM计算。在图8至图11中示出相应的计算结果。
图8是示出对橡胶材料的变形行为的分析结果的截面图像图。注意，图8示出使用三维模型数据在将整体三维模型沿着z方向拉伸了15％的情况下对应变状态和应力分布的分析结果。对于应力分布，浓度越大的部分表示应力值越高。
图9是示出通过FEM计算所获得的橡胶材料的应力-应变曲线的图，并且图10是通过将图9所示的关系转换成杨氏模量和应变之间的关系所获得的图。根据这些结果可知：三层的三维模型示出与对橡胶部分和填料部分进行了二值转换的传统的三维模型不同的变形行为。
上述结果以极高的准确度从微观层面表现了橡胶材料的变形行为。因此，基于上述分析结果，可以精确地预测实际的橡胶材料的变形行为，由此可以有效地控制橡胶材料的诸如断裂强度和损耗等的材料特性。
附图标记说明
1 橡胶材料变形行为预测系统
3 计算机
4 线缆
5 键盘
6 计算机本体
7 显示器
8 FD
9 FDU
13 连接器

Claims (5)

1.一种橡胶材料变形行为预测方法，包括：
生成通过向橡胶添加填料所形成的橡胶材料的三维模型，所述三维模型形成为包括橡胶部分、填料部分和吸附有填料的橡胶层部分的三层；
对构成所述三维模型的所述橡胶层部分应用用于根据基于分子动力学方法所获得的厚度信息和温度信息来指定应力和应变之间的关系的结构条件；以及
分析所述橡胶材料的变形行为。
2.根据权利要求1所述的橡胶材料变形行为预测方法，其特征在于，
通过对应用了所述结构条件的三维模型应用有限元方法来分析所述橡胶材料的变形行为。
3.一种橡胶材料变形行为预测装置，包括：
用于获得各自示出通过向橡胶添加填料所形成的橡胶材料的截面形状的多个切片图像的部件；
用于将各个所述切片图像转换成二值转换图像以区别所述橡胶材料中所包含的橡胶部分和填料部分的部件；
用于在所述二值转换图像中设置吸附有填料的橡胶层部分并且将所述二值转换图像转换成三层图像的部件；
用于层叠所述三层图像以生成所述橡胶材料的三维模型的部件；
用于对所述三维模型应用基于分子动力学方法所指定的结构条件的部件，其中，所述基于分子动力学方法所指定的结构条件是用于根据基于分子动力学方法所获得的所述橡胶层部分的厚度信息和温度信息来指定应力和应变之间的关系的结构条件；
用于使用应用了所述结构条件的三维模型来分析所述橡胶材料的变形行为的部件；以及
用于呈现对所述橡胶材料的变形行为的分析结果的部件。
4.根据权利要求3所述的橡胶材料变形行为预测装置，其特征在于，
用于分析所述橡胶材料的变形行为的部件通过对应用了所述结构条件的三维模型应用有限元方法来进行分析。
5.根据权利要求3所述的橡胶材料变形行为预测装置，其特征在于，
用于呈现对所述橡胶材料的变形行为的分析结果的部件基于所述分析结果计算应变分布或应力分布，区别应变分布区域或应力分布区域，并且指定各个区域的位置。
CN2010800065557A 2009-02-03 2010-02-02 橡胶材料变形行为预测装置和橡胶材料变形行为预测方法 Active CN102308195B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2009022923A JP4603082B2 (ja) 2009-02-03 2009-02-03 ゴム材料の変形挙動予測装置及びゴム材料の変形挙動予測方法
JP2009-022923 2009-02-03
PCT/JP2010/051733 WO2010090295A1 (ja) 2009-02-03 2010-02-02 ゴム材料の変形挙動予測装置及びゴム材料の変形挙動予測方法

Publications (2)

Publication Number Publication Date
CN102308195A CN102308195A (zh) 2012-01-04
CN102308195B true CN102308195B (zh) 2013-11-20

Family

ID=42542185

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010800065557A Active CN102308195B (zh) 2009-02-03 2010-02-02 橡胶材料变形行为预测装置和橡胶材料变形行为预测方法

Country Status (6)

Country Link
US (1) US9002687B2 (zh)
EP (1) EP2395339A4 (zh)
JP (1) JP4603082B2 (zh)
KR (1) KR101296502B1 (zh)
CN (1) CN102308195B (zh)
WO (1) WO2010090295A1 (zh)

Families Citing this family (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5269732B2 (ja) * 2009-09-28 2013-08-21 株式会社ブリヂストン ゴム材料の変形挙動予測方法およびそれに用いられる装置
JP5559594B2 (ja) * 2010-05-20 2014-07-23 住友ゴム工業株式会社 ゴム材料のシミュレーション方法
EP2562534A3 (en) * 2011-08-23 2014-05-21 Sumitomo Rubber Industries, Ltd. Method for analyzing rubber compound with filler particles
JP2013057638A (ja) * 2011-09-09 2013-03-28 Sumitomo Rubber Ind Ltd ゴム材料のシミュレーション方法
JP5395864B2 (ja) * 2011-09-14 2014-01-22 住友ゴム工業株式会社 ゴム材料のシミュレーション方法
JP5503618B2 (ja) * 2011-10-03 2014-05-28 住友ゴム工業株式会社 ゴム材料のシミュレーション方法
JP2013108800A (ja) * 2011-11-18 2013-06-06 Sumitomo Rubber Ind Ltd ゴム材料のシミュレーション方法
JP5548180B2 (ja) * 2011-11-24 2014-07-16 住友ゴム工業株式会社 高分子材料のシミュレーション方法
JP6072525B2 (ja) * 2012-12-04 2017-02-01 東洋ゴム工業株式会社 フィラー充填ポリマーモデルの作成装置、その方法及びプログラム
JP6082303B2 (ja) * 2013-04-04 2017-02-15 住友ゴム工業株式会社 高分子材料のシミュレーション方法
JP6055359B2 (ja) * 2013-04-15 2016-12-27 住友ゴム工業株式会社 高分子材料のシミュレーション方法
JP6200193B2 (ja) * 2013-04-15 2017-09-20 住友ゴム工業株式会社 高分子材料のシミュレーション方法
KR102262622B1 (ko) * 2013-10-07 2021-06-08 스미토모 고무 코교 카부시키카이샤 필러 배합 고무의 유한 요소 모델의 작성 방법
JP5913260B2 (ja) * 2013-11-14 2016-04-27 住友ゴム工業株式会社 高分子材料のシミュレーション方法
US10520394B2 (en) 2014-09-29 2019-12-31 Lord Corporation Elastomeric degradation indicator devices, systems, and methods
KR20160131278A (ko) * 2015-05-06 2016-11-16 엘지전자 주식회사 이동 단말기
KR102028018B1 (ko) 2017-03-27 2019-10-02 주식회사 엘지화학 고무 보강용 유-무기 복합체의 고무 보강 효과를 예측하는 방법
EP3447717A1 (en) * 2017-08-24 2019-02-27 Tata Consultancy Services Limited Systems and methods for determining properties of composite materials for predicting behaviour of structures
CN107782608B (zh) * 2017-09-08 2020-12-29 吉林大学 一种材料孔挤压实验方法
US20210097211A1 (en) * 2018-04-09 2021-04-01 Covestro Llc Techniques to custom design products
CN108709801B (zh) * 2018-04-16 2020-09-22 四川理工学院 一种检测载荷下高聚物复合材料中填料分散状态的方法
US10823655B2 (en) * 2018-10-10 2020-11-03 Deere & Company Method and apparatus to measure strains in adhesives of joints
CN112595634B (zh) * 2020-12-14 2021-09-24 青岛理工大学 一种三维颗粒材料的内部变形分析实验装置及方法

Family Cites Families (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0919941B1 (en) * 1997-11-25 2005-09-21 Sumitomo Rubber Industries Limited Method of and apparatus for simulating rolling tyre
EP1293917B1 (en) * 2001-09-18 2005-07-27 Sumitomo Rubber Industries Ltd. A method for tire rolling simulation
EP1526468B1 (en) * 2003-10-17 2009-09-30 Sumitomo Rubber Industries Limited Method of simulating viscoelastic material
US7373284B2 (en) * 2004-05-11 2008-05-13 Kimberly-Clark Worldwide, Inc. Method of evaluating the performance of a product using a virtual environment
JP4299735B2 (ja) * 2004-06-07 2009-07-22 株式会社ブリヂストン ゴム製品の弾性応答性能の予測方法及び設計方法
JP4594043B2 (ja) * 2004-11-15 2010-12-08 住友ゴム工業株式会社 ゴム材料のシミュレーション方法
JP4695399B2 (ja) * 2005-01-11 2011-06-08 住友ゴム工業株式会社 フィラー配合材料のシミュレーションモデル生成方法
JP4694847B2 (ja) * 2005-01-18 2011-06-08 株式会社ブリヂストン ゴム材料の形態表示方法
JP4602776B2 (ja) * 2005-01-18 2010-12-22 株式会社ブリヂストン ゴム材料の変形挙動予測方法及びゴム材料の変形挙動予測装置
JP4681417B2 (ja) * 2005-10-12 2011-05-11 住友ゴム工業株式会社 高分子材料のシミュレーション方法
JP2007193560A (ja) * 2006-01-19 2007-08-02 Hitachi Ltd 制御システム向け分散データベースシステム
JP2007272416A (ja) * 2006-03-30 2007-10-18 Bridgestone Corp ゴム製品の弾性応答性能の予測方法、設計方法、及び弾性応答性能予測装置
JP4685747B2 (ja) * 2006-11-09 2011-05-18 住友ゴム工業株式会社 ゴム材料解析モデルの作成方法
JP5012094B2 (ja) * 2007-03-06 2012-08-29 横浜ゴム株式会社 粘弾性特性決定方法、シミュレーションプログラムおよびシミュレーション装置、および、粘弾性特性を備えたシミュレーションモデルを用いるシミュレーション方法
EP2164694B1 (en) * 2007-07-02 2012-05-30 MAGMA Giessereitechnologie GmbH Method for describing the statistical orientation distribution of particles in a simulation of a mould filling process.
JP4253822B2 (ja) * 2007-07-18 2009-04-15 株式会社テクニカルスルー 粒子法の解析プログラム

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
JP特开2006-193560A 2006.07.27
JP特开2006-200938A 2006.08.03
JP特开2007-107968A 2007.04.26
JP特开2007-272416A 2007.10.18
JP特开2008-122154A 2008.05.29
JP特开2008-213724A 2008.09.18

Also Published As

Publication number Publication date
KR20110111458A (ko) 2011-10-11
US20110288838A1 (en) 2011-11-24
CN102308195A (zh) 2012-01-04
JP4603082B2 (ja) 2010-12-22
EP2395339A4 (en) 2015-10-21
JP2010181194A (ja) 2010-08-19
WO2010090295A1 (ja) 2010-08-12
US9002687B2 (en) 2015-04-07
KR101296502B1 (ko) 2013-08-13
EP2395339A1 (en) 2011-12-14

Similar Documents

Publication Publication Date Title
CN102308195B (zh) 橡胶材料变形行为预测装置和橡胶材料变形行为预测方法
Gadelrab et al. Densification modeling of fused silica under nanoindentation
WO2010090294A1 (ja) ゴム材料の変形挙動予測装置及びゴム材料の変形挙動予測方法
JP5269732B2 (ja) ゴム材料の変形挙動予測方法およびそれに用いられる装置
JP4602776B2 (ja) ゴム材料の変形挙動予測方法及びゴム材料の変形挙動予測装置
Yun et al. A self-optimizing inverse analysis method for estimation of cyclic elasto-plasticity model parameters
Pilch et al. Measurement of whole-field surface displacements and strain using a genetic algorithm based intelligent image correlation method
Amirrahmat et al. Micro shear bands: Precursor for strain localization in sheared granular materials
Čech et al. Examination of Berkovich indenter tip bluntness
US20080294397A1 (en) Analytical method, recording medium, and analyzing apparatus
Hamar et al. State-of-health estimation using a neural network trained on vehicle data
Lee et al. System identification of spatial distribution of structural parameters using modified transitional Markov Chain Monte Carlo method
Chen et al. Full-field prediction of stress and fracture patterns in composites using deep learning and self-attention
Krier et al. Introduction of the real tip defect of Berkovich indenter to reproduce with FEM nanoindentation test at shallow penetration depth
Zeng et al. Bayesian model updating for structural dynamic applications combing differential evolution adaptive metropolis and kriging model
US8160848B2 (en) Apparatus for generating coarse-grained simulation image of sample to be measured with a probe of a scanning probe microscope
JP4694847B2 (ja) ゴム材料の形態表示方法
Aswal et al. Strain‐based joint damage estimation approach robust to unknown non‐stationary input force
Landauer et al. A materials data framework and dataset for elastomeric foam impact mitigating materials
Lyngdoh et al. Experimental monitoring and modeling of fatigue damage for 3D-printed polymeric beams under irregular loading
Smyl et al. Stacked elasticity imaging approach for visualizing defects in the presence of background inhomogeneity
JP2009110228A (ja) 固体粒子を分散させた高分子化合物の分子運動解析方法および解析プログラム
Boubaker et al. Non-destructive testing of baked anodes based on modal analysis and principal component analysis
Heindel et al. Fatigue monitoring and maneuver identification for vehicle fleets using a virtual sensing approach
JP2007271369A (ja) ゴム材料の変形挙動予測装置及びゴム材料の変形挙動予測方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant