怎样用单单元模拟单轴压缩试验

ilovechu

我现在在做混凝土的单轴压缩试验的模拟来验证材料特性。请问在建立一个10x10x10mm的solid单元后，是用BOUNDARY_PRESCRIBED_MOTION控制节点的位移进行加载还是使用load segment控制加载压力？
我现在使用BOUNDARY_PRESCRIBED_MOTION来实现匀速加载，但是最后的结果得到的峰值压应力是0.6e-3MPa，比试验小了几个数量级，不知道是什么原因？我的单位是mm-g-ms，下面是部分k文件：
*SECTION_SOLID
$#   secid    elform       aet   
         1         1         0
*BOUNDARY_PRESCRIBED_MOTION_SET
$#    nsid       dof       vad      lcid        sf       vid     death     birth
         2         3         0         1       1.0         01.00000E28       0.0
*SET_NODE_LIST_TITLE
Node set 2
$#     sid       da1       da2       da3       da4    solver      
         2       0.0       0.0       0.0       0.0MECH
$#    nid1      nid2      nid3      nid4      nid5      nid6      nid7      nid8
         5         6         8         7         0         0         0         0
*DEFINE_CURVE
$#    lcid      sidr       sfa       sfo      offa      offo    dattyp     lcint
         1         0       1.0       1.0       0.0       0.0         0         0
$#                a1                  o1  
                 0.0              -0.001
                10.0              -0.001
              1000.0              -0.001