2 Comments
  1. David Mašín
    David Mašín 12 months ago

    Hi Sahar, please also check the initial void ratio for triaxial tests indicated by excalibre. This is the only issue I can think about, otherwise you have some error in Abaqus input.

  2. Sahar Saab Author
    Sahar Saab 12 months ago

    Hi David, this is the input file for my abaqus model. This input file is for the sand example “lapProtocalTemplate-Sand” shared on Excalibre automatic calibration page. I got the hypoplastic parameters from Excalibre and input them into Abaqus. The peak of deviatoric stress-strain curve from Abaqus is much higher than the experimental curve and the Excalibre curve. I can’t find the mistake in my input.

    *Heading
    ** Job name: CD-50-Excalibre-Sample-Sand Model name: CD-50-Excalibre-Sample-Sand
    ** Generated by: Abaqus/CAE 2016
    *Preprint, echo=NO, model=NO, history=NO, contact=NO
    **
    ** PARTS
    **
    *Part, name=Part-1
    *Node
    1, 0., 0.
    2, 0.100000001, 0.
    3, 0., 0.100000001
    4, 0.100000001, 0.100000001
    5, 0.0500000007, 0.
    6, 0.100000001, 0.0500000007
    7, 0.0500000007, 0.100000001
    8, 0., 0.0500000007
    *Element, type=CAX8
    1, 1, 2, 4, 3, 5, 6, 7, 8
    *Nset, nset=Set-1, generate
    1, 8, 1
    *Elset, elset=Set-1
    1,
    *Nset, nset=top-n
    3, 4, 7
    *Nset, nset=bot-n
    1, 2, 5
    *Nset, nset=axis
    1, 3, 8
    ** Section: Section-1
    *Solid Section, elset=Set-1, material=Material-1
    ,
    *End Part
    **
    **
    ** ASSEMBLY
    **
    *Assembly, name=Assembly
    **
    *Instance, name=Part-1-1, part=Part-1
    *End Instance
    **
    *Elset, elset=EL, instance=Part-1-1
    1,
    *Elset, elset=_Surf-1_S2, internal, instance=Part-1-1
    1,
    *Elset, elset=_Surf-1_S3, internal, instance=Part-1-1
    1,
    *Surface, type=ELEMENT, name=Surf-1
    _Surf-1_S2, S2
    _Surf-1_S3, S3
    *End Assembly
    **
    ** MATERIALS
    **
    *Material, name=Material-1
    *Density
    1.0927,
    *Depvar
    13,
    *User Material, constants=16, unsymm
    36.5, 1.,49846., 0.243, 0.572, 1.144, 1.372, 0.09
    4.5, 0., 0., 0., 0., 0., 0., 0.51
    **
    ** BOUNDARY CONDITIONS
    **
    ** Name: axis-bc Type: Displacement/Rotation
    *Boundary
    Part-1-1.axis, 1, 1
    ** Name: bot-bc Type: Displacement/Rotation
    *Boundary
    Part-1-1.bot-n, 2, 2
    ** Name: top-bc Type: Displacement/Rotation
    *Boundary
    Part-1-1.top-n, 2, 2
    **
    ** PREDEFINED FIELDS
    **
    ** Name: Predefined Field-1 Type: Geostatic stress
    *Initial Conditions, type=STRESS, GEOSTATIC
    Part-1-1.Set-1, -50., 0.1, -50., 0., 1., 1.
    ** Name: Predefined Field-2 Type: Void ratio
    *Initial Conditions, TYPE=RATIO
    Part-1-1.Set-1, 0.51
    ** —————————————————————-
    **
    ** STEP: Step-1
    **
    *Step, name=Step-1, nlgeom=NO
    *Geostatic
    **
    ** LOADS
    **
    ** Name: Load-1 Type: Pressure
    *Dsload
    Surf-1, P, 50.
    **
    ** OUTPUT REQUESTS
    **
    *Restart, write, frequency=0
    **
    ** FIELD OUTPUT: F-Output-1
    **
    *Output, field
    *Node Output
    U,
    *Element Output, directions=YES
    E, S
    **
    ** HISTORY OUTPUT: H-Output-1
    **
    *Output, history, variable=PRESELECT
    *End Step
    ** —————————————————————-
    **
    ** STEP: Step-2
    **
    *Step, name=Step-2, nlgeom=NO, inc=1000
    *Static
    0.001, 1., 1e-15, 1.
    **
    ** BOUNDARY CONDITIONS
    **
    ** Name: top-bc Type: Displacement/Rotation
    *Boundary
    Part-1-1.top-n, 2, 2, -0.05
    **
    ** LOADS
    **
    ** Name: Load-1 Type: Pressure
    *Dsload
    Surf-1, P, 50.
    **
    ** OUTPUT REQUESTS
    **
    *Restart, write, frequency=0
    **
    ** FIELD OUTPUT: F-Output-1
    **
    *Output, field
    *Node Output
    U,
    *Element Output, directions=YES
    E, S
    **
    ** HISTORY OUTPUT: H-Output-1
    **
    *Output, history, variable=PRESELECT
    *End Step

Leave a reply

©2020 SoilModels

Log in with your credentials

Forgot your details?