Ansys APDL / Example —— 法兰连接

目 录

  • 〇. 模型说明
  • 一. 前处理
    • 1. 设置参数
      • 1.1 螺栓垫片参数
      • 1.2 法兰参数
      • 1.3 加载点参数
      • 1.4 输出文件
    • 2. 创建单元
      • 2.1 创建上法兰外圈通孔环带单元
        • 2.1.1 设置单元属性及材料性质
        • 2.1.2 创建局部坐标系
        • 2.1.3 创建几何图形
        • 2.1.4 生成网格
      • 2.2 创建上法兰内圈通孔环带单元
      • 2.3 创建上法兰环带单元 (内圈+外圈 绑定)
        • 2.3.3 定义实常数及接触单元类型
        • 2.3.4 创建接触单元
      • 2.4 创建下法兰外圈盲孔环带单元
        • 2.4.3 体的创建、切分及删除
        • 2.4.3 划分网格
      • 2.5 创建下法兰内圈通孔环带单元
      • 2.6 创建下法兰环带单元 (内圈+外圈 绑定)
      • 2.7 创建外圈盲孔螺栓连接
        • 2.7.1 创建垫片单元
        • 2.7.2 创建螺栓杆单元
        • 2.7.3 创建螺栓间接触
        • 2.7.4 创建螺栓垫片间接触
        • 2.7.2.5 创建预紧单元
      • 2.8 创建内圈通孔螺栓连接
      • 2.9 创建垫片与法兰间的接触
    • 3. 装配与约束创建
    • 4. 施加边界条件
  • 二. 加载与求解
  • 三. 后处理

〇. 模型说明

简单网格由ANSYS生成,复杂网格由HyperMesh处理。

Mechanical APDL Product Launcher

一. 前处理

1
2
3
4
5
6
finish                      ! 退出当前处理器
/clear,all                  ! 清除所有
/prep7                      ! 进人前处理器
*afun,rad                   ! 指定角度单位为弧度
pi = acos(-1)               ! 获取圆周率pi的数值
*afun,deg                   ! 将角度单位切换为度

1. 设置参数

单位:mm 。

1.1 螺栓垫片参数

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
Ds = 2500                                   ! 螺栓分度圆直径
num = 80                                    ! 螺栓个数

d_bolt = 42                                 ! 螺栓公称直径 ―― d
d_hole = 45                                 ! 螺栓孔直径 ―― d0
l_bolt = 260                                ! 螺栓长度 ―― l
l_b = 20                                    ! 旋合长度

de_bolt = sqrt(4*1120/pi)                   ! 内外圈螺栓的有效直径 ―― de
ds_bolt = 38                                ! 螺栓光杆截面直径 ―― ds


d_washer_inside=45                          ! 垫片内径
d_washer_outside=78                         ! 垫片外径
t_washer=8                                  ! 垫片厚度

1.2 法兰参数

1
2
3
4
5
6
D_flange = 2800                             ! 法兰外径
t_flange = 80                               ! 法兰厚度
a = 150                                     ! 螺栓边距

t = 20                                      ! 筒体壁厚
width = 100                                 ! 螺栓孔附近法兰环带宽度

1.3 加载点参数

1
L = 2000                                    ! 加载点与法兰连接面的距离

1.4 输出文件

1
PARSAV,SCALAR ,'parameters','txt',' '       ! 将上述定义的参数以ANSYS能识别的方式存储到parameters.txt文件中

??ANSYS中不区分大小写,通常情况下,为了提高命令流的可读性,约定:命令流中小写表示用户自定义的变量,大写为ANSYS系统固有变量。

??运行上述命令流后,定义的参数可在ANSYS界面中,依次点击 菜单栏 >> Parameters >> Scalar Parameters 中查看,如下所示:

在这里插入图片描述

2. 创建单元

1
2
3
4
5
finish                                 ! 退出当前处理器
/clear,all                             ! 清除所有
/prep7                                 ! 进人前处理器
*afun,deg                              ! 指定角度单位为度
PARRES,NEW,'parameters','txt',' '      ! 从文件para1.txt中读入参数

2.1 创建上法兰外圈通孔环带单元

2.1.1 设置单元属性及材料性质

  • SOLID186 —— 3-D 20-Node Structural Solid,DOF: UX, UY, UZ。

??KEYOPT(2)
????Element technology:
????0 – Uniform reduced integration (default)
????1 – Full integration

??分别对应C3D8R和C3D8单元。

  • MESH200 —— Meshing Facet 2-20 nodes 2-D/3-D space,DOF: None。

????MESH200 is a “mesh-only” element, contributing nothing to the solution

????不参与求解计算,生成高维度网格时所需要的低维度网格。

????部分KEYOPT(1)的取值,如下:

在这里插入图片描述

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
*GET,etmax,ETYP,,NUM,MAX     ! 获取当前定义的最大单元类型编号并将其赋予变量etmax

! (1). 定义单元类型 —— 3D 20节点 完全积分单元

ET,etmax+1,SOLID186          ! 定义单元类型,单元类型参考号为当前最大单元类型号加1,SOLID186单元。
KEYOPT,etmax+1,2,1           ! 设置单元关键选项,将单元类型参考号为etmax+1的单元设置为完全积分单元。

! (2). 定义线弹性各项同性材料

MPTEMP,,,,,,,,               ! 为材料属性定义温度表
MPTEMP,1,0                   ! 为材料属性定义温度表
MPDATA,EX,etmax+1,,1.69e5    ! 定义材性,材料参考编号设置为etmax+1,弹性模量EX=1.69e5。
MPDATA,PRXY,etmax+1,,0.275   ! 定义材性,材料参考编号为etmax+1,泊松比设置为PRXY=0.275。

! (3). 定义网格单元

ET,etmax+2,200               ! 定义单元类型
KEYOPT,etmax+2,1,7           ! 设置单元关键选项

! (4). 为单元赋属性

EMODIF,ALL,MAT,etmax+1       ! 对已存在单元进行修改,实际上,实现的是为所有单元赋材料属性。

??上述命令流定义的单元类型及材料属性,在Ansys界面中,按如下方式查看:

在这里插入图片描述

单元类型的查看

在这里插入图片描述

材料属性的查看

2.1.2 创建局部坐标系

1
2
3
CSYS,0                       ! 激活总体直角坐标系统
LOCAL,15,0,L,0,0,0,0,90      ! 建立以(L,0,0)为原点,局部z轴与总体x轴平行,坐标系号为15的局部直角坐标。
LOCAL,16,1,L,0,0,0,0,90      ! 建立以(L,0,0)为原点,局部z轴与总体x轴平行,坐标系号为16的局部柱角坐标。

坐标系统的有关命令

??ANSYS中的坐标轴:x、y、z —— 白、绿、蓝

  • CSYS

??使用功能:激活一个预先已定义的坐标系统。
??使用格式:CSYS,KCN,其中,KCN为将要激活的坐标系统参考号。

KCN 说明
0 直角坐标 (默认)
1 用z轴作为旋转轴的柱坐标系
2 球坐标系
4 或 WP 工作平面
5 用y轴作为旋转轴的柱坐标系
11及以上 由用户已预先定义的局部坐标系统
  • LOCAL

??使用功能:通过位置和方向指定一个局部坐标系。
??使用格式:LOCAL,KCNKCSXCYCZCTHXYTHYZTHZXPAR1PAR2

其中,KCN为用户给定的大于等于11的坐标系统参考号,作为识别该局部坐标系的标志;
???KCS:坐标系的类型;
???XCYCZC:新坐标系原点位置 (在直角坐标系统里) Location (in the global Cartesian coordinate system) of the origin of the new coordinate system;
???THXYTHYZTHZX:分别绕z轴、x轴和y轴旋转的角度,其正向为:xy,yz,zx。

KCS 说明
0 / CART 直角坐标 (Cartesian)
1 / CYLIN 柱坐标 (Cylindrical (circular or elliptical))
2 / SPHE 球坐标 (Spherical (or spheroidal))
3 / TORO Toroidal
  • CLOCAL

??使用功能:Defines a local coordinate system relative to the active coordinate system。
??使用格式:LOCAL,KCNKCSXLYLZLTHXYTHYZTHZXPAR1PAR2

其中,KCN为用户给定的大于等于11的坐标系统参考号,作为识别该局部坐标系的标志;
???XLYLZL:Location (in the active coordinate system) of the origin of the new coordinate system (R, θ, Z for cylindrical, R, θ,Φ for spherical or toroidal).;
???THXY:First rotation about local Z (positive X toward Y);
???THYZ:Second rotation about local X (positive Y toward Z);
???THZX:Third rotation about local Y (positive Z toward X)。

??GUI中局部坐标系统的创建/删除如下:

在这里插入图片描述

局部坐标系统的创建与删除

??GUI中模型坐标系统的查看如下:

在这里插入图片描述

坐标系统的查看

??GUI中当前激活的坐标系统的查看如下:

在这里插入图片描述

当前被激活坐标系统的查看

2.1.3 创建几何图形

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
*GET,kmax,KP,,NUM,MAX, , , ,         ! 获取最大关键点编号将其赋予变量kmax

!(1). 定义局部坐标系40 —— 方便螺栓孔的创建

CSYS,15                              ! 激活局部坐标系15
CLOCAL,40,1,Ds/2,,0,0,0,0            ! 根据激活的局部坐标系15定义新的局部柱坐标系40;
                                     ! 原点为(Ds/2,0,0),各轴方向与坐标系15一致;
                                     ! 局部坐标系40随之被激活。

! (2). 创建螺栓孔的半圆轮廓线

WPCSYS,-1,                            ! 根据当前坐标系定义工作平面的位置
K,1,d_hole/2,0                        ! 创建关键点1,其在激活的局部坐标系40中的坐标为(d_hole/2,0,0)。
K,2,d_hole/2,45                       ! 创建关键点2,其在激活的局部坐标系40中的坐标为(d_hole/2,45,0)。
K,3,d_hole/2,90                       ! 创建关键点3,其在激活的局部坐标系40中的坐标为(d_hole/2,90,0)。
K,4,d_hole/2,135                      ! 创建关键点4,其在激活的局部坐标系40中的坐标为(d_hole/2,135,0)。
K,5,d_hole/2,180                      ! 创建关键点5,其在激活的局部坐标系40中的坐标为(d_hole/2,180,0)。

! K, 关键点编号, x坐标, y坐标, z坐标
! 缺省时为0;创建关键点时,可指定关键点的编号。

L,1,2                                 ! 由两个关键点创建一条线(可以是直线或是弧线取决于坐标系)
*REPEAT,4,1,1                         ! 上一条命令重复执行4次,始末关键点编号增量为1
                                      ! eg. 创建线,线的起点、终点编号分别为2、3。
! 创建线时,系统自动为线分配编号。


! L,1,2
! *REPEAT,4,1,1
! 以上两行命令相当于:
! L,1,2
! L,2,3
! L,3,4
! L,4,5

! (3). 创建法兰环带轮廓线
CSYS,15                              ! 激活局部坐标系15
CLOCAL,42,1,0,0,0                    ! 根据激活的局部坐标系16定义新的局部柱坐标系42
CSYS,16                              ! 激活局部坐标系16
CLOCAL,42,1,0,0,0                    ! 根据激活的局部坐标系16定义新的局部柱坐标系42
                                     ! 局部柱坐标系42随之被激活
                                                                       
WPCSYS,-1,                                      ! 根据当前坐标系定义工作平面的位置
k,6,Ds/2+width/2                                ! 创建关键点
k,7,Ds/2+width/2,180/num                        ! 创建关键点
k,8,Ds/2,180/num                                ! 创建关键点
k,9,Ds/2-width/2,180/num                        ! 创建关键点
k,10,Ds/2-width/2,0                             ! 创建关键点

L,6,7
*REPEAT,4,1,1
L,1,6
*REPEAT,5,1,1

! (4). 由轮廓线生成面 Area from lines ---> AL

AL,1,10,5,9                   ! 创建平面,平面由边界线围成,边界线的编号分别为1、10、5、9。
*REPEAT,4,1,1,1,1             ! 上一条命令重复执行4次,各边界线编号增量均为1

??以上命令流运行后,将得到:

在这里插入图片描述

点、线和面的创建

  • WPCSYS

??使用功能:根据当前坐标系来定义工作平面的位置
??使用格式:WPCSYS,WNKCN

其中:
??WN : Window number whose viewing direction will be modified to be normal to the working plane (defaults to1). If WN is a negative value, the viewing direction will not be modified.
??KCN为坐标系统参考号,默认值是当前激活坐标系。

??几何在工作平面上绘制,工作平面只能为当前坐标系的xy平面,z方向为其法线方向,因此 WN 确定视线来自z轴正方向还是负方向。

??GUI中工作平面的定义如下:

在这里插入图片描述

工作平面的定义

  • *REPEAT

??使用功能:重复以前的命令,主要是重复执行上一条命令。
??使用格式:*REPEAT,NTOT,VINC1,VINC2,VINC3,VINC4,VINC5,VINC6,VINC7,VINC8,VINC9,VINC10,VINC11

其中,NTOT : 包括第1次在内,前述命令被执行的次数,若NTOT=4,则上一条命令还会被重复执行3次;
???VINC1,VINC2,…,VINC11 : 施加到前述命令上第1个~第11个数据域的值增量。

  • AL

??使用功能:通过已定义的边界线生成一个面。
??使用格式:AL,L1,L2,L3,L4,L5,L6,L7,L8,L9,L10

其中,L1,L2,…,L10 为线的编号,至少需要3条线才能生成一个面。

2.1.4 生成网格

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
! (1) 由几何面生成面网格

LESIZE,9, , ,2, , , , ,0        ! 设置线的网格大小,编号为9的线被分成2段。
*REPEAT,5,1,,,,                 ! 上一条命令重复执行5次
LESIZE,1, , ,2, , , , ,0        ! 设置线的网格大小,编号为1的线被分成2段。
*REPEAT,4,1,,,,                 ! 上一条命令重复执行4次
LESIZE,8, , ,2, , , , ,0        ! 设置线的网格大小,编号为8的线被分成2段。

MSHAPE,0,2D                     ! 指定划分单元的形状,生成四边形单元。
MSHKEY,1                        ! 采用映射网格划分方式划分网格。
AMESH,ALL                       ! 将所有面划分网格

! (2) 由面网格拉伸成体网格

TYPE,etmax+1                      ! 给随后生成的单元激活一个单元类型号
EXTOPT,ESIZE,8,1,                 ! EXTOPT:由面单元生成体单元的控制选项;
                                  ! 8表示在体生成或体扫略方向上单元分割数量为8;
                                  ! 1表示在体生成或体扫略方向上的间隔率为1(默认)。
EXTOPT,ACLEAR,1                   ! 体单元网格生成后清除面单元网格。
VEXT,ALL, , ,0,0,t_flange         ! 通过给的偏移量由面生体


! (3) 网格镜像及阵列生成法兰环带全部单元

! ① 对称单元
CSYS,15                                           ! 激活局部直角坐标系15
Alls                                              ! 选择所有实体
VSYMM,Y,ALL, , , ,0,0                             ! 对体进行镜像操作

! ② 阵列单元
CSYS,16                                           ! 激活局部柱坐标系16(原点位于螺栓群形心处)
VGEN,num,ALL, , , ,360/num, , ,0                  ! 复制体,重复生成次数num,复制全部体,
                                                  ! 当前激活坐标系下,关键点坐标值的偏移量为(--,360/num,0)
                                                  ! 单个偏转角度为360/num

NUMMRG,ALL                                        ! 对实体进行合并
NUMCMP,ALL                                        ! 压缩所定义项的编号

CM,e_flange_throughHole,elem                      ! 由所选实体生成一个元件
                                                  ! 元件名称:e_flange_throughHole
                                                  ! 元件中的数据类型:elem
ALLSel,ALL                                        ! 选择所有实体
CDWRITE,DB,'Topf_TH_Outside','cdb',,'',''         ! 输出文件,文件名:Topf_TH_Outside,格式:.cdb。

! .db 文件 —— 二进制数据库文件
! .cdb 文件 —— 有限元求解文件

??以上命令流将实现,由面生成体网格,在法兰厚度方向上有8层单元,如下图所示:

在这里插入图片描述

创建的法兰环带单元 (局部)

  • LESIZE

??使用功能:对所选择的线设置网格单元大小。
??使用格式:LESIZE,NL1SIZEANGSIZNDIVSPACEKFORCLAYER1LAYER1KYNDIV

其中:

??NL1 : 线的编号,也可以为 ALL、P 或元件名;
??SIZE : 若 NDIV 为空,SIZE 指单元的边长,分段数将根据线长自动计算,圆整得到;
??NDIV : 若为正,NDIV 表示每条线的分段数;
??KYNDIV : 当 KYNDIV 为0时,表示智能化网格划分无效,如果等分不匹配,映射网格失败;当 KYNDIV 为1时,对于大曲率或邻接区域,智能化网格优先使用。

  • VSYMM

??使用功能:通过对坐标轴对称镜像生成体
??使用格式:VSYMM,NcompNV1NV2NINCKINCNOELEMIMOVE

其中:

??Ncomp : 对称控制选项,Ncomp = X,Y 或 Z,镜像平面的法线方向;
??NV1NV2NINC: 按增量 NINCNV1NV2 将要进行对称镜像的体编号的范围,其中 NV1 也可以为 ALL、P 或元件名;
??NOELEM: 是否也生成单元和节点的控制键,若为0,如果存在相关的节点和单元,也生成相关的节点和单元,若为1,不生成节点和单元。
??IMOVE: 设置是否删除/新定义体,若为0,保留原来的体,生成额外的体,若为1,删除原来的体,创建新的体。

??相类似的操作命令有:KSYMM、LSYMM、ARSYM。

  • VGEN

??使用功能:对体进行复制操作,在不同的坐标系中可实现平移阵列和环形阵列;
??使用格式:VGEN,ITIMENV1NV2NINCDXDYDZKINCNOELEMIMOVE

其中:

??ITIME : 重复生成的次数,默认值为2;
??NV1NV2NINC: 按增量 NINCNV1NV2 将要进行对称镜像的体编号的范围,其中 NV1 也可以为 ALL、P 或元件名;
??DXDYDZ: 在当前激活坐标系中,体中各关键点坐标值的偏移量;
??KINC: 要生成关键点编号的增量,如果其值为0,由系统自动确定;
??NOELEM: 是否也生成单元和节点的控制键,若为0,如果存在相关的节点和单元,也生成相关的节点和单元,若为1,不生成节点和单元。
??IMOVE: 设置是否删除/新定义体,若为0,保留原来的体,生成额外的体,若为1,删除原来的体,创建新的体。

  • TYPE

??使用功能:设置单元类型属性指示器,为随后生成的单元指定一个单元类型号。
??使用格式:TYPE,ITYPE,其中,ITYPE 为单元类型号。

  • EXTOPT

??使用功能:由面单元生成体单元的相关控制选项。
??使用格式:EXTOPT,LabVal1Val2Val3

其中,Lab 为识别控制选项的标签;Val1Val2Val3将根据Lab 的不同而变化。

  • VEXT

??使用功能:通过给定偏移量由面生体。
??使用格式:VEXT,NA1NA2NINCDXDYDZRXRYRZ

其中:

??NA1NA2NINC:设置将要被拖拉面的范围,即按增量 NINCNA1 增大到 NA2,其中,NA1 也可以为 ALL、P 或元件名;
??DXDYDZ:为在激活坐标系中,作用于关键点坐标值在X,Y和Z方向的增量 (其他坐标系下类似);
??RXRYRZ:为在激活坐标系中,作用于关键点坐标值在X,Y和Z方向的缩放因子。

  • NUMMRG

??使用功能:对实体进行合并 / Merges coincident or equivalently defined items 。
??使用格式:NUMCMP,LabelsTOLERGTOLERActionSwitch

其中,Labels : 合并操作选项,它可以为:NODE(节点)、ELEM(单元)、KP(关键点)、MAT(材料编号)、TYPE(单元类型)、REAL(实常数)、CP(耦合)、CE(约束方程) 或 ALL(以上全部)。

  • NUMCMP

??使用功能:压缩所定义项的编号。
??使用格式:NUMCMP,Labels

  • CM

??使用功能:由所选实体生成一个组件 / Groups geometry items into a component.
??使用格式:CM,CnameEntity

其中:
??Cname: 创建的组件名称;
??Entity: 组件中的数据类型,它的值可以为:VOLU、AREA、LINE、KP、ELEM 或 NODE;

  • ALLSEL

??使用功能:选择所有实体 / Selects all entities with a single command.
??使用格式:ALLSEL,LabTEntity

其中:
??LabT: 选择的类型,若为 ALL,选择所有指定的实体类型及其下层附属实体(默认方式),若为 BELOW,选择所有直接相关的实体项及其下层的实体类型;
??Entity: 选择所依据的实体类型,它的值可以为:ALL(所有的实体类型)(默认方式)、VOLU、AREA、LINE、KP、ELEM 或 NODE;

2.2 创建上法兰内圈通孔环带单元

??按照上节步骤,可创建上法兰内圈通孔环带单元。

2.3 创建上法兰环带单元 (内圈+外圈 绑定)

??通过在上法兰外圈环带和内圈环带间建立约束,形成整体的上法兰环带单元。

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
FINISH                             ! 退出当前处理器
/CLEAR,ALL                         ! 清除所有
/prep7                             ! 进人前处理器
*AFUN,DEG                          ! 指定角度单位为度
PARRES,NEW,'parameters','txt',' '  ! 从文件para1.txt中读入参数

! 1. 导入cdb文件

CDread,db,'Topf_TH_Inside','cdb',,'',''   ! 读取并导入文件Topf_TH_Inside.cdb
CDread,db,'Topf_TH_Outside','cdb',,'',''  ! 读取并导入文件Topf_TH_Outside.cdb

! 2. 定义坐标系

CSYS,0                         ! 激活全局直角坐标系
LOCAL,15,0,L,0,0,0,0,90        ! 局部直角坐标系15
LOCAL,16,1,L,0,0,0,0,90        ! 局部柱坐标系16

2.3.3 定义实常数及接触单元类型

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
CSYS,16                                  ! 激活局部坐标系16
*GET,etmax,ETYP, ,NUM,MAX                ! 获取当前定义的最大单元类型编号将其赋予参数etmax
*GET,rlmax,RCON, ,NUM,MAX                ! 获取最大实常数将其赋予参数rlmax

! (1). 定义实常数
R,rlmax+1                                ! 定的单元的实常数,实常数编号为RLMAX+1。

! (2). 定义目标面单元
ET,etmax+1,TARGE170                      ! 定义目标面单元类型,单元类型:TARGE170,编号:etmax+1。

! (3). 定义接触面单元
ET,etmax+2,CONTA174                      ! 定义接触面单元类型,单元类型:CONTA174,编号:etmax+2。  
KEYOPT,etmax+2,2,2                       ! 设置单元的关键选项
KEYOPT,etmax+2,4,1                       ! 设置单元的关键选项
KEYOPT,etmax+2,5,3                       ! 设置单元的关键选项
KEYOPT,etmax+2,9,1                       ! 设置单元的关键选项
KEYOPT,etmax+2,12,5                      ! 设置单元的关键选项

??ANSYS 支持刚体/柔体的面-面的接触单元,刚性面被当做“目标”面(Abaqus中的主面,粗网格硬材料为主面),分别用 Target169 和 Target170 来模拟 2D 和 3D 的目标面,柔性体的表面被当做“接触”面(Abaqus中的从面),用 Conta171,Conta172,Conta173,,Conta174 来模拟。一个目标单元和一个接触单元叫做一个“接触对”,程序通过一个共享的实常数来识别“接触对”,给目标单元和接触单元指定相同的实常的号。

??在ANSYS GUI 中,目标面显示为粉色,接触面显示为绿色,如下所示:

在这里插入图片描述

左侧主面 右侧从面

??Abaqus使用单纯的主-从接触算法:从面上的节点不能侵入主面的任何部分。该算法对主面没有做限制,主面可以在从面的节点之间侵入从面。

??ANSYS 中几种接触单元如下:

在这里插入图片描述

在这里插入图片描述

  • CONTA174 KEYOPT

??设置关键选项的格式:KEYOPT,ITYPEKNUMVALUE

??(2). KEYOPT(2) —— Contact algorithm

0 1 2 3 4
Augmented Lagrangian (default) Penalty function Multipoint constraint (MPC) Lagrange multiplier on contact normal and penalty on tangent Pure Lagrange multiplier on contact normal and tangent

??(4). KEYOPT(4) —— Location of contact detection point

0 1 2 3
On Gauss point (for general cases) On nodal point - normal from contact surface On nodal point - normal to target surface On nodal point - normal from contact surface (projection-based method)

??(5). KEYOPT(5) —— CNOF/ICONT Automated adjustment

0 1 2 3 4
No automated adjustment Close gap with auto CNOF Reduce penetration with auto CNOF Close gap/reduce penetration with auto CNOF Auto ICONT

??(9). KEYOPT(9) —— Effect of initial penetration or gap

0 1 2 3 4 5
Include both initial geometrical penetration or gap and offset Exclude both initial geometrical penetration or gap and offset Include both initial geometrical penetration or gap and offset, but with ramped effects Include offset only (exclude initial geometrical penetration or gap) Include offset only (exclude initial geometrical penetration or gap), but with ramped effects Include offset only (exclude initial geometrical penetration or gap) regardless of the initial contact status (near-field or closed)

??(10). KEYOPT(10) —— Contact stiffness update

0 1 2
Each iteration based on the current mean stress of underlying elements. The actual elastic slip does not to exceed the maximum allowable limit (SLTO) within a substep Each load step if FKN is redefined during the load step. Each iteration based on the current mean stress of underlying elements. The actual elastic slip never exceeds the maximum allowable limit (SLTO) during the entire solution.

??(12). KEYOPT(12) —— Behavior of contact surface

0 1 2 3 4 5 6
Standard Rough No separation (sliding permitted) Bonded No separation (always) Bonded (always) Bonded (initial contact)
  • TARGE170 KEYOPT

??设置关键选项的格式:KEYOPT,ITYPEKNUMVALUE

??(2). KEYOPT(2) —— Boundary conditions for rigid target nodes

??0 – Automatically constrained by the program

??1 – Specified by user

  • TSHAP

??使用功能:为目标单元 TARGE169 和 TARGE170 指定 2D 和 3D 几何表面。
??使用格式:TSHAP,shape

其中:
??shape:为目标单元 TARGE169 和 TARGE170 指定几何形状。其有效的形状有:PILO(2D 或3D 引导节点)。

??使用提示:使用该命令可以生产面与面接触时的刚性目标面,对于 2D 是 TARGE169 、CONTA171 和 CONTA172,对于 3D 是 TARGE170 、CONTA173 和 CONTA174 。执行该命令后,其后生成的单元具有同样的形状,直到用另一个形状来代替。

2.3.4 创建接触单元

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
REAL,rlmax+1                       ! 为随后生成的目标面和接触面单元指定实常数

! (1). 创建目标面单元

TYPE,etmax+1                       ! 为随后生成的单元激活一个单元类型号

! 定义选择集

CSYS,16
 
CMSEL,S,ele_Topf_TH_Outside        ! 选择组件,组件名称为ele_Topf_TH_Outside。

NSEL,,                                                   ! 选中依附于所选单元上的节点
NSEL,R,EXT                                               ! 选择所选单元的外部(表面)节点
NSEL,R,LOC,X,Ds/2-0.5*width-0.5,Ds/2-0.5*width+0.5       ! 选择x坐标值位于给定区间的所有节点
NSEL,R,LOC,Z,-0.5,t_flange+0.5                           ! 选择z坐标值位于给定区间的

ESLN                                                     ! 选中依附于所选节点上的单元
ESURF,,TOP                                               ! 在已选单元的自由面上生成重叠单元

! (2). 创建接触面单元

TYPE,ETMAX+2  

! 定义选择集

CSYS,16
                                                      ! 给随后生成的单元激活一个单元类型号
CMSEL,S,ele_Topf_TH_Inside                            ! 创建新选择集
                                                      ! 将组件ele_Topf_TH_Inside放入到新选组集中。
NSEL,,                                                ! 选中依附于上一步选择单元的节点
NSEL,R,EXT                                            ! 选择所选单元的外部(表面)节点
NSEL,R,LOC,X,Ds/2-0.5*width-0.5,Ds/2-0.5*width+0.5    ! 选择上一步所选节点x坐标值位于给定区间的节点
NSEL,R,LOC,Z,-0.5,t_flange+0.5                        ! 选择上一步所选节点z坐标值位于给定区间的节点

ESLN                                                  ! 选中依附于所选节点上的单元
ESURF,,TOP                                            ! 在已选单元的自由面上生成重叠单元

! (3). 输出文件

ALLSEL,ALL                                            ! 全选
CDWRITE,db,'Topflange_TH','cdb',,'',''                ! 输出文件Topflange_TH.cdb

! ele_Topflange_TH.cdb文件中包含约束信息

??由上所创建的上法兰内外圈环带单元如下:

在这里插入图片描述

2.4 创建下法兰外圈盲孔环带单元

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
FINISH                             ! 退出当前处理器
/CLEAR,ALL                         ! 清除所有
/prep7                             ! 进人前处理器
*AFUN,DEG                          ! 指定角度单位为度

PARRES,NEW,'parameters','txt',' '  ! 从文件para1.txt中读入参数
NUMCMP,ALL                         ! 压缩所定义项的编号

! 1. 设置单元属性及材料性质

! (1). 定义SOLID186单元

*GET,etmax,ETYP,,NUM,MAX
ET,etmax+1,SOLID186
KEYOPT,etmax+1,2,1

ET,etmax+3,SOLID186  ! 旋合段螺栓的单元
KEYOPT,etmax+3,2,1
 
! (2). 定义球墨铸铁材料属性及钢材材料属性

MPTEMP,1,0  
MPDATA,EX,etmax+1,,1.69e5      ! 指定弹性模量      
MPDATA,PRXY,etmax+1,,0.275     ! 指定泊松比

MPTEMP,1,0  
MPDATA,EX,etmax+3,,2e5
MPDATA,PRXY,etmax+3,,0.3

! (3). 定义MESH200单元

ET,etmax+2,200
KEYOPT,etmax+2,1,7

! (4). 为单元赋材料属性

EMODIF,ALL,MAT,etmax+1

! 2. 创建局部坐标系

CSYS,0
LOCAL,15,0,L,0,0,0,0,90           ! 建立局部直角坐标系15
LOCAL,16,1,L,0,0,0,0,90           ! 建立局部柱坐标系16
CLOCAL,32,0,0,0,-t_flange,0,0,0   ! 局部直角坐标系32被激活

2.4.3 体的创建、切分及删除

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
! (1). 创建螺栓扇区体

WPCSYS,-1,,  
CYL4,0,0,Ds/2-width/2,0,Ds/2+width/2,180/num,t_flange
        ! 生成圆柱体
        ! 圆心坐标(0,0);内径:Ds/2-width/2,起始角度:0;
        ! 外径:Ds/2+width/2,结束角度:360/num/2;
        ! 体的高度:t_flange。
       
CM,v1,VOLU        ! 创建元件,元件名称:v1,元件类型:volu;
                  ! 此时组件v1中有1个体

CYL4,Ds/2,0,d_hole/2,0,de_bolt/2,90,t_flange  
CYL4,Ds/2,0,d_hole/2,90,de_bolt/2,180,t_flange
CM,v2,VOLU ! 此时组件v2中有3个体

! (2). 布尔操作 Overlap

VOVLAP,ALL     ! 体搭接,实际上将体v1用v2切割成部分,原来的v1整体被删除。


CSYS,32                             ! 激活局部直角坐标系
CLOCAL,33,1,Ds/2,0,0,0,0,0          ! 创建局部柱坐标系,原点在分度圆上。

! (3). 线拖拉成面

K,17,d_hole/2,90,-50                          
K,18,d_hole/2,90,t_flange+50            
L,17,18                             ! 创建直线,直线的编号为L1,此为被拖拉直线。                    

CSYS,16   ! 激活局部柱坐标系16
K,21,Ds/2-width/2,90/num,,
K,24,Ds/2-width/2,360/num,,
L,21,24                             ! 创建直线,直线的编号为L8,此为路径线。

ADRAG,1, , , , , ,8             ! 由线沿路径线拖拉成面。
                                ! 将线L1沿路径L8拖拉生成面,生成的面编号为A1

! (4). 切分体
VSBA,6,1                        ! 用面分割体,被分割的体号为V6,切分面号为A1。
WPOF,,,l_b                      ! 偏移工作平面,l_b螺栓在盲孔中的旋合长度。
VSBW,ALL                        ! 用工作平面分割体

! (5).  删除体

VDELE,7, , ,1                   ! 删除没有画网格的体,删除7号体。
                                ! 1表示删除体及其以下的图元
VDELE,11, , ,1                  ! 删除没有画网格的体
VDELE,9, , ,1

! (6). 再次切分体

KWPLAN,-1,18,29,14
VSBW,6                          ! 用工作平面分割体

??以上命令运行后将得到如下几何体:

在这里插入图片描述

  • CYL4

??使用功能:在工作平面上生成一个圆面或圆柱体。
??使用格式:CYL4,XCENTERYCENTERRAD1THETA1RAD2THETA2,DEPTH

其中:
??XCENTERYCENTER: 圆心在工作平面x和y方向的坐标值;
??RAD1RAD2: 内、外半径;
??THETA1THETA2: 开始、结束角度;
??DEPTH: 柱体高度。

?&emsp在 GUI 中,操作如下:

在这里插入图片描述

Create hollow cylinder

  • KWPLAN

??使用功能:三个关键点定义一个工作平面
??使用格式:KWPLAN,WNKORIGKXAXKPLAN

其中:
??KORIG: Keypoint number defining the origin of the working plane coordinate system;
??KORIG: Keypoint number defining the x-axis orientation;
??KORIG: Keypoint number defining the working plane。
即分别为:工作平面坐标原点,工作平面x轴上的一点,位于工作平面上的任意一点。

2.4.3 划分网格

??(1). 划分面网格 —— 2D 网格

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
! 为线布种子

ALLSEL,ALL
LESIZE,2, , ,2, , , , ,0  
LESIZE,47, , ,2, , , , ,0  
LESIZE,29, , ,2, , , , ,0  
LESIZE,4, , ,2, , , , ,0  
LESIZE,43, , ,2, , , , ,0  
LESIZE,45, , ,2, , , , ,0  
LESIZE,34, , ,2, , , , ,0

LCCAT,2,47   ! 粘接线以方便映射分网
LCCAT,4,29   ! 粘接线以方便映射分网

CSYS,32
LSEL,S,RADIUS,,d_hole/2
LSEL,R,LOC,Z,0  
LESIZE,ALL, , ,4, , , , ,0  ! 螺栓孔半圆轮廓线上划分4个网格

 
LSEL,S,LENGTH,,d_hole/2-de_bolt/2  
LSEL,R,LOC,Z,0
LESIZE,ALL, , ,1, , , , ,0  ! z=0平面内,螺栓孔与有效直径间的线划分为1个网格


! 为面划分网格

ASEL,R,LOC,Z,0

TYPE,etmax+2
MSHAPE,0,2D   ! 网格形状为2D四边形
MSHKEY,1      ! 映射网格
AMESH,ALL

??以上命令运行后,将得到如下位置处的面网格:

在这里插入图片描述

2D 网格的生成

在这里插入图片描述

单元的材料属性

  • LCCAT

??使用功能:由多条线连成一条连接线,以便于面的映射网格划分。
??使用格式:LCCAT,NL1NL2

其中: NL1NL2 : 要连接的线号,也可以为 ALL、P或元件名。

??对于一个形状不规则且包含的边界线多余4条的面,在没有进行任何处理之前,是不能采用映射网格划分的。

  • LSEL

??使用功能:选择一组线的子集。
??使用格式:LSEL,TypeItemCompVMINVMAXVINCKSWP

其中: Type: Label identifying the type of select / 选择类型的有效标签,它的值有:S、R、A、U、ALL、NONE、INVE、STAT,具体解释如下:

Type Description Type Description
S Select a new set (default). ALL Restore the full set.
R Reselect a set from the current set. NONE Unselect the full set.
A Additionally select a set and extend the current set. INVE Invert the current set (selected becomes unselected and vice versa).
U Unselect a set from the current set. STAT Display the current select status.

在这里插入图片描述

ANSYS 中基本选择功能图示

  • VATT

??使用功能:给所选的体设置划分网格单元属性。/ Associates element attributes with the selected, unmeshed volumes.
??使用格式:VATT,MATREALTYPEESYSSECNUM

其中: MATREALTYPEESYSSECNUM,具体解释如下:

??Material number, real constant set number, type number, coordinate system number, and section number to be associated with selected, unmeshed volumes.

??(2). 体网格的生成 —— 3D 网格

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
! 布种子

LSEL,S,LENGTH,,l_b  ! 按线的长度选择线
! LSEL,U,LOC,Z,0               ! 上一选择集中去除指定位置的线(显然此处是多余的)
! LSEL,U,LOC,Z,l_b
LESIZE,ALL, , ,2, , , , ,0

! 为线划设置分网原则
LSEL,S,LENGTH,,t_flange-l_b           ! 按线长选择直线
LESIZE,ALL, , ,6, , , , ,0

! 为线划设置分网原则  
LSEL,,LENGTH,,,d_hole/2-de_bolt/2
LESIZE,ALL, , ,1, , , , ,0  

! 指定随后生成单元的类型及材料属性

VSEL, , , ,1,2
VATT,etmax+3, ,etmax+3,0  ! 为旋合段螺栓赋予材料属性

TYPE,etmax+1
MAT,etmax+1   ! 赋球墨铸铁材料

! 生成单元即执行分网

ALLSEL,ALL
VSWEEP,ALL  ! 生成网格

! 合并节点及压缩编号

ALLSEL,BELOW,VOLU                        ! 选择体及其附属几何(如:生成体所必须的点、线、面等)
NUMMRG,ALL                               ! 对实体进行合并
NUMCMP,ALL                               ! 压缩所定义项的编号

CM,ele_Bottomf_BH_Outside,ELEM                    ! 由所选单元生成一个元件

! 文件输出

ALLSEL,ALL                                        ! 选择所有实体
CDWRITE,db,'Bottomf_BH_Outside','cdb',,'',''      ! 输出文件,文件名:Bottomf_BH_Outside,格式:.cdb。

EPLOT                                             ! Plot Elements

??以上命令运行后,将得到整个下法兰外圈环带单元,如下:

在这里插入图片描述

2.5 创建下法兰内圈通孔环带单元

??按照本文 2.1 节步骤,可创建下法兰内圈通孔环带单元。

2.6 创建下法兰环带单元 (内圈+外圈 绑定)

??通过在下法兰外圈环带和内圈环带间建立约束,形成整体的下法兰环带单元,与本文 2.3 节类似,所创建的下法兰单元如下所示:

在这里插入图片描述

2.7 创建外圈盲孔螺栓连接

2.7.1 创建垫片单元

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
FINISH                             ! 退出当前处理器
/CLEAR,ALL                         ! 清除所有
/prep7                             ! 进人前处理器
*AFUN,DEG                          ! 指定角度单位为度
PARRES,NEW,'parameters','txt',' '  ! 从文件para1.txt中读入参数

CSYS,0
LOCAL,11,0,L,0,0,0,0,-90  ! 建立以(L,0,0)为原点,局部z轴与总体x轴平行,坐标系号为11的局部直角坐标。
CLOCAL,12,1,0,0,0,0,0,0   ! 以局部坐标系11为参照,建立局部柱坐标系12。
CSYS,12                   ! 激活局部柱坐标系12

 ! (1) 螺栓单元及材料指定

*GET,etmax,ETYP,0,NUM,MAX
ET,etmax+1,185  
KEYOPT,etmax+1,2,2  
MPTEMP,1,0  
MPDATA,EX,etmax+1,,2e5
MPDATA,PRXY,etmax+1,,0.28
MPDATA,DENS,etmax+1,,7.85e-9  

 ! (2) 垫片单元及材料指定

ET,etmax+2,185  
KEYOPT,etmax+2,2,2  
MPTEMP,1,0  
MPDATA,EX,etmax+2,,2e5
MPDATA,PRXY,etmax+2,,0.28
MPDATA,DENS,etmax+2,,7.85e-9  

! ET,20,MESH200            ! 单元类型号为20
! KEYOPT,20,1,6            ! 截面三角形单元4或者6

 ! (3) 创建垫片

! vsel,none

*GET,kvmax,VOLU,,NUM,MAX, , , ,
CLOCAL,21,0,0,0,-t_flange-t_washer     ! 创建局部直角坐标系21

WPCSYS,-1, ! 定义工作平面

CYL4,Ds/2,0,d_washer_inside/2,0,d_washer_outside/2,90,t_washer    
CYL4,Ds/2,0,d_washer_inside/2,90,d_washer_outside/2,180,t_washer    
CYL4,Ds/2,0,d_washer_inside/2,180,d_washer_outside/2,270,t_washer    
CYL4,Ds/2,0,d_washer_inside/2,270,d_washer_outside/2,360,t_washer

! 四个半圆柱体拼成一个垫片,方便网格划分

TYPE,etmax+2
MAT,etmax+2

VGLUE,ALL ! 提粘接

 ! (4) 划分网格

LSEL,S,LENGTH,,t_washer      ! 选择线
LESIZE,ALL, , ,3, , , , ,0   ! 分三段

LSEL,S,LENGTH,,0,d_washer_outside/2-d_washer_inside/2
LESIZE,ALL, , ,3, , , , ,0

VSWEEP,ALL    ! 选择线

CSYS,12
VGEN,num,ALL, , , 0, 360/num, 0, ,0  ! 为划分网格创建的连接线不能复制和阵列

CM,washer_outside,VOLU
ALLSEL,BELOW,VOLU
NUMMRG,NODE
CM,e_washer_outside,ELEM
ALLSEL
NUMCMP,ALL

??由此创建的单个垫片网格如下所示:

在这里插入图片描述

  • VGLUE

??使用功能:体粘接。

??使用格式:VGLUE,NV1,NV2,NV3,NV4,NV5,NV6,NV7,NV8,NV9

其中: NV1NV9 : 要粘接的体号,也可以为 ALL、P或元件名。

??合并各体的公共点,公共边,实际上就是布尔操作,合并体,保留边界线,即 Abaqus 中 Merge >> Retain

在这里插入图片描述

  • VSWEEP

??使用功能:利用与邻近的面单元采用扫略方式对体进行网格划分。

??使用格式:VSWEEP,VNUM,SRCA,TRGA,LSMO

其中:
??VNUM: 体的编号,也可以为 ALL、P或元件名;
??SRCA: 源面的编号;
??TRGA: 目标面的编号。

2.7.2 创建螺栓杆单元

??螺栓杆采用梁单元建模,在ANSYS中,梁单元的截面显示:

??PlotCtrls >> Style >> Size and Shape >> 勾选 Display of element

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
 ! (1) 创建螺栓单元
 
CSYS,0
LOCAL,11,0,L,0,0,0,0,-90  ! 建立以(L,0,0)为原点,局部z轴与总体x轴平行,坐标系号为11的局部直角坐标。
CLOCAL,12,1,0,0,0,0,0,0   ! 以局部坐标系11为参照,建立局部柱坐标系12。
CSYS,12                   ! 激活局部柱坐标系12

*GET,etmax,ETYP,,NUM,MAX
*GET,lmax,LINE,,NUM,MAX
*GET,kpmax,KP,,NUM,MAX
*GET,smax,SECP,0,NUM,MAX

ET,etmax+1,BEAM189  ! 定义梁单元
 
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,etmax+1,,2e5
MPDATA,PRXY,etmax+1,,0.3

SECTYPE,smax+1,BEAM,CSOLID, ,0
SECOFFSET,CENT
SECDATA,de_bolt/2,10,3,            !应力截面直径,圆周方向10份,半径方向3份。

SECTYPE,smax+2,BEAM,CSOLID, ,0
SECOFFSET,CENT
SECDATA,de_bolt/2,10,3,            !可设置螺栓缩进

CSYS,12  
CLOCAL,21,0,0,0,0

zmin=t_flange-l_b                               ! 梁单元最下端点
zmid=(-t_flange-t_washer+l_b)/2                 ! 梁单元预紧位置点

K, ,Ds/2,0,zmin,
K, ,Ds/2,0,zmin-5,
K, ,Ds/2,0,zmid,
K, ,Ds/2,0,-t_flange-t_washer+5,
K, ,Ds/2,0,-t_flange-t_washer,

LSTR,kpmax+1,kpmax+2  

! 创建直线 Defines a straight line irrespective of the active coordinate system.

*REPEAT,4,1,1

LESIZE,lmax+1, , ,2,1, , , ,0     ! 单元数量2
LESIZE,lmax+2, , ,6,1, , , ,0     ! 单元数量4
LESIZE,lmax+3, , ,3,1, , , ,0     ! 单元数量4
LESIZE,lmax+4, , ,2,1, , , ,0     ! 单元数量2

TYPE,etmax+1
MAT,etmax+1
SECNUM,smax+1
LMESH,lmax+1
LMESH,lmax+2

TYPE,etmax+1
MAT,etmax+1
SECNUM,smax+2
LMESH,lmax+3
LMESH,lmax+4

! 显示梁截面
 
/ESHAPE,1    
/REPLOT

LSEL,S,,,lmax+1,lmax+4
CSYS,12  
LGEN,num,ALL, , , ,360/num, , ,0

CM,line_shanks_Outside,LINE
CDREAD,DB,'BottomFlange','cdb',,'',''

2.7.3 创建螺栓间接触

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
! 创建螺栓间的接触(与旋合螺栓的接触)

CSYS,0
LOCAL,11,0,L,0,0,0,0,-90
CLOCAL,12,1,0,0,0,0,0,0

/com, 这是备注

CSYS,12

*GET,etmax,ETYP, ,NUM,MAX  
*GET,rmax,REAL, ,NUM,MAX
 
ET,etmax+1,TARGE170  ! 8 nodes
KEYOPT,etmax+1,2,1   ! 8 nodes  Specified by user

ET,etmax+2,CONTA174
KEYOPT,etmax+2,2,2
KEYOPT,etmax+2,4,1
KEYOPT,etmax+2,12,5

! Multipoint constraint (MPC)
! On nodal point - normal from contact surface
! Bonded (always)


*do,i,1,num

! 创建目标单元

CLOCAL,201,1,Ds/2,360/num*(i-1),0  ! 创建局部柱坐标系201

REAL,rmax+i                        ! 定义实常数并激活
TYPE,etmax+1                       ! 目标单元

TSHAP,PILO                         ! 目标单元由各节点张成,以下为选节点操作

CMSEL,S,line_shanks_Outside        ! 选择螺栓杆轴线
ALLSEL,BELOW,LINE                  ! 选择其下附属节点
E,NODE(0,0,zmin)                   ! 在最下节点处生成单元


ALLSEL,ALL


! 创建接触单元

REAL,rmax+i                        ! 激活实常数
TYPE,etmax+2                       ! 接触单元


CMSEL,S,ele_BottomBolts            ! 选择旋合段螺栓单元
ALLSEL,BELOW,ELEM                  ! 选择单元下属节点

NSEL,R,EXT                         ! 选择表面节点
NSEL,R,LOC,x,0,de_bolt/2+0.1      
NSEL,R,LOC,z,zmin-0.5,zmin+0.5     ! 选择螺栓单元上表面节点
ESURF, ,TOP,                       ! 节点张成单元


ALLSEL,ALL
CSYS,12

*enddo

2.7.4 创建螺栓垫片间接触

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
! 创建螺栓与垫片间的接触(垫片表面与螺栓头间的接触)

*GET,etmax,ETYP, ,NUM,MAX  
*GET,rmax,REAL, ,NUM,MAX
 
ET,etmax+1,TARGE170
KEYOPT,etmax+1,2,1

ET,etmax+2,CONTA174
KEYOPT,etmax+2,2,2
KEYOPT,etmax+2,4,1
KEYOPT,etmax+2,12,5
CSYS,12

*do,i,1,num

REAL,rlmax+i                        !实常数
CLOCAL,202,1,Ds/2,360/num*(i-1),0

TYPE,etmax+1                        !目标单元
TSHAPE,PILO

CMSEL,S,line_shanks_Outside    
ALLSEL,BELOW,LINE
E,NODE(0,0,-t_flange-t_washer)
ALLSEL,ALL

TYPE,etmax+2                        !接触单元
CMSEL,S,ele_washer_outside
ALLSEL,BELOW,ELEM
NSEL,R,EXT
NSEL,R,LOC,x,d_washer_inside/2-0.1,d_washer_outside/2+0.1
NSEL,R,LOC,z,-t_flange-t_washer-0.1,-t_flange-t_washer+0.1

ESLN,S
ESURF, ,TOP,

ALLSEL,ALL
CSYS,12

*enddo

2.7.2.5 创建预紧单元

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
/com,(7) 建立预紧单元   !!!

CLOCAL,21,1,0,0,0

*do,i,1,num

CMSEL,S,line_shanks_Outside
ALLSEL,BELOW,LINE


NSEL,R,LOC,Y,(360/num)*i-(360/num/2+1),(360/num)*i+(360/num/2+1)   ! 选中螺栓上的全部节点
ESLN,,1                                                            ! 选中螺栓单元
ESEL,R,ENAME,,189      

PSMESH,1000+i, , ,ALL, , 11,Z, ,node(0,0,zmid),, , ,               !注意这个命令,创建预紧单元
 
*enddo

ALLSEL,ALL
CDWRITE,db,'All_Bolts_Outside','cdb',,'',''    ! 只含全部下法兰与外圈螺栓
  • ESURF

??使用功能:在存在已选择单元的自由表面上生成重叠单元。/ Generates elements overlaid on the free faces of selected nodes.

??使用格式:ESURF,XNODE,Tlab,Shape

其中:
??Tlab: 用来在法线的合适方向生成目标和接触单元,其值可以为: TOP、BOTTOM 和 REVERS,TOP 表示在梁和壳单元的上表面生成沿法线方向相同的目标单元和接触单元(默认设置)。

  • E

??使用功能:通过节点相连生成一个单元 / Defines an element by node connectivity.

??使用格式:E,I,J,K,L,M,N,O,P

其中:I,J,K,L,M,N,O,P 分别为组成单元的节点编号。

??使用提示:注意节点顺序要符合ANSYS的规定,当前默认的 MAT、TYPE 等属性也将赋给生成的单元。

  • ESLN

??使用功能:选中依附于所选节点上的单元。

??使用格式:ESLN,Type,EKEY,NodeType

其中:Type 确定节点选择类型的标签,其值有S、R、A 和 U;

???EKEY 为节点选择控制键,若为0,选中至少有一个节点被选择的单元(默认);若为1,只有单元中所有节点被选择,该单元才选中。

  • PSMESH

??使用功能:生成预拉伸剖面网格。

??使用格式:PSMESH,SECID,Name,P0,Egroup,NUM,KCN,KDIR,VALUE,NDPLANE,PSTOL,PSTYPE,ECOMP,NCOMP

其中:

??SECID: 唯一的剖面编号,并确定没有被指定给某个剖面;
??Name : 唯一的不超过8个字符的描述性名称;
??P0: 预拉伸节点编号,如果不存在,将生成一个,默认的编号是最大节点编号+1;
??Egroup,NUM : 将要处理的单元组,若 Egroup 为ALL,表示对所有选择的单元进行操作,并忽略 NUM
??KCN :对分离表面及其法向的坐标系参考编号;
??KDIR :在 KCN 坐标系统中与分离面相垂直的方向(x、y 或 z);
??VALUE :沿 KDIR 轴线的点,用于分离表面的定位;
??NDPLANE :用来确定分离表面且已存在的节点编号;

2.8 创建内圈通孔螺栓连接

??与 2.7 节类似,可创建内圈通孔垫片、螺栓以及螺栓与垫片间的接触。

2.9 创建垫片与法兰间的接触

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
FINISH                             ! 退出当前处理器
/CLEAR,ALL                         ! 清除所有
/prep7                             ! 进人前处理器
*AFUN,DEG                          ! 指定角度单位为度
PARRES,NEW,'parameters','txt',' '  ! 从文件para1.txt中读入参数

 ! ***************************** ---- 1、建法兰与垫片间的接触关系 ---- **************************** !

! (1) 导入 .cdb 文件

! CDREAD,DB,'TowerFlange','cdb',,'',''  ! 特别注意导入次序

CDREAD,DB,'All_Bolts_Outside','cdb',,'',''
CDREAD,DB,'All_Bolts_Inside','cdb',,'',''
CDREAD,DB,'TopFlange','cdb',,'',''

CSYS,0
LOCAL,11,0,L,0,0,0,0,90 ! 建立以(LL,0,0)为原点,局部z轴正向与总体x轴正向相同,坐标系号为15的局部直角坐标。
LOCAL,12,1,L,0,0,0,0,90 ! 建立以(LL,0,0)为原点,局部z轴正向与总体x轴正向相同,坐标系号为16的局部柱坐标。
CLOCAL,13,0
CLOCAL,14,1

! (2) 建立法兰与垫片间接触(标准接触)

! 法兰186单元、垫片是185单元 建接触时可以用带中间节点的高阶单元?。

! 定义单元类型

CSYS,12

*GET,rlmax,RCON,0,NUM,MAX  ! 获取最大实常数
*GET,etmax,ETYP, ,NUM,MAX  ! 获取最大单元类型号

REAL,rlmax+1
ET,etmax+1,170   ! Targe170 8 nodes

ET,etmax+2,174   ! Conta174 8 nodes
KEYOPT,etmax+2,4,0
KEYOPT,etmax+2,5,4  
KEYOPT,etmax+2,9,2
KEYOPT,etmax+2,10,2

! On Gauss point
! Auto ICONT
! Include both initial geometrical penetration or gap and offset, but with ramped effects

! Each iteration based on the current mean stress of underlying elements. The actual elastic slip never
! exceeds the maximum allowable limit (SLTO) during the entire solution.

! 定义材性

MP,MU,etmax+1,0.3 ! 指定摩擦系数为0.3
MAT,etmax+1

! 创建接触单元

TYPE,etmax+2

CSYS,12
CMSEL,S,ele_washer_Inside
CMSEL,A,ele_washer_Outside

NSLE
NSEL,R,EXT  
NSEL,U,LOC,Z,t_flange+0.5,t_flange+t_washer+0.5
NSEL,U,LOC,Z,-(t_flange+t_washer+0.5),-(t_flange+0.5)  ! 选节点

ESLN,S,0
ESURF   ! 创建单元

! 创建目标单元

TYPE,etmax+1

CSYS,12
CMSEL,S,ele_BottomF_TH_Inside
CMSEL,A,ele_TopF_TH_Outside
CMSEL,A,ele_TopF_TH_Inside

NSLE
NSEL,R,EXT
NSEL,U,LOC,Z,-(t_flange-0.5),t_flange-0.5 ! 选节点

ESLN,S,0
ESURF  ! 创建单元

ALLSEL,ALL
CDWRITE,db,'FlangeWithBolts','cdb',,'',''

??由此创建的接触单元如下所示:

在这里插入图片描述

3. 装配与约束创建

4. 施加边界条件

二. 加载与求解
三. 后处理