ANSYS渡槽槽身动水压力的施加（2）——U型渡槽

U型渡槽动水压力荷载施加命令及说明

• 程序中需要用到ANSYS重启动，因为需提取前一步加速度结果以施加部分动水压力；
• 默认Y方向为重力方向，X方向为横槽向，Z方向为纵槽向；
• 需准备地震波文件；
• 需先将槽底面和槽壁面编为组元：bottom_areas,left_areas,right_areas，U型渡槽下部曲线在水平向既是槽壁，在竖向又是槽底；

pi=3.1415926

H_water=3.6   !槽内水深

R_water=2.3    !半槽宽

rou_water=1000   !水体密度

M_water=rou_water*(2*R_water*(H_water-R_water)+0.5*pi*R_water**2)

 !控制点------截面槽底中心处节点
node_control=8450
coord_y_control=ny(node_control)
coord_x_control=nx(node_control)
coord_z_control=nz(node_control)
coord_y_water=coord_y_control+H_water

/prep7
!定义mass21单元
et,99,mass21
!定义弹簧单元
et,98,combin14
!初始实常数
nr=100

!读入地震波
*SET,NT,200      !地震加速度步数
*SET,DT,0.01      !地震加速度步长
*DIM,accel_x,,NT
*DIM,accel_y,,NT
*DIM,accel_z,,NT
*VREAD,accel_x,'accel_x','txt'
(F12.9)
*VREAD,accel_y,'accel_y','txt'
(F12.9)
*VREAD,accel_z,'accel_z','txt'
(F12.9)

!***********************************附加质量（开始）**********************************

 

!将附加质量的渡槽底面的所有节点编为一个数组

allsel

CMSEL,S,bottom_areas,AREA

NSLA,S,1

*get,n_bottom,node,,count

*dim,node_bottom,array,n_bottom

*get,nmin,node,,num,min

node_bottom(1)=nmin

*do,i,2,n_bottom

  *get,nnum,node,nmin,nxth

  nmin=nnum

  node_bottom(i)=nmin

*enddo

allsel

 

!将附加质量的渡槽左壁面的所有节点编为一个数组

CMSEL,S,left_areas,AREA

NSLA,S,1

nsel,U,loc,x,coord_x_control-0.001,coord_x_control+0.001

*get,n_left,node,,count

*dim,node_left,array,n_left

*get,nmin,node,,num,min

node_left(1)=nmin

*do,i,2,n_left

  *get,nnum,node,nmin,nxth

  nmin=nnum

  node_left(i)=nmin

*enddo

allsel

<