你在这里
使用ABAQUS有限元分析程序求解塑性接触问题。
星期四，2008-04-10 12:41 -lazghabt
大家好:
1.我已经建立了有限元模型二维平面应变使用ABAQUS/标准(CPE4R)元素，
2.我正在用二维模型模拟圆柱体和平板之间的接触。
3.给定的圆柱体材料为无塑性的弹性各向同性材料E=310 GPa
4.该板具有Ramberg-Osgood塑性，杨氏模量E=210 GPa。
5.接触点处的网格是非常精细的元素，特征长度约为1微米。(我的模型单位是(mm,N,sec))。
6.我用赫兹接触理论验证了我的模型，效果很好。
7.我尝试用mindlin理论(压力和剪切)来验证，当我在接触区绘制压力剖面时，我得到了一条表现良好的曲线。然而，接触区剪切表现不佳，表现出较大的波动和噪声响应。
8.我试图修改接触区的网格尺寸(实际上增加了)，但响应仍然是相同的。
注:如果我从模型中去除塑性，它在正常压力和剪切下都能正常工作。
问题:
1.如何才能纠正这种情况?
2.是否有特定的塑料模型我应该使用(我的材料是金属)?
3.接触控制可以用来补救吗?如果是，怎么做?
»
- lazghabt的博客
- 登录或注册发表评manbetx体育论
- 24139年读
评论
塑性接触问题
你好,
在深入讨论细节之前，我想问几个问题:
-你如何应用加载?
-你使用了哪些接触属性?法向和切向?
-我建议使用其他具有各向同性硬化等的塑性模型(只是为了检查)
-最重要的是，你的正应力值和剪应力值有什么不同?可能剪切部分根据载荷类型可以忽略不计?
-此外，如果您期望有限滑动，请确保在ABAQUS中使用有限滑动算法。
这是我目前能想到的几点。如果你能解释一下，我可能会给你更多的信息。
祝一切顺利。
阿米尔
aslam alikum, 2D接触问题
aslam alikum。我需要你的帮助，我的兄弟。我面临着接触问题代码的问题。这个想法是使用ALE在2D中接触问题
(任意拉格朗日-欧拉算法)。如果你能帮我，就给我传个信。我的电子邮箱是mimo-shishi@hotmail.comm
附加的输入文件。
你好阿米尔:
感谢您抽出时间来回复。为了回答你所有的问题，我在这里包括了我正在工作的输入文件，也许这可以对我遇到的困难有所帮助。
**+--------------------------------+ #
**单位| #
**+--------------------------------+ #
**长度:mm (1e - 3m) #
**力:牛顿(1n) #
**应力:MPa #
**角度:degrees #
* * #
**---------------------------------------------------------------------------------------------------------------#
**
*预印本,模型=是的,历史=不,联系= NO
*标题
圆柱/平面接触nowell & hill模型
**
**--------------------------------------+
* *定义领域的种子节点 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**--------------------------------------+
**
*节点,NSET = NCenterS
1, 0。, 0。
*节点,NSET = LINE1S
2, 12.0， -28.6478897565
162, 12.0， -14.3239448783
234,12.0， -7.16197243914
274, 12.0， -3.58098621957
290, 12.0， -2.44461992589
1314,12.0, 2.44461992589
1330,12.0, 3.58098621957
1370,12.0, 7.16197243914
1442, 12.0, 14.3239448783
1602, 12.0, 28.6478897565
**
**--------------------------------------+
* *生成LINE1S的节点 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**--------------------------------------+
**
* NGEN NSET = LINE1S
2, 162年,1
162、234、1
234、274、1
274、290、1
290、1314、1
1314、1330、1
1330、1370、1
1370、1442、1
1442、1602、1
**
**---------------------------------------------+
使用* * *生成么,3和4 NCOPY命令 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**---------------------------------------------+
**
* ncopy, shift，旧set = line1s，新set = line2s，更改number =3202
-0.002、0、0
0, 0, 0, 0, 0, 1, 0。
**
* ncopy, shift，旧set = line1s，新set = line3s，更改number =38424
-0.03、0、0
0, 0, 0, 0, 0, 1, 0。
**
* ncopy, shift，旧set = line1s，新set = line4s，更改number =115272
-0.12、0、0
0, 0, 0, 0, 0, 1, 0。
**
* ncopy, shift，旧set = line1，新set = line5，更改number =281776
-0.48、0、0
0, 0, 0, 0, 0, 1, 0。
**
* ncopy, shift，旧set = line1s，新set = line6s，更改number =666016
-9.0、0、0
0, 0, 0, 0, 0, 1, 0。
**
* ncopy, shift，旧set = line6s，新set = line7s，更改号码=512320
-2.4、0、0
0, 0, 0, 0, 0, 1, 0。
**
**+----------------------------------------------------+
**填充球体的节点。|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+----------------------------------------------------+
**
* NFILL NSET =球体
LINE1S么2,1601年
么LINE3S 11, 3202年
LINE3S LINE4S 12, 6404年
LINE4S LINE5S 13, 12808年
LINE5S LINE6S 15, 25616
LINE6S LINE7S 10, 51232年
**
**+----------------------------------------------------+
* *的地图节点从圆柱坐标范围 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+----------------------------------------------------+
**
* NSET NSET =球体
NCenterS
* NMAP, NSET =球体,TYPE =圆柱形
0。,12.0,0,0,12.0,1。
0。1。0。
1。1。1。
**
**---------------------------------------------+
**生成球体元素。|
**---------------------------------------------+
**
*元素,ELSET = ZONE1S = CPE4类型
1290291、1892、1891
* ELGEN ELSET = ZONE1S
1 1024年1 1 2,1601,1024
**
*元素,ELSET = ZONE2S = CPE4类型
2049, 274276, 3478, 3476
2585、3492、3494、6696、6694
2569、1314、1316、4518、4516
* ELGEN ELSET = ZONE2S
2049、8、2、1,12,3202,528
2585、512、2、1,11日,3202年,528年
2569、8、2、1,12,3202,528
**
*元素,ELSET = ZONE3S = CPE4类型
8385, 234238, 6642, 6638
10099、38698、38702、45106、45102
8659、1330、1334、7738、7734
* ELGEN ELSET = ZONE3S
8385、10、4、1、18、6404、284
10099、264、4、1、6404、284
8659、10、4、1、18、6404、284
**
*元素,ELSET = ZONE4S = CPE4类型
13497, 162170, 12978, 12970
14946、115506、115514、128322、128314
13648、1370、1378、14186、14178
* ELGEN ELSET = ZONE4S
1 13497 9, 8日,22日,12808年,160年
14946、142、8、13、12808、160
1 13648 9, 8日,22日,12808年,160年
**
*元素,ELSET = ZONE5S = CPE4类型
17017、2、18、25634、25618
18127、281938、281954、307570、307554
17107、1442、1458、27074、27058
* ELGEN ELSET = ZONE5S
1 17017年10,16日,26日,25616年,100年
1 18127年,80年,16日,15日,25616年,100年
1 17107年10,16日,26日,25616年,100年
**
*元素,ELSET = ZONE6S = CPE4类型
19617、666018、666050、717282、717250
* ELGEN ELSET = ZONE6S
19617年,50岁,32岁,1、10、51232、50
**
**
**+-------------------------------------------+
为SPHERE定义网格细化MPCs。|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+-------------------------------------------+
**
ZONE1S和ZONE2S之间线性连接的**节点集
**
* NSET NSET = NSET1LnrMPCaZ1Z2S,生成
290290年,3202年
* NSET NSET = NSET1LnrMPCbZ1Z2S,生成
3492、3492、3202
* NSET NSET = NSET1LnrMPCpZ1Z2S,生成
1891、1891、3202
** 1 2
* NSET NSET = NSET2LnrMPCaZ1Z2S,生成
3492、4514、2
* NSET NSET = NSET2LnrMPCbZ1Z2S,生成
3494、4516、2
* NSET NSET = NSET2LnrMPCpZ1Z2S,生成
3493、4515、2
** 1 2
* NSET NSET = NSET3LnrMPCaZ1Z2S,生成
1314、1314、3202
* NSET NSET = NSET3LnrMPCbZ1Z2S,生成
4516、4516、3202
* NSET NSET = NSET3LnrMPCpZ1Z2S,生成
2915、2915、3202
**
ZONE2S和ZONE3S之间线性连接的**节点集
**
* NSET NSET = NSET1LnrMPCaZ2Z3S,生成
274、32294、6404
* NSET NSET = NSET1LnrMPCbZ2Z3S,生成
6678、38698、6404
* NSET NSET = NSET1LnrMPCpZ2Z3S,生成
3476、35496、6404
** 2 3
* NSET NSET = NSET2LnrMPCaZ2Z3S,生成
38698、39750、4
* NSET NSET = NSET2LnrMPCbZ2Z3S,生成
38702、39754、4
* NSET NSET = NSET2LnrMPCpZ2Z3S,生成
38700、39752、4
** 2 3
* NSET NSET = NSET3LnrMPCaZ2Z3S,生成
1330、33350、6404
* NSET NSET = NSET3LnrMPCbZ2Z3S,生成
7734、39754、6404
* NSET NSET = NSET3LnrMPCpZ2Z3S,生成
4532、36552、6404
**
ZONE3S和ZONE4S之间线性连接的**节点集
**
* NSET NSET = NSET1LnrMPCaZ3Z4S,生成
234、102698、12808
* NSET NSET = NSET1LnrMPCbZ3Z4S,生成
13042、115506、12808
* NSET NSET = NSET1LnrMPCpZ3Z4S,生成
6638、109102、12808
** 3 4
* NSET NSET = NSET2LnrMPCaZ3Z4S,生成
115506、116634、8
* NSET NSET = NSET2LnrMPCbZ3Z4S,生成
115514、116642、8
* NSET NSET = NSET2LnrMPCpZ3Z4S,生成
115510、116638、8
** 3 4
* NSET NSET = NSET3LnrMPCaZ3Z4S,生成
1370、103834、12808
* NSET NSET = NSET3LnrMPCbZ3Z4S,生成
14178、116642、12808
* NSET NSET = NSET3LnrMPCpZ3Z4S,生成
7774、110238、12808
**
ZONE4S和ZONE5S之间线性连接的**节点集
**
* NSET NSET = NSET1LnrMPCaZ4Z5S,生成
162、256322、25616
* NSET NSET = NSET1LnrMPCbZ4Z5S,生成
25778、281938、25616
* NSET NSET = NSET1LnrMPCpZ4Z5S,生成
12970、269130、25616
** 4 5
* NSET NSET = NSET2LnrMPCaZ4Z5S,生成
281938年,283202年,16岁
* NSET NSET = NSET2LnrMPCbZ4Z5S,生成
281954年,283218年,16岁
* NSET NSET = NSET2LnrMPCpZ4Z5S,生成
281946年,283210年,16岁
** 4 5
* NSET NSET = NSET3LnrMPCaZ4Z5S,生成
1442、257602、25616
* NSET NSET = NSET3LnrMPCbZ4Z5S,生成
27058、283218、25616
* NSET NSET = NSET3LnrMPCpZ4Z5S,生成
14250、270410、25616
**
ZONE5S和ZONE6S之间线性连接的**节点集
**
* NSET NSET = NSET1LnrMPCaZ5Z6S,生成
666018年,667586年,32岁
* NSET NSET = NSET1LnrMPCbZ5Z6S,生成
666050年,667618年,32岁
* NSET NSET = NSET1LnrMPCpZ5Z6S,生成
666034年,667602年,32岁
**
* MPC
线性、NSET1LnrMPCpZ1Z2S NSET1LnrMPCaZ1Z2S NSET1LnrMPCbZ1Z2S
线性、NSET2LnrMPCpZ1Z2S NSET2LnrMPCaZ1Z2S NSET2LnrMPCbZ1Z2S
线性、NSET3LnrMPCpZ1Z2S NSET3LnrMPCaZ1Z2S NSET3LnrMPCbZ1Z2S
**
线性、NSET1LnrMPCpZ2Z3S NSET1LnrMPCaZ2Z3S NSET1LnrMPCbZ2Z3S
线性、NSET2LnrMPCpZ2Z3S NSET2LnrMPCaZ2Z3S NSET2LnrMPCbZ2Z3S
线性、NSET3LnrMPCpZ2Z3S NSET3LnrMPCaZ2Z3S NSET3LnrMPCbZ2Z3S
**
线性、NSET1LnrMPCpZ3Z4S NSET1LnrMPCaZ3Z4S NSET1LnrMPCbZ3Z4S
线性、NSET2LnrMPCpZ3Z4S NSET2LnrMPCaZ3Z4S NSET2LnrMPCbZ3Z4S
线性、NSET3LnrMPCpZ3Z4S NSET3LnrMPCaZ3Z4S NSET3LnrMPCbZ3Z4S
**
线性、NSET1LnrMPCpZ4Z5S NSET1LnrMPCaZ4Z5S NSET1LnrMPCbZ4Z5S
线性、NSET2LnrMPCpZ4Z5S NSET2LnrMPCaZ4Z5S NSET2LnrMPCbZ4Z5S
线性、NSET3LnrMPCpZ4Z5S NSET3LnrMPCaZ4Z5S NSET3LnrMPCbZ4Z5S
**
线性、NSET1LnrMPCpZ5Z6S NSET1LnrMPCaZ5Z6S NSET1LnrMPCbZ5Z6S
**
**--------------------------------------------------------------------------+
* *接触表面的定义范围 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**--------------------------------------------------------------------------+
**
* ELSET ELSET = ESurfSphere,生成
1, 1024, 1
2049、2056、1
8385、8394、1
13497、13505、1
17017、17026、1
2569、2576、1
8659、8668、1
13648、13656、1
17107、17116、1
**
*表面,名称= SurfSphere类型=元素
ESurfSphere, S1
**
**+--------------------------------------------+
* *定义材料属性的球体 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+--------------------------------------------+
**
*材料,NAME =球体
*弹性、类型=各向同性
310000.0, 0.27
*密度
3900.0
**
*固体部分,ELSET = Zone1S材料=球体
1.0
*固体部分,ELSET = Zone2S材料=球体
1.0
*固体部分,ELSET = Zone3S材料=球体
1.0
*固体部分,ELSET = Zone4S材料=球体
1.0
*固体部分,ELSET = Zone5S材料=球体
1.0
*固体部分,ELSET = Zone6S材料=球体
1.0
*刚体，REF NODE=NCenterS,POSITION=INPUT,ELSET=Zone6S,DENSITY=2700。
**
**~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**~~~~~~~~~~~~~~~~~~~~~~~~~~~~~########################################~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**~~~~~~~~~~~~~~~~~~~~~~~~~~~~~# 有限元分析的输入命令 #~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**~~~~~~~~~~~~~~~~~~~~~~~~~~~~~########################################~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**
**--------------------------------------+
* *定义种子节点板 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**--------------------------------------+
**
**
*节点,NSET = LINE1P
1181539， -8.0, 0.0
1181859， -6.0, 0.0
1182019， -3.0, 0.0
1182091， -1.5, 0.0
1182131， -0.75, 0.0
1182147， -0.512, 0.0
1183171, 0.512, 0.0
1183187, 0.75, 0.0
1183227,1.5, 0.0
1183299, 3.0, 0.0
1183459, 6.0, 0.0
1183779, 8.0, 0.0
**
**----------------------------------------+
* *生成LINE1P的节点板 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**----------------------------------------+
**
* NGEN NSET = LINE1P
1181539, 1181859, 1
1181859, 1182019, 1
1182019, 1182091, 1
1182091, 1182131, 1
1182131, 1182147, 1
**
1182147, 1183171, 1
1183171, 1183187, 1
**
1183187, 1183227, 1
1183227, 1183299, 1
1183299, 1183459, 1
1183459, 1183779, 1
**
**---------------------------------------------------------------------+
使用* * *生成LINE2P LINE3P & LINE4P NCOPY板的命令 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**---------------------------------------------------------------------+
**
* ncopy, shift，旧set = line1p，新set = line2p，更改number =4482
0。,-0.002,0。
0, 0, 0, 0, 0, 1, 0。
**
* ncopy, shift，旧set = line1p，新set = line3p，更改number =53784
0。,-0.03,0。
0, 0, 0, 0, 0, 1, 0。
**
* ncopy, shift，旧set = line1p，新set = line4p，更改number =161352
0。,-0.12,0。
0, 0, 0, 0, 0, 1, 0。
**
* ncopy, shift，旧set = line1p，新set = line5p，更改number =394416
0。,-0.48,0。
0, 0, 0, 0, 0, 1, 0。
**
* ncopy, shift，旧set = line1p，新set = line6p，更改number =788832
0。,-1.92,0。
0, 0, 0, 0, 0, 1, 0。
**
* ncopy, shift, old set = line1p, new set = line7p, change number =1075680
0。,-7.68,0。
0, 0, 0, 0, 0, 1, 0。
**
**+----------------------------------------------------+
**填充PLATE的节点。|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+----------------------------------------------------+
**
* NFILL NSET =板
Line1p, line2p, 2,2241
LINE2P LINE3P 11, 4482年
LINE3P LINE4P 12, 8964年
LINE4P LINE5P 13, 17928年
LINE5P LINE6P 11, 35856年
LINE6P LINE7P 4, 71712年
**
**+----------------------------------------------------+
* *从直角坐标映射节点板的 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+----------------------------------------------------+
**
* NMAP, NSET =板类型=矩形
0, 0, 0, 1, 0, 0。
1。1。0。
1。1。1。
**-----------------------------------------------+
* *生成元素的网格板的区域 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**-----------------------------------------------+
** Zone1P的要素
**
*元素,ELSET = ZONE1P = CPE4类型
20117, 1182147, 1182148, 1184389, 1184388
* ELGEN ELSET = ZONE1P
20117、1024、1、1、2、2241、1024
**
** Zone2P元素
**
*元素,ELSET = ZONE2P = CPE4类型
22165, 1182131, 1182133, 1186615, 1186613
22701, 1186629, 1186631, 1191113, 1191111
22685, 1183171, 1183173, 1187655, 1187653
* ELGEN ELSET = ZONE2P
22165、8、2、1,12,4482,528
22701、512、2、1,11日,4482年,528年
22685、8、2、1,12,4482,528
**
** Zone3P元素
**
*元素,ELSET = ZONE3P = CPE4类型
28501, 1182091, 1182095, 1191059, 1191055
30215, 1235915, 1235919, 1244883, 1244879
28775, 1183187, 1183191, 1192155, 1192151
* ELGEN ELSET = ZONE3P
28501、10、4、1、18、8964、284
30215、264、4、1、8964、284
28775、10、4、1、18、8964、284
**
** Zone4P的元素
**
*元素,ELSET = ZONE4P = CPE4类型
33613, 1182019, 1182027, 1199955, 1199947
35062, 1343443, 1343451, 1361379, 1361371
33764, 1183227, 1183235, 1201163, 1201155
* ELGEN ELSET = ZONE4P
1 33613 9, 8日,22日,17928年,160年
35062、142、8、13、17928、160
1 33764 9, 8日,22日,17928年,160年
**
** Zone5P的元素
**
*元素,ELSET = ZONE5P = CPE4类型
37133, 1181859, 1181875, 1217731, 1217715
38243, 1576435, 1576451, 1612307, 1612291
37223, 1183299, 1183315, 1219171, 1219155
* ELGEN ELSET = ZONE5P
1 37133年10,16日,22日,35856年,100年
1 38243年,80年,16日,11日,35856年,100年
1 37223年10,16日,22日,35856年,100年
**
** Zone6P的元素
**
*元素,ELSET = ZONE6P = CPE4类型
39333, 1181539, 1181571, 1253283, 1253251
40113, 1970691, 1970723, 2042435, 2042403
39393, 1183459, 1183491, 1255203, 1255171
* ELGEN ELSET = ZONE6P
39333年10,32岁,1,15日,71712年,70
40113年,50岁,32岁,1,4,71712,70
39393年10,32岁,1,15日,71712年,70
**
**+-------------------------------------------------------------+
* *沿对称轴为数据输出定义节点集Y |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+-------------------------------------------------------------+
**
* NSET NSET =对称,生成
1182659,1187141, 2241
1187141,1236443, 4482
1236443,1344011, 8964
1344011,1577075, 17928
1577075,1971491, 35856
1971491, 2258339, 71712
**
**+-------------------------------------------------------------+
**用于PLATE网格区域的网格细化MPCs。|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+-------------------------------------------------------------+
**
ZONE1P和ZONE2P之间线性连接的**节点集
**
* NSET NSET = NSET1LnrMPCaZ1Z2P,生成
1182147, 1182147, 4482
* NSET NSET = NSET1LnrMPCbZ1Z2P,生成
1186629, 1186629, 4482
* NSET NSET = NSET1LnrMPCpZ1Z2P,生成
1184388, 1184388, 4482
** 1 2
* NSET NSET = NSET2LnrMPCaZ1Z2P,生成
1186629, 1187651, 2
* NSET NSET = NSET2LnrMPCbZ1Z2P,生成
1186631, 1187653, 2
* NSET NSET = NSET2LnrMPCpZ1Z2P,生成
1186630, 1187652, 2
** 1 2
* NSET NSET = NSET3LnrMPCaZ1Z2P,生成
1183171, 1183171, 4482
* NSET NSET = NSET3LnrMPCbZ1Z2P,生成
1187653, 1187653, 4482
* NSET NSET = NSET3LnrMPCpZ1Z2P,生成
1185412, 1185412, 4482
**
ZONE2P和ZONE3P之间的线性连接的**节点集
**
* NSET NSET = NSET1LnrMPCaZ2Z3P,生成
1182131, 1226951, 8964
* NSET NSET = NSET1LnrMPCbZ2Z3P,生成
1191095, 1235915, 8964
* NSET NSET = NSET1LnrMPCpZ2Z3P,生成
1186613, 1231433, 8964
** 2 3
* NSET NSET = NSET2LnrMPCaZ2Z3P,生成
1235915, 1236967, 4
* NSET NSET = NSET2LnrMPCbZ2Z3P,生成
1235919, 1236971, 4
* NSET NSET = NSET2LnrMPCpZ2Z3P,生成
1235917, 1236969, 4
** 2 3
* NSET NSET = NSET3LnrMPCaZ2Z3P,生成
1183187, 1228007, 8964
* NSET NSET = NSET3LnrMPCbZ2Z3P,生成
1192151, 1236971, 8964
* NSET NSET = NSET3LnrMPCpZ2Z3P,生成
1187669, 1232489, 8964
**
ZONE3P和ZONE4P之间的线性连接的**节点集
**
* NSET NSET = NSET1LnrMPCaZ3Z4P,生成
1182091, 1325515, 17928
* NSET NSET = NSET1LnrMPCbZ3Z4P,生成
1200019, 1343443, 17928
* NSET NSET = NSET1LnrMPCpZ3Z4P,生成
1191055, 1334479, 17928
** 3 4
* NSET NSET = NSET2LnrMPCaZ3Z4P,生成
1343443, 1344571, 8
* NSET NSET = NSET2LnrMPCbZ3Z4P,生成
1343451, 1344579, 8
* NSET NSET = NSET2LnrMPCpZ3Z4P,生成
1343447, 1344575, 8
** 3 4
* NSET NSET = NSET3LnrMPCaZ3Z4P,生成
1183227, 1326651, 17928
* NSET NSET = NSET3LnrMPCbZ3Z4P,生成
1201155, 1344579, 17928
* NSET NSET = NSET3LnrMPCpZ3Z4P,生成
1192191, 1335615, 17928
**
ZONE4P和ZONE5P之间线性连接的**节点集
**
* NSET NSET = NSET1LnrMPCaZ4Z5P,生成
1182019, 1540579, 35856
* NSET NSET = NSET1LnrMPCbZ4Z5P,生成
1217875, 1576435, 35856
* NSET NSET = NSET1LnrMPCpZ4Z5P,生成
1199947, 1558507, 35856
** 4 5
* NSET NSET = NSET2LnrMPCaZ4Z5P,生成
1576435, 1577699, 16岁
* NSET NSET = NSET2LnrMPCbZ4Z5P,生成
1576451, 1577715, 16岁
* NSET NSET = NSET2LnrMPCpZ4Z5P,生成
1576443, 1577707, 16岁
** 4 5
* NSET NSET = NSET3LnrMPCaZ4Z5P,生成
1183299, 1541859, 35856
* NSET NSET = NSET3LnrMPCbZ4Z5P,生成
1219155, 1577715, 35856
* NSET NSET = NSET3LnrMPCpZ4Z5P,生成
1201227, 1559787, 35856
* * * *
ZONE5P和ZONE6P之间线性连接的**节点集
**
* NSET NSET = NSET1LnrMPCaZ5Z6P,生成
1181859, 1898979, 71712
* NSET NSET = NSET1LnrMPCbZ5Z6P,生成
1253571, 1970691, 71712
* NSET NSET = NSET1LnrMPCpZ5Z6P,生成
1217715, 1934835, 71712
** 5 6
* NSET NSET = NSET2LnrMPCaZ5Z6P,生成
1970691, 1972259, 32
* NSET NSET = NSET2LnrMPCbZ5Z6P,生成
1970723, 1972291, 32
* NSET NSET = NSET2LnrMPCpZ5Z6P,生成
1970707, 1972275, 32
** 5 6
* NSET NSET = NSET3LnrMPCaZ5Z6P,生成
1183459, 1900579, 71712
* NSET NSET = NSET3LnrMPCbZ5Z6P,生成
1255171, 1972291, 71712
* NSET NSET = NSET3LnrMPCpZ5Z6P,生成
1219315, 1936435, 71712
**
**--------------------------------------------------------------------------+
* *“线性”类型之间的网状transitons mpc网格板的区域 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**--------------------------------------------------------------------------+
* MPC
**
线性、NSET1LnrMPCpZ1Z2P NSET1LnrMPCaZ1Z2P NSET1LnrMPCbZ1Z2P
线性、NSET2LnrMPCpZ1Z2P NSET2LnrMPCaZ1Z2P NSET2LnrMPCbZ1Z2P
线性、NSET3LnrMPCpZ1Z2P NSET3LnrMPCaZ1Z2P NSET3LnrMPCbZ1Z2P
**
线性、NSET1LnrMPCpZ2Z3P NSET1LnrMPCaZ2Z3P NSET1LnrMPCbZ2Z3P
线性、NSET2LnrMPCpZ2Z3P NSET2LnrMPCaZ2Z3P NSET2LnrMPCbZ2Z3P
线性、NSET3LnrMPCpZ2Z3P NSET3LnrMPCaZ2Z3P NSET3LnrMPCbZ2Z3P
**
线性、NSET1LnrMPCpZ3Z4P NSET1LnrMPCaZ3Z4P NSET1LnrMPCbZ3Z4P
线性、NSET2LnrMPCpZ3Z4P NSET2LnrMPCaZ3Z4P NSET2LnrMPCbZ3Z4P
线性、NSET3LnrMPCpZ3Z4P NSET3LnrMPCaZ3Z4P NSET3LnrMPCbZ3Z4P
**
线性、NSET1LnrMPCpZ4Z5P NSET1LnrMPCaZ4Z5P NSET1LnrMPCbZ4Z5P
线性、NSET2LnrMPCpZ4Z5P NSET2LnrMPCaZ4Z5P NSET2LnrMPCbZ4Z5P
线性、NSET3LnrMPCpZ4Z5P NSET3LnrMPCaZ4Z5P NSET3LnrMPCbZ4Z5P
**
线性、NSET1LnrMPCpZ5Z6P NSET1LnrMPCaZ5Z6P NSET1LnrMPCbZ5Z6P
线性、NSET2LnrMPCpZ5Z6P NSET2LnrMPCaZ5Z6P NSET2LnrMPCbZ5Z6P
线性、NSET3LnrMPCpZ5Z6P NSET3LnrMPCaZ5Z6P NSET3LnrMPCbZ5Z6P
**----------------------------------------------------------------------------------+
**生成节点集的应用程序的负载到左边的球体。|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**----------------------------------------------------------------------------------+
* NSET NSET = LFTEDGE_S,生成
2、640402、25616
* NSET NSET = RGTEDGE_S,生成
1602、642002、25616
*方程
2
NCenterS LFTEDGE_S, 1, 1.0, -1.0
2
NCenterS RGTEDGE_S, 1, 1.0, -1.0
**
2
LFTEDGE_S 2 1.0, NCenterS 2 -1.0
2
RGTEDGE_S 2 1.0, NCenterS 2 -1.0
**-----------------------------------------------------------------------------------+
**生成节点和元素集合，用于加载和BC在板上的应用。|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**-----------------------------------------------------------------------------------+
**创建node_set左边缘
* NSET NSET = LFTEDGE_P,生成
1183779, 2187747, 71712
* NSET NSET = RGTEDGE_P,生成
1253251, 2185507, 71712
**
**使用*等式将体张力施加到板的右边缘。
*方程
2
RGTEDGE_P 1 1.0, 1181539, 1, -1.0
**
* NSET NSET = NRgtMasterP
1181539
**
**--------------------------------------------------------------------------+
* *接触表面的定义 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**--------------------------------------------------------------------------+
* ELSET ELSET = ESurfPlate,生成
20117、21140、1
22165、22172、1
28501、28510、1
33613、33621、1
37133、37142、1
39333、39342、1
**
22685、22692、1
28775、28784、1
33764、33772、1
37223、37232、1
39393、39402、1
**
*表面,名称= SurfPlate类型=元素
ESurfPlate, S1
* * Fayza
**
* ELSET ELSET = ESurfBulk,生成
39333、40313、70
*表面,名称= ESurfBulk类型=元素
ESurfBulk, S4
**+------------------------------------------------------+
* *板衬底的定义材料属性 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+------------------------------------------------------+
* * *材料,NAME =板
* * *弹性、类型=各向同性
* * 210000.0, 0.3
* * *密度
* * 7800.0
*材料,NAME =板
*密度
7800.0
*弹性
210000.0, 0.3
* * *弹塑性变形
,170。e3 * * 206.8, 0.3, 4.76, 7.42, 0
**
* * *塑料硬化=约翰逊厨师
* * 382.0499392, 2441.447443, 0.655129073, 1 0
**
*塑料
422.66349, 0.00000
430.81336, 0.00100
438.80393, 0.00200
446.64406, 0.00300
454.34181, 0.00400
461.90454, 0.00500
469.33900, 0.00600
476.65141, 0.00700
483.84749, 0.00800
490.93254, 0.00900
497.91146, 0.01000
504.78881, 0.01100
511.56883, 0.01200
518.25547, 0.01300
524.85243, 0.01400
531.36316, 0.01500
537.79090, 0.01600
544.13871, 0.01700
550.40943, 0.01800
556.60577, 0.01900
562.73028, 0.02000
568.78534, 0.02100
574.77324, 0.02200
580.69613, 0.02300
586.55605, 0.02400
592.35492, 0.02500
598.09459, 0.02600
603.77680, 0.02700
609.40323, 0.02800
614.97544, 0.02900
620.49497, 0.03000
625.96324, 0.03100
631.38164, 0.03200
636.75149, 0.03300
642.07404, 0.03400
647.35051, 0.03500
652.58205, 0.03600
657.76977, 0.03700
662.91473, 0.03800
668.01797, 0.03900
673.08046, 0.04000
678.10314, 0.04100
683.08694, 0.04200
688.03271, 0.04300
692.94130, 0.04400
697.81354, 0.04500
702.65019, 0.04600
707.45201, 0.04700
712.21974, 0.04800
716.95408, 0.04900
721.65571, 0.05000
726.32528, 0.05100
730.96344, 0.05200
735.57079, 0.05300
740.14794, 0.05400
744.69546, 0.05500
749.21390, 0.05600
753.70382, 0.05700
758.16573, 0.05800
762.60014, 0.05900
767.00756, 0.06000
771.38845, 0.06100
775.74328, 0.06200
780.07250, 0.06300
784.37656, 0.06400
788.65588, 0.06500
792.91087, 0.06600
797.14194, 0.06700
*固体部分,ELSET = Zone1P材料=板
1.0
*固体部分,ELSET = Zone2P材料=板
1.0
*固体部分,ELSET = Zone3P材料=板
1.0
*固体部分,ELSET = Zone4P材料=板
1.0
*固体部分,ELSET = Zone5P材料=板
1.0
*固体部分,ELSET = Zone6P材料=板
1.0
**
**-----------------------+
**振幅定义|
**-----------------------+
**
*振幅,name = LoadRamp
0., 0。, 1.0, 1。
**
*振幅,name = BattmRamp
0.00 0.00 0.25 1.00 0.50 0.00
**
**+------------------------------------------------------+
定义SURFACE CONTACT属性。|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+------------------------------------------------------+
*接触副，相互作用=钢-钢，有限滑动
SurfPlate, SurfSphere
**
**
*表面相互作用,NAME = Steel-Steel
*摩擦
**，滑移公差= 0.005
0.6,
*表面行为，压力过闭=指数
1 e - 005721 - 5
**
**+------------------------------------------------------+
定义边界条件。|~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+------------------------------------------------------+
*边界,TYPE =位移
Line7P 2 2
LFTEDGE_P 2 2
**
**+------------------------------------------------------+
* *步骤1:应用正常负载 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+------------------------------------------------------+
**
*一步,名称= APPLY-NORMAL NLGEOM = NO
*静态
* *,直接
* *,稳定= 0.0002
0.25, 1。
*边界,TYPE =位移
* * Line7P 2 2
NCenterS 1 1
**
* CLOAD
NCenterS 2 -571.5
**
* *的交互
**
**接触控件的交互:圆柱planeconcontact
*接触控制，从=表面，主=表面，自动公差，稳定=1。
**
**
**步骤1的输出请求
**
*输出、历史、频率= 99
*联系输出,奴隶= SurfPlate = SurfSphere大师
CPRESS, CSHEAR1
*联系文件、奴隶= SurfPlate主= SurfSphere频率= 1
应该注意
*接触印刷,奴隶= SurfPlate主= SurfSphere总数=不,总结=是的,频率= 99
* * CPRESS、CSHEAR1 CDISP
应该注意
*节点打印，nset = line1p，合计= no，汇总= yes，频率=99
* * COORD
*节点打印，NSET= NCenterS，合计=NO，汇总=YES，频率=99
射频、CF、U
*输出,频率= 10,变量=预选
*最后一步
**
**+------------------------------------------------------+
* *步骤2:切向载荷应用到左边 |~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~#
**+------------------------------------------------------+
**
* step, name = apply q, nlgeom = no
施加切向载荷。
*静态
* *,直接
* *,稳定= 0.0002
0.1, 1。
*边界、OP = NEW TYPE = DISP
Line7P 2 2
LFTEDGE_P 1 1
NCenterS 1 1
* * *边界、OP = NEW TYPE = DISP,振幅= BattmRamp
* * Line7P, 1, 1, 0.02
* * LFTEDGE_P, 1, 1, 0.02
* * 1181539, 1, 1, 0.02
**
*边界、OP = NEW TYPE = DISP
Line7P 2 2
LFTEDGE_P 1 1
**
* CLOAD, OP = MOD,振幅= BattmRamp
NCenterS 1 -171.45
* DLOAD,振幅= BattmRamp
ESurfBulk、P4 -431.6087917
**
**
* *的交互
**
**接触控件的交互:圆柱planeconcontact
*接触控制，从=表面，主=表面，自动公差，稳定=1。
**
**步骤2的输出请求
**
*输出、历史、频率= 1
*联系输出,奴隶= SurfPlate = SurfSphere大师
* * CPRESS CSHEAR1
*节点输出，NSET= NCenterS
射频、CF、U
*节点输出，NSET= NRgtMasterP
U
*联系文件、奴隶= SurfPlate主= SurfSphere频率= 1
应该注意
* * CDISP
*接触印刷,奴隶= SurfPlate主= SurfSphere总数=不,总结=是的,频率= 99
* * CPRESS、CSHEAR1 CDISP
应该注意
*节点打印，nset = line1p，合计= no，汇总= yes，频率=99
COORD
*节点打印，NSET= NCenterS，合计=NO，汇总=YES，频率=99
射频、CF、U
*NODE PRINT, NSET= NRgtMasterP
U
*输出、字段、变量=预选,频率= 1
*最后一步
关于你的模型接触问题的一些意见/建议
你好,
我简单看了一下你的输入文件。我发现没什么可提的。你可能已经知道或者错过了。
-你不能在接触对中使用有限滑动选项。你能用它成功地运行作业吗?
-没有接触对定义。我想你没有在接触对选项中定义主表面和从表面名称。
-你已经使用了接触稳定选项和模型稳定，这是不推荐的。使用它们有点棘手，因为必须将接触阻尼应力与接触法向应力和切向应力进行比较，以检查接触阻尼应力是否足够高。同样，如果你在步骤定义中使用稳定选项，那么模型的静态耗散必须与系统的总应变能进行比较，以避免系统过阻尼(这会导致虚假的结果)。只是一个提示!根据我的经验，必须通过运行几次来对这个参数进行微调，以达到稳定参数的合理值。
-在ABAQUS分析手册中有一节讨论了与接触/交互模型相关的常见困难。请看一看。关于波动剪应力，它说明如下
—默认情况下，一旦建立接触就应用切向接触约束。在大多数情况下，这将产生令人满意的结果和合理的收敛。然而，经验表明，当接触建立时，在增量中施加法向约束，并在随后的增量中施加切向约束，有时可以改善收敛性，特别是在摩擦应力对接触应力有强烈影响的情况下。
在这种情况下，您可以更改默认行为，将摩擦延迟到接触点关闭的增量之后的增量。如果接触区随着分析的进展而迅速变化，则不建议这样做;在这种情况下，关闭后立即没有摩擦会导致摩擦力的快速、非物理振荡。此设置适用于各种摩擦，包括粗糙摩擦;但是它对用户子程序FRIC没有影响，它在每次增量结束时发生接触时被调用。您可以根据需要恢复默认行为。
在这一刻，在我短暂的一瞥中，我只想到了上面提到的几点。我会在周末的时候看看细节，回来的时候你可以试试上面提到的事情。也许他们会解决你的问题。
愿一切都好！
一个