用188单元建的桁架模型,在进入后处理的时候,通过建立单元表提取的列表显示的最大等效应力值,跟云图上显示的最大等效应力值
问题:用188单元建的桁架模型,在进入后处理的时候,通过建立单元表提取的列表显示的最大等效应力值,跟云图上显示的最大等效应力值不一样???
命令流如下:
Finish
/Clear,all
/UNITS,SI
/Filename,beamopt
/PreP7
ET,1, BEAM188
ET,2,MASS21
KEYOPT,2,3,2
R,1,3664
R,2,784.2325
MP,EX,1,2.06e11
MP,PRXY,1,0.3
MP,DENS,l,7850
SECTYPE,1,BEAM,I
SECOFFSET, CENT
SECDATA,0.25,0.25,0.34,0.014,0.014,0.009
SECTYPE,2,BEAM,I
SECOFFSET, CENT
SECDATA,0.25,0.25,0.25,0.014,0.014,0.009
SECTYPE,3,BEAM,L
SECOFFSET, CENT
SECDATA,0.063,0.1,0.007,0.007
SECTYPE,4,BEAM,CTUBE
SECOFFSET, CENT
SECDATA,0.067,0.073
SECTYPE,5,BEAM,HREC
SECOFFSET, CENT
SECDATA,0.1,0.1,0.006,0.006,0.006,0.006
SECTYPE,6,BEAM,I
SECOFFSET, CENT
SECDATA,0.138,0.138,0.36,0.0158,0.0158,0.012
CSYS,0
X=86!井架倾斜角变量
*AFUN,DEG !设置角度的输入输出单位为度
K24X=-4.555+47.500/TAN(X)
K24Y=47.500
K,l,-4.555,0,0
K,24,K24X,K24Y
KBETW,1,24,4,DIST,8.275
KBETW,4,24,5,DIST,1.800
KBETW,5,24,6,DIST,0.625
KBETW,6,24,10,DIST,9.600
KBETW,10,24,11,DIST,0.450
KBETW,11,24,15,DIST,9.600
KBETW,15,24,16,DIST,0.450
KBETW,16,24,20,DIST,9.600
KBETW,20,24,21,DIST,0.450
KFILL,1,4
KFILL,6,10
KFILL,11,15
KFILL,16,20
KFILL,21,24
K27X=-4.555+8.275*COS(X)
K27Y=8.275*SIN(X)
K27Z=-2.423
K47X=K24X
K47Y=K24Y
K47Z=-2.050
K,27,K27X,K27Y,K27Z
K,47,K47X,K47Y,K47Z
KBETW,27,47,28,DIST,1.800
KBETW,28,47,29,DIST,0.625
KBETW,29,47,33,DIST,9.600
KBETW,33,47,34,DIST,0.450
KBETW,34,47,38,DIST,9.600
KBETW,38,47,39,DIST,0.450
KBETW,39,47,43,DIST,9.600
KBETW,43,47,44,DIST,0.450
KFILL,1,27,2,25,1
KFILL,29,33
KFILL,34,38
KfILL,39,43
KFILL,44,47
*Do,I,1,23,1
L,I,I+1
*ENDDO
*DO,I,25,46,1
L,I,I+1
*ENDDO
L,1,25
*Do,I,2,24,1
L,I,I+23
*ENDDO
L,3,25
L,4,26
L,7,29
L,7,31
L,31,9
L,9,33
L,12,34
L,12,36
L,14,38
L,14,36
L,17,39
L,17,41
L,19,41
L,19,43
L,22,44
L,22,46
L,24,46
LSYMM,X,ALL
KBETW,47,94,95,,0.2
KBETW,47,94,96,,0.8
KBETW,46,93,97,,0.153
KBETW,46,93,98,,0.847
KBETW,37,84,99
KBETW,35,82,100
KBETW,32,79,101
KBETW,30,77,102
KBETW,28,75,,,0.0656
KBETW,28,75,,,0.36
KBETW,28,75,,,0.64
KBETW,28,75,,,0.9344
KBETW,27,74,,,0.0656
KBETW,27,74,,,0.36
KBETW,27,74,,,0.64
KBETW,27,74,,,0.9344
L,47,95
L,95,96
L,96,94
L,46,97
L,97,98
L,98,93
L,97,96
L,95,98
L,45,92
L,42,89
L,40,87
L,37,99
L,99,84
L,35,100
L,100,82
L,32,101
L,101,79
L,30,102
L,102,77
L,28,103
L,103,104
L,104,105
L,105,106
L,106,75
L,27,107
L,107,108
L,108,109
L,109,110
L,110,74
L,103,107
L,107,104
L,104,108
L,104,109
L,109,105
L,109,106
L,106,110
L,92,42
L,42,87
L,99,40
L,99,87
L,100,37
L,100,84
L,101,35
L,101,82
L,102,32
L,102,79
L,30,104
L,77,105
k,111,-4.555,0,-5.195
k,112,4.555,0,-5.195
L,27,111
L,74,112
L,24,71
esize,,1!每根杆划分为一个单元
!H340X250X9X14工字钢
LATT,,,,,,,1
LMESH,1,5
LMESH,24,27
LMESH,46,50
LMESH,70,71
LMESH,87,91
LMESH,132,136
LMESH,110,113
LMESH,156,157
LMESH,192,201
LMESH,202,208,6
!H250X250X9X14工字钢
LATT,,,,,,,2
LMESH,6,23
LMESH,28,45
LMESH,52,54,2
LMESH,57,59,2
LMESH,62,64,2
LMESH,67,69,2
LMESH,92,109
LMESH,114,131
LMESH,138,140,2
LMESH,143,145,2
LMESH,148,150,2
LMESH,153,155,2
LMESH,181,191
LATT,,,,,,,3
LMESH,51,55,2
LMESH,56,60,2
LMESH,61,65,2
LMESH,66,68,2
LMESH,72,86
LMESH,137,141,2
LMESH,142,146,2
LMESH,147,151,2
LMESH,152,154,2
LMESH,158,172
LMESH,203,207
LATT,,,,,,,4
LMESH,173,178
LMESH,209,220
LATT,,,,,,,5
LMESH,179,180
LATT,,,,,,,6
LMESH,221,222
LATT,,,,,,,1
LMESH,223!最上面四个顶点处
TYPE,2
REAL,1
E,48
E,66
E,84
E,102
TYPE,2
REAL,2
E,39
E,38
E,75
E,74
NUMMRG,NODE, , , ,LOW ! 合并节点
NUMCMP,NODE ! 压缩节点编号
/solu
!静力分析
antype,static
acel,0,9.8,!定义重力加速度,系统自动算出井架本身重量
!绳载荷施加到顶端四个点
fk,24,fy,-4346
fk,47,fy,-4346
fk,71,fy,-4346
fk,94,fy,-4346
!风荷载(侧面来风,从-x轴往+x轴吹风)
*do,i,1,5,1!(侧前一段:17.989KN/9)
fk,i,fx,1999
*enddo
*do,i,25,28,1
fk,i,fx,1999
*enddo
*do,i,48,52,1!(侧后一段:4.499KN/9)
fk,i,fx,500
*enddo
*do,i,72,75,1
fk,i,fx,500
*enddo
*do,i,6,10,1!(侧前二段:27.134KN/10)
fk,i,fx,2713
*enddo
*do,i,29,33,1
fk,i,fx,2713
*enddo
*do,i,53,57,1!(侧后二段:10.854KN/10)
fk,i,fx,1085
*enddo
*do,i,76,80,1
fk,i,fx,1085
*enddo
*do,i,11,15,1!(侧前三段:28.822KN/10)
fk,i,fx,2882
*enddo
*do,i,34,38,1
fk,i,fx,2882
*enddo
*do,i,58,62,1!(侧后三段:10.944KN/10)
fk,i,fx,1094
*enddo
*do,i,81,85,1
fk,i,fx,1094
*enddo
*do,i,16,20,1!(侧前四段:29.767KN/10)
fk,i,fx,2977
*enddo
*do,i,39,43,1
fk,i,fx,2977
*enddo
*do,i,63,67,1!(侧后四段:10.417KN/10)
fk,i,fx,1042
*enddo
*do,i,86,90,1
fk,i,fx,1042
*enddo
*do,i,21,24,1!(侧前五段:21.424KN/8)
fk,i,fx,2678
*enddo
*do,i,44,47,1
fk,i,fx,2678
*enddo
*do,i,68,71,1!(侧后五段:6.641KN/8)
fk,i,fx,830
*enddo
*do,i,91,94,1
fk,i,fx,830
*enddo
!定义边界条件
dk,1,all
dk,48,all
dk,112,all
dk,111,all
Solve
到这一步,结果如下图所示:
然后进入后处理,查看最大等效应力云图:
显示最大等效应力为248MPa
然后通过建立单元表,列表显示出单元的轴向应力、弯曲应力、最大等效应力,最小等效应力,进而找出最大等效应力单元
其中 :最大等效应力=轴向应力+弯曲应力 最小等效应力=轴向应力-弯曲应力
ETABLE,SDIR,SMISC,31
ETABLE,SBZT,SMISC,34
SADD,YL_MAX,SDIR,SBZT
SADD,YL_MIN,SDIR,SBZT,,-1
/page,99999,132,99999,240!输出小技巧
PRETAB,SDIR,SBZT,YL_MAX,YL_MIN
问题:由列表可知最大等效应力分布在198单元,为130.2MPa,与云图上显示的最大等效应力248MPa不一样