有没有人做过桥墩的稳定分析
如题！本人现在在做一个独柱墩稳定分析，桥是半径45m的曲线箱梁，固结一直径150cm的圆柱独柱墩，要分析在风载和离心力作用下的墩柱稳定性和强度，望高手给于指点
!互通桥，3孔16米现浇箱梁，双幅桥，半径分别为49.37m和42.62m位于曲线上，墩柱采用1.5m接1.8m桩基，分析桥墩的稳定性和强度；含活载偏载和风压、离心力。离心力采用实际轮子位置添加。
!分析了左幅和右幅添加了横系梁后情况
!独圆柱墩稳定性和强度分析
!******************************************************
/COM, Structural
!----------定义单元及材料等---------------------
/PREP7
et,2,shell63 !定义shell63单元
et,1,solid45 !定义solid45单元
r,2,0.1 !定义shell63的实常数号
r,1,0.1 !定义solid45的实常数号
mp,ex,1,3e10 !定义solid45单元的弹性模量
mp,prxy,1,0.3 !定义solid45单元的泊松系数
mp,dens,1,2600 !定义混凝土的密度
!定义参数*******************************************************
WH1=35 !外侧桥墩1高度
WH2=34.5 !外侧桥墩2高度
NH1=32 !内侧桥墩1高度
NH2=30 !内侧桥墩2高度
WR=49.37 !外侧设计线半径
NR=40.62 !内侧设计线半径
DUNR=0.75 !墩柱半径
ZUGR=0.9 !桩基半径
ZUGH=20 !桩长
lxlxs1=0.1954 !定义外侧离心力系数1，设计速度35KM/H
lxlxs2=0.2375 !定义内侧离心力系数2
wind=990 !定义35m处设计风压（Pa）
!***************************************************************
!----------建立几何模型-------------------------
!*******************外侧桥建模**********************************
wprota,,90,, !绕Y轴旋转工作平面
cswpla,11,0
csys,wp
!系梁关键点系梁1*1.5m
k,200001,45.9005927,16.6101671,15.0
k,200002,45.5507009,17.5469573,15.0
k,200003,38.4011505,14.8765950,15.0
k,200004,38.7510423,13.9398049,15.0
k,200101,45.9005927,16.6101671,16.5
k,200102,45.5507009,17.5469573,16.5
k,200103,38.4011505,14.8765950,16.5
k,200104,38.7510423,13.9398049,16.5
k,300001,37.1874608,31.6204811,15.0
k,300002,36.5319104,32.3756325,15.0
k,300003,30.7686201,27.3724946,15.0
k,300004,31.4241704,26.6173431,15.0
k,300101,37.1874608,31.6204811,16.5
k,300102,36.5319104,32.3756325,16.5
k,300103,30.7686201,27.3724946,16.5
k,300104,31.4241704,26.6173431,16.5
!活载加载关键节点（全局笛卡儿坐标系）
n,30001,41.284,8.6970,-1.000
n,30002,43.045,9.0680,-1.000
n,30003,49.846,10.501,-1.000
n,30004,51.607,10.872,-1.000
n,30005,40.572,11.570,-1.000
n,30006,42.303,12.064,-1.000
n,30007,48.987,13.970,-1.000
n,30008,50.718,14.463,-1.000
n,30009,38.223,17.860,-1.000
n,30010,39.854,18.622,-1.000
n,30011,46.151,21.564,-1.000
n,30012,47.781,22.326,-1.000
n,30013,33.196,26.038,-1.000
n,30014,34.613,27.149,-1.000
n,30015,40.081,31.439,-1.000
n,30016,41.498,32.549,-1.000
n,30017,31.128,28.479,-1.000
n,30018,32.456,29.694,-1.000
n,30019,37.583,34.386,-1.000
n,30020,38.911,35.601,-1.000
n,30021,21.851,36.091,-1.000
n,30022,22.783,37.630,-1.000
n,30023,26.383,43.576,-1.000
n,30024,27.315,45.115,-1.000
!
csys,wp
cswpla,11,1 !以当前工作平面为中心定义局部柱坐标
*AFUN,DEG !修改角度输入采用" 度"
k,1,WR,,,
k,2,WR,20.48,,
k,3,WR,40.96,,
k,4,WR,61.44,,
k,102,WR,20.48,WH1
k,103,WR,40.96,WH2
L,2,102,WH1/0.5 !墩柱分段0.5M一段
L,3,103,WH2/0.5
!
k,100001,NR,,, !编号比外侧大100000
k,100002,NR,20.48,,
k,100003,NR,40.96,,
k,100004,NR,61.44,,
k,100102,NR,20.48,NH1
k,100103,NR,40.96,NH2
L,100002,100102,NH1/0.5 !墩柱分段0.5M一段
L,100003,100103,NH2/0.5
!
kwpave,2
CYL4,0,0,DUNR, , , , !定义外侧第一个桥墩顶面
MSHAPE,1,2D
MSHKEY,0
amesh,1
kwpave,3
CYL4,0,0,DUNR, , , , !定义外侧第二个桥墩顶面
MSHAPE,1,2D
MSHKEY,0
amesh,2
!
kwpave,100002
CYL4,0,0,DUNR, , , , !定义内侧第一个桥墩顶面
MSHAPE,1,2D
MSHKEY,0
amesh,3
kwpave,100003
CYL4,0,0,DUNR, , , , !定义内侧第二个桥墩顶面
MSHAPE,1,2D
MSHKEY,0
amesh,4
kwpave,2
CSWPLA,11,1,1,1, !使局部坐标移到工作平面原点，采用柱坐标
vdrag,1,,,,,,1 !拖拉生成第一个墩柱网格
kwpave,3
CSWPLA,11,1,1,1, !使局部坐标移到工作平面原点，采用柱坐标
vdrag,2,,,,,,2 !拖拉生成第一个墩柱网格
kwpave,100002
CSWPLA,11,1,1,1, !使局部坐标移到工作平面原点，采用柱坐标
vdrag,3,,,,,,3 !拖拉生成第一个墩柱网格
kwpave,100003
CSWPLA,11,1,1,1, !使局部坐标移到工作平面原点，采用柱坐标
vdrag,4,,,,,,4 !拖拉生成第一个墩柱网格
!*********************************************************
!******************墩柱建模分网结束***********************
!*********************************************************
!建立箱梁截面
kwpave,1 !工作平面原点移到1号关键点
wprot,,-90,0 !旋转工作平面
CSWPLA,11,0,,, !使局部坐标移到工作平面原点,并采用笛卡儿坐标
!生成外箱梁断面外缘点
k,1001,2.120,0.000,0.000
k,1002,2.120,0.600,0.000
k,1003,4.370,0.820,0.000
k,1004,4.370,1.000,0.000
k,1005,-4.370,1.000,0.000
k,1006,-4.370,0.820,0.000
k,1007,-2.120,0.600,0.000
k,1008,-2.120,0.000,0.000
!生成外箱梁断面内部淘空点
k,2001,1.440,0.200,0.000
k,2002,1.620,0.380,0.000
k,2003,1.620,0.600,0.000
k,2004,0.870,0.780,0.000
k,2005,-0.870,0.780,0.000
k,2006,-1.620,0.600,0.000
k,2007,-1.620,0.380,0.000
k,2008,-1.440,0.200,0.000
a,1001,1002,1003,1004,1005,1006,1007,1008
a,2001,2002,2003,2004,2005,2006,2007,2008
asba,25,26,,delete,delete !通过淘空形成箱梁断面
ldiv,56,,,18 !控制网格划分精度
ldiv,60,,,10
MSHAPE,1,2D
MSHKEY,0
asel,all
esize,0.500
amesh,27
agen,2,27,,,-8.75,,,100000,0,0 !复制箱梁截面
!************************************
!外侧箱梁拖拉生成
k,9999,
kmove,9999,0,,,, !将点移到总体坐标系原点
larc,1,4,9999,WR !生成截面扫略路径线
esize,0.8,, !全联分层0.8m一个网格面
vdrag,27,,,,,,137 !划分箱梁网格,沿着137号线
!内侧箱梁拖拉生成
larc,100001,100004,9999,NR !生成截面扫略路径线222号
esize,0.8,, !全联分层0.8m一个网格面
vdrag,25,,,,,,222 !划分箱梁网格,沿着222号线
v,200001,200002,200003,200004,200101,200102,200103,200104
v,300001,300002,300003,300004,300101,300102,300103,300104
MSHAPE,1,3D
MSHKEY,0
esize,0.500
vmesh,7
vmesh,8
!******耦合系梁和墩柱*********
kwpave,2
wprota,,90
cswpla,11,1
csys,wp
wpstyl,,,,,,1 !定义工作平面坐标类型为柱坐标
nsel,,loc,z,15,16.5,0.5
nsel,r,loc,x,DUNR,DUNR+0.5
CP,100,ALL,ALL
kwpave,3
cswpla,11,1
csys,wp
nsel,,loc,z,15,16.5,0.5
nsel,r,loc,x,DUNR,DUNR+0.5
CP,200,ALL,ALL
kwpave,100002
cswpla,11,1
csys,wp
nsel,,loc,z,15,16.5,0.5
nsel,r,loc,x,DUNR,DUNR+0.5
CP,300,ALL,ALL
kwpave,100003
cswpla,11,1
csys,wp
nsel,,loc,z,15,16.5,0.5
nsel,r,loc,x,DUNR,DUNR+0.5
CP,400,ALL,ALL
ALLSEL,ALL
!箱梁和墩柱交界处节点耦合
kwpave,2
NSEL,,LOC,Z,0
NSEL,R,LOC,X,0,DUNR
CP,500,ALL,ALL
kwpave,3
NSEL,,LOC,Z,0
NSEL,R,LOC,X,0,DUNR
CP,600,ALL,ALL
kwpave,100002
NSEL,,LOC,Z,0
NSEL,R,LOC,X,0,DUNR
CP,700,ALL,ALL
kwpave,100003
NSEL,,LOC,Z,0
NSEL,R,LOC,X,0,DUNR
CP,800,ALL,ALL
ALLSEL,ALL
!施加桥墩约束(墩柱底部固结)
cswpla,11,1
csys,wp
kwpave,2
nsel,,loc,z,WH1
nsel,r,loc,x,0,DUNR
D,ALL,ALL
kwpave,3
nsel,,loc,z,WH2
nsel,r,loc,x,0,DUNR
D,ALL,ALL
kwpave,100002
nsel,,loc,z,NH1
nsel,r,loc,x,0,DUNR
D,ALL,ALL
kwpave,100003
nsel,,loc,z,NH2
nsel,r,loc,x,0,DUNR
D,ALL,ALL
!桥台约束（竖向约束）
KWPAVE,9999
NSEL,S,LOC,Z,0
NSEL,R,LOC,Y,0
D,ALL,UY,0
NSEL,S,LOC,Z,0
NSEL,R,LOC,Y,61.44
D,ALL,UY,0
ALLSEL,ALL
!
!/SOLU
ACEL,0,9.8,0, !重力
!桥面铺装、护栏和风压
cswpla,11,1
csys,wp
asel,s,loc,z,-1
sfa,all,1,pres,2860 !施加桥面铺装2860和护栏（1263.6）
ALLSEL,ALL
SFA,93,1,PRES,wind*2 !施加风压（考虑护栏迎风面）
!wpstyle,0.2,0.5,-1,60,0.1,1,1,,0.2 !设置格栅
!施加活载和离心力********************************************************
cswpla,11,1
csys,11,1
wprota,11.9
csys,wp
!nrotat,30001,30004
!nsel,s,node,,30001,30004
nrotat,49750
F,49750,FX,15000*lxlxs2
F,49750,FZ,15000
nrotat,49534
F,49534,FX,15000*lxlxs2
F,49534,FZ,15000
nrotat,45256
F,45256,FX,15000*lxlxs1
F,45256,FZ,15000
nrotat,44992
F,44992,FX,15000*lxlxs1
F,44992,FZ,15000
wprota,-11.9
wprota,15.92
!nsel,s,node,,30005,30008,1
!nrotat,30005,30008,1
nrotat,49753
F,49753,FX,120000*lxlxs2
F,49753,FZ,120000
nrotat,49537
F,49537,FX,120000*lxlxs2
F,49537,FZ,120000
nrotat,45194
F,45194,FX,120000*lxlxs1
F,45194,FZ,120000
nrotat,44996
F,44996,FX,120000*lxlxs1
F,44996,FZ,120000
wprota,-15.92
wprota,25.04
!nsel,s,node,,30009,30012,1
!nrotat,30009,30012,1
nrotat,49762
F,49762,FX,140000*lxlxs2
F,49762,FZ,140000
nrotat,49546
F,49546,FX,140000*lxlxs2
F,49546,FZ,140000
nrotat,45204
F,45204,FX,140000*lxlxs1
F,45204,FZ,140000
nrotat,45006
F,45006,FX,140000*lxlxs1
F,45006,FZ,140000
wprota,-25.04
wprota,38.11
!nsel,s,node,,30013,30016,1
!nrotat,30013,30016,1
nrotat,49773
F,49773,FX,35000*lxlxs2
F,49773,FZ,35000
nrotat,49557
F,49557,FX,35000*lxlxs2
F,49557,FZ,35000
nrotat,45285
F,45285,FX,35000*lxlxs1
F,45285,FZ,35000
nrotat,45021
F,45021,FX,35000*lxlxs1
F,45021,FZ,35000
wprota,-38.11
wprota,42.46
!nsel,s,node,,30017,30020,1
!nrotat,30017,30020,1
nrotat,49777
F,49777,FX,65000*lxlxs2
F,49777,FZ,65000
nrotat,49561
F,49561,FX,65000*lxlxs2
F,49561,FZ,65000
nrotat,45289
F,45289,FX,65000*lxlxs1
F,45289,FZ,65000
nrotat,45025
F,45025,FX,65000*lxlxs1
F,45025,FZ,65000
wprota,-42.46
wprota,58.81
!nsel,s,node,,30021,30024,1
!nrotat,30021,30024,1
nrotat,49792
F,49792,FX,35000*lxlxs2
F,49792,FZ,35000
nrotat,49576
F,49576,FX,35000*lxlxs2
F,49576,FZ,35000
nrotat,45307
F,45307,FX,35000*lxlxs1
F,45307,FZ,35000
nrotat,45043
F,45043,FX,35000*lxlxs1
F,45043,FZ,35000
wprota,-58.81
allsel,all
!***************************
由于采用ANSYS9建模，而且直接选节点号等，可能在不同系统上运行不了，望见谅 最近经常做些稳定方面的分析,不过大多都是一类的分支点失稳,极值点失稳倒是可以共同研究,共同进步.
你可以先看看
http://www.bridgehead.com.cn/viewthread.php?tid=197&extra=page%3D1关于稳定的大讨论 贴两篇关于稳定方面的论文----转贴至okok.org
115_大跨度桥梁静风稳定性分析及其在ANSYS中的实现.rar 102_4000t全回转浮吊下部回转支撑结构稳定性分析.rar
非常感谢版主
你发的内容正在研究中 对于高墩的悬臂施工，稳定性是一个比较大的问题！有没有一种通过现场监控的方法来进行施工过程中的安全控制呢？ 谢谢楼主，太感动了。向你多多学习了。
呵呵 学习
学习中 大家共同进步
页:
[1]