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)模型。

5.方茴说：“那时候我们不说爱，爱是多么遥远、多么沉重的字眼啊。

我们只说喜欢，就算喜欢也是偷偷摸摸的。

”6.方茴说：“我觉得之所以说相见不如怀念，是因为相见只能让人在现实面前无奈地哀悼伤痛，而怀念却可以把已经注定的谎言变成童话。

”7.在村头有一截巨大的雷击木，直径十几米，此时主干上唯一的柳条已经在朝霞中掩去了莹光，变得普普通通了。

8.这些孩子都很活泼与好动，即便吃饭时也都不太老实，不少人抱着陶碗从自家出来，凑到了一起。

9.石村周围草木丰茂，猛兽众多，可守着大山，村人的食物相对来说却算不上丰盛，只是一些粗麦饼、野果以及孩子们碗中少量的肉食。

遗憾，每个遗憾都有它的青春美。

4.方茴说：“可能人总有点什么事，是想忘也忘不了的。

”5.方茴说：“那时候我们不说爱，爱是多么遥远、多么沉重的字眼啊。

我们只说喜欢，就算喜欢也是偷偷摸摸的。

”6.方茴说：“我觉得之所以说相见不如怀念，是因为相见只能让人在现实面前无奈地哀悼伤痛，而怀念却可以把已经注定的谎言变成童话。

”7.在村头有一截巨大的雷击木，直径十几米，此时主干上唯一的柳条已经在朝霞中掩去了莹光，变得普普通通了。

8.这些孩子都很活泼与好动，即便吃饭时也都不太老实，不少人抱着陶碗从自家出来，凑到了一起。

9.石村周围草木丰茂，猛兽众多，可守着大山，村人的食物相对来说却算不上丰盛，只是一些粗麦饼、野果以及孩子们碗中少量的肉食。

混凝土是目前应用最为广泛的建筑材料之一。

为了解混凝土结构的受力机理和破坏过程，在大型有限元软件ANSYS中，专门设置了Sdid65单元来模拟混凝土或钢筋混凝土结构，提供了很多缺省参数，从而为使用者提供了很大的方便。

1 Solid65单元Sdid65单元是专为混凝土、岩石等抗压能力远大于抗拉能力的非均匀材料开发的单元。

它可以模拟混凝土中的加强钢筋(或玻璃纤维、型钢等)，以及材料的拉裂和压溃现象。

1.1 几点假设1)只允许在每个积分点正交的方向开裂。

第一章ANSYS与结构分析1.1 ANSYS功能与软件结构工程和制造业的生命力在于产品的创新，而计算机的发展和广泛应用大大提高了产品开发、设计、分析和制造的效率和产品性能，用计算机对设计产品实时或进行随后的分析称为计算机辅助工程。

即CAE（Computer Aided Engineering）。

该技术是由计算机技术和工程分析技术相结合形成的新兴技术，它涉及计算力学、计算数学、结构动力学、数字仿真技术、工程管理学与计算机技术等学科。

随着有限元理论和计算机硬件的发展，CAJ软件和技术越来越成熟，已逐渐成为工程师实现工程创新和产品创新的得力助手和有效工具。

大型通用CAE软件可对多种类型功能和产品物理力学性能进行分析，其应用范围及其广泛，如ANSYS、ADINA、NASTRAN、MARC、ABAQUS、ADAMS、I-DEAS、SAP 等。

ANSYS软件是融结构、流体、电磁场、声场和热场分析于一体的大型大型通用有限元分析软件，可广泛应用于土木、地质、矿业、材料、机械、仪器仪表、热工电子、水利、生物医学和原子能等工程的分析和科学研究。

它可在大多数计算机和操作系统（如Windows、UNIX、Linux、HP-UX等）中运行，可与大多数CAD软件接口。

1970年，Dr.John.Swanson成立了Swanson Analysis System，Inc，后来重组后改称AN-SYS公司，总部设在美国宾西法尼亚州的匹兹堡。

近几年来，ANSYS 软件发展迅速，功能不断增强，目前最高版本为11.0beta。

1.1.1 ANSYS软件的技术特点ANSYS的主要技术特点如下：（1）强大的建模能力：仅靠ANSYS本身就可建立各种复杂的几何模型，可采用自底向上、自顶向下或两者混合建模方法，通过各种布尔运算和操作建立所需几何实体。

（2）强大的求解能力：ANSYS提供了数种求解器，主要类型有迭代求解器（预条件共轭梯度、雅可比共轭梯度、不完全共轭梯度），直接求解器（波前、稀疏矩阵）、特征值求解法（分块Lanczos法、子空间法、凝聚发、QR阻尼法）、并行求解器(分布式并行、代数多重网格)等,用户可根据问题类型选择合适的求解器。

用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结构分析基本原理1 应力-应变关系本文将介绍结构分析中材料线性理论，在后续再介绍材料非线性的理论。

在线弹性理论中应力-应变关系：(1) 其中：{σ}：应力分量，即在ANSYS软件里以S代替σ形式出现。

[D]：弹性矩阵或弹性刚度矩阵或应力-应变矩阵。

利用(14)～(19)给出了其具体表达式。

(4)给出了其逆矩阵的表达式。

通过给出完整的[D]可以定义少数的各向异性单元。

在ANSYS 中利用命令：TB,ANEL来输入具体数值。

：弹性应变矢量。

在ANSY中以EPEL形式输出。

{ε}：总的应变矢量，即{εth}：热应变矢量，(3)给出了其定义式，在ANSYS中以EPTH形式给出。

注意：{εel}：是由应力引起的应变。

软件中的剪切应变( εxy、εyz和εxz)是工程应变，他们是拉伸应变的两倍。

ε通常用来表示拉伸应变，但为了简化输出而采用此表示。

将在材料的非线性分析中说明总应变的分量，以EPTO形式输出。

图1 单元的应力矢量图如图1给出了单元应力矢量图。

ANSYS程序中规定正应力和正应变拉伸是为正，压缩时为负。

(1)式还可以被写作以下形式：(2)三维情况下，热应变矢量为：(3) 其中：：方向的正割热膨胀系数。

ΔT=T-T refT：问题中节点当前温度。

T ref：参考温度也就是应变自由时的温度。

用TREF或MP命令输入。

柔度矩阵的定义：(4) 其中：E x: 方向上的杨氏模量，在MP命令中用EX输入。

v xy:主泊松比，在MP命令中用PRXY输入。

v yx:次泊松比，在MP命令中用NUXY输入。

G xy: 平面上的剪切模量，在MP命令中用GXY输入。

此外，[D]-1是对称矩阵，因此(5)(6)(7)由(5)～(7)，可知νxy、νyz、νxz、νyx、νzy和νzx是不独立的，因此程序中必须输入νxy、νyz和νxz(以PRXY, PRYZ, and PRXZ标记输入)或νyx、νzy和νzx(以NUXY, NUYZ, and NUXZ标记输入)。

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结构分析的一些细节一将命令流代码录入文本文件中（如xx.txt），通过Utility>File>Read Input From命令读入文本文件，可迅速解题。

网格细化网格细化的命令在 Main Menu>reprocessor>Meshing>Modify Mesh菜单中（对应于Mesh Tool面板下方的网格细化操作按钮Refine）。

操作中有下拉列表框Level ofrefinement用来选择网格细化的精度。

细化单元操作仅仅对于平面的三角形、四边形网格以及体的四面体网格适用，对于使用六面体进行划分的三维几何模型不能进行网格细化。

耦合约束前处理器中的Coupling/Ceqn 选项用来定义耦合约束。

ANSYS中可以设置一种特殊的称为耦合的加载方式，一个耦合设置是一组被约束在一起，有着相同大小但数值未知的自由度，在铰的处理、接触分析等问题中往往需要用到耦合约束。

选择Main Menu>reprocessor>Couping/Ceqn>Couple DOFs命令。

生成耦合约束对应的命令CP，操作格式如下：CP，N（耦合编号），耦合自由度，节点1，节点2，节点3……施加荷载荷载可以加在实体模型或有限元模型（划分好网格的模型）上。

不管采取什么加载方式，ANSYS求解前都将荷载转化到有限元模型。

加载到实体模型上的荷载将自动转化到其所属的节点、单元上。

当删除实体模型时，ANSYS将自动删除其上所有荷载。

但建议用户直接将荷载加在实体模型上，因为在实体模型上加载操作简单（尤其是在直接拾取图形时），另一方面重新划分网格或局部修改网格时不对荷载产生影响。

在位移约束中有两个非常重要的约束方式：施加在对称面上的对称约束和反对称约束。

对称约束指限制对称面内所有节点的旋转自由度，同时限制垂直于对称面的位移自由度，对应着Main Menu>Solution>DefineLoads>Apply>Structural>Displacement>Symmetry B.C命令。

ANSYS结构分析基础篇一、总体介绍进行有限元分析的基本流程:1.分析前的思考1)采用哪种分析(静态,模态,动态...)2)模型是零件还是装配件(零件可以form a part形成装配件,有时为了划分六面体网格采用零件,但零件间需定义bond接触)3)单元类型选择(线单元,面单元还是实体单元)4)是否可以简化模型(如镜像对称,轴对称)2.预处理1)建立模型2)定义材料3)划分网格4)施加载荷及边界条件3.求解4.后处理1)查看结果(位移,应力,应变,支反力)2)根据标准规范评估结构的可靠性3)优化结构设计高阶篇:一、结构的离散化将结构或弹性体人为地划分成由有限个单元,并通过有限个节点相互连接的离散系统。

这一步要解决以下几个方面的问题:1、选择一个适当的参考系,既要考虑到工程设计习惯,又要照顾到建立模型的方便。

2、根据结构的特点,选择不同类型的单元。

对复合结构可能同时用到多种类型的单元,此时还需要考虑不同类型单元的连接处理等问题。

3、根据计算分析的精度、周期及费用等方面的要求,合理确定单元的尺寸和阶次。

4、根据工程需要,确定分析类型和计算工况。

要考虑参数区间及确定最危险工况等问题。

5、根据结构的实际支撑情况及受载状态,确定各工况的边界约束和有效计算载荷。

二、选择位移插值函数1、位移插值函数的要求在有限元法中通常选择多项式函数作为单元位移插值函数,并利用节点处的位移连续性条件,将位移插值函数整理成以下形函数矩阵与单元节点位移向量的乘积形式。

位移插值函数需要满足相容(协调)条件,采用多项式形式的位移插值函数,这一条件始终可以满足。

但近年来有人提出了一些新的位移插值函数,如:三角函数、样条函数及双曲函数等,此时需要检查是否满足相容条件。

2、位移插值函数的收敛性(完备性)要求:1) 位移插值函数必须包含常应变状态。

2)位移插值函数必须包含刚体位移。

3、复杂单元形函数的构造对于高阶复杂单元,利用节点处的位移连续性条件求解形函数,实际上是不可行的。

ANSYS结构分析报告教程篇ANSYS结构分析基础篇⼀、总体介绍进⾏有限元分析的基本流程：1.分析前的思考1)采⽤哪种分析(静态，模态，动态...）2)模型是零件还是装配件(零件可以form a part形成装配件，有时为了划分六⾯体⽹格采⽤零件，但零件间需定义bond接触)3)单元类型选择(线单元，⾯单元还是实体单元)4)是否可以简化模型(如镜像对称，轴对称)2.预处理1)建⽴模型2)定义材料3)划分⽹格4)施加载荷及边界条件3.求解4.后处理1)查看结果(位移，应⼒，应变，⽀反⼒)2)根据标准规范评估结构的可靠性3)优化结构设计⾼阶篇：⼀、结构的离散化将结构或弹性体⼈为地划分成由有限个单元，并通过有限个节点相互连接的离散系统。

这⼀步要解决以下⼏个⽅⾯的问题:1、选择⼀个适当的参考系，既要考虑到⼯程设计习惯，⼜要照顾到建⽴模型的⽅便。

2、根据结构的特点，选择不同类型的单元。

对复合结构可能同时⽤到多种类型的单元，此时还需要考虑不同类型单元的连接处理等问题。

3、根据计算分析的精度、周期及费⽤等⽅⾯的要求，合理确定单元的尺⼨和阶次。

4、根据⼯程需要，确定分析类型和计算⼯况。

要考虑参数区间及确定最危险⼯况等问题。

5、根据结构的实际⽀撑情况及受载状态，确定各⼯况的边界约束和有效计算载荷。

⼆、选择位移插值函数1、位移插值函数的要求在有限元法中通常选择多项式函数作为单元位移插值函数，并利⽤节点处的位移连续性条件，将位移插值函数整理成以下形函数矩阵与单元节点位移向量的乘积形式。

位移插值函数需要满⾜相容（协调）条件，采⽤多项式形式的位移插值函数，这⼀条件始终可以满⾜。

但近年来有⼈提出了⼀些新的位移插值函数，如：三⾓函数、样条函数及双曲函数等，此时需要检查是否满⾜相容条件。

2、位移插值函数的收敛性（完备性）要求：1）位移插值函数必须包含常应变状态。

2）位移插值函数必须包含刚体位移。

3、复杂单元形函数的构造对于⾼阶复杂单元，利⽤节点处的位移连续性条件求解形函数，实际上是不可⾏的。

ANSYS单元结构解释ANSYS单元结构解释结构线单元2-D LINK1 3-D LINK8 LINK10 LINK11 LINK180结构梁单元2-D BEAM3 BEAM23 BEAM543-D BEAM4 BEAM24 BEAM44 BEAM188 BEAM189结构线单元2-DLINK1单元可用于不同的工程应用中，依具体的应用，该单元可模拟桁架、链杆及弹簧等。

该二维杆单元每个节点的自由度只考虑x,y 两个方向的线位移，是一种可承受单轴拉压的单元。

因为只用于铰接结构，故本单元不能承受弯矩作用。

结构线单元3-DLINK8单元是这种单元LINK1的三维情况。

LINK10 元素描述 LINK10 是一种 3 D 轴向拉或压的杆单元，有双线性刚度阵。

在仅处于拉状态时，如果进入压状态（在模拟松弛的电缆或松弛的链情况) 时，刚化被取消。

它也能被用于动态分析 ( 由于惯性或阻尼作用)，此时希望使用松弛单元的特性但是单元的运动不是主要的。

This element is a line version of SHELL41 with KEYOPT(1) = 2, the “cloth” option.如果分析的目的要研究元素(由于没有松弛元素) 的运动, 应该用一种相似不会松弛的单元, 像 LINK8 或 PIPE59。

LINK10 也不能在最后结果是拉状态并且该处是集中力的情况时使用，但是迭代收敛结果可以使用在松弛条件。

（LINK10 should also not be used for static convergence applications where the final solution is known to be a taut structure but a slack condition is possible while iterating to a final converged solution.）如果要使用LINK10，那么因为这种原因，应当使用一种不同的单元或缓慢运动技术。

ANSYS结构分析教程篇ANSYS结构分析基础篇⼀、总体介绍进⾏有限元分析的基本流程：1.分析前的思考1)采⽤哪种分析(静态，模态，动态...）2)模型是零件还是装配件(零件可以form a part形成装配件，有时为了划分六⾯体⽹格采⽤零件，但零件间需定义bond接触)3)单元类型选择(线单元，⾯单元还是实体单元)4)是否可以简化模型(如镜像对称，轴对称)2.预处理1)建⽴模型2)定义材料3)划分⽹格4)施加载荷及边界条件3.求解4.后处理1)查看结果(位移，应⼒，应变，⽀反⼒)2)根据标准规范评估结构的可靠性3)优化结构设计⾼阶篇：⼀、结构的离散化将结构或弹性体⼈为地划分成由有限个单元，并通过有限个节点相互连接的离散系统。

这⼀步要解决以下⼏个⽅⾯的问题:1、选择⼀个适当的参考系，既要考虑到⼯程设计习惯，⼜要照顾到建⽴模型的⽅便。

2、根据结构的特点，选择不同类型的单元。

对复合结构可能同时⽤到多种类型的单元，此时还需要考虑不同类型单元的连接处理等问题。

3、根据计算分析的精度、周期及费⽤等⽅⾯的要求，合理确定单元的尺⼨和阶次。

4、根据⼯程需要，确定分析类型和计算⼯况。

要考虑参数区间及确定最危险⼯况等问题。

5、根据结构的实际⽀撑情况及受载状态，确定各⼯况的边界约束和有效计算载荷。

⼆、选择位移插值函数?1、位移插值函数的要求在有限元法中通常选择多项式函数作为单元位移插值函数，并利⽤节点处的位移连续性条件，将位移插值函数整理成以下形函数矩阵与单元节点位移向量的乘积形式。

位移插值函数需要满⾜相容（协调）条件，采⽤多项式形式的位移插值函数，这⼀条件始终可以满⾜。

但近年来有⼈提出了⼀些新的位移插值函数，如：三⾓函数、样条函数及双曲函数等，此时需要检查是否满⾜相容条件。

2、位移插值函数的收敛性（完备性）要求：?1）位移插值函数必须包含常应变状态。

2）位移插值函数必须包含刚体位移。

3、复杂单元形函数的构造对于⾼阶复杂单元，利⽤节点处的位移连续性条件求解形函数，实际上是不可⾏的。

一、求解分析（结构分析）（一）求解设置（二）边界条件l对称与反对称边界条件——实体和单元1）针对对称边界条件下实体结构的分析，可利用ANSYS对称边界条件设置，求解半个或者1/4实体结构，将所得结果对称/循环，得到整体结果分析；2）针对反对称边界条件下实体结构的分析，可利用ANSYS反对称边界条件设置，求解半个实体结构，将所得结果按180度CYCLIC循环对称定义，注意反对称要求如下因素亦满足反对称条件：材料、约束方程、载荷、外形。

l位移边界条件——实体和单元1. 位移约束与强制位移位移约束（displacement constraint）是在节点、或关键点（自由点）上施加某种条件以限制其沿某一自由度方向的运动强制位移（enforced displacement）是在约束点（节点或关键点）上施加某种条件以促使其沿某一自由度方向运动。

2. 限制刚体位移问题一：分析中有时会遇到这样一种情况：即外加载荷是整体平衡的，从理论上来说不会引起刚体位移，只会引起结构变形。

但在进行静力分析时，如果不施加任何约束却会由于刚度矩阵的奇异无法计算，这是怎么回事？这种情况下约束应该如何施加？答1：这种情况叫做Pure Neumann boundary value problem。

这种情况下所得到的位移都是相对位移加上一个常数，常数即为刚体位移。

一个很简单的例子就是一根一维杆两端加大小相等方向相反的力，杆内任意两点之间有相对位移，但每一点的绝对位移却是整个杆的刚体位移加上相对位移。

但是固定杆上的一个点，就会使这个常数即刚体位移为零。

对于Pure Neumann boundary value problem，讨论位移或者温度没有意义，有意义的量是位移和温度的导数的函数。

梁，杆，壳单元可以通过固定任意一个节点，如固定刚体，刚体转动。

对于体单元或者二维平面单元，固定一个点，会导致应力奇异。

应该固定一个面或一条线，这样就不会发生应力奇异了。

1 ANSYS概述1.1 ANSYS简介ANSYS是一种广泛的商业套装工程分析软件。

所谓工程分析软件，主要是在机械结构系统受到外力负载所出现的反应，例如应力、位移、温度等，根据该反应可知道机械结构系统受到外力负载后的状态，进而判断是否符合设计要求。

一般机械结构系统的几何结构相当复杂，受的负载也相当多，理论分析往往无法进行。

想要解答，必须先简化结构，采用数值模拟方法分析。

由于计算机行业的发展，相应的软件也应运而生，ANSYS软件在工程上应用相当广泛，在机械、电机、土木、电子及航空等领域的使用，都能达到某种程度的可信度，颇获各界好评。

使用该软件，能够降低设计成本，缩短设计时间。

到80年代初期，国际上较大型的面向工程的有限元通用软件主要有：ANSYS, NASTRAN, ASKA, ADINA, SAP等。

以ANSYS为代表的工程数值模拟软件，是一个多用途的有限元法分析软件，它从1971年的2.0版本与今天的5.7版本已有很大的不同，起初它仅提供结构线性分析和热分析，现在可用来求结构、流体、电力、电磁场及碰撞等问题的解答。

它包含了前置处理、解题程序以及后置处理，将有限元分析、计算机图形学和优化技术相结合，已成为现代工程学问题必不可少的有力工具。

1.2ANSYS软件主要功能ANSYS软件是融结构、热、流体、电磁、声学于一体的大型通用有限元软件，可广泛的用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、生物医学、水利、日用家电等一般工业及科学研究。

该软件提供了不断改进的功能清单，具体包括：结构高度非线性分析、电磁分析、计算流体力学分析、设计优化、接触分析、自适应网格划分及利用ANSYS参数设计语言扩展宏命令功能。

1.3ANSYS软件主要特点主要技术特点：∙唯一能实现多场及多场耦合分析的软件∙唯一实现前后处理、求解及多场分析统一数据库的一体化大型FEA软件∙唯一具有多物理场优化功能的FEA软件∙唯一具有中文界面的大型通用有限元软件∙强大的非线性分析功能∙多种求解器分别适用于不同的问题及不同的硬件配置∙支持异种、异构平台的网络浮动，在异种、异构平台上用户界面统一、数据文件全部兼容∙强大的并行计算功能支持分布式并行及共享内存式并行∙多种自动网格划分技术∙良好的用户开发环境支持的图形传递标准：∙SAT∙Parasolid∙STEP与CAD软件的接口∙Unigraphics∙Pro/ENGINEER∙I-Deas∙Catia∙CADDS∙SolidEdge∙SolidWorks1.4 运行环境（ANSYS5.7）Computer: Pentium-class systemMemory (RAM): 64 MB以上Hard Disk: 500MB以上自由空间.Operating System: Microsoft Windows 2000, Windows NT 4.0 (SP 5 or higher) or Windows 98 Graphics: A Windows 2000/NT 4.0 or 98 supported Graphics Card, capable of 1024x768 in High Color (16-bit).A 17 inch monitor (or larger) compatible with the above mentioned card is recommended.2ANSYS 的基本使用2.1 ANSYS环境简介ANSYS有两种模式：一种是交互模式（Interactive Mode），另一个是非交互模式（Batch Mode）。

ansys分析混凝土的若干问题1. 讨论两种Ansys求极限荷载的方法（1）力加载可以通过对应的方法（比如说特征值屈曲）估计结构的极限荷载的大致范围，然后给结构施加一个稍大的荷载，打开自动荷载步二分法进行非线性静力分析，最后计算会因不收敛终止，则倒数第二个子步对应的就是结构的极限荷载；另外，也可以选择弧长法，采用足够的子步（弧长法可以一直分析到极限承载力之后的过程）同样可以从绘制的荷载位移曲线或计算结果中找出结构的极限荷载。

（2）位移加载给结构施加一个比较大的位移，打开自动荷载步二分法进行非线性分析，保证足够的子步数，这样也可以分析到极限荷载以后，通过绘制荷载位移曲线或查看相应结果文件也可知道结构的极限荷载。

希望众高手讨论一下（1）弧长法求极限荷载的收敛性问题，如何画到荷载位移曲线的下降段？（2）位移法求极限荷载的具体步骤？2. 需要注意的问题1. 由于SOLID 65单元本身是基于弥散裂缝模型和最大拉应力开裂判据，因此在很多情况下会因为应力集中而使混凝土提前破坏，从而和试验结果不相吻合，因此，在实际应用过程中应该对单元分划进行有效控制，根据作者经验，当最小单元尺寸大于5cm 时，就可以有效避免应力集中带来的问题；2. 支座是另一个需要注意的问题。

在有限元分析中，很多时候约束是直接加在混凝土节点上，这样很可能在支座位置产生很大的应力集中，从而使支座附近的混凝土突然破坏，造成求解失败。

因此，在实际应用过程中，应该适当加大支座附近单元的尺寸或者在支座上加一些弹性垫块，避免支座的应力集中；3. 六面体的SOLID 65 单元一般比四面体的单元计算要稳定且收敛性好，因此，只要条件允许，应该尽量使用六面体单元；4. 正确选择收敛标准，一般位移控制加载最好用位移的无穷范数控制收敛，而用力控制加载时可以用残余力的二范数控制收敛。

在裂缝刚刚出现和接近破坏的阶段，可以适当放松收敛标准，保证计算的连续性；3. 关于下降段的问题1）在实际混凝土中都有下降段，但是在计算的时候要特别小心下降段的问题。