《JTG3363-2019》附录L算例计算问题

    1、由于实行了新基础规范，原07规范附录P调整为了附录L，桩基础效应计算原理也未变。
; O  h/ Y. o% o+ O2 c' J- M    2、由于新规范的发行，在阅读附录L条文说明（0.7规范一样）P185-186页时发现个算例。
    3、问题所在：本人通过计算发现精确解x0=3.913mm，Mmax=344.962kN*m，根据规范解x0=4.1203mm，Mmax=354.43kN*m
& L1 |, j8 j8 s2 ?# p5 h2 `         均与提供的算例数值不符，算例精确解x0=4.35mm，Mmax=335.49kN*m，本规范解x0=4.44mm，Mmax=340.31kN*m。
         如下图：
3 ^+ l" Y- f; N+ d* p         

雪狼儿88 发表于 2020-5-11 09:35
1 |8 M7 {* z8 x* {3 T! ]8 _楼主这自己计算的精确解不是参照规范吗，还有两个规范解值是分别按07和19规范？？？

    1、我计算的精确解是采用微分方程数值解法直接用电脑解微分方程计算出来的。
) |0 @- K9 ^& U( {& k  Z: K    2、“本规范解x0=4.44”是规范原话（07、19规范都是一样的，本附录未变动）。
    3、“根据规范解x0=4.1203”是我采用规范附录L（或附录P）计算过程，但是无量纲系数采用的幂级数通解递推公式得出。
7 |. ~2 h) t; ?  M, r" j7 v    4、3楼是我完全按照规范计算过程，无量纲系数采用规范无量纲系数用表查得。a*h>4.0时取a*h=4.0无量纲系数，
3 k" n( o  y0 i' k' B4 ?         无量纲系数略微小差异，计算结果x0=4.1392，其实“4.1392”与“4.1203”得差异主要还是因为a*h>4.0按a*h=4.0时照成的。
    5、结论：也就是说我完全按照规范计算过程得出得x0=4.1392与规范算例x0=4.44差异是明显的；我采用计算机数值解法解出得结果
         x0=3.913与算例提供得精确解x0=4.35相差也是明显，
+ B! X! Y# t* x$ a( m

tautau 发表于 2020-5-12 15:45
我知道你的为什么对不上，我的是严格按照算例的参数取值。

    你都没有原始数据，如果是凑数的话，我也可以完全凑到的。我可只是说可能规范不严谨。
/ q0 Z0 B0 G; J( e  n1 |    微分方程数值求解的结果理论上是约等于精确值的。我的有限元计算结果接近，你应该是模拟有问题。
" ?3 R. M# ^8 Z2 h$ J$ |: C

ntsjy2013 发表于 2020-5-12 17:22
/ e7 ], l+ Y  b* a就是对不上，midas结果小了很多M=215kNm，x=3.1mm，excel按规范计算相差无几。
" e) O  y! f+ \$ r, B3 x% cmidas边界条件为底部约束z ...

有限元计算是不需要考虑换算地基比例系数的，就像数值解法一样，可以建n层土基，得出的结果不需要修正。用换算比例系数就有问题，需要修正。

这点差别可以忽略吧。

ntsjy2013 发表于 2020-5-9 08:17
这点差别可以忽略吧。

) o+ }0 e2 `+ t' i关键规范上说精确解啊，理论上精确解至少应该非常相近。即使采用规范系数用表计算，水平位移也是不同的哦
如下图

! |5 \2 ]% D$ Q$ `+ |2 r$ v$ q

楼主这自己计算的精确解不是参照规范吗，还有两个规范解值是分别按07和19规范？？？

这个算例我也算过，excel截图是规范算例算的。另外我用midas也验算了一下，结果也是吻合的。

. l) Y7 ^8 V; o: }  u, \

tautau 发表于 2020-5-11 13:31
这个算例我也算过，excel截图是规范算例算的。另外我用midas也验算了一下，结果也是吻合的。

你的桩变形系数0.4113都对不上，不知道你计算的原始数据多少

完全按照规范计算计算过程用表格

计算表格本来向上传，结果8楼没有原始计算表格

midas计算结果不同，excel按规范计算结果接近。

loupoo 发表于 2020-5-11 15:56
" q% I9 u' a, p! [6 I6 u你的桩变形系数0.4113都对不上，不知道你计算的原始数据多少

" Q+ z( K3 ?- w5 s! n我知道你的为什么对不上，我的是严格按照算例的参数取值。
5 W6 V  \" U. m

    我用有限元软件计算的结果与计算机微分数值解法得出的结果是接近，单元划分采用0.2m。
( b+ V3 }  o9 |3 J; @    计算结果
    有限元水平位移x0=3.96mm，弯矩Mmax=343.33kN*m
) [6 {8 g1 T2 V$ Y8 o* @$ ?    数值解法水平位移x0=3.913mm，弯矩Mmax=344.962kN*m
' q0 [2 p7 n; H) S! r
& O  U  G+ e* ?+ D

ntsjy2013，midas算的挺准的啊。

tautau 发表于 2020-5-12 15:45
我知道你的为什么对不上，我的是严格按照算例的参数取值。

$ b# z( y) `. ~7 k5 c: k# ~凑数与规范一样，桩的弹性模量需要调整2.7E+7kPa
$ ^, {0 Y* x1 K2 J9 P, V6 W