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

[分享]在ANSYS里施加地震惯性力的方法[b][size=16px]在ANSYS里施加地震惯性力的方法[/size][/b]在ANSYS里做地震分析时，需要对结构施加地震惯性荷载，地震惯性力是通过加速度的方式输入进结构的，然后与结构的质量一起形成动力计算时的惯性荷载，下面说一下在ANSYS里施加地震惯性力的方法。首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：[font=宋体]-0.227109E-02-0.209046E+000.467072E+01-0.413893E-02-0.168195E+000.261523E+01-0.574753E-02-0.157890E+000.809014E-01-0.731227E-02-0.152996E+000.119975E+01-0.876865E-02-0.138102E+000.130902E+01-0.101067E-01-0.131582E+000.143611E+00[/font].......................然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：[font=宋体]0.100000E-010.200000E-010.300000E-010.400000E-010.500000E-010.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,transbtime=0.01!定义计算起始时间etime=15.00!定义计算结束时间dtime=0.01!定义计算时间步长*DO,itime,btime,etime,dtimetime,itime[/font][font=宋体][font=宋体]AUTOTS,0NSUBST,1,,,1KBC,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.28761.23332.9938”这样的格式，此行的数据格式也就修改为“(3f10.4)”这样的格式。04行：从数据文件time.txt里读时间值到数组ACCEXYZ的第零列06行：将数组ACCEXYZ的第零行赋值，如果不对行插值的话也可以不赋值2、ANSYS中施加地震动荷载的讨论本贴中所说的在结构上施加惯性荷载的方法仅适用于对刚性基础上的结构的地震分析，对于考虑结构-地基相互作用的分析，由于需要考虑地基对结构的影响，其中最重要的因素就是地基辐射阻尼的影响，这在ANSYS里还不太好做。目前用得比较多的，也是最常用的方法就是用无质量地基模型，即在建立结构和地基的模型后，对地基周围的截断边界加固定约束，地基的密度设置为零，但地基弹性模量不为零，考虑地基的刚度影响，然后用上面的方法加整体惯性力。但是，这种方法由于没有能够考虑地基的辐射阻尼，往往导致计算结构偏于保守（计算结果往往较考虑地基辐射阻尼时要大）。ANSYS地震分析算例!一个四层弹簧－质点模型的地震ANSYS分析!ExampleofseismicanalysisinANSYS!包括频率分析，谱分析和时程分析!Withmodel&frequencyanalysis,spectrumanalysisandtransientanalysis!作者：陆新征，清华大学土木工程系!Author:LuXinzhengDept.CivilEngrg.ofTsinghuaUniversity[Money=20]FINI/CLEAR/UNITS,SI/PREP7DAMPRATIO=0.02!振型阻尼系数都为0.02ET,1,BEAM188KEYOPT,1,7,1KEYOPT,1,8,1ET,2,MASS21,,,4MP,EX,1,210E9MP,NUXY,1,0.3MP,DENS,1,7800MPDATA,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,0K,2,0,3,0K,3,0,6,0K,4,0,9,0K,5,0,12,0K,6,0,0,10000L,1,2L,2,3L,3,4L,4,5LSEL,,,,1,4LATT,1,0,1,,6,,3ALLSLESIZE,ALL,,,3,,,,,1LMESH,ALLNUMMRG,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)/SOLUD,1,ALL,ANTYPE,0SOLVEFINISH/SOLU!模态分析ANTYPE,2MODOPT,SUBSP,8MXPAND,8,,,1SOLVEFINI!得到自振频率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!按7度多遇地震，地震影响系数0.08第一组III场地TG＝0.45SRSS,0.01,DISPSOLVEFINISH/POST1/INP,,MCOMNSORT,U,SUM*GET,DMAX_SPE,SORT,,MAX!反应谱最大变形!开始地震计算FINISH/CONFIG,NRES,20000/SOLUANTYPE,TRANSTRNOPT,FULLALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)*DO,I,1,1001ACEL,AC(I),0,0TIME,I*0.02OUTRES,ALL,ALLSOLVE*ENDDOFINISH/POST26FILENSOL,2,20,U,X,PLVAR,2,,,,,,,,,,STORE,MERGEXVAR,1PLVAR,2,*GET,DMAX_TRANS,VARI,2,EXTREM,VMAX,,!时程分析最大位移数据文件record.txt-0.0014-0.0108-0.0101-0.0088-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.0042-0.0066-0.0131-0.019-0.0196-0.00660.0030.0141-0.0049-0.0128-0.0144-0.0203-0.026-0.0325-0.0306-0.0172-0.0197-0.0163-0.0164-0.00670.00250.0150.02360.02520.03360.04630.04920.04190.03590.02710.02350.03390.04120.0530.06390.07320.06520.05990.040.040.0063-0.0515-0.0787-0.0603-0.0484-0.025-0.00590.01340.03080.04990.0710.09950.12190.15290.14490.11550.09350.08920.09260.08390.09010.09930.12090.0328-0.1475-0.2066-0.1989-0.2034-0.1816-0.1725-0.1752-0.1753-0.1805-0.163-0.1347-0.1087-0.0782-0.0429-0.00170.0360.07850.11640.15980.1960.24120.27290.30360.320.34170.28210.2324-0.1198-0.2373-0.164-0.1865-0.1095-0.0753-0.01730.01130.05330.08950.11860.17570.0576-0.2631-0.1547-0.1729-0.1012-0.05790.0237-0.067-0.198-0.1641-0.1685-0.1481-0.1231-0.1001-0.0751-0.0523-0.0271-0.00440.0188-0.0095-0.0433-0.0838-0.0951-0.0716-0.0599-0.0334-0.01080.01850.0420.0673-0.0097-0.0372-0.0040.00110.03440.05650.08830.1130.13630.02190.02410.06830.06890.13180.13530.204-0.0931-0.1308-0.0692-0.05460.00720.0675-0.1067-0.1488-0.1071-0.1162-0.0762-0.0559-0.0215-0.0126-0.0674-0.0324-0.0337-0