2 Comments
  1. David Mašín
    David Mašín 2 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 2 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

©2019 SoilModels

Log in with your credentials

Forgot your details?