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

loupoo 发表于 2020-5-8 11:22:10 | 显示全部楼层 | 阅读模式
    1、由于实行了新基础规范，原07规范附录P调整为了附录L，桩基础效应计算原理也未变。# ^. L5 t; t+ D
    2、由于新规范的发行，在阅读附录L条文说明（0.7规范一样）P185-186页时发现个算例。
. D8 r8 f/ y" c0 j$ Z4 I    3、问题所在：本人通过计算发现精确解x0=3.913mm，Mmax=344.962kN*m，根据规范解x0=4.1203mm，Mmax=354.43kN*m2 V9 {) G1 n) h& D
         均与提供的算例数值不符，算例精确解x0=4.35mm，Mmax=335.49kN*m，本规范解x0=4.44mm，Mmax=340.31kN*m。1 A9 d$ ]0 P* X/ r- ^& K7 K& n
         如下图：- m. o  S! K1 B1 @+ j. G
         
更多图片 小图 大图
组图打开中，请稍候......
使用道具 举报 回复
精彩评论倒序浏览

2.6万查看21评论

loupoo 发表于 2020-5-11 10:44:45
雪狼儿88 发表于 2020-5-11 09:35" T. g" |& z! j6 H! x
楼主这自己计算的精确解不是参照规范吗，还有两个规范解值是分别按07和19规范？？？

( e9 |6 ~1 X2 s8 }7 h$ W/ a; v    1、我计算的精确解是采用微分方程数值解法直接用电脑解微分方程计算出来的。) T7 J4 ^1 B! d3 L% Y( I* S( S
    2、“本规范解x0=4.44”是规范原话（07、19规范都是一样的，本附录未变动）。. m/ c- y: R8 N
    3、“根据规范解x0=4.1203”是我采用规范附录L（或附录P）计算过程，但是无量纲系数采用的幂级数通解递推公式得出。
. B, j6 w* ^8 {; `; F. }    4、3楼是我完全按照规范计算过程，无量纲系数采用规范无量纲系数用表查得。a*h>4.0时取a*h=4.0无量纲系数，
6 {7 V' S7 P- W' m' R! r0 N: L         无量纲系数略微小差异，计算结果x0=4.1392，其实“4.1392”与“4.1203”得差异主要还是因为a*h>4.0按a*h=4.0时照成的。$ U0 ^6 ~3 w0 B+ X& G
    5、结论：也就是说我完全按照规范计算过程得出得x0=4.1392与规范算例x0=4.44差异是明显的；我采用计算机数值解法解出得结果1 s3 u" v" E# v- [
         x0=3.913与算例提供得精确解x0=4.35相差也是明显，! E7 c) e5 |  A. G6 |, ]
举报 回复
loupoo 发表于 2020-5-12 16:09:55
tautau 发表于 2020-5-12 15:45
6 u/ I3 h# v2 u8 z" \$ Q我知道你的为什么对不上，我的是严格按照算例的参数取值。

+ n) q* s" b: a( Y7 {% Q( o/ S$ L4 [    你都没有原始数据，如果是凑数的话，我也可以完全凑到的。我可只是说可能规范不严谨。8 h& C/ w4 [$ k" G+ B: B
    微分方程数值求解的结果理论上是约等于精确值的。我的有限元计算结果接近，你应该是模拟有问题。
1 `4 T1 M2 V; n+ U, k" D0 a, O$ v, d& o5 s; {9 o6 E2 X; e
截图202005121607317591..png
举报 回复
loupoo 发表于 2020-5-13 08:51:23
ntsjy2013 发表于 2020-5-12 17:22
( ]& z+ p0 ^0 K1 l" z就是对不上，midas结果小了很多M=215kNm，x=3.1mm，excel按规范计算相差无几。
$ @! E7 \4 h5 E# m9 d8 Q3 N. Tmidas边界条件为底部约束z ...
( ~8 P+ H4 a0 y
有限元计算是不需要考虑换算地基比例系数的，就像数值解法一样，可以建n层土基，得出的结果不需要修正。用换算比例系数就有问题，需要修正。
; k. y( E3 q1 V1 m8 n: p4 z
举报 回复
ntsjy2013 发表于 2020-5-9 08:17:20
这点差别可以忽略吧。

点评

loupoo
关键规范上说精确解啊，理论上精确解至少应该非常相近。即使采用规范系数用表计算，水平位移也是不同的哦 如下图  详情 回复 发表于 2020-5-9 16:56
举报 回复
loupoo 发表于 2020-5-9 16:56:37
ntsjy2013 发表于 2020-5-9 08:17
3 y3 y5 r" @& Q) I* H9 g7 J这点差别可以忽略吧。

& j6 {' ~& |1 E7 v3 ?9 h/ r关键规范上说精确解啊，理论上精确解至少应该非常相近。即使采用规范系数用表计算，水平位移也是不同的哦. U# C3 [- o4 D
如下图1 u% S2 V! B. ^) s' ^2 e5 a
+ T. [9 \! F; q# p) h2 r0 l' }
截图202005091656297438..png
举报 回复
雪狼儿88... 发表于 2020-5-11 09:35:34
楼主这自己计算的精确解不是参照规范吗，还有两个规范解值是分别按07和19规范？？？
1111.png

点评

loupoo
1、我计算的精确解是采用微分方程数值解法直接用电脑解微分方程计算出来的。 2、“本规范解x0=4.44”是规范原话（07、19规范都是一样的，本附录未变动）。 3、“根据规范解x0=4.1203”是我采用规范附  详情 回复 发表于 2020-5-11 10:44
举报 回复
tautau 发表于 2020-5-11 13:31:50
这个算例我也算过，excel截图是规范算例算的。另外我用midas也验算了一下，结果也是吻合的。
更多图片 小图 大图
组图打开中，请稍候......

点评

loupoo
你的桩变形系数0.4113都对不上，不知道你计算的原始数据多少  详情 回复 发表于 2020-5-11 15:56
举报 回复
loupoo 发表于 2020-5-11 15:56:04
本帖最后由 loupoo 于 2020-5-11 15:59 编辑 3 u. W# u, G6 D: e& j
, B" w6 ]" d7 P% m
, ^! t% _. x5 N( K
截图202005111552355560..png tautau 发表于 2020-5-11 13:31
3 [- q) O1 Q: M, S2 r/ v1 z/ y这个算例我也算过，excel截图是规范算例算的。另外我用midas也验算了一下，结果也是吻合的。
. {2 a, v( z8 f1 @2 E) e* d/ p- z4 Y
你的桩变形系数0.4113都对不上，不知道你计算的原始数据多少. j: {2 M. L; d# p/ J/ c( ?' L5 W

点评

tautau
我知道你的为什么对不上，我的是严格按照算例的参数取值。  详情 回复 发表于 2020-5-12 15:45
举报 回复
loupoo 发表于 2020-5-11 16:02:43
完全按照规范计算计算过程用表格
举报 回复
loupoo 发表于 2020-5-11 16:38:51
计算表格本来向上传，结果8楼没有原始计算表格
举报 回复
ntsjy2013 发表于 2020-5-12 15:02:26
midas计算结果不同，excel按规范计算结果接近。
举报 回复
tautau 发表于 2020-5-12 15:45:03
loupoo 发表于 2020-5-11 15:56
- K# y. b0 I  w你的桩变形系数0.4113都对不上，不知道你计算的原始数据多少
) B5 `  F  \4 R6 r
我知道你的为什么对不上，我的是严格按照算例的参数取值。1 E3 P) B0 ]; x' p

点评

loupoo
凑数与规范一样，桩的弹性模量需要调整2.7E+7kPa  详情 回复 发表于 2020-5-12 16:24
举报 回复
loupoo 发表于 2020-5-12 15:51:33
    我用有限元软件计算的结果与计算机微分数值解法得出的结果是接近，单元划分采用0.2m。
- z& Z0 h& X% r    计算结果
; p% e4 ?. }+ n0 Z' Q' i( F    有限元水平位移x0=3.96mm，弯矩Mmax=343.33kN*m6 o1 u7 c; X. h  l( v
    数值解法水平位移x0=3.913mm，弯矩Mmax=344.962kN*m
. G- j# o4 c  s- i% T0 Z7 ^* I5 M5 y* t7 g
截图202005121546582959..png
截图202005121547411035..png
举报 回复
tautau 发表于 2020-5-12 15:51:53
ntsjy2013，midas算的挺准的啊。
截图202005121548587360..png
举报 回复
loupoo 发表于 2020-5-12 16:24:30
tautau 发表于 2020-5-12 15:45
- W7 D  V0 S! U4 x; d我知道你的为什么对不上，我的是严格按照算例的参数取值。

8 ]0 x1 n# @! P+ q  L凑数与规范一样，桩的弹性模量需要调整2.7E+7kPa8 B7 R$ J3 T' O$ E' D; z
截图202005121622593694..png
举报 回复
12下一页
发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

loupoo 实名认证实名认证

高级工程师

  • 0

    关注

  • 6

    粉丝

  • 8

    主题

TA 的主题

 
 
  • QQ：56984982
    		•
    点击这里给我发消息
    电话：13527553862
    站务咨询群桥头堡站务咨询桥梁专业交流群：
    中国桥梁专业领袖群
    工作时间
    8:00-18:00
    网站地图小黑屋手机版Archiver桥梁论坛专业领导者 渝ICP备11004164号-1|重庆公安备案:2012004611
    Copyright © 2001-2021 Qiaotoubao Inc.