ANSYS中混凝土的计算问题
基于ANSYS的大体积混凝土的水化热模拟研究
基于ANSYS的大体积混凝土的水化热模拟研究一、概览话说这大家伙儿混凝土,可是咱们建筑行业里头的顶梁柱呢!它不仅结实耐用,而且造型多样,满足了咱们各种建筑需求。
然而混凝土的诞生可不是一蹴而就的,它是经过无数科学家和工程师的努力研究、实验、改进才逐渐发展起来的。
这其中大体积混凝土作为一种特殊的混凝土形式,因其施工难度大、质量要求高等特点,一直是建筑工程领域的研究热点。
那么大体积混凝土的水化热问题又是个啥情况呢?别着急接下来咱就来详细说说这个话题。
1.1 研究背景和意义随着社会的发展和科技的进步,大体积混凝土在建筑领域的应用越来越广泛。
然而大体积混凝土的水化热问题一直是困扰工程界的一个难题。
水化热是指水泥与水反应产生的热量,这种热量在一定程度上会影响混凝土的性能和使用寿命。
因此研究大体积混凝土的水化热规律,对于提高混凝土结构的抗裂性、耐久性和安全性具有重要意义。
ANSYS是一款广泛应用于工程领域的有限元分析软件,可以模拟各种物理现象和过程。
利用ANSYS对大体积混凝土的水化热进行模拟研究,可以更直观地了解其内部发生的热力学过程,为实际工程提供有力的理论支持。
此外这种方法还可以避免因现场条件限制而导致的实际试验结果与理论预测之间的误差,提高工程质量。
1.2 国内外研究现状嗨,亲爱的读者朋友们!让我来和大家分享一下关于大体积混凝土的水化热模拟研究的最新进展。
你知道吗这个课题在国内外的研究现状中一直备受关注,研究成果也是五花八门,各有千秋。
首先让我们看看国内的研究现状,近年来随着科技的发展和社会需求的变化,大体积混凝土的应用越来越广泛。
然而其水化热效应如何,能否通过模拟进行预测,一直是困扰我们的难题。
国内的一些学者们对此进行了深入研究,提出了一系列的观点和方法,为我们理解和利用大体积混凝土提供了新的视角。
然后我们再来看看国外的研究情况,由于文化背景和科研环境的不同,国外的研究风格和方法也有所不同。
一些国外的研究者更倾向于直接实验验证,他们设计了各种实验方案,通过对比实验结果来分析大体积混凝土的水化热效应。
ANSYS计算大体积混凝土温度场的关键技术
ANSYS计算大体积混凝土温度场的关键技术在大型混凝土结构建设的过程中,温度场分析对保障混凝土结构的安全性至关重要。
ANSYS作为工程领域中常用的数值模拟工具,能够对混凝土温度场进行准确的计算,为混凝土结构的设计和施工提供科学依据。
但是,对于大体积混凝土的温度场计算,存在一些关键技术需要考虑,下面将进行详细介绍。
1. 混凝土的物理性质混凝土的物理性质是温度场模拟中的关键因素之一。
混凝土在浇筑后的初凝期、成型期、硬化期、老化期等各个阶段的物理性质都存在巨大的变化。
因此,在进行混凝土的温度场计算前,需要准确地测量混凝土在不同时间点的物理性质,如热导率、比热容、密度等。
2. 热源的模拟混凝土的温度场计算需要考虑混凝土内部的各种热源对温度场的影响。
建筑中的热源包括太阳辐射、室内外温度差、人体热辐射等,需要对这些热源进行准确的模拟。
3. 初始条件的设置混凝土温度场计算的初始条件设置直接影响计算结果的准确性。
混凝土在浇筑后的初始温度值、初始变形状态的设置等都需要进行准确、合理的处理。
4. 传热模型的选择对于大体积混凝土的温度场计算,需要选择合适的传热模型。
传热模型可以根据混凝土的物理性质和热源的模拟情况,选择适用于不同情况下的传热模型,如动态传热模型、静态传热模型等。
5. 计算方法的选择针对大体积混凝土温度场的计算,需要选择合适的计算方法。
常用的方法有有限元法、有限差分法等,需要根据混凝土内部温度场、变形场等的变化情况,选择合适的计算方法。
6. 数值模拟使用ANSYS进行混凝土温度场计算,需要进行数值模拟。
数值模拟是对真实物理系统的数学模拟,通过建立数学模型,利用计算机运算获得物理系统的各种行为特性,如温度场、应力场、变形场等。
7. 计算结果的验证在进行混凝土温度场计算后,需要对计算结果进行验证。
验证结果通常采用实验测试的方式进行验证,如温度测试、原位应力测试、变形测试等。
,针对大体积混凝土温度场计算,需要考虑混凝土的物理性质、热源模拟、初始条件设置、传热模型选择、计算方法选择、数值模拟和计算结果验证等方面的关键技术,以保证计算结果的准确性和可靠性。
ANSYS分析钢筋混凝土结构技巧及实例详解
0 前言利用ANSYS分析钢筋混凝土结构时,其有限元模型主要有分离式和整体式两种模型。
这里结合钢筋混凝土材料的工作特性,从模型建立到非线性计算再到结果分析的全过程讲述了利用ANSYS进行钢筋混凝土结构分析的方法与技巧,并以钢筋混凝土简支梁为例,采用分离式有限元模型,说明其具体应用。
1 单元选取与材料性质1. 1 混凝土单元ANSYS中提供了上百种计算单元类型,其中Solid65单元是专门用于模拟混凝土材料的三维实体单元。
该单元是八节点六面体单元,每个节点具有三个方向的自由度( UX , UY , UZ) 。
在普通八节点线弹性单元Solid45 的基础上,该单元增加了针对于混凝土的材性参数和组合式钢筋模型,可以综合考虑包括塑性和徐变引起的材料非线性、大位移引起的几何非线性、混凝土开裂和压碎引起的非线性等多种混凝土的材料特性。
使用Solid65 单元时,一般需要为其提供如下数据:1)、实常数(Real Constants) :定义弥散在混凝土中的最多三种钢筋的材料属性,配筋率和配筋角度。
对于墙板等配筋较密集且均匀的构件,一般使用这种整体式钢筋混凝土模型。
如果采用分离式配筋,那么此处则不需要填写钢筋实常数。
2)、材料模型(Material Model) :在输入钢筋和混凝土的非线性材料属性之前,首先必须定义钢筋和混凝土材料在线弹性阶段分析所需的基本材料信息,如:弹性模量,泊松比和密度。
3)、数据表(Data Table) :利用数据表进一步定义钢筋和混凝土的本构关系。
对于钢筋材料,一般只需要给定一个应力应变关系的数据表就可以了,譬如双折线等强硬化(bilinear isotropic hardening)或随动硬化模型( kinematic hardening plasticity)等。
而对于混凝土模型,除需要定义混凝土的本构关系外,还需要定义混凝土材料的破坏准则。
在ANSYS中,常用于定义混凝土本构关系的模型有:1)多线性等效强化模型(Multilinear isotropic hardening plas2ticity ,MISO模型),MISO模型可包括20条不同温度曲线,每条曲线可以有最多100个不同的应力-应变点;2)多线性随动强化模型(Multilinear kinematic hardening plas2ticity ,MKIN 模型),MKIN 模型最多允许5个应力-应变数据点;3)Drucker2Prager plasticity(DP)模型。
ansys 钢筋混凝土建模
ansys 钢筋混凝土建模Ansys 钢筋混凝土建模在现代工程领域中,钢筋混凝土结构的应用极为广泛,从高楼大厦到桥梁隧道,从水利设施到工业厂房,无一不见其身影。
为了确保这些结构的安全性、可靠性和经济性,对其进行准确的力学分析至关重要。
Ansys 作为一款功能强大的有限元分析软件,为钢筋混凝土建模提供了高效且精确的解决方案。
钢筋混凝土是一种由钢筋和混凝土两种材料共同作用的复合材料。
混凝土具有较高的抗压强度,但抗拉强度较低;而钢筋则具有良好的抗拉性能。
在实际结构中,两者协同工作,共同承受外力。
因此,在Ansys 中进行钢筋混凝土建模时,需要准确地模拟这两种材料的特性以及它们之间的相互作用。
首先,我们来谈谈混凝土的建模。
在 Ansys 中,混凝土通常可以采用实体单元进行模拟。
对于混凝土的本构关系,我们可以选择合适的模型,如经典的混凝土损伤塑性模型(Concrete Damaged Plasticity Model)。
这个模型能够较好地考虑混凝土在受压和受拉时的非线性行为,包括混凝土的开裂、压碎等现象。
在定义混凝土的材料参数时,需要输入诸如弹性模量、泊松比、抗压强度、抗拉强度等参数。
这些参数的准确取值对于模型的准确性至关重要。
一般来说,可以通过实验测试或者参考相关的规范和标准来获取这些参数。
接下来是钢筋的建模。
钢筋在 Ansys 中有多种建模方法,常见的有两种:一种是使用杆单元(Link Element)来模拟钢筋,另一种是将钢筋嵌入到混凝土实体单元中(Embedded Element)。
使用杆单元模拟钢筋时,需要定义钢筋的截面积、弹性模量、屈服强度等参数。
这种方法计算效率较高,但对于钢筋与混凝土之间的粘结滑移行为模拟不够精确。
将钢筋嵌入到混凝土实体单元中的方法能够更准确地考虑钢筋与混凝土之间的相互作用,但计算量相对较大。
在这种方法中,需要确保钢筋单元与混凝土单元之间的节点协调。
在钢筋混凝土建模中,还需要考虑钢筋与混凝土之间的粘结滑移。
基于ANSYS的钢筋混凝土梁非线性分析若干问题研究
20 1 1 4月 年
水利 与建 筑工程 学报
0 fⅥ trReo re n c e tr ] e s u c sa d Ar ̄tcu a
V0 . o 2 19 N . p ., r 2011
基 于 A S S的 钢 筋 混 凝 土 梁 非 线 性 NY 分 析 若 干 问题 研 究
Re e r h O lSe e a o l ms a o tNo ln a s a c i v r lPr b e b u n i e r Ana y i f l ss o RC a s Ba e n ANS Be m s d o YS
GA0 ,Z Li ENG a . u 2 W ANG n 2 Xio y n Yo g
单元选择 、 材料本构关系 、 加速 收敛等若干 问题 , A S S对一钢 筋混凝土简 支梁进 行非线 性分析 , 用 NY 研 究梁荷载与跨中挠度关系 、 刚度变化及裂缝 发展形态 , 为类似的模 拟分析提供可靠 的理论依据 。
关 键 词 :A S S 钢 筋 混 凝 土 梁 ; 线 性 ;有 限元 NY; 非 中 图分 类 号 : U 7 . T 35 1 文 献 标识 码 : A 文 章 编 号 :17— 14 (0 1o—0 9—0 62 14 2 1)2 02 3
0 前
言
1 非线性分析时应注 意的若干问题
1 1 关 于模 型 .
钢筋混 凝土 是土木 行业 中应用 最广泛 的工程 材 料, 其性 质 复 杂 , 料 非 线 性 和几 何 非 线 性 同 时存 材
钢筋混 凝土 有限元 模 型根据 钢筋 的处理方 式主
在 。用传 统方法 分 析 往往 难 度 较 大 , 有 限元 方 法 而 在钢筋 混凝 土非线性 分析 中显示 出越来 越 大的实用
基于ANSYS的混凝土板的仿真计算
化 , 强化 切线模 量取为 0 0 E, 其 .1 其应力一应变关系如图 2所示。
1 模型 的 建立 1 1 模 型 的主要 参 数 .
混 凝 土 弹 性 模 量 E =3×1 m , 松 比 =0 1 7 混 0 N/ m2 泊 .6 ,
凝土单轴抗拉强度为 3 125MP , . 1 a单轴抗拉 强度 为 2 .6MP 。 4 0 a 模型的尺寸为 60 0ln×600ll, 0 l. 0 fl q3 qT 板的厚度为 10Fl, 、 l 5 lt 下表 f 上 T 层 的厚度 均为 3 考 虑作用在 双 向板上 的荷载 为 1 N/ 2 0mm, 0k m 。 横 向和纵 向均配筋 HR 3 5 2 5 , B 3  ̄1 @10 泊松 比 =0 2 , .5 钢筋 的屈 服应力设计值 为 3 1 6MP 。 3 . a
14 混 凝 土 本 构 关 系及 强 度 模 型 .
1 混 凝 土 的本 构 关 系 。混 凝 土 材 料 性 质 复 杂 , 不 仅 成 分 多 ) 它
图 3 有限元加载模型
S S计算 结 果分 析 样, 而且混凝 土的应力一 应变关系是高 度非线性 的 , 受其组成 、 3 AN Y 且 通过采用通用有限元分析软件 AN Y S S对现浇混凝 土空心板 成型工艺和使用环境 的严 重影 响 。特别 是在 复杂 的应 力状 态和
其有 限元加 载模 型如 图 3 所示 , 载荷大小为 1 N/ 2 0k m 。 i5 d 单元有八 个 节点 。每 个节 点 都具 有 沿 3, 方 向的 自由 6 5 Y, " 载,
度, 在三维等参元 S l4 oi 5的基础上 ,oi 5单元 增加 了针 对混凝 d Sl6 d 土的材料参数和组合 式钢筋 模型 。Sl 6 许 每个单 元有 4种 oi 5允 d
基于ANSYS的混凝土泵车上车结构非线性计算
截面积 ,考虑到油缸长度随臂架 的转动而变化 ,为
处理。
保证油缸质量不变 ,模型中将油缸密度进行了等效 其 中
( 端 部 软管 简化 。以质量 点形 式 模拟 ,采 用 4)
一结构重力冲击系数 ,取 1 . 2;
庀 混凝 土重 力 冲击 系数 ,取 1 。 . 3
5节 臂 风 载
臂 节所 受风载 为面载荷 ,模 型 中以均布载 荷
,
^ 一
x
1 -C n en 0 e , J ol cl t n d s d Y
基于 A Y NS S的混凝土泵 车上 车结构非线性计算
孙武 和 ,周振平 1 ,冯
( . 林大 学 1吉
敏z ,赵 二飞 -
长春 10 2 3 0 5: 徐州 2 10 2 0 4)
机 械科 学 与工程 学 院力 学实验 中心 ,吉 林
2 徐 工集 团工程 机械 股份 有 限公 司建 设机械 分 公司 ,江苏 .
[ 要 ] 以混 凝土 泵 车上 车结构 为研 究对 象 ,应用 A S S软件 对其 有 限元模 型进 行线 性 与非线 性计 摘 NY 算 ,通 过 两者 计 算 结果 对 比 ,给 出结 构 最大 位 移 的定 性 分析 ,其 过 程及 结 论对 臂 架设 计 及 改 型具 有 参考 价值 和指 导 意义 。 [ 关键 词 ]混 凝土 泵车 ; 车结构 ; 限元 ; 上 有 几何非 线性 [ 中图分 类号 ] U 4 T 66 [ 文献标 识码 ] B [ 文章编 号 ]10 — 5X ( 0 1 0— 04 0 0 154 2 1 ) 7 0 8— 4
IE 6管 NS 利用有限元多种单元类型的特点 ,对 Po r E模 PP 1 单 元 是 A YS中 具 有 拉 、压 、扭 、弯 / 节点组成 ,每节点具有 6自由度 , 型进行 中面提取 和几 何处 理 ,划分 网格 ,添加 有 限 的管单元 ,由 2
ANSYS中混凝土的计算问题
发信人: rubors (宝马), 信区: FEA标题: 混凝土单元的应用(solid65)[转载]发信站: 同舟共济站(2002年09月08日17:16:34 星期天), 站内信件ANSYS中混凝土的计算问题【精华】最近做了点计算分析,结合各论坛关于这方面的讨论,就一些问题探讨如下,不当之处,敬请指正。
一、关于模型钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种,即分离式、分布式和组合式模型。
考虑钢筋和混凝土之间的粘结和滑移,则采用引入粘结单元的分离式模型;假定混凝土和钢筋粘结很好,不考虑二者之间的滑移,则三种模型都可以;分离式和分布式模型适用于二维和三维结构分析,后者对杆系结构分析比较适用。
裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型,后者目前尚处研究之中,主要应用的是前两种。
离散裂缝模型和分布裂缝模型各有特点,可根据不同的分析目的选择使用。
随着计算速度和网格自动划分的快速实现,离散裂缝模型又有被推广使用的趋势。
就ANSYS而言,她可以考虑分离式模型(solid65+link8,认为混凝土和钢筋粘结很好,如要考虑粘结和滑移,则可引入弹簧单元进行模拟,比较困难!),也可采用分布式模型(带筋的solid65)。
而其裂缝的处理方式则为分布裂缝模型。
二、关于本构关系混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类,其中研究最多的是非线性弹性和弹塑性本构关系,其中不乏实用者。
混凝土破坏准则从单参数到五参数模型达数十个模型,或借用古典强度理论或基于试验结果等,各个破坏准则的表达方式和繁简程度各异,适用范围和计算精度差别也比较大,给使用带来了一定的困难。
就ANSYS而言,其问题比较复杂些。
1.ANSYS混凝土的破坏准则与屈服准则是如何定义的?采用tb,concr,matnum则定义了W-W破坏准则(failure criterion),而非屈服准则(yield criterion)。
ANSYS,ADINA在钢筋混凝土结构非线性分析中的应用与算例分析
表1跨中截面的应变和应力计算 结果对比
结论与展望
结论与展望
通过以上分析和算例对比,我们可以得出以下结论:ANSYS和ADINA在钢筋混 凝土结构非线性分析中均具有较高的准确性和可靠性,二者的计算结果非常接近。 这些软件提供了丰富的建模功能、强大的求解器和灵活的材料本构模型选择,使 得它们能够广泛应用于各种复杂的钢筋混凝土结构非线性分析问题中。
ANSYS和ADINA在钢筋混凝土结构非线性分析中的应用步骤和常见问题
钢筋混凝土结构非线性分析基本原理和方法
ANSYS和ADINA都是广泛用于工程模拟的有限元软件,它们在钢筋混凝土结构 非线性分析中具有相似的应用步骤。具体而言,可以分为以下几个步骤:
钢筋混凝土结构非线性分析基本原理和方法
1、建立模型:根据实际结构,建立相应的计算模型,包括几何形状、材料属 性、边界条件等。
谢谢观看
ANSYS,ADINA在钢筋混凝土结构 非线性分析中的应用与算例分
析
目录
01 引言
03 表1跨中截面的应变 和应力计算结果对比
02 算例分析 04 结论与展望 Nhomakorabea 内容摘要
ANSYS和ADINA在钢筋混凝土结构非线性分析中的应用与算例分析
引言
引言
钢筋混凝土结构非线性分析在工程实践中具有重要意义,它能够帮助工程师 们更准确地预测结构的性能和行为。在众多非线性分析软件中,ANSYS和ADINA是 最常用的两个工具。本次演示将介绍这两个软件在钢筋混凝土结构非线性分析中 的应用原理、方法和算例,并对其优势和不足进行分析,以期为工程师们提供有 益的参考。
整体式钢筋混凝土模型算例 (ANSYS)
在土木工程结构中,最为常用的一种结构形式就是钢筋混凝土结构,在各类房屋、水坝、桥梁、道路中都有广泛应用。
ANSYS软件提供了专门的钢筋混凝土单元和材料模型。
本算例将介绍ANSYS软件分析混凝土一些基本应用。
(1) 首先建立有限元模型,这里我们选用ANSYS软件自带的专门针对混凝土的单元类型Solid 65,进入ANSYS主菜单Preprocessor->Element Type->Add/Edit/Delete,选择添加Solid 65号混凝土单元。
(2) 点击Element types窗口中的Options,设定Stress relax after cracking 为Include,即考虑混凝土开裂后的应力软化行为,这样在很多时候都可以提高计算的收敛效率。
(3) 下面我们要通过实参数来设置Solid 65单元中的配筋情况。
进入ANSYS主菜单Preprocessor-> Real Constants->Add/Edit/Delete,添加实参数类型1与Solid 65单元相关,输入钢筋的材料属性为2号材料,但不输入钢筋面积,即这类实参数是素混凝土的配筋情况。
(4) 再添加第二个实参数,输入X方向配筋为0.05,即X方向的体积配筋率为5%。
(5) 下面输入混凝土的材料属性。
混凝土的材料属性比较复杂,其力学属性部分一般由以下3部分组成:基本属性,包括弹性模量和泊松比;本构关系,定义等效应力应变行为;破坏准则,定义开裂强度和压碎强度。
下面分别介绍如下。
(6) 首先进入ANSYS主菜单Preprocessor-> Material Props-> MaterialModels,在Define Material Model Behavior 窗口中选择Structural-> Linear -> Elastic-> Isotropic,输入弹性模量和泊松比分别为30e9和0.2 (7) 下面输入混凝土的等效应力应变关系,这里我们选择von Mises屈服面,该屈服面对于二维受力的混凝土而言精度还是可以接受的。
用Ansys或Abaqus分析钢管混凝土结构或构件
用Ansys或Abaqus分析钢管混凝土结构或构件用Ansys或Abaqus分析钢管混凝土结构或构件以上两个软件国外都有人用来分析钢管混凝土结构,但建模的方法不尽相同。
关键在于钢管和混凝土本构关系的选取以及两者之间的界面处理方法,各位有没有这方面的经验能向我们大家介绍一下。
==========程序中大概只有Drucker-Prager比较适合描述受约束混凝土的本构关系,因为这个模型可以考虑 hydrostatic stress (流体静应力)的影响。
在程序中,需要输入cohesion, angle of internal friction,(one more for ANSYS is the angle of dilatancy)。
值得注意的是,两个软件确定这几个参数的公式各不相同,很是令人头疼。
其实user manuals不可能给出明确的表达式,因为到目前为止,好像没有研究把钢管的强度,混凝土的强度,含钢率等等因素(i.e. the confinement)全部在Drucker-Prager 中考虑进去。
至于两种材料的界面,日本的 Hanbin Ge曾用link element来模拟,但在他的文章中,没有详细的描述。
轴压状况下,好像可以忽略滑移。
偏压可能情况有所不同。
==========韩教授书上的混凝土应力-应变关系,可以简单理解为单向受力的混凝土本构关系(考虑了钢管的约束),因此不能用于多向应力状态下混凝土的有限元分析。
材料非线性有限元分析,需要定义材料的屈服面,流动准则,强化准则,等等。
对受约束的混凝土,还要考虑体积膨胀,钢管对它的约束等因素。
显然,不是一个简单的应力-应变曲线所能概括的。
==========三向有限元分析,需要定义屈服面、流动准则和强化准则等等,而考虑钢管约束的混凝土本构关系,只是应力-应变关系。
对钢管混凝土的有限元分析,主要困难是如何定义屈服面,和模拟两个材料之间的滑移,我曾经用过接触分析(contact analysis)来求轴压构件的承载力,发现最大承载力能够比较精确地求得,但是精确的荷载-位移曲线很难获得,因为商用软件(Ansys\Marc)里面的D-P模型是塑性模型。
简支箱梁ansys分析例题
简支箱梁约束扭转算例一、工程背景已知某预应力混凝土简支箱梁,计算跨径为40m,沿梁长等截面。
截面尺寸如图1所示。
采用C40混凝土,剪切模量G=1.445×104MPa,弹性模量E=3.40×104MPa。
荷载为跨中作用一偏心荷载P=451.0kN,偏心距e=2.35m(计算约束扭转时,可以简化为集中力矩M k=1060.0kN)。
具体分两个工况进行:(1)跨中截面腹板位置作用一对对称集中竖向荷载,荷载大小为P/2=225.5kN;(2)跨中截面腹板位置作用一对反对称集中竖向荷载,荷载大小为P/2=225.5kN。
分别计算跨中截面、1/4跨位置截面上的正应力与剪应力分布,并绘制相应的正应力和剪应力分布曲线。
图1 箱梁截面尺寸(尺寸单位:cm)二、Ansys计算分析采用壳单元计算,荷载采用一对称荷载和反对称荷载加载。
1、Ansys命令流命令流见附录。
2、计算结果图a.工况一(a) 1/2跨正应力云图和应力曲线图(b) 1/2跨剪应力云图和应力曲线图b.工况二(a) 1/4跨正应力云图和应力曲线图(b) 1/4跨剪应力云图和应力曲线图注:由于路径选择时重复了一段,因此右翼缘由应力重叠现象附录命令流finish/clear/PREP7 ANTYPE,STATIC !定义箱梁厚度(单位:米)t1=0.22t2=0.30t3=0.34ET,1,SHELL63!定义单元R,1,t1R,2,t2R,3,t3!定义实常数跟所选单元有关本单元采用壳单元只需要厚度MP,EX,1,3.4E10!弹性模量1为材料参考号MP,PRXY,1,0.3!泊松比K,1,-4.75,0.955K,2,-2.35,0.955K,3,-2.35,-1.165K,4,2.35,-1.165K,5,2.35,0.955K,6,4.75,0.955K,7,2.35,0.955,40!创建关键点L,1,2$L,2,3$L,3,4L,4,5$L,5,6$L,5,2$L,5,7!连接各点ADRAG,1,6,5,,,,7ADRAG,2,3,4,,,,7!创建面根据直线编号2,3 路径线7ASEL,S,LOC,Y,0.955$AATT,1,1,1!选择面S选择新面作为子集loc按坐标值选择ASEL,S,LOC,X,-2.35$AA TT,1,2,1!Y中心坐标y值为0.955的面ASEL,S,LOC,X,2.35$AATT,1,2,1!设置面的单元属性相应定义点线体的命令不同ASEL,S,LOC,Y,-1.165$AATT,1,3,1!材料号,实常数,单元类型ALLSEL,ALL!选择所有图素NUMMRG,ALL !合并实体单元NUMCMP,ALL!压缩定义编号/PNUM,LINE,1/PNUM,AREA,1!显示面和线的编号LSEL,S,LENGTH,,40!选择长度为40的线段LESIZE,ALL,,,100!将所有的线每根分成100份即上述长度为40的线段LSEL,S,LENGTH,,2.4,4.7!选择长度在2.4-4.7长度单位内的线段LESIZE,ALL,,,5ALLSEL,ALLAMESH,ALL!划分网格FINISH!施加边界条件/SOLU!工况1对称集中竖向荷载DL,3,,UX$DL,3,,UY$DL,3,,UZDL,16,,UY!对线进行约束!F,662,FY,-225500!F,57,FY,-225500!加上集中荷载节点力F,662,FY,-225500!工况2反对称集中竖向荷载F,57,FY,225500SOLVEFINISH/POST1PATH,STRESS1,7!由于这个箱梁的特点,使得选路径时重复了一段,在应力图中有反映PPA TH,1,,-4.75,0.955,20$PPATH,2,,-2.35,0.9 55,20PPA TH,3,,-2.35,-1.165,20$PPATH,4,,2.35,-1.1 65,20PPA TH,5,,2.35,0.955,20$PPATH,6,,4.75,0.955 ,20PPA TH,7,,-2.35,0.955,20!七个路径点PDEF,STRESS1,S,ZPLPAGM,STRESS1,5.0!跨中截面正应力云图PLPATH,STRESS1!跨中截面正应力曲线PDEF,STRESS1,S,XYPLPAGM,STRESS1,5.0!跨中截面剪应力云图PLPATH,STRESS1!跨中截面剪应力曲线PATH,STRESS2,7PPA TH,1,,-4.75,0.955,10$PPATH,2,,-2.35,0.9 55,10PPA TH,3,,-2.35,-1.165,10$PPATH,4,,2.35,-1.1 65,10PPA TH,5,,2.35,0.955,10$PPATH,6,,4.75,0.955 ,10PPA TH,7,,-2.35,0.955,10PDEF,STRESS2,S,ZPLPAGM,STRESS2,5.0!1/4跨截面正应力云图PLPATH,STRESS2!1/4跨截面正应力曲线PDEF,STRESS2,S,XYPLPAGM,STRESS2,5.0!1/4跨截面剪应力云图PLPATH,STRESS2!1/4跨截面剪应力曲线。
基于ANSYS平台的混凝土徐变计算
摘
要: 通过 比较 , 选取 国 内规 范使 用徐 变 系数 的 3种拟 合 方法 中最好 的一 种 , 据 叠加 根
原 理 , 用 A S S二 次开发接 口编制 了相 应 的 A S S用 户子 程 序 , A S S平 台上 实现 了 利 NY NY 在 NY
混凝土的三维徐 变计算。在此基础上, 通过实例验证 了本文方法在 实际工程计算中的有效性。
维普资讯
第 3 卷第 3 2 期
20 0 6年 9月
湖
南
交
通
科
技
V 13 . 0 . 2 No 3 S p. 0 e 2 06
HU NAN C OMMUNI A I N CI NC C T O S E E AND T C E HNOL GY O
厶
【 、,川 l+r ’ l , r
() 、 一 1
警
单 位 时间 内的徐 变应 变增量 为 :
A = 占 t £r )一占 tr ) = e ( +△ ,0 ( ,0 【 £ tT )一 ( ,O ( +A ,O £T )】+ 维普资讯 Nhomakorabea3 期
何继访 , 基于 A S S 台的混凝土徐变计算 等: NY 平
表 3 l (l ^ 20r 拟合结果 比较 f t (。= 0 a / m}
段明德方法 钟 一颜方法 赵新 生方法
l9 0
S S 是通过求解蠕变的本构微分方程 , Y ̄ 得出单位时 间间隔内蠕变应变增量来实现蠕 变计算 。因此 , 利 用 A S S 蠕变计算功能计算混 凝土徐变 , NY ̄ 需要求 解单位时间间隔内混凝土的徐变应变增量 。 当混 凝 土 的应 力小 于 其强 度 的 5 % 时 , 0 混凝 土
ANSYS中混凝土的计算问题 本构模型
ANSYS中混凝土的计算问题最近做了点计算分析,结合各论坛关于这方面的讨论,就一些问题探讨如下,不当之处敬请指正。
一、关于模型钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种,即分离式、分布式和组合式模型。
考虑钢筋和混凝土之间的粘结和滑移,则采用引入粘结单元的分离式模型;假定混凝土和钢筋粘结很好,不考虑二者之间的滑移,则三种模型都可以;分离式和分布式模型适用于二维和三维结构分析,后者对杆系结构分析比较适用。
裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型,后者目前尚处研究之中,主要应用的是前两种。
离散裂缝模型和分布裂缝模型各有特点,可根据不同的分析目的选择使用。
随着计算速度和网格自动划分的快速实现,离散裂缝模型又有被推广使用的趋势。
就ANSYS而言,她可以考虑分离式模型(solid65+link8,认为混凝土和钢筋粘结很好,如要考虑粘结和滑移,则可引入弹簧单元进行模拟,比较困难!),也可采用分布式模型(带筋的solid65)。
而其裂缝的处理方式则为分布裂缝模型。
二、关于本构关系混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类,其中研究最多的是非线性弹性和弹塑性本构关系,其中不乏实用者。
混凝土破坏准则从单参数到五参数模型达数十个模型,或借用古典强度理论或基于试验结果等,各个破坏准则的表达方式和繁简程度各异,适用范围和计算精度差别也比较大,给使用带来了一定的困难。
就ANSYS而言,其问题比较复杂些。
1 ANSYS混凝土的破坏准则与屈服准则是如何定义的?采用tb,concr,matnum则定义了W-W破坏准则(failure criterion),而非屈服准则(yi eld criterion)。
W-W破坏准则是用于检查混凝土开裂和压碎用的,而混凝土的塑性可以另外考虑(当然是在开裂和压碎之前)。
理论上破坏准则(failure criterion)和屈服准则(y ield criterion)是不同的,例如在高静水压力下会发生相当的塑性变形,表现为屈服,但没有破坏。
ansys中混凝土结构的计算
(1)转贴一个例题,来自中国有限元联盟论坛 -运行没有问题/title, fixed - fixed concrete beam example/prep7et,1,65mp,ex,1,3e7 ! steel rebar (units are pounds, inches)MP,PRXY, 1,0.3mp,ex,2,1e6, ! concreteMP,PRXY, 2,0.3mp,dens,2,.00025tb,concr,2tbdata,1,.3,.5,200,4000 ! shear coeffs, tensile and compress strength r,1,1,.03,0,0 ! mat 1 (steel), 3 percent reinforcement in x dirr,2,1,.01,0,0 ! mat 1 (steel), 1 percentr,3,1,.04,0,0 ! mat 1 (steel), 4 percentblock,,100,,5,,5block,,100,5,10,,5block,,100,10,15,,5!vovlap,allNUMMRG,KP, , , ,LOWnumcmp,voluesize,5mat,2 ! concrete materialreal,1 ! rebarvmesh,1real,2vmesh,2real,3vmesh,3nsel,s,loc,xd,all,allnsel,s,loc,x,100d,all,allnsel,allfini/solunsel,s,loc,y,15sf,all,pres,100nsel,allOUTRES,ALL,ALL,nsub,10solvefini/post1/DEVICE,VECTOR,1set,lastplcrackfini(2)! ANSYS SOLID65环向布置钢筋的例子――运行没有问题! 一个管道,环向配筋率为1%,纵向配筋率为0.5%,径向配筋率为0.1%! 作者: 陆新征清华大学土木系!FINISH/CLEAR/PREP7!*! 单元属性ET,1,SOLID65!*KEYOPT,1,1,0 !表示考虑大变形KEYOPT,1,5,0 !表示只打印质心的线性解KEYOPT,1,6,0 !表示只打印质心的线性解KEYOPT,1,7,1 !表示开裂后考虑应力松弛,有助于计算收敛!*!实参数1:不同方向配筋R,1,2,.001, , ,2, .01,RMORE, 90, ,2,.005 ,90 ,90 ,!材料属性!混凝土基本材料属性MPTEMP,,,,,,,,MPTEMP,1,0MPDATA,EX,1,,30E2MPDATA,PRXY,1,,.2!屈服准则TB,MISO,1,1,5, !多线性等向强化,材料号为1,数据的温度数为1,对某一给定温度数据的点数为5。
ansys作业
高速弹丸侵彻混凝土靶板一,问题描述一个直径D=15mm,长度L=43mm的圆柱弹丸,材料为钢,以1500m/s的速度侵彻两层混凝土靶板。
混凝土靶板尺寸为400mm╳400mm╳60mm,四周边界固定,试模拟混凝土靶板的破坏现象。
二,建模分析混凝土靶板采用*MAT_SOLID-CONCRETE,采用四分之一模型,三维实体单元划分网格。
对称面上施加对称边界约束,靶板边界施加固定约束。
弹丸和混凝土靶板之间采用侵蚀接触。
采用cm-g-µm单位制。
计算时间为50µms,每1µs输出一个结果数据文件。
三,求解步骤第一步:设置工作目录和工作文件1,D盘新建文件夹penetration3d-concrete.2,启动Ansys 10.0,弹出10.0Launcher窗口,如下图在simulation enviroment 下拉框中选择ANSYS, 在License下拉框中选择ANSYS LS—DYNA;3,单击File Management 选项卡,弹出工作目录和工作文件设置窗口,单击Browse 按钮,在Working Directory框中选择D:\penetration3d-concrete,在Job Name 框中输入penetration3d-concrete,单击Run按钮,进入ANSYS 操作界面。
如图第二步:选择单元类型1,选择Main Menu>Preprocessor>Element Type>Add/Edit/Delete命令,弹出Element Type对话框;如图2,单击Add按钮,弹出Library of Element Type对话框;3,在Library of Element Types选择栏中选择LS—DYNA Explicit 3D Solid 164,单击OK按钮,关闭对话框;4,在Element Types对话框中单击Options按钮,弹出Solid 164 Element Types Options对话框;5,在Solid Element Formulation 栏中激活Const.Stress(def)选项,在Material Continuum栏中激活Lagrangian选项,单击OK按钮,关闭对话框;6,在Element Types对话框中单击Close按钮,关闭对话框。
ansys分析混凝土的若干问题
ansys分析混凝土的若干问题1. 讨论两种Ansys求极限荷载的方法(1)力加载可以通过对应的方法(比如说特征值屈曲)估计结构的极限荷载的大致范围,然后给结构施加一个稍大的荷载,打开自动荷载步二分法进行非线性静力分析,最后计算会因不收敛终止,则倒数第二个子步对应的就是结构的极限荷载;另外,也可以选择弧长法,采用足够的子步(弧长法可以一直分析到极限承载力之后的过程)同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。
(2)位移加载给结构施加一个比较大的位移,打开自动荷载步二分法进行非线性分析,保证足够的子步数,这样也可以分析到极限荷载以后,通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。
希望众高手讨论一下(1)弧长法求极限荷载的收敛性问题,如何画到荷载位移曲线的下降段?(2)位移法求极限荷载的具体步骤?2. 需要注意的问题1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据,因此在很多情况下会因为应力集中而使混凝土提前破坏,从而和试验结果不相吻合,因此,在实际应用过程中应该对单元分划进行有效控制,根据作者经验,当最小单元尺寸大于5cm 时,就可以有效避免应力集中带来的问题;2. 支座是另一个需要注意的问题。
在有限元分析中,很多时候约束是直接加在混凝土节点上,这样很可能在支座位置产生很大的应力集中,从而使支座附近的混凝土突然破坏,造成求解失败。
因此,在实际应用过程中,应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块,避免支座的应力集中;3. 六面体的SOLID 65 单元一般比四面体的单元计算要稳定且收敛性好,因此,只要条件允许,应该尽量使用六面体单元;4. 正确选择收敛标准,一般位移控制加载最好用位移的无穷范数控制收敛,而用力控制加载时可以用残余力的二范数控制收敛。
在裂缝刚刚出现和接近破坏的阶段,可以适当放松收敛标准,保证计算的连续性;3. 关于下降段的问题1)在实际混凝土中都有下降段,但是在计算的时候要特别小心下降段的问题。
ANSYS混凝土HJC模型
HJC 混凝土模型介绍最近应用在计算机模型上的混凝土模型已经投入了特别多的努力。
大部分模型大部分应用在相对小应变、低应变率和低压力情况下即支持民用土木工程。
而在混凝土冲击计算中遭受大应变高应变率和高压力的混凝土模型方面投入的很少。
这篇论文为遭受大应变高应变率和高压力的混凝土提出了一个连续计算模型。
适合在计算机中用拉格朗日和欧拉算法。
这个模型类似于Osborn 【1】发展的模型,但是提出的模型扩展到包括材料损伤、应变率效果、和把压力和空隙率作为参数的永久压缩函数。
论文的末尾对模型进行了描述,测定了非限制压缩强度为0.048GPa 的混凝土的参数。
侵彻计算与试验结果进行了比较。
模型的描述 模型的大体描述见图一。
模型的强度部分见图一的最顶。
归一化强度被定义为/*c f σσ=这里σ*为实际等效应力,f ’C 是准静态异轴向压缩强度,具体表达式为: ()[]⎥⎥⎦⎤⎢⎢⎣⎡⎪⎪⎭⎫ ⎝⎛++-=*.**ln 11εσC BP D A N (1) 此处D 是损伤(10≤≤D ),P *=P/f ‘C 是归一化强度(此处P 是实际的压力)。
0,,*,εεε=是无参数的应变率(,ε是实际应变率,0,ε=1.0S -1是参考应变率)。
归一化最大静水压力Cf T T ,*=,此处T 是材料能够抵抗的最大的静水压力(maximum tensile hydrostatic pressure ) 材料常数为A 、B 、N 、C 、SMAX 。
A 为归一化粘结强度,B 为归一化压力的强化系数,N 是压力强化指数,C 是应变率系数。
SMAX 是即能发展的归一化最大强度。
损伤的的累积方式类似于johnson-cook 模型,johnson-cook 模型从等效塑性应变来积累损伤。
此处讨论模型的累积损伤来自等效塑性应变和塑性体积应变。
用下式来表达: f Pf P P P D μεμε+∆+∆∑= (2) 此处P ε∆和P μ∆是等效塑性应变和塑性体积应变,相应的在一次积分循环中:f P f P με+=f(P)是在常压力下塑性损伤应变。
ANSYS混凝土问题分析
ANSYS混凝土问题分析1.关于模型钢筋混凝土有限元模型根据钢筋的处理方式分为三种:分离式、整体式和组合式模型◆分离式模型:把混凝土和钢筋作为不同的单元来处理,即混凝土和钢筋各自被划分为足够小的单元,两者的刚度矩阵是是分开来求解的,考虑到钢筋是一种细长的材料,通常可以忽略起横向抗剪强度,因此可以将钢筋作为线单元处理。
钢筋和混凝土之间可以插入粘结单元来模拟钢筋与混凝土之间的粘结和滑移。
一般钢筋混凝土是存在裂缝的,而开裂必然导致钢筋和混凝土变形的不协调,也就是说要发生粘结的失效与滑移,所以此种模型的应用最为广泛.◆整体式模型:将钢筋分布与整个单元中,假定混凝土和钢筋粘结很好,并把单元视为连续均匀材料,与分离式模型不同的是,它求出的是综合了混凝土与钢筋单元的整体刚度矩阵;与组合式不同之点在于它不是先分别求出混凝土与钢筋对单元刚度的贡献然后再组合,而是一次求得综合的刚度矩阵.◆组合式模型组合式模型分为两种:一种是分层组合式,在横截面上分成许多混凝土层和若干钢筋层,并对截面的应变作出某些假设,这种组合方式在钢筋混凝土板、壳结构中应用较广;另一种组合方法是采用带钢筋膜的等参单元。
当不考虑混凝土和钢筋二者之间的滑移,三种模型都可以。
分离式和整体式模型使用于二维和三维结构分析。
就ANSYS而言,可以考虑分离式模型:混凝土(SOLID65)+钢筋(LINK单元或PIPE 单元),认为混凝土和钢筋粘结很好。
如要考虑粘结和滑移,则可引入弹簧单元进行模拟,如果比较困难也可以采用整体式模型(带筋的SOLID65)。
2.本构关系及破坏准则◆本构关系混凝土本构关系的模型对钢筋混凝土结构的非线性分析有重大影响。
混凝土的本构就是表示在各种外荷载作用下的混凝土应力应变的响应关系。
在建立混凝土本构关系时一般都是基于现有的连续介质力学的本构理论,在结合混凝土的力学特性,确定甚至调整本构关系中各种所需的材料参数.通常,混凝土的本构关系可以分为线性弹性、非线性弹性、弹塑性及其他力学理论等四类。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
ANSYS中混凝土的计算问题【精华】最近做了点计算分析,结合各论坛关于这方面的讨论,就一些问题探讨如下,不当之处敬请指正。
一、关于模型钢筋混凝土有限元模型根据钢筋的处理方式主要分为三种,即分离式、分布式和组合式模型。
考虑钢筋和混凝土之间的粘结和滑移,则采用引入粘结单元的分离式模型;假定混凝土和钢筋粘结很好,不考虑二者之间的滑移,则三种模型都可以;分离式和分布式模型适用于二维和三维结构分析,后者对杆系结构分析比较适用。
裂缝的处理方式有离散裂缝模型、分布裂缝模型和断裂力学模型,后者目前尚处研究之中,主要应用的是前两种。
离散裂缝模型和分布裂缝模型各有特点,可根据不同的分析目的选择使用。
随着计算速度和网格自动划分的快速实现,离散裂缝模型又有被推广使用的趋势。
就ANSYS而言,她可以考虑分离式模型(solid65+link8,认为混凝土和钢筋粘结很好,如要考虑粘结和滑移,则可引入弹簧单元进行模拟,比较困难!),也可采用分布式模型(带筋的solid65)。
而其裂缝的处理方式则为分布裂缝模型。
二、关于本构关系混凝土的本构关系可以分为线弹性、非线性弹性、弹塑性及其它力学理论等四类,其中研究最多的是非线性弹性和弹塑性本构关系,其中不乏实用者。
混凝土破坏准则从单参数到五参数模型达数十个模型,或借用古典强度理论或基于试验结果等,各个破坏准则的表达方式和繁简程度各异,适用范围和计算精度差别也比较大,给使用带来了一定的困难。
就ANSYS而言,其问题比较复杂些。
1 ANSYS混凝土的破坏准则与屈服准则是如何定义的?采用tb,concr,matnum则定义了W-W破坏准则(failure criterion),而非屈服准则(yield criterion)。
W-W破坏准则是用于检查混凝土开裂和压碎用的,而混凝土的塑性可以另外考虑(当然是在开裂和压碎之前)。
理论上破坏准则(failure criterion)和屈服准则(yield criterion)是不同的,例如在高静水压力下会发生相当的塑性变形,表现为屈服,但没有破坏。
而工程上又常将二者等同,其原因是工程结构不容许有很大的塑性变形,且混凝土等材料的屈服点不够明确,但破坏点非常明确。
定义tb,concr matnum后仅仅是定义了混凝土的破坏准则和缺省的本构关系,即W—W破坏准则、混凝土开裂和压碎前均为线性的应力应变关系,而开裂和压碎后采用其给出的本构关系。
但屈服准则尚可另外定义(随材料的应力应变关系,如tb,MKIN,则定义的屈服准则是V on Mises,流动法则、硬化法则也就确定了)。
2 定义tb,concr后可否定义其它的应力应变关系当然是可以的,并且只有在定义tb,concr后,有些问题才好解决。
例如可以定义tb,miso,输入混凝土的应力应变关系曲线(多折线实现),这样也就将屈服准则、流动法则、硬化法则等确定了。
这里可能存在一点疑问,即ANSYS中的应力应变关系是拉压相等的,而混凝土材料显然不是这样的。
是的,因为混凝土受拉段非常短,认为拉压相同影响很小,且由于定义的tb,concr中确定了开裂强度,所以尽管定义的是一条大曲线,但应用于受拉部分的很小。
三、具体的系数及公式1 定义tb,concr时候的两个系数如何确定?一般的参考书中,其值建议先取为0.3~0.5(江见鲸),原话是“在没有更仔细的数据时,不妨先取0.3~0.5进行计算”,足见此0.3~0.5值的可用程度。
根据我的经验和理由,建议此值取大些,即开裂的剪力传递系数取0.5,(定要>0.2)闭合的剪力传递系数取1.0。
支持此说法的还有现行铁路桥规的抗剪计算理论,以及原公路桥规的容许应力法的抗计剪计算。
2 定义混凝土的应力应变曲线单向应力应变曲线很多,常用的可参考国标混凝土结构规范,其中给出的应力应变曲线是二次曲线+直线的下降段,其参数的设置按规范确定即可。
当然如有实测的应力应变曲线更好了。
四、关于收敛的问题ANSYS混凝土计算收敛(数值)是比较困难的,主要影响因素是网格密度、子步数、收敛准则等,这里讨论如下。
1 网格密度:网格密度适当能够收敛。
不是网格越密越好,当然太稀也不行,这仅仅是就收敛而言的,不考虑计算费用问题。
但是究竟多少合适,没有找到规律,只能靠自己针对情况慢慢试算。
2 子步数:NSUBST的设置很重要,设置太大或太小都不能达到正常收敛。
这点可以从收敛过程图看出,如果F范数曲线在[F]曲线上面走形的很长,可考虑增大nsubst。
或者根据经验慢慢调正试算。
3 收敛精度:实际上收敛精度的调正并不能彻底解决收敛的问题,但可以放宽收敛条件以加速吧。
一般不超过5%(缺省是0.5%),且使用力收敛条件即可。
4 混凝土压碎的设置:不考虑压碎时,计算相对容易收敛;而考虑压碎则比较难收敛,即便是没有达到压碎应力时。
如果是正常使用情况下的计算,建议关掉压碎选项;如果是极限计算,建议使用concr+MISO且关闭压碎检查;如果必设压碎检查,则要通过大量的试算(设置不同的网格密度、NSUBST)以达到目的,但也很困难。
5 其他选项:如线性搜索、预测等项也可以打开,以加速收敛,但不能根本解决问题。
6 计算结果:仅设置concr,不管是否设置压碎,其一般P-F曲线接近二折线;采用concr+miso则P-F曲线与二折线有差别,其曲线形状明显是曲线的。
例题1!题目:钢筋混凝土简支梁模拟计算!方法:分离式;solid65和link8!材料:混凝土采用concr和钢筋为弹性材料,但不考虑压碎,为方便,假定钢筋置于梁底两侧./config,nres,2000/prep7!定义单元及其材料特性等rd0=20.0 !钢筋直径et,1,solid65et,2,link8mp,ex,1,33e3mp,prxy,1,0.20r,1hntra=28hntrl=2.6tb,concr,1tbdata,,0.7,1.0,hntrl,-1mp,ex,2,2.1e5mp,prxy,2,0.30r,2,acos(-1)*0.25*rd0*rd0!定义梁体即单元划分blc4, , ,100,200,3000/view,1,1,1,1/ang,1gplot!定义网分时边长控制lsel,s,loc,z,1,2999lsel,r,loc,y,0latt,2,2,2lesize,all,,,20 !钢筋网格数目lmesh,alllsel,s,loc,z,0lesize,all,,,4 !截面上的网格数目4x4vsel,allvatt,1,1,1 mshape,0,3d mshkey,1 vmesh,all allsel,allfinish/solu!施加约束lsel,s,loc,z,0 lsel,r,loc,y,0dl,all,,uydl,all,,uzlsel,alllsel,s,loc,z,3000 lsel,r,loc,y,0dl,all,,uylsel,allksel,s,loc,x,0 ksel,r,loc,y,0 dk,all,ux allsel,all!施加荷载qmz=0.3asel,s,loc,y,200 sfa,all,1,pres,qmz allsel,all nsubst,40 outres,all,all time,qmz*10 neqit,40solvefinish/post1pldisp,1 etable,zxyl,ls,1 plls,zxyl,zxyl,1 finish/post26nsol,2,33,u,y prod,3,1,,,,,,1/100 prod,4,2,,,,,,-1 xvar,4plvar,3例题2!题目:钢筋混凝土简支梁模拟计算!方法:分离式;solid65和link8!材料:混凝土采用concr+Miso和钢筋为弹性材料,但不考虑压碎,增加网格密度,为方便,假定钢筋置于梁底两侧./config,nres,2000/prep7!定义单元及其材料特性等rd0=20.0 !钢筋直径et,1,solid65et,2,link8mp,ex,1,33e3mp,prxy,1,0.20r,1hntra=28hntrl=2.6tb,concr,1tbdata,,0.7,1.0,hntrl,-1tb,miso,1,,14tbpt,,0.0002,hntra*0.19tbpt,,0.0004,hntra*0.36tbpt,,0.0006,hntra*0.51tbpt,,0.0008,hntra*0.64tbpt,,0.0010,hntra*0.75tbpt,,0.0012,hntra*0.84tbpt,,0.0014,hntra*0.91tbpt,,0.0016,hntra*0.96tbpt,,0.0018,hntra*0.99tbpt,,0.0020,hntra*1.00tbpt,,0.0025,hntra*0.95tbpt,,0.0030,hntra*0.90tbpt,,0.0035,hntra*0.85tbpt,,0.0040,hntra*0.80mp,ex,2,2.1e5mp,prxy,2,0.30r,2,acos(-1)*0.25*rd0*rd0!定义梁体即单元划分blc4, , ,100,200,3000/view,1,1,1,1/ang,1gplot!定义网分时边长控制lsel,s,loc,z,1,2999 lsel,r,loc,y,0latt,2,2,2lesize,all,,,20 lmesh,alllsel,s,loc,z,0 lesize,all,,,4vsel,allvatt,1,1,1 mshape,0,3d mshkey,1vmesh,allallsel,allfinish/solu!施加约束lsel,s,loc,z,0lsel,r,loc,y,0dl,all,,uydl,all,,uzlsel,alllsel,s,loc,z,3000 lsel,r,loc,y,0dl,all,,uylsel,allksel,s,loc,x,0ksel,r,loc,y,0dk,all,uxallsel,all!施加荷载qmz=0.3asel,s,loc,y,200sfa,all,1,pres,qmz allsel,alloutres,all,alltime,qmz*10 nsubst,40neqit,40solvefinish/post1pldisp,1 etable,zxyl,ls,1 plls,zxyl,zxyl,1 finish/post26nsol,2,33,u,y prod,3,1,,,,,,1/100 prod,4,2,,,,,,-1 xvar,4plvar,3。