急寻：钢管混凝土拱桥的建模命令流！！！

楼上提醒的是啊！
只是做的太粗糙，希望大家多给些建议！
/prep7
/title,The analyse of steel arch bridge--chenglin 2006.06.25

!*********************************************************************************
! 下面定义材料特性
!*********************************************************************************
et,1,beam188 ! 钢管特性
mp,ex,1,2.1e11
mp,dens,1,7800
mp,prxy,1,0.3
!n,90000,0,0,30 ! 参考点

et,2,beam188 ! 拱肋内混凝土特性 50#混凝土
mp,ex,2,3.5e10
mp,dens,2,2600
mp,prxy,2,0.1667

et,3,beam4 ! 纵梁特性 30#混凝土
mp,ex,3,3.0e10
mp,dens,3,2600
mp,prxy,3,0.1667

et,4,beam4 ! 横梁特性 30#混凝土
mp,ex,4,3.0e10
mp,dens,4,2600
mp,prxy,4,0.1667

et,5,beam188 ! 风撑特性
mp,ex,5,2.1e11
mp,dens,5,7800
mp,prxy,5,0.3

et,6,link10 ! 吊杆特性 （钢绞线）
mp,ex,6,1.90e11
mp,dens,6,7800
mp,prxy,6,0.3
keyopt,6,3,0 ! 只拉吊杆

et,7,beam4 ! 端横梁特性 30#混凝土
mp,ex,7,3.0e10
mp,dens,7,2600
mp,prxy,7,0.1667

et,8,link10 !系杆特性（钢绞线）
mp,ex,8,2.0e11
mp,dens,8,7800
mp,prxy,8,0.3
keyopt,8,3,0 !只拉系杆

!**********************************************************************
! 下面是建立模型
!**********************************************************************
! 以下将主梁划分为14个单元，各单元尺寸为：0.3,0.35,0.45,……，0.45，0.35，0.3。
n,1
n,2,0.30
n,3,0.65
n,13,5.15
n,14,5.5
n,15,5.8
fill,3,13,9 !表示在前片梁3～13号节点间插入9个节点

n,16,0,0,-0.94
n,17,0.30,0,-0.94
n,18,0.65,0,-0.94
n,28,5.15,0,-0.94
n,29,5.5,0,-0.94
n,30,5.8,0,-0.94
fill,18,28,9 !表示在后片梁18～28号节点间插入9个节点

type,3
mat,3
r,3,0.032,1.067e-4,6.827e-5,0.16,0.2,,
real,3
*do,i,1,14
e,i,i+1
e,i+15,i+16
*enddo
/view,1,1,1,1
gplot
!**********************************************
!0# 横梁(左右端横梁)
type,7
r,7,0.22*0.2,0.2*0.2*0.2*0.22/12,0.22*0.22*0.22*0.2/12,0.22,0.20,,
real,7
mat,7
e,1,16
e,15,30
!**********************************************
!1～13# 横梁(端横梁)
type,4
mat,4
r,11,0.0128,2.7307e-5,6.827-6,0.08,0.16,,
real,11
*do,i,2,14
e,i,i+15
*enddo
!**********************************************
! 吊杆
*do,i,3,13
n,100+i,nx(i),(-0.13793)*(nx(i)-2.9)*(nx(i)-2.9)+1.16,0 ! 前吊杆节点
n,200+i,nx(i),(-0.13793)*(nx(i)-2.9)*(nx(i)-2.9)+1.16,-0.94 ! 后吊杆节点
*enddo
type,6
mat,6
r,6,0.016*0.016*3.14/4
real,6
!SECTYPE,6, BEAM, CSOLID, , 0 !定义截面形状
!SECOFFSET, CENT
!SECDATA,0.008,8,8 !赋予截面形状参数
!secnum,6 !调用截面
*do,i,3,13
e,100+i,i
e,200+i,i+15
*enddo
!*****************************************************************************************
! 风撑
type,5
mat,5
!r,5,0.001205,0.139e-5,0.139e-5,0.1,0.1,,
!real,5
sectype,1,beam,ctube,,1
secoffset,cent
secdata,0.0460,0.05,16
e,106,206
e,108,208
e,110,210

!**********************************************
! 拱肋建模
type,2
mat,2
!r,2,0.017658,0.248e-4,0.248e-4,0.15,0.15,,
!real,2 !生成前拱肋单元
sectype,2,beam,csolid,,1
secoffset,cent
secdata,0.075,0.0705,16

e,1,103
e,113,15
e,16,203
e,213,30
*do,i,3,12
e,i+100,i+101
e,i+200,i+201
*enddo
!**********************************************
/eshape,1
/rep,fast
eplot
!**********************************************
/solu
antype,2
modopt,lanb,15
mxpand,15,,,0
modopt,lanb,15,0,100000000,,off
d,1,ux,0
d,1,uy,0
d,16,ux,0
d,16,uy,0
d,15,uy
d,30,uy
solve
finish

/post1
set,list
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
以上命令流只是建立了一个大致的体形，还有以下的问题没有解决好：
1）钢管拱的模拟。钢管拱中填充有混凝土，在做模态分析的时候选用什么单元来模拟？是应该用梁单元还是实体单元？
2)吊杆和钢管拱之间的节点联结怎么处理？是否应该耦合节点坐标？
3）这是一个混凝土实体模型拱桥，混凝土梁该采用梁单元还是实体单元？
4）桥两端是简支撑于墩台上，由于支撑接触面相对于构件或单元尺寸较大，在有限元模拟时该如何约束？
5）纵梁中的预应力如何施加？（上命令中还没有考虑）施工中考虑的是在梁中心预留孔道待混凝土强度达到后进行后张，那么实际建模是采用分离法还是降温法？
6）桥面板的处理，个人感觉由于模型的尺寸效应，板（厚度5cm）对整个结构的刚度影响是非常大的，如何模拟较好？
……

感觉有限元模拟实际模型要考虑的因素很多，如何能尽可能的准确建模符合实际受力还是很不搞的。
希望做过类似的同仁不吝赐教。