求助：接触分析中出现的问题

清舞飞扬

我做的是一个钢针刺进皮肤的数值模拟。我暂时把皮肤处理成弹塑性材料，考虑了shear failure
但是在计算过程中老师出现下面的错误，导致计算无法进行下去
The elements contained in element set ErrElemExcessDistortion-Step1 have distorted excessively.

There is only one excessively distorted element

The elements contained in element set ErrElemExcessDistortion-Step1 have distorted excessively.
哪位帮我看看，这问题怎么解决啊？
还有ABAQUS里面带的Adaptive mesh control,怎么用啊？我试了一下，好像没什么用啊
选参数时我都是采取的默认，我应该改一下什么参数，能使效果明显一点啊？
拜托了，救命啊

清舞飞扬

这是我的inp文件
*Heading
** Job name: Job-3-29 Model name: Model-1
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=Part-1
*End Part
**  
*Part, name=Part-2
*End Part
**  
**
** ASSEMBLY
**
*Assembly, name=Assembly
**  
*Instance, name=Part-1-1, part=Part-1
*Node
      1,           0.,  0.100000001
      2,           0.,           0.

..............
      92, 0.0149999997, 0.0914999992
*Element, type=RAX2
1,  1,  5
2,  5,  6
3,  6,  7
4,  7,  8
5,  8,  9
6,  9, 10
7, 10, 11
8, 11, 12
9, 12, 13
10, 13, 14
11, 14, 15
12, 15, 16
13, 16, 17
14, 17, 18
15, 18, 19
16, 19, 20
17, 20, 21
18, 21, 22
19, 22, 23
20, 23, 24
21, 24, 25
22, 25, 26
23, 26, 27
24, 27, 28
25, 28, 29
26, 29, 30
27, 30, 31
28, 31, 32
29, 32, 33
30, 33, 34
31, 34, 35
32, 35, 36
33, 36, 37
34, 37, 38
35, 38, 39
36, 39, 40
37, 40, 41
38, 41, 42
39, 42, 43
40, 43, 44
41, 44, 45
42, 45, 46
43, 46, 47
44, 47, 48
45, 48, 49
46, 49, 50
47, 50, 51
48, 51, 52
49, 52, 53
50, 53,  2
51,  2, 54
52, 54, 55
53, 55, 56
54, 56, 57
55, 57, 58
56, 58, 59
57, 59, 60
58, 60, 61
59, 61, 62
60, 62, 63
61, 63, 64
62, 64, 65
63, 65, 66
64, 66, 67
65, 67, 68
66, 68, 69
67, 69, 70
68, 70, 71
69, 71, 72
70, 72, 73
71, 73, 74
72, 74, 75
73, 75, 76
74, 76, 77
75, 77, 78
76, 78, 79
77, 79, 80
78, 80, 81
79, 81, 82
80, 82,  3
81,  4, 83
82, 83,  1
83,  3, 84
84, 84, 85
85, 85, 86
86, 86, 87
87, 87, 88
88, 88, 89
89, 89, 90
90, 90, 91
91, 91, 92
92, 92,  4
*Node
     93,           0.,  0.100000001,           0.
*Nset, nset=Part-1-1-RefPt_, internal
93,
*Elset, elset=Part-1-1, generate
  1,  92,   1
*End Instance
**  
*Instance, name=Part-2-1, part=Part-2
*Node
      1,           0., -0.0425000004
      2, 0.0175000001, -0.0425000004
      3, 0.0175000001,           0.
      4,           0.,           0.
      5, 0.0949999988,           0.
      6,           0., -0.0850000009
      7, 0.0949999988, -0.0850000009
      8,  0.150000006, -0.150000006
      9,  0.150000006,           0.
     10,           0., -0.150000006
     11, 0.000682924525, -0.0425000004
     12, 0.00136513286, -0.0425000004
     13, 0.0020475504, -0.0425000004
     14, 0.00273080496, -0.0425000004
..............
         
   5696, 5850,  148,    6,  278
*Nset, nset=_PickedSet2, internal, generate
    1,  5850,     1
*Elset, elset=_PickedSet2, internal, generate
    1,  5696,     1
** Region: (Section-1:Picked), (Controls:EC-1)
*Elset, elset=_PickedSet2, internal, generate
    1,  5696,     1
** Section: Section-1
*Solid Section, elset=_PickedSet2, controls=EC-1, material=Material-1
1.,
*End Instance
**  
*Nset, nset=_PickedSet9, internal, instance=Part-2-1
   8,  10, 310, 311, 312, 313, 314, 315, 316, 317, 318, 319, 320, 321, 322, 323
324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 336, 337, 338, 339
340, 341, 342, 343, 344, 345, 346, 347, 348, 349, 350, 351, 352, 353, 354, 355
356, 357, 358, 359, 360, 361, 362, 363, 364, 365, 366, 367, 368, 369, 370, 371
372,
*Elset, elset=_PickedSet9, internal, instance=Part-2-1
2913, 2914, 2915, 2916, 2917, 2918, 2919, 2920, 2921, 2922, 2923, 2924, 3981, 4014, 4047, 4080
4113, 4146, 4179, 4212, 4245, 4278, 4311, 4344, 4377, 4410, 4443, 4476, 4509, 4542, 4575, 4608
4641, 4674, 4707, 4740, 4773, 4806, 4839, 4872, 4905, 4938, 4971, 5004, 5037, 5070, 5103, 5136
5169, 5202, 5235, 5268, 5301, 5334, 5367, 5400, 5433, 5466, 5499, 5532, 5565, 5598, 5631, 5664
*Nset, nset=_PickedSet10, internal, instance=Part-1-1
93,
*Nset, nset=_PickedSet68, internal, instance=Part-2-1
    1,    2,  ............

, ......................3218
*Elset, elset=_PickedSet68, internal, instance=Part-2-1, generate
    1,  2912,     1
*Nset, nset=_PickedSet78, internal, instance=Part-2-1
    1,   .........

, 3218
*Elset, elset=_PickedSet78, internal, instance=Part-2-1, generate
    1,  2912,     1
*Nset, nset=_PickedSet81, internal, instance=Part-2-1, generate
    1,  5850,     1
*Elset, elset=_PickedSet81, internal, instance=Part-2-1, generate
    1,  5696,     1
*Nset, nset=_PickedSet82, internal, instance=Part-2-1, generate
    1,  5850,     1
*Elset, elset=_PickedSet82, internal, instance=Part-2-1, generate
    1,  5696,     1
*Nset, nset=_PickedSet87, internal, instance=Part-2-1, generate
    1,  5850,     1
*Elset, elset=_PickedSet87, internal, instance=Part-2-1, generate
    1,  5696,     1
*Nset, nset=_PickedSet88, internal, instance=Part-2-1, generate
    1,  5850,     1
*Elset, elset=_PickedSet88, internal, instance=Part-2-1, generate
    1,  5696,     1
*Nset, nset=_PickedSet106, internal, instance=Part-2-1, generate
    1,  5850,     1
*Elset, elset=_PickedSet106, internal, instance=Part-2-1, generate
    1,  5696,     1
*Nset, nset=_PickedSet123, internal, instance=Part-2-1
    1,    2,    3,   .........................

     ........................3218
*Elset, elset=_PickedSet123, internal, instance=Part-2-1, generate
    1,  2912,     1
*Nset, nset=_PickedSet138, internal, instance=Part-2-1, generate
    1,  5850,     1
*Elset, elset=_PickedSet138, internal, instance=Part-2-1, generate
    1,  5696,     1
*Elset, elset=__PickedSurf122_SNEG, internal, instance=Part-1-1
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, 83, 84
85, 86, 87, 88, 89, 90, 91, 92
*Surface, type=ELEMENT, name=_PickedSurf122, internal
__PickedSurf122_SNEG, SNEG
*Surface, type=NODE, name=_PickedSet123_CNS_, internal
_PickedSet123, 1.
*Rigid Body, ref node=Part-1-1.Part-1-1-RefPt_, elset=Part-1-1.Part-1-1
*End Assembly
**
** ELEMENT CONTROLS
**
*Section Controls, name=EC-1, DISTORTION CONTROL, length ratio=0.05, hourglass=ENHANCED
1., 1., 1.
*Amplitude, name=Amp-1
0., 0., 0.01, 1.
**
** MATERIALS
**
*Material, name=Material-1
*Density
1.2e-09,
*Elastic
15., 0.4
*Plastic
  8.68,   0.
  9.41, 0.04
   10., 0.19
10.05,  0.25
*Shear Failure
0.25
**
** INTERACTION PROPERTIES
**
*Surface Interaction, name=IntProp-1
**
** BOUNDARY CONDITIONS
**
** Name: BC-1 Type: Displacement/Rotation
*Boundary
_PickedSet9, 1, 1
_PickedSet9, 2, 2
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1
*Dynamic, Explicit
, 0.01
*Bulk Viscosity
0.06, 1.2
** Mass Scaling: Semi-Automatic
**               Whole Model
*Fixed Mass Scaling, factor=10000.
**
** BOUNDARY CONDITIONS
**
** Name: BC-2 Type: Displacement/Rotation
*Boundary, amplitude=Amp-1
_PickedSet10, 1, 1
_PickedSet10, 2, 2, -0.03
_PickedSet10, 6, 6
**
** INTERACTIONS
**
** Interaction: Int-1
*Contact Pair, interaction=IntProp-1, mechanical constraint=KINEMATIC, cpset=Int-1
_PickedSurf122, _PickedSet123_CNS_
**
** OUTPUT REQUESTS
**
*Restart, write, number interval=1, time marks=NO
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*Node Output
A, RF, U, V
*Element Output, directions=YES
LE, PE, PEEQ, PS, S, STATUS, VE, VEEQ, VS
*Contact Output
CSTRESS,
**
** HISTORY OUTPUT: H-Output-1
**
*Output, history, variable=PRESELECT
*End Step

没什么

楼主问题怎么解决的？我现在也出现相同的问题

ydliu_zn

1) 将单元换成对应的减速积分单元（如果有的话）。
2) 使用网格重划分技术，比如说你提到的ALE法，对应的关键字为*ADAPTIVE MESH, ELSET=ADAPT看看帮助文件的12.2.2就可以了。你还可以用其他重划分网格的方法。
3) 细化网格，一可以防止沙漏，二可以减少distortion。

另外，钢针刺入皮肤，尺寸貌似很小，针头太尖的话本身就不容易收敛，可以适当修改修改模型。