混凝土非线性分析-应力应变曲线输出概要
钢筋混凝土非线性分析
五、基本概念 1、本构关系：材料力学性质的数学表达式 2、屈服极限：由弹性变形变为非弹性变形的转折点的应力 屈服条件：某一点出现塑性变形时应力状态应满足的条件 屈服函数：表示屈服条件的函数 屈服面： 屈服函数在应力空间中表示的曲面
3、强化：屈服极限提高的现象 软化：应力降低、应变增大的现象 拉伸强化：混凝土受拉构件中主裂缝之间混凝土仍承担 一部分拉应力的现象 4、反复加载：周期性静力荷载作用下交替产生拉、压应力 重复加载：周期性静力荷载作用下仅产生单向应力
b)应力随时间变化，且σ＜0.5fc： •粘弹性流变模型——弹性继承理论
•老化理论
【朱】Page25 式1.41
假定：混凝土各向同性 瞬时应力、应变具有线性关系 徐变变形与应力之间具有线性关系 应力变化，徐变变形值按相应应力增量引起的徐变变形迭加 变形的绝对值与应力符号无关 在相同应力条件下，不同龄期的徐变曲线可以成为“平行曲 线”，徐变变形可按应力值和混凝土的龄期加以推算。 如已知龄期τi，较早龄期τ1，加载时刻t， 则：
第一章：绪论
一、学习非线性分析的意义 （当前混凝土结构设计存在的问题） 1、混凝土材料工作状态的非线性 2、钢筋和混凝土共同工作条件——变形协调 3、结构内力计算和截面设计不协调 4、节点的理想化（刚接、铰接）与实际状态不符 5、长期荷载下徐变、应力松弛引起的结构内力重分布 6、动力荷载作用下的材料特性与静力下不同
a)应力不变，且σ＜0.5fc （线性徐变或有限徐变）： 幂表达式 指数表达式 双曲线表达式 对数表达式
其中各常数可以调整，用以考虑 时间和不同因素的影响
在此基础上，另加调整参数，对表达式进行修正 【朱】Page24 式1.37 考虑自由收缩、水泥水化程度 式1.38、1.39 考虑湿度、尺寸、龄期 式1.40 考虑湿度、尺寸、龄期、配合比、其它
混凝土受压应力-应变全曲线方程(描述)
混凝土受压应力－应变全曲线方程混凝土受压应力－应变全曲线方程混凝土的应力－应变关系是钢筋混凝土构件强度计算、超静定结构内力分析、结构延性计算和钢筋混凝土有限元分析的基础，几十年来，人们作了广泛的努力，研究混凝土受压应力－应变关系的非线性性质，探讨应力与应变之间合理的数学表达式，1942年，Whitney 通过混凝土圆柱体轴压试验，提出了混凝土受压完整的应力应变全曲线数学表达式，得出了混凝土脆性破坏主要是由于试验机刚度不足造成的重要结论，这一结论于1948年由Ramaley 和Mchenry 的试验研究再次证实，1962年，Barnard 在专门设计的具有较好刚性且能控制应变速度的试验机上，试验了一批棱柱体试件以及试件两靖被放大的圆柱体试件，试验再次证明，混凝土的突然破坏并非混凝土固有特性，而是试验条件的结果，即混凝土的脆性破坏可用刚性试验机予以防止，后来由很多学者（如M.Sagin ，P.T.Wang ，过镇海等）所进行的试验，都证明混凝土受压应力－应变曲线确实有下降段存在，那么混凝土受压应力与应变间的数学关系在下降段也必然存在，研究这一数学关系的工作一刻也没有停止。
钢筋混凝土结构是目前使用最为广泛的一种结构形式。
但是，对钢筋混凝土的力学性能还不能说已经有了全面的掌握。
近年来，随着有限元数值方法的发展和计算机技术的进步，人们已经可以利用钢筋混凝土有限元分析方法对混凝土结构作比较精确的分析了。
由于混凝土材料性质的复杂性，对混凝土结构进行有限元分析还存在不少困难，其中符合实际的混凝土应力应变全曲线的确定就是一个重要的方面。
1、混凝土单轴受压全曲线的几何特点经过对混凝土单轴受压变形的大量试验大家一致公认混凝土单轴受压变过程的应力应变全曲线的形状有一定的特征。
典型的曲线如图1所示，图中采用无量纲坐标。
sc c E E N f y x 0,,===σεε 式中，c f 为混凝土抗压强度；c ε为与c f 对应的峰值应变；0E 为混凝土的初始弹性模量；s E 为峰值应力处的割线模量。
基于ABAQUS的混凝土破坏实验非线性分析
基于ABAQUS的混凝土破坏实验非线性分析一、内容综述近年来，随着建筑工程技术的飞速发展，现代建筑结构对于材料性能和结构安全性的要求也日益提高。
作为一种极具代表性的建筑材料，在各种工程领域中得到了广泛应用。
混凝土在复杂荷载下的破坏过程通常表现出强烈的非线性特性，这给结构的分析和设计带来了极大的挑战。
传统的线性分析方法在处理这类非线性问题时存在诸多局限性，基于非线性理论的混凝土破坏实验分析方法成为了研究的热点。
ABAQUS，作为目前国际上广泛应用的有限元分析软件，具备高效、精确的非线性分析能力，为混凝土破坏实验提供了有力的工具。
1.1 研究背景和意义随着现代建筑事业的不断发展，高层建筑和大跨度结构越来越普遍，这对混凝土结构的承载能力和抗震性能要求越来越高。
混凝土破坏实验是研究混凝土性能的重要手段之一，其中非线性分析方法在混凝土破坏实验中具有重要的应用价值。
ABAQUS作为国际上广泛应用的有限元分析软件，在混凝土破坏实验非线性分析中具有很高的精度和可靠性。
本文通过基于ABAQUS的混凝土破坏实验非线性分析，旨在研究混凝土在不同条件下的破坏模式和力学行为，为混凝土结构和结构的抗震设计提供理论依据。
本文的研究背景是考虑到现代建筑物对混凝土性能的高要求以及在地震、暴风等自然灾害中对结构安全性的挑战。
研究混凝土破坏实验非线性分析有助于深入了解混凝土内部的应力分布和破坏机制，为实际工程提供更加准确的计算和分析方法。
研究ABAQUS在混凝土破坏实验非线性分析中的应用，可以为相关领域的研究提供借鉴和参考，推动有限元分析技术在混凝土结构研究中的进一步发展。
1.2 ABAQUS软件简介ABAQUS（阿巴克斯）是一款国际知名的非线性有限元分析软件，广泛应用于工程领域的结构分析和模拟。
它具有全面、准确、可靠的特点，能有效地处理各种复杂材料和结构问题。
在混凝土破坏实验非线性分析中，ABAQUS软件发挥着重要作用，为研究者提供了便捷的工具来模拟混凝土在受力状态下的破坏过程。
钢筋混凝土板的非线性分析
钢筋混凝土板的非线性分析钢筋混凝土板的非线性分析钢筋混凝土板是一种常用的结构构件，在建筑和桥梁中广泛应用。
由于其在使用过程中会受到各种荷载的作用，因此需要对其进行非线性分析，以确保其安全可靠。
非线性分析是指在分析过程中考虑材料和结构的非线性特性，包括材料的本构关系、几何非线性和接触非线性等因素。
在钢筋混凝土板的非线性分析中，需要考虑以下几个方面。
1. 材料的本构关系钢筋混凝土板的材料包括混凝土和钢筋两部分，它们的本构关系是非线性的。
混凝土的本构关系可以采用双曲正切模型或Drucker-Prager 模型等进行描述，而钢筋的本构关系则可以采用弹塑性模型或Ramberg-Osgood模型等进行描述。
在进行非线性分析时，需要考虑这些材料的本构关系对结构的影响。
2. 几何非线性钢筋混凝土板在受到荷载作用后会发生变形，这种变形会导致结构的几何非线性。
几何非线性包括平面内的弯曲变形和平面外的扭转变形等。
在进行非线性分析时，需要考虑这些几何非线性因素对结构的影响。
3. 接触非线性钢筋混凝土板在使用过程中会受到多种荷载的作用，其中包括接触荷载。
接触非线性是指结构中两个或多个体之间的接触面会发生变形，从而影响结构的力学性能。
在进行非线性分析时，需要考虑接触非线性对结构的影响。
以上三个方面是钢筋混凝土板非线性分析的关键因素，下面将对其进行详细介绍。
1. 材料的本构关系混凝土的本构关系可以用双曲正切模型或Drucker-Prager模型等进行描述。
其中，双曲正切模型是一种常用的混凝土本构模型，其本构方程如下：σ = f(ε) = σc + α(ε-εc) + β(ε-εc)/(1+(ε-εc)/γ)其中，σ为混凝土的应力，ε为混凝土的应变，σc和εc分别为混凝土的极限应力和极限应变，α、β和γ为模型参数。
该模型可以较好地描述混凝土的非线性本构关系。
钢筋的本构关系可以采用弹塑性模型或Ramberg-Osgood模型等进行描述。
混凝土受压应力-应变全曲线方程(描述)
混凝土受压应力－应变全曲线方程混凝土受压应力－应变全曲线方程混凝土的应力－应变关系是钢筋混凝土构件强度计算、超静定结构内力分析、结构延性计算和钢筋混凝土有限元分析的基础，几十年来，人们作了广泛的努力，研究混凝土受压应力－应变关系的非线性性质，探讨应力与应变之间合理的数学表达式，1942年，Whitney通过混凝土圆柱体轴压试验，提出了混凝土受压完整的应力应变全曲线数学表达式，得出了混凝土脆性破坏主要是由于试验机刚度不足造成的重要结论，这一结论于1948年由Ramaley和Mchenry的试验研究再次证实，1962年，Barnard在专门设计的具有较好刚性且能控制应变速度的试验机上，试验了一批棱柱体试件以及试件两靖被放大的圆柱体试件，试验再次证明，混凝土的突然破坏并非混凝土固有特性，而是试验条件的结果，即混凝土的脆性破坏可用刚性试验机予以防止，后来由很多学者（如M.Sagin，P.T.Wang，过镇海等）所进行的试验，都证明混凝土受压应力－应变曲线确实有下降段存在，那么混凝土受压应力与应变间的数学关系在下降段也必然存在，研究这一数学关系的工作一刻也没有停止。
钢筋混凝土结构是目前使用最为广泛的一种结构形式。
但是，对钢筋混凝土的力学性能还不能说已经有了全面的掌握。
近年来，随着有限元数值方法的发展和计算机技术的进步，人们已经可以利用钢筋混凝土有限元分析方法对混凝土结构作比较精确的分析了。
由于混凝土材料性质的复杂性，对混凝土结构进行有限元分析还存在不少困难，其中符合实际的混凝土应力应变全曲线的确定就是一个重要的方面。
1、混凝土单轴受压全曲线的几何特点经过对混凝土单轴受压变形的大量试验大家一致公认混凝土单轴受压变过程的应力应变全曲线的形状有一定的特征。
典型的曲线如图1所示，图中采用无量纲坐标。
sc c E E N f y x 0,,===σεε 式中，c f 为混凝土抗压强度；c ε为与c f 对应的峰值应变；0E 为混凝土的初始弹性模量；s E 为峰值应力处的割线模量。
混凝土结构非线性分析与设计
混凝土结构非线性分析与设计混凝土结构是工业建筑和民用建筑中最常见的建筑结构之一。
在设计和分析混凝土结构时，通常需要考虑结构的稳定性、刚度、承载能力等方面，以保证其能够满足工程使用的要求。
然而，混凝土结构是一个典型的非线性结构，这意味着结构在受外力作用下的变形和应力状态是非线性的，而这种非线性往往会导致结构的预测性能和可靠性变差。
因此，对于混凝土结构的非线性分析与设计来说，更具有挑战性。
本文将详细探讨混凝土结构非线性分析与设计的相关问题。
混凝土结构的非线性特性混凝土结构存在多种非线性性质，这些非线性特性主要包括以下几点：1. 弹塑性行为混凝土结构在小变形范围内表现出弹性行为，但随着外部载荷的增加，混凝土就开始表现出一定程度的塑性行为。
这是因为混凝土的本身是一个复合材料，内部含有许多孔隙和缺陷，难以表现出完全的线性弹性特性。
2. 非线性材料特性与金属材料相比，混凝土是一种非常脆弱的材料，其受压性能强于其受拉性能。
在混凝土的受拉过程中，裂缝会逐渐扩展。
这种非线性行为会对混凝土结构的整体性能产生很大影响。
3. 多项式应力应变关系当混凝土受到剪切力或扭矩时，其应变和应力之间的关系并不是线性的，而是一个多项式函数。
这种关系导致混凝土结构的非线性行为更加复杂。
混凝土结构的非线性分析混凝土结构的非线性分析方法有多种，主要包括有限元法、弹塑性分析法、极限等效塑性法等。
有限元法是一种非线性分析方法中应用最为广泛的方法之一。
它通过将结构划分成很多的有限元单元，在各个单元上建立力学方程，求解整个结构的应力和位移场。
这种方法可以考虑裂缝的产生和扩展等非线性因素，因而能够比较准确地模拟混凝土结构的非线性行为。
弹塑性分析法是在有限元法的基础上发展起来的一种方法。
它将结构划分为弹性区和塑性区，并同时考虑两个区域的力学行为。
弹性区的行为遵循线性弹性理论，而塑性区的行为则遵循材料的应力-应变曲线。
这种方法适用于不同类型的混凝土结构，并可以将结构的裂缝和塑性变形等非线性因素考虑在内。
混凝土受压应力-应变全曲线方程(描述)
混凝土受压应力－应变全曲线方程混凝土受压应力－应变全曲线方程混凝土的应力－应变关系是钢筋混凝土构件强度计算、超静定结构内力分析、结构延性计算和钢筋混凝土有限元分析的基础，几十年来，人们作了广泛的努力，研究混凝土受压应力－应变关系的非线性性质，探讨应力与应变之间合理的数学表达式，1942年，Whitney 通过混凝土圆柱体轴压试验，提出了混凝土受压完整的应力应变全曲线数学表达式，得出了混凝土脆性破坏主要是由于试验机刚度不足造成的重要结论，这一结论于1948年由Ramaley 和Mchenry 的试验研究再次证实，1962年，Barnard 在专门设计的具有较好刚性且能控制应变速度的试验机上，试验了一批棱柱体试件以及试件两靖被放大的圆柱体试件，试验再次证明，混凝土的突然破坏并非混凝土固有特性，而是试验条件的结果，即混凝土的脆性破坏可用刚性试验机予以防止，后来由很多学者（如M.Sagin ，P.T.Wang ，过镇海等）所进行的试验，都证明混凝土受压应力－应变曲线确实有下降段存在，那么混凝土受压应力与应变间的数学关系在下降段也必然存在，研究这一数学关系的工作一刻也没有停止。
钢筋混凝土结构是目前使用最为广泛的一种结构形式。
但是，对钢筋混凝土的力学性能还不能说已经有了全面的掌握。
近年来，随着有限元数值方法的发展和计算机技术的进步，人们已经可以利用钢筋混凝土有限元分析方法对混凝土结构作比较精确的分析了。
由于混凝土材料性质的复杂性，对混凝土结构进行有限元分析还存在不少困难，其中符合实际的混凝土应力应变全曲线的确定就是一个重要的方面。
1、混凝土单轴受压全曲线的几何特点经过对混凝土单轴受压变形的大量试验大家一致公认混凝土单轴受压变过程的应力应变全曲线的形状有一定的特征。
典型的曲线如图1所示，图中采用无量纲坐标。
sc c E E N f y x 0,,===σεε 式中，c f 为混凝土抗压强度；c ε为与c f 对应的峰值应变；0E 为混凝土的初始弹性模量；s E 为峰值应力处的割线模量。
混凝土的应力应变关系及其分析方法
混凝土的应力应变关系及其分析方法混凝土是一种广泛应用于建筑工程中的材料，其工程性能与强度密切相关。
了解混凝土在受力下的应变变化特征，可以有效地指导混凝土结构的设计和施工过程。
本文将就混凝土的应力应变关系及其分析方法进行探讨和介绍。
一、混凝土的应力应变关系混凝土在受力下的应变变化特征与其材料性质、构造和外部荷载等因素密切相关。
在混凝土受力过程中，其应力应变关系通常分为弹性阶段和塑性阶段两个阶段来进行研究。
1. 弹性阶段在混凝土受力时，施加在其表面的应力随之产生应变。
当荷载较小时，混凝土会在受力后立即回弹并恢复初始状态，这一阶段称为弹性阶段。
在弹性阶段，混凝土的应变与应力成正比，即应力-应变曲线为一条直线。
这种情况下，混凝土的弹性模量可以用来表征其弹性性能。
弹性模量取决于混凝土的配合比、孔隙率、龄期等因素，其值一般在30~40GPa之间。
2. 塑性阶段当混凝土受到更大的荷载时，超过了其弹性极限，就会进入塑性阶段。
在这个阶段中，混凝土会先出现一定程度的塑性变形，然后在荷载升高的情况下继续变形，最后极限荷载达到时发生破坏。
在塑性阶段中，混凝土的应力-应变曲线不再是一条直线，而呈现出拐点和曲线段落。
混凝土的应变变化主要表现为体积变化和剪切变形。
这时，我们需要使用一些塑性力学理论来分析混凝土在受力过程中的变形特征。
二、混凝土应力应变关系的分析方法了解混凝土在受力下的应力应变关系对于工程设计和施工至关重要。
下面我们将介绍一些目前常用的分析方法。
1. 材料试验法材料试验法是通过试验的方式确定混凝土的应力应变特性。
通过制作不同尺寸规格的混凝土试样，在规定的试验条件下进行荷载试验，并记录荷载与应变的关系。
在试验中，我们可以得到混凝土的应力-应变曲线。
通过分析应力-应变曲线，我们可以知道混凝土的弹性模量、弹性极限、屈服强度、极限强度等指标，从而为工程设计提供数据支持。
2. 数值模拟法数值模拟法基于有限元分析原理，将复杂的结构体系离散化成若干个单元，进而分析其应力应变特性。
混凝土的应力-应变关系分析
混凝土的应力-应变关系分析一、引言混凝土是一种广泛使用的建筑材料，用于各种类型的建筑和基础工程。
混凝土的应力-应变关系是混凝土工程设计和结构分析中非常重要的一个因素。
本文将详细分析混凝土的应力-应变关系，包括混凝土的力学性质、应力-应变曲线的形状和特点、影响应力-应变关系的因素以及实验方法。
二、混凝土的力学性质混凝土是一种复合材料，由水泥、骨料、砂和水等组成。
混凝土的力学性质受到其组成和制备方法的影响。
混凝土的力学性质包括弹性模量、抗拉强度、抗压强度、剪切强度等。
1. 弹性模量混凝土的弹性模量是指在弹性阶段，混凝土的应变与应力之比。
弹性模量是混凝土的刚度指标，通常用于计算混凝土结构的变形和挠度。
混凝土的弹性模量通常介于20-40 GPa之间，取决于混凝土的成分和强度等级。
2. 抗拉强度混凝土的抗拉强度通常比抗压强度低很多。
这是因为混凝土的骨料在混凝土中的分布不均匀，导致混凝土在拉伸过程中难以传递应力。
混凝土的抗拉强度通常介于2-10 MPa之间。
3. 抗压强度混凝土的抗压强度是指混凝土在压缩过程中的最大承载能力。
混凝土的抗压强度通常是设计混凝土结构时最关键的性质之一。
混凝土的抗压强度通常介于10-50 MPa之间。
4. 剪切强度混凝土的剪切强度通常比抗压强度低很多。
这是因为混凝土在剪切过程中容易出现裂缝，导致混凝土的强度降低。
混凝土的剪切强度通常介于0.2-0.5 MPa之间。
三、应力-应变曲线的形状和特点混凝土的应力-应变曲线通常具有非线性的形状。
在应力较小的情况下，混凝土的应变与应力呈线性关系。
然而，随着应力的增加，混凝土开始发生非线性变形。
在一定应力范围内，混凝土的应力-应变曲线呈现出一个明显的拐点，称为峰值点。
在峰值点之后，混凝土开始出现裂缝和破坏，应力开始降低。
在应变较大的情况下，混凝土的应力与应变之间呈现出一个平台，称为残余强度。
混凝土的应力-应变曲线的形状和特点受到许多因素的影响，包括混凝土的强度等级、骨料类型和分布、水胶比、养护条件等。
混凝土材料的应力-应变特性原理
混凝土材料的应力-应变特性原理一、前言混凝土是一种常用的建筑材料，在现代建筑中得到广泛的应用。
混凝土的应力-应变特性是混凝土材料的重要性能之一，是混凝土结构设计的基础。
本文将对混凝土材料的应力-应变特性进行详细介绍。
二、混凝土的应力-应变曲线混凝土材料的应力-应变特性通常是用应力-应变曲线来表示。
应力-应变曲线可以反映混凝土材料的强度、韧性和变形性能等特性。
1. 应力-应变曲线的基本形态应力-应变曲线的基本形态如图1所示。
曲线的第一段是线性段，称为弹性阶段；第二段是非线性段，称为塑性阶段；第三段是断裂阶段，称为破坏阶段。
图1 应力-应变曲线的基本形态2. 弹性阶段弹性阶段是应力-应变曲线的线性段，其斜率称为弹性模量。
在弹性阶段，混凝土材料的应变与应力成正比，而且在去除载荷后，混凝土材料完全恢复原来的形态。
3. 塑性阶段塑性阶段是应力-应变曲线的非线性段，也称为屈服阶段。
在这个阶段，混凝土材料开始发生塑性变形，应力-应变曲线的斜率开始减小。
在这个阶段，混凝土材料的应变增加，但应力增加的速率减慢。
4. 破坏阶段破坏阶段是应力-应变曲线的最后一段，也称为断裂阶段。
在这个阶段，混凝土材料的应力急剧下降，出现明显的裂纹和破坏。
在这个阶段，混凝土材料已经失去了承载能力。
三、混凝土的应力-应变特性的影响因素混凝土的应力-应变特性受到许多因素的影响，包括混凝土材料的成分、制备工艺、试验条件等。
1. 混凝土材料的成分混凝土材料的成分是影响其应力-应变特性的重要因素之一。
常见的混凝土材料成分包括水泥、骨料、粉煤灰、膨胀剂等。
其中，水泥的种类、含量和水灰比对混凝土的强度和变形性能有很大的影响。
2. 制备工艺混凝土的制备工艺也会影响其应力-应变特性。
制备工艺包括搅拌时间、搅拌方式、养护方式等。
其中，搅拌时间和搅拌方式对混凝土的均匀性和孔隙度有影响，养护方式对混凝土的强度和变形性能有影响。
3. 试验条件试验条件也会影响混凝土的应力-应变特性。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
各位同学，这里几乎已经是最完整的命令流了，关于这方面工作，初步建议如下：（1）初步的模拟最好先不要考虑混凝土下降段，不要考虑混凝土的压碎，同时单元尺寸不要过小，以免影响收敛。
（2）反复试算结果表明，混凝土的MISO和KINH模型、DP和MISES准则、剪力传递系数对模拟结果影响均不大，但在一定程度上影响了收敛性的好坏。
（3）当解决了收敛问题之后，最大的难点是下降段，对于一个桥墩（柱）的反复荷载试验，引起承载力下降的原因大体可归结为：塑性铰区混凝土压碎脱落、纵筋屈曲，对于以剪切破坏为主的试件，还包括剪切滑移的影响。
对于以黏结破坏为主的试件，还包括纵筋－混凝土黏结滑移破坏。
因此，能否正确模拟这些将是模型正确与否的关键。
众所周知，以link8单元模拟钢筋自然无法模拟纵筋屈曲，这样就把一个因素忽略了。
（4）作者们最后采用了扩大的破坏面，坦率的讲，这也是为了模拟混凝土的压碎破坏而采取的没有办法的办法，鉴于作者水平有限，至今也无法对此做出更为合理的解释。
破坏面到底扩大多少，仍无规律可循。
且一旦打开压碎开关，会出现诸如难以收敛、“假压碎”等一系列问题。
所以我们建议大家起初还是先不要考虑压碎为好。
定义应力应变曲线1，定义变量：拾取主菜单：Main Menu>Time Hist postproc>Define Variables>在随之弹出的对话框中点击Add键，定义第一个变量序号为2，选取第一个变量stress,确定与之对应的下一级选项（如Y-direction SY等）；返回定义变量对话框，再点击add键，定义第二个变量序号为3，选取第二个变量strain-elastic及以及对赢得下一级选项（如Y-dir'n EPEL Y等，在应力-应变图中，其向量的取向应相同）。
同理再定义变量4，选取变量strain-plastic及与之对应的下一级选项如Y-dir'n EPEL Y等），在应力-应变图中，应变是弹性应变和塑性应变累加的总应变。
为使其实现相加，还需进行以下操作：拾取主菜单：Main Menu>Time Hist postproc>math operation>add,定义计算变量序号为5，同时在相应交互框内输入3和4。
点击确认键，则由变量3，4代表的应变之和就存在变量5中。
2，绘制应力-应变曲线：拾取主菜单：Main Menu>Time Hist postproc>setting>graph.设置x轴向变量为单变量，并将其变量序号定义为5。
点击确定键退出退化框。
拾取应用菜单：Utility Menu>plot ctrls>styles>Graphs>Modify axis.将x,y坐标轴分别命名为Y-strains,Y-stress,拾取主菜单：Main Menu>Time Hist postproc>graph variables. 在对话框上"the first variable"对应的交互框中输入2。
点击确定键，则预想的应力-应变曲线就显示在屏幕上。
ok!试试看！欢迎大家继续批评指正！/prep7et,1,solid65et,2,link8et,3,solid45r,1r,2,236e-6r,3,75.4e-6r,4,151.2e-6!混凝土本构关系mp,ex,1,2.522e10mp,prxy,1,0.2mp,dens,1,2500tb,concr,1,1tbdata,,0.5,0.95,5e6,-1tb,miso,1,1,7tbtemp,0tbpt,,0.00068,17.15e6tbpt,,0.00136,28.03e6tbpt,,0.00204,33.62e6tbpt,,0.00272,35.96e6tbpt,,0.0034,36.53e6tbpt,,0.00408,36.53e6!纵筋的本构关系mp,ex,2,2.0e11mp,prxy,2,0.28tb,bkin,2tbtemp,0tbdata,,362e6,2e9！箍筋的本构关系mp,EX,3,2.1E11mp,PRXY,3,.27tb,bkin,3tbtemp,0tbdata,,272.5e6,2.1e9!建立模型local,11,1,0,0,0,0,-90,0 ！建立局部坐标系，采用柱坐标，绕y轴顺时针转90度csys,11*do,j,1,8n,j,0.185,j*360/8,-0.6*enddongen,3,100,1,8,1,0,0,0.3 ！复制节点，3次，编号增量，源节点1~8，0,0,0.3 ngen,4,100,201,208,1,0,0,0.15ngen,7,100,501,508,1,0,0,0.3type,2 ！指定单元属性，建立基础箍筋单元real,3mat,3*do,i,1,6*do,j,1,7e,(i-1)*100+j,(i-1)*100+j+1*enddoe,(i-1)*100+1,(i-1)*100+8*enddotype,2 ！建立柱上箍筋单元real,4mat,3*do,i,1,5*do,j,1,7e,(i+5)*100+j,(i+5)*100+j+1*enddoe,(i+5)*100+1,(i+5)*100+8*enddotype,2 ！建立纵向钢筋real,2mat,2*do,i,1,11*do,j,1,8e,(i-1)*100+j,i*100+j*enddo*enddoallsel,all ！合并节点，压缩节点nummrg,allnumcmp,allwpcsys,1,11cylind,0,0.185,-0.6,0,0,360 ！建立圆柱，内外径i-or，高度z1-z2，起终角度cylind,0,0.185,0,0.45,0,360cylind,0.185,0.3,-0.6,0,0,360cylind,0,0.185,0.45,2.25,0,360cylind,0,0.185,2.25,2.55,0,360cylind,0.185,0.20,2.25,2.55,0,360cylind,0,0.185,2.55,2.565,0,360wprota,,90 ！绕y轴逆时针转90度vsel,s,,,1,7,1 ！选择体，选择类型，item，comb，最小编号，最大编号，增量vsbw,all ！体相减numcmp,allwprota,,,90 ！绕z轴逆时针转90度vsel,s,,,1,14,1vsbw,allnumcmp,allvsel,s,,,1,28,1vglue,allnumcmp,allwprota,,,-90wprota,,-90allsel,allvsel,s,,,1,16,1vsel,a,,,21,24,1vatt,1,1,1vsel,s,,,17,20,1vsel,a,,,25,28,1vatt,2,1,3allsel,alllsel,s,,,1,8,1 ！选择线lsel,a,,,11,22,1 ！补选（a）线lsel,a,,,25,40,1lsel,a,,,43,46,1lesize,all,,,2 ！线尺寸为2lsel,s,,,47,50,1lsel,a,,,53,56,1lsel,a,,,58,73,1lesize,all,,,2lsel,s,,,98,107,1lsel,a,,,110,112,2lesize,all,,,1lsel,s,,,23,24,1lsel,a,,,9,10,1lesize,all,,,2lsel,s,,,78,82,1lesize,all,,,3lsel,s,,,83,87,1lesize,all,,,6lsel,s,,,41,42,1lsel,a,,,88,92,1lesize,all,,,2lsel,s,,,93,97,1lesize,all,,,1mshkey,1mshape,0,3dvsel,s,,,1,28vmesh,allallsel,allnummrg,allnumcmp,allwprota,,90csys,0allsel,all/solu ！进入加载、求解csys,11nsel,s,loc,x,0.3 ！选择半径在x=3处的节点，并施加约束d,all,allallsel,allOUTRES,ALL,2 ！输出获得的解csys,0pred,on ！打开预应力选项cnvtol,f,,0.05,2 !定义收敛条件,使用缺省的VALUE autots,1 !打开自动时间步控制lnsrch,1 !打开线性搜索ncnv,2 !如果不收敛时结束而不退出neqit,50 !每一子步中方程的迭代次数限值savetime,1 ！定义第1载荷步nsubst,10asel,s,,,29,32,1sfa,all,,pres,5468699 ！在已选择的面上施加压力荷载allsel,allacel,,9.8 ！重力加速度y=9.8kbc,1 ！阶跃加载方式lswrite,1time,10 ！定义第2荷载步cp,1,ux,298,300,303 !将节点298,300,303ux向位移耦合d,298,ux,0.015 ！再借点298处施加ux位移0.015 allsel,allKBC,0allsel,allNSUBST,20,1000000,10lswrite,2time,20 ！定义第3载荷步d,298,ux,-0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,3time,30 ！定义第4载荷步d,298,ux,0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,4time,40 ！定义第5载荷步d,298,ux,-0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,5time,50 ！定义第6载荷步d,298,ux,0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,6time,60 ！定义第7载荷步d,298,ux,-0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,7time,70 ！定义第8载荷步d,298,ux,0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,8time,80d,298,ux,-0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,9time,90d,298,ux,0.015allsel,allKBC,0allsel,allNSUBST,100,1000000,10lswrite,10time,100d,298,ux,-0.015 allsel,allKBC,0allsel,allNSUBST,100,1000000,10 lswrite,11time,110d,298,ux,0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10 lswrite,12time,120d,298,ux,-0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10 lswrite,13time,130d,298,ux,0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10 lswrite,14time,140d,298,ux,-0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10 lswrite,15time,150d,298,ux,0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10 lswrite,16time,160d,298,ux,-0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10 lswrite,17time,170d,298,ux,0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10 lswrite,18time,180d,298,ux,-0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10 lswrite,19time,190d,298,ux,0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10 lswrite,20time,200d,298,ux,-0.03allsel,allKBC,0allsel,allNSUBST,200,1000000,10 lswrite,21time,210d,298,ux,0.045allsel,allKBC,0allsel,allNSUBST,300,1000000,10lswrite,22lssolve,1,22,1命令流：1、单调水平荷载作用下的命令流：建模：!mono brick wall!solid65-整体式建模!*Finish/clear/units,si/com,Structural/config,nres,400000/prep7/title,mono-brick-wall!*depth_sup=0.24 !弹性支垫240mm!单元属性!for brick, mortar and leveling coatet,1,solid65!*KEYOPT,1,1,0KEYOPT,1,5,0KEYOPT,1,6,0KEYOPT,1,7,1!for elastic supportet,2,solid45!*--------------------------------砌体材料-------------------------------------- !MU10,M5砌体材料属性mp,prxy,1,0.15mp,dens,1,1700mp,ex,1,0.24e10!砌体屈服准则TB,MKIN,1TBTEMP,,strainTBdata,1,0.2E-3,1.0E-3,2.0E-3,3.0E-3,4.8E-3TBTEMP,,,!TBDATA,1,0.48E6,0.82E6,1.32E6,1.5E6,1.32E6TBDATA,1,0.48E6,1.31E6,2.11E6,2.4E6,2.11E6!/XRANGE,0,0.02!TBPLOT,MKIN,1!砌体破坏准则!hntrl=0.13E6hntrl=0.21E6tb,concr,1,,4,tbdata,,0.2,0.9,hntrl,-1 !不考虑压碎!*--------------------------------弹性支撑-------------------------------------- !弹性支撑暂取钢筋材料属性mp,ex,2,2.0e11mp,dens,2,7.600mp,prxy,2,0.30!real constantsr,1!*--------------------------------modeling-------------------------------------- !geometry modelblock,0,3.12,0,0.24,0,3.12block,-0.24,3.12,0,0.24,3.12,3.36Vsel,s,volu,,1,2,1$aslv,s$lsla,slesize,all,0.24Vsel,s,volu,,1$vatt,1,1,1, !brick wallVsel,s,volu,,2$vatt,2,1,2, !brick wallvsel,all !必须再次选择所有单元VSWEEP,ALLsavenumcmp,allallsel!*/device,vector,1/eshape,1!*-----------------------------------耦合--------------------------------------- allsel,all!elastic support nodes setesel,s,type,,2$nsle,s,all$nsel,r,loc,z,3.12nsel,u,loc,x,-0.24$cm,slavenode,node!other node setesel,s,type,,1$nsle,s,all$cm,masternode,node/input,ucouple,macfinish耦合文件ucouple,mac：!*-----------------------------------开始--------------------------------------- allsel !最好保留这句命令!*******将从属节点编号依次存入数组**************** cmsel,s,slavenode*get,count1_node,node,0,count*del,slave_node*dim,slave_node,array,count1_node*get,slave_node(1),node,0,num,min*do,i,2,count1_nodeslave_node(i)=ndnext(slave_node(i-1))*enddo!*******将主节点编号依次存入数组****************allselcmsel,s,masternode*get,count2_node,node,0,count*del,master_node*dim,master_node,array,count2_node*get,master_node(1),node,0,num,min*do,i,2,count2_nodemaster_node(i)=ndnext(master_node(i-1))*enddo!********将与从属节点耦合的节点数组初始化**************** *del,cp_node*dim,cp_node,array,count1_node*do,i,1,count1_nodecp_node(i)=0*enddo!*********开始选择程序****************allselcmsel,s,masternode*do,i,1,count1_nodekk=1k=1*dowhile,kkk=nnear(slave_node(i))nsel,s,cp,,allkk=nsel(k)+0.001cmsel,s,masternodensel,u,node,,kcm,masternode,node*enddocp_node(i)=k*enddo!*******选择完毕****************!*******开始耦合****************allsel,all/prep7*do,i,1,count1_nodecp,next,UX,slave_node(i),cp_node(i)cp,next,UY,slave_node(i),cp_node(i)cp,next,UZ,slave_node(i),cp_node(i)*enddo!*******耦合完毕****************求解方法NR法：brick wall NR solution!displacement load!用NR法!*---------------------------------solution------------------------------------- /soluantype,0!-----------constraint----------------nsel,s,loc,z,0d,all,all!*esel,s,type,,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uy!*!esel,s,type,,2$nsle,s,all$nsel,r,loc,x,-0.24!d,all,uzlumpm,onsolcontrol,onpred,onlnsrch,onautots,1 !打开自动时间步控制ncnv,2 !如果不收敛时结束而不退出compression=0.24E7!*------------------------------------load steps-------------------------------- !*------------------------------------1st--------------------------------------- acel,,,9.8nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.119esln,s,0,allsfe,all,6,pres,0,0.03*compressionallsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,2,200,1neqit,15KBC,0nlgeom,1cnvtol,f,,0.01,2,0.05allsellswrite,1!*--------------------------------2nd step-------------------------------------- nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6!d,all,ux,2.79E-3d,all,ux,1.0E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,100neqit,15KBC,0nlgeom,1cnvtol,f,,0.05,2,0.05allsellswrite,2!*--------------------------------3rd step-------------------------------------- nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6!d,all,ux,2.79E-3d,all,ux,2.5E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,100neqit,15KBC,0nlgeom,1cnvtol,f,,0.09,2,0.25allsellswrite,3!*--------------------------------4th step-------------------------------------- nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6!d,all,ux,2.79E-3d,all,ux,4.0E-3allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,100neqit,15KBC,0nlgeom,1cnvtol,f,,0.1,2,0.5allsellswrite,4!*allsellssolve,1,4,1FINISH也可以用弧长法，不过没有一定经验这比较难控制：!mono brick wall-1 arclength solution!displacement load!*---------------------------------solution------------------------------------- /soluantype,0!-----------constraint----------------nsel,s,loc,z,0d,all,all!*esel,s,type,,2$nsle,s,all$nsel,r,loc,x,-0.24 d,all,uy!*lumpm,onsolcontrol,onncnv,2 !如果不收敛时结束而不退出compression=0.24E7!*--------------load steps----------------!*--------------first step----------------- acel,,,9.8nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.12 esln,s,0,allsfe,all,6,pres,0,0.03*compression allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,,1E-5, !使用弧长法NSUBST,200neqit,15KBC,0nlgeom,1cnvtol,f,,0.1,2,0.5allsel,alllswrite,1!*-------------second step---------------- esel,s,type,,2$nsle,s,all$nsel,r,loc,x,-0.24 d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6d,all,ux,2.50E-3!d,all,ux,3.51E-3allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,,1E-5, !使用弧长法kbc,0NSUBST,250!arctrm,u,0.0028,313,uxneqit,15nlgeom,1!cnvtol,f,,0.005,2,0.01!cnvtol,f,,0.05,2,0.5cnvtol,f,,0.1,2,0.5allsellswrite,2!*-------------third step--------------- esel,s,type,,2$nsle,s,all$nsel,r,loc,x,-0.24 d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6!d,all,ux,2.79E-3!d,all,ux,3.18E-3d,all,ux,3.00E-3allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,,1E-7, !使用弧长法kbc,0NSUBST,2000!arctrm,u,0.0028,313,uxneqit,15nlgeom,1!cnvtol,f,,0.005,2,0.01!cnvtol,f,,0.05,2,0.5cnvtol,f,,0.5,2,0.8allsellswrite,3allsel!*-------------fourth step--------------- esel,s,type,,2$nsle,s,all$nsel,r,loc,x,-0.24 d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6d,all,ux,3.51E-3allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLarclen,on,,1E-7, !使用弧长法NSUBST,300!arctrm,u,0.0028,313,uxneqit,15nlgeom,1!cnvtol,f,,0.005,2,0.01!cnvtol,f,,0.05,2,0.5cnvtol,f,,0.6,2,0.9allsellswrite,4allsel!*lssolve,1,4,1FINISH后处理：/post26esel,s,type,,1eplotlines,5000nsol,2,326,u,x,ux!nsol,2,313,u,x,ux*do,i,1,14jj=node((i-1)*0.24,0.24,0)kk=node((i-1)*0.24,0,0)! ll=node(225.5+(i-1)*225.5,3300,0) rforce,3,jj,f,x,fxadd,6,3,6rforce,4,kk,f,x,fxadd,6,4,6! rforce,5,ll,f,x,fx!add,6,5,6*enddoprod,7,2,,,,,,1000prod,8,6,,,,,,-1/1000/axlab,x,Displacement(mm)/axlab,y,F(kN)/XRANGE,0,4/YRANGE,0,100xvar,7!rforce,5,340,f,x,fx!esol,5,87,221,f,x,fx2、循环水平荷载作用下的命令流：出求解文件不同外，其他都相同；!cyclic NR solution!displacement load!*---------------------------------solution------------------------------------- /soluantype,0!-----------constraint----------------nsel,s,loc,z,0d,all,all!*esel,s,type,,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uy!*esel,s,type,,2$nsle,s,all$nsel,r,loc,x,3.36d,all,uy!*lumpm,onsolcontrol,onpred,onlnsrch,onautots,1 !打开自动时间步控制ncnv,2 !如果不收敛时结束而不退出compression=0.24E7!*----------------------------load steps---------------------------------------- !*--------------------------------1st step-------------------------------------- acel,,,9.8nsel,s,loc,z,3.36$nsel,r,loc,x,0.001,3.119esln,s,0,allsfe,all,6,pres,0,0.03*compressionallsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,10neqit,15nlgeom,1cnvtol,f,,0.1,2,0.5allsel,alllswrite,1!*----------------------------初裂荷载20%做1次循环------------------------------ !*---------------------------------用NR法---------------------------------------!*--------------------------------2nd step-------------------------------------- esel,s,type,,2$nsle,s,all$nsel,r,loc,x,-0.24d,all,uz!*esel,s,type,,2$nsle,s,all$nsel,r,loc,x,3.36d,all,uz!*nsel,s,loc,x,-0.24!esln,s,0,all!sfe,all,4,pres,0,60*0.0593E6d,all,ux,1.48E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,20,2000,10neqit,15NROPT ,FULL,,OFFKBC,0nlgeom,1cnvtol,f,,0.1,2,0.5allsellswrite,2!*--------------------------------3rd step-------------------------------------- nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-1.48E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,40,2000,20neqit,15NROPT ,FULL,,OFFKBC,0nlgeom,1cnvtol,f,,0.1,2,0.5allsellswrite,3!*--------------------------------4th step-------------------------------------- nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,1.48E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,40,2000,20neqit,15NROPT ,FULL,,OFFKBC,0nlgeom,1cnvtol,f,,0.1,2,0.5allsellswrite,4!*-----------------------------------one cycle----------------------------------!*----------------------------初裂荷载50%做1次循环------------------------------ !*---------------------------------用NR法---------------------------------------!*--------------------------------5th step-------------------------------------- nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,2.85E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,100,2000,10neqit,15NROPT ,FULL,,OFFKBC,0nlgeom,1cnvtol,f,,0.1,2,0.5allsellswrite,5!*--------------------------------6th step-------------------------------------- nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-2.85E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,30neqit,15NROPT ,FULL,,OFFKBC,0nlgeom,1cnvtol,f,,0.1,2,0.5allsellswrite,6!*--------------------------------7th step-------------------------------------- nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,2.85E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALLNSUBST,200,2000,30neqit,15NROPT ,FULL,,OFFKBC,0nlgeom,1cnvtol,f,,0.1,2,0.5allsellswrite,7!*------------------------------------two cycle---------------------------------!*----------------------------初裂荷载80%做1次循环------------------------------ !*---------------------------------用NR法---------------------------------------!*---------------------------------8th step------------------------------------- nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,3.5E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,100,2000,20neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,8!*---------------------------------9th step------------------------------------- nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-3.5E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,9!*---------------------------------10th step------------------------------------- nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,3.5E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,10!*----------------------------------three cycle---------------------------------!*----------------------------初裂荷载100%做1次循环----------------------------- !*---------------------------------用NR法---------------------------------------!*--------------------------------11th step------------------------------------- nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,4.42E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,100,2000,20neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,11!*--------------------------------12th step------------------------------------- nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-4.42E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,12!*--------------------------------13th step------------------------------------- nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,4.42E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,13!*-----------------------------------four cycle---------------------------------!*----------------------------61.4KN的荷载做1次循环----------------------------- !*---------------------------------用NR法---------------------------------------!*--------------------------------14th step------------------------------------- nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,6.72E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,100,2000,20neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,14!*--------------------------------15th step------------------------------------- nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-6.72E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,15!*--------------------------------16th step------------------------------------- nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,6.72E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,16!*-----------------------------------five cycle---------------------------------!*----------------------------74.1KN的荷载做1次循环----------------------------- !*---------------------------------用NR法---------------------------------------!*--------------------------------17th step------------------------------------- nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,9.01E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,100,2000,20neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,17!*--------------------------------18th step------------------------------------- nsel,s,loc,x,-0.24ddele,all,uxnsel,s,loc,x,3.36d,all,ux,-9.01E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,18!*--------------------------------19th step------------------------------------- nsel,s,loc,x,3.36ddele,all,uxnsel,s,loc,x,-0.24d,all,ux,9.01E-4allsel,all!*OUTRES,ERASEOUTRES,ALL,ALL!arclen,on,,1E-7, !使用弧长法kbc,0NSUBST,200,2000,30neqit,15KBC,0nlgeom,1cnvtol,f,,0.5,2,0.9allsellswrite,19!*------------------------------------six cycle--------------------------------- allsellssolve,1,19,1FINISH混凝土也好，砌体也好，开裂之后ANSYS的弧长法很难计算，特别麻烦，这时候只要不是计算下降段最好还是用NR法。