本发明专利技术公开了一种各向同性不可压缩超弹性体本构模型及其应用，属于固体力学、计算力学和实验力学技术领域。本发明专利技术提出的本构模型只采用单一类型实验数据就能获得较准确的橡胶类不可压缩超弹性材料的模型，比当前已有模型精度和可靠性都高，且可以仅仅进行简单的单轴拉伸实验测试获得较准确、全面的橡胶材料特性模型，不必采用当前国内尚难以进行的等双轴拉伸和平面拉伸测试，具有极高的工程应用价值和计算力学使用前景。
【技术实现步骤摘要】
一种各向同性不可压缩超弹性体本构模型及其应用
本专利技术涉及工程和科学研究中对各类橡胶材料、生物组织材料的力学性能进行描述和建模，以及在此基础上的工程和科学应用，特别是计算力学和实验力学的相关领域。
技术介绍
诸如橡胶、肌肉、韧带等材料可以用不可压缩超弹性模型来描述其力学特性，建立一个能完整描述该类材料力学的本构模型是工程和科学研究的重大问题。当前虽然已经提出来基于应变的不变量模型、基于三个主伸长率的模型，以及基于分子链的模型，但上述橡胶、生物组织类材料缺乏一个完备而且简洁的模型，具体表现在：或者已有模型的描述能力不够，不能完全描述材料在各种变形状态下的力学特性；或者模型过于复杂，在实际的使用中比较困难。此外，通常情况下，建立橡胶类不可压缩超弹性材料的模型至少需要单轴拉伸、平面拉伸和等双轴拉伸三类实验，但当前只有第一种实验即单轴拉伸实验在实践和理论上比较成熟，接近理想状态，后两种实验理论上还存在一定的缺陷，同时由于实验设备与条件的限制，实验非常困难，国内基本上没有实验室能较好的进行该两类实验，基本上要依托国外如北美AXEL专业实验室才能进行，实验成本和费用都非常高。随着工业的发展，对橡胶类不可压缩超弹性材料力学性能模型的要求越来越高，经检索，中国专利申请号201610303626.9，申请日为2016年5月10日，专利技术创造名称为：橡胶材料的超弹性本构模型建模方法；中国专利申请号201610519914.8，申请日为2016年7月4日，专利技术创造名称为：一种考虑相关效应的橡胶材料粘超弹本构模型的建模方法；上述申请案均提出了构建材料本构模型的方法(本构模型参数识别方法)，并未提出本构模型本身，目前国内尚未有公开的本构关系申请案。
技术实现思路
1.专利技术要解决的技术问题为了解决上述现有技术存在的问题，本专利技术提出了一种各向同性不可压缩超弹性体本构模型及其应用，本专利技术提出的模型不仅能很好的拟合所有实验曲线，且只用一种类型的材料实验就可获得较准确地任意状态下可靠的材料力学特性，可有效解决工程实践中橡胶类材料平面拉伸和等双轴拉伸实验困难，及已有模型拟合精度低，不可靠的问题。2.技术方案为达到上述目的，本专利技术提供的技术方案为：本专利技术的一种各向同性不可压缩超弹性体本构模型，为任意平面变形状态的应力-伸长率函数与平面拉伸应力-伸长率函数的关系模型，具体为如下形式：式中，Ti表示在任意三轴应力状态(λi,λj,λk)下λi方向的名义应力；Ti(λi,λj)表示在平面应力状态(λi,λj,1/λiλj)下λi方向的名义应力；λi,λj,λk分别为三个受力方向的主伸长率，下标i,j,k为数字1、2、3的组合排列；k为一个待识别的材料参数，k在λj大于1或者小于1时可取不同的两个值；Tplanar表示为平面拉伸材料力学性能试验中测得的在伸长率为λ时的名义应力函数。更进一步地，在三轴应力状态下，所述本构模型为：式中，p为静水压力，需要结合材料的边界条件确定。更进一步地，在单轴拉伸实验中，单轴拉伸名义应力函数Tuniaxial(λ)与平面拉伸应力函数Tplanar(λ)用以下关系转换：更进一步地，在等双轴拉伸实验中，等双轴拉伸的名义应力函数Tbiaxial(λ)和平面拉伸应力函数Tplanar(λ)用以下关系转换：式中，Tbiaxial(λ)为等双轴拉伸实验中在伸长率为λ时的名义应力函数，g为一个材料常数。本专利技术的一种各向同性不可压缩超弹性体本构模型的应用，使用任何一种实验数据，拟合材料参数，借助上述模型建立各向同性不可压缩超弹性体本构模型用于工程计算、数值仿真或者材料实验。3.有益效果采用本专利技术提供的技术方案，与已有的公知技术相比，具有如下显著效果：(1)本专利技术以两个主伸长率为模型变量，建立了各种实验类型数据之间的关系，在此基础上构建了一个全新的本构模型，该模型是一种不可压缩超弹性模型，可直接用于各向同性不可压缩超弹性材料，也可以作为基本模型嵌入到粘超弹性模型中使用，同时可描述材料的各向异性，可靠度大于已有模型；(2)本专利技术的一种各向同性不可压缩超弹性体本构模型，使用一种实验类型实验数据就能可靠地预测橡胶类材料所有变形状态下的力学响应，不必采用当前国内尚难以进行的等双轴拉伸和平面拉伸测试，节省了昂贵的材料多轴测试费用，具有极高的工程应用价值和计算力学使用前景，是橡胶类不可压缩超弹性材料本构特性的一个新发现。附图说明图1为本专利技术的各向同性不可压缩超弹性体本构模型拟合结果示意图；图2为实施例1获得的名义应力-伸长率曲线(Treloar)。具体实施方式为进一步了解本专利技术的内容，结合附图和实施例对本专利技术作详细描述。针对已有的超弹性本构模型或者不能很好地同时拟合所有类型的实验曲线(即拟合精度低)；或者不能可靠地由一种类型实验数据(譬如仅用单轴拉伸实验数据)准确预测其它类型的实验数据(如等双轴和平面拉伸实验数据)，模型可靠度低的问题，本专利技术提出了一种新的各向同性不可压缩超弹性体本构模型及其参数识别方法。实施例1本实施例提出了任意平面变形状态的应力-伸长率函数与平面拉伸应力-伸长率函数的关系模型，通过该关系模型，可直接建立橡胶类不可压缩超弹性材料本构关系并识别模型参数，具体为如下形式：式中，Ti表示在任意三轴应力状态(λi,λj,λk)下λi方向的名义应力；Ti(λi,λj)表示在平面应力状态(λi,λj,1/λiλj)下λi方向的名义应力(力/初始面积)；λi,λj,λk分别为三个受力方向的主伸长率(λi,λj在平面应力状态下即为第一和第二主伸长率)，下标i,j,k为数字1、2、3的一个排列；k为一个待识别的材料参数，在一般情况可以令其等于0，在伸长率λj大于1和小于1时可取同样的值或者不一样的值；Tplanar表示为平面拉伸材料力学性能试验中测得的在伸长率为λ时的名义应力函数。式(1)为本实施例提出的本构模型的最一般形式，描述了在平面应力状态下的橡胶类不可压缩超弹性材料本构关系。该模型以平面拉伸应力函数Tplanar为基函数，该函数可以采取任何合适的表达形式，包括多项式函数、指数函数、幂函数等等。对于更一般的三轴应力状态，所述本构关系为：式中，p为静水压力，需要结合材料的边界条件才能确定。而在单轴拉伸实验中，单轴拉伸名义应力函数Tuniaxial(λ)与平面拉伸应力函数Tplanar(λ)用以下关系转换：在等双轴拉伸实验中，等双轴拉伸的名义应力和平面拉伸的应力关系为：式中，Tbiaxial(λ)为等双轴拉伸实验中在伸长率为λ时的名义应力函数，g为一个材料常数。本实施例提出的本构模型可以单独用任何一种实验数据(单轴拉伸、平面拉伸或等双轴拉伸)拟合参数，然后利用该模型推出任意变形状态下的应力响应，具有模型精度高、可靠性好的特点，由于本实施例提出了包括单轴拉伸实验、平面拉伸实验和等双轴拉伸实验这三个实验状态在内的任意双轴拉伸实验状态中应力函数相互转换的模型，只要得到一种实验数据，利用所述模型就可推知另外实验状态下的材料应力响应，不必像现有方法中需要分别进行上述三种类型的材料力学实验，然后同时用这三种实验数据来拟合模型。通过对实验数据的拟合显示，只采用一类实验数据如单轴拉伸实验数据，就能较准确地得出材料在其它变形状态下的应力。且本实施例提出的本本文档来自技高网...
【技术保护点】
一种各向同性不可压缩超弹性体本构模型，其特征在于：该本构模型为任意平面变形状态的应力‑伸长率函数与平面拉伸应力‑伸长率函数的关系模型，具体为如下形式：
【技术特征摘要】
1.一种各向同性不可压缩超弹性体本构模型，其特征在于：该本构模型为任意平面变形状态的应力-伸长率函数与平面拉伸应力-伸长率函数的关系模型，具体为如下形式：式中，Ti表示在任意三轴应力状态(λi,λj,λk)下λi方向的名义应力；Ti(λi,λj)表示在平面应力状态(λi,λj,1/λiλj)下λi方向的名义应力；λi,λj,λk分别为三个受力方向的主伸长率，下标i,j,k为数字1、2、3的组合排列；k为一个待识别的材料参数，k在λj大于1或者小于1时可取不同的两个值；Tplanar表示为平面拉伸材料力学性能试验中测得的在伸长率为λ时的名义应力函数。2.根据权利要求1所述的一种各向同性不可压缩超弹性体本构模型，其特征在于：在三轴应力状态下，所述本构模型为：式中，p为静水压力，需要结合材料...
【专利技术属性】
技术研发人员：魏志刚，王孝义，
申请(专利权)人：安徽工业大学，
类型：发明
国别省市：安徽,34
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。