CN103154703A - 橡胶制品的弹性响应性能的预测方法、设计方法以及弹性响应性能预测装置 - Google Patents
橡胶制品的弹性响应性能的预测方法、设计方法以及弹性响应性能预测装置 Download PDFInfo
- Publication number
- CN103154703A CN103154703A CN2011800483590A CN201180048359A CN103154703A CN 103154703 A CN103154703 A CN 103154703A CN 2011800483590 A CN2011800483590 A CN 2011800483590A CN 201180048359 A CN201180048359 A CN 201180048359A CN 103154703 A CN103154703 A CN 103154703A
- Authority
- CN
- China
- Prior art keywords
- centerdot
- prime
- beta
- response performance
- rubber
- 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.)
- Granted
Links
- 229920001971 elastomer Polymers 0.000 title claims abstract description 108
- 239000005060 rubber Substances 0.000 title claims abstract description 108
- 230000004044 response Effects 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000013461 design Methods 0.000 title claims description 6
- 230000003993 interaction Effects 0.000 claims abstract description 13
- 238000004458 analytical method Methods 0.000 claims abstract description 12
- 230000001419 dependent effect Effects 0.000 claims description 25
- 230000004001 molecular interaction Effects 0.000 claims description 20
- 238000004132 cross linking Methods 0.000 claims description 11
- 230000009477 glass transition Effects 0.000 claims description 6
- 238000010008 shearing Methods 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 2
- 238000007906 compression Methods 0.000 claims description 2
- 238000002425 crystallisation Methods 0.000 description 15
- 230000008025 crystallization Effects 0.000 description 15
- 238000005192 partition Methods 0.000 description 12
- 239000013536 elastomeric material Substances 0.000 description 8
- 244000043261 Hevea brasiliensis Species 0.000 description 5
- 229920003052 natural elastomer Polymers 0.000 description 5
- 229920001194 natural rubber Polymers 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000009878 intermolecular interaction Effects 0.000 description 4
- 230000008863 intramolecular interaction Effects 0.000 description 4
- 229920003048 styrene butadiene rubber Polymers 0.000 description 4
- RRHGJUQNOFWUDK-UHFFFAOYSA-N Isoprene Chemical compound CC(=C)C=C RRHGJUQNOFWUDK-UHFFFAOYSA-N 0.000 description 3
- 230000000052 comparative effect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 229920001195 polyisoprene Polymers 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000000945 filler Substances 0.000 description 2
- 230000008676 import Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 206010013710 Drug interaction Diseases 0.000 description 1
- 244000287680 Garcinia dulcis Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000004913 activation Effects 0.000 description 1
- 230000003542 behavioural effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- 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
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/44—Resins; Plastics; Rubber; Leather
- G01N33/445—Rubber
-
- 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/0058—Kind of property studied
- G01N2203/0092—Visco-elasticity, solidification, curing, cross-linking degree, vulcanisation or strength properties of semi-solid materials
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Biochemistry (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Analytical Chemistry (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Tires In General (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
在使用有限元分析法来预测表示橡胶制品的变形行为的弹性响应性能的弹性响应性能预测方法中,利用使用表示分子间的相互作用的参数来表示橡胶制品的应变能量的温度和应变依赖性的结构方程式,来预测橡胶制品的弹性响应性能。
Description
技术领域
本发明涉及一种橡胶制品的弹性响应性能的预测方法、设计方法以及弹性响应性能预测装置,特别是涉及一种使用有限元分析法(FEA)来预测橡胶制品的弹性响应性能的方法、使用了该预测方法的橡胶制品的设计方法以及弹性响应性能预测装置。
背景技术
以下方法是已经经过了数十年的实例:在设计橡胶制品时,使用三维有限元分析法(FEA)来预测其弹性响应性能,应用其分析或者仿真结果。使用于该FEA计算中的、使橡胶材料的应力-应变的关系反映到分析中的能量的结构方程式从线性弹性方程式向Mooney-Rivlin方程式过渡,近来,在该能量的结构方程式中导入了大变形区域内的非线性结构方程式。
另外,近年来,关于这些橡胶材料结构方程式,提出了以下结构方程式:以从基于橡胶的分子链拉伸的分子统计热力学发展出的网格变形理论为基础,能够使用具有物理意义的参数来表示在以轮胎为代表的很多工业用橡胶材料的设计中为重要因素的温度依赖性(日本专利4299735号公报)。
在此,下面示出该结构方程式的概要。
首先,使用剪切弹性率G和拉伸方向的拉伸比λ来如以下式(1)那样表示橡胶的拉伸应力。
[式1]
另一方面,如以下式(2)所示,根据统计热力学使用I1对亥姆霍兹自由能A进行微分来计算橡胶的剪切弹性率G。
[式2]
[式3]
在此,当k表示波尔兹曼常数、ΔT表示从高分子的玻璃化转变温度Tg至测量温度T为止的差时,β等于1/(kΔT)。另外,Aα表示能级。Z是使系统的能量标准化的分配函数,使用以下式(4)来表示。
[式4]
在此,如以下式(5)所示,根据统计热力学的想法,假设Uα与哈密顿算符H相等。
Uα=H(Temp,Constraints) ···(5)
其中,哈密顿算符H被表示为规定统计热力学中的微观状态的温度条件与限制条件的函数。
作为表示橡胶分子的温度依赖性的模型,设定低温状态与高温状态这两个状态,在I1a与I1b这种应变能量场中分配r1与r2这种橡胶分子的状态数。此时,当使用函数表示所有状态时,能够使用以下式(6)来表示。
[式5]
在此,κ表示分子的回弹能量。使用r1=-1、r2=1表示橡胶分子的能量状态不同的两种状态,因此作为统计热力学平均,表示为<r1×r2>。该积表示两个两极端的能量状态,在该积为+1的情况下表示橡胶分子为同一状态,在该积为-1的情况下表示橡胶分子为不同的状态。在此,当考虑为对所有的分子施加相同的应变能量场时,设为I1=I1a=I1b,由此能够如以下式(7)那样表示上述式(6)。
[式6]
将上述式(7)代入到上述式(3),由此如以下式(9)那样表示橡胶弹性率G。
[式7]
在此,在考虑分子的拉伸的情况下,使用郎之万函数如以下式(10)那样表示上述式(9)的第二项的熵项。
[式8]
在此,n被定义为交联点之间的统计分子链的链节数。
这样,在日本专利4299735号公报中,提出了能够使用具有物理意义的参数来表示温度依赖性的结构方程式。
另外,如图6所示,日本专利4299735号公报所记载的应变能量函数(实线)能够以宽幅的拉伸比与温度区来再现实验值(标绘点)。
发明内容
发明要解决的问题
然而,在日本专利4299735号公报所记载的应变能量函数中存在以下问题:在低温条件下的低拉伸比区域中,实验值与理论值偏离,因此无法高精度地预测低温条件下的橡胶制品的弹性响应性能。
本发明是用于解决上述问题而完成的,目的在于提供一种能够提高低温条件下的橡胶制品的弹性响应性能的预测精度的弹性响应性能预测方法、橡胶制品设计方法以及弹性响应性能预测装置。
用于解决问题的方案
为了达到上述目的,本发明的第一方式是,在预测表示橡胶制品的变形行为的弹性响应性能的弹性响应性能预测方法中,利用使用表示分子间的相互作用的参数来表示上述橡胶制品的应变能量的温度和应变依赖性的结构方程式,来预测该橡胶制品的弹性响应性能。
另外,本发明的第二方式是,在预测表示橡胶制品的变形行为的弹性响应性能的弹性响应性能预测装置中,利用使用表示分子间的相互作用的参数来表示上述橡胶制品的应变能量的温度和应变依赖性的结构方程式,来预测该橡胶制品的弹性响应性能。
根据本发明的第一方式和第二方式,在包括以微观级别观察填充类橡胶时的橡胶制品的低拉伸比区域(例如,拉伸比为100%以下)的广泛的变形区域中,利用使用表示分子间的相互作用的参数来表示应变能量的温度和应变依赖性的结构方程式,来预测橡胶制品的弹性响应性能。由此,能够高精度地预测低温条件下的低拉伸比区域的弹性响应性能。由此,特别是,对于实用上重要的由填充类的橡胶材料构成的橡胶复合制品,使用具有物理性意义的参数来表示在很多工业用橡胶制品的设计中重要的温度和应变依赖性,并且在低温条件下能够高精度地预测橡胶制品的弹性响应性能。
根据本发明的第三方式,上述结构方程式为下面示出的式(I),
[式9]
ΔA=(U1-TS1)+p(V1-V0)-(U0-TS0)… (I)
其中,A表示亥姆霍兹自由能,U0表示没有变形的状态下的内能,U1表示变形状态下的内能。p表示压力,V0表示没有变形的状态下的体积,V1表示变形状态下的体积,T表示绝对温度。S0表示没有变形的状态下的熵,S1表示变形状态下的熵。使用下面示出的式(II)~式(IV)来表示式(I)的各项。
[式10]
其中,I1、I2、I3作为橡胶的三维轴上的xyz方向的变形λ1、λ2、λ3这三个拉伸比的函数,使用I1=λ1 2+λ2 2+λ3 2、I2=λ1 2×λ2 2+λ2 2×λ3 2+λ3 2×λ1 2、I3=λ1 2×λ2 2×λ3 2来表示。n表示上述交联点之间的统计分子链的链节数,κ表示分子间相互作用能量系数。当将R设为气体常数、将Tg设为玻璃化转变温度时,β=1/RT,β’=1/R(T-Tg)。使用作为表示分子间的相互作用的参数的局部相互作用系数λmicro,以下面示出的式(V)来表示I1’。
[式11]
在本发明的第四方式中,在预测表示橡胶制品的变形行为的弹性响应性能的弹性响应性能预测方法中,利用使用表示分子间的相互作用的参数来表示上述橡胶制品的弹性率的温度和应变依赖性的结构方程式,来预测该橡胶制品的弹性响应性能。
在本发明的第五方式中,在预测表示橡胶制品的变形行为的弹性响应性能的弹性响应性能预测装置中,利用使用表示分子间的相互作用的参数来表示上述橡胶制品的弹性率的温度和应变依赖性的结构方程式,来预测该橡胶制品的弹性响应性能。
根据本发明的第四方式和第五方式,在包括以微观级别观察填充类橡胶时的橡胶制品的低拉伸比区域(例如,拉伸比为100%以下)的广泛的变形区域中,利用使用表示分子间的相互作用的参数来表示弹性率的温度和应变依赖性的结构方程式,来预测橡胶制品的弹性响应性能。由此,能够高精度地预测低温条件下的低拉伸比区域的弹性响应性能。由此,特别是,关于实用上重要的由填充类的橡胶材料构成的橡胶复合制品,使用具有物理性意义的参数来表示在很多工业用橡胶制品的设计中重要的温度和应变依赖性,并且在低温条件下能够高精度地预测橡胶制品的弹性响应性能。
根据本发明的第六方式,上述结构方程式为下面示出的式(VI)。
[式12]
其中,G表示剪切弹性率,W表示应变能量函数,A表示亥姆霍兹自由能。U表示内能,T表示绝对温度,S表示熵。I1作为橡胶的三维轴上的xyz方向的变形λ1、λ2、λ3这三个拉伸比的函数,使用I1=λ1 2+λ2 2+λ3 2来表示。使用下面示出的式(VII)、式(VIII)、式(IX)来表示式(VI)的各项。
[式13]
其中,n表示交联点之间的统计分子链的链节数,κ表示分子间相互作用能量系数,v表示交联密度,λ表示拉伸比或者压缩比。当将R设为气体常数、将Tg设为玻璃化转变温度时,β’=1/R(T-Tg)。使用I2=λ1 2×λ2 2+λ2 2×λ3 2+λ3 2×λ1 2来表示I2。使用作为表示分子间的相互作用的参数的局部相互作用系数λmicro,以下面示出的式(X)来表示I1’。
[式14]
根据本发明的第七方式,使用有限元分析法来预测表示上述橡胶制品的变形行为的弹性响应性能。
以上,通过详细说明的本发明的表示橡胶制品的应变能量的温度和应变依赖性的结构方程式或者表示橡胶制品的弹性率的温度和应变依赖性的结构方程式来预测该橡胶制品的弹性响应性能。由此,经过需要的仿真、特别是经过橡胶的微观级别的仿真,能够以微观级别来设计用于达到该橡胶制品的期望性能的最佳橡胶材料。另外,能够提供一种有效且精度良好的有用的橡胶制品的设计方法。
发明的效果
如以上说明的那样,根据本发明的弹性响应性能预测方法、橡胶制品设计方法以及弹性响应性能预测装置,利用使用表示分子间的相互作用的参数来表示应变能量的温度和应变依赖性的结构方程式或者表示弹性率的温度和应变依赖性的结构方程式,来预测橡胶制品的弹性响应性能。由此,得到以下效果:能够高精度地预测低温条件下的弹性响应性能。
附图说明
图1是表示本发明的实施方式所涉及的弹性响应性能的预测装置的概要图。
图2是用于导出本结构方程式的微结构模型的图。
图3是表示交联的丁苯橡胶的应变应力曲线的与温度依赖性有关的预测值和实测值的图表。
图4是表示交联的丁苯橡胶在室温时的应变应力曲线的、与能量弹性所引起的应力和熵弹性所引起的应力有关的计测值和计算值的图表。
图5是表示交联的天然橡胶在室温时的应变应力曲线的、与能量弹性所引起的应力和熵弹性所引起的应力有关的计测值和计算值的图表。
图6是表示以往技术中的橡胶的应变应力曲线的与温度依赖性有关的预测值和实测值的图表。
具体实施方式
下面,参照附图详细说明本发明的实施方式。
如图1所示,本发明的第一实施方式所涉及的执行弹性响应性能的预测的弹性响应性能的预测装置50由计算机运算处理系统构成,该计算机运算处理系统通过用于执行弹性响应性能预测的弹性响应性能的预测程序来执行后述的处理。此外,这种计算机系统例如具备CPU、ROM、RAM、硬盘、输入输出终端、其它需要的单元等。上述弹性响应性能的预测程序预先存储于硬盘等。
在本发明的第一实施方式的弹性响应性能的预测方法中,使用表示构成橡胶制品的橡胶材料的应变能量的温度和应变依赖性的结构方程式、特别是优选使用上述式(I)~(IV)来预测该橡胶制品的弹性响应性能。
在此,说明表示橡胶材料的应变能量的温度和应变依赖性的结构方程式的导出原理。下面,明确用于进行统计热力学计算的微结构模型并导出基于该模型的结构方程式。
图2表示微结构模型。Hint为表示分子间或者分子内相互作用的哈密顿算符,Hrot为表示分子旋转的哈密顿算符,Htrans为表示分子平移的哈密顿算符。
在此,能够使用以下式(11)来示出表示该微结构模型的热力学方程式。
[式15]
在此,A表示亥姆霍兹自由能,U表示内能,p表示压力。V表示体积,T表示绝对温度,S表示熵。另外,为表示分子间或者分子内相互作用的分配函数,为表示分子旋转的分配函数,为表示分子移动的分配函数。如以下式(12)~式(14)那样表示与各分子运动模式对应的哈密顿算符。
[式16]
在此,rl和r2表示分子间或者分子内相互作用的两个分子的能量状态。当设为I1’表示上述图2的微结构模型中的分子间的应变的不变量时,rl和r2与以+1或者-1表示的、I1’所表示的变形引起的应变能量的四个组合(-1,-1),(+1,+1),(+1,-1),(-1,+1)对应。κ为表示相互作用能量的大小的系数。δ为表示相互作用能量所干预的状态的系数,在两个分子处于相同能量状态时成为-1。这是由于,在两个分子处于相同能量状态时,κ×δ意味着相互作用能量相当稳定而为负的系数。另一方面,在处于不同的能量状态的情况下,假设不产生相互作用而κ×δ变为0。m表示分子的质量,Pθ表示相对于连结两个分子的轴的角度的、表示分子旋转的运动量。表示相对于以连结两个分子的轴为中心的旋转角度的分子旋转的运动量,P表示相对于质量m的分子的运动量,IM表示分子旋转的惯性力矩。
如以下式(15)那样表示哈密顿算符与分配函数的关系。
[式17]
通过将与各分子运动模式对应的哈密顿算符代入到使用上述(15)表示的关系式,由此如以下式(16)~式(18)那样表示各个分配函数。
[式18]
通过求解上述关系来求出以下式(19)~式(21)。
[式19]
这些分配函数为图2示出的局部模型,因此当拉伸到分子数N的宏观模型时,得到以下式(22)~式(24)。
[式20]
另一方面,如以下式(25)~式(27)那样表示分配函数与热力学的状态函数的关系。
[式21]
在此,β表示1/RT。当将Tg设为玻璃化转变温度时,β’表示以下式。
[式22]
n表示交联点之间的统计分子链的链节数,E表示分子间相互作用能量系数。Ω(F)表示保持分子的形态的数量,F表示变形的分子所产生的力量,a表示统计分子的长度,k表示波尔兹曼常数。通过将分子数N的分配函数代入到这些式,来求出上述式(I)。于是,由橡胶的变形引起的亥姆霍兹自由能的变化被表示为没有变形的状态与变形状态的差。
在上述式(I)中,在各项中,适用上述式(II)~(IV)。另外,考虑局部相互作用系数λmicro,使用上述式(V)来定义I1’。
在本实施方式中的弹性响应性能的预测方法中,使用以上述那样导出的式(I)~(V)表示的结构方程式来预测该橡胶制品的弹性响应性能。
另外,在本实施方式中,使用以下式(28)来定义交联点之间的统计分子链的数量n,并表示温度依赖性。由此,即使在应用于拉伸结晶性橡胶的微观级别的应变-应力行为分析技术的情况下,也能够考虑温度和应变依赖性。
[式23]
n=Δ·(1-Xc)·exp(-ε·β)… (28)
在此,α表示统计链段运动的频率因子,ε表示统计链段运动的活化能,Xc表示晶化率。仅在作为对象的材料示出拉伸结晶性的情况下,以内能的变化ΔU=U1-U0的函数用以下式(29)来表示Xc。
[式24]
其中,ΔH相当于结晶融化时的融解焓。在天然橡胶的情况下,作为ΔH,使用6MJ/m3。此外,在作为对象的材料不示出拉伸结晶性的情况下,将Xc设为0。另外,不会由于拉伸而100%结晶化,因而设为Xc不会变为1以上。
通过将上述式(28)、式(29)应用于式(I),内能项U在拉伸结晶性的橡胶中随着变形加大而示出负值(发热)。另外,熵S的分子的拉伸现象受到拉伸结晶化引起的交联点之间的分子运动性的限制,n的值与晶化成比例地变小,S增加。由此,得到表示在拉伸结晶化所引起的大变形区域中观察到的分子的拉伸或者非线性行为的结构方程式。
接着,说明第一实施方式所涉及的进行橡胶材料的变形行为的预测时的动作。
首先,通过计算机50来生成例如作为轮胎的橡胶制品的三维模型,三维模型数据被存储到HDD。
接着,用户操作输入输出终端来设定分析对象的三维模型。接着,用户操作输入输出终端来设定分析处理中的分析条件。分析条件的设定包括设定三维模型的橡胶部分的结构条件和填充剂部分的结构条件。作为三维模型的橡胶部分的结构条件,设定上述式(I)~(V)、式(28)的结构方程式(橡胶具有拉伸结晶化的情况)或者式(I)~(V)、式(28)、式(29)的结构方程式(橡胶不具有拉伸结晶化的情况)。
然后,计算机50从HDD读取所设定的分析对象的三维模型的三维模型数据,对由读取的三维模型数据表示的三维模型的橡胶部分和填充剂部分的结构条件附加所设定的分析条件来重新构成三维模型数据。
接着,计算机50使用重新构成的三维模型数据,根据加载到轮胎的车辆的负载、轮胎的内压、轮胎的滚动等,通过有限元法对在所设定的结构条件下改变三维模型时的三维模型的应变、内部应力分布以及三维模型整体的应力值进行分析。
接着,计算机50将通过分析求得的三维模型的应变状态、内部应力分布以及三维模型整体的应力值显示在显示器中而结束处理。
接着,使用作为上述结构方程式的式(I)~(V),预测橡胶的应变应力曲线的温度依赖性,并且使用交联的丁苯橡胶来测量应变应力曲线的温度依赖性。将通过上述预测得到的预测值(实线)与测量值(标绘点)进行比较。图3示出比较结果。可以说使用式(I)~(V)得到的曲线特别是在低拉伸比区域中示出与测量值的标绘点的良好的一致性。
另外,使用交联的丁苯橡胶,对于室温时的应变应力曲线,计测能量弹性所引起的应力和熵弹性所引起的应力,并且进行与上述结构方程式的式(I)、(II)相对应部分(U1、-TS1)的计算。将上述计测值(标绘点)与计算值(实线)进行比较。图4示出比较结果。可知通过式(I)、(II)的对应部分得到的曲线示出与计测值的标绘点的良好的一致性。
另外,使用交联的天然橡胶,对于室温时的应变应力曲线,计测能量弹性所引起的应力和熵弹性所引起的应力,并且应用式(28)、式(29)进行与上述结构方程式的式(I)、(II)相对应部分(U1、-TS1)的计算。将上述计测值(标绘点)与计算值(实线)进行比较。图5示出比较结果。可知通过式(I)、(II)的对应部分得到的曲线示出与计测值的标绘点的良好的一致性。由于天然橡胶示出拉伸结晶化,因此当拉伸比为3以上时存在能量弹性变为负的趋势。由于也良好地示出了该现象,因此可以说即使对于拉伸结晶化也能够实现良好的模型化。
如上所述,根据第一实施方式所涉及的橡胶制品的弹性响应性能的预测装置,在以微观级别观察填充类橡胶时的橡胶制品的广泛的变形区域中,通过使用加入了低应变时的临接分子的影响的局部相互作用系数λmicro来表示应变能量的温度和应变依赖性的结构方程式,来预测橡胶制品的弹性响应性能。由此,能够高精度地预测低温条件下的、特别是低拉伸比区域的弹性响应性能。
另外,定义微结构模型,加入看起来微小的相互作用的部分的误差、特别是低变形区域的能量弹性的贡献、拉伸结晶化这种现象,能够减少与实测值的偏离。因此,能够高精度地预测橡胶制品的弹性响应性能。
另外,除了表示以往的橡胶弹性的熵弹性的贡献以外,特征在于考虑能量弹性。在日本专利4299735号公报中提出的应变能量函数也同样地考虑了能量弹性,但是没有明确地定义其微结构模型。因此,无法表示看起来微小的相互作用的部分的误差、特别是低变形区域的能量弹性的贡献、拉伸结晶化这种现象。在本发明中,明确用于统计热力学计算的微结构模型,导出基于该模型的结构方程式,由此导出不依赖于变形区域、橡胶种类的普遍的橡胶弹性的结构方程式。
在以往的橡胶的结构方程式中,无法表现出上述图4、图5示出的橡胶的热力学行为。在本发明中,将表现上述图4、图5所示那样橡胶的热力学行为的结构方程式导入到以往的计算科学仿真。由此,能够在反映了现实中使用的原材料的力学行为的状态下进行预测计算。
另外,明确了对于使用通过统计热力学计算得到的结构方程式来表示的应力应变曲线的温度依赖性,包括具有拉伸结晶性的天然橡胶在内与通过实测得到的结果良好地一致。由此,通过将所提出的结构方程式导入到有限元法的应力应变曲线,在拉伸结晶性以及非结晶性的橡胶两者的纳米级的橡胶的应力应变行为中,能够表示包括其弹性贡献、粘性贡献、塑性贡献在内的行为。
接着,说明第二实施方式。
在第二实施方式中,使用表示橡胶制品的弹性率的温度和应变依赖性的结构方程式来预测橡胶制品的弹性响应性能这一点与第一实施方式不同。
在第二实施方式的弹性响应性能的预测方法中,将与上述第一实施方式同样地导出的、使用上述式(VI)来表示的结构方程式应用于橡胶的微观级别的变形行为分析技术,来预测该橡胶制品的弹性响应性能。
使用上述式(VII)、(VIII)、(IX)来表示上述式(VI)中的各项。其中,使用局部相互作用系数λmicro,以上述式(X)来表示I1’。
另外,在本实施方式中,使用上述式(28)、式(29)来定义交联点之间的统计分子链的数量n,并表示拉伸结晶化。由此,即使应用于拉伸结晶性橡胶的微观级别的变形行为分析技术的情况下,也能够考虑温度和应变依赖性。
此外,第二发明的实施方式所涉及的弹性响应性能的预测装置的其它结构和作用与第一实施方式相同,因此省略说明。
这样,根据第二实施方式所涉及的橡胶制品的弹性响应性能的预测装置,在以微观级别观察填充类橡胶时的橡胶制品的广泛的变形区域中,通过使用加入了低应变时的临接分子的影响的局部相互作用系数λmicro来表示弹性率的温度和应变依赖性的结构方程式,来预测橡胶制品的弹性响应性能。由此,能够高精度地预测低温条件下的、特别是低拉伸比区域的弹性响应性能。
另外,定义微结构模型,加入看起来微小的相互作用的部分的误差、特别是低变形区域的能量弹性的贡献、拉伸结晶化这种现象,能够减少与实测值的偏离。因此,能够高精度地预测橡胶制品的弹性响应性能。
日本申请2010-225985的公开整体通过参照取入到本说明书中。
关于本说明书所记载的所有文献、专利申请以及技术标准,各个文献、专利申请以及技术标准通过参照被取入的情况与具体且分别记载的情况相同程度地,通过参照被取入到本说明书中。
Claims (8)
1.一种弹性响应性能预测方法,预测表示橡胶制品的变形行为的弹性响应性能,在该弹性响应性能预测方法中,
利用使用表示分子间的相互作用的参数来表示上述橡胶制品的应变能量的温度和应变依赖性的结构方程式,来预测该橡胶制品的弹性响应性能。
2.根据权利要求1所述的弹性响应性能预测方法,其特征在于,
上述结构方程式为下面示出的式(I),
[式1]
ΔA=(U1-TS1)+p(V1-V0)-(U0-TS0)… (I)
其中,A表示亥姆霍兹自由能,U0表示没有变形的状态下的内能,U1表示变形状态下的内能,p表示压力,V0表示没有变形的状态下的体积,V1表示变形状态下的体积,T表示绝对温度,S0表示没有变形的状态下的熵,S1表示变形状态下的熵,使用下面示出的式(II)~式(IV)来表示式(I)的各项,
[式2]
其中,I1、I2、I3作为橡胶的三维轴上的xyz方向的变形λ1、λ2、λ3这三个拉伸比的函数,以I1=λ1 2+λ2 2+λ3 2、I2=λ1 2×λ2 2+λ2 2×λ3 2+λ3 2×λ1 2、I3=λ1 2×λ2 2×λ3 2来表示,n表示交联点之间的统计分子链的链节数,κ表示分子间相互作用能量系数,当将R设为气体常数、将Tg设为玻璃化转变温度时,β=1/RT,β’=1/R(T-Tg),使用作为表示上述分子间的相互作用的参数的局部相互作用系数λmicro,以下面示出的式(V)来表示I1’,
[式3]
3.一种弹性响应性能预测方法,预测表示橡胶制品的变形行为的弹性响应性能,在该弹性响应性能预测方法中,
利用使用表示分子间的相互作用的参数来表示上述橡胶制品的弹性率的温度和应变依赖性的结构方程式,来预测该橡胶制品的弹性响应性能。
4.根据权利要求3所述的弹性响应性能预测方法,其特征在于,
上述结构方程式为下面示出的式(VI),
[式4]
其中,G表示剪切弹性率,W表示应变能量函数,A表示亥姆霍兹自由能,U表示内能,T表示绝对温度,S表示熵,I1作为橡胶的三维轴上的xyz方向的变形λ1、λ2、λ3这三个拉伸比的函数,以I1=λ1 2+λ2 2+λ3 2来表示,以下面示出的式(VII)、式(VIII)、式(IX)来表示式(VI)的各项,
[式5]
其中,n表示交联点之间的统计分子链的链节数,κ表示分子间相互作用能量系数,v表示交联密度,λ表示拉伸比或者压缩比,当将R设为气体常数、将Tg设为玻璃化转变温度时,β’=1/R(T-Tg),使用I2=λ1 2×λ2 2+λ2 2×λ3 2+λ3 2×λ1 2来表示I2,使用作为表示上述分子间的相互作用的参数的局部相互作用系数λmicro,以下面示出的式(X)来表示I1’,
[式6]
5.根据权利要求1~4中的任一项所述的弹性响应性能预测方法,其特征在于,
使用有限元分析法来预测表示上述橡胶制品的变形行为的弹性响应性能。
6.一种橡胶制品设计方法,使用根据权利要求1~5中的任一项所述的弹性响应性能预测方法来设计橡胶制品。
7.一种弹性响应性能预测装置,预测表示橡胶制品的变形行为的弹性响应性能,在该弹性响应性能预测装置中,
利用使用表示分子间的相互作用的参数来表示上述橡胶制品的应变能量的温度和应变依赖性的结构方程式,来预测该橡胶制品的弹性响应性能。
8.一种弹性响应性能预测装置,预测表示橡胶制品的变形行为的弹性响应性能,在该弹性响应性能预测装置中,
利用使用表示分子间的相互作用的参数来表示上述橡胶制品的弹性率的温度和应变依赖性的结构方程式,来预测该橡胶制品的弹性响应性能。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010225985A JP5474726B2 (ja) | 2010-10-05 | 2010-10-05 | ゴム製品の弾性応答性能の予測方法、設計方法、及び弾性応答性能予測装置 |
JP2010-225985 | 2010-10-05 | ||
PCT/JP2011/072887 WO2012046740A1 (ja) | 2010-10-05 | 2011-10-04 | ゴム製品の弾性応答性能の予測方法、設計方法、及び弾性応答性能予測装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103154703A true CN103154703A (zh) | 2013-06-12 |
CN103154703B CN103154703B (zh) | 2015-04-01 |
Family
ID=45927735
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180048359.0A Expired - Fee Related CN103154703B (zh) | 2010-10-05 | 2011-10-04 | 橡胶制品的弹性响应性能的预测方法、设计方法以及弹性响应性能预测装置 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9405869B2 (zh) |
EP (1) | EP2626688B1 (zh) |
JP (1) | JP5474726B2 (zh) |
CN (1) | CN103154703B (zh) |
WO (1) | WO2012046740A1 (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105989244A (zh) * | 2016-05-10 | 2016-10-05 | 中国人民解放军国防科学技术大学 | 橡胶材料的超弹性本构模型建模方法 |
CN106295021A (zh) * | 2016-08-15 | 2017-01-04 | 浙江大学 | 一种临氢厚壁圆柱壳弹性应力应变的预测方法 |
CN106372273A (zh) * | 2016-08-15 | 2017-02-01 | 浙江大学 | 一种临氢厚壁圆柱壳弹性极限载荷预测方法 |
CN106612760A (zh) * | 2016-12-12 | 2017-05-10 | 贵州省烟草科学研究院 | 便携式螺旋背负打孔机 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5913260B2 (ja) * | 2013-11-14 | 2016-04-27 | 住友ゴム工業株式会社 | 高分子材料のシミュレーション方法 |
JP7088748B2 (ja) | 2018-05-29 | 2022-06-21 | 住友重機械工業株式会社 | 溶接装置及び溶接方法 |
CN110376054A (zh) * | 2019-06-06 | 2019-10-25 | 中国船舶重工集团公司第七一九研究所 | 一种基于化学反应的材料力学性能统一表征方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2007265266A (ja) * | 2006-03-29 | 2007-10-11 | Bridgestone Corp | 粘弾性応答性能予測方法、ゴム製品設計方法、及び粘弾性応答性能予測装置 |
JP4299735B2 (ja) * | 2004-06-07 | 2009-07-22 | 株式会社ブリヂストン | ゴム製品の弾性応答性能の予測方法及び設計方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010225985A (ja) | 2009-03-25 | 2010-10-07 | Fujifilm Corp | 光電変換半導体層とその製造方法、光電変換素子、及び太陽電池 |
US8155933B2 (en) * | 2009-05-26 | 2012-04-10 | King Fahd University Of Petroleum & Minerals | Method of modeling residual stresses during laser cutting |
-
2010
- 2010-10-05 JP JP2010225985A patent/JP5474726B2/ja not_active Expired - Fee Related
-
2011
- 2011-10-04 EP EP11830672.9A patent/EP2626688B1/en not_active Not-in-force
- 2011-10-04 CN CN201180048359.0A patent/CN103154703B/zh not_active Expired - Fee Related
- 2011-10-04 WO PCT/JP2011/072887 patent/WO2012046740A1/ja active Application Filing
- 2011-10-04 US US13/877,607 patent/US9405869B2/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4299735B2 (ja) * | 2004-06-07 | 2009-07-22 | 株式会社ブリヂストン | ゴム製品の弾性応答性能の予測方法及び設計方法 |
JP2007265266A (ja) * | 2006-03-29 | 2007-10-11 | Bridgestone Corp | 粘弾性応答性能予測方法、ゴム製品設計方法、及び粘弾性応答性能予測装置 |
Non-Patent Citations (4)
Title |
---|
A.D.DROZDOV: ""A model for the mechanical response of composites with thermoplastic-elastomer matrices"", 《COMPOSITES SCIENCE AND TECHNOLOGY》, vol. 66, 31 December 2006 (2006-12-31), pages 2648 - 2663, XP 025128662, DOI: doi:10.1016/j.compscitech.2006.03.017 * |
张赟等: ""剪切变形下非晶态高聚物的力学行为"", 《应用数学和力学》, vol. 25, no. 10, 31 October 2004 (2004-10-31), pages 997 - 1006 * |
杨小震: ""链分子的构象弹性理论"", 《中国科学(B辑)》, vol. 31, no. 1, 28 February 2001 (2001-02-28), pages 78 - 87 * |
王利荣等: ""橡胶隔振器有限元建模技术及静态弹性特性分析"", 《汽车工程》, vol. 24, no. 6, 31 December 2002 (2002-12-31), pages 480 - 485 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105989244A (zh) * | 2016-05-10 | 2016-10-05 | 中国人民解放军国防科学技术大学 | 橡胶材料的超弹性本构模型建模方法 |
CN106295021A (zh) * | 2016-08-15 | 2017-01-04 | 浙江大学 | 一种临氢厚壁圆柱壳弹性应力应变的预测方法 |
CN106372273A (zh) * | 2016-08-15 | 2017-02-01 | 浙江大学 | 一种临氢厚壁圆柱壳弹性极限载荷预测方法 |
CN106295021B (zh) * | 2016-08-15 | 2018-05-29 | 浙江大学 | 一种临氢厚壁圆柱壳弹性应力应变的预测方法 |
CN106372273B (zh) * | 2016-08-15 | 2018-05-29 | 浙江大学 | 一种临氢厚壁圆柱壳弹性极限载荷预测方法 |
CN106612760A (zh) * | 2016-12-12 | 2017-05-10 | 贵州省烟草科学研究院 | 便携式螺旋背负打孔机 |
CN106612760B (zh) * | 2016-12-12 | 2023-03-21 | 贵州省烟草科学研究院 | 便携式螺旋背负打孔机 |
Also Published As
Publication number | Publication date |
---|---|
EP2626688B1 (en) | 2019-02-20 |
WO2012046740A1 (ja) | 2012-04-12 |
EP2626688A1 (en) | 2013-08-14 |
CN103154703B (zh) | 2015-04-01 |
EP2626688A4 (en) | 2016-05-25 |
JP5474726B2 (ja) | 2014-04-16 |
US9405869B2 (en) | 2016-08-02 |
US20130211802A1 (en) | 2013-08-15 |
JP2012078296A (ja) | 2012-04-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103154703B (zh) | 橡胶制品的弹性响应性能的预测方法、设计方法以及弹性响应性能预测装置 | |
Kamarian et al. | Free vibration analysis of conical shells reinforced with agglomerated Carbon Nanotubes | |
Malekzadeh et al. | Low velocity impact analysis of functionally graded carbon nanotubes reinforced composite skew plates | |
Tan et al. | Thermal mechanical constitutive model of fiber reinforced shape memory polymer composite: based on bridging model | |
Huisman et al. | Monte Carlo study of multiply crosslinked semiflexible polymer networks | |
Österlöf et al. | A finite strain viscoplastic constitutive model for rubber with reinforcing fillers | |
Assidi et al. | Equivalent mechanical properties of biological membranes from lattice homogenization | |
CN103154702B (zh) | 橡胶制品的弹性响应性能的预测方法、设计方法以及弹性响应性能预测装置 | |
JP6048358B2 (ja) | 解析装置 | |
Bazkiaei et al. | A framework for model base hyper-elastic material simulation | |
Yerezhep et al. | A multifaceted approach for cryogenic waste tire recycling | |
Wang et al. | Size effects of elastic properties for auxetic cellular structures: bending energy-based method | |
JP4299735B2 (ja) | ゴム製品の弾性応答性能の予測方法及び設計方法 | |
Yamashita et al. | Calculation of Strain Energy Density Function Using Ogden Model and Mooney–Rivlin Model Based on Biaxial Elongation Experiments of Silicone Rubber | |
JP2007265266A (ja) | 粘弾性応答性能予測方法、ゴム製品設計方法、及び粘弾性応答性能予測装置 | |
CN106096258B (zh) | 一种含氟聚氨酯表面疏水性的分析方法 | |
Deng et al. | Moment Lyapunov exponents and stochastic stability of coupled viscoelastic systems driven by white noise | |
Pierce et al. | Implementation of a non-orthogonal constitutive model for the finite element simulation of textile composite draping | |
Yang et al. | On the energy conservation during the active deformation in molecular dynamics simulations | |
CN112287465B (zh) | 一种面内轮胎柔性环模型实时仿真的实现方法 | |
af Segerstad et al. | A constitutive equation for open-cell cellular solids, including viscoplasticity, damage and deformation induced anisotropy | |
Ntetsika et al. | A hybrid immersed boundary/coarse-graining method for modeling inextensible semi-flexible filaments in thermally fluctuating fluids | |
Das et al. | On Mechanical Properties of Graphene Sheet Estimated Using Molecular Dynamics Simulations | |
Shepherd et al. | Review of Modelling and Additive Manufacturing of Auxetic Materials for Application in Sport, 8th Asia-Pacific Congress on Sports Technology | |
Akutagawa | Nano-mechanical modeling for rubbery materials |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150401 |