求教：弥散裂缝模型

hym-0387

请问有没有做过三维的混凝土弥散裂缝模型的？按默认的参数设置，模型不收敛。然后改动参数后，能收敛，但是得出来的应力-应变曲线如下图，没有下降段。不知道是参数的问题，还是有其他需要注意的地方没考虑到，烦请知道的高手谈一谈。由于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 编辑 ]