编者按:
熟悉ANSYS Mechanical的朋友知道,早年ANSYS经典界面风行一时,后来随着2000年后ANSYS Workbench平台的推出,经过十多年的发展完善,其易用性、功能性进入了一个非常强健稳定的状态,现在用Workbench平台进行分析的工程师越来越多,毋容置疑的是其易用性远超ANSYS经典界面,在功能角度也实现了相当的水平。早年学习ANSYS的朋友会使用一些经典的练习案例进行学习,熟悉软件的操作及基本特性,那这些经典案例是非常有学习意义的,不过这些官方的经典案例并没有Workbench的版本,所以南京安世亚太集中资源对一些经典的ANSYS学习算例进行了梳理,在workbench中进行了一些复现的尝试,并将以连载的方式与爱好者们分享,希望能对大家的学习工作有所帮助。
案例六:某钢筋混凝土结构分析
一工程背景
钢筋混凝土结构在建筑和部分机械结构中经常被用到,钢筋布置不合理会使得结构在载荷影响下(例如恒定静载,冲击载荷,地震载荷等)发生脆性断裂。所以,对于混凝土和钢筋的力学性能研究非常重要。本案例以某简化的跨海桥桥墩结构为例,介绍如何在 ANSYS Workbench环境中定义混凝土材料,同时在ANSYS Workbench环境中如何考虑加强筋进行结构仿真。本案例中的结构并非实际真实结构,但是,通过该案例可以让广大用户了解如何在ANSYS Workbench环境中定义高级材料非线性模型中的微平面材料模型,同时如何在ANSYS Workbench环境中定义结构加强筋。二 问题描述
本案例中的结构模型取材自某跨海大桥的桥墩。所有尺寸均非实际尺寸,且不考虑跨海大桥的拉索结构,假设所有载荷全部施加桥墩横梁上。同时,桥墩横梁上布置了上下两层钢筋,每层九条,每条间距0.5m。仿真中,钢筋材料定义为默认结构钢材料:弹性模量2e5MPa;泊松比为0.3。图1 几何模型
图2 加强筋布置
其中,混凝土材料采用ANSYS中的微平面材料模型,加强筋采用REINF264单元。常规实体单元和加强筋单元在节点位置处连接。 三 边界条件和载荷 本案例主要探究加强筋对混凝土桥墩结构变形的影响。将桥墩底部固定,在横梁上表面施加100000N的沿重力方向的载荷。 四 微平面材料模型简介 微平面材料模型是由Bazant和Gambarova研究提出的,由于结构刚度的衰减退化与方向有关,所以每一个平面都有潜在的失效风险,并导致宏观上的结构的破坏。该材料模型认为在每个平面上,材料都满足应力-应变法则。通过对每个微平面的评估来考虑这种与方向相关的刚度变化。这种材料模型适合模拟由不同属性的成分聚集而成的工程材料,比如混凝土材料是由沙石填充到水泥中聚合而成的。 微平面理论包含以下三个重要步骤: 1 通过一个动态约束将宏观应变张量和其相应的微平面进行关联(映射); 2 在微平面定义连续性法则,在每一个微平面定义单轴连续方程(例如应力和应变分量); 3 基于能量守恒准则对整体材料进行均化,得到材料的整体响应。 微平面材料模型假设在每一个微平面都存在微观自由能,且所有微平面上的微观自由能的积分等于宏观Helmholtz能量:![7687132ae40986edba126d4d5c018c5b.png](https://i-blog.csdnimg.cn/blog_migrate/cef356c213e1e90b2667f696290a80e7.png)
下图所示为相同损伤发展速率下,不同最大损伤的应力应变曲线
图3 相同损伤发展速率下,不同最大损伤的应力应变曲线
下图所示为相同最大损伤下,不同损伤发展速率的应力应变曲线
图4 相同最大损伤下,不同损伤发展速率的应力应变曲线
通过以上两张图,可以很好得理解和是如何影响材料特性的。 等效应变 用于描述损伤发展准则,其公式如下:图5 Discrete加强筋模型和Smeared加强筋模型
ANSYS中建立加强筋单元需要两个步骤: 1 用网格独立法或标准法定义加强筋模型的材料、几何和对齐方式; 2 建立加强筋单元并嵌入至基础单元里。 网格独立法对基础单元的形状没有特定的要求,这为加强筋单元的建立提供了更多的可能性与灵活性。其具体思路是先用Mesh200单元代替加强筋单元,对几何进行网格划分,然后再使用加强筋单元替代Mesh200单元。诸如加强筋材料、加强筋横截面以及加强筋对齐方式等单元信息,用户可以在加强筋单元横截面形式(SECDATA)中进行定义。 标准法建立加强筋单元对于基础网格为梁单元、管单元和分层单元等单元类型较为便利。所有的加强筋单元的单元属性都在加强筋单元横截面形式中进行定义,而加强筋的位置由基础单元网格的节点位置控制。 本案例中采用网格独立法建立加强筋单元。按照上述理论,分别在Workbench Mechanical(本案例采用静力学分析)的Geometry和Static Structural中插入命令: 1 针对geometry的所有线体,插入如下命令: sectype,2,reinf,discrete !该线体对应的材料号为2,为该线体定义加强筋横截面 secdata,2,2827.3,MESH !2827.3 为加强筋的横截面积,mesh 表示由mesh200 单元确定加强筋的位置 et,2,200,2 keyopt,2,1,3 secnum,2 上述命令的含义是用MESH200单元对线体进行网格划分,并将加强筋横截面属性赋予对应的MESH200单元。 2在Static Structural中插入如下命令: Ereinf !建立加强筋单元 需要指出的是,加强筋单元无法通过命令选择单元类型的方式定义,只能通过该方法生成。此外,该命令并没有生成新的节点,加强筋单元和基础单元共享节点信息。六 计算结果
如下图所示,不考虑加强筋的结构最大变形为0.21mm,最大拉应力为0.08MPa;考虑加强筋后结构最大变形为0.14mm,最大拉应力为0.05MPa。不考虑加强筋结构变形云图
考虑加强筋结构变形云
图6 结构变形结果
不考虑加强筋结构最大主应力云图
考虑加强筋结构最大主应力云图
图7 结构最大主应力结果
由于WorkbenchMechanical中暂时不支持直接输出加强筋的结果,用户可以通过插入命令流的方式,或者直接将结果导入ANSYS经典截面查看加强筋的应力分布。注意,此时线体几何已经不是由Beam单元进行网格划分的,所以Beam Tool在此无法输出加强筋的结果。图8 加强筋等效应力云图
我们 | 南京安世亚太
▼