1、CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN初级培训ANSYS/STRUCTURAL基础基础安世亚太西安办事处安世亚太西安办事处2022-7-30第1页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册有限元法基础有限元法基础有限单元法在50年代起源于航空工程飞机结构的矩阵分析，结构矩阵分析认为，一个结构可以看作是由有限个力学小单元相互连接而成的集合体，表征单元力学特性的刚度矩阵可比喻作建筑中的砖瓦。装配在一起就能提供整个
2、结构的力学特性。应用有限元法求解连续体时，把求解区域分为有限个单元，并在每个单元上指定有限个节点。一般可以认为相邻单元在节点连接成一组单元的集合体，用以模拟或逼近求解区域进行分析，同时选定场函数的节点值。例如取节点位移作为基本未知量。假设一个插值函数近似地表示位移分布规律。再利用变分原理或其它方法建立单元节点力和位移之间的力学特性关系。得到一组以节点位移为未知量的代数方程组。从而求解节点位移分量。2000，4第2页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册下面以小变形弹性静力问题为例，加以详细
3、介绍。下面以小变形弹性静力问题为例，加以详细介绍。几何方程：几何方程：eij=1/2(ui,j+uj,i)物理方程：物理方程：s sij=aijklekl平衡方程：平衡方程：s sij,j+fi=0边界条件：边界条件：位移已知边界条件位移已知边界条件 ui=ui （在边界（在边界u上位移已知）上位移已知）外力已知边界条件外力已知边界条件 s sij,j+pi=0（在边界（在边界p上外力已知）上外力已知）简例（续）简例（续）2000，4第3页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册简例（续）简
4、例（续）上标表示单元号码，下标表示节点号码。因铰接不存在力矩，每节点力和位移各有两个分量，共四个自由度，用四个方程来描述力位移关系：F1x1=K11U11+K12V11+K13U21+K14V21 F1y1=K21U11+K22V11+K23U21+K24V21 （1-1）F1x2=K31U11+K32V11+K33U21+K34V21 F1y2=K41U11+K42V11+K43U21+K44V212000，4第4页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册简例（续）简例（续）即 F 1=K
5、 11在（1-1）中令U11=1，V11=U21=V21=0则F1x1=K11，F1y1=K21，F1x2=K31，F1y2=K41 上式表明，当节点1沿X方向产生一单位位移，而单元1的其余节点位移为零时，各节点施于单元1上的力将组成一平衡力系，表示单元1抵抗位移U11的刚度。根据节点间位移协调关系。U11=U22，V11=V22 又根据各节点的平衡条件有 F=K2000，4第5页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册有限元分析步骤有限元分析步骤有限元法可分为几步：结构的离散化 选择位移模
6、式 即假定位移是坐标的某种简单的函数这种函数称为位移模式或插值函数通常选多项式作为位移模式一般来说，多项式的项数应等于单元的自由度数。f=Nf=Ne e (1-9)(1-9)f单元内任一点的位移列阵 e 单元的节点位移列阵 N形函数矩阵 2000，4第6页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册有限元分析步骤（续）有限元分析步骤（续）分析单元的力学特性由（1-9）导出节点位移表示单元应变的关系式 =B=Be e (1-10)(1-10)单元内任一点的应变列阵 B单元应变矩阵由（1-10）导出
7、 =DB=DBe e 单元内任一点的应力矩阵 D与单元材料有关的弹性矩阵利用变分原理，建立作用于单元上的节点力和位移之间的关系式 Fe=Kee 2000，4第7页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册有限元分析步骤（续）有限元分析步骤（续）集合所有单元的平衡方程，集合依据的是所有相邻 单元在公共节点处的位移相等；建立总体的有限元方程组。引入边界条件 求解有限元方程组，得到未知节点位移 计算单元应力，对不同的单元，对应力的处理还有不同的方法2000，4第8页，共68页。CFD Analysi
8、s with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS文件结构文件结构二进制文件二进制文件Jobname.db （数据库文件）Jobname.dbb （备份文件）Jobname.rst （结构分析结果文件）Jobname.rth （热分析结果文件）Jobname.rmg （电磁场分析结果文件）Jobname.rfl （流体分析结果文件）Jobname.tri （三角化刚度矩阵文件）Jobname.emat（单元矩阵文件）Jobname.esav （单元保存文件）2000，4第9页，共68页。CFD Analysis with A
9、NSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS文件结构（续）文件结构（续）文本文件Jobname.log（命令日志文件）Jobname.err（错误及警告信息文件）ANSYS6.0可以改变 Jobname Work directory2000，4第10页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS内存管理内存管理2000，4第11页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis
10、 with ANSYS/FLOTRAN培训手册ANSYS内存管理内存管理(续续)Workspace2000，4第12页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS 分析过程分析过程 建模(/prep7)加载及求解(/solu)观察结果(/post1,/post26)2000，4第13页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册建模建模 确定工作名和分析标题（/TITLE）定义单位 选单
11、元 定义实常数 定义材料性质 创建几何模型 建立有限元模型2000，4第14页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册 ANSYS/Structural求解功能求解功能Static-结构静力问题（包括线性和非线性问题）结构静力问题（包括线性和非线性问题）Modal-模态振动特性计算分析（结构固有频率和振型）模态振动特性计算分析（结构固有频率和振型）Harmonic-谐波分析谐波分析Transient-瞬态分析瞬态分析Spectrum-谱分析谱分析Eigen Buckling-特征值屈曲分析（
12、线性）特征值屈曲分析（线性）Substructural-子结构分析子结构分析。ANSYS/Structural求解功能求解功能2000，4第15页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS非线性非线性材料非线性材料非线性几何非线性几何非线性单元非线性单元非线性2000，4第16页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册几何非线性几何非线性大应变大应变大挠度大挠度应力刚化应力刚化旋转
13、软化旋转软化2000，4第17页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册材料非线性材料非线性速率无关的塑性速率无关的塑性速率相关的塑性速率相关的塑性超弹性超弹性粘弹性粘弹性混凝土混凝土非线弹性非线弹性蠕变蠕变膨胀膨胀2000，4第18页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册材料非线性材料非线性 流动准则流动准则 屈服准则屈服准则 强化准则强化准则2000，4第19页，共68页。CFD A
14、nalysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册单元非线性单元非线性接触接触 点点-点点 点点-线线 点点-面面 面面-面面 刚刚-柔柔 柔柔-柔柔单元生死单元生死(焊接焊接,土壤开挖土壤开挖)非线性单元非线性单元 LINK10,COMBINE40.2000，4第20页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册几何模型与有限元模型几何模型与有限元模型关键点线面体单元节点几何模型有限元模型2000，4第21页，共68页。C
15、FD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册加载加载载荷的含义（六大类）载荷的含义（六大类）l 自由度约束（位移，对称边界条件，温度等）自由度约束（位移，对称边界条件，温度等）l 节点力（力，力矩，热流量）节点力（力，力矩，热流量）l 表面载荷（压力，对流，热流量）表面载荷（压力，对流，热流量）l 体载荷（结构分析中温度，热分析中热生成率，磁体载荷（结构分析中温度，热分析中热生成率，磁 场场分析中电流密度）分析中电流密度）l 惯性载荷（重力加速度，角加速度）惯性载荷（重力加速度，角加速度）l 耦合场载荷（上述
16、载荷的一种特殊情况，一种分析的结耦合场载荷（上述载荷的一种特殊情况，一种分析的结果作为载荷加于另一种分析中）果作为载荷加于另一种分析中）2000，4第22页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册实体几何模型载荷实体几何模型载荷优点优点改变网格不影响载荷改变网格不影响载荷涉及到的加载实体少涉及到的加载实体少缺点缺点生成的单元在当前激活的单元座标下，节点生成的单元在当前激活的单元座标下，节点为总体直角座标，因此实体与有限元模型可为总体直角座标，因此实体与有限元模型可能有不同座标系统和载荷方向能
17、有不同座标系统和载荷方向实体载荷在凝聚分析中不方便，因载荷加实体载荷在凝聚分析中不方便，因载荷加在主自由度上施加关键点约束较繁锁在主自由度上施加关键点约束较繁锁 不能显示所有实体载荷不能显示所有实体载荷2000，4第23页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册有限元载荷有限元载荷优点优点凝聚分析加载方便凝聚分析加载方便不用考虑约束扩展不用考虑约束扩展缺点缺点网格改变后需重新加载网格改变后需重新加载图形拾取加载不方便图形拾取加载不方便2000，4第24页，共68页。CFD Analysis
18、with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS坐标系种类坐标系种类总体坐标系总体坐标系局部坐标系局部坐标系节点坐标系节点坐标系单元坐标系单元坐标系结果坐标系结果坐标系显示坐标系显示坐标系2000，4第25页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS坐标系类型坐标系类型直角坐标系直角坐标系圆柱坐标系圆柱坐标系球坐标系球坐标系环形坐标系环形坐标系2000，4第26页，共68页。CFD Analysis with ANS
19、YS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册误差估计误差估计SERR SERR 单元能量范数误差单元能量范数误差 e ei i e ei i=1/2=1/2volvolDDT T-1-1 d(vol)d(vol)PRERRPRERR能量范数百分比误差能量范数百分比误差 E E E=100e/E=100e/（U+eU+e）1/21/2 U U整个模型应变能整个模型应变能SDSG SDSG i i连接单元的节点中，最大应力差值的绝对值连接单元的节点中，最大应力差值的绝对值 i in n=a an n-i in n 2000，4第27页，共68页。CFD
20、 Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册选择单元策略选择单元策略 对线性结构（应力）分析，建议采用高阶单元对线性结构（应力）分析，建议采用高阶单元 对非线性应力分析，用低阶单元采用较密网格，对非线性应力分析，用低阶单元采用较密网格，而不用较粗网格高阶单元而不用较粗网格高阶单元 对场分析，低阶单元通常与高阶单元一样可以取得好的对场分析，低阶单元通常与高阶单元一样可以取得好的结果结果 对应力分析，四边形比三角形结果要好对应力分析，四边形比三角形结果要好2000，4第28页，共68页。CFD Analysis wi
21、th ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS单元介绍单元介绍ANSYS单元分类单元分类*ANSYS常用单元常用单元 介绍介绍 ANSYS单元介绍单元介绍2000，4第29页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS单元分类单元分类*ANSYS常用单元常用单元 介绍介绍*ANSYS常用单元常用单元 推荐推荐 ANSYS单元介绍单元介绍第30页，共68页。CFD Analysis with ANSYS/FLOTRANCF
22、D Analysis with ANSYS/FLOTRAN培训手册ANSYSANSYS的单元按形状可分为以下三类：的单元按形状可分为以下三类：质点单元（如质量单元）质点单元（如质量单元）线单元（如杆、梁、管单元）线单元（如杆、梁、管单元）面单元（二维平面面单元（二维平面 、壳等）、壳等）体单元（如体单元（如SOLIDSOLID等）等）其它辅助单元其它辅助单元下面就这些类型中常用的单元作简要介绍。下面就这些类型中常用的单元作简要介绍。ANSYS单元分类单元分类2000，4第31页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSY
23、S/FLOTRAN培训手册 ANSYS单元分类单元分类1.杆单元，包括二维杆单元和三维杆单元，线性调节元杆单元，包括二维杆单元和三维杆单元，线性调节元，主要包括：，主要包括：LINK1，LINK8，LINK10，LINK11，LINK180等。等。2.弹簧阻尼单元，包括弹簧阻尼单元，包括COMBIN系列：系列：COMBIN7,COMBIN14,COMBIN37,COMBIN40等。等。3.质量元，质量元，MASS21。第32页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN初级培训 ANSYS单元分类单元分
24、类4.梁单元，分为二维和三维单元，弹性和塑性；可定梁单元，分为二维和三维单元，弹性和塑性；可定 义义各种截面。主要包括：各种截面。主要包括：BEAM3，BEAM4，BEAM23，BEAM24，BEAM44，BEAM188，BEAM189。5.管单元，分为弹性和塑性管元，直管、弯管、管单元，分为弹性和塑性管元，直管、弯管、T 形管。主形管。主要包括：要包括：PIPE系列。系列。6.二维二维/三维预拉单元三维预拉单元PRET179，可以处理预紧螺栓。可以处理预紧螺栓。7.结构实体单元，分为二维和三维，弹性和塑性，结构实体单元，分为二维和三维，弹性和塑性，低阶低阶 单元和高阶单元。单元和高阶单元。主
25、要包括：主要包括：PLANE系列，系列，SOLID系列。系列。第33页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN初级培训 ANSYS单元分类单元分类8.壳单元，分为薄壳、中厚壳，弹性壳元，塑性壳单壳单元，分为薄壳、中厚壳，弹性壳元，塑性壳单 元，元，膜单元。主要包括：膜单元。主要包括：SHELL系列。系列。9.流体单元，流体单元，FLUID系列。系列。10.超弹单元，超弹单元，HYPER系列。系列。11.粘弹性单元，粘弹性单元，VISCO系列。系列。12.二维、三维表面效应单元。二维、三维表面效应单元
26、。热表面效应单元：热表面效应单元：SURF151，SURF152；结构表面效应单元：结构表面效应单元：SURF153，SURF154。第34页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN初级培训 ANSYS单元分类单元分类13.接触单元，接触单元，分为二维和三维点分为二维和三维点-点接触、点点接触、点-面接触、面接触、面面-面接触。接触单元和目标单元在某些情况下成对使用。面接触。接触单元和目标单元在某些情况下成对使用。主要包括：主要包括：接触单元：接触单元：CONTA系列；目标单元：系列；目标单元：TR
27、AGE系列。系列。14.网格划分辅助单元，网格划分辅助单元，MESH200。与求解无关，与求解无关，不不 影响计算结果。在用低级单元创建高级单元等影响计算结果。在用低级单元创建高级单元等 状态时使用，例如拖拉。状态时使用，例如拖拉。15.其它单元其它单元 超单元：超单元：MATRIX50，主要用于子结构分析等方面；主要用于子结构分析等方面；刚度、阻尼、质量单元：刚度、阻尼、质量单元：MATRIX27，可用于弹性运动学响应分析等方面。可用于弹性运动学响应分析等方面。第35页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/F
28、LOTRAN初级培训ANSYS单元介绍单元介绍ANSYS单元分类常用单元 介绍*ANSYS常用单元推荐第36页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS常用单元介绍常用单元介绍Link8 3维杆单元，可用来模拟桁架、缆索、连杆、弹簧等。维杆单元，可用来模拟桁架、缆索、连杆、弹簧等。沿杠的轴向承受拉压，不承受弯矩。每个节点有沿杠的轴向承受拉压，不承受弯矩。每个节点有3个自由度（沿节个自由度（沿节点坐标系点坐标系x、y、z的平动）。具有塑性、蠕变、膨胀、应力刚的平动）。具有塑性、蠕变、膨
29、胀、应力刚化、大变形、大应变等功能。化、大变形、大应变等功能。单元无选项要求。单元无选项要求。所需实常数：截面积和初始应变。所需实常数：截面积和初始应变。同类型单元有同类型单元有Link1(2D)、Link10、Link11等。等。2000，4第37页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册Combin14 弹簧阻尼单元，可用于一维、二维、三维中的轴向弹簧阻尼单元，可用于一维、二维、三维中的轴向拉压和扭转。拉压和扭转。具有轴向拉压、扭转能力。具有轴向拉压、扭转能力。单元选项：单元选项：K1-
30、线性、非线性求解选项；线性、非线性求解选项；K21D自由度选项（包括平动自由自由度选项（包括平动自由 度和转动自由度）度和转动自由度）K3-2D、3D自由度选项自由度选项 （包括平动自由度和转动自由度）（包括平动自由度和转动自由度）实常数：实常数：K-弹簧刚度；弹簧刚度；CV1-阻尼系数；阻尼系数；CV2-非线性阻尼系数。非线性阻尼系数。ANSYS常用单元介绍常用单元介绍2000，4第38页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS常用单元介绍常用单元介绍BEAM4 三维弹性梁单元三
31、维弹性梁单元 可承受轴向拉压、扭转和弯曲，每个节点六个自由度可承受轴向拉压、扭转和弯曲，每个节点六个自由度 （三个平动（三个平动和三个转动自由度）和三个转动自由度）具有应力刚化和大变形功能。具有应力刚化和大变形功能。单元选项：单元选项：k2刚化矩阵选项（主矩阵或一致矩阵）；刚化矩阵选项（主矩阵或一致矩阵）；k6是否输出单元力和力矩；是否输出单元力和力矩；k7是否计算回转矩阵；是否计算回转矩阵；k9额外中间选项输出；额外中间选项输出；k10按长度单位还是按长度比例进行载按长度单位还是按长度比例进行载 荷偏移处理。荷偏移处理。实常数定义：可以通过截面参数按截面图形输入和实常数定义：可以通过截面参数
32、按截面图形输入和 定义截面尺寸、截面类型。定义截面尺寸、截面类型。同类型的单元：同类型的单元：BEAM3、BEAM23、BEAM24、BEAM44、BEAM188、BEAM189等。等。2000，4第39页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS常用单元介绍常用单元介绍PLANE42 2D结构实体单元。结构实体单元。可用于平面应力、平面应变、轴对称。有四个节点，每个节可用于平面应力、平面应变、轴对称。有四个节点，每个节点两个平动自由度。点两个平动自由度。具有塑性、蠕变、膨胀、应力刚
33、化、大变形和大应变功能。具有塑性、蠕变、膨胀、应力刚化、大变形和大应变功能。单元选项：单元选项：K1单元坐标系定义（平行于总体坐标系或以单元坐标系定义（平行于总体坐标系或以IJ边为基边为基准）；准）；K2是否包括额外的位移形状；是否包括额外的位移形状；K3-单元行为（平面应力、应变、轴对称等）单元行为（平面应力、应变、轴对称等）K4额外的应力输出（没有、积分点、节点）额外的应力输出（没有、积分点、节点）K5额外的表面输出选项。额外的表面输出选项。实常数：无实常数要求（单元厚度按单位厚度考虑）实常数：无实常数要求（单元厚度按单位厚度考虑）对应的高阶单元是对应的高阶单元是PLANE82，接受非轴对
34、称载荷的接受非轴对称载荷的 轴对称单元是轴对称单元是PLANE25。同类型的单元还有同类型的单元还有PLANE182、PLANE183等。等。2000，4第40页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN初级培训ANSYS常用单元介绍SOLID95 D20节点结构实体单元。节点结构实体单元。是是3D八节点实体单元八节点实体单元SOLID45的高阶单元。该单元在保证的高阶单元。该单元在保证精度的同时允许使用不规则的形状。具有相容的位移形式，精度的同时允许使用不规则的形状。具有相容的位移形式，适用于曲线边
35、界的建模。适用于曲线边界的建模。具有塑性、蠕变、膨胀、应力刚化、大变形和大应变功能。具有塑性、蠕变、膨胀、应力刚化、大变形和大应变功能。单元选项：单元坐标系；额外单元输出选项；额外表面输出选项；单元选项：单元坐标系；额外单元输出选项；额外表面输出选项；积分原则。积分原则。同类型的单元还有：同类型的单元还有：SOLID185、SOLID186等。等。第41页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN初级培训ANSYS常用单元介绍SHELL93 八节点结构壳单元。八节点结构壳单元。每个节点有六个自由度（三
36、个平动和三个转动自每个节点有六个自由度（三个平动和三个转动自 由由 度）。在平面的两个方向上，变形的形状都是度）。在平面的两个方向上，变形的形状都是 二次的。二次的。具有塑性、应力刚化、大变形和大应变功能。具有塑性、应力刚化、大变形和大应变功能。单元选项：单元坐标系；额外应力输出选项；单元选项：单元坐标系；额外应力输出选项；非线性积分点输出选项等。非线性积分点输出选项等。实常数：单元各点的厚度；单元实常数：单元各点的厚度；单元X轴的转角；轴的转角；增加的单位质量面等。增加的单位质量面等。同类型的单元：同类型的单元：第42页，共68页。CFD Analysis with ANSYS/FLOTRA
37、NCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS常用单元介绍常用单元介绍SHELL63 线性单元（线性单元（4节点），节点），具有应力刚化、大变形功能。具有应力刚化、大变形功能。SHELL143 四节点塑性壳单元，具有塑性、蠕变、四节点塑性壳单元，具有塑性、蠕变、膨胀、应力刚化、大变形和大应变功能。膨胀、应力刚化、大变形和大应变功能。可处理薄壳和中厚度的壳单元。可处理薄壳和中厚度的壳单元。SHELL181 有限应变壳单元。特别适用于具有线性、有限应变壳单元。特别适用于具有线性、大角度转动大角度转动/或非线性大应变特性的应用问题。或非线性大应变特性的应用问题。非
38、线性分析中考虑了壳厚度的变化。非线性分析中考虑了壳厚度的变化。单元内积分可用完全积分和缩减积分。单元内积分可用完全积分和缩减积分。采用采用SHELL43有收敛困难的问题可用有收敛困难的问题可用 SHELL181代替。代替。SHELL43 塑性大应变壳单元。塑性大应变壳单元。2000，4第43页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN初级培训ANSYS常用单元介绍SURF153 2D表面效应单元表面效应单元SURF154 3D表面效应单元表面效应单元 单元选项：单元特点（平面应力、应变、轴对称等），单
39、元选项：单元特点（平面应力、应变、轴对称等），是否有中间节点、压力载荷法向等。是否有中间节点、压力载荷法向等。实常数：基础刚度、表面张力、厚度实常数：基础刚度、表面张力、厚度第44页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册ANSYS常用单元介绍常用单元介绍接触单元：接触单元：目标面；接触面；目标面；接触面；可以通过接触导向来定义接触对；可以通过接触导向来定义接触对；单元选项：单元自由度、接触算法、接触检查点、单元选项：单元自由度、接触算法、接触检查点、CNOF的自动调节、接触时间（载荷）预报
40、、的自动调节、接触时间（载荷）预报、伪接触预报、接触初期渗透伪接触预报、接触初期渗透/间隙考虑否、间隙考虑否、接触刚度修正。接触刚度修正。实常数：接触刚度（实常数：接触刚度（KFN）、）、穿透容差值（穿透容差值（FTOLN）、）、摩擦系数等。摩擦系数等。2000，4第45页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN初级培训表面效应单元表面效应单元表面效应单元（表面效应单元（Surf154Surf154）应用简介应用简介单元定义：单元定义：4-84-8个节点及材料特性个节点及材料特性单元退化单元退化实常
41、数实常数 厚度等，厚度缺损为单位厚度。厚度等，厚度缺损为单位厚度。质量和体积的计算要使用单元厚度。质量和体积的计算要使用单元厚度。ADMSUADMSU，附加质量（单位面积）附加质量（单位面积）EFSEFS，基础刚度基础刚度 .第46页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册表面效应单元表面效应单元表面效应单元（表面效应单元（Surf154)Surf154)的加载的加载节点和单元载荷节点和单元载荷对于对于faceface为为1 1，2 2，3 3，4 4，5 5压力的正值作用在单元坐标系的正面
42、（除压力的正值作用在单元坐标系的正面（除非在非在Z Z的负方向作用法向压力）；的负方向作用法向压力）；正或负根据正或负根据Keyopt(6)Keyopt(6)可以移动到流体约束的自由表面；可以移动到流体约束的自由表面；face=1face=1，压力作用在表面的法向；，压力作用在表面的法向；face=2face=2，压力作用在表面的切向（，压力作用在表面的切向（X)X)；face=3face=3，压力作用在表面的切向（，压力作用在表面的切向（Y)Y)；2000，4第47页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FL
43、OTRAN培训手册表面效应单元表面效应单元Face=4Face=4，法向压力法向压力 在每个积分点，压力的大小是在每个积分点，压力的大小是Pi+xPj+yPk+ZPlPi+xPj+yPk+ZPl 其中，其中，Pi,Pj,Pk,PlPi,Pj,Pk,Pl是用是用SFESFE命令输入的命令输入的 VAL1,VAL2,VAL3,VAL4 VAL1,VAL2,VAL3,VAL4，X,Y,ZX,Y,Z是当前是当前 位置的笛卡儿系下的坐标值位置的笛卡儿系下的坐标值2000，4第48页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/F
44、LOTRAN培训手册表面效应单元表面效应单元Face=5Face=5，与表面成一定的角度的压力，与表面成一定的角度的压力 压力值的大小是压力值的大小是Pi Pi；压力的方向是：压力的方向是：（P PjjX+PX+Pk kY+PlZ)/(PjY+PlZ)/(Pj2 2+Pk+Pk2 2+Pl+Pl2 2)1/21/2压力值可以根据压力值可以根据Keyopts(11)Keyopts(11)和和Keyopts(12)Keyopts(12)调整，调整，但，当使用但，当使用SffunSffun或或SfgradSfgrad命令时，方向是不会调整的。命令时，方向是不会调整的。2000，4第49页，共68页。
45、CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册表面效应单元表面效应单元Suuf154Suuf154的输入的输入单元名单元名 SURF154SURF154节点：节点：i,j,k,l if keyopt(4)=1 i,j,k,l if keyopt(4)=1 （无中间节点）（无中间节点）I,j,k,l,m,n,o,p if keyopt(4)=0 I,j,k,l,m,n,o,p if keyopt(4)=0 （有中间节点）（有中间节点）自由度：自由度：ux uy uzux uy uz实常数：实常数：材料特性：材料
46、特性：Dens,Visc,DampDens,Visc,Damp2000，4第50页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册表面效应单元表面效应单元 Surface LoadsSurface Loads Pressures:Face=1 (I,J,K,L)Pressures:Face=1 (I,J,K,L)表面法向表面法向 Face=2 (I,J,K,L)Face=2 (I,J,K,L)切向（切向（+X)+X)Face=3 (I,J,K,L)Face=3 (I,J,K,L)切向（切向（+Y)+
47、Y)Face=4 (I,J,K,L)Face=4 (I,J,K,L)法向法向 （锥形（锥形)Face=Face=5 5 (I,J,K,L)(I,J,K,L)根据输入的矢量根据输入的矢量2000，4第51页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册*ANSYS单元分类单元分类*ANSYS常用单元常用单元 介绍介绍 ANSYS常用单元推荐常用单元推荐ANSYS单元介绍单元介绍2000，4第52页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis wi
48、th ANSYS/FLOTRAN培训手册实体单元推荐实体单元推荐几乎不可压缩材料几乎不可压缩材料(塑性塑性)可忽略弯曲的体结构变形采用可忽略弯曲的体结构变形采用 Plane182,Solid185Plane182,Solid185 选择选择缩减积分缩减积分(B-Bar)B-Bar)。对于小应变应用采用非协调模式单元对于小应变应用采用非协调模式单元 Plane42,Solid45Plane42,Solid45。对于大应变应用采用具有对于大应变应用采用具有URI(URI(特别对于大模型特别对于大模型)的的Plane182Plane182 和和 Solid185Solid185 或具有或具有URIU
49、RI的的 Solid95Solid95 。也可以采用也可以采用 Visco106,Visco106,Visco107 Visco107 和和 Visco108Visco108 单元单元(甚至对于与速率无关的塑性甚至对于与速率无关的塑性)。2000，4第53页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册实体单元推荐实体单元推荐不可压缩材料不可压缩材料(超弹性超弹性)类似橡胶材料的不可压缩条件约束需要类似橡胶材料的不可压缩条件约束需要 U-P U-P 混合公式，可采用混合公式，可采用单元单元 Hyp
50、er56,Hyper58,Hyper158 Hyper56,Hyper58,Hyper158 或或 Hyper74 Hyper74。2000，4第54页，共68页。CFD Analysis with ANSYS/FLOTRANCFD Analysis with ANSYS/FLOTRAN培训手册壳单元壳单元-概述概述 当结构的总体厚度相对于典型长度很小时可使用壳单元，长度当结构的总体厚度相对于典型长度很小时可使用壳单元，长度比厚度大十倍以上的问题可决定使用壳单元。比厚度大十倍以上的问题可决定使用壳单元。公开的文献中有各种壳理论，这源于对壳位移的不同近似描述公开的文献中有各种壳理论，这源于对壳位
