用于ABAQUS显式分析梁单元的混凝土单轴本构模型_王强

2011年7月第27卷第4期沈阳建筑大学学报（自然科学版）

Journal of Shenyang Jianzhu University （Natural Science ）Jul ．2011Vol ．27，No ．4

收稿日期：2010－12－31

基金项目：住房和城乡建设部科技基金项目（2008－K1－15）作者简介：王强（1971—），男，副教授，博士，主要从事工程结构抗震研究．

文章编号：2095－1922（2011）04－0679－06

用于ABAQUS 显式分析梁单元的

混凝土单轴本构模型

王

强，潘天林，刘

明，李

哲

（沈阳建筑大学土木工程学院，辽宁沈阳110168）

摘

要：目的为实现采用梁单元进行钢筋混凝土杆系结构的弹塑性响应分析，对其混凝土本

构关系进行二次开发，使ABAQUS 软件提供的混凝土材料模型能用于三维梁单元．方法利用

ABAQUS 用户自定义材料程序VUMAT 接口，开发用于显式动力分析的梁单元混凝土单轴本构模型，并编制相应的计算程序，对低周往复加载下的钢筋混凝土柱进行数值模拟计算．结果数值模拟结果能够较好地反映轴力对钢筋混凝土构件滞回性能的影响以及钢筋混凝土柱的双向弯曲耦合性能．结论笔者所开发的混凝土本构模型能够用于多维受力状态下钢筋混凝土梁柱构件的受力行为分析，

满足钢筋混凝土杆系结构动力弹塑性分析的需求．关键词：混凝土；滞回性能；本构模型；ABAQUS ；VUMAT 中图分类号：TU375.3

文献标志码：A

Study on a Uniaxial Constitutive Model of Concrete for Explicit Dynamic Beam Elements of ABAQUS

WANG Qiang ，PAN Tianlin ，LIU Ming ，LI Zhe

（School of Civil Engineering ，Shenyang Jianzhu University ，Shenyang ，China ，110168）

Abstract ：In order to use the beam element of FEM software ABAQUS for analyzing the elastic-plastic dy-namic response of RC truss structures ，

it is necessary to carry out a secondary development of the concrete constitutive for spatial beam element．In this paper ，a uniaxial constitutive model of concrete is established．The material subroutine of this model is successfully developed and applied to explicit dynamic module of

ABAQUS by means of user-defined subroutine interface VUMAT．Afterwards ，

the hysteretic performance of RC columns under cyclic loading is numerically simulated and compared with experiment results．The results show that the uniaxial constitutive model can rightly simulate the influence on the hysteretic performance of RC columns under varies axial load ，as well as the bi-axes bending coupling performance．The established model can meet the demand of analyzing the elastic-plastic dynamic response of RC frame structures．

Key words ：concrete ；hysteretic performance ；constitutive model ；ABAQUS ；VUMAT

混凝土结构在大震作用下通常会进入塑性状态，采用弹性分析方法进行结构的受力分析不能真实反映结构实际受力情况．进行结构的动力弹

塑性响应分析，特别是基于构件材料层次分析模型的弹塑性响应分析，能够较为准确地把握结构在大震作用下的非线性形态，对于评估结构的抗

680沈阳建筑大学学报（自然科学版）第27卷

震安全性具有重要意义．目前各国学者及工程界

已开始致力于此方面的研究

［1－4］

．通用有限元软件ABAQUS 具有较好的计算

稳定性、丰富的单元材料模型以及强大的前后处理功能，目前已在结构构件的非线性分析中得到

了广泛的应用［5－7］

，特别是其显式分析模块（ABAQUS /EXPLICIT ），由于其采用中心差分法

求解动力平衡方程，计算中无需形成结构的整体刚度矩阵，具有计算收敛性好的特点，更适于结构动力弹塑性响应分析．但在ABAQUS 显式分析模块中，软件提供的混凝土材料模型不能用于三维梁单元．若采用实体单元进行高层建筑等杆系结构的整体分析，

则计算工作量较大，难以满足工程计算需求．笔者基于纤维模型［8］

，利用ABAQUS 显式分析模块的用户自定义材料子程序VU-MAT ，对梁单元的混凝土材料模型进行二次开发，以满足结构动力弹塑性响应分析的需求．

1纤维梁单元模型

基于材料单轴本构关系的纤维模型是将构件

沿纵向划分为若干子段，再沿构件横截面划分成纤维束．每个纤维只考虑它的轴向本构关系，且可

定义不同的本构关系．柱横截面变形符合平截面

假定．对截面纤维的当前状态积分就可以得到截面的双向抗弯刚度、双向抵抗矩以及轴力，进而沿杆长进行积分，就可以得到精确的杆件单元刚度矩阵．纤维模型可以自然、简单地描述构件的双向弯曲－轴力耦合效应．1.1

基本假定

（1）构件截面变形满足平截面假定；（2）不考虑钢筋与混凝土之间的相对滑移；（3）不考虑构

件的剪切非线性及与其他变形的耦合关系．1.2单元截面刚度矩阵

梁单元类型为ABAQUS 显式分析模块中的B31梁单元［9］．该单元是基于铁摩辛柯（Timosh-enko ）梁理论构建的，可以考虑剪切变形．B31梁单元具有两个节点，一个积分点，转角和位移采用线性插值，如图1（a ）所示．采用GREEN 应变计算公式，可考虑大应变．单元质量阵为对角阵形式．采用矩形梁截面描述构件截面中的混凝土部

分，

将其划分为25个积分点或更多，如图1（b ）所示；同时采用箱型截面按等面积原则、等位置代替截面中的钢筋，划分为16个积分点或更多，如图1（c ）所示．每个积分点即为一个纤维

．

图1B31梁单元的积分点设置

Fig.1

Integration points of B31beam element

假设梁单元的横截面坐标轴分别为y 、z 轴，纵向坐标轴为x 轴．由单元节点位移通过插值函数可以得到轴向积分点处变形向量

d （x ）=｛Φz （x ）Φy （x ）ε0（x ）｝T ．（1）

根据截面积分点的位置，由轴向积分点处变

形向量可以得到纤维的应变向量

ε（x ）25?1=H 25?3d （x ）．

（2）

其中截面纤维几何位置转换矩阵H =［

H 1H 2…H 25］T

，H i =［－y i

z i

1］，i =1，2， (25)

由纤维的应变向量与材料的本构关系可得截面应力向量σ=E ε，其中E 为纤维切线刚度对角

阵．截面恢复力向量

F （x ）=｛M z

M y

N ｝T =H T A σ=H T AE ε=

H T AEH d （x ）．

（3）

式中：M z ，

M y ，N 分别为截面上绕y 、z 轴的弯矩及轴向力；A 为纤维面积对角阵．

整理可得单元截面的刚度矩阵为K sec =H T AEH ．

（4）

运用单元形函数矩阵，可以从截面刚度矩阵推得单元刚度矩阵K e =

∫l

B T K

sec

B d x．式中，B 为

单元形函数矩阵，

l 为单元长度．

第27卷王强等：用于ABAQUS 显式分析梁单元的混凝土单轴本构模型681

2

材料的本构模型

2.1

钢筋的本构模型

钢筋在反复荷载作用下本构模型采用ABAQUS 中自带的随动强化模型［9］，并考虑钢筋屈服硬化，钢筋屈服后刚度取E =0.01E 0，对应的单轴本构模型如图2所示．其中E 0为初始弹性模

量，E 为屈服后弹性模量，f y 为屈服应力，εy 为屈服应变

．

图2钢筋的本构模型

Fig.2

Constitutive model of steel

2.2

混凝土的本构模型

笔者采用基于文献［

10］提出的混凝土本构模型，如图3所示．其中E c 0为原点切线模量；E cr 为损伤后弹性模量；εcm 为混凝土所经历的最大压应变；f c 为混凝土抗压强度；ε0为混凝土峰值应力所对应的应变，ε0=0.002；εu 为混凝土的极限压应变，

εu =4ε0．混凝土受压骨架曲线采用Kent 和Park 所提出并由B．D．Scott 改进的混凝土应

力－应变曲线［11］

．由于混凝土的抗拉强度很低，且在滞回过程中一旦开裂，混凝土就不能再承受

拉力，

因此抗拉强度对混凝土构件滞回性能影响较小

［12－13］

．故在本构模型中忽略混凝土的抗拉强度，并忽略裂面效应影响．

混凝土卸载及再加载曲线均取为直线形式．

卸载时考虑刚度的退化，卸载模量按式（7）确定：

E cr =

E c0εc ≤ε0，

E c0

ε0

ε()

cm

0.9

εc ＞ε0{

．

（7）

当混凝土卸载至零压应力时，如继续卸载则材料应力保持为零．若混凝土卸载至零压应力之前又开始加载，则沿原路径返回．再加载时加载曲线始终指向骨架曲线上所经历的最大应变点．若应力超过骨架曲线上的相应点，则沿骨架曲线加

载．若混凝土应力在达到骨架曲线之前开始卸载，则按照所经历的最大压应变来确定卸载刚度

．

图3混凝土的本构模型

Fig.3

Constitutive model of concrete

3

用户自定义材料子程序（VU-MAT ）实现

依据前述的混凝土本构模型，笔者基于用户

自定义材料子程序VUMAT 接口，编制了计算程

序，并嵌入ABAQUS /EXPLICIT 模块中［9］．主程

序通过ABAQUS 输入文件中的关键字

“*USER-MATERIAL ”来判断是否使用了用户自定义材

料，并提供混凝土本构模型所需的材料参数［11］

．在ABAQUS 中对编制的VUMAT 子程序进行调

试，来跟踪每一步调用子程序时变量的更新情况，从而及时发现所产生的错误．调试时要在com-mand 窗口中输入“abaqus －j 文件名．inp －user

程序名．for －debug －explicit ”，在VISUAL STU-DIO 开发环境中打开子程序，然后设置断点进行

调试．在VUMAT 中只有程序中定义的数组和变量能够进行新旧变量更替，如果另定义更新变量必须特别声明存储特性，否则子程序不会保存上一步变量数值．编程中还应避免除零问题．

为保证程序编制思路的可靠性，笔者在进行混凝土本构模型开发之前，首先编制了理想弹塑性材料的VUMAT ，并与ABAQUS 自带理想弹塑性模型进行对比，得到的结果基本一致．

4算例验证

为充分验证模型的有效性，笔者分别对不同

加载制度下的钢筋混凝土柱滞回性能进行计算分

析．试件情况见文献［14］，构造和配筋如图4所

682沈阳建筑大学学报（自然科学版）第27卷

示．各试件的加载规则见表1，其中试件SP1与SP2

为笔者构造的加载制度，

SP3与SP4则为文献［14］中的试件TP74和TP77．材料参数取值见表

2．

图4

钢筋混凝土柱试件Fig.4

Reinforced concrete column specimens

表1

算例加载制度

Table 1

Loading pattern of example

试件编号

加载图示

加载规则

轴力/kN

备注

SP1－轴向往复加载

SP20

无轴压单向往复侧推

SP3160

有轴压单向往复侧推

SP4160

有轴压双向往复侧推

表2

材料基本参数

Table 2

Basic parameters of materials

参数项屈服强度/MPa 屈服应变泊松比弹性模量/104

MPa 混凝土29．660．0020．252钢筋

357

0．0017

0．3

20

由于采用显式动力方法进行拟静力分析，必

须减小惯性力对整个构件的影响．采取的措施是降低加载速率和减小计算时步，这样可以使加速度趋近于很小，从而忽略惯性力影响．图5为计算所得试件SP1在轴向往复拉压

时，

ABAQUS 计算输出的角部混凝土纤维受压应力应变曲线（压为负）．该曲线符合笔者所给出混凝土的本构模型，表明笔者编制的材料本构子程序是正确的．图6为试件SP2计算所得的水平加载滞回曲线．可以看出无轴压时构件的滞回曲线呈梭形，且较为饱满，属于典型的受弯构件滞回性能［15］

．而且对试验的“超前指向”现象也有所表现，

即加载曲线并不指向前一循环的开始卸载点，而是指向前一循环的开始卸载点位移更大的一点

．

图5试件SP1角部混凝土纤维的应力应变关系Fig.5

Stress-strain relationship of corner concrete fiber of specimen

SP1

图6试件SP2计算所得滞回曲线

Fig.6

Calculated hysteresis curve of specimen SP2

图7、

8分别为试件SP3的试验实测与计算所得滞回曲线，二者对比可以发现在加载初期0 20mm 时实验曲线与计算曲线基本一致，只是峰值点处计算值略小，这可能是对于混凝土受箍筋

第27卷王强等：用于ABAQUS 显式分析梁单元的混凝土单轴本构模型683

约束使得强度提高考虑不足．在后期加载20 60mm 时，计算所得滞回曲线较为丰满．造成此现象的原因主要是没有考虑钢筋的滑移，特别是加载后期实际构件已产生滑移，而计算模型并没有表现出来．而且采用的钢筋本构模型为线性强化模型，与真实钢筋的本构关系有一定误差，耗能能力更强一些，

所以导致计算所得的滞回曲线比试验所得的曲线要饱满一些．对于试验结果中的“超前指向”现象，计算结果同样能够予以较好的描述．此外由图8与图6对比可以看出轴压力的存在使得构件极限承载力略有提高，而滞回曲线产生捏拢现象

．

图7试件SP3实测滞回曲线

Fig.7

Hysteresis curve of specimen

SP3

图8

试件SP3计算所得滞回曲线

Fig.8

Calculated hysteresis curve of specimen SP3

图9、

10分别为试件SP4的实验与计算结果．由计算结果可以看出，当方向1保持位移恒定，方向2的加载使得方向1产生荷载跌落现象，反之亦然，这在试验曲线中有相应的体现．可以认为计算模型能够较好地反映钢筋混凝土柱的双向弯曲耦合性能．计算所得滞回曲线仍较试验曲线丰满，计算峰值略低于实验值

．

图9SP4试验滞回曲线Fig.9

Hysteresis curve of

SP4

图10

SP4计算滞回曲线

Fig.10

Calculated hysteresis curve of SP4

684沈阳建筑大学学报（自然科学版）第27卷

5结论

（1）笔者建立的模型可以正确反映轴力对钢筋混凝土构件滞回性能的影响，能够较好地模拟钢筋混凝土柱的双向弯曲耦合性能以及滞回曲线中的超前指向与捏拢现象，可以用于多维受力状态下钢筋混凝土梁柱构件的受力行为分析，能够满足空间框架结构动力弹塑性分析的需求．（2）采用箱型截面等效代替考虑钢筋混凝土杆件中的钢筋，有效地解决了杆件采用梁单元模型时难以考虑钢筋作用的问题．

（3）由于采用的模型未考虑钢筋的滑移，对整个结果的精确性有一定的影响，有待于进一步研究．

参考文献：

［1］秦从律，张爱晖．基于截面纤维模型的弹塑性时程分析方法［J］．浙江大学学报，2005（7）：1003－

1008．

（Qin Conglü，Zhang Aihui．Non linear time history

analysis based on section fiber model［J］．Journal of

Zhejiang University（Engineering Science），2005

（7）：1003－1008．）

［2］Li Kangning．3-D Analysis of RC f rame-w al l building damaged in the1995hyogoken-nanbu earth-

quake［C］//Process of the12th World Conference on

Earthquake Engineering，New Zealand：Auckland，

2000．

［3］Mazzoni S，Mckenna F，Scott M H，et al．Opensees command language manual［R］．PEER，Berkeley：U-

niversity of California，2004．

［4］汪训流，陆新征，叶列平．往复荷载下钢筋混凝土柱受力性能的数值模拟［J］．工程力学，2007（12）：

76－81．

（Wang Xunliu，Lu Xinzheng，Ye Lieping．Numerical

simulation for the hysteresis behavior of RC columns

under cyclic loads［J］Engineering Mechanics，2007，

24（12）：76－81．）

［5］王金昌，陈页开．ABAQUS在土木工程中的应用［M］．杭州：浙江大学出版社，2006．

（Wang Jinchang，Chen Yekai．Application of

ABAQUS in civil engineering［M］．Hangzhou：Zhe-

jiang University Press，2006．）

［6］Navalurkar R K，Hsu C T．Fracture analysis of high strength concrete members［J］．Journal of Materials

in Civil Engineering，2001，13（3）：185－193．

［7］Chung W，Sotelino E D．Nonlinear finite-element a-nalysis of composite steel girder bridges［J］．Journal

of Structural Engineering，2005，131（2）：304－313．［8］Spacone E，Fillippou F C，Taucer F．Fiber beam-col-umn model for non-liner analysis of R/C frames

［J］．Journal of Earthquake Engineering and Structur-

al Dynamics，1996，25：711－725．

［9］庄茁，由小川，廖剑辉，等．基于ABAQUS的有限元分析和应用［M］．北京：清华大学出版社，2009．

（Zhuang Zhuo，You Xiaochuan，Liao Jianhui，et al．

FEM analysis and application based on ABAQUS

［M］．Beijing：Tsinghua University Press，2009．）［10］王强，吕西林，雷淑忠．离散单元法在钢筋混凝土构件非线性分析的应用［J］．沈阳建筑大学学报：

自然科学版，2005（2）：91－95．

（Wang Qiang，Lu Xilin，Lei Shuzhong．Application

of the DEM on the nonlinear analysis of reinforced

concrete members［J］．Journal of Shenyang Jianzhu

University：Natural Science，2005（2）：91－95．）［11］Scott B D，Park R，Priestley M J N．Stress-strain be-havior of concrete confined overlapping hoops at low

and high strain rates［J］．ACI Journal，1982，79（1）：

13－27．

［12］Légeron F，Paultre P．Uniaxial confinement model for normal and high-strength concrete columns［J］．Jour-

nal of Structural Engineering，2003，129（2）：241－

252．

［13］Mander J B，Priestley M J N，Park R．Theoretical stress-strain model for confined concrete［J］．Journal

of Structural Engineering，1988，114（8）：1804－

1825．

［14］Kawashima K，Watanabe G，Hayakawa R．Seismic performance of RC bridge columns subjected to bi-

lateral excitation［C］//Proc．35th joint meeting，pan-

el on wind and seismic effects，Japan：Tsukuba Sci-

ence City，2003．

［15］张新培．钢筋混凝土抗震结构非线性分析［M］．北京：科学出版社，2003．

（Zhang Xinpei．Nonlinear seismic study on rein-

forced concrete structures［M］．Beijing：Science

Press，2003．）

ABAQUS中Cohesive单元建模方法

复合材料模型建模与分析 1. Cohesive单元建模方法 几何模型 使用内聚力模型（cohesive zone）模拟裂纹的产生和扩展，需要在预计产生裂纹的区域加入cohesive层。建立cohesive层的方法主要有： 方法一、建立完整的结构（如图1（a）所示），然后在上面切割出一个薄层来模拟cohesive 单元，用这种方法建立的cohesive单元与其他单元公用节点，并以此传递力和位移。 方法二、分别建立cohesive层和其他结构部件的实体模型，通过“tie”绑定约束，使得cohesive单元两侧的单元位移和应力协调，如图1（b）所示。 （a）cohesive单元与其他单元公用节点（b）独立的网格通过“tie”绑定 图1.建模方法 上述两种方法都可以用来模拟复合材料的分层失效，第一种方法划分网格比较复杂；第二种方法赋材料属性简单，划分网格也方便，但是装配及“tie”很繁琐；因此在实际建模中我们应根据实际结构选取较简单的方法。 材料属性 应用cohesive单元模拟复合材料失效，包括两种模型：一种是基于traction-separation描述；另一种是基于连续体描述。其中基于traction-separation描述的方法应用更加广泛。 而在基于traction-separation描述的方法中，最常用的本构模型为图2所示的双线性本构模型。它给出了材料达到强度极限前的线弹性段和材料达到强度极限后的刚度线性降低软化阶段。注意图中纵坐标为应力，而横坐标为位移，因此线弹性段的斜率代表的实际是cohesive单元的刚度。曲线下的面积即为材料断裂时的能量释放率。因此在定义cohesive

ABAQUS中的三种混凝土本构模型(20200706140516)

ABAQUS用连续介质的方法建立描述混凝土模型不采用宏观离散裂纹的方法描述裂纹的水平的在每一个积分点上单独计算其中。 低压力混凝土的本构关系包括： Con crete Smeared cracki ng model (ABAQUS/Sta ndard) Concrete Brittle cracki ng model (ABAQUS/Explicit) Con crete Damage plasticity model 高压力混凝土的本构关系： Cap model 1、ABAQUS/Standard 中的弥散裂缝模型Concrete Smeared cracking model (ABAQUS/Standard): 只能用于ABAQUS/Standard 中 裂纹是影响材料行为的最关键因素，它将导致开裂以及开裂后的材料的各向异性 用于描述：单调应变、在材料中表现出拉伸裂纹或者压缩时破碎的行为 在进行参数定义式的Keywords: *CONCRETE *TENSION STIFFENING *SHEAR RETENTION *FAILURE RATIOS 2、ABAQUS/Explicit 中脆性破裂模型Concrete Brittle cracking model (ABAQUS/Explicit): 适用于拉伸裂纹控制材料行为的应用或压缩失效不重要，此模型考虑了由于裂纹引起的材料 各向异性性质，材料压缩的行为假定为线弹性，脆性断裂准则可以使得材料在拉伸应力过大 时失效。 在进行参数定义式的Keywords *BRITTLE CRACKING, *BRITTLE FAILURE, *BRITTLE SHEAR 3、塑性损伤模型Concrete Damage plasticity model : 适用于混凝土的各种荷载分析，单调应变，循环荷载，动力载荷，包含拉伸开裂(cracking)和压缩破碎(crushing)，此模型可以模拟硬度退化机制以及反向加载刚度恢复的混凝土力学特性 在进行参数定义式的Keywords: *CONCRETE DAMAGED PLASTICITY *CONCRETE TENSION STIFFENING *CONCRETE COMPRESSION HARDENING *CONCRETE TENSION DAMAGE

abaqus 中梁板壳单元的弯曲问题beamplateshell

ABAQUS中梁板壳单元的弯曲问题 曲哲 2007-4-3 一、Euler-Bernoulli梁与Timoshenko梁 在ABAQUS的单元库中，所有三次插值的梁单元（如B23，B33等），均为Euler-Bernoulli梁，而所 有线性和二次插值的梁单元（如B21，B22，B31，B32等），均为Timoshenko梁。 （1）细长梁与深梁 B23为2结点三次插值的Euler-Bernoulli梁。由于在形成单元刚度矩阵时等效载荷项的被积函数至少 是3次的，所以至少需要2个积分点才能达到完全的高斯积分。而在ABAQUS中，B23有3个积分点， 这意味着被积函数可以达到5次。总之B23是完全积分的单元。而B21和B22分别为2结点线性插值和3 结点二次插值的Timoshenko梁，并且默认的采用减缩积分来避免剪切锁死。B22只有2个积分点，B21 只有1个积分点，它们都只能达到1次的插值精度。 表1：集中力作用下悬臂梁的自由端挠度（mm） 细长梁（l/h=10）深梁（l/h=3） 材力解 1 2 4 材力解 4 2 单元个数 1 0.1080 0.1080 4.000 4.000 4.000 0.10800.1080 B23（E-B梁） 4.000 B21（Timoshenko梁） 3.734 3.955 4.010 4.000 0.10860.1145 0.1160 0.1080 B22（Timoshenko梁） 4.028 4.028 4.028 4.000 0.11650.1165 0.1165 0.1080 表1比较了上述三种梁单元在应用于细长梁和深梁受弯时的表现。问题描述如图1所示，为端部受集 中载荷的悬臂梁。E-B梁B23完全忠实于材料力解的解答，不考虑剪切应变的影响，并且只用1个B23单 元就可以得到与材力解一致的结果。B21和B22考虑了梁的剪切变形，其分析得到的挠度略大于材力解。 同时可以看出，B21和B22用于细长梁时并没有发生剪切自锁。 图1：悬臂梁的构型图与截面图 图2：网格划分（2个单元）

ABAQUS应用梁单元计算简支梁

ABAQUS应用梁单元计算简支梁 对于梁的分析可以使用梁单元、壳单元或是固体单元。Abaqus的梁单元需要设定线的方向，用选中所需要的线后，输入该线梁截面的主轴1方向单位矢量(x,y,z)，截面的主轴方向在截面Profile设定中有规定。 注意： 因为ABAQUS软件没有UNDO功能，在建模过程中，应不时地将本题的CAE模型（阶段结果）保存，以免丢失已完成的工作。 简支梁，三点弯曲，工字钢构件，结构钢材质，E=210GPa，μ=0.28，ρ=7850kg/m3（在不计重力的静力学分析中可以不要）。F=10kN，不计重力。计算中点挠度，两端转角。理论解：I=2.239×10-5m4，w中=2.769×10-3m，θ边=2.077×10-3。 文件与路径： 顶部下拉菜单File, Save As ExpAbq00。 一部件 1 创建部件：Module，Part，Create Part， 命名为Prat-1；3D，可变形模型，线，图形大约范围10(程序默认长度单位为m)。 2 绘模型图：选用折线，从(0,0)→(2,0)→(4,0)绘出梁的轴线。 3 退出：Done。 二性质 1 创建截面几何形状：Module，Property，Create Profile， 命名为Profile-1，选I型截面，按图输入数据，l=0.1，h=0.2，b l=0.1，b2=0.1，t l=0.01，t2=0.01，t3=0.01，关闭。 2 定义梁方向：Module，Property，Assign Beam Orientation， 选中两段线段，输入主轴1方向单位矢量(0,0,1)或(0,0,-1)，关闭。 3 定义截面力学性质：Module，Property，Create Section， 命名为Section-1，梁，梁，截面几何形状选Profile-1，输入E=210e9（程序默认单位为N/m2，GPa=109 N/m2）， G=82.03e9，ν=0.28，关闭。

ABQUS中的三种混凝土本构模型

. ABQUS中的三种混凝土本构模型 ABAQUS 用连续介质的方法建立描述混凝土模型不采用宏观离散裂纹的方法描述裂纹的水平的在每一个积分点上单独计算其中。 低压力混凝土的本构关系包括： Concrete Smeared cracking model (ABAQUS/Standard) Concrete Brittle cracking model (ABAQUS/Explicit) Concrete Damage plasticity model 高压力混凝土的本构关系： Cap model 1、ABAQUS/Standard中的弥散裂缝模型Concrete Smeared cracking model (ABAQUS/Standard)：——只能用于ABAQUS/Standard中 裂纹是影响材料行为的最关键因素，它将导致开裂以及开裂后的材料的各向异性 用于描述：单调应变、在材料中表现出拉伸裂纹或者压缩时破碎的行为 在进行参数定义式的Keywords： *CONCRETE *TENSION STIFFENING *SHEAR RETENTION *FAILURE RATIOS 2、ABAQUS/Explicit中脆性破裂模型Concrete Brittle cracking model (ABAQUS/Explicit) ： 适用于拉伸裂纹控制材料行为的应用或压缩失效不重要，此模型考虑了由于裂纹引起的材料各向异性性质，材料压缩的行为假定为线弹性，脆性断裂准则可以使得材料在拉伸应力过大时失效。 在进行参数定义式的Keywords *BRITTLE CRACKING, *BRITTLE FAILURE, *BRITTLE SHEAR 3、塑性损伤模型Concrete Damage plasticity model： 适用于混凝土的各种荷载分析，单调应变，循环荷载，动力载荷，包含拉伸开裂(cracking)和压缩破碎(crushing)，此模型可以模拟硬度退化机制以及反向加载刚度恢复的混凝土力学特性 在进行参数定义式的Keywords： *CONCRETE DAMAGED PLASTICITY *CONCRETE TENSION STIFFENING *CONCRETE COMPRESSION HARDENING *CONCRETE TENSION DAMAGE *CONCRETE COMPRESSION DAMAGE 1 / 1'.

ABAQUS中的三种混凝土本构模型

ABQUS中的三種混凝土本構模型 ABAQUS?用連續介質的方法建立描述混凝土模型不采用宏觀離散裂紋的方法描述裂紋的水平的在每一個積分點上單獨計算其中。 低壓力混凝土的本構關系包括： Concrete Smeared cracking model (ABAQUS/Standard) Concrete Brittle cracking model (ABAQUS/Explicit) Concrete Damage plasticity model 高壓力混凝土的本構關系： Cap model 1、ABAQUS/Standard中的彌散裂縫模型Concrete Smeared cracking model (ABAQUS/Standard)： ——只能用于ABAQUS/Standard中 裂紋是影響材料行為的最關鍵因素，它將導致開裂以及開裂后的材料的各向異性 GAGGAGAGGAFFFFAFAF

用于描述?：單調應變?、在材料中表現出拉伸裂紋或者壓縮時破碎的行為 在進行參數定義式的Keywords： *CONCRETE *TENSION STIFFENING *SHEAR RETENTION *FAILURE RATIOS 2、ABAQUS/Explicit中脆性破裂模型Concrete Brittle cracking model (ABAQUS/Explicit) ： 适用于拉伸裂纹控制材料行为的应用或压缩失效不重要，此模型考虑了由于裂纹引起的材料各向异性性质，材料压缩的行为假定为线弹性，脆性断裂准则可以使得材料在拉伸应力过大时失效。 在进行参数定义式的Keywords *BRITTLE CRACKING, *BRITTLE FAILURE, *BRITTLE SHEAR GAGGAGAGGAFFFFAFAF

ABAQUS简支梁分析(梁单元和实体单元)

基于ABAQUS简支梁受力和弯矩的相关分析 （梁单元和实体单元） 对于简支梁，基于 ABAQUS2016，首先用梁单元分析了梁受力作用下的应力，变形，剪力和力矩；对同一模型，并用实体单元进行了相应的分析。另 外，还分析了梁结构受力和弯矩作用下的剪力及力矩分析。 对于CAE仿真分析具体细节操作并没有给出详细的操作，不过在后面上传了对应的cae，odb，inp文件。不过要注意的是本文采用的是ABAQUS2016进行计算，低版本可能打不开，可以自己提交inp文件自己计算即可。可以到小木虫搜索：“基于ABAQUS简支梁受力和弯矩的相关分析”进行相应文件下载。 对于一简支梁，其结构简图如下所示，梁的一段受固支，一段受简支，在梁的两端受集中载荷，梁的大直径D=180mm，小直径d=150mm，a=200mm，b=300mm，l=1600mm，F=300000N。现通过梁单元和实体单元分析简支梁的受力情况，变形情况，以及分析其剪力和弯矩等。材料采用45#钢，弹性模量 E=2.1e6MPa，泊松比v=0.28。 图1 简支梁结构简图 1.梁单元分析 ABAQUS2016中对应的文件为beam-shaft.cae ，beam-shaft.odb，beam-shaft.inp。 在建立梁part的时候，采用三维线性实体，按照图1所示尺寸建立，然后在台阶及支撑梁处进行分割，结果如图2所示。

图2 建立part并分割 接下来为梁结构分配材料，创建材料，定义弹性模量和泊松比，创建梁截面形状，如图3，非别定义两个圆，圆的直接分别为180和150mm。然后创建两个截面，截面选择梁截面，再选择图2中的所有梁，定义梁的方向矢量为（0，0，-1）（点击图3中的n2，n1，t那个图标即可创建梁的方向矢量），最后把创建好的梁赋给梁结构。 图3 创建梁截面形状 接下来装配实体，再创建分析步，在创建分析步的时候，点击主菜单栏的Output，编辑Edit Field Output Request，在SF前面打钩，这样就可以在结果后处理中输出截面剪力和力矩，如图4所示。在Load加载中，在固支处剪力边界条件，约束x，y，z，及绕x和y轴的转动，如图5所示，同理，在固支另一处约束y，z，及绕x和y轴的转动。在梁的两端添加集中力，集中力的大小为300000N。最后对实体部件进行分网，采用B32梁单元，网格尺寸为10。完成

ABAQUS钢筋混凝土损饬塑性模型有限元分析

ABAQUS钢筋混凝土损饬塑性模型有限元分析 发表时间：2009-10-12 刘劲松刘红军来源：万方数据 钢筋混凝土材料，是一种非匀质的力学性能复杂的建筑材料。随着计算机和有限元方法的发展，有限元法已经成为研究混凝土结构的一个重要的手段。由于数值计算具有快速、代价低和易于实现等诸多优点，这种分析方法已经广泛用于实际工程中。然而，要在有限元软件中尽可能准确地模拟混凝土这种材料，是不容易的，国内外学者提出了基于各种理论的混凝土本构模型。但是迄今为止，还没有一种理论被公认为可以完全描述混凝土的本构关系。 ABAQUS是大型通用的有限元分析软件，其在非线性分析方面的巨大优势，获得了广大用户的认可，在结构分析领域的应用趋于广泛。本文把规范建议的混凝土本构关系，应用到损伤塑性模型，对一悬臂梁进行了精细的有限元建模计算和探讨。 1 混凝土损伤塑性模型 ABAQUS在钢筋混凝土分析上有很强的能力。它提供了三种混凝土本构模型：混凝土损伤塑性模型，混凝土弥散裂缝模型和ABAQUS/Explicit中的混凝土开裂模型。其中混凝土损伤塑性模型可以用于单向加载、循环加载以及动态加载等场合，它使用非关联多硬化塑性和各向同性损伤弹性相结合的方式描述了混凝土破碎过程中发生的不可恢复的损伤。这一特性使得损伤塑性模型具有更好的收敛性。 2 模型材料的定义 2.1 混凝土的单轴拉压应力-应变曲线 本模型中选用的混凝土本构关系是《混凝土结构设计规范》所建议的曲线，其应力应变关系可由函数表达式定义。 2.2 钢筋的本构关系 钢筋采用本构关系为强化的二折线模型，无刚度退化。折线第一上升段的斜率，为钢筋本身的弹性模量，第二上升段为钢筋强化段，此时的斜率大致可取为第一段的1/100。 2.3 损伤的定义 损伤是指在单调加载或重复加载下，材料性质所产生的一种劣化现象，损伤在宏观方面的表现就是(微)裂纹的产生。材料的损伤状态，可以用损伤因子来描述。根据前面确定的混凝土非弹性阶段的应力一应变关系。可求得损伤因子的数值。 2.4混凝土塑性数值的计算 混凝土在单向拉伸，压缩试验中得到的数据，通常是以名义应变和名义应力表示的，为了准确地描述大变形过程中截面积的改变，需要使用真实应变和真实应力，可通过它们之间的换算公式计算。真实应变是由塑性应变和弹性应变两部分构成的。在ABAQUS中定义塑性材料参数时，需要使用塑性应变。 3 钢筋混凝土悬臂梁实例分析 3.1 模型设计 该悬臂梁的具体情况如图1所示，梁截面尺寸为200mm×300mm，梁长1500mm；纵筋为HRB335钢筋，箍筋为HPB235钢筋，混凝土强度等级为C30。混凝土和钢筋的各力学参数均取自《混凝土结构设计规范》的标准值。

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

ABQUS中的三种混凝土本构模型 ABAQUS 用连续介质的方法建立描述混凝土模型不采用宏观离散裂纹的方法描述裂纹的水平的在每一个积分点上单独计算其中。 低压力混凝土的本构关系包括： Concrete Smeared cracking model (ABAQUS/Standard) Concrete Brittle cracking model (ABAQUS/Explicit) Concrete Damage plasticity model 高压力混凝土的本构关系： Cap model 1、ABAQUS/Standard中的弥散裂缝模型Concrete Smeared cracking model (ABAQUS/Standard)：——只能用于ABAQUS/Standard中 裂纹是影响材料行为的最关键因素，它将导致开裂以及开裂后的材料的各向异性 用于描述：单调应变、在材料中表现出拉伸裂纹或者压缩时破碎的行为 在进行参数定义式的Keywords： *CONCRETE *TENSION STIFFENING *SHEAR RETENTION *FAILURE RATIOS 2、ABAQUS/Explicit中脆性破裂模型Concrete Brittle cracking model (ABAQUS/Explicit) ： 适用于拉伸裂纹控制材料行为的应用或压缩失效不重要，此模型考虑了由于裂纹引起的材料各向异性性质，材料压缩的行为假定为线弹性，脆性断裂准则可以使得材料在拉伸应力过大时失效。 在进行参数定义式的Keywords *BRITTLE CRACKING, *BRITTLE FAILURE, *BRITTLE SHEAR 3、塑性损伤模型Concrete Damage plasticity model： 适用于混凝土的各种荷载分析，单调应变，循环荷载，动力载荷，包含拉伸开裂(cracking)和压缩破碎(crushing)，此模型可以模拟硬度退化机制以及反向加载刚度恢复的混凝土力学特性 在进行参数定义式的Keywords： *CONCRETE DAMAGED PLASTICITY *CONCRETE TENSION STIFFENING *CONCRETE COMPRESSION HARDENING *CONCRETE TENSION DAMAGE *CONCRETE COMPRESSION DAMAGE

(仅供参考)Abaqus混凝土损伤塑性模型的参数标定

Abaqus 混凝土损伤塑性模型的参数标定 1. 塑性参数（Plasticity ） 1) 剪胀角（Dilation Angle ） = 30° 2) 流动势偏移量（Eccentricity ） 3) 双轴受压与单轴受压极限强度比 = 1.16 4) 不变量应力比 = 0.667 5) 粘滞系数（Visosity Parameter ） = 0.0005 2. 受压本构关系 应力-Yield Stress ：第一行应输入本构模型刚进入非弹性段非弹性应变为0时所对应的应力。 非弹性应变-Inelastic Strain （受拉时为开裂应变-Cracking Strain ）：根据应力按混凝土本构模型得出对应的应变值，并通过 , 和 ，得出非弹性应变。 3. 受压损伤因子（Damage Parameter ）计算 根据《Abaqus Analysis User's Manual (6.10)》 - 20.6.3 “Concrete damaged plasticity ”中公式： 假设非弹性应变 in c ε中塑性应变 pl c ε所占的比例为c β，通过转换可得损伤因子c d 的计算公式： () () 0 011in c c in c c c c E E d βεσβε-=+- 根据《ABAQUS 混凝土损伤塑性模型参数验证》规定，混凝土受压时c β的取值范围为0.35 ~ 0.7。

4. 受拉损伤因子（Damage Parameter ）计算 受拉损伤因子的计算与受压损伤因子的计算方法基本相同，只需将对应受压变量更换为受拉即可： () () 0011in t t in t t t t E E d βεσβε-=+- 而根据参考文献混凝土受拉时t β的取值范围为0.5 ~ 0.95。 5. 损伤恢复因子 受拉损伤恢复因子（Tension Recovery ）：缺省值0t w =。 受压损伤恢复因子（Compression Recovery ）：缺省值1c w =。

abaqus有限元分析报告开裂梁要点

Abaqus梁的开裂模拟计算报告 1.问题描述 利用ABAQUS有限元软件分析如图1.1所示的钢筋混凝土梁的裂缝开展。参考文献Brena et al.（2003）得到梁的基本数据： 图1.1 Brena et al.（2003）中梁C尺寸 几何尺寸：跨度3000mm，截面宽203mm，高406mm的钢筋混凝土梁 由文献Chen et al. 2011得材料特性： 1.混凝土：抗压强度f c’=35.1MPa，抗拉强度f t= 2.721MPa，泊松比ν=0.2，弹性模量 E c=28020MPa； 2.钢筋：弹性模量为E c=200GPa，屈服强度f ys=f yc=440MPa，f yv=596MPa 3.混凝土垫块：弹性模量为E c=28020MPa，泊松比ν=0.2 2.建模过程 1)Part 打开ABAQUS使用功能模块，弹出窗口Create Part，参数为：Name：beam；Modeling Space：2D；Type：Deformable；Base Feature─Shell；Approximate size：2000。点击Continue 进入Sketch二维绘图区。由于该梁关于Y轴对称，建模的时候取沿X轴的一半作为模拟对象。 使用功能模块，分别键入独立点(0,0),(1600,0),(1600,406),(406,0),(0,0)并按下下方提 示区的Done，完成草图。 图2.1 beam 部件二维几何模型

相同的方法建立混凝土垫块： 图2.2 plate 部件二维几何模型 所选用的点有(0,0),(40,0),(40,10),(0,10) 受压区钢筋： 在选择钢筋的base feature的时候选择wire，即线模型。 图2.3 compression bar 部件二维几何模型 选取的点(0,0),(1575,0) 受拉区钢筋： 图2.4 tension bar 部件二维几何模型 选取的点(0,0),(1575,0) 箍筋： 图2.5 stirrup 部件二维几何模型 选取的点为(0,0),(0,330) 另外，此文里面为了作对比，部分的模型输入尺寸的时候为m，下面无特别说明尺寸都为mm。

ABAQUS简支梁分析(梁单元和实体单元)

基于ABAQUS 简支梁受力和弯矩的相关分析 （梁单元和实体单元） 对于简支梁，基于ABAQUS2016，首先用梁单元分析了梁受力作用下的应 力，变形，剪力和力矩；对同一模型，并用实体单元进行了相应的分析。另 外，还分析了梁结构受力和弯矩作用下的剪力及力矩分析。 对于CAE 仿真分析具体细节操作并没有给出详细的操作，不过在后面上 传了对应的cae, odb , inp 文件。不过要注意的是本文采用的是 ABAQUS2016 进行计算，低版本可能打不开，可以自己提交 inp 文件自己计算即可。可以到 小木虫搜索：“基于ABAQUS 简支梁受力和弯矩的相关分析”进行相应文件 下载。 对于一简支梁，其结构简图如下所示，梁的一段受固支，一段受简支，在 梁的两端受集中载荷，梁的大直径 D=180mm ，小直径d=150mm ，a=200mm ， b=300mm , l=1600mm , F=300000N 。现通过梁单元和实体单元分析简支梁的受 力情况，变形情况，以及分析其剪力和弯矩等。材料采用 45#钢，弹性模量 E=2.1e6MPa,泊松比 v=0.28。 1.梁单元分析 ABAQUS2016 中对应的文件为 beam-shaft.cae , beam-shaft.odb , beam- shaft.inp 。 在建立梁part 的时候，采用三维线性实体，按照图1所示尺寸建立，然后 在台阶及支撑梁处进行分割，结果如图 2所示 l b b a a A A C B A 图1简支梁结构简图

图2建立part并分割 接下来为梁结构分配材料，创建材料，定义弹性模量和泊松比，创建梁截面形状，如图3,非别定义两个圆，圆的直接分别为180和150mm。然后创建两个截面，截面选择梁截面，再选择图2中的所有梁，定义梁的方向矢量为 （0，0，-1）（点击图3中的n2, n 1，t那个图标即可创建梁的方向矢量），最后把创建好的梁赋给梁结构。 图3创建梁截面形状 接下来装配实体，再创建分析步，在创建分析步的时候，点击主菜单栏的Output，编辑Edit Field Output Request,在SF前面打钩，这样就可以在结果后处理中输出截面剪力和力矩，如图4所示。在Load加载中，在固支处剪力边界条件，约束x，y，z，及绕x和y轴的转动，如图5所示，同理，在固支另一处约束y，z，及绕x和y轴的转动。在梁的两端添加集中力，集中力的大小为300000N。最后对实体部件进行分网，采用B32梁单元，网格尺寸为10。完成

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

Abaqus混凝土材料塑性损伤模型浅析与参数设置 【壹讲壹插件】欢迎转载，作者：星辰-北极星，QQ群：431603427 https://www.wendangku.net/doc/5616111117.html, Abaqus混凝土材料塑性损伤模型浅析与参数设置 (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梁结构经典计算 一榀轻钢结构库房框架，结构钢方管构件，材质E=210GPa，μ=， ρ=7850kg/m3（在不计重力的静力学分 析中可以不要）。F=1000N，此题要计入重力。计算水平梁中点下降位移。 文件与路径 顶部下拉菜单File, Save As ExpAbq02。 一部件 创建部件，命名为Prat-1。 3D，可变形模型，线，图形大约范围20(m)。 选用折线绘出整个图示屋架。 退出Done。 二性质 1 创建截面几何形状：Module，Property，Create Profile， 将截面(1)命名为Profile-1，选Box型截面，按图输入数据，关闭。直至完成截面(3)。 2 定义各段梁的方向： 选中所有立杆，输入截面主惯性轴1方向单位矢量(1,0,0)，选中横梁和斜杆，输入截面主轴1方向单位矢量(0,1,0)，关闭。还有好办法，请大家自己捉摸。

3 定义截面力学性质： 将截面(1) Profile-1命名为Section-1，梁，梁，截面几何形状选 Profile-1，输入E=210GPa，G=，ν=，ρ=7850，关闭。直至完成截面(3) Section-3。 4 将截面的几何、力学性质附加到部件上： 选中左右立柱和横梁，将各Section-1~3信息注入Part-1的各个杆件上，要对号入座。 5 保存模型： 将本题的CAE模型保存为。 三组装 创建计算实体，以Prat-1为原形，用Independent方式或Dependent生成实体。 四分析步 创建分析步，命名为Step-1，静态Static，通用General。 注释：无，时间：不变，非线性开关：关。 五载荷 1 施加位移边界条件： 命名为BC-1，在分析步Step-1中，性质：力学，针对位移和转角。 选中立柱两脚，约束全部自由度。 2 创建载荷： 命名为Load-1，在分析步Step-1中，性质：力学，选择集中力。 选中顶点，施加Fy=F2=-1000(N)。 六网格 对部件Prat-1进行。 1 撒种子： 针对部件，全局种子大约间距。 2 划网格： 针对部件，OK。 3 保存你的模型： 将本题的CAE模型保存为。

(完整版)Abaqus分析实例(梁单元计算简支梁的挠度)精讲

Abaqus分析实例（梁单元计算简支梁的挠度）精讲 对于梁的分析可以使用梁单元、壳单元或是固体单元。Abaqus的梁单元需要设定线的方向，用选中所需要的线后，输入该线梁截面的主轴1方向单位矢量(x,y,z)，截面的主轴方向在截面Profile设定中有规定。 注意： 因为ABAQUS软件没有UNDO功能，在建模过程中，应不时地将本题的CAE模型（阶段结果）保存，以免丢失已完成的工作。 简支梁，三点弯曲，工字钢构件，结构钢材质，E=210GPa，μ=0.28，ρ=7850kg/m3（在不计重力的静力学分析中可以不要）。F=10k N，不计重力。计算中点挠度，两端转角。理论解：I=2.239×10-5m4，w中=2.769×10-3m，θ边=2.077×10-3。 文件与路径： 顶部下拉菜单File, Save As ExpAbq00。 一部件 1 创建部件：Module，Part，Create Part， 命名为Prat-1；3D，可变形模型，线，图形大约范围10(程序默认长度单位为m)。 2 绘模型图：选用折线，从(0,0)→(2,0)→(4,0)绘出梁的轴线。 3 退出：Done。 二性质 1 创建截面几何形状：Module，Property，Create Profile， 命名为Profile-1，选I型截面，按图输入数据，l=0.1，h=0.2，b l=0.1，b2=0.1，t l=0.01，t2=0.01，t3=0.01，关闭。 2 定义梁方向：Module，Property，Assign Beam Orientation，

选中两段线段，输入主轴1方向单位矢量(0,0,1)或(0,0,-1)，关闭。 3 定义截面力学性质：Module，Property，Create Section， 命名为Section-1，梁，梁，截面几何形状选Profile-1，输入E=210e9（程序默认单位为N/m2，GPa=109N/m2）， G=82.03e9，ν=0.28，关闭。 4 将截面的几何、力学性质附加到部件上：Module，Property，Assign Section， 选中两段线段，将Section-1信息注入Part-1。 三组装 创建计算实体：Module，Assembly，顶部下拉菜单Instance，Create， Create Instance，以Prat-1为原形，用Independent方式生成实体。 四分析步 创建分析步：Module，Step， Create Step，命名为Step-1，静态Static，通用General。注释：无，时间：不变，非线性 开关：关。 五载荷 1 施加位移边界条件：Module，Load，Create Boundary Condition， 命名为BC-1，在分析步Step-1中，性质：力学，针对位移和转角，Continue。选中梁左端，Done，约束u1、u2、u3、u R1、u R2各自由度。 命名为BC-2，在分析步Step-1中，性质：力学，针对位移和转角，Continue。选中梁右端，Done，约束u2、u3、u R1、u R2各自由度。 2 创建载荷：Module，Load，Create Load， 命名为Load-1，在分析步Step-1中，性质：力学，选择集中力Concentrated Force，Continue。选中梁中点，Done，施加F y(CF2)=-10000（程序默认单位为N）。 六网格 对实体Instance进行。 1 撒种子：Module，Mesh，顶部下拉菜单Seed，Instance， Global Seeds，Approximate g lobal size 0.2全局种子大约间距0.2。 2 划网格：Module，Mesh，顶部下拉菜单Mesh，Instance，yes。 七建立项目 1 建立项目：Module，Job，Create Job，Instance，

三种混凝土本构模型

ABAQUS中的三种混凝土本构模型 2010-05-12 22:19:14| 分类：ABAQUS | 标签：|字号大中小订阅 资料来自SIMWE论坛shanhuimin923，特表示感谢！ ABAQUS 用连续介质的方法建立描述混凝土模型不采用宏观离散裂纹的方法描述裂纹 的水平的在每一个积分点上单独计算其中。 低压力混凝土的本构关系包括： Concrete Smeared cracking model (ABAQUS/Standard) Concrete Brittle cracking model (ABAQUS/Explicit) Concrete Damage plasticity model 高压力混凝土的本构关系： Cap model 1、ABAQUS/Standard中的弥散裂缝模型Concrete Smeared cracking model (ABAQUS/Standard)： ——只能用于ABAQUS/Standard中 裂纹是影响材料行为的最关键因素，它将导致开裂以及开裂后的材料的各向异性用于描述：单调应变、在材料中表现出拉伸裂纹或者压缩时破碎的行为 在进行参数定义式的Keywords： *CONCRETE *TENSION STIFFENING *SHEAR RETENTION *FAILURE RATIOS 2、ABAQUS/Explicit中脆性破裂模型Concrete Brittle cracking model (ABAQUS/Explicit) ：适用于拉伸裂纹控制材料行为的应用或压缩失效不重要，此模型考虑了由于裂纹引起的材料各向异性性质，材料压缩的行为假定为线弹性，脆性断裂准则可以使得材料在拉伸应力过大 时失效。 在进行参数定义式的Keywords *BRITTLE CRACKING, *BRITTLE FAILURE, *BRITTLE SHEAR 3、塑性损伤模型Concrete Damage plasticity model： 适用于混凝土的各种荷载分析，单调应变，循环荷载，动力载荷，包含拉伸开裂(cracking)和压缩破碎(crushing)，此模型可以模拟硬度退化机制以及反向加载刚度恢复的混凝土力学 特性 在进行参数定义式的Keywords： *CONCRETE DAMAGED PLASTICITY *CONCRETE TENSION STIFFENING *CONCRETE COMPRESSION HARDENING *CONCRETE TENSION DAMAGE *CONCRETE COMPRESSION DAMAGE

悬臂梁—有限元ABAQUS线性静力学分析实例

线性静力学分析实例——以悬臂梁为例 线性静力学问题是简单且常见的有限元分析类型，不涉及任何非线性（材料非线性、几何非线性、接触等），也不考虑惯性及时间相关的材料属性。在ABAQUS中，该类问题通常采用静态通用（Static，General）分析步或静态线性摄动（Static，Linear perturbation）分析步进行分析。 线性静力学问题很容易求解，往往用户更关系的是计算效率和求解效率，希望在获得较高精度的前提下尽量缩短计算时间，特别是大型模型。这主要取决于网格的划分，包括种子的设置、网格控制和单元类型的选取。在一般的分析中，应尽量选用精度和效率都较高的二次四边形/六面体单元，在主要的分析部位设置较密的种子；若主要分析部位的网格没有大的扭曲，使用非协调单元（如CPS4I、C3D8I）的性价比很高。对于复杂模型，可以采用分割模型的方法划分二次四边形/六面体单元；有时分割过程过于繁琐，用户可以采用精度较高的二次三角形/四面体单元进行网格划分。 悬臂梁的线性静力学分析 1.1 问题的描述 一悬臂梁左端受固定约束，右端自由，结构尺寸如图1-1所示，求梁受载后的Mises应力、位移分布。 ν 材料性质：弹性模量3 = E=，泊松比3.0 2e 均布载荷：F=103N 图1-1 悬臂梁受均布载荷图 1.2 启动ABAQUS 启动ABAQUS有两种方法，用户可以任选一种。 （1）在Windows操作系统中单击“开始”--“程序”--ABAQUS 6.10 --

ABAQUS/CAE。 （2）在操作系统的DOS窗口中输入命令：abaqus cae。 启动ABAQUS/CAE后，在出现的Start Section（开始任务）对话框中选择Create Model Database。 1.3 创建部件 在ABAQUS/CAE顶部的环境栏中，可以看到模块列表：Module：Part，这表示当前处在Part（部件）模块，在这个模块中可以定义模型各部分的几何形体。可以参照下面步骤创建悬臂梁的几何模型。 （1）创建部件。对于如图1-1所示的悬臂梁模型，可以先画出梁结构的二维截面（矩形），再通过拉伸得到。 单击左侧工具区中的（Create Part）按钮，或者在主菜单里面选择Part--Create，弹出如图1-2所示的Create Part对话框。 图1-2 Create Part对话框 在Name（部件名称）后面输入Beam，Modeling Space（模型所在空间）设

ABAQUS(Explict)混凝土开裂模型翻译

混凝土开裂模型 适用模块：Abaqus/Explicit Abaqus/CAE 参考 ●“Material library: overview,” Section 18.1.1 ●“Inelastic behavior,” Section 20.1.1 ●*BRITTLE CRACKING ●*BRITTLE FAILURE ●*BRITTLE SHEAR ●“Defining brittle cracking” in “Defining other mechanical models,” Section 12.9.4 of the Abaqus/CAE User's Manual 概述 Abaqus/Explicit模块中脆性断裂模型： ●提供一种通用模型来模拟包括梁单元、桁架单元、壳单元以及实体单元在内的所有单元 形式； ●也可以用来模拟诸如陶瓷及脆性岩石的其他材料； ●用于模拟受拉开裂占主导地位的材料本构行为； ●假设受压行为是线弹性的； ●必须与线弹性模型(“Linear elastic behavior,” Section 19.2.1)同时使用，它也 定义了材料开裂前的本构行为； ●用于模拟脆性行为占主导地位的本构关系是十分准确的，基于此，假设受压行为是线弹 性的是合理的； ●该模型主要是用于钢筋混凝土结构的分析，同时也适用于素混凝土； ●基于脆性失效准则，将失效单元删除； 关于失效单元删除的内容详见“A cracking model for concrete and other brittle materials,” Section 4.5.3 of the Abaqus Theory Manual. 关于ABAQUS中可用混凝土本构模型的相关讨论参见“Inelastic behavior,” Section 20.1.1。 钢筋 ABAQUS中，混凝土结构中的钢筋是通过指定Rebar单元实现的。Rebar单元是一维应变理论单元（杆单元），既可以单独定义，也可以镶嵌在有向曲面上。关于Rebar单元的相关讨论参见“Defining rebar as an element property,” Section 2.2.4。Rebar单元特别地用来模拟钢筋的弹塑性行为，并且可叠加在用于模拟素混凝土标准单元的网格上。基于这种模拟方法，混凝土的开裂行为与Rebar是没有关系的。混凝土和钢筋之间的相互作用，例如粘结滑移、销栓作用，可以引入拉神硬化（强化）的概念来近似模拟混凝土裂缝处荷载向钢筋转移的特点。 开裂 Abaqus/Explicit中使用弥散裂缝模型来表征混凝土脆性行为的非连续性。这种方法并不关注于单个宏观裂缝，相反地，只是独立地计算有限元模型质点处的本构关系。裂缝的存在对于计算的影响在于：裂缝的存在影响质点处的应力和材料刚度。 为简化本部分讨论内容，“开裂”一词实质上代表的是一个方向——所考虑单个质点处裂缝的方向，与其最相近的物理解释为：在一个质点附近出现一系列连续的微裂缝，其方向