ABAQUS-材料本构模型及编程

合集下载

橡胶材料在ABAQUS的材料参数设定

橡胶材料在ABAQUS的材料参数设定

橡胶材料在ABAQUS的材料参数设定ABAQUS是一款常用的有限元分析软件,能够进行多种工程问题的模拟和分析。

在ABAQUS中,要设定橡胶材料的材料参数,首先需要选择适当的材料模型,并根据实验数据来确定材料参数的具体数值。

橡胶材料的性质是非线性的,所以在ABAQUS中通常使用Hyperelastic材料模型。

下面将详细介绍橡胶材料在ABAQUS中的材料参数设定。

橡胶材料的本构模型由于橡胶材料的高度可压缩性和非线性行为,经典的线性弹性模型不能准确地描述橡胶的力学性能。

在ABAQUS中,默认的橡胶材料模型是非线性的Hyperelastic材料模型,可选的模型包括:Mooney-Rivlin模型、Neo-Hookean模型、Ogden模型等。

这些模型的主要区别在于其形式和需要确定的参数数量。

在选择合适的模型时,需要根据实验数据的特点和需求来进行选择。

材料参数的确定确定橡胶材料的材料参数是非常重要的,这些参数直接影响到模拟结果的准确性。

通常,可以通过实验测试来测量材料的拉伸或压缩行为,以及其它的力学性能,例如剪切刚度和各个方向上的应变能函数。

利用这些实验数据,可以利用ABAQUS提供的拟合工具进行参数拟合,从而得到合理的材料参数。

拟合工具ABAQUS提供了多种实验数据拟合工具,用于确定材料模型的参数。

其中最常用的是通过拉伸实验数据进行拟合来确定材料的应变能函数。

该方法基于ABAQUS的材料模型来计算应变能函数,然后将实验数据拟合到计算结果得到最佳拟合参数。

在ABAQUS中,可以通过以下步骤进行材料参数设定:1. 创建材料模型:选择合适的Hyperelastic材料模型,并为其分配一个名称。

2.确定材料参数:根据实验数据的特点和要求,选择适当的材料参数。

3.输入材料参数:将确定的材料参数输入到ABAQUS中,可以通过输入文件或者ABAQUS/CAE图形界面进行设定。

4.材料测试:使用所设定的材料参数进行模拟测试,验证材料模型的准确性。

基于Abaqus子程序的高分子材料本构关系实现

基于Abaqus子程序的高分子材料本构关系实现

基于Abaqus子程序的高分子材料本构关系实现作者:汪品红来源:《计算机辅助工程》2013年第05期摘要:对于高分子材料的仿真,业界一般使用经典的弹塑性本构模型来描述其应力应变关系,但其真实的应力应变关系与经典的弹塑性本构模型存在一定差异,从而导致仿真与实际测试之间的差异.Abaqus提供UMAT/VUMAT子程序接口,让用户可以自己构建新的材料本构模型.通过撰写新的材料本构子程序实现高分子材料应力应变关系在仿真中的准确描述,减少仿真与测试之间的差异.同时,在卸载段可以通过卸载标志符的选择定义不同的卸载路径,方便用户使用.关键词:高分子材料;本构关系; Abaqus; UMAT; VUMAT中图分类号: TB324; TB115.1文献标志码: B引言高分子材料在日常生活中有着广泛的应用,因此其不可避免地出现在仿真分析中.当前没有一种商业软件具有适合高分子材料的材料本构模型.Abaqus是一款优秀的商业软件,其提供的子程序接口UMAT/VUMAT允许用户根据使用需求自定义材料本构.[1]使用该方法,可有效解决在仿真中由于材料本构不适用而导致的仿真与实际测试差异过大的问题.1高分子材料本构一般描述方法业界通常使用弹塑性本构定义高分子材料的材料属性.屈服强度一般取材料曲线上第一个峰值点.弹性模量的取法有2种不同的方式:对于应力应变关系曲线有明显直线段的,以第一段直线的斜率作为材料的弹性模量(切线法);对于曲线没有明显直线段的材料,则使用原点与屈服点连成的直线的斜率作为弹性模量(割线法).2种方式与真实应力应变曲线的比较见图1.图 1高分子材料测试材料曲线与仿真曲线比较由图1可知,无论使用何种方式,仿真使用的应力应变曲线都与实际材料的应力应变曲线有较大差异.将切线法获得材料数据代入到手机电池盖三点弯曲中进行仿真,见图2,其仿真与测试力位移曲线在最高点的差异约为23%,见图3.对于手机等一些电子类产品,高分子材料的仿真非常重要.在跌落或弯折测试中,高分子材料的应力应变关系与弹塑性本构的差异造成仿真预测不准确,必须定义正确的高分子材料本构.2Abaqus VUMAT子程序Abaqus提供丰富的材料本构模型库,能够满足绝大多数仿真材料模型的需要;同时,还提供UMAT/VUMAT子程序接口,让用户可以用FORTRAN语言编程,自己定义需要的材料本构模型,对Abaqus材料库中没有包含的材料进行计算.几乎可以把用户材料属性赋予Abaqus 中的任何单元,其中UMAT用在隐式仿真计算中,VUMAT用在显式仿真计算中.由于隐式计算与显式计算的差别,导致UMAT与VUMAT也有一定的差异,但是经过简单的改写即可完成它们之间的转换.本文使用准静态仿真分析方法,属于显式求解,所以只介绍VUMAT.3高分子材料VUMAT本构介绍由图1可知,高分子材料的本构与弹塑性本构最大的差异在于弹性段是直线还是曲线.弹性段的路径也直接影响到卸载的路径.因此,对高分子材料本构的定义关键在于非线性弹性段的实现,即要根据当前的应力值实时获取下一增量步所用的弹性模量值.程序整体流程见图4.图 4程序整体流程3.1弹性段多段线性的实现在弹性段,程序根据弹性模量和泊松比计算应力增量.由于弹性段为非线性,需要根据应力或应变更新用于计算的弹性模量值,直至达到屈服点,因此需要在输入文件中输入材料真实应力应变曲线,通过查表计算的函数,根据当前应力σ所在的位置,计算当前的弹性模量.应力应变曲线输入时,输入格式为:用查表的方法,直到σn3.2卸载路径的选择屈服发生后,需要选择弹性模量参与相关计算,有2个作用:一是用来计算屈服后加载段的应力试探值(不对该增量步真实应力产生影响,只起对比判断的作用);二是用来作为屈服后卸载的路径(为实现不同卸载路径,在程序中设置一个flag位,其值由用户自己输入),用户可以根据实际的需要选择卸载的路径.如图4中,共设置3种卸载路径:沿切线卸载、沿割线卸载以及沿曲线卸载等.用户也可以根据需要增加其他的卸载方式.4子程序的验证为验证子程序是否能实现设计的功能,取一个1/8的网格模型进行单轴拉伸仿真,单元类型为C3D8R.输出其应力应变曲线,与材料真实应力应变曲线比较,见图5.图 5使用VUMAT后加载应力应变曲线与材料曲线对比使用VUMAT后,加载的应力应变曲线与材料测试得到的真实应力应变曲线完全重合,说明VUMAT可以完全反映材料在加载过程中的力学行为.在卸载过程中,分别实现沿弹性段的切线、割线以及曲线卸载.为进一步验证,将VUMAT用于图2所示的手机电池盖三点弯模型中进行仿真与试验对比.在使用弹塑性本构模型时,仿真与测试力位移曲线的最大差异约为23%,而引入使用VUMAT 编写的高分子材料本构后,其仿真与测试的差异减少到4.5%,见图6.从实际项目的验证结果看,使用VUMAT后电池盖测试的力位移曲线与仿真的力位移曲线基本重合,仿真与测试的差异也明显减小.将该本构应用于其他高分子材料和实际案例,其仿真精度均明显改善,也说明该子程序在实际工程中的适用性.图 6使用VUMAT后电池盖力位移曲线对比5结束语使用VUMAT子程序后,高分子材料在加载段的力学特性与测试的真实应力应变曲线一致,同时将其应用在工程实际问题上,也与测试曲线基本一致,验证该程序的适用性.由于高分子材料的卸载特性较为复杂,还需进一步研究,所以程序只给出3种方式供用户按照实际需求进行选择.参考文献:[1]庄茁,张帆,岑松,等. 基于Abaqus的有限元分析和应用[M]. 北京:清华大学出版社, 2009: 509512.(编辑武晓英)。

abaqus的多晶硅本构模型

abaqus的多晶硅本构模型

abaqus的多晶硅本构模型多晶硅是一种常见的半导体材料,广泛应用于电子元器件、太阳能电池等领域。

在使用ABAQUS软件进行多晶硅材料仿真时,需要选择合适的本构模型来描述其力学行为。

本文将介绍ABAQUS中多晶硅的本构模型及其特点。

在ABAQUS中,多晶硅的本构模型主要有两种:弹塑性模型和弹性模型。

弹塑性模型考虑了多晶硅材料的塑性变形,而弹性模型则仅考虑弹性变形。

我们来介绍多晶硅的弹塑性模型。

ABAQUS中常用的弹塑性模型有Hypoelastic模型和Von Mises塑性模型。

Hypoelastic模型将多晶硅的塑性行为近似为弹性行为,通过引入有效应力和等效塑性应变来描述材料的塑性变形。

Von Mises塑性模型则根据材料的屈服准则来确定塑性区域,并引入等效应力和等效塑性应变来描述塑性行为。

我们来介绍多晶硅的弹性模型。

ABAQUS中常用的弹性模型有线性弹性模型和非线性弹性模型。

线性弹性模型假设多晶硅的应力和应变之间存在线性关系,即满足胡克定律。

非线性弹性模型则考虑了多晶硅材料的非线性行为,通过引入材料的刚度矩阵和应变能函数来描述材料的应力应变关系。

根据实际需求,选择合适的本构模型对多晶硅材料进行仿真。

如果需要考虑多晶硅的塑性行为,可以选择弹塑性模型;如果只需考虑多晶硅的弹性行为,可以选择弹性模型。

需要注意的是,弹塑性模型需要输入多个材料参数,如弹性模量、泊松比、屈服应力等;而弹性模型只需要输入弹性模量和泊松比即可。

在使用ABAQUS软件进行多晶硅材料仿真时,除了选择合适的本构模型外,还需要根据实际情况设置加载条件和边界条件。

通过在ABAQUS中定义加载和边界条件,可以模拟多晶硅在不同应力和应变条件下的力学行为。

ABAQUS软件提供了多种多晶硅本构模型供用户选择。

根据实际需求,选择合适的本构模型可以准确描述多晶硅材料的力学行为,并为工程设计和优化提供有力支持。

基于Abaqus子程序的高分子材料本构关系实现

基于Abaqus子程序的高分子材料本构关系实现
mo d e l t o s i mul a t e t h e p o l y me r s .Th e UMAT/ VUMAT s ub r o u t i n e i s pr o v i d e d b y Aba q us t o he l p u s e r t o d e in f e ma t e r i a l c o n s t i t u t i v e mo d e l a c c o r di ng t o t he a c t u a l n e e d s . Th e s t r e s s - s t r a i n r e l a t i o n s h i p o f p o l y me r s l c a n b e i mp l e me n t e d p r e c i s e l y b y de in f e d VUMAT,a n d S O t h e d i f f e r e n c e be t we e n t he s i mu l a t i o n a n d t h e t e s t c a n b e r e d u c e d. Me a n whi l e,t h e d i f f e r e n t u n l o a d pa t h s c a n b e d e in f e d b y di f f e r e n t u n l o a d la f g s . Ke y wo r d s:p o l y me r s;c o ns t i t u t i v e r e l a t i o n;Ab a q us ;UMAT;VUMAT
基于 A b a q u s 子程 序 的高 分 子 材 料本 构 关 系实 现
汪 品红

于ABAQUS有限元建模材料本构分析

于ABAQUS有限元建模材料本构分析

于 ABAQUS 有限元建模材料本构分析摘要:在ABAQUS建模时应选取合理材料本构才能更好的进行分析,本文将对钢筋本构选用模型,混凝土本构选用的损伤模型、受拉损伤因子及受压损伤因子通过规范确定,从而为ABAQUS有限元建模分析提供有力分析。

关键词:钢筋本构;混凝土本构一、钢筋本构ABAQUS中可用Embed将膜或链杆单元嵌入混凝土中,结构自由度由软件自发耦合。

对于钢筋的定义方式,ABAQUS包含多方面的定义。

其中包括定义钢筋的截面积、方向、间距、钢筋对应的单元体的边界编号和其在该边的相对位置。

用户在使用ABAQUS建模时可灵活选用。

常用的钢筋本构有三种。

有双折线模型,双斜线模型及三折线模型三种,考虑计算精度及计算方便,一般选取双斜线模型进行计算。

二、混凝土本构ABAQUS中混凝土本构的模型主要有两种,一是弹塑性断裂(Smearde Crack Model),主用于受压,而受拉用固定弥散裂缝模型来表述。

二是弹塑性断裂和损伤的混凝土模型。

它针对第一种改进了三点。

1、导入损伤参数,折减弹性刚度,以此模拟损伤积聚的过程。

2、导入非关联硬化。

3、手动操控裂缝闭合表现行为,可更真实反映工程实况[2]。

弹塑性损伤模型的原理可以用以下方程来概括。

(1)(2)(3)(4)式1规定了参考损伤效应条件下的有效应力,式2规定了弹性应变和有效应力的数值函数,式3、4则规定了材料的塑性行为。

如混凝土单轴受力,混凝土压、拉时由损伤而起的弹性刚度退化,用Dc与Dt量化阐述。

(5)(6)在ABAQUS中,用以下公式来模拟混凝土受往复荷载作用的损伤指标。

(7)(8)(9)(10)公式中的ωc 和ωt作为参数,ABAQUS中默认俩参数分别为1和0。

在混凝土弹塑性损伤模型中混凝土的弹塑性屈服面如以下公式所表示。

(11)其中:的σc 和σt为承受压、拉时,混凝土材料的有效黏聚应力;σb0和σc0为双、单轴受压的初始屈服应力。

材料参数Kc定值2/3。

abaqus 中的r o 本构

abaqus 中的r o 本构

一、abaqus 中的r o 本构是什么?abaqus 是一种常用的有限元分析软件,可用于模拟和分析复杂的工程结构和材料行为。

在abaqus 中,r o 本构是一种用于描述材料应力-应变行为的数学模型,它通过一些参数来描述材料的变形和强度特性。

r o 本构可以用于模拟各种材料的行为,包括金属、塑料、混凝土等。

二、r o 本构的基本原理r o 本构是基于线弹性理论和塑性理论的,它假设材料在弹性阶段遵循胡克定律,而在塑性阶段则遵循流变规律。

r o 本构通过描述应力和应变之间的关系来表示材料的力学行为,通常采用多种数学公式来拟合不同材料的实际行为。

三、abaqus 中的r o 本构的参数在abaqus 中,r o 本构通常由一些参数来描述,这些参数包括弹性模量、屈服强度、硬化模量、屈服准则等。

这些参数需要根据材料的实际性质进行调整,以精确描述材料的行为。

四、abaqus 中的r o 本构的应用r o 本构在abaqus 中被广泛应用于模拟和分析各种材料的行为。

它可以用于模拟材料在受力过程中的变形、应力分布、破裂行为等各种力学性质,为工程设计和分析提供可靠的依据。

五、abaqus 中的r o 本构的优缺点r o 本构具有较高的通用性和适用性,可以模拟各种材料的行为,并且可以通过调整参数来精确描述材料的力学特性。

然而,r o 本构也存在一些局限性,例如在模拟非线性、温度效应和动态加载等方面的表现可能不够准确。

六、结语abaqus 中的r o 本构是一种重要的数学模型,它可以用于描述材料的应力-应变行为,为工程设计和分析提供可靠的数值模拟结果。

然而,使用r o 本构时需要充分了解材料的实际特性,并且需要进行合理的参数调整,以获得准确可靠的模拟结果。

希望本文可以帮助读者对abaqus 中的r o 本构有更深入的了解。

根据上面的内容,我们可以进一步深入探讨abaqus 中的r o 本构的应用和调整参数的方法,以及其在工程设计和材料分析中的实际意义。

ABAQUS中Cohesive单元建模方法讲解

ABAQUS中Cohesive单元建模方法讲解

复合材料模型建模与分析1。

Cohesive单元建模方法1。

1 几何模型使用内聚力模型(cohesive zone)模拟裂纹的产生和扩展,需要在预计产生裂纹的区域加入cohesive层。

建立cohesive层的方法主要有:方法一、建立完整的结构(如图1(a)所示),然后在上面切割出一个薄层来模拟cohesive 单元,用这种方法建立的cohesive单元与其他单元公用节点,并以此传递力和位移.方法二、分别建立cohesive层和其他结构部件的实体模型,通过“tie”绑定约束,使得cohesive单元两侧的单元位移和应力协调,如图1(b)所示。

(a)cohesive单元与其他单元公用节点(b)独立的网格通过“tie”绑定图1。

建模方法上述两种方法都可以用来模拟复合材料的分层失效,第一种方法划分网格比较复杂;第二种方法赋材料属性简单,划分网格也方便,但是装配及“tie”很繁琐;因此在实际建模中我们应根据实际结构选取较简单的方法。

1.2 材料属性应用cohesive单元模拟复合材料失效,包括两种模型:一种是基于traction-separation描述;另一种是基于连续体描述。

其中基于traction—separation描述的方法应用更加广泛.而在基于traction-separation描述的方法中,最常用的本构模型为图2所示的双线性本构模型。

它给出了材料达到强度极限前的线弹性段和材料达到强度极限后的刚度线性降低软化阶段. 注意图中纵坐标为应力,而横坐标为位移,因此线弹性段的斜率代表的实际是cohesive 单元的刚度.曲线下的面积即为材料断裂时的能量释放率.因此在定义cohesive的力学性能时,实际就是要确定上述本构模型的具体形状:包括刚度、极限强度、以及临界断裂能量释放率,或者最终失效时单元的位移。

常用的定义方法是给定上述参数中的前三项,也就确定了cohesive的本构模型。

Cohesive单元可理解为一种准二维单元,可以将它看作被一个厚度隔开的两个面,这两个面分别和其他实体单元连接。

abaqus中johnson-cook本构模型理解 -回复

abaqus中johnson-cook本构模型理解 -回复

abaqus中johnson-cook本构模型理解-回复Abaqus中Johnson-Cook本构模型理解引言:材料的本构模型是描述材料力学行为的数学方程。

在有限元分析中,本构模型可以用于模拟材料的变形和损伤行为,从而预测材料在不同加载条件下的响应。

Johnson-Cook本构模型是一种常用的本构模型,广泛应用于材料科学和工程领域。

本文将从基本原理开始,逐步解释和理解Abaqus 中Johnson-Cook本构模型。

1. 弹塑性本构模型首先需要了解的是,弹塑性本构模型是最基本的材料模型之一。

它基于线弹性理论,假设材料在小应变范围内具有弹性行为,而在大应变范围内表现出塑性行为。

弹塑性本构模型可以描述材料的应力-应变关系,并预测材料的弹性变形和塑性变形。

2. 材料的温度效应在考虑Johnson-Cook本构模型之前,还需要考虑材料的温度效应。

温度对材料力学行为的影响是复杂而重要的。

温度的增加可以引起材料的软化、蠕变和断裂等现象。

因此,在模拟材料行为时,必须考虑材料的温度效应,并选择适当的本构模型来描述。

3. Johnson-Cook本构模型的基本原理Johnson-Cook本构模型是一种经验模型,用于描述材料的塑性行为和温度效应。

它采用以下形式的应力-应变关系:σ= (A + B ε^n) (1 + C ln(ε˙/ε˙_0))^m (1 - T/T_m)^p其中,σ是材料的应力,ε是应变,ε˙是应变速率,T是材料的温度,A、B、C、n、m、p和T_m是需要通过实验来确定的材料参数。

4. 材料参数的确定为了使用Johnson-Cook本构模型,需要通过实验来确定材料参数。

这些参数通常由材料的拉伸实验和冲击实验等得到。

拉伸实验可以提供材料的应力-应变曲线,以及材料的屈服强度和断裂应变等信息。

冲击实验可以提供材料的应变率敏感性和断裂韧性等信息。

根据实验数据,可以使用不同的方法来确定Johnson-Cook本构模型的参数。

ABAQUS-材料本构模型及编程

ABAQUS-材料本构模型及编程

材料本构模型及编程-ABAQUS-UMAT材料本构模型及编程实现:简介1、什么时候用用户定义材料(User-defined material, UMAT)?很简单,当ABAQUS没有提供我们需要的材料模型时。

所以,在决定自己定义一种新的材料模型之前,最好对ABAQUS已经提供的模型心中有数,并且尽量使用现有的模型,因为这些模型已经经过详细的验证,并被广泛接受。

2、好学吗?需要哪些基础知识?先看一下ABAQUS手册(ABAQUS Analysis User's Manual)里的一段话:Warning: The use of this option generally requires considerable expertise. The user is cautioned that the implementation of any realis tic constitutive model requires extensive development and testing. Initial testing on a single element model with prescribed traction l oading is strongly recommended.但这并不意味着非力学专业,或者力学基础知识不很丰富者就只能望洋兴叹,因为我们的任务不是开发一套完整的有限元软件,而只是提供一个描述材料力学性能的本构方程(Constitutive equation)而已。

当然,最基本的一些概念和知识还是要具备的,比如应力(stress),应变(strain)及其分量; volumetric part和deviatoric part;模量(modulus)、泊松比(Poisson’s ratio)、拉美常数(Lame constant);矩阵的加减乘除甚至求逆;还有一些高等数学知识如积分、微分等。

基于abaqus的木材本构关系数值模拟方法

基于abaqus的木材本构关系数值模拟方法

基于abaqus的木材本构关系数值模拟方法1.引言1.1 概述概述部分的内容示例:在工程领域中,对木材的力学性能和行为的研究一直存在着广泛的关注。

木材作为一种常见的工程材料,其强度、刚度和耐久性对于工程结构的设计和安全具有重要的影响。

研究木材本构关系是了解木材力学行为及其材料特性的基础,对于优化木材结构的设计和评估具有重要意义。

随着计算机技术的进步和有限元分析方法的广泛应用,通过数值模拟方法来研究木材本构关系变得越来越普遍和有效。

其中,ABAQUS作为一款常用的有限元分析软件,具有强大的建模和分析能力,被广泛应用于各个领域的工程问题求解中。

基于ABAQUS的数值模拟方法可以通过建立木材的有限元模型,模拟木材受力过程中的变形、应力分布和破坏行为。

本文旨在介绍基于ABAQUS的数值模拟方法,用于研究木材的本构关系。

首先,将会概述木材本构关系的研究现状,包括已有的实验数据和模型。

然后,介绍基于ABAQUS的数值模拟方法,包括材料参数的输入、模型的建立和加载条件的设置。

通过数值模拟可以得到木材在各种力学载荷下的力学性能,进而分析其行为特性。

最后,通过结果分析,我们将对该数值模拟方法的有效性和实用性进行评价,并指出研究的局限性以及未来的发展方向。

本文的研究旨在提供一种基于ABAQUS的数值模拟方法,为工程师和研究人员提供一个有效的工具,用于理解和预测木材的力学行为和性能,以支持木材结构的合理设计和优化。

希望本文的内容能够对读者在木材本构关系的研究和应用方面提供一定的指导和参考。

文章结构部分的内容可以按照以下方式进行编写:本篇文章的结构如下:一、引言1.1 概述1.2 文章结构1.3 目的二、正文2.1 木材本构关系的研究现状2.2 基于abaqus的数值模拟方法三、结论3.1 结果分析3.2 研究的局限性和展望在引言部分,我们将对本文的研究背景、意义和目的进行介绍。

首先,我们会提出木材本构关系的研究现状,包括目前已有的研究成果以及存在的问题。

(完整版)ABAQUS中的三种混凝土本构模型

(完整版)ABAQUS中的三种混凝土本构模型

ABQUS中的三种混凝土本构模型ABAQUS 用连续介质的方法建立描述混凝土模型不采用宏观离散裂纹的方法描述裂纹的水平的在每一个积分点上单独计算其中。

低压力混凝土的本构关系包括:Concrete Smeared cracking model (ABAQUS/Standard)Concrete Brittle cracking model (ABAQUS/Explicit)Concrete Damage plasticity model高压力混凝土的本构关系:Cap model1、ABAQUS/Standard中的弥散裂缝模型Concrete Smeared cracking model (ABAQUS/Standard):——只能用于ABAQUS/Standard中裂纹是影响材料行为的最关键因素,它将导致开裂以及开裂后的材料的各向异性用于描述:单调应变、在材料中表现出拉伸裂纹或者压缩时破碎的行为在进行参数定义式的Keywords:*CONCRETE*TENSION STIFFENING*SHEAR RETENTION*FAILURE RATIOS2、ABAQUS/Explicit中脆性破裂模型Concrete Brittle cracking model (ABAQUS/Explicit) :适用于拉伸裂纹控制材料行为的应用或压缩失效不重要,此模型考虑了由于裂纹引起的材料各向异性性质,材料压缩的行为假定为线弹性,脆性断裂准则可以使得材料在拉伸应力过大时失效。

在进行参数定义式的Keywords*BRITTLE CRACKING,*BRITTLE FAILURE,*BRITTLE SHEAR3、塑性损伤模型Concrete Damage plasticity model:适用于混凝土的各种荷载分析,单调应变,循环荷载,动力载荷,包含拉伸开裂(cracking)和压缩破碎(crushing),此模型可以模拟硬度退化机制以及反向加载刚度恢复的混凝土力学特性在进行参数定义式的Keywords:*CONCRETE DAMAGED PLASTICITY*CONCRETE TENSION STIFFENING*CONCRETE COMPRESSION HARDENING*CONCRETE TENSION DAMAGE*CONCRETE COMPRESSION DAMAGE。

abaqus应变软化本构模型

abaqus应变软化本构模型

abaqus应变软化本构模型abaqus应变软化本构模型是一种常用的材料本构模型,用于描述材料在加载过程中的应变软化行为。

本文将介绍abaqus应变软化本构模型的原理和应用。

应变软化是指材料在受到加载时,随着应变的增加,材料的刚度逐渐降低的现象。

这种现象在很多材料中都存在,特别是在一些脆性材料中,如混凝土、岩石等。

在模拟这些材料的力学行为时,使用abaqus应变软化本构模型可以更准确地描述材料的实际行为。

abaqus应变软化本构模型是基于塑性力学理论的,它假设材料的本构关系在弹性阶段和塑性阶段是不同的。

在弹性阶段,材料遵循胡克定律,即应力与应变成线性关系。

而在塑性阶段,材料的刚度会随着应变的增加而降低,这种现象可以通过引入软化函数来描述。

在abaqus中,常用的应变软化本构模型包括弹塑性本构模型和本构模型。

弹塑性本构模型适用于强度较高的材料,如钢材。

而本构模型适用于较脆性的材料,如混凝土、岩石等。

在abaqus中,应变软化本构模型的参数可以通过试验数据进行确定。

常用的试验包括压缩试验、拉伸试验和剪切试验等。

通过对试验数据的拟合,可以得到材料的本构参数,进而进行数值模拟。

应变软化本构模型在工程实践中有着广泛的应用。

例如,在土木工程中,模拟混凝土的破坏过程是一项重要的任务。

混凝土在受到加载时会发生应变软化现象,使用abaqus应变软化本构模型可以更准确地模拟混凝土的破坏过程,为工程设计提供可靠的依据。

除了土木工程,应变软化本构模型还可以应用于岩石力学、金属材料力学等领域。

在岩石力学中,岩石在受到加载时会发生应变软化现象,使用abaqus应变软化本构模型可以更好地模拟岩石的破坏行为,为岩石工程提供可靠的分析结果。

在金属材料力学中,金属材料的应变软化行为对于模拟金属的变形和破坏过程至关重要,使用abaqus应变软化本构模型可以更准确地描述金属材料的力学行为。

abaqus应变软化本构模型是一种常用的材料本构模型,可以很好地描述材料在加载过程中的应变软化行为。

umat

umat

材料本构模型及编程实现:简介一起学习UMATZHANG chunyu1、什么时候用用户定义材料(User-defined material, UMAT)?很简单,当ABAQUS没有提供我们需要的材料模型时。

所以,在决定自己定义一种新的材料模型之前,最好对ABAQUS已经提供的模型心中有数,并且尽量使用现有的模型,因为这些模型已经经过详细的验证,并被广泛接受。

2、好学吗?需要哪些基础知识?先看一下ABAQUS手册(ABAQUS Analysis User's Manual)里的一段话:Warning: The use of this option generally requires consider able expertise. The user is cautioned that the implementationof any realistic constitutive model requires extensive devel opment and testing. Initial testing on a single element model with prescribed traction loading is strongly recommended.但这并不意味着非力学专业,或者力学基础知识不很丰富者就只能望洋兴叹,因为我们的任务不是开发一套完整的有限元软件,而只是提供一个描述材料力学性能的本构方程(Constitutive equation)而已。

当然,最基本的一些概念和知识还是要具备的,比如应力(stress),应变(strain)及其分量;volumetric part和deviatoric part;模量(modulus)、泊松比(Poisson’s ratio)、拉美常数(Lame constant);矩阵的加减乘除甚至求逆;还有一些高等数学知识如积分、微分等。

ABAQUS钢筋混凝土本构模型

ABAQUS钢筋混凝土本构模型

ABAQUS钢筋混凝土本构模型钢是各向同性材料,其本构关系理论比较成熟,考虑了其弹性、弹塑性、强化、断裂和包辛格效应并得到充分验证。

基本参数:密度:ρ=7800kg/m^3弹性模量:E_s=2.07×10^5泊松比:ν =0.31.2 混凝土混凝土在拉压方向上的力学性能不同,存在着强化、软化、开裂、损伤等复杂的力学行为。

如何在有限元程序中准确模拟混凝土的本构关系,对于后续有限元计算结构的合理性和准确性尤为重要。

基本参数:密度:ρ=2200~2400kg/m^3弹性模量:E_c(与强度有关)泊松比:ν =0.18~0.22(建议取0.2)Ψe fb0/fc Kυ30°0.11.160.6677.5e-042 混凝土单轴应力-应变关系2.1 混凝土单轴受压应力-应变关系混凝土材料在单轴压缩下的应力-应变关系由弹性段、强化段和软化段组成。

图1 混凝土单轴应力-应变关系ε_{c0}^{el}——未损伤或者未考虑损伤的混凝土受压弹性应变,材料无损时的弹性应变ε_c^{el}——考虑损伤的混凝土受压弹性应变(损伤导致刚度减小,相应的弹性应变就增大了)ε_c^{pl}——混凝土受压塑性应变(总应变减去考虑损伤的受压弹性应变)ε_c^{in}——混凝土受压非弹性应变(包括了一部分塑性应变和受压损伤导致的刚度变小产生的应变等)1.弹性段定义——确定初始切线模量E0(1) 确定弹性极限点(ε_{c,e0},σ_{c,e0}) \\建议一般取σ_{c,e0}=f_c/3 \\则初始切线弹性模量为E_0=ε_{c,e0}/σ_{c,e0} \\(2) 混凝土的弹性模量Ec(3) 也可以采用如下方法进行确定:首先计算混凝土拉伸开裂时的割线模量,并按此割线模量取值确定混凝土压缩应力-应变关系曲线上升段中割线模量的等值点,以此作为混凝土受压受力阶段的弹塑性分界点,通过这样的方法可以保证混凝土的压缩弹性模量和拉伸弹性模量取值保持一致。

abaqus的多晶硅本构模型

abaqus的多晶硅本构模型

abaqus的多晶硅本构模型多晶硅是一种常用的材料,具有广泛的应用领域,如电子器件、太阳能电池等。

为了准确描述多晶硅的力学行为,在abaqus中使用了多晶硅本构模型。

多晶硅本构模型是建立在晶体塑性理论基础上的。

晶体塑性理论认为,材料的塑性行为是由晶体中的位错运动引起的。

位错是晶体中的一种缺陷,它可以沿晶体晶格滑移面移动,导致材料发生塑性变形。

多晶硅的本构模型考虑了这种位错滑移的影响。

多晶硅的本构模型可以分为两个部分:晶体塑性本构模型和晶粒间转变本构模型。

晶体塑性本构模型描述了单个晶体的力学行为,晶粒间转变本构模型描述了晶粒之间的相互作用。

在abaqus中,晶体塑性本构模型通常使用von Mises屈服准则来描述晶体的塑性行为。

von Mises屈服准则假设材料的塑性变形发生在位错滑移面上,并且滑移面上的应力满足一定的准则。

根据这个准则,abaqus可以计算出材料的应力和应变。

晶粒间转变本构模型描述了晶粒之间的相互作用。

多晶硅中存在大量的晶粒,它们之间的相互作用会影响材料的力学行为。

在abaqus 中,晶粒间转变本构模型通常采用屈服面法来描述晶粒间的相互作用。

屈服面法假设晶粒的屈服行为可以用一个屈服面来描述,当晶体的应力达到这个屈服面时,晶粒会发生转变。

使用abaqus的多晶硅本构模型可以准确地模拟多晶硅的力学行为。

通过对多晶硅材料进行拉伸、压缩等加载条件的仿真分析,可以得到材料的应力、应变等力学特性。

这些数据对于材料的设计和优化非常重要。

abaqus的多晶硅本构模型是一种有效的工具,可以帮助工程师和科研人员更好地理解和预测多晶硅材料的力学行为。

通过对多晶硅的力学行为进行仿真分析,可以为多晶硅相关领域的研究和应用提供有力的支持。

在ABAQUS中开发实现Duncan_Chang本构模型

在ABAQUS中开发实现Duncan_Chang本构模型
降低的影响。
1034
岩土力学
2004 年
对于卸载情况的处理方法,不同文献所考虑的 准则有所不同,但卸载弹性模量计算公式应统一表 达为
Eur
= kur Pa
⎜⎛ ⎝
−σ 1 Pa
⎟⎞ nur ⎠
(8)
至 此 , 由 式 (2) ~ (8) 可 以 全 部 确 定 DuncanChang E-B 模型刚度系数矩阵的元素,然后,可按 照三维问题或平面应变问题的数学表达式形成相应 的刚度系数矩阵 [D] ,如对于 Duncan-Chang E-B 三 维模型问题有
修改稿收到日期:2003-09-04
作者简介:徐远杰,男,1956 年生,教授,结构破坏与数值仿真研究所所长,从事岩土工程数值仿真研究。
第7期
徐远杰等:在 ABAQUS 中开发实现 Duncan-Chang 本构模型
1033
一不足,本文作者利用 ABAQUS 提供的二次开发 用户子程序,完成了 Duncan-Chang 本构模型[1]的接 口开发工作,并且对若干模型算例进行了测试,计 算结果表明,模型算例的求解效率和精度令人满意, 有望成为强大的解决大型复杂土工问题的非线性数 值分析工具。同时,开发与维护的工作量将大大降 低,是解决三维复杂土工数值分析问题的快捷方案 和有效途径。
XU Yuan-jie, WANG Guan-qi, LI Jian, TANG Bi-hua
(School of Civil and Architectural Engineering, Wuhan University, Wuhan 430072, China)
Abstract: The software of ABAQUS, which is developed and supported by HKS Inc., is currently one of the most powerful nonlinear finite element analysis implement in the world. It includes many different kinds of material models. However, the release of ABAQUS currently absents Duncan-Chang constitutive model that has been widely used to simulate some deformation problems in the field of geotechnical engineering. The shortage limits the software’s application scopes in some extents. In this paper, the procedure and outline to develop Duncan-Chang constitutive model in ABAQUS by using User Subroutine were proposed and implemented. Two numerical examples were also solved and compared with conventional triaxial compression test data. All numerical results exhibit that it will be possible to accomplish some complex stress-strain analysis for soil engineering problems by fully making use of the powerful solving platform after ABAQUS is added Duncan-Chang constitutive model. Meanwhile, there are some advantages of faster convergence, higher efficiency, higher precision and more conveniently preprocessing-postprocessing. The programming difficulties to develop soil engineering and workloads of maintenance are also greatly decreased. Key words: finite element; Duncan-Chang constitutive model; ABAQUS; redevelopment

abaqus的多晶硅本构模型

abaqus的多晶硅本构模型

abaqus的多晶硅本构模型abaqus是一种常用的有限元分析软件,可用于模拟材料的力学行为。

本文将介绍abaqus中的多晶硅本构模型,该模型可以用来描述多晶硅在不同载荷下的力学响应。

多晶硅是一种常见的半导体材料,由多个晶粒组成。

在力学加载下,晶粒之间会发生相对位移和变形,从而影响整体材料的力学性能。

为了准确地模拟多晶硅的力学行为,需要使用合适的本构模型来描述晶粒间的相互作用。

多晶硅的本构模型通常基于弹塑性理论,考虑了晶粒的变形、滑移和断裂等现象。

其中一个常用的本构模型是晶体塑性有限元法(CPFE),它将晶粒视为一个个离散的小单元,并考虑了晶粒内部的位错运动和晶界的滑移。

在abaqus中,可以通过定义材料的本构行为来使用多晶硅本构模型。

首先,需要指定多晶硅的晶体学方向和晶粒的排列方式。

然后,可以选择合适的本构模型来描述晶粒的力学行为。

常用的多晶硅本构模型包括晶体塑性有限元法、晶体塑性本构模型和粘弹性本构模型等。

晶体塑性有限元法将晶粒视为一个个小单元,通过求解位错运动方程和晶界滑移方程来描述晶粒的变形行为。

晶体塑性本构模型则将晶粒内部的变形行为表示为一个连续的本构关系,常用的模型有von Mises模型和Hill模型等。

粘弹性本构模型则考虑了晶粒间的粘滞行为,将晶粒视为具有弹性和粘性特性的材料。

在abaqus中,可以通过定义材料的本构行为来选择合适的多晶硅本构模型。

可以根据实际情况选择合适的模型参数,并进行力学分析。

通过对多晶硅材料进行模拟和分析,可以得到材料的应力、应变分布以及其他力学性能指标,为工程设计和材料优化提供重要参考。

总结起来,abaqus的多晶硅本构模型可以用来模拟多晶硅在不同载荷下的力学行为。

通过选择合适的本构模型和参数,可以准确地描述多晶硅的变形和破坏行为。

这对于研究多晶硅材料的力学性能、优化材料设计以及预测材料的寿命和可靠性具有重要意义。

abaqus材料子程序

abaqus材料子程序

各向同性材料损伤本构模型SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,+ RPL,DDSDDT,DRPLDE,DRPLDT,+ STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,CMNAME, + NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,PNEWDT, + CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)INCLUDE 'ABA_PARAM.INC'CHARACTER*80 CMNAMEDIMENSION STRESS(NTENS),STATEV(NSTATV),+ DDSDDE(NTENS,NTENS),DDSDDT(NTENS),+ DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),+ TIME(2),PREDEF(1),DPRED(1),PROPS(NPROPS),+ COORDS(3),DROT(3,3),DFGRD0(3,3),DFGRD1(3,3)DIMENSION STRANT(6),TSTRANT(4),PT(1)DIMENSION OLD_STRESS(6)DIMENSION DOLD_STRESS(6),D_STRESS(6)DIMENSION C(6,6),CD(6,6),DSTRESS(6),BSTRESS(6),ROOT(3),+ DFMNDE(6),DDMDE(6),DCDDM(6,6),ATEMP1(6), ATEMP2(6)PARAMETER (ZERO=0.D0,ONE=1.D0,TWO=2.D0,FOUR=4.D0,HALF = 0.5D0) C startC IF (NPROPS.LT.2) THENC WRITE(7,*) '** ERROR: UMAT REQUIRES *NPROPS=2'C STOPC END IFE11 =PROPS(1)V12 =PROPS(2)G12 =PROPS(1)/TWO/(ONE+PROPS(2))C Critical values of streSSEsXT=PROPS(3)XC=PROPS(4)XS=PROPS(5)GX=PROPS(6) !Fracture energy in matrixETA=0.001C Current strainDO I = 1, NTENSSTRANT(I) = STRAN(I) + DSTRAN(I) END DOC StiffnessDO I = 1, 6DO J = 1, 6C(I,J)=ZEROEND DOEND DOATEMP = (1+V12)*(1-TWO*V12)C(1,1) = E11*(1-V12)/ATEMPC(2,2) = E11*(1-V12)/ATEMPC(3,3) = E11*(1-V12)/ATEMPC(1,2) = E11*V12/ATEMPC(1,3) = E11*V12/ATEMPC(2,3) = E11*V12/ATEMPC(4,4) = G12C(5,5) = G12C(6,6) = G12DO I = 2, 6DO J = 1, I-1C(I,J) = C(J,I)END DOEND DOC Critical values of strainsXET=XT/(C(1,1)-2*V12*C(1,2))XEC=XC/(C(1,1)-2*V12*C(1,2))XES=XS/C(4,4)DMOLD = STATEV(1)C Strain initiation criterionA11 = STRANT(1)**TWO+STRANT(2)**TWO+STRANT(3)**TWO A12 = A11 / XET / XECA21 = STRANT(1)+STRANT(2)+STRANT(3)A22 = (XEC - XET) / XEC / XET * A21A31 = STRANT(4)**TWO+STRANT(5)**TWO+STRANT(6)**TWO A32 = A31 / XES**TWOA1= A12 + A22 + A32C B11 = STRANT(2)**TWOC B12 = B11 / XET / XECC B21 = STRANT(2)C B22 = (XEC - XET) / XEC / XET * B21C B31 = STRANT(5)**TWOC B32 = B31 / XES**TWOC B1= B12 + B22 + B32C C11 = STRANT(3)**TWOC C12 = C11 / XET / XECC C21 = STRANT(3)C C22 = (XEC - XET) / XEC / XET * C21C C31 = STRANT(6)**TWOC C32 = C31 / XES**TWOC C1= C12 + C22 + C32STATEV(2)=A1C STATEV(3)=B1C STATEV(4)=C1FMN = ZEROIF (A1.GT.ZERO) THEN FMN =SQRT(A1)C IF (B1.GT.ONE) THEN C FMN =FMN+SQRT(B1) C IF(C1.GT.ONE) THEN C FMN =FMN+SQRT(C1)C END IFC END IFEND IFSTATEV(5)=FMNC write(*,*) FMNDM = ZERODDMDFMN = ZERODO I = 1, 6DFMNDE(I) = ZERODDMDE(I) = ZEROEND DOIF (FMN .GT. ONE) THENC CALCULATE DM, DDMDFMNC WRITE(6,*)FMNT1 = (C(1,1)-2*V12*C(1,2)) * XET**2 * CELENT / GXT2 = (ONE - FMN) * T1DM = ONE - EXP(T2)/FMNC WRITE(6,*)'T1 ',T1,' T2', T2, ' DM', DMC write(*,*) DMC CALCULATE THE DERIVATIVE OF DAMAGE VARIABLE WITH RESPECT TO FAILURE C RITERIONDDMDFMN = (ONE / FMN + T1) * (ONE - DM)C CALCULATE DFMNDEIF (DM .GT. DMOLD) THENDFMNDE(1) = HALF/FMN*(TWO*STRANT(1)+XEC-XET)/XET/XEC DFMNDE(2) = HALF/FMN*(TWO*STRANT(2)+XEC-XET)/XET/XEC DFMNDE(3) = HALF/FMN*(TWO*STRANT(3)+XEC-XET)/XET/XEC DFMNDE(4) = ONE/FMN*TWO*STRANT(4)/XES**TWODFMNDE(5) = ONE/FMN*TWO*STRANT(5)/XES**TWODFMNDE(6) = ONE/FMN*TWO*STRANT(6)/XES**TWODO I = 1, 6DDMDE(I) = DFMNDE(I) * DDMDFMNEND DOEND IFEND IFDM = MAX (DM, DMOLD)C write(6,*) DMC SAVE THE OLD STRESS TO OLD_STRESSDO I = 1, NTENSOLD_STRESS(I) = STRESS(I) END DOC Effective stiffnessDO I = 1, 6DO J = 1, 6CD(I,J)=C(I,J)END DOEND DOIF (DM.NE.ZERO) THENCD(1,1) = (ONE - DM)*C(1,1) CD(1,2) = (ONE - DM)*C(1,2) CD(2,1) = CD(1,2)CD(2,2) = (ONE - DM)*C(2,2) CD(1,3) = (ONE - DM)*C(1,3) CD(3,1) = CD(1,3)CD(2,3) = (ONE- DM)*C(2,3)CD(3,2) = CD(2,3)CD(4,4) = (ONE - DM)*C(4,4) CD(5,5) = (ONE - DM)*C(5,5) CD(6,6) = (ONE - DM)*C(5,5) END IFC Elastic derivativeDO I = 1, 6DO J = 1, 6DCDDM(I,J) = ZEROEND DOEND DOCC CALCULATE DC/DDMCDCDDM(1,1) = -C(1,1) DCDDM(1,2) = -C(1,2) DCDDM(2,1) = -C(2,1)DCDDM(2,2) = -C(2,2)DCDDM(2,3) = -C(2,3)DCDDM(3,2) = -C(3,2)DCDDM(3,3) = -C(3,3)DCDDM(3,1) = -C(3,1)DCDDM(1,3) = -C(1,3)DCDDM(4,4) = -C(4,4)DCDDM(5,5) = -C(5,5)DCDDM(6,6) = -C(6,6)C UPDATE THE JACOBIANDO I = 1, NTENSATEMP1(I) = ZERODO J = 1, NTENSATEMP1(I) = ATEMP1(I) + DCDDM(I,J) * STRANT(J) END DOEND DODDSDDE=0DO I = 1, NTENSDO J = 1, NTENSDDSDDE(I,J)=CD(I,J)+(ATEMP1(I)*DDMDE(J))*DTIME/(DTIME+ETA) END DO END DOC Update streSSEsDO I = 1, NTENSSTRESS(I)=ZERODO J = 1, NTENSC IF(DM.LT.0.5) THENSTRESS(I)=STRESS(I)+ CD(I,J) * STRANT(J)C ELSEC STRESS(I)=STRESS(I)+ CD(I,J) * STRANT(J) * (1-DM)C ENDIFEND DOEND DOC EnergyDO I = 1, NDISSE = SSE + HALF * (STRESS(I) + OLD_STRESS(I)) * DSTRAN(I) END DO DO I = NDI+1, NTENSSSE = SSE + (STRESS(I) + OLD_STRESS(I)) * DSTRAN(I)END DOSTATEV(1)=DMRETURNEND。

Abaqus混凝土材料模型解读与参数设置

Abaqus混凝土材料模型解读与参数设置

Abaqus混凝土材料塑性损伤模型浅析与参数设置【壹讲壹插件】欢迎转载,作者:星辰-北极星,QQ群:431603427Abaqus混凝土材料塑性损伤模型浅析与参数设置 (1)第一部分:Abaqus自带混凝土材料的塑性损伤模型 (2)1.1概要 (2)1.2学习笔记 (2)1.3 参数定义与说明 (3)1.3.1材料模型选择:Concrete Damaged Plasticity (3)1.3.2 混凝土塑性参数定义 (3)1.3.3 混凝土损伤参数定义: (4)1.3.4 损伤参数定义与输出损伤之间的关系 (4)1.3.5 输出参数: (4)第二部分:根据GB50010-2010定义材料损伤值 (5)第三部分:星辰-北极星插件介绍:POLARIS-CONCRETE (6)3.1 概要 (6)3.2 插件的主要功能 (6)3.3 插件使用方法: (6)3.3.1 插件界面: (6)3.3.2 生成结果 (7)3.4、算例: (9)3.4.1三维实体简支梁模型说明 (9)3.4.2 计算结果: (9)第一部分:Abaqus自带混凝土材料的塑性损伤模型1.1概要首先我要了解Abaqus内自带的参数模型是怎样的,了解其塑性模型,进而了解其损伤模型,其帮助文档Abaqus Theory Manual 4.5.1 An inelastic constitutive model for concrete讲述的是其非弹性本构,4.5.2 Damaged plasticity model for concrete and other quasi-brittle materials则讲述的塑性损伤模型,同时在Abaqus Analysis User's Manual 22.6 Concrete也讲述了相应的内容。

1.2学习笔记1、混凝土塑性损伤本构模型中的损伤是一标量值,数值范围为(0无损伤~1完全失效[对于混凝土塑性损伤一般不存在]);2、仅适用于脆性材料在中等围压条件(为围压小于轴抗压强度1/4);3、拉压强度可设置成不同数值;4、可实现交变载荷下的刚度恢复;默认条件下,由拉转压刚度恢复,由压转拉刚度不变;5、强度与应变率相关;6、使用的是非相关联流动法则,刚度矩阵为非对称,因此在隐式分析步设置时,需在分析定义other-》Matrix storate-》Unsymmetric。

abaqus导入材料本构方程

abaqus导入材料本构方程

abaqus导入材料本构方程Abaqus导入材料本构方程材料本构方程是描述材料力学性能的数学模型,它可以用来预测材料在不同应力和应变条件下的力学响应。

Abaqus是一种常见的有限元分析软件,它可以用来模拟和分析各种材料的力学行为。

在Abaqus中,导入材料本构方程可以帮助工程师更准确地模拟材料的行为,从而更好地预测结构的性能。

在Abaqus中,导入材料本构方程的步骤如下:1. 定义材料类型:在Abaqus中,需要首先定义材料的类型。

常见的材料类型包括线性弹性材料、非线性弹性材料和塑性材料等。

根据材料的力学性质,选择合适的材料类型。

2. 定义材料参数:根据材料本构方程的形式和参数,输入材料的力学性质参数。

这些参数可以是材料的弹性模量、泊松比、屈服强度等。

根据具体材料的测试数据或者文献资料,确定合适的参数值。

3. 编写材料本构方程:根据材料的力学行为和材料本构方程的形式,编写表示材料本构方程的数学表达式。

在Abaqus中,可以使用用户子程序(User Subroutine)来定义材料本构方程。

用户子程序是一种自定义的程序,可以将用户自己编写的材料本构方程导入Abaqus 中。

4. 导入材料本构方程:将编写好的用户子程序导入Abaqus中,并指定该材料本构方程对应的材料类型和材料参数。

通过这一步骤,Abaqus就可以根据用户定义的材料本构方程来模拟和分析材料的力学行为了。

导入材料本构方程可以提高模拟和分析的准确性,使得仿真结果更符合实际情况。

通过使用合适的材料本构方程,可以考虑材料的非线性行为、塑性变形以及破坏等情况,从而更全面地评估结构的性能。

例如,在模拟金属材料的变形行为时,常常需要考虑材料的塑性本质和应力应变曲线的非线性特性。

通过导入合适的材料本构方程,可以更准确地预测材料的变形和破坏过程,从而指导结构设计和优化。

Abaqus导入材料本构方程是一种重要的工具,可以帮助工程师更好地模拟和分析材料的力学行为。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

材料本构模型及编程-ABAQUS-UMAT材料本构模型及编程实现:简介1、什么时候用用户定义材料(User-defined material,UMAT)?很简单,当ABAQUS没有提供我们需要的材料模型时。

所以,在决定自己定义一种新的材料模型之前,最好对ABAQUS已经提供的模型心中有数,并且尽量使用现有的模型,因为这些模型已经经过详细的验证,并被广泛接受。

2、好学吗?需要哪些基础知识?先看一下ABAQUS手册(ABAQUS Analysis User's Manual)里的一段话:Warning:The use of this option generally requires considerable expertise.The user is cautioned that the implementation of any realistic constit utive model requires extensive development and testing.Initial testing on a single element model with prescribed traction loading is strongly rec ommended.但这并不意味着非力学专业,或者力学基础知识不很丰富者就只能望洋兴叹,因为我们的任务不是开发一套完整的有限元软件,而只是提供一个描述材料力学性能的本构方程(Constitutive equation)而已。

当然,最基本的一些概念和知识还是要具备的,比如应力(stress),应变(strain)及其分量;volumetric part和deviatoric part;模量(modulus)、泊松比(Poisson’s ratio)、拉美常数(Lame constant);矩阵的加减乘除甚至求逆;还有一些高等数学知识如积分、微分等。

3、UMAT的基本任务?我们知道,有限元计算(增量方法)的基本问题是:已知第n步的结果(应力,应变等),;然后给出一个应变增量,计算新的应力。

UMAT要完成这一计算,并要计算Jacobian矩阵DDSDDE(I,J)=。

是应力增量矩阵(张量或许更合适),是应变增量矩阵。

DDSDDE(I,J)定义了第J个应变分量的微小变化对第I个应力分量带来的变化。

该矩阵只影响收敛速度,不影响计算结果的准确性(当然,不收敛自然得不到结果)。

4、怎样建立自己的材料模型?本构方程就是描述材料应力应变(增量)关系的数学公式,不是凭空想象出来的,而是根据实验结果作出的合理归纳。

比如对弹性材料,实验发现应力和应变同步线性增长,所以用一个简单的数学公式描述。

为了解释弹塑性材料的实验现象,又提出了一些弹塑性模型,并用数学公式表示出来。

对各向同性材料(Isotropic material),经常采用的办法是先研究材料单向应力-应变规律(如单向拉伸、压缩试验),并用一数学公式加以描述,然后把讲该规律推广到各应力分量。

这叫做“泛化“(generalization)。

5、一个完整的例子及解释下面这个UMAT取自ABAQUS手册,是一个用于大变形下的弹塑性材料模型。

希望我的注释能帮助初学者理解。

需要了解J2理论。

SUBROUTINE UMAT(STRESS,STATEV,DDSDDE,SSE,SPD,SCD,RPL,DDSDDT,1DRPLDE,DRPLDT,STRAN,DSTRAN,TIME,DTIME,TEMP,DTEMP,PREDEF,DPRED,2CMNAME,NDI,NSHR,NTENS,NSTATV,PROPS,NPROPS,COORDS,DROT,3PNEWDT,CELENT,DFGRD0,DFGRD1,NOEL,NPT,LAYER,KSPT,KSTEP,KINC)STRESS--应力矩阵,在增量步的开始,保存并作为已知量传入UMAT;在增量步的结束应该保存更新的应力;STRAN--当前应变,已知。

DSTRAN—应变增量,已知。

STATEV--状态变量矩阵,用来保存用户自己定义的一些变量,如累计塑性应变,粘弹性应变等等。

增量步开始时作为已知量传入,增量步结束应该更新;DDSDDE=。

需要更新DTIME—时间增量dt。

已知。

NDI—正应力、应变个数,对三维问题、轴对称问题自然是3(11,22,33),平面问题是2(11,22);已知。

NSHR—剪应力、应变个数,三维问题时3(12,13,23),轴对称问题是1(12);已知。

NTENS=NTENS NSHR,已知。

PROPS材料常数矩阵,如模量啊,粘度系数啊等等;作为已知量传入,已知。

DROT—对finite strain问题,应变应该排除旋转部分,该矩阵提供了旋转矩阵,详见下面的解释。

已知。

PNEWDT—可用来控制时间步的变化。

如果设置为小于1的数,则程序放弃当前计算,并用新的时间增量DTIME X PNEWDT作为新的时间增量计算;这对时间相关的材料如聚合物等有用;如果设为大余1的数,则下一个增量步加大DTIME为DTIME X PNEWDT。

可以更新。

其他变量含义可参看手册,暂时用不到。

CINCLUDE'ABA_PARAM.INC'定义了一些参数,变量什么的,不用管CCHARACTER*8CMNAMECDIMENSION STRESS(NTENS),STATEV(NSTATV),DDSDDE(NTENS,NTENS),1DDSDDT(NTENS),DRPLDE(NTENS),STRAN(NTENS),DSTRAN(NTENS),2PREDEF(1),DPRED(1),PROPS(NPROPS),COORDS(3),DROT(3,3),3DFGRD0(3,3),DFGRD1(3,3)矩阵的尺寸声明CC LOCAL ARRAYSC----------------------------------------------------------------C EELAS-ELASTIC STRAINSC EPLAS-PLASTIC STRAINSC FLOW-DIRECTION OF PLASTIC FLOWC----------------------------------------------------------------C局部变量,用来暂时保存弹性应变、塑性应变分量以及流动方向DIMENSION EELAS(6),EPLAS(6),FLOW(6)CPARAMETER(ZERO=0.D0,ONE=1.D0,TWO=2.D0,THREE=3.D0,SIX=6.D0,1ENUMAX=.4999D0,NEWTON=10,TOLER=1.0D-6)CC----------------------------------------------------------------C UMAT FOR ISOTROPIC ELASTICITY AND ISOTROPIC MISES PLASTICITYC CANNOT BE USED FOR PLANE STRESSC----------------------------------------------------------------C PROPS(1)-EC PROPS(2)-NUC PROPS(3..)-SYIELD AN HARDENING DATAC CALLS HARDSUB FOR CURVE OF YIELD STRESS VS.PLASTIC STRAINC----------------------------------------------------------------CC ELASTIC PROPERTIESC获取杨氏模量,泊松比,作为已知量由PROPS向量传入EMOD=PROPS(1)EENU=PROPS(2)νEBULK3=EMOD/(ONE-TWO*ENU)3KEG2=EMOD/(ONE ENU)2GEG=EG2/TWO GEG3=THREE*EG3GELAM=(EBULK3-EG2)/THREEλDO K1=1,NTENSDO K2=1,NTENSDDSDDE(K1,K2)=ZEROEND DOEND DO弹性部分,Jacobian矩阵很容易计算注意,在ABAQUS中,剪切应变采用工程剪切应变的定义,所以剪切部分模量是G而不是2G!CC ELASTIC STIFFNESSCDO K1=1,NDIDO K2=1,NDIDDSDDE(K2,K1)=ELAMEND DODDSDDE(K1,K1)=EG2 ELAMEND DODO K1=NDI 1,NTENSDDSDDE(K1,K1)=EGEND DOCC RECOVER ELASTIC AND PLASTIC STRAINS AND ROTATE FORWARDC ALSO RECOVER EQUIVALENT PLASTIC STRAINC读取弹性应变分量,塑性应变分量,并旋转(调用了ROTSIG),分别保存在EELAS和EPLAS中;CALL ROTSIG(STATEV(1),DROT,EELAS,2,NDI,NSHR)CALL ROTSIG(STATEV(NTENS 1),DROT,EPLAS,2,NDI,NSHR)读取等效塑性应变EQPLAS=STATEV(1 2*NTENS)先假设没有发生塑性流动,按完全弹性变形计算试算应力CC CALCULATE PREDICTOR STRESS AND ELASTIC STRAINCDO K1=1,NTENSDO K2=1,NTENSSTRESS(K2)=STRESS(K2) DDSDDE(K2,K1)*DSTRAN(K1)END DOEELAS(K1)=EELAS(K1) DSTRAN(K1)END DOC计算Mises应力C CALCULATE EQUIVALENT VON MISES STRESSCSMISES=(STRESS(1)-STRESS(2))**2 (STRESS(2)-STRESS(3))**21(STRESS(3)-STRESS(1))**2DO K1=NDI 1,NTENSSMISES=SMISES SIX*STRESS(K1)**2END DOSMISES=SQRT(SMISES/TWO)C根据当前等效塑性应变,调用HARDSUB得到对应的屈服应力C GET YIELD STRESS FROM THE SPECIFIED HARDENING CURVE CNVALUE=NPROPS/2-1CALL HARDSUB(SYIEL0,HARD,EQPLAS,PROPS(3),NVALUE)CC DETERMINE IF ACTIVELY YIELDINGC如果Mises应力大余屈服应力,屈服发生,计算流动方向IF(SMISES.GT.(ONE TOLER)*SYIEL0)THENCC ACTIVELY YIELDINGC SEPARATE THE HYDROSTATIC FROM THE DEVIATORIC STRESS C CALCULATE THE FLOW DIRECTIONCSHYDRO=(STRESS(1) STRESS(2) STRESS(3))/THREEDO K1=1,NDIFLOW(K1)=(STRESS(K1)-SHYDRO)/SMISESEND DODO K1=NDI 1,NTENSFLOW(K1)=STRESS(K1)/SMISESEND DOC根据J2理论并应用Newton-Rampson方法求得等效塑性应变增量C SOLVE FOR EQUIVALENT VON MISES STRESSC AND EQUIVALENT PLASTIC STRAIN INCREMENT USING NEWTON ITERATION CSYIELD=SYIEL0DEQPL=ZERODO KEWTON=1,NEWTONRHS=SMISES-EG3*DEQPL-SYIELDDEQPL=DEQPL RHS/(EG3 HARD)CALL HARDSUB(SYIELD,HARD,EQPLAS DEQPL,PROPS(3),NVALUE)IF(ABS(RHS).LT.TOLER*SYIEL0)GOTO10END DOCC WRITE WARNING MESSAGE TO THE.MSG FILECWRITE(7,2)NEWTON2FORMAT(//,30X,'***WARNING-PLASTICITY ALGORITHM DID NOT',1'CONVERGE AFTER',I3,'ITERATIONS')10CONTINUEC更新应力,应变分量C UPDATE STRESS,ELASTIC AND PLASTIC STRAINS ANDC EQUIVALENT PLASTIC STRAINCDO K1=1,NDISTRESS(K1)=FLOW(K1)*SYIELD SHYDROEPLAS(K1)=EPLAS(K1) THREE/TWO*FLOW(K1)*DEQPLEELAS(K1)=EELAS(K1)-THREE/TWO*FLOW(K1)*DEQPLEND DODO K1=NDI 1,NTENSSTRESS(K1)=FLOW(K1)*SYIELDEPLAS(K1)=EPLAS(K1) THREE*FLOW(K1)*DEQPLEELAS(K1)=EELAS(K1)-THREE*FLOW(K1)*DEQPLEND DOEQPLAS=EQPLAS DEQPLCC CALCULATE PLASTIC DISSIPATIONCSPD=DEQPL*(SYIEL0 SYIELD)/TWOCC计算塑性变形下的Jacobian矩阵FORMULATE THE JACOBIAN(MATERIAL TANGENT)C FIRST CALCULATE EFFECTIVE MODULICEFFG=EG*SYIELD/SMISESEFFG2=TWO*EFFGEFFG3=THREE/TWO*EFFG2EFFLAM=(EBULK3-EFFG2)/THREEEFFHRD=EG3*HARD/(EG3 HARD)-EFFG3c...if(props(7).lt..001)go to99c...DO K1=1,NDIDO K2=1,NDIDDSDDE(K2,K1)=EFFLAMEND DODDSDDE(K1,K1)=EFFG2 EFFLAMEND DODO K1=NDI 1,NTENSDDSDDE(K1,K1)=EFFGEND DODO K1=1,NTENSDO K2=1,NTENSDDSDDE(K2,K1)=DDSDDE(K2,K1) EFFHRD*FLOW(K2)*FLOW(K1) END DOEND DOc...99continuec...ENDIFC将弹性应变,塑性应变分量保存到状态变量中,并传到下一个增量步C STORE ELASTIC AND(EQUIVALENT)PLASTIC STRAINSC IN STATE VARIABLE ARRAYCDO K1=1,NTENSSTATEV(K1)=EELAS(K1)STATEV(K1 NTENS)=EPLAS(K1)END DOSTATEV(1 2*NTENS)=EQPLASCRETURNENDc...c...子程序,根据等效塑性应变,利用插值的方法得到对应的屈服应力SUBROUTINE HARDSUB(SYIELD,HARD,EQPLAS,TABLE,NVALUE)CINCLUDE'ABA_PARAM.INC'CDIMENSION TABLE(2,NVALUE)CPARAMETER(ZERO=0.D0)CC SET YIELD STRESS TO LAST VALUE OF TABLE,HARDENING TO ZERO CSYIELD=TABLE(1,NVALUE)HARD=ZEROC IF MORE THAN ONE ENTRY,SEARCH TABLECIF(NVALUE.GT.1)THENDO K1=1,NVALUE-1EQPL1=TABLE(2,K1 1)IF(EQPLAS.LT.EQPL1)THENEQPL0=TABLE(2,K1)IF(EQPL1.LE.EQPL0)THENWRITE(7,1)1FORMAT(//,30X,'***ERROR-PLASTIC STRAIN MUST BE`,1`ENTERED IN ASCENDING ORDER')CALL XITENDIFCC CURRENT YIELD STRESS AND HARDENINGCDEQPL=EQPL1-EQPL0SYIEL0=TABLE(1,K1)SYIEL1=TABLE(1,K1 1)DSYIEL=SYIEL1-SYIEL0HARD=DSYIEL/DEQPLSYIELD=SYIEL0 (EQPLAS-EQPL0)*HARD GOTO10ENDIFEND DO10CONTINUEENDIFRETURNENDWelcome To Download !!!欢迎您的下载,资料仅供参考!。

相关文档
最新文档