【原创】框架结构的时程分析!
:)偶作的框架结构的时程分析,望斑竹可以给俺加一分,鼓励一下,到现偶还是0分:(,/PREP7
/color,pbak,1
/rep
ET,1,BEAM188
mp,ex,1,25.5e9
mp,ex,2,28e9
mp,gxy,1,9.8e9
mp,gxy,2,10.77e9
mp,prxy,1,0.2
mp,dens,1,2500
!difine sections
SECTYPE,1,BEAM, RECT
SECOFFSET,CENT
SECDATA,.45,.45,0,0,0,0,0,0,0,0
SECTYPE,2,BEAM, RECT
SECOFFSET,CENT
SECDATA,0.25,.5,0,0,0,0,0,0,0,0
SECTYPE,3, BEAM,RECT
SECOFFSET,CENT
SECDATA,0.25,0.7,0,0,0,0,0,0,0,0
SECTYPE,4, BEAM, RECT, , 0
SECOFFSET, CENT
SECDATA,0.45,0.55,0,0,0,0,0,0,0,0
k,1,
k,2,5,
k,3,12.2
!define the direction points of element beam188
k,99993,1,
k,99994,1,,6
k,99995,1,,12
k,99996,1,,18
k,99997,1,,24
k,99998,1,,30
k,100000,,1
k,100001,,1,6
k,100002,,1,12
k,100003,,1,18
k,100004,,1,24
k,100005,,1,30
k,100006,,-1
k,100007,5,
k,100008,12.2,
!creat keypoints
kgen,2,1,3,1,,4,,3
kgen,5,4,6,1,,3.6,,3
*do,j,1,13,3
*do,i,j,j+2,1
l,i,i+3
*enddo
*enddo
*do,i,4,16,3
l,i,i+1
l,i+1,i+2
*enddo
/VIEW,1,1,1,1
/ANG,1
/REP,FAST
lsel,all
lgen,6,all,,,,,6
!define sections of columns
lsel,s,line,,1,3,1
latt,2,1,1,,99993,,4
lsel,s,line,,4,15,1
latt,2,1,1,,99993,,1
lsel,s,line,,26,28,1
latt,2,1,1,,99994,,4
lsel,s,line,,29,40,1
latt,2,1,1,,99994,,1
lsel,s,line,,51,53,1
latt,2,1,1,,99995,,4
lsel,s,line,,54,65,1
latt,2,1,1,,99995,,1
lsel,s,line,,76,78,1
latt,2,1,1,,99996,,4
lsel,s,line,,79,90,1
latt,2,1,1,,99996,,1
lsel,s,line,,101,103,1
latt,2,1,1,,99997,,4
lsel,s,line,,104,115,1
latt,2,1,1,,99997,,1
lsel,s,line,,126,128,1
latt,2,1,1,,99998,,4
lsel,s,line,,129,140,1
latt,2,1,1,,99998,,1
!define the section of beams
lsel,s,line,,17,25,2
lsel,a,line,,16
latt,1,1,1,,100000,,3
lsel,s,line,,18,24,2
latt,1,1,1,,100000,,2
lsel,s,line,,42,50,2
lsel,a,line,,41
latt,1,1,1,,100001,,3
lsel,s,line,,43,49,2
latt,1,1,1,,100001,,2
lsel,s,line,,67,75,2
lsel,a,line,,66
latt,1,1,1,,100002,,3
lsel,s,line,,68,74,2
latt,1,1,1,,100002,,2
lsel,s,line,,92,100,2
lsel,a,line,,91
latt,1,1,1,,100003,,3
lsel,s,line,,93,99,2
latt,1,1,1,,100003,,2
lsel,s,line,,117,125,2
lsel,a,line,,116
latt,1,1,1,,100004,,3
lsel,s,line,,118,124,2
latt,1,1,1,,100004,,2
lsel,s,line,,142,150,2
lsel,a,line,,141
latt,1,1,1,,100005,,3
lsel,s,line,,143,149,2
latt,1,1,1,,100005,,2
!creat beams of the long
*do,k,0,2,1
*do,i,3,6,1
*do,j,0,4,1
l,3*i-k+18*j,3*i-k+18*(j+1)
*enddo
*enddo
*enddo
*do,j,0,4,1
*SET,i,6+18*j
l,i,i+18
*enddo
*do,j,1,4,1
*SET,i,2+18*j
l,i,i+18
*enddo
l,20,4
*do,j,1,4,1
*SET,i,4+18*j
l,i,i+18
*enddo
l,22,5
!define sections
lsel,s,line,,191,210,1
lsel,a,line,,216,220,1
latt,1,1,1,,100006,,3
lsel,s,line,,171,190,1
lsel,a,line,,221,225,1
latt,1,1,1,,100007,,3
lsel,s,line,,151,170,1
lsel,a,line,,211,215,1
latt,1,1,1,,100008,,3
!creat elememts
lsel,all
lesize,all,,,20
lmesh,all
/eshape,1
eplot
!delete the direction points
kdele,99993,99998,1
kdele,100000,100008,1
/solu
ksel,s,loc,y,0,,
dk,all,all,0
/eshape,0
/rep
save
AnType,Modal
ModOpt,Lanb,15
Solve
Finish
/post1
set,list
set,1,15
pldisp0
ancntr10,0.5
set,next
pldisp,2
fini
/solu
NT=100!时程曲线有NT个点
DT=0.02 !时间间隔
*dim,ac,,NT
/input,tianjin,txt !天津波东西向共2秒
NSUBST,1, , ,1 !1个子步
OUTRES,ALL,1 !输出每个子步的结果
ANTYPE,TRANS !时程分析
*do,i,1,NT
ACEL,ac(i),0,0
TIME,i*DT
solve
*enddo
save
/post26
NSOL,2,484,U,x, UX_2
STORE,MERGE
NSOL,3,484,U,Y, UY_3
STORE,MERGE
NSOL,4,484,U,Z, UZ_4
STORE,MERGE
NSOL,5,8191,U,x, UX_5
STORE,MERGE
NSOL,6,8191,U,Y, UY_6
STORE,MERGE
NSOL,7,8191,U,Z, UZ_7
STORE,MERGE
prvar,2,3,4,5,6,7 !列表显示位移
XVAR,1
PLVAR,2,5 !绘制node484和8191的x方向位移曲线
XVAR,1
PLVAR,3,6
XVAR,1
PLVAR,4,7
! /EXIT,ALL
回复: 【原创】框架结构的时程分析!
附一图:回复: 【原创】框架结构的时程分析!
注意:/input,tianjin,txt
需要在working directory下建立一个文本文件:tianjin.txt, 我修改后的天津波已附后)。
回复: 【原创】框架结构的时程分析!
useful. these days I am working this problem.thanks a lot!
回复: 【原创】框架结构的时程分析!
太感谢了回复: 【原创】框架结构的时程分析!
将“NSOL,7,8191,U,Z, UZ_7 ”作适当修改即可的,偶得回去试试,
回复: 【原创】框架结构的时程分析!
做时程分析的时候,要用到恢复力模型,请问怎么样才能做出来。你这个好像没有用到这个。如果会,可不可以贴出来Re:回复: 【原创】框架结构的时程分析!
不好意思,偶现在不用ANSYS了,恢复力模型在ANSYS里面偶没试过,目前改用ABAQUS了做了,研究中ING 好东西哦 楼主辛苦 不知道我什么时候能做成谢谢!!!
蛮实用的哦 好帖子,顶 好东西,学习学习 Good files ,thanks you! 好东西,谢谢楼主 不错 2# jacobi 先学习一下 先学习一下,谢谢啊 谢了,非常感谢。
页:
[1]
2