在ANSYS里施加地震惯性力的方法Word版.doc

传播优秀Word版文档 ，希望对您有帮助，可双击去除！[分享]在ANSYS里施加地震惯性力的方法[b][size=16px]在ANSYS里施加地震惯性力的方法[/size][/b]在ANSYS里做地震分析时，需要对结构施加地震惯性荷载，地震惯性力是通过加速度的方式输入进结构的，然后与结构的质量一起形成动力计算时的惯性荷载，下面说一下在ANSYS里施加地震惯性力的方法首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：[font=宋体] -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[/font].......................然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：[font=宋体] 0.100000E-01 0.200000E-01 0.300000E-01 0.400000E-01 0.500000E-01 0.600000E-01[/font].......................编写如下的命令流文件，并命名为acce.inp[font=宋体]*dim,ACCEXYZ,TABLE,2000,3 !01行*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000 !02行[/font][font=宋体](3e16.6) !03行[/font][font=宋体]*vread,ACCEXYZ(1,0),time,txt !04行(e16.6) !05行[/font][font=宋体]ACCEXYZ(0,1)=1 !06行ACCEXYZ(0,2)=2 !07行，同上ACCEXYZ(0,3)=3 !08行，同上[/font]finish[font=宋体]/SOLUANTYPE,trans btime=0.01 !定义计算起始时间etime=15.00 !定义计算结束时间dtime=0.01 !定义计算时间步长*DO,itime,btime,etime,dtimetime,itime[/font][font=宋体][font=宋体]AUTOTS,0传播优秀Word版文档 ，希望对您有帮助，可双击去除！NSUBST,1, , ,1 KBC,1[/font]acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3) !施加三个方向的地震加速度SOLVE*ENDDO[/font]最后，在命令窗口里输入/input,acce,inp即可对结构进行地震动力分析。

[b]说明和讨论：[/b]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地震分析算例传播优秀Word版文档 ，希望对您有帮助，可双击去除！ !一个四层弹簧－质点模型的地震ANSYS分析!Example of seismic analysis in ANSYS!包括频率分析，谱分析和时程分析!With model & frequency analysis, spectrum analysis and transient analysis!作者：陆新征，清华大学土木工程系!Author: Lu Xinzheng Dept. Civil Engrg. of Tsinghua University [Money=20] FINI/CLEAR/UNITS,SI /PREP7 DAMPRATIO=0.02 !振型阻尼系数都为0.02ET,1,BEAM188KEYOPT,1,7,1KEYOPT,1,8,1ET,2,MASS21,,,4MP,EX,1,210E9 MP,NUXY,1,0.3 MP,DENS,1,7800 MPDATA,DAMP,1,,DAMPRATIOSECTYPE,3,BEAM,ISECDATA,0.2,0.2,0.5,0.01,0.01,0.008R,1,1.6E2,1.6E2R,2,1.2E2,1.2E2K,1,0,0,0 K,2,0,3,0 K,3,0,6,0 K,4,0,9,0 K,5,0,12,0 K,6,0,0,10000L,1,2L,2,3L,3,4L,4,5LSEL,,,,1,4传播优秀Word版文档 ，希望对您有帮助，可双击去除！LATT,1,0,1,,6,,3 ALLSLESIZE,ALL, , ,3 , , , , ,1 LMESH,ALL NUMMRG,ALLNUMCMP,ALLTYPE,2REAL,1E,2TYPE,2REAL,2E,8E,14E,20*SET,NT,1001*SET,DT,0.02 *DIM,AC,,NT *VREAD,AC(1),RECORD,TXT(F8.3)/SOLU D,1,ALL,ANTYPE,0SOLVEFINISH/SOLU !模态分析ANTYPE,2MODOPT,SUBSP,8 MXPAND,8, , ,1 SOLVEFINI!得到自振频率1*GET,FREQ1,MODE,1,FREQ/ESHAPE,1.0 /SOLUANTYPE,SPECTR !谱分析SPOPT,SPRSMDAMP,1,DAMPRATIO, , , , , , SED,1SVTYPE,2,1FREQ,0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667SV,,0.154350625,0.191590625,0.210210625,0.224175625,0.250716714,0.329514922,0.484352764,0.598723486,0.786897371FREQ,2.222,10,100000SV,,1.034212766,1.034212766,0.3528传播优秀Word版文档 ，希望对您有帮助，可双击去除！!按7度多遇地震，地震影响系数0.08 第一组III场地TG＝0.45SRSS,0.01,DISPSOLVEFINISH/POST1/INP,,MCOMNSORT,U,SUM *GET,DMAX_SPE,SORT, ,MAX !反应谱最大变形!开始地震计算FINISH/CONFIG,NRES,20000/SOLU ANTYPE,TRANS TRNOPT,FULL ALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)*DO,I,1,1001ACEL,AC(I),0,0 TIME,I*0.02 OUTRES,ALL,ALLSOLVE *ENDDO FINISH /POST26 FILENSOL,2,20,U,X, PLVAR,2, , , , , , , , , , STORE,MERGE XVAR,1 PLVAR,2,*GET,DMAX_TRANS,VARI,2,EXTREM,VMAX, , ! 时程分析最大位移数据文件record.txt-0.0014-0.0108-0.0101-0.0088传播优秀Word版文档 ，希望对您有帮助，可双击去除！-0.0095-0.012-0.0142-0.0128-0.011-0.0085-0.0085-0.0131-0.0176-0.0194-0.0162-0.0144-0.0108-0.0082-0.。