冰荷载作用下海洋平台结构响应分析命令流，并请教一问题
小弟对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,BEAM4ET,3,SHELL63 !定义实常数R,1,1.2,0.03R,2,0.8,0.02R,3,0.5,0.02R,4,0.3,0.02R,5,0.06,0.0002,0.00045,0.3,0.2R,6,0.02,0.02,0.02,0.02 !定义BEAM4单元截面SECTYPE,1,BEAM,RECT,,0SECOFFSET,USER,,-0.1 SECDATA,0.2,0.3,0,0,0,0,0,0,0,0 !定义材料参数MP,EX,1,2.0E11MP,NUXY,1,0.3MP,DENS,1,7800 !导管架的关键点K,1,-10,-7.5,0K,2,0,-9.5,0K,3,10,-7.5,0K,4,10,7.5,0K,5,0,9.5,0K,6,-10,7.5,0K,7,-10,0,0K,8,-5.588,0,0K,9,0,0,0K,10,5.588,0,0K,11,10,0,0K,12,-9.454,-3.438,10.91K,13,-9.454,3.438,10.91K,14,0,-4.313,11.01K,15,0,4.313,11.01K,16,9.454,-3.438,10.91K,17,9.454,3.438,10.91K,18,-4.75,-7.75,10.53K,19,4.75,-7.75,10.53K,20,-4.75,7.75,10.53K,32,4.75,7.75,10.53 K,21,-9,-6.25,20K,22,0,-7.75,20K,23,9,-6.25,20K,24,9,6.25,20K,25,0,7.75,20K,26,-9,6.25,20K,27,-9,0,20K,28,-4.982,0,20K,29,0,0,20K,30,4.982,0,20K,31,9,0,20 K,41,-8,-5,40K,42,0,-6,40K,43,8,-5,40K,44,8,5,40K,45,0,6,40K,46,-8,5,40K,47,-8,0,40K,48,-4.364,0,40K,49,0,0,40K,50,4.364,0,40K,51,8,0,40 K,61,-7,-3.75,60K,62,0,-4.25,60K,63,7,-3.75,60K,64,7,3.75,60K,65,0,4.25,60K,66,-7,3.75,60K,68,-3.720,0,60K,70,3.720,0,60 K,81,-6,-2.5,80K,82,0,-2.5,80K,83,6,-2.5,80K,84,6,2.5,80K,85,0,2.5,80K,86,-6,2.5,80K,88,-3,0,80K,90,3,0,80 !基础约束关键点K,71,-10.1,-7.625,-2K,72,0,-9.675,-2K,73,10.1,-7.625,-2K,74,10.1,7.625,-2K,75,0,9.675,-2K,76,-10.1,7.625,-2 !导管架－甲板连接处关键点K,91,-6,-2.5,82K,92,0,-2.5,82K,93,6,-2.5,82K,94,6,2.5,82K,95,0,2.5,82K,96,-6,2.5,82 !甲板部分关键点K,101,-6,-2.5,90K,102,0,-2.5,90K,103,6,-2.5,90K,104,6,2.5,90K,105,0,2.5,90K,106,-6,2.5,90K,107,-8,-4.5,90K,108,-6,-4.5,90K,109,0,-4.5,90K,110,6,-4.5,90K,111,8,-4.5,90K,112,-8,0,90K,113,-6,0,90K,114,0,0,90K,115,6,0,90K,116,8,0,90K,117,-8,-2.5,90K,118,8,-2.5,90K,119,-8,2.5,90K,120,8,2.5,90K,121,-8,4.5,90K,122,-6,4.5,90K,123,0,4.5,90K,124,6,4.5,90K,125,8,4.5,90 K,151,-6,-2.5,98K,152,0,-2.5,98K,153,6,-2.5,98K,154,6,2.5,98K,155,0,2.5,98K,156,-6,2.5,98K,157,-9,-5.5,98K,158,-6,-5.5,98K,159,0,-5.5,98K,160,6,-5.5,98K,161,9,-5.5,98K,162,-9,0,98K,163,-6,0,98K,164,0,0,98K,165,6,0,98K,166,9,0,98K,167,-9,-2.5,98K,168,9,-2.5,98K,169,-9,2.5,98K,170,9,2.5,98K,171,-9,5.5,98K,172,-6,5.5,98K,173,0,5.5,98K,174,6,5.5,98K,175,9,5.5,98 !结构框架连线L,1,2
!平面L,2,3L,3,11L,11,4L,4,5L,5,6L,6,7L,7,1L,9,2L,9,5L,8,1L,8,2L,8,9L,8,5L,8,6L,8,7L,10,2L,10,3L,10,11L,10,4L,10,5L,10,9 L,21,22L,22,23L,23,31L,31,24L,24,25L,25,26L,26,27L,27,21L,29,22L,29,25L,28,21L,28,22L,28,29L,28,25L,28,26L,28,27L,30,22L,30,23L,30,31L,30,24L,30,25L,30,29 L,41,42L,42,43L,43,51L,51,44L,44,45L,45,46L,46,47L,47,41L,49,42L,49,45L,48,41L,48,42L,48,49L,48,45L,48,46L,50,42L,50,43L,50,44L,50,45L,50,49 L,61,62L,62,63L,63,64L,64,65L,65,66L,66,61L,62,65L,68,61L,68,62L,68,65L,68,66L,70,62L,70,63L,70,64L,70,65 L,81,82L,82,83L,83,84L,84,85L,85,86L,86,81L,82,85L,88,81L,88,82L,88,85L,88,86L,90,82L,90,83L,90,84L,90,85 L,1,21
!侧面L,21,41L,41,61L,61,81L,6,26L,26,46L,46,66L,66,86L,13,6L,13,7L,13,27L,13,26L,12,7L,12,1L,12,21L,12,27L,7,27L,27,47L,26,47L,21,47L,47,66L,47,61L,61,86 L,2,22L,22,42L,42,62L,62,82L,5,25L,25,45L,45,65L,65,85L,15,5L,15,9L,15,29L,15,25L,9,29L,14,9L,14,2L,14,22L,14,29L,29,49L,25,49L,22,49L,49,65L,49,62L,62,85 L,3,23L,23,43L,43,63L,63,83L,4,24L,24,44L,44,64L,64,84L,17,4L,17,11L,17,31L,17,24L,11,31L,16,11L,16,3L,16,23L,16,31L,31,51L,24,51L,23,51L,51,64L,51,63L,63,84 L,18,1
!立面L,18,2L,18,22L,18,21L,19,2L,19,3L,19,23L,19,22L,21,42L,23,42L,42,61L,42,63L,62,81L,62,83 L,20,6L,20,5L,20,25L,20,26L,32,5L,32,4L,32,24L,32,25L,26,45L,24,45L,45,66L,45,64L,65,86L,65,84 L,1,71
!桩基础连线L,2,72L,3,73L,4,74L,5,75L,6,76 L,81,91
!导管架—甲板连接处连线L,82,92L,83,93L,84,94L,85,95L,86,96 L,101,102
!甲板平面L,102,103L,103,115L,115,104L,104,105L,105,106L,106,113L,113,101L,107,108L,108,109L,109,110L,110,111L,112,113L,113,114L,114,115L,115,116L,117,101L,103,118L,119,106L,104,120L,121,122L,122,123L,123,124L,124,125L,107,117L,117,112L,112,119L,119,121L,108,101L,106,122L,109,102L,102,114L,114,105L,105,123L,110,103L,104,124L,111,118L,118,116L,116,120L,120,125 L,151,152L,152,153L,153,165L,165,154L,154,155L,155,156L,156,163L,163,151L,157,158L,158,159L,159,160L,160,161L,162,163L,163,164L,164,165L,165,166L,167,151L,153,168L,169,156L,154,170L,171,172L,172,173L,173,174L,174,175L,157,167L,167,162L,162,169L,169,171L,158,151L,156,172L,159,152L,152,164L,164,155L,155,173L,160,153L,154,174L,161,168L,168,166L,166,170L,170,175 L,91,101
!甲板撑杆L,101,151L,92,102L,102,152L,93,103L,103,153L,94,104L,104,154L,95,105L,105,155L,96,106L,106,156 L,91,102L,102,93L,94,105L,105,96L,91,106L,95,102L,93,104 L,151,102L,102,153L,154,105L,105,156L,106,151L,104,153L,102,155 !甲板面A,107,108,101,117A,108,109,102,101A,109,110,103,102A,110,111,118,103A,117,101,113,112A,101,102,114,113A,102,103,115,114A,103,118,116,115A,112,113,106,119A,113,114,105,106A,114,115,104,105A,115,116,120,104A,119,106,122,121A,106,105,123,122A,105,104,124,123A,104,120,125,124 A,157,158,151,167A,158,159,152,151A,159,160,153,152A,160,161,168,153A,167,151,163,162A,151,152,164,163A,152,153,165,164A,153,168,166,165A,162,163,156,169A,163,164,155,156A,164,165,154,155A,165,166,170,154A,169,156,172,171A,156,155,173,172A,155,154,174,173A,154,170,175,174 !框架网格划分LSEL,S,,,95,102LSEL,A,,,118,125LSEL,A,,,141,148LSEL,A,,,192,203LATT,1,1,1LESIZE,ALL,,,1LMESH,ALL LSEL,S,,,1,94LSEL,A,,,103,117LSEL,A,,,126,140LSEL,A,,,149,191LATT,1,3,1LESIZE,ALL,,,1LMESH,ALL LSEL,S,,,284,295LATT,1,2,1LESIZE,ALL,,,1LMESH,ALL LSEL,S,,,296,319LATT,1,4,1LESIZE,ALL,,,1LMESH,ALL !甲板梁网格划分LSEL,S,,,192,283LATT,1,5,2,,,,1LESIZE,ALL,,,3LMESH,ALL !甲板面网格划分ASEL,A,,,ALLAATT,1,6,3AMESH,ALL D,31,ALLD,32,ALLD,33,ALLD,34,ALLD,35,ALLD,36,ALLallsel,allfini 本帖最后由 西风独自凉 于 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
http://www.bridgehead.com.cn/data/attachment/album/201005/31/084853v6vc6q63x8vvk0vy.jpg
《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
页:
[1]