冰荷载作用下海洋平台结构响应分析命令流，并请教一问题

西风独自凉... 发表于 2010-5-31 08:41:15 | 显示全部楼层 | 阅读模式
小弟对ansys不是太熟练，做完硕士课题后，感觉学好ansys还是必要的，再说离开学还有一段时间，小弟就较为系统地学了下ansys，所用参考教材为：《Algor、Ansys在桥梁工程中的应用方法与实例》，《ANSYS.结构有限元高级分析方法与范例应用（第一版）》，《ANSYS.结构有限元高级分析方法与范例应用（第二版）》，《ANSYS土木工程应用实例》，《ansys高级工程有限元分析范例精选——祝效华》。这些教材，基本上学习ansys最好的教材了，但遗憾的是这些书中都有错误。如能将这些错误改正，对自己和别人都大有益处。
这些书中的错误有些地方小弟可以自己看出来，有些地方就无能为力了，在此，特向各位师长，学长，同仁请教。
首先，将自己发现并改正的例子发出来，供大家讨论。让请大家帮忙解决一下本人解决不了的问题。
《ANSYS.结构有限元高级分析方法与范例应用（第一版）》第14章冰荷载作用下海洋平台结构响应分析模型没有加约束，造成分析错误，本人提供完整的模型命令流如下，按此命令流和书中的分析及后处理可以得到书中的分析效果。
模型如下：
FINISH
/CLEAR
/FILENAME,PLATFORM
/TITLE, PLATFORM-STUDY
/PREP7
!定义单元类型
ET,1,PIPE16
ET,2,BEAM4
ET,3,SHELL63
!定义实常数
R,1,1.2,0.03
R,2,0.8,0.02
R,3,0.5,0.02
R,4,0.3,0.02
R,5,0.06,0.0002,0.00045,0.3,0.2
R,6,0.02,0.02,0.02,0.02
!定义BEAM4单元截面
SECTYPE,1,BEAM,RECT,,0
SECOFFSET,USER,,-0.1
SECDATA,0.2,0.3,0,0,0,0,0,0,0,0
!定义材料参数
MP,EX,1,2.0E11
MP,NUXY,1,0.3
MP,DENS,1,7800
!导管架的关键点
K,1,-10,-7.5,0
K,2,0,-9.5,0
K,3,10,-7.5,0
K,4,10,7.5,0
K,5,0,9.5,0
K,6,-10,7.5,0
K,7,-10,0,0
K,8,-5.588,0,0
K,9,0,0,0
K,10,5.588,0,0
K,11,10,0,0
K,12,-9.454,-3.438,10.91
K,13,-9.454,3.438,10.91
K,14,0,-4.313,11.01
K,15,0,4.313,11.01
K,16,9.454,-3.438,10.91
K,17,9.454,3.438,10.91
K,18,-4.75,-7.75,10.53
K,19,4.75,-7.75,10.53
K,20,-4.75,7.75,10.53
K,32,4.75,7.75,10.53
K,21,-9,-6.25,20
K,22,0,-7.75,20
K,23,9,-6.25,20
K,24,9,6.25,20
K,25,0,7.75,20
K,26,-9,6.25,20
K,27,-9,0,20
K,28,-4.982,0,20
K,29,0,0,20
K,30,4.982,0,20
K,31,9,0,20
K,41,-8,-5,40
K,42,0,-6,40
K,43,8,-5,40
K,44,8,5,40
K,45,0,6,40
K,46,-8,5,40
K,47,-8,0,40
K,48,-4.364,0,40
K,49,0,0,40
K,50,4.364,0,40
K,51,8,0,40
K,61,-7,-3.75,60
K,62,0,-4.25,60
K,63,7,-3.75,60
K,64,7,3.75,60
K,65,0,4.25,60
K,66,-7,3.75,60
K,68,-3.720,0,60
K,70,3.720,0,60
K,81,-6,-2.5,80
K,82,0,-2.5,80
K,83,6,-2.5,80
K,84,6,2.5,80
K,85,0,2.5,80
K,86,-6,2.5,80
K,88,-3,0,80
K,90,3,0,80
!基础约束关键点
K,71,-10.1,-7.625,-2
K,72,0,-9.675,-2
K,73,10.1,-7.625,-2
K,74,10.1,7.625,-2
K,75,0,9.675,-2
K,76,-10.1,7.625,-2
!导管架－甲板连接处关键点
K,91,-6,-2.5,82
K,92,0,-2.5,82
K,93,6,-2.5,82
K,94,6,2.5,82
K,95,0,2.5,82
K,96,-6,2.5,82
!甲板部分关键点
K,101,-6,-2.5,90
K,102,0,-2.5,90
K,103,6,-2.5,90
K,104,6,2.5,90
K,105,0,2.5,90
K,106,-6,2.5,90
K,107,-8,-4.5,90
K,108,-6,-4.5,90
K,109,0,-4.5,90
K,110,6,-4.5,90
K,111,8,-4.5,90
K,112,-8,0,90
K,113,-6,0,90
K,114,0,0,90
K,115,6,0,90
K,116,8,0,90
K,117,-8,-2.5,90
K,118,8,-2.5,90
K,119,-8,2.5,90
K,120,8,2.5,90
K,121,-8,4.5,90
K,122,-6,4.5,90
K,123,0,4.5,90
K,124,6,4.5,90
K,125,8,4.5,90
K,151,-6,-2.5,98
K,152,0,-2.5,98
K,153,6,-2.5,98
K,154,6,2.5,98
K,155,0,2.5,98
K,156,-6,2.5,98
K,157,-9,-5.5,98
K,158,-6,-5.5,98
K,159,0,-5.5,98
K,160,6,-5.5,98
K,161,9,-5.5,98
K,162,-9,0,98
K,163,-6,0,98
K,164,0,0,98
K,165,6,0,98
K,166,9,0,98
K,167,-9,-2.5,98
K,168,9,-2.5,98
K,169,-9,2.5,98
K,170,9,2.5,98
K,171,-9,5.5,98
K,172,-6,5.5,98
K,173,0,5.5,98
K,174,6,5.5,98
K,175,9,5.5,98
!结构框架连线
L,1,2
!平面
L,2,3
L,3,11
L,11,4
L,4,5
L,5,6
L,6,7
L,7,1
L,9,2
L,9,5
L,8,1
L,8,2
L,8,9
L,8,5
L,8,6
L,8,7
L,10,2
L,10,3
L,10,11
L,10,4
L,10,5
L,10,9
L,21,22
L,22,23
L,23,31
L,31,24
L,24,25
L,25,26
L,26,27
L,27,21
L,29,22
L,29,25
L,28,21
L,28,22
L,28,29
L,28,25
L,28,26
L,28,27
L,30,22
L,30,23
L,30,31
L,30,24
L,30,25
L,30,29
L,41,42
L,42,43
L,43,51
L,51,44
L,44,45
L,45,46
L,46,47
L,47,41
L,49,42
L,49,45
L,48,41
L,48,42
L,48,49
L,48,45
L,48,46
L,50,42
L,50,43
L,50,44
L,50,45
L,50,49
L,61,62
L,62,63
L,63,64
L,64,65
L,65,66
L,66,61
L,62,65
L,68,61
L,68,62
L,68,65
L,68,66
L,70,62
L,70,63
L,70,64
L,70,65
L,81,82
L,82,83
L,83,84
L,84,85
L,85,86
L,86,81
L,82,85
L,88,81
L,88,82
L,88,85
L,88,86
L,90,82
L,90,83
L,90,84
L,90,85
L,1,21
!侧面
L,21,41
L,41,61
L,61,81
L,6,26
L,26,46
L,46,66
L,66,86
L,13,6
L,13,7
L,13,27
L,13,26
L,12,7
L,12,1
L,12,21
L,12,27
L,7,27
L,27,47
L,26,47
L,21,47
L,47,66
L,47,61
L,61,86
L,2,22
L,22,42
L,42,62
L,62,82
L,5,25
L,25,45
L,45,65
L,65,85
L,15,5
L,15,9
L,15,29
L,15,25
L,9,29
L,14,9
L,14,2
L,14,22
L,14,29
L,29,49
L,25,49
L,22,49
L,49,65
L,49,62
L,62,85
L,3,23
L,23,43
L,43,63
L,63,83
L,4,24
L,24,44
L,44,64
L,64,84
L,17,4
L,17,11
L,17,31
L,17,24
L,11,31
L,16,11
L,16,3
L,16,23
L,16,31
L,31,51
L,24,51
L,23,51
L,51,64
L,51,63
L,63,84
L,18,1
!立面
L,18,2
L,18,22
L,18,21
L,19,2
L,19,3
L,19,23
L,19,22
L,21,42
L,23,42
L,42,61
L,42,63
L,62,81
L,62,83
L,20,6
L,20,5
L,20,25
L,20,26
L,32,5
L,32,4
L,32,24
L,32,25
L,26,45
L,24,45
L,45,66
L,45,64
L,65,86
L,65,84
L,1,71
!桩基础连线
L,2,72
L,3,73
L,4,74
L,5,75
L,6,76
L,81,91
!导管架—甲板连接处连线
L,82,92
L,83,93
L,84,94
L,85,95
L,86,96
L,101,102
!甲板平面
L,102,103
L,103,115
L,115,104
L,104,105
L,105,106
L,106,113
L,113,101
L,107,108
L,108,109
L,109,110
L,110,111
L,112,113
L,113,114
L,114,115
L,115,116
L,117,101
L,103,118
L,119,106
L,104,120
L,121,122
L,122,123
L,123,124
L,124,125
L,107,117
L,117,112
L,112,119
L,119,121
L,108,101
L,106,122
L,109,102
L,102,114
L,114,105
L,105,123
L,110,103
L,104,124
L,111,118
L,118,116
L,116,120
L,120,125
L,151,152
L,152,153
L,153,165
L,165,154
L,154,155
L,155,156
L,156,163
L,163,151
L,157,158
L,158,159
L,159,160
L,160,161
L,162,163
L,163,164
L,164,165
L,165,166
L,167,151
L,153,168
L,169,156
L,154,170
L,171,172
L,172,173
L,173,174
L,174,175
L,157,167
L,167,162
L,162,169
L,169,171
L,158,151
L,156,172
L,159,152
L,152,164
L,164,155
L,155,173
L,160,153
L,154,174
L,161,168
L,168,166
L,166,170
L,170,175
L,91,101
!甲板撑杆
L,101,151
L,92,102
L,102,152
L,93,103
L,103,153
L,94,104
L,104,154
L,95,105
L,105,155
L,96,106
L,106,156
L,91,102
L,102,93
L,94,105
L,105,96
L,91,106
L,95,102
L,93,104
L,151,102
L,102,153
L,154,105
L,105,156
L,106,151
L,104,153
L,102,155
!甲板面
A,107,108,101,117
A,108,109,102,101
A,109,110,103,102
A,110,111,118,103
A,117,101,113,112
A,101,102,114,113
A,102,103,115,114
A,103,118,116,115
A,112,113,106,119
A,113,114,105,106
A,114,115,104,105
A,115,116,120,104
A,119,106,122,121
A,106,105,123,122
A,105,104,124,123
A,104,120,125,124
A,157,158,151,167
A,158,159,152,151
A,159,160,153,152
A,160,161,168,153
A,167,151,163,162
A,151,152,164,163
A,152,153,165,164
A,153,168,166,165
A,162,163,156,169
A,163,164,155,156
A,164,165,154,155
A,165,166,170,154
A,169,156,172,171
A,156,155,173,172
A,155,154,174,173
A,154,170,175,174
!框架网格划分
LSEL,S,,,95,102
LSEL,A,,,118,125
LSEL,A,,,141,148
LSEL,A,,,192,203
LATT,1,1,1
LESIZE,ALL,,,1
LMESH,ALL
LSEL,S,,,1,94
LSEL,A,,,103,117
LSEL,A,,,126,140
LSEL,A,,,149,191
LATT,1,3,1
LESIZE,ALL,,,1
LMESH,ALL
LSEL,S,,,284,295
LATT,1,2,1
LESIZE,ALL,,,1
LMESH,ALL
LSEL,S,,,296,319
LATT,1,4,1
LESIZE,ALL,,,1
LMESH,ALL
!甲板梁网格划分
LSEL,S,,,192,283
LATT,1,5,2,,,,1
LESIZE,ALL,,,3
LMESH,ALL
!甲板面网格划分
ASEL,A,,,ALL
AATT,1,6,3
AMESH,ALL
D,31,ALL
D,32,ALL
D,33,ALL
D,34,ALL
D,35,ALL
D,36,ALL
allsel,all
fini
使用道具 举报 回复
精彩评论倒序浏览

2520查看2评论

西风独自凉... 发表于 2010-5-31 08:45:46
本帖最后由 西风独自凉 于 2010-5-31 08:52 编辑

加载及后处理
/sol

antype,transient

prnopt,full

nsubst,10,0,0

time,1.9

alphad,0.05

betad,0.05

outres,all,1

F,5,FY,0

F,25,FY,0

lswrite,1,



time,2.0

F,5,FY,1e5

F,25,FY,1e5

lswrite,2,



time,3.0

F,5,FY,0

F,25,FY,0

lswrite,3,



time,10

nsubst,50,0,0

lswrite,4,

lssolve,1,4,1,

fini





/post26

nsol,2,73,u,y,uy

/axlab,x,time

/axlab,y,uy

plvar,2

fini


                               
登录/注册后可看大图




举报 回复
西风独自凉... 发表于 2010-5-31 08:53:04

《ANSYS土木工程应用实例》李围4.4节钢筋混凝土梁分析有错误，本人始终找不出来，希望大家人尽其才，帮忙看看。
!采用均匀网格划分，同时取消两种垫片,DMX=6.415,最大第一主应力263.629，28个子部计算完毕
FINISH
/CLEAR
/clear,nostart
/config,nres,5000
/prep7
!定义单元类型
et,1,solid65            !混凝土单元
et,2,pipe20                !钢筋单元
et,3,plane42                !用于拉伸成体单元之用

!定义实常数
r,1,18,8.99                !纵向受拉钢筋实常数       
r,2,8,3.99                !纵向受压钢筋，箍筋实常数
r,3                        !混凝土
                       
!定义材料属性
mp,ex,1,24000                !混凝土材料属性
mp,prxy,1,0.2               
tb,conc,1,1,9
tbdata,,0.4,1,3,-1

mp,ex,2,2e5                !纵向受拉钢筋
mp,prxy,2,0.3
tb,bkin,2,1,2,1
tbdata,,350

mp,ex,3,2e5                !横向箍筋，受压钢筋材料属性
mp,prxy,3,0.25
tb,bkin,3,1,2,1
tbdata,,200

/pnum,node,1
/pnum,elem,1

!!产生所有的节点**********!!!
n,1
n,9,150
fill,1,9
ngen,11,9,1,9,1,,30
ngen,19,1000,1,99,1,,,-75
/view,1,1,1,1       

!!!!!!!箍筋，受压钢筋*******!!!!
type,2
real,2
mat,3
!水平箍筋
*do,ii,11,16,1
  e,ii,ii+1
*enddo
*do,ii,83,88,1
  e,ii,ii+1
*enddo
!竖直箍筋
*do,ii,11,74,9
  e,ii,ii+9
*enddo
*do,ii,17,80,9
  e,ii,ii+9
*enddo
!产生整个模型的箍筋
egen,19,1000,1,28,1
!此时最大的单元编号为532

!纵向受压钢筋
*do,ii,83,17083,1000
  e,ii,ii+1000
*enddo
*do,ii,89,17089,1000
  e,ii,ii+1000
*enddo


!!!!!纵向受拉钢筋************!!!!!!!!!!
type,2
real,1
mat,2
*do,ii,11,17011,1000
  e,ii,ii+1000
*enddo
*do,ii,17,17017,1000
  e,ii,ii+1000
*enddo
!此时最大单元编号为604
/view,1,1,1,1
/pnum,elem,0
/pnum,node,0
/eshape,1
eplot

!!!**********生成面单元，以便拉伸成体单元!!!!!!

k,1
k,2,150,
k,3,150,300
k,4,,300
a,1,2,3,4

lsel,s,loc,y,0
lsel,a,loc,y,300
lesize,all,,,8
lsel,all
lsel,s,loc,x,0
lsel,a,loc,x,150
lesize,all,,,10
type,3
amesh,all

type,1                        !拉伸成混凝土单元
real,3
mat,1
extopt,esize,18
extopt,aclear,1
vext,all,,,,,-1350


/pnum,mat,1
/pnum,node,0
/pnum,elem,0
/view,1,1,1,1
eplot
allsel
nummrg,all
numcmp,all
eplot


/solu
nsel,s,loc,y,-0                !右边界约束
nsel,r,loc,z,-1275
d,all,uy
d,all,ux


asel,s,loc,z,0                !对称面约束
da,all,symm

nsel,all                !施加外部载荷
nsel,s,loc,y,300
nsel,r,loc,z,-600
d,all,uy,-5

allsel
gplot

nlgeom,on
NSUBST,200   
OUTRES,ALL,ALL  
NEQIT,50
pred,on
CNVTOL,F, ,0.05,2,0.5
allsel
SOLVE   
FINISH

/post1
allsel
/DEVICE,VECTOR,1
PLCRACK,0,1                !显示第一开裂位置
PLCRACK,0,2                !显示第二开裂位置
/DEVICE,VECTOR,0
PLNSOL,U,Y,0,1

ESEL,S,TYPE,,2
PLNSOL,S,1,0,1
举报 回复
发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

西风独自凉 实名认证

高级工程师

  • 0

    关注

  • 27

    粉丝

  • 99

    主题

TA 的主题

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