求教：弥散裂缝模型
请问有没有做过三维的混凝土弥散裂缝模型的？按默认的参数设置，模型不收敛。然后改动参数后，能收敛，但是得出来的应力-应变曲线如下图，没有下降段。不知道是参数的问题，还是有其他需要注意的地方没考虑到，烦请知道的高手谈一谈。由于INP文件大过附件的要求，故粘贴出来。（节点坐标有删减以适应字数要求）*Heading
** Job name: concrete-smeared-cracking-model Model name: Model-1
** Generated by: Abaqus/CAE Version 6.8-1
*Preprint, echo=NO, model=NO, history=NO, contact=NO
**
** PARTS
**
*Part, name=concrete
*Node
1, 200., 200., 50.
2, 200., 190., 50.
3, 200., 180., 50.
4, 200., 170., 50.
2440, 10., 170., 30.
2441, 10., 160., 30.
2442, 10., 150., 30.
2443, 10., 140., 30.
2444, 10., 130., 30.
2445, 10., 120., 30.
2446, 10., 110., 30.
2447, 10., 100., 30.
2462, 10., 160., 20.
2463, 10., 150., 20.
2464, 10., 140., 20.
2465, 10., 130., 20.
2466, 10., 120., 20.
2644, 0., 20., 0.
2645, 0., 10., 0.
2646, 0., 0., 0.
*Element, type=C3D8R
1,127,128,149,148, 1, 2, 23, 22
2,128,129,150,149, 2, 3, 24, 23
3,129,130,151,150, 3, 4, 25, 24
4,130,131,152,151, 4, 5, 26, 25
1999, 2623, 2624, 2645, 2644, 2497, 2498, 2519, 2518
2000, 2624, 2625, 2646, 2645, 2498, 2499, 2520, 2519
*Nset, nset=_PickedSet2, internal, generate
1,2646, 1
*Elset, elset=_PickedSet2, internal, generate
1,2000, 1
** Section: concrete
*Solid Section, elset=_PickedSet2, material=concrete
*End Part
**
*Part, name="end plate"
*Node
1, 10., 100., 50.
2, 10., 80., 50.
3, 10., 60., 50.
4, 10., 40., 50.
5, 10., 20., 50.
6, 10., 0., 50.
7, 10., -20., 50.
8, 10., -40., 50.
9, 10., -60., 50.
10, 10., -80., 50.
11, 10., -100., 50.
12, 10., 100., 33.3333321
13, 10., 80., 33.3333321
14, 10., 60., 33.3333321
15, 10., 40., 33.3333321
16, 10., 20., 33.3333321
17, 10., 0., 33.3333321
18, 10., -20., 33.3333321
19, 10., -40., 33.3333321
20, 10., -60., 33.3333321
21, 10., -80., 33.3333321
22, 10., -100., 33.3333321
23, 10., 100., 16.666666
24, 10., 80., 16.666666
25, 10., 60., 16.666666
26, 10., 40., 16.666666
27, 10., 20., 16.666666
28, 10., 0., 16.666666
29, 10., -20., 16.666666
30, 10., -40., 16.666666
31, 10., -60., 16.666666
32, 10., -80., 16.666666
33, 10., -100., 16.666666
34, 10., 100., 0.
35, 10., 80., 0.
36, 10., 60., 0.
37, 10., 40., 0.
38, 10., 20., 0.
39, 10., 0., 0.
40, 10., -20., 0.
41, 10., -40., 0.
42, 10., -60., 0.
43, 10., -80., 0.
44, 10., -100., 0.
45, -10., 100., 50.
46, -10., 80., 50.
47, -10., 60., 50.
48, -10., 40., 50.
49, -10., 20., 50.
50, -10., 0., 50.
51, -10., -20., 50.
52, -10., -40., 50.
53, -10., -60., 50.
54, -10., -80., 50.
55, -10., -100., 50.
56, -10., 100., 33.3333321
57, -10., 80., 33.3333321
58, -10., 60., 33.3333321
59, -10., 40., 33.3333321
60, -10., 20., 33.3333321
61, -10., 0., 33.3333321
62, -10., -20., 33.3333321
63, -10., -40., 33.3333321
64, -10., -60., 33.3333321
65, -10., -80., 33.3333321
66, -10., -100., 33.3333321
67, -10., 100., 16.666666
68, -10., 80., 16.666666
69, -10., 60., 16.666666
70, -10., 40., 16.666666
71, -10., 20., 16.666666
72, -10., 0., 16.666666
73, -10., -20., 16.666666
74, -10., -40., 16.666666
75, -10., -60., 16.666666
76, -10., -80., 16.666666
77, -10., -100., 16.666666
78, -10., 100., 0.
79, -10., 80., 0.
80, -10., 60., 0.
81, -10., 40., 0.
82, -10., 20., 0.
83, -10., 0., 0.
84, -10., -20., 0.
85, -10., -40., 0.
86, -10., -60., 0.
87, -10., -80., 0.
88, -10., -100., 0.
*Element, type=C3D8R
1, 45, 46, 57, 56,1,2, 13, 12
2, 46, 47, 58, 57,2,3, 14, 13
3, 47, 48, 59, 58,3,4, 15, 14
4, 48, 49, 60, 59,4,5, 16, 15
5, 49, 50, 61, 60,5,6, 17, 16
6, 50, 51, 62, 61,6,7, 18, 17
7, 51, 52, 63, 62,7,8, 19, 18
8, 52, 53, 64, 63,8,9, 20, 19
9, 53, 54, 65, 64,9, 10, 21, 20
10, 54, 55, 66, 65, 10, 11, 22, 21
11, 56, 57, 68, 67, 12, 13, 24, 23
12, 57, 58, 69, 68, 13, 14, 25, 24
13, 58, 59, 70, 69, 14, 15, 26, 25
14, 59, 60, 71, 70, 15, 16, 27, 26
15, 60, 61, 72, 71, 16, 17, 28, 27
16, 61, 62, 73, 72, 17, 18, 29, 28
17, 62, 63, 74, 73, 18, 19, 30, 29
18, 63, 64, 75, 74, 19, 20, 31, 30
19, 64, 65, 76, 75, 20, 21, 32, 31
20, 65, 66, 77, 76, 21, 22, 33, 32
21, 67, 68, 79, 78, 23, 24, 35, 34
22, 68, 69, 80, 79, 24, 25, 36, 35
23, 69, 70, 81, 80, 25, 26, 37, 36
24, 70, 71, 82, 81, 26, 27, 38, 37
25, 71, 72, 83, 82, 27, 28, 39, 38
26, 72, 73, 84, 83, 28, 29, 40, 39
27, 73, 74, 85, 84, 29, 30, 41, 40
28, 74, 75, 86, 85, 30, 31, 42, 41
29, 75, 76, 87, 86, 31, 32, 43, 42
30, 76, 77, 88, 87, 32, 33, 44, 43
*Nset, nset=_PickedSet2, internal, generate
1,88, 1
*Elset, elset=_PickedSet2, internal, generate
1,30, 1
** Section: steel
*Solid Section, elset=_PickedSet2, material=steel
*End Part
**
**
** ASSEMBLY
**
*Assembly, name=Assembly
**
*Instance, name=concrete-1, part=concrete
*End Instance
**
*Instance, name="end plate-1", part="end plate"
210., 100., 0.
*End Instance
**
*Instance, name="end plate-2", part="end plate"
-10., 100., 0.
*End Instance
**
*Nset, nset="Set-plot end plane", instance="end plate-1", generate
1,44, 1
*Elset, elset="Set-plot end plane", instance="end plate-1", generate
1,30, 1
*Nset, nset=Set-X, instance="end plate-1", generate
1,44, 1
*Elset, elset=Set-X, instance="end plate-1", generate
1,30, 1
*Nset, nset=concrete, instance=concrete-1, generate
1,2646, 1
*Elset, elset=concrete, instance=concrete-1, generate
1,2000, 1
*Nset, nset=Set-nX, instance="end plate-2", generate
45,88, 1
*Elset, elset=Set-nX, instance="end plate-2", generate
1,30, 1
*Elset, elset="Set-plot element", instance=concrete-1
403,
*Elset, elset="_Surf-end plate-x_S2", internal, instance="end plate-2", generate
1,30, 1
*Surface, type=ELEMENT, name="Surf-end plate-x"
"_Surf-end plate-x_S2", S2
*Elset, elset="_Surf-end plate-nx_S1", internal, instance="end plate-1", generate
1,30, 1
*Surface, type=ELEMENT, name="Surf-end plate-nx"
"_Surf-end plate-nx_S1", S1
*Elset, elset=_Surf-concrete-x_S2, internal, instance=concrete-1, generate
1, 100, 1
*Surface, type=ELEMENT, name=Surf-concrete-x
_Surf-concrete-x_S2, S2
*Elset, elset=_Surf-concrete-nx_S1, internal, instance=concrete-1, generate
1901,2000, 1
*Surface, type=ELEMENT, name=Surf-concrete-nx
_Surf-concrete-nx_S1, S1
** Constraint: Constraint-1
*Tie, name=Constraint-1, adjust=yes
Surf-concrete-x, "Surf-end plate-nx"
** Constraint: Constraint-2
*Tie, name=Constraint-2, adjust=yes
Surf-concrete-nx, "Surf-end plate-x"
*End Assembly
**
** MATERIALS
**
*Material, name=concrete
*Concrete
11.38, 0.
14.3011, 8.61e-06
17.5, 3.88e-05
20.5108, 7.48e-05
22.5806,0.00014
24.6505, 0.000189
26.5323, 0.000261
28.0376, 0.000344
29.543, 0.000427
30.4839, 0.000511
31.4247, 0.000628
31.8011, 0.000714
32.5538, 0.000837
32.9301,0.00107
*Failure Ratios
1.16, 0.1,1.28, 0.333
*Shear Retention
1., 0.01
*Tension Stiffening, type=DISPLACEMENT
0.05,
*Density
2.4e-09,
*Elastic
32500., 0.2
*Material, name=steel
*Elastic
210000., 0.3
**
** BOUNDARY CONDITIONS
**
** Name: BC-X Type: Displacement/Rotation
*Boundary
Set-X, 1, 1
Set-X, 2, 2
Set-X, 3, 3
Set-X, 4, 4
Set-X, 5, 5
Set-X, 6, 6
** Name: BC-nX Type: Displacement/Rotation
*Boundary
Set-nX, 1, 1
Set-nX, 2, 2
Set-nX, 3, 3
Set-nX, 4, 4
Set-nX, 5, 5
Set-nX, 6, 6
** ----------------------------------------------------------------
**
** STEP: Step-1
**
*Step, name=Step-1, nlgeom=YES, inc=100000
*Static
0.001, 1., 1e-15, 0.1
**
** BOUNDARY CONDITIONS
**
** Name: BC-X Type: Displacement/Rotation
*Boundary
Set-X, 1, 1, -2.
**
** OUTPUT REQUESTS
**
*Restart, write, frequency=0
**
** FIELD OUTPUT: F-Output-1
**
*Output, field
*Node Output
CF, RF, U
*Element Output, directions=YES
LE, PE, PEEQ, PEEQT, S
**
** HISTORY OUTPUT: H-Output-element
**
*Output, history
*Element Output, elset="Set-plot element"
E11, MISES
**
** HISTORY OUTPUT: H-Output-end plane
**
*Node Output, nset="Set-plot end plane"
RF1, U1
*End Step
[ 本帖最后由 hym-0387 于 2009-5-21 10:13 编辑 ]
页:
[1]