PQ-Fiber

概述

PQ-Fiber是清华大学土木工程系结构工程研究所潘鹏副教授和曲哲博士基于大型通用有限元程序ABAQUS开发的一组材料单轴滞回本构模型的集合。主要用于在钢筋混凝土结构、钢结构等的弹塑性时程分析中定义杆系结构的材料本构,同时可用于任何需要定义材料单轴滞回本构模型的情况。

作者以FORTRAN编译文件.obj的形式在网上免费发布PQ-Fiber的最新版本。为适用于不同的系统平台和需求,PQ-Fiber发布如下版本:

PQ-Fiber_Std20_32:适用于32位系统平台和ABAQUS/Standard

PQ-Fiber_Std20_64:适用于64位系统平台和ABAQUS/Standard

作者提供上述版本以供广大科研与工程设计人员使用,发布的版本没有功能限制,请使用者尊重知识产权。

PQFiber版本信息:V2.0

推荐平台:ABAQUS 6.9.1Visual Studio 2005Intel Fortran 9.1

包含的材料模型(详细介绍见第三节):

       UConcrete01UConcrete02Uconcrete03, USteel01USteel02USteel03

使用过程中如有问题,请与作者联系:

潘鹏(Email: panpeng@tsinghua.edu.cn

通信地址:北京清华大学土木工程系,100084

 

版本信息:v2.0 下载-单击右键另存为
         

包含的材料模型(详细介绍见第三节):

       UConcrete01UConcrete02USteel01USteel02USteel03

使ABAQUS示例inp文件下载
【安装ABAQUSFortran成功后,把PQFiber.obj和示例文件放在同一文件夹下,双击abarun.bat

 


 

相关下载:

陆新征-曲哲塑性铰模型

MARC 自定义塑性铰示例源程序

MSC.MARC THUFIBER纤维模型程序

ATC-63报告建议的远场地震波库

在MARC 2005R2中加入简单钢筋混凝土纤维模型源程序

在MARC中实现基于多点位移控制的推覆分析算法源程序

 

ABAQUS中使用

1. ABAQUS中使用本模型

1)定义材料

Properties模块中定义User Material,如图1所示。材料名的前几个字母必须与第三节中定义的某一个材料名相一致。需要分别选择General选项卡中的User MaterialDepvar两个选项。

User Material选项中定义该材料所需要的所有材料属性,如图2所示。在Depvar选项中定义该材料所需的状态变量的个数,如图3所示。

也可以在.inp文件中直接添加用户自定义材料,下面给出了一个例子。

*Material, name=UConcrete01

*Depvar

      5,

*User Material, constants=4

30., 0.002,   10., 0.005

1 定义用户自定义材料

 

2 定义材料属性

 

3 定义状态变量的个数

注意:

如果使用ABAQUS中的铁木辛科梁单元(如B21B22B31B32等),除了上述定义外,还需要在.inp文件为梁单元截面额外定义横向剪切刚度,定义方法详见ABAQUS keyword reference manual中的*Transverse Shear Stiffness关键字。可以直接在.inp文件中添加,也可以利用ABAQUS/CAE提供的Keyword editor,如图4所示。剪切刚度可以是一个大数,它对计算结果的影响不大。

 

4 Keyword editor中为截面定义剪切刚度

 

2)调用用户子程序

Job模块中,在Edit Job对话框的General选顶卡中选择本文提供的.obj文件作为User subroutine file,如图5所示。

由于ABAQUS的一个分析模型只能接受一个用户子程序文件,如果使用了本文编译后的材料模型文件,将无法使用其它用户子程序。用户如确有需要,可与本模型作者联系。

 

5 Job中调用用户子程序

 

2. ABAQUS/Standard中定义钢筋混凝土梁单元

可以使用*rebar关键字在梁单元中定义钢筋,使用方法详见ABAQUS keyword reference manual。下面给出了定义钢筋混凝土梁单元一个例子,其中分别用UConcrete02USteel02来定义梁单元中的混凝土和钢筋的材料本构,并在梁截面的四角分别定义了四根钢筋。

 

*Beam Section, elset=Pier, material=UCONCRETE02, section=RECT

300., 300.

0.,0.,-1.

*Transverse shear stiffness

 1.0e16, 1.0e16

*rebar, element=beam, material=USTEEL02-235, name=rebar1

Pier, 201, 107, 107

*rebar, element=beam, material=USTEEL02-235, name=rebar2

Pier,201,-107,107

*rebar, element=beam, material=USTEEL02-235, name=rebar3

Pier,201,107,-107

*rebar, element=beam, material=USTEEL02-235, name=rebar4

Pier,201,-107,-107

 

在同一模型中可以使用不同自定义钢筋与混凝土模型,也可以采用同一种自定义钢筋或混凝土模型定义不同等级的钢筋或者混凝土。如在图6所示的Material Manager中,用USteel02模型定义了屈服强度分别为235MPa335MPa的两种钢筋,用UConcrete02模型定义了轴压强度分别为30MPa60MPa的两种混凝土。只要材料名的前几个字母与第三节中定义的材料名相同即可完成调用,材料名的后半部分可以随意定义。

6 在一个模型中使用不同的自定义材料

 

3. 结果输出

自定义材料结果的输入通过状态变量实现。如果想得到弹性应变以外的输出,用户需要在Step模块中定义对状态变量(SDV)的输出,如图7所示,或者在.inp文件中直接添加对SDV的输出,如下所示。另外,ABAQUS目前不支持梁单元中*rebar定义的钢筋的任何输出。各材料模型中状态变量的物理意义详见第三节。

 

*Element Output, directions=YES

E, S, SDV, SE, SF

 

7 输出状态变量SDV

 

材料模型

PQ-Fiber v1.2提供如下材料滞回模型。

调用名

USteel01

描述

弹塑性随动硬化单轴本构模型。

材料参数

props(1)

弹性模量

props(2)

屈服强度

props(3)

硬化刚度系数,等于第二阶段刚度与弹性模量之比。

状态变量

SDV(1)

塑性应变

SDV(2)

Back stress

图例

 

往复加载时的单轴应力应变关系

 


 


调用名

USteel02

描述

再加载刚度按Clough本构退化的随动硬化单轴本构模型

材料参数

props(1)

弹性模量

props(2)

屈服强度

props(3)

硬化刚度系数,等于第二阶段刚度与弹性模量之比。

状态变量

SDV(1)

历史最大拉应变

SDV(2)

对应于历史最大拉应变的应力

SDV(3)

历史最大压应变

SDV(4)

对应于历史最大压应变的应力

SDV(5)

屈服记号

图例

 

往复加载时的单轴应力应变关系

 


 

调用名

USteel03

描述

拉压不等强的弹塑性随动硬化单轴本构模型

材料参数

props(1)

弹性模量

props(2)

受拉时的屈服强度

props(3)

受压时的屈服强度

props(4)

硬化刚度系数,等于第二阶段刚度与弹性模量之比。

状态变量

SDV(1)

塑性应变

SDV(2)

Back stress

图例

 

往复加载时的单轴应力应变关系

 


 

调用名

UConcrete01

描述

忽略抗拉强度的混凝土模型。

(1) 拉应力始终为零;

(1) 受压骨架线上升段采用Hognested曲线,下降段为直线;

(2) 受卸载刚度随历史最大压应变的增大而减小;

(3) 受拉后反向加载时,直至达到上次受压卸载的残余应变时再开始反向承载。

材料参数

props(1)

轴心受压强度

props(2)

峰值压应变,即达到轴心受压强度时的应变

props(3)

极限受压强度

props(4)

极限压应变

props(5)

截面钢筋屈服的临界应变,定义为混凝土受拉边缘的应变。

状态变量

SDV(1)

初始化变量,无实际意义

SDV(2)

历史最大压应变

SDV(3)

受压残余应变,即卸载至应力为零时的压应变

SDV(4)

卸载/再加载刚度

SDV(5)

截面屈服标志

图例

 

往复加载时的单轴应力应变关系

图例中fu = props(3)eu = props(4)

 


 

调用名

UConcrete02

描述

考虑抗拉强度的混凝土模型。

(1) 受拉骨架线由线性上升段和线性下降段组合;

(2) 受拉卸载时指向原点;

(3) 受压骨架线上升段采用Hognested曲线,下降段为直线;

(4) 受卸载刚度随历史最大压应变的增大而减小,且不小于达到极限压应变时的卸载刚度。

材料参数

props(1)

轴心受压强度

props(2)

峰值压应变,即达到轴心受压强度时的应变

props(3)

极限受压强度

props(4)

极限压应变

props(5)

达到极限压应变时的卸载刚度与初始弹性模量之比

props(6)

轴心受拉强度

props(7)

受拉软化模量,即受拉骨架线下降段的刚度,输入其绝对值即可。

props(8)

截面钢筋屈服的临界应变,定义为混凝土受拉边缘的应变。

状态变量

SDV(1)

初始化变量,无实际意义

SDV(2)

历史最大压应变

SDV(3)

受压残余应变,即卸载至应力为零时的压应变

SDV(4)

卸载/再加载刚度

SDV(5)

截面屈服标志

图例

 

往复加载时的单轴应力应变关系

 

个人信息

研究工作

实际工程

论文工作

教学工作

资料下载

专题

其他

 

我们的实验室

抗倒塌专业委员会