24小时热门版块排行榜

返回列表

难过先生

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 3
散金: 55
红花: 2
帖子: 84
在线: 45.8小时
虫号: 4250802
注册: 2015-11-27
性别: GG
专业: 防灾工程

[求助] ansys中螺栓连接的接触分析已有1人参与

本人建了一个两块板之间通过螺栓连接的有限元模型，做接触分析。模型建好之后不能计算，在计算的时候总是闪退。请各位亲帮忙看看是怎么回事，谢谢。命令流如下：!螺栓连接
FINISH
/CLEAR
/PREP7
ET,1,SOLID95  !单元号1
MP,EX,1,206000  !材料号1
MP,PRXY,1,0.3
MP,MU,1,0.44 !抗滑移摩擦系数
TB,BKIN,1  !随动强化准则
TBDATA,,420,0,,,,
MP,EX,2,206000  !材料号2
MP,PRXY,2,0.3
MP,MU,2,0.44
EY=960/206000$FY=960  !定义10.9级螺栓屈服强度及屈服应变
TB,MISO,2,,2 !随动强化准则
TBPT,,EY,FY
TBPT,,0.011,1040

ET,2,TARGE170    !目标单元
ET,3,CONTA174       !接触单元

KEYOPT,3,2,0 !增广的拉格朗日法
KEYOPT,3,5,1 !闭合间隙
KEYOPT,3,9,1 !忽略几何穿透/间隙和CNOF
KEYOPT,3,12,0  !标准接触
!!!!!!!!!!!!定义参数,单位为MM，MP,N
L=500    !板长500mm
B=200       !板宽200mm
T=40       !板厚40
R=10       !螺栓孔及螺杆半径
D=420       !螺栓连接中心到板边的距离
L1=80       !螺栓中心到板宽方向的距离
L2=60       !螺栓中心到板长方向的距离
LRD=80    !螺栓中心之间的距离
RLM=20    !螺帽的半径
DLM=20    !螺帽的厚度
YJL=500
!!!!!!!!!!!!建模开始!!!!!!!!!!!!!!
!!!!!!!!!!!!板1的几何模型!!!!!!

*GET,KSTART,KP,0,NUM,MAX
K,KSTART+1,0,0,0
K,KSTART+2,500,0,0
K,KSTART+3,500,200,0
K,KSTART+4,0,200,0
*GET,LSTART,LINE,0,NUM,MAX
NUMSTR,LINE,LSTART       !按顺序生成连线
L,KSTART+1,KSTART+2
L,KSTART+2,KSTART+3
L,KSTART+3,KSTART+4
L,KSTART+4,KSTART+1
AL,LSTART+1,LSTART+2,LSTART+3,LSTART+4  !由封闭线段生成面

VEXT,ALL,,,,,T ! 在Z方向延伸，将面延伸成体
CM,BAN1,VOLU !定义板1为元件
VSEL,NONE
ASEL,NONE

!!!!!!!!!!!!板2的几何模型!!!!!!
WPOFFS,L-2*L1,0,T    !!!!!!!!!!!!!移动工作平面
CSYS,4          !以工作平面为坐标系
*GET,KSTART,KP,0,NUM,MAX
K,KSTART+1,0,0,0
K,KSTART+2,500,0,0
K,KSTART+3,500,200,0
K,KSTART+4,0,200,0
*GET,LSTART,LINE,0,NUM,MAX
NUMSTR,LINE,LSTART       !按顺序生成连线
L,KSTART+1,KSTART+2
L,KSTART+2,KSTART+3
L,KSTART+3,KSTART+4
L,KSTART+4,KSTART+1
AL,LSTART+1,LSTART+2,LSTART+3,LSTART+4  !由封闭线段生成面
VEXT,ALL,,,,,T ! 在Z方向延伸，将面延伸成体
CM,BAN2,VOLU !定义板2为元件
VSEL,NONE
ASEL,NONE
!!!!!!!!!!!!!!!螺杆建模!!!!!!!!!!!!!!
WPCSYS,-1,0 !!!!!!!工作平面回到原坐标系下
WPOFFS,0,B/2,0    !移动工作平面
WPOFFS,D,LRD/2,0    !移动工作平面
CSYS,4
!!!!!!!!!!!!!!!!上面螺栓杆建模开始!!!!!!!!!!!
CYLIND,R,,-DLM,0,0,360
CYLIND,R,,0,T,0,360
CYLIND,R,,T,2*T,0,360
CYLIND,R,,2*T,2*T+DLM,0,360
!!!!!!!!!!!!!!!!下面螺栓杆建模开始!!!!!!!!!!!
WPCSYS,-1,0
WPOFFS,0,B/2,0    !移动工作平面
WPOFFS,D,-LRD/2,0    !移动工作平面
CSYS,4
CYLIND,R,,-DLM,0,0,360
CYLIND,R,,0,T,0,360
CYLIND,R,,T,2*T,0,360
CYLIND,R,,2*T,2*T+DLM,0,360
VADD,ALL                !这里的体与体的相加是有相互接触的地方相加，所以只有螺杆的几部分加成了一个整体。
CM,LG,VOLU                      !!!!!!!!!!!!!定义螺杆为元件LG
CMSEL,S,BAN1,VOLU
CMSEL,A,LG,VOLU
VSBV,BAN1,LG,SEPO,DELETE,KEEP    !板1减螺杆（板中与螺栓的重合部分被删除，对应部分在螺栓中不删除。）
CMSEL,U,LG,VOLU                   !去掉螺杆
CM,BAN1,VOLU                      !!!!!!!!!!!!!定义带螺栓孔板为元件BAN1
VSEL,NONE
ASEL,NONE
CMSEL,S,BAN2,VOLU
CMSEL,A,LG,VOLU
VSBV,BAN2,LG,SEPO,DELETE,KEEP    !板2减螺杆（板中与螺栓的重合部分被删除，对应部分在螺栓中不删除。）
CMSEL,U,LG,VOLU                   !去掉螺杆
CM,BAN2,VOLU                      !!!!!!!!!!!!!定义带螺栓孔板为元件BAN1
VSEL,NONE
ASEL,NONE
!!!!!!!!!!!!!!!!螺帽部分建模!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!上面螺杆对应螺帽建模!!!!!!!!!!!!!!!!!
WPCSYS,-1,0
WPOFFS,0,B/2,0    !移动工作平面
WPOFFS,D,LRD/2,0    !移动工作平面
CSYS,4
CYLIND,RLM,,-DLM,0,0,360 !板低螺帽
CYLIND,RLM,,2*T,2*T+DLM,0,360 !板顶螺帽
!!!!!!!!!!!!!!!!下面螺杆对应螺帽建模!!!!!!!!!!!!!!!!!
WPCSYS,-1,0
WPOFFS,0,B/2,0    !移动工作平面
WPOFFS,D,-LRD/2,0    !移动工作平面
CSYS,4
CYLIND,RLM,,-DLM,0,0,360 !板低螺帽
CYLIND,RLM,,2*T,2*T+DLM,0,360 !板顶螺帽
CM,LM,VOLU
CMSEL,S,LM,VOLU
CMSEL,A,LG,VOLU
VSBV,LM,LG,,DELETE,KEEP    !!!在螺帽上挖螺栓孔
VGLUE,ALL                   !!!!进行粘结运算，母体个数不变，但是边界是共享的
CMSEL,U,LG,VOLU          !去掉螺杆部分
CM,LM,VOLU          !!!!!!!!!!!!重新定义挖了孔的螺帽为元件LM
VSEL,NONE
ASEL,NONE
!!!!!!!!!!!!!!!!!!!!网格划分，建立有限元模型!!!!!!!!!!!!!!!!!!!!!（可以尝试将节点部分剖分后再画网格）
ALLSEL,ALL
CMSEL,S,LG,VOLU
CMSEL,A,LM,VOLU
VATT,2,,1,0             !定义属性
!LSEL,S,LENGTH,,15.71
!LESIZE,ALL,,,3
!LSEL,A,LENGTH,,20
!LESIZE,ALL,,,4
!LSEL,A,LENGTH,,40
!LESIZE,ALL,,,8
!LSEL,A,LENGTH,,31.42
!LESIZE,ALL,,,6
!VSWEEP,ALL
!ALLSEL,ALL
ESIZE,5
MSHAPE,0,2D
MSHKEY,0
VSWEEP,ALL

CMSEL,S,BAN1,VOLU
WPCSYS,-1,0
WPOFFS,L-2*L1
WPROTA,,,90
CSYS,4
VSBW,ALL    !利用工作面切分图素（利用XY平面切分，切而不分）
CM,BAN1,VOLU
CMSEL,S,BAN1,VOLU
VATT,1,,1,0
ESIZE,15
MSHAPE,0,2D
MSHKEY,0
VSWEEP,ALL

CMSEL,S,BAN2,VOLU
!CMPLOT
WPCSYS,-1,0
WPOFFS,L,0,T
WPROTA,,,90
CSYS,4
VSBW,ALL    !利用工作面切分图素（利用XY平面切分，切而不分）
CM,BAN2,VOLU
CMSEL,S,BAN2,VOLU
VATT,1,,1,0
ESIZE,15
MSHAPE,0,2D
MSHKEY,0
VSWEEP,ALL

!!!!!!!!!!!!!!!!!!!!!!!!!!建立接触!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

!!!!!!!!!!!!!!!!!!BAN1下表面与对应螺母上表明的接触关系!!!!!!!!!!!!!!!!!!
ALLSEL,ALL
WPCSYS,-1,0
CMSEL,S,BAN1,VOLU
ASLV,S
ASEL,R,LOC,Z,0    !BAN1底面
ASEL,U,,,11
NSLA,S,1
R,5
REAL,5
MAT,1
TYPE,2
ESURF,ALL !!!!!!!!!目标面生成

CMSEL,S,LM,VOLU
ASLV,S !选择与所选体相关的面
ASEL,R,LOC,Z,0  !从当前子集中再选择一组面，形成新的当前子集  板下螺帽的上表明
NSLA,S,1
R,5
REAL,5
MAT,2
TYPE,3
ESURF,ALL !!!!!!接触面生成

!!!!!!!!!!!!!!!!!!BAN2上表面与对应螺母下表明的接触关系!!!!!!!!!!!!!!!!!!

CMSEL,S,BAN2,VOLU
ASLV,S
ASEL,R,LOC,Z,2*T    !BAN1底面
ASEL,U,,,24
NSLA,S,1
R,6
REAL,6
MAT,1
TYPE,2
ESURF,ALL !!!!!!!!!目标面生成

CMSEL,S,LM,VOLU
ASLV,S !选择与所选体相关的面
ASEL,R,LOC,Z,2*T  !从当前子集中再选择一组面，形成新的当前子集  板下螺帽的上表明
NSLA,S,1
R,6
REAL,6
MAT,2
TYPE,3
ESURF,ALL !!!!!!接触面生成

!!!!!!!!!!!!!!!!!!BAN1侧面与螺杆侧面的接触关系!!!!!!!!!!!!!!!!!!

CMSEL,S,BAN1,VOLU
ASLV,S
ASEL,R,LOC,Z,0,T
ASEL,U,LOC,X,0
ASEL,U,LOC,X,L
ASEL,U,LOC,X,L-2*L1
ASEL,U,LOC,Y,0
ASEL,U,LOC,Y,B
ASEL,U,LOC,Z,0
ASEL,U,LOC,Z,T
NSLA,S,1
R,7
REAL,7
MAT,1
TYPE,2
ESURF,ALL !!!!!!目标面生成

CMSEL,S,LG,VOLU
ASLV,S
ASEL,R,LOC,Z,0,T
NSLA,S,1
R,7
REAL,7
MAT,2
TYPE,3
ESURF,ALL !!!!!!!!!接触面生成

!!!!!!!!!!!!!!!!!!BAN2侧面与螺杆侧面的接触关系!!!!!!!!!!!!!!!!!!

CMSEL,S,BAN2,VOLU
ASLV,S
ASEL,R,LOC,Z,T,2*T
ASEL,U,LOC,X,L-2*L1
ASEL,U,LOC,X,L
ASEL,U,LOC,X,2*L-2*L1
ASEL,U,LOC,Y,0
ASEL,U,LOC,Y,B
ASEL,U,LOC,Z,T
ASEL,U,LOC,Z,2*T
NSLA,S,1
R,8
REAL,8
MAT,1
TYPE,2
ESURF,ALL !!!!!!目标面生成

CMSEL,S,LG,VOLU
ASLV,S
ASEL,R,LOC,Z,T,2*T
NSLA,S,1
R,8
REAL,8
MAT,2
TYPE,3
ESURF,ALL !!!!!!!!!接触面生成

!!!!!!!!!!!!!!!!!!BAN2和BAN1的接触关系!!!!!!!!!!!!!!!!!!
CMSEL,S,BAN1,VOLU
ASLV,S
ASEL,R,LOC,Z,T    !BAN1底面
ASEL,U,,,12
NSLA,S,1
R,9
REAL,9
MAT,1
TYPE,2
ESURF,ALL !!!!!!!!!目标面生成

CMSEL,S,BAN2,VOLU
ASLV,S
ASEL,R,LOC,Z,T    !BAN1底面
ASEL,U,,,18
NSLA,S,1
R,9
REAL,9
MAT,1
TYPE,3
ESURF,ALL !!!!!!!!!接触面生成

allsel,all
VSEL,S,,,11
VSEL,A,,,12
PSMESH,10000,YJM1,,VOLU,11,0,z,60,,,,elems    !创建预紧面并划分网格
PSMESH,100001,YJM2,,VOLU,12,0,z,60,,,,elems    !创建预紧面并划分网格

!!!!!!!!!!!!!!!!!!!!施加边界条件!!!!!!!!!!!!!!!
allsel,all
FINI
/SOLU
ASEL,S,LOC,X,0
NSLA,S,1
D,ALL,ALL
!!!!!!!!!施加预紧荷载!!!!!!!!!!!!!

SLOAD,10000,PL01,,FORCE,500,1,2
SLOAD,100001,PL01,,FORCE,500,1,2
nsubst,25,100,1
SOLVE@wuming524

回复此楼