发布时间 : 星期二 文章abaqus有限元分析报告开裂梁 - 图文更新完毕开始阅读c0ffa340e718964bcf84b9d528ea81c758f52e7a
Abaqus梁的开裂模拟计算报告
1. 问题描述
利用ABAQUS有限元软件分析如图1.1所示的钢筋混凝土梁的裂缝开展。参考文献Brena et al.（2003）得到梁的基本数据：
图1.1 Brena et al.（2003）中梁C尺寸
几何尺寸：跨度3000mm，截面宽203mm，高406mm的钢筋混凝土梁 由文献Chen et al. 2011得材料特性：
1.混凝土：抗压强度fc’=35.1MPa，抗拉强度ft=2.721MPa，泊松比ν=0.2，弹性模量Ec=28020MPa；
2.钢筋：弹性模量为Ec=200GPa，屈服强度fys=fyc=440MPa，fyv=596MPa 3.混凝土垫块：弹性模量为Ec=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。
2) Property
Module中选择property，然后选择功能模块
对不同的材料进行赋值，下面是各种
材料输入时候的数据：
① 混凝土本构关系： 模型一：（该模型的尺寸单位为m）在弹出的对话框中命名为beam，在Mechanical选项中点击Elastic，在Young’s Modulus中输入28020000000，Poisson’s Ratio为0.2；类似的方法找到Concrete Damaged Plasticity，按混凝土结构设计规范对受压取点8个，受拉取点7个。下面是计算过程：
由规范中附录C中C2 混凝土本构关系： C.2.3混凝土单轴受拉的应力-应变曲线公式：
??(1?dt)Ec?，所用参数可以参考规范（混凝土结构设计规范GB50010-2010）
C.2.4混凝土单轴受压的应力-应变曲线公式：
??(1?dc)Ec?，所用参数可以参考规范
以及塑性应变与总应变的关系：
???el??pl，其中?el??E
以及塑性应变与总应变的关系：
???el??cr，?el与上式相同
借助matlab软件计算受压、拉时的本构关系方程，其中塑性应变分别取0、0.0005、0.001、0.002、0.003、0.004、0.005、0.008得对应应力为23171225.88、35072364.4、32320980.31、24236963.88、17763125.96、13340141.18、10334750.06、5892086.123；开裂应变分别取0、0.001、0.003、0.005、0.008、0.01、0.05得对应应力2721000、2648625.968、875761.085、511922.937、334706.391、50489.31、16333.482。
由于水平限制，算得的受拉损伤因子未能通过abaqus的算前检测，这里没有输入。 模型二：（该模型单元尺寸为m）Elastic中的设置与模型一相同，计算其本构关系的时候按照文献Finite-Element Modeling of Intermediate Crack Debonding in FRP-Plated RC Beams，G. M. Chen中的公式：
单轴受压：????1?[(??p/?p)?2](?/?p)?(?/?p)2，所用参数可见于引用的文献；
受压的时候为了寻找输入塑性应变的起始点，令abaqus中提供的应力--应变输入方式的名义值与真实值关系公式为0，有：
ln(1??nom)??nom(1??nom)E?0
?nom??nom?1?[(???/??)?2](?nom/??)?(?nom/??)2pl解之得：?nom?22.76652048MPa，?nom?0.00081284，?ln?0。
?nom从点0.00081284开始以步长0.0002得到83个点。在输入到abaqus前使用公式
pl?ln?ln(1??nom)??nom(1??nom)E
?true??nom(1??nom)进行转换，输入的具体值可见附录。
wtwt3(?c2wcr)wt单轴受拉： ?[1?(c1)]e?(1?c13)e(?c2)，所用参数可见于引用文献；
ftwcrwcr由于hc为裂缝带宽（上面提到的文献中也命名为平面四节点单元的特征裂缝长度），取
?t2e，其中e为单元网格长度，这里取10mm，即hc为14mm。
在清华江见鲸或者J.G. ROTS撰写的文献中都能找到其他单元类型的hc与e的换算公
式，这里提到的hc也可称为Lcr，这提供了把本文中的开裂位移转化为开裂应变的方法：
wt??crLcr
对wt进行离散是采用前密后疏的方法，一开始的10个点步长为0.00012，中间,28个点步长为0.0012，最后13个点步长为0.006，加上零点一共53个点。
把位移转换成应变以后，同样地，使用名义值--真实值的转换公式得到数据，具体数值可见于附录。由于水平问题，算得的受拉损伤因子未能通过abaqus的算前检测，这里没有输入。
模型三：该模型使用的尺寸单位为mm，受压时的取值只需在模型二的取值基础上进行单位换算即可。本构关系的输入方式为应力--位移的方法，在输入类型中选择Displacement
单轴受拉：对wt进行离散，取等步长0.0012，共51个点。这里还进行了受拉损伤因子dt的计算，按公式：dt?wt，具体数据可见附录。
[wt?(hc?t)/Ec]② 钢筋本构关系：
类似的，把钢筋的杨氏模量输入到Elastic中，三个模型的受拉受压钢筋都是200GPa，在Plastic中按照名义值--真实值得方法，取得两个点：屈服强度440MPa，塑性应变为0；屈服强度448.8MPa，塑性应变为0.02。箍筋的屈服强度取596MPa，塑性应变为0。
③ 垫块的本构关系：
这里垫块的本构可以按混凝土的本构输入，也可以按钢筋的输入，为了方便计算取钢筋的本构关系作为其材料属性。
建立完本构关系后需对混凝土等截面属性进行赋值，点击
，弹出Create Section对
话框，将Category设为solid，Type设为Homogeneous，其余参数保持默认，点击Continue，material为beam，thickness为203mm，垫块也是类似的输入方式；对于钢筋，Category设为beam，Type设为truss，如受拉钢筋，由文献得直径为16mm，面积为200.96（mm2），由于同一水平面上有两根，取值402；受压钢筋和箍筋分别按直径9.5mm和7mm计算。
然后再同一个环境栏中使用
，提示区要求用户选择赋予截面的部件，分别对上述
创建过的部件赋予材料属性，完成操作。
3) Assembly