在ANSYS里施加地震惯性力的方法

在ANSYS里做地震分析时,需要对结构施加地震惯性荷载,地震惯性力是通过加速度的方式输入进结构的,然后与结构的质量一起形成动力计算时的惯性荷载,下面说一下在ANSYS里施加地震惯性力的方法。

首先,将三个方向的地震加速度放到一个文本文件里,如accexyz.txt,在这个数据文件里共放三列数据,每列为一个方向的地震加速度值,这里仅给出数据文件中前几行的数据: -0.227109E-02 -0.209046E+00 0.467072E+01

-0.413893E-02 -0.168195E+00 0.261523E+01

-0.574753E-02 -0.157890E+00 0.809014E-01

-0.731227E-02 -0.152996E+00 0.119975E+01

-0.876865E-02 -0.138102E+00 0.130902E+01

-0.101067E-01 -0.131582E+00 0.143611E+00

.......................

然后,再建一个文本文件用来存放三个方向的地震加速度时间点,如time.txt,在这个数据文件里仅一列数据,对应于加速度数据文件里每一行的时间点,这里给出数据文件中前几行数据:

0.100000E-01

0.200000E-01

0.300000E-01

0.400000E-01

0.500000E-01

0.600000E-01

.......................

编写如下的命令流文件,并命名为acce.inp

*dim,ACCEXYZ,TABLE,2000,3 !01行

*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行

(3e16.6) !03行

*vread,ACCEXYZ(1,0),time,txt !04行

(e16.6) !05行

ACCEXYZ(0,1)=1 !06行

ACCEXYZ(0,2)=2 !07行,同上

ACCEXYZ(0,3)=3 !08行,同上

finish

/SOLU

ANTYPE,trans

btime=0.01 !定义计算起始时间

etime=15.00 !定义计算结束时间

dtime=0.01 !定义计算时间步长

*DO,itime,btime,etime,dtime

time,itime

AUTOTS,0

NSUBST,1, , ,1

KBC,1

acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度

SOLVE

*ENDDO

最后,在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。

说明和讨论:

1、命令流中各行说明:

01行:定义2000行,3列的数组,(行数根据数据文件里加速度点数来定)ACCEXYZ用来存放三个方向的加速度值

02行:从数据文件accexyz.txt里读加速度值到数组ACCEXYZ,2000为行数,可根据情况修改。其中的JIK,3,2000非常重要,它决定着将从加速度数据文件中数据输入到ACCEXYZ数组时的存放格式,这里用到按JIK方式变化,即读进来的数据依次放入ACCEXYZ(I,J)(J从1到3,I从1到2000,K默认从1到1)。根据这行命令下面(3e16.6)的格式,每次从数据文件里读一行三个数据,分别放入ACCEXYZ(1,1)、ACCEXYZ(1,2)、ACCEXYZ(1,3),

接着再读下一行的三个数据,分别放入ACCEXYZ(2,1)、ACCEXYZ(2,2)、ACCEXYZ(2,3),依次类推。当然依据数据文件的格式,也可以选用IJK, IKJ, JIK, JKI, KIJ, KJI等格式,其中IJK为默认。

03行:读数据的格式,每行三个数值。由于数据文件中的数据是用科学记数年法表示的,因此,这里也用相应FORTRAN的科学记数法的格式。如果数据文件里的数值是如

“0.2876 1.2333 2.9938”这样的格式,此行的数据格式也就修改为“(3f10.4)”这样的格式。

04行:从数据文件time.txt里读时间值到数组ACCEXYZ的第零列

06行:将数组ACCEXYZ的第零行赋值,如果不对行插值的话也可以不赋值

2、ANSYS中施加地震动荷载的讨论

本贴中所说的在结构上施加惯性荷载的方法仅适用于对刚性基础上的结构的地震分析,对于考虑结构-地基相互作用的分析,由于需要考虑地基对结构的影响,其中最重要的因素就是地基辐射阻尼的影响,这在ANSYS里还不太好做。目前用得比较多的,也是最常用的方法就是用无质量地基模型,即在建立结构和地基的模型后,对地基周围的截断边界加固定约束,地基的密度设置为零,但地基弹性模量不为零,考虑地基的刚度影响,然后用上面的方法加整体惯性力。但是,这种方法由于没有能够考虑地基的辐射阻尼,往往导致计算结构偏于保守(计算结果往往较考虑地基辐射阻尼时要大)。

在ANSYS里做地震分析时,需要对结构施加地震惯性荷载,地震惯性力是通过加速度的方式输入进结构的,然后与结构的质量一起形成动力计算时的惯性荷载,下面说一下在ANSYS里施加地震惯性力的方法。

首先,将三个方向的地震加速度放到一个文本文件里,如accexyz.txt,在这个数据文件里共放三列数据,每列为一个方向的地震加速度值,这里仅给出数据文件中前几行的数据: -0.227109E-02 -0.209046E+00 0.467072E+01

-0.413893E-02 -0.168195E+00 0.261523E+01

-0.574753E-02 -0.157890E+00 0.809014E-01

-0.731227E-02 -0.152996E+00 0.119975E+01

-0.876865E-02 -0.138102E+00 0.130902E+01

-0.101067E-01 -0.131582E+00 0.143611E+00

.......................

然后,再建一个文本文件用来存放三个方向的地震加速度时间点,如time.txt,在这个数据文件里仅一列数据,对应于加速度数据文件里每一行的时间点,这里给出数据文件中前几行数据:

0.100000E-01

0.200000E-01

0.300000E-01

0.400000E-01

0.500000E-01

0.600000E-01

.......................

编写如下的命令流文件,并命名为acce.inp

*dim,ACCEXYZ,TABLE,2000,3 !01行

*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行

(3e16.6) !03行

*vread,ACCEXYZ(1,0),time,txt !04行

(e16.6) !05行

ACCEXYZ(0,1)=1 !06行

ACCEXYZ(0,2)=2 !07行,同上

ACCEXYZ(0,3)=3 !08行,同上

finish

/SOLU

ANTYPE,trans

btime=0.01 !定义计算起始时间

etime=15.00 !定义计算结束时间

dtime=0.01 !定义计算时间步长

*DO,itime,btime,etime,dtime

time,itime

AUTOTS,0

NSUBST,1, , ,1

KBC,1

acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度

SOLVE

*ENDDO

最后,在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。

说明和讨论:

1、命令流中各行说明:

01行:定义2000行,3列的数组,(行数根据数据文件里加速度点数来定)ACCEXYZ用来存放三个方向的加速度值

02行:从数据文件accexyz.txt里读加速度值到数组ACCEXYZ,2000为行数,可根据情况修改。其中的JIK,3,2000非常重要,它决定着将从加速度数据文件中数据输入到ACCEXYZ数组时的存放格式,这里用到按JIK方式变化,即读进来的数据依次放入ACCEXYZ(I,J)(J从1到3,I从1到2000,K默认从1到1)。根据这行命令下面(3e16.6)的格式,每次从数据文件里读一行三个数据,分别放入ACCEXYZ(1,1)、ACCEXYZ(1,2)、ACCEXYZ(1,3),

接着再读下一行的三个数据,分别放入ACCEXYZ(2,1)、ACCEXYZ(2,2)、ACCEXYZ(2,3),依次类推。当然依据数据文件的格式,也可以选用IJK, IKJ, JIK, JKI, KIJ, KJI等格式,其中IJK为默认。

03行:读数据的格式,每行三个数值。由于数据文件中的数据是用科学记数年法表示的,因此,这里也用相应FORTRAN的科学记数法的格式。如果数据文件里的数值是如

“0.2876 1.2333 2.9938”这样的格式,此行的数据格式也就修改为“(3f10.4)”这样的格式。

04行:从数据文件time.txt里读时间值到数组ACCEXYZ的第零列

06行:将数组ACCEXYZ的第零行赋值,如果不对行插值的话也可以不赋值

2、ANSYS中施加地震动荷载的讨论

本贴中所说的在结构上施加惯性荷载的方法仅适用于对刚性基础上的结构的地震分析,对于考虑结构-地基相互作用的分析,由于需要考虑地基对结构的影响,其中最重要的因素就是地基辐射阻尼的影响,这在ANSYS里还不太好做。目前用得比较多的,也是最常用的方法就是用无质量地基模型,即在建立结构和地基的模型后,对地基周围的截断边界加固定约束,地基的密度设置为零,但地基弹性模量不为零,考虑地基的刚度影响,然后用上面的方法加整体惯性力。但是,这种方法由于没有能够考虑地基的辐射阻尼,往往导致计算结构偏于保守(计算结果往往较考虑地基辐射阻尼时要大)。


相关内容

  • 地震波输入
  • 目前, 地震波输入以加速度时程输入为主. 在不考虑土-结构共同作用问题中, 该加速度时程是固定基底的加速度-时间曲线; 考虑土-结构共同作用问题中, 这一加速度时程是基岩运动的加速度-时间关系. 不同的地震波输入位置都可以称为地震波输入面. 从实质上讲, 在地震动力响应问题中, 地震波的时程曲线是地 ...

  • 基于ANSYS对工字钢梁结构的静力学分析
  • 基于ANSYS对工字钢梁结构的静力学分析 基于ANSYS对工字钢梁结构的静力学分析 李建莉 (兰州资源环境职业技术学院,甘肃 兰州 730021) 摘要:基于ANSYS对工字钢梁结构的静力学特性进行理论计算和研究.建立了工字梁的有限元模型,得出了在给定参数下工字梁的一些静力学结果,并对计算结果进行静 ...

  • 基于ANSYS二次开发的埋地管道应力分析软件及实现
  • the Winkler elastic footing beam theory. The model is used in re-development of ANSYS software by VC to highlight calculation equation for soil spring ...

  • 毕业论文--联轴器的有限元分析
  • 摘要 本文首先扼要得为大家介绍了联轴器的性能.功用.分类和有限元方法.接下来分别简要介绍凸缘联轴器和十字轴式万向联轴器的优缺点,接着完成开始有限元分析的前期准备,准备好两种联轴器的二维图形的尺寸,根据准备好的尺寸图开始在Proe中建立两种联轴器的三维建模.然后准备用Ansys开始有限元分析,再把三维 ...

  • ANSYS悬臂梁的自由端受力的有限元计算[1]
  • 悬臂梁自由端受力的有限元计算 任柳杰 [1**********] 一. 计算目的 1. 掌握ANSYS软件的基本几何形体构造.网格划分.边界条件施加等方法. 2. 熟悉有限元建模.求解及结果分析步骤和方法. 3. 利用ANSYS软件对梁结构进行有限元计算. 4. 梁的变形.挠曲线等情况的分析. 5. ...

  • ansys轴对称结构的静力分析
  • 第八章 轴对称结构的静力分析 在工程实践所应用的结构中,有许多结构是可以由一个截面绕某固定轴旋转而生成的,如果这种结构所受的外载荷和边界条件也沿此轴对称,则称此结构为轴对称结构.在有限元理论中对于此类结构有专门的简化方法,在ANSYS中也可以通过结构的轴对称性简化模型,减少模型规模.缩短计算时间,提 ...

  • ansys桁架结构分析实例
  • 平面桁架的静力分析 摘要:近些年来,ANSYS 工程软件在工程领域内运用的很多,在分析线性有限元模型上比其他软件更具有优势.而在ANSYS 软件中最常用的是线性静力分析,尽管很多的材料不一样,但结果确基本一致.本文主要是要对平面桁架进行静力分析. 关键字:线性:桁架:有限元:结构 The plane ...

  • 有限元题库
  • 有限元考试复习 1. 性与非线性问题 平面应力问题 (1) 均匀薄板(2)载荷平行于板面且沿厚度方向均匀分布 在六个应力分量中,只需要研究剩下的平行于XOY平面的三个应力分量,即σx.σy.τxy=τyx (σz=0,τzx=τxz=0,τzy=τyz=0). 一般σz=0,εz并不一定等于零,但可 ...

  • 半潜式钻井平台甲板结构极限强度分析
  • 江苏科技大学 科生毕业设计(论 学 院 船舶与海洋工程 专 业 船舶与海洋工程 学生姓名 王杰 班级学号 1040101128 指导老师 施兴华 二零一四年六月 文) 本 江苏科技大学本科生毕业论文 半潜式钻井平台甲板结构极限强度分析 Ultimate strength analysis of se ...