请问大家关于桩单元的弹簧刚度的问题

winni

关于桩单元的弹簧刚度问题，一直很迷惑。

我做了个抗拔桩的例子
1.当cs_sk=1e11 cs_scoh=0 cs_sfric=0 cs_nk=0 cs_ncoh=0 cs_nfric=0

cycle 0
print sel nod stiff range id 2 3 ；2为桩端点，3为桩身点

node-id      x-comp       y-comp       z-comp   
---------   -----------  -----------  -----------
        2 (  5.277e+011,  0.000e+000,  0.000e+000)
        3 (  1.005e+012,  0.000e+000,  0.000e+000)

2、当cs_sk=0 cs_scoh=0 cs_sfric=0 cs_nk=1e11 cs_ncoh=0 cs_nfric=0

node-id      x-comp       y-comp       z-comp   
---------   -----------  -----------  -----------
        2 (  5.027e+011,  2.500e+010,  2.500e+010)
        3 (  1.005e+012,  0.000e+000,  0.000e+000)

3、当cs_sk=1e11 cs_scoh=0 cs_sfric=0 cs_nk=1e11 cs_ncoh=0 cs_nfric=0

node-id      x-comp       y-comp       z-comp   
---------   -----------  -----------  -----------
        2 (  5.277e+011,  2.500e+010,  2.500e+010)
        3 (  1.005e+012,  0.000e+000,  0.000e+000)


从上面看来，1、得到的节点刚度与设置的刚度值相比大了1个数量级。
2、设置法向刚度nk为什么对桩端节点还会产生x-comp.
3、无论设置法向还是切向刚度，对桩身中间的节点刚度影响是一样的。


=================================================================

另外，根据FLAC的理论公式，Fs/L=CSstiff(up-um),可是我算出来的刚度，跟设置的刚度相比，少了两个数量级。

这样就出现了3个刚度，并且我怎么也没看出他们之间有什么关系。而模型的实际弹簧刚度到底以哪个才是正确的呢？？？我很迷惑啊。大家来讨论一下吧。

eagleming

x-comp 应该是力的单位，是刚度乘面积，和你桩的面积有关。大家在说说，我理解的也不一定对

winni

2# eagleming

我的模型跟彭文斌书上关于桩结构单元那一章的模型基本上一样，只是改成了抗拔桩。桩直径是1米，面积是0.785平方米，桩端弹簧面积1平方米。
书上p184,sk=nk=1.3e11,sfric=10,用print sel nod stiff range id 2获得5.4e11N/m，弹簧刚度为5.4e11N/m3,与设置的刚度1.3e11还是有距离。

capture

期待解决中，楼主根据Flac理论公式Fs/L=CSstiff(up-um)具体是怎样计算的呢？

lhy407

学习下。。。。。。

饿狼鼎

想问问楼主，壳单元和实体单元不考虑滑移作用，是通过接触面来实现的，还是通过merge，或者布尔运算实现的呢~？