ANSYS耦合与约束方程
ansys中如何进行刚性连接
ansys中如何进行刚性连接一般来说，按“杆梁壳体”单元顺序，只要后一种单元的自由度完全包含前一种单元的自由度，则只要有公共节点即可，不需要约束方程，否则需要耦合自由度与约事方程。
例如：（1）杆与梁、壳、体单元有公共节点即可，不需要约束方程。
（2）梁与壳有公共节点即可，也不需要约束写约束方程；壳梁自由度数目相同，自由度也相同，尽管壳的rotz是虚的自由度，也不妨碍二者之间的关系，这有点类同于梁与杆的关系。
（3）梁与体则要在相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
（4）壳与体则也要相同位置建立不同的节点，然后在节点处耦合自由度与施加约束方程。
举例：有一长为100mm的矩形截面梁，截面为10X1mm，与一规格为20mmX7mmX10mm的实体连接，约束实体的端面，在梁端施加大小为3N的y方向的压力，梁与实体都为一材料，弹性模量为30Gpa，泊松比为0.3。
本例主要讲解梁与实体连接处如何利用耦合及约束方程进行处理。
命令流如下：FINI/CLE/FILNAME,BEAM_AND_SOLID_ELEMENTS_CONNECTION !定义工作文件名/TITLE,COUPLE_AND_CONSTRAINT_EQUATION !定义工作名/PREP7ET,1,SOLID95 !定义实体单元类型为SOLID95ET,2,BEAM4 !定义梁单元类型为BEAM4MP,EX,1,3E4 !定义材料的弹性模量MP,PRXY,1,0.3 !定义泊松比R,1 !定义实体单元实常数R,2,10.0,10/12.0,1000/12.0,10.0,1.0 !定义梁单元实常数BLC4,,,20,7,10 !创建矩形块为实体模型WPOFFS,0,3.5 !将工作平面向Y方向移动3.5WPROTA,0,90 !将工作平面绕X轴旋转90度VSBW,ALL !将实体沿工作平面剖开WPOFFS,0,5 !将工作平面向Y方向移动5WPROTA,0,90 !将工作平面绕X轴旋转90度VSBW,ALL !将实体沿工作平面剖开WPCSYS,-1 !将工作平面设为与总体笛卡儿坐标一致K,100,20,3.5,5 !创建关键点K,101,120,3.5,5 !创建关键点L,100,101 !连接关键点生成梁的线实体LSEL,S,LOC,X,21,130 !选择梁线LATT,1,2,2 !指定梁的单元属性LESIZE,ALL,,,10 !指定梁上的单元份数LMESH,ALL !划分梁单元VSEL,ALL !选择所有实体VATT,1,1,1 !设置实体的单元属性ESIZE,1 !指定实体单元尺寸MSHAPE,0,2D !设置实体单元为2DMSHKEY,1 !设置为映射网格划分方法VMESH,ALL !划分实体单元ALLS !全选FINI !退出前处理!------------------------/SOLU !进入求解器ASEL,S,LOC,X,0 !选择实体的端面DA,ALL,ALL !约束实体端面ALLS !全选FK,101,FY,-3.0 !在两端施加Y向压力CP,1,UX,1,21 !耦合节点1和节点21X方向自由度CP,2,UY,1,21 !耦合节点1和节点21Y方向自由度CP,3,UZ,1,21 !耦合节点1和节点21Z方向自由度CE,1,0,626,UX,1,2328,UX,-1,1,ROTY,-ABS(NZ(626)-NZ(2328)) !设置约束方程CE,2,0,67,UX,1,4283,UX,-1,1,ROTZ,-ABS(NY(67)-NY(4283)) !设置约束方程CE,3,0,67,UZ,1,4283,UZ,-1,1,ROTX,-ABS(NY(67)-NY(4283)) !设置约束方程ALLS !全选SOLVE !保存FINI !退出求解器!------------------------/POST1 !进入通用后处理PLNSOL, U,Y, 0,1.0 !显示Y方向位移PLNSOL, S,EQV, 0,1.0 !显示等效应力ETABLE,ZL1,SMISC,1 !读取梁单元上I节点X方向的力ETABLE,ZL2,SMISC,7 !读取梁单元上J节点X方向的力ETABLE,MZ1,SMISC,6 !读取梁单元上I节点Z方向的力矩ETABLE,MZ2,SMISC,12 !读取梁单元上J节点Z方向的力矩PLETAB,ZL1 !显示梁单元X方向的力PLETAB,MZ1 !显示梁单元Z方向力矩!**********************************************上面所述的不同单元之间的接连方法主要是用耦合自由度和约束方程来实现的，有一定的局限性，只适用于小位移，下面介绍一种支持大位移算法的方法，MPC 法。
ANSYS_高级接触问题
图3-3
• 例：Seal.dat (图3-3) • Step 5. 建立接触面单元 • 设置接触单元属性、选择可变形体表面节点，并在可 变形体上建立接触单元（过程与在可变形体上建立目标单 元相同） Main Menu > Preprocessor > Modeling > Create > Elements > Surf/Contact > Surf to Surf (ESURF) • 这些接触单元与基体有同样的阶数（低阶或高阶）。 • 注意，在壳或梁单元上建立目标单元或接触单元时， 可以选择要在梁或壳单元的顶层还是底层建立单元。
· 接触分析中自动时间步的其它注意事项：
• 与所有其它非线性分析一样，对接触问题，时间步长是非 常有力的提高收敛性的工具。 • 采用足够小的时间步长以获得收敛。 • 对于冲击瞬态分析，必须使用足够数量的计算步以描 述表面间的动量转移。 • 对于路径相关现象（如接触摩擦），相对较小的最大 时间步长对计算精度是必须的。
ANSYS 高级接触问题
• • 接触问题概述 在工程中会遇到大量的接触问题，如齿轮的啮合、法兰 联接、机电轴承接触、卡头与卡座、密封、板成形、冲 击等等。接触是典型的状态非线性问题，它是一种高度 非线性行为。接触例子如图1：
ANSYS 高级接触问题
• 分析中常常需要确定两个或多个相互接触物体的 位移、接触区域的大小和接触面上的应力分布。 • 接触分析存在两大难点： • 在求解之前，你不知道接触区域的范围；表面之 间是接触还是分开是未知的；表面之间突然接触 或突然不接触会导致系统刚度的突然变化。 • 大多数接触问题需要计算摩擦。摩擦是与路径有 关的现象，摩擦响应还可能是杂乱的，使问题求 解难以收敛。
§3 自动时间步、控制
ANSYS中杆单元和壳单元的单元耦合问题
ANSYS中杆单元和壳单元的单元耦合问题在比较复杂的结构的有限元分析中，不同的结构部件通常使用不同类型的单元来模拟。
通常情况下，不同类型的单元的各个节点的自由度数目是不同的，不同类型单元的连接节点处的自由度的耦合问题，是一个比较令人头疼的问题。
在ANSYS中通常可以用耦合命令CP来耦合不同类型单元在连接节点处的自由度(DOF)。
也可以用CE命令来认为添加自由度之间的约束方程来达到耦合的目的。
下面是一个简单的算例，使用了CE命令来耦合连接节点处的自由度。
模型是航天器的机翼的一个Section的某一个隔框。
上下表皮是薄壳结构，用Shell63单元来模拟，在上下表皮之间有起支撑作用的杆件，用link8单元来模拟。
建模的时候，link8单元和shell63单元在连接有各自独立的节点。
即：link8单元和shell63单元的节点在连接处是重合的，但是，节点编号是各自独立的。
link8单元在每个节点有ux,uy,uz3个平动自由度；shell63在每个节点有ux,uy,uz这3个平动自由度和rotx,roty,rotz这3个转个自由，共6个自由度。
在耦合节点处，两个耦合节点的ux,uy,uz自由度应该是相等的。
这个等式可以用CE命令来描述。
完整的命令流如下：finish/clear,start/prep7!定义第一种材料属性；mp,ex,1,30e6mp,prxy,1,0.3!定义shell63单元和实常数；et,1,shell63r,1,1e-3!建立几何模型；rectng,31.8,33.2,0,0.3556agen,2,1,1,1,0,0,1a,1,4,8,5a,6,7,3,2KL,7,0.5, ,KL,3,0.5, ,在关键点处生成节点；nkpt,100,4 !与编号为117的节点耦合nkpt,101,9 !与编号为169的节点耦合nkpt,102,10 !与编号为120的节点耦合nkpt,103,7 !与编号为160的节点耦合mat,1type,1real,1lesize,1,,,6lesize,3,,,6lesize,5,,,6lesize,7,,,6lesize,9,,,6lesize,10,,,6lesize,11,,,6lesize,12,,,6lesize,2,,,6lesize,4,,,6lesize,6,,,6lesize,8,,,6MSHAPE, 0, 2DMSHKEY, 1allselamesh,all!定义第二种材料属性；mp,ex,2,30e4mp,prxy,2,0.3!定义link8单元和实常数；et,2,link8r,2,28.26e-6mat,2type,2real,2e,101,102e,100,101e,102,103!CE, NEQN, CONST, NODE1, Lab1, C1, NODE2, Lab2, C2, NODE3, !Lab3, C3ce,1,0,100,ux,1,117,ux,-1 !节点100的ux＝节点117的ux;ce,2,0,100,uy,1,117,uy,-1 !节点100的uy＝节点117的uy;ce,3,0,100,uz,1,117,uz,-1 !节点100的uz＝节点117的uz;ce,4,0,101,ux,1,169,ux,-1 !同上；ce,5,0,101,uy,1,169,uy,-1ce,6,0,101,uz,1,169,uz,-1ce,7,0,102,ux,1,120,ux,-1ce,8,0,102,uy,1,120,uy,-1ce,9,0,102,uz,1,120,uz,-1ce,10,0,103,ux,1,160,ux,-1ce,11,0,103,uy,1,160,uy,-1ce,12,0,103,uz,1,160,uz,-1!施加约束，底面约束所有的自由度；DA,1,all!施加分布载荷；SFA,2,1,PRES,-1e4finish/solusolve !求解；finish/post1PLNSOL, S,X, 0,1.0finishPLNSOL, S,X, 0,1.0 对应的结果云图PLNSOL, S,EQV, 0,1.0对应的结果云图。
桥梁计算(常用的计算方法)
＊＊＊桥梁仿真单元类型 (1)一、建议选用的单元类型 (1)二、常见桥梁连接部位 (2)三、桥梁基础的处理方式 (2)＊＊＊桥梁常见模型处理 (2)一、桥梁中常用的模型可以用相应的单元 (2)二、桥梁建模要综合运用各种合适的单元 (3)三、选用合适的分析方法 (3)施加预应力的方式 (3)一、预应力的模拟方式 (3)二、建立预应力的模型 (5)＊＊＊土弹簧的模拟 (5)＊＊＊混凝土的模拟 (5)工况组合 (6)一、典型的荷载工况步骤 (6)二、存储组合后的荷载工况 (6)风荷载的确定 (7)地震波的输入 (7)初应力荷载 (8)Ansys可采用两种方法来实现铰接： (8)AUTOCAD模型输入 (9)用ANSYS作桥梁计算十三（其他文件网格划分） (12)(一)时间选项 (13)(二)子步数和时间步大小 (13)(三)自动时间步长 (14)(四)阶跃或递增载荷 (14)关于阶跃载荷和逐渐递增载荷的说明： (14)一、用于动态和瞬态分析的命令 (14)二、非线性选项 (14)三、输出控制 (15)重新启动一个分析 (16)一、重启动条件 (16)二、一般重启动的步骤 (17)三、边界条件重建 (17)在Ansys单元库中，有近200种单元类型，在本章中将讨论一些在桥梁工程中常用到的单元，包括一些单元的输人参数，如单元名称、节点、自由度、实常数、材料特性、表面荷载、体荷载、专用特性、关键选项KEYOPl等。
＊＊＊关于单元选择问题这是一个大问题，方方面面很多，主要是掌握有限元的理论知识。
首先当然是由问题类型选择不同单元，二维还是三维，梁，板壳，体，细梁，粗梁，薄壳，厚壳，膜等等，再定义你的材料：各向同性或各向异性，混凝土的各项’参数，粘弹性等等。
接下来是单元的划分与网格、精度与求解时间的要求等选择，要对各种单元的专有特性有个大概了解。
使用Ansys，还要了解Ansys的一个特点是笼统与通用，因此很多东西被掩盖到背后去了。
耦合和约束方程82页PPT
36、如果我们国家的法律中只有某种 神灵， 而不是 殚精竭 虑将神 灵揉进 宪法， 总体上 来说， 法律就 会更好 。—— 马克·吐 温 37、纲纪废弃之日，便是暴政兴起之 时。— —威·皮 物特
38、若是没有公众舆论的支持，法律 是丝毫 没有力 量的。 ——菲 力普斯 39、一个判例造出另一个判例，它们 迅速累 聚，进 而变成 法律。 ——朱 尼厄斯
40、人类法律，事物有规律，这是不 ——德 谟克利 特 67、今天应做的事没有做，明天再早也 是耽误 了。——裴斯 泰洛齐 68、决定一个人的一生，以及整个命运 的，只 是一瞬 之间。 ——歌 德 69、懒人无法享受休息之乐。——拉布 克 70、浪费时间是一桩大罪过。——卢梭
ANSYS多物理耦合场有限元分析详细步骤操作
ANSYS热分析
传导
• 传导引起的热通量流由传导的傅立叶定律决定：
q*
Knn
T n
heat flow rate per unit area in
direction n
Where,
Knn thermalconductivity in direction n
T temperature
• 对流一般作为面边界条件施加
TB
Ts
ANSYS热分析
热力学第一定律
• 能量守恒要求系统的能量改变与系统边 界处传递的热和功数值相等。
• 能量守恒在一个微小的时间增量下可以 表示为方程形式
Estored Ein thru the boundary Eout thru the boundary Egenerated 0
SHELL57 SHELL131
LINK31,32,33,34
Quadratic
PLANE77 PLANE35
SOLID90 SOLID87
SHELL132
• 材料特性
– 至少需要 Kxx — 稳态分析热传导系数。 – 如果是瞬态分析，则需要比热 (C) 。 – 优先设置为 “thermal” （热分析），在 GUI 方式中只显示热材料特性。
T
Dt
如果时间步长 太大, 就不能 得到足够的温度梯度。
t
一种方法是先指定一个相对较保守的初始时间步长，然后使用自动时间步长 按需要增加时间步。下面说明使用自动时间步长大致估计初始时间步长的方 法。
ANSYS热分析
时间步大小说明 (续)
在瞬态热分析中大致估计初始时间步长，可以使用Biot和Fourier数。 Biot 数 是无量纲的对流和传导热阻的比率:
ansys单元划分技巧[指南]
![ansys单元划分技巧[指南]](https://img.taocdn.com/s3/m/0eed4cc0b8f3f90f76c66137ee06eff9aff84952.png)
ansys单元划分技巧众所周知，对于有限元分析来说，网格划分是其中最关键的一个步骤，网格划分的好坏直接影响到解算的精度和速度。
在ANSYS中，大家知道，网格划分有三个步骤：定义单元属性（包括实常数）、在几何模型上定义网格属性、划分网格。
在这里，我们仅对网格划分这个步骤所涉及到的一些问题，尤其是与复杂模型相关的一些问题作简要阐述。
一、自由网格划分自由网格划分是自动化程度最高的网格划分技术之一，它在面上（平面、曲面）可以自动生成三角形或四边形网格，在体上自动生成四面体网格。
通常情况下，可利用ANSYS的智能尺寸控制技术（SMARTSIZE 命令）来自动控制网格的大小和疏密分布，也可进行人工设置网格的大小（AESIZE、LESIZE、KESIZE、ESIZE等系列命令）并控制疏密分布以及选择分网算法等（MOPT命令）。
对于复杂几何模型而言，这种分网方法省时省力，但缺点是单元数量通常会很大，计算效率降低。
同时，由于这种方法对于三维复杂模型只能生成四面体单元，为了获得较好的计算精度，建议采用二次四面体单元（92号单元）。
如果选用的是六面体单元，则此方法自动将六面体单元退化为阶次一致的四面体单元，因此，最好不要选用线性的六面体单元（没有中间节点，比如45号单元），因为该单元退化后为线性的四面体单元，具有过刚的刚度，计算精度较差；如果选用二次的六面体单元（比如95号单元），由于其是退化形式，节点数与其六面体原型单元一致，只是有多个节点在同一位置而已，因此，可以利用TCHG命令将模型中的退化形式的四面体单元变化为非退化的四面体单元，减少每个单元的节点数量，提高求解效率。
在有些情况下，必须要用六面体单元的退化形式来进行自由网格划分，比如，在进行混合网格划分（后面详述）时，只有用六面体单元才能形成金字塔过渡单元。
对于计算流体力学和考虑集肤效应的电磁场分析而言，自由网格划分中的层网格功能（由LESIZE命令的LAYER1和LAYER2域控制）是非常有用的。
ansys 热 结构耦合分析
第21章热－结构耦合分析热－结构耦合问题是结构分析中通常遇到的一类耦合分析问题。
由于结构温度场的分布不均会引起结构的热应力，或者结构部件在高温环境中工作，材料受到温度的影响会发生性能的改变，这些都是进行结构分析时需要考虑的因素。
为此需要先进行相应的热分析，然后在进行结构分析。
热分析用于计算一个系统或部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度（热通量）等。
本章主要介绍在ANSYS中进行稳态、瞬态热分析的基本过程，并讲解如何完整的进行热－结构耦合分析。
21.1 热－结构耦合分析简介热－结构耦合分析是指求解温度场对结构中应力、应变和位移等物理量影响的分析类型。
对于热－结构耦合分析，在ANSYS中通常采用顺序耦合分析方法，即先进行热分析求得结构的温度场，然后再进行结构分析。
且将前面得到的温度场作为体载荷加到结构中，求解结构的应力分布。
为此，首先需要了解热分析的基本知识，然后再学习耦合分析方法。
21.1.1 热分析基本知识ANSYS热分析基于能量守恒原理的热平衡方程，用有限元法计算各节点的温度，并导出其它热物理参数。
ANSYS热分析包括热传导、热对流及热辐射三种热传递方式。
此外，还可以分析相变、有内热源、接触热阻等问题。
热传导可以定义为完全接触的两个物体之间或一个物体的不同部分之间由于温度梯度而引起的内能的交换。
热对流是指固体的表面和与它周围接触的流体之间，由于温差的存在引起的热量的交换。
热辐射指物体发射电磁能，并被其它物体吸收转变为热的热量交换过程。
如果系统的净热流率为0，即流入系统的热量加上系统自身产生的热量等于流出系统的热量：q流入+q生成-q流出=0，则系统处于热稳态。
在稳态热分析中任一节点的温度不随时间变化。
瞬态传热过程是指一个系统的加热或冷却过程。
在这个过程中系统的温度、热流率、热边界条件以及系统内能随时间都有明显变化。
ANSYS热分析的边界条件或初始条件可分为七种：温度、热流率、热流密度、对流、辐射、绝热、生热。