Calibrate hypoplastic model

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