钢筋混凝土管的动力响应特性分析实例

发布日期:[08-12-04 12:37:31] 浏览人次:[]

信号持续时间,采用正弦信号激震,也可采用RICKER信号,但要从数组文件中读取,可考虑设定一FLAG控制
XHCXSJ=NINT(1/ZYPL/DT)+1 !信号持续的时间步长
T0I=NINT(XHCXSJ/2)
F=ZYPL!经MATLAB实际输出后获得的信号实际频率

!######信号记录步长处理#########
*IF,DST/DT,LT,NINT(DST/DT),THEN !计算信号记录的时间步长间隔,实际记录的时间步长与设定值可能不同
RECORDSTEP=NINT(DST/DT)-1
*ELSE
RECORDSTEP=NINT(DST/DT)
*ENDIF
!#################时间步计算结束#########################
!##################单元剖分################
ET,1,PLANE183

*DO,I,1,JZFQ+BJCSM+1
LOCAL, 11, 1, 0, 0!定义局部坐标系,柱体
CSYS,11
LSEL,S,LOC,X,JHRB(I),
LESIZE,ALL,,R1ANGLE
ALLSEL,ALL
*ENDDO
CSYS,0

AMESH,ALL
AREFINE,1,9,,1!精细剖分
AREFINE,1,4,,1!精细剖分

!#################赋予材料属性#################
*DO,I,1,TCCS
MP,DENS,I,TC(I,1)
MP,NUXY,I,TC(I,2)
MP,EX,I,TC(I,3)
*ENDDO

LOCAL, 11, 1, 0, 0!定义局部坐标系,柱体
CSYS,11
*DO,I,1,JZFQ+BJCSM+1
NSEL,S,LOC,X,JHRA(I,1),JHRA(I,2)
NSEL,R,LOC,Y,JHRA(I,3),JHRA(I,4)
ESLN,S !边界节点单元
MAT,AMAT(I)
EMODIFY,ALL!确认材料参数修改
ALLSEL,ALL !重新选择ALL,准备下一步操作
*ENDDO

CSYS,0
/PNUM,MAT,1
/REPLOT
EPLOT

LOCAL, 11, 1, 0, 0!定义局部坐标系,柱体
CSYS,11

LSEL,S,LOC,X,R1
NSEL,ALL
NSLL,S,1
*GET, R1NUM, NODE, 0, COUNT
*GET, R1NMIN, NODE, 0,NUM, MIN!获取最小节点编号

*DIM,R1NODE,ARRAY,R1NUM,1!定义内圈节点数组
*DIM,R1NODEXY,ARRAY,R1NUM,2!定义内圈节点坐标数组,存储极坐标

R1NODE(1,1)=R1NMIN !赋内圈节点编号的最小值
*GET,RNX, NODE, R1NODE(1,1), LOC, X!获取内圈的节点的坐标,
*GET,RNY, NODE, R1NODE(1,1), LOC, Y!在此为极坐标
R1NODEXY(1,1)=RNX
R1NODEXY(1,2)=RNY

*DO,I,2,R1NUM
R1NODE(I,1)=NDNEXT(R1NODE(I-1))
!对内圈节点赋值,NDNEXT的输出为选定的节点中编号大于且最接近该节点的节点编号
*GET,RNX, NODE, R1NODE(I,1), LOC, X!获取内圈的节点的坐标,
*GET,RNY, NODE, R1NODE(I,1), LOC, Y!在此为极坐标
R1NODEXY(I,1)=RNX
R1NODEXY(I,2)=RNY
*ENDDO

!###############按极坐标中角度,由小到大排列节点编号和对应坐标
*DO,I,1,R1NUM-1
*DO,J,I+1,R1NUM
*IF,R1NODEXY(I,2),GT,R1NODEXY(J,2),THEN
RY=R1NODEXY(I,2)
R1NODEXY(I,2)=R1NODEXY(J,2)
R1NODEXY(J,2)=RY
RX=R1NODEXY(I,1)
R1NODEXY(I,1)=R1NODEXY(J,1)
R1NODEXY(J,1)=RX
RNODE=R1NODE(I)
R1NODE(I)=R1NODE(J)
R1NODE(J)=RNODE
*ENDIF
*ENDDO
*ENDDO
!####################排序OVER#####################
*CFOPEN,R1NODEXY,DAT!输出井壁内表面节点编号和坐标
*DO,K,1,R1NUM
NUM_NODE=R1NODE(K)
RX=R1NODEXY(K,1)
RY=R1NODEXY(K,2)
*VWRITE,NUM_NODE,RX,RY
(F10.0,F13.2,F13.2)
*ENDDO
*CFCLOSE

ALLSEL,ALL
FINISH

!#################求解过程##########################
!##与竖直方向采用直角坐标的不同:
!##1、震源位置变换方式
!##2、记录点位置
!##3、边界条件设定,不作专门的吸收边界,而以空间扩大来处理

/SOLU
*CFOPEN,USUM,DAT
!采用REIKER子波,在初始时段施加
NSUBST,1, , ,1 !1个子步
ANTYPE,TRANS !时程分析
LUMPM,0
OUTPR,BASIC,SCBCJG, !每隔一定时间步长输出基本值

NSEL,S,LOC,
|<< << < 1 2 3 > >> >>|
www.mapeng.net 马棚网
www.mapeng.net
文章作者:未知 | 文章来源:网络 | 责任编辑:admin | 发送至邮箱: | 加入收藏:
本文关键字:钢筋  混凝土    动力响应  特性  分析  实例
相关资讯
热点资讯
推荐资讯

关于我们 | 站点导航 | 使用帮助 | 友情链接 | 广告服务 | 免责声明 | 新手上路
设为首页 | 加入收藏 | 在线留言 | 马棚网QQ群:{92562572}{102901272}{333259257} | 交流QQ: 客户服务 客户服务 客户服务