24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 8814  |  回复: 3

难过先生

新虫 (小有名气)

[求助] 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
回复此楼

» 收录本帖的淘帖专辑推荐

ANSYS命令流

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
难过先生: 金币+40, ★★★★★最佳答案 2017-04-03 15:14:09
在求解之前,你选择了部分节点施加约束,然后直接求解了,这时不对的。
  在求解之前执行命令:
      allsel,all
  然后再求解,看是否还有其它问题。
2楼2017-04-03 08:08:28
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

littlepoint

木虫 (职业作家)

楼上听说是个问题。另外,接触问题不容易收敛,需要平稳、分步建立接触关系。

发自小木虫Android客户端
3楼2017-04-03 12:11:05
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

难过先生

新虫 (小有名气)

引用回帖:
2楼: Originally posted by shgao20 at 2017-04-03 08:08:28
在求解之前,你选择了部分节点施加约束,然后直接求解了,这时不对的。
  在求解之前执行命令:
      allsel,all
  然后再求解,看是否还有其它问题。

谢谢
4楼2017-04-03 15:14:38
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 难过先生 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 262求调剂 +4 励志一定发文章 2026-04-02 5/250 2026-04-02 11:07 by 励志一定发文章
[考研] 化工调剂303分,过四级 +3 栖梧待风 2026-04-02 3/150 2026-04-02 10:29 by 王保杰33
[考研] 材料科学与工程调剂 +18 深V宿舍吧 2026-03-30 19/950 2026-04-02 10:28 by sanrepian
[基金申请] esi高被引论文是不是能对中标有所加分和帮助呢 +4 redcom 2026-04-01 4/200 2026-04-02 10:16 by wqy8563
[考研] 266求调剂 +4 学员97LZgn 2026-04-02 4/200 2026-04-02 09:52 by yulian1987
[考研] 材料调剂 +12 一样YWY 2026-04-01 12/600 2026-04-02 09:15 by olim
[考研] 材料工程322分 +7 哈哈哈吼吼吼哈 2026-04-01 7/350 2026-04-02 07:05 by yjolah
[考研] 318求调剂,计算材料方向 +7 吸喵有害笙命 2026-04-01 8/400 2026-04-01 23:54 by 义文wang
[考研] 085600 295分求调剂 +15 W55j 2026-03-30 18/900 2026-04-01 23:00 by 来看流星雨10
[考研] 材料专硕306英一数二 +7 z1z2z3879 2026-03-31 7/350 2026-04-01 14:50 by ZXlzxl0425
[考研] 330分求调剂 +11 qzenlc 2026-03-29 11/550 2026-04-01 14:32 by chenqifeng666
[考研] 267求调剂 +13 uiybh 2026-03-31 13/650 2026-04-01 10:25 by 探123
[考研] 311(085601)求调剂 +12 liziyeyeye 2026-03-28 13/650 2026-04-01 00:34 by fmesaito
[考研] 求调剂 生物学 377分 +6 zzll03 2026-03-31 6/300 2026-03-31 17:33 by 唐沐儿
[考研] 求调剂 +9 akdhjs 2026-03-31 10/500 2026-03-31 17:14 by zhyzzh
[考研] 福建理工大学材料学院先进合金团队招收考研调剂学生 +3 大华金商都 2026-03-30 4/200 2026-03-31 01:04 by 方英俊602
[考研] 一志愿南开大学0710生物学359求调剂 +5 兔兔兔111223314 2026-03-29 7/350 2026-03-30 18:29 by 兔兔兔111223314
[考研] 071010 323 分求调剂 +3 Baekzhy 2026-03-27 3/150 2026-03-30 14:24 by andresqi
[考研] 求佛 +7 迷人的哈哈 2026-03-28 7/350 2026-03-28 16:47 by 催化大白
[考研] 考研调剂 +9 小蜡新笔 2026-03-26 9/450 2026-03-27 11:10 by 不吃魚的貓
信息提示
请填处理意见