|
编者按：
熟悉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能量：
为了考虑材料退化和破坏，通过参数对微观自由能方程进行修正，和分别表示微平面应变的体积项和偏离项，由微平面应变拆分而来，由宏观概念的应变通过投影的方式转化而来。
下图所示为相同损伤发展速率下，不同最大损伤的应力应变曲线 
图3 相同损伤发展速率下，不同最大损伤的应力应变曲线 下图所示为相同最大损伤下，不同损伤发展速率的应力应变曲线 
图4 相同最大损伤下，不同损伤发展速率的应力应变曲线 通过以上两张图，可以很好得理解和是如何影响材料特性的。
注意：这里的等效应变 特指微平面材料模型中的上述方程式，与ANSYS Mechanical中的等效应变不是同一个概念。 在ANSYS中，可以通过命令流定义微平面材料模型，具体命令如下：tbdata,1,k0,k1,k2,gamma,alpha,beta在ANSYS中，加强筋指纤维长度远大于其横截面的纤维。对于结构分析而言，一般梁单元的弯曲、扭转和剪切刚度在加强筋单元中均不予考虑，故指存在其轴向刚度。ANSYS将加强筋单元与其基础单元进行绑定，故在整个计算过程中，二者的相对位移也是不存在的。换言之，加强筋单元的运动，取决于也仅取决于其基础单元的运动。在ANSYS中，存在两种加强筋模型：discrete和smeared。Discrete选项用于建立材料、横截面以及对齐方向不一致的加强筋。程序分别将不同的加强筋处理成不同的只考虑单轴刚度的梁。Smeared选项用于以层形式的加强筋。每一层中，单向刚度的加强筋的材料、横截面以及对齐方式均相同。 
图5 Discrete加强筋模型和Smeared加强筋模型 1 用网格独立法或标准法定义加强筋模型的材料、几何和对齐方式；网格独立法对基础单元的形状没有特定的要求，这为加强筋单元的建立提供了更多的可能性与灵活性。其具体思路是先用Mesh200单元代替加强筋单元，对几何进行网格划分，然后再使用加强筋单元替代Mesh200单元。诸如加强筋材料、加强筋横截面以及加强筋对齐方式等单元信息，用户可以在加强筋单元横截面形式（SECDATA）中进行定义。标准法建立加强筋单元对于基础网格为梁单元、管单元和分层单元等单元类型较为便利。所有的加强筋单元的单元属性都在加强筋单元横截面形式中进行定义，而加强筋的位置由基础单元网格的节点位置控制。本案例中采用网格独立法建立加强筋单元。按照上述理论，分别在Workbench Mechanical（本案例采用静力学分析）的Geometry和Static Structural中插入命令：1 针对geometry的所有线体，插入如下命令：sectype,2,reinf,discrete !该线体对应的材料号为2，为该线体定义加强筋横截面 secdata,2,2827.3,MESH !2827.3为加强筋的横截面积，mesh表示由mesh200单元确定加强筋的位置上述命令的含义是用MESH200单元对线体进行网格划分，并将加强筋横截面属性赋予对应的MESH200单元。2在Static Structural中插入如下命令：需要指出的是，加强筋单元无法通过命令选择单元类型的方式定义，只能通过该方法生成。此外，该命令并没有生成新的节点，加强筋单元和基础单元共享节点信息。六 计算结果 如下图所示，不考虑加强筋的结构最大变形为0.21mm，最大拉应力为0.08MPa；考虑加强筋后结构最大变形为0.14mm，最大拉应力为0.05MPa。
不考虑加强筋结构变形云图
考虑加强筋结构变形云 图6 结构变形结果
不考虑加强筋结构最大主应力云图
考虑加强筋结构最大主应力云图 图7 结构最大主应力结果 由于WorkbenchMechanical中暂时不支持直接输出加强筋的结果，用户可以通过插入命令流的方式，或者直接将结果导入ANSYS经典截面查看加强筋的应力分布。注意，此时线体几何已经不是由Beam单元进行网格划分的，所以Beam Tool在此无法输出加强筋的结果。
图8 加强筋等效应力云图
|
