Skip to main content

Help Required for cam clay soil model

Submitted by rahul_umkc on

Hello everyone,

I am trying to simulate the force vs displacement response when cam clay soil model containing rigid foundation is acted by the concentrated load at the interface of rigid body with the soil. I don't know why i am receiving "time increment error".

Input file

** ASSEMBLY

**

*Assembly, name=Assembly

** 

*Instance, name=Part-1-1, part=Part-1

*End Instance

** 

*Nset, nset=Set-3, instance=Part-1-1, generate

   1,  661,    1

*Elset, elset=Set-3, instance=Part-1-1, generate

   1,  200,    1

*Instance, name=Part-1-2, part=Part-1

*End Instance

**

*Nset, nset=Set-3, instance=Part-1-2, generate

   1,  661,    1

*Elset, elset=Set-3, instance=Part-1-2, generate

   1,  200,    1

*Instance, name=Part-2-1, part=Part-2

        8.25,           4.,           0.

*End Instance

** 

*Nset, nset=Set-1, instance=Part-2-1

 2,

*Nset, nset=Set-2, instance=Part-1-2

 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223, 224, 225, 226

 227, 228, 229, 230, 231, 622, 625, 627, 629, 631, 633, 635, 637, 639, 641, 643

 645, 647, 649, 651, 653, 655, 657, 659, 661

*Elset, elset=Set-2, instance=Part-1-2, generate

 181,  200,    1

*Nset, nset=Set-3, instance=Part-1-2, generate

   1,  661,    1

*Elset, elset=Set-3, instance=Part-1-2, generate

   1,  200,    1

*Nset, nset=Set-4, instance=Part-2-1

 12,

*Nset, nset=Set-5, instance=Part-1-2

  21,  42,  63,  84, 105, 126, 147, 168, 189, 210, 231, 291, 332, 373, 414, 455

 496, 537, 578, 619, 660

*Elset, elset=Set-5, instance=Part-1-2, generate

  20,  200,   20

*Nset, nset=Set-6, instance=Part-1-2

   1,  22,  43,  64,  85, 106, 127, 148, 169, 190, 211, 235, 295, 336, 377, 418

 459, 500, 541, 582, 623

*Elset, elset=Set-6, instance=Part-1-2, generate

   1,  181,   20

*Rigid Body, ref node=Part-2-1.Part-2-RefPt_, elset=Part-2-1.Part-2

*End Assembly

*Amplitude, name=Amp-1

             0.,              1.,            10.,              1.

**

** MATERIALS

**

*Material, name=Material-1

*Clay Plasticity

 0.1174, 1.,15., 1., 1.

*Permeability, specific=0.0361

1, 0.889

*Porous Elastic

 0.026, 0.28,   0.

*Initial Conditions, TYPE=RATIO

Set-3, 0.889

*Initial Conditions, type=STRESS, GEOSTATIC

Set-3, -30., 4, -30., -1, 1., 1.

**

** INTERACTION PROPERTIES

**

*Surface Interaction, name=IntProp-1

1.,

*Friction

0.,

**

** BOUNDARY CONDITIONS

**

** Name: BC-1 Type: Displacement/Rotation

*Boundary

Set-2, 1, 1

Set-2, 2, 2

** Name: BC-2 Type: Displacement/Rotation

*Boundary

Set-5, 1, 1

** Name: BC-3 Type: Displacement/Rotation

*Boundary

Set-6, 1, 1

**

 *constraint controls, print=yes

** ----------------------------------------------------------------

** ----------------------------------------------------------------

** **

* STEP: Step-1

**

*Step, name=Step-1, nlgeom=YES

*Geostatic, utol=1.

1., 1., 1e-05, 1.

**

** LOADS

**

** Name: Load-1   Type: Concentrated force

*Cload, amplitude=Amp-1

Set-4, 2, -30.

**

** OUTPUT REQUESTS

**

*Restart, write, frequency=0

**

** FIELD OUTPUT: F-Output-1

**

*Output, field, variable=PRESELECT

**

** HISTORY OUTPUT: H-Output-1

**

*Output, history, variable=PRESELECT

*End Step

** STEP: Step-2

**

*Step, name=Step-2

*Static, direct

10.,100.,

**

** OUTPUT REQUESTS

**

*Restart, write, frequency=0

**

** FIELD OUTPUT: F-Output-1

**

*Output, field, variable=PRESELECT

**

** HISTORY OUTPUT: H-Output-1

**

*Output, history, variable=PRESELECT

*End Step

Any help would be highly appreciated......

Regards

Rahul