地震分析算例Word下载.docx
- 文档编号:19929147
- 上传时间:2023-01-12
- 格式:DOCX
- 页数:44
- 大小:95.09KB
地震分析算例Word下载.docx
《地震分析算例Word下载.docx》由会员分享,可在线阅读,更多相关《地震分析算例Word下载.docx(44页珍藏版)》请在冰豆网上搜索。
K,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,10000
需要说明的是,关键点6为后面建立梁单元所需的截面方向控制点,在上一个例子中已经做过介绍。
(13)完成关键点输入后下面建立直线模型。
在ANSYS主菜单中选择Preprocessor->
Lines->
StraightLine,依次连接关键点1~5。
(14)下面给建立完的几何模型赋予材料属性,在ANSYS菜单中选择Meshing->
MeshAttributes->
PickedLines,选中所有的直线,进入LineAttributes窗口,选择相关选项材料属性,实参数,单元类型和截面类型都为1,并选择勾上参考点方向选项,然后点击OK后输入关键点6作为截面方向控制点。
(15)在ANSYS主菜单中选择Preprocessor->
Meshing->
SizeCntrls->
ManualSize->
Lines->
PickedLines,在ElementSizesonPickedLines窗口中设定NDIVNo.ofelementdivisions为3,即将每条直线分为3段
(16)在ANSYS主菜单中选择Preprocessor->
Mesh->
Lines,选择所有的直线,完成直线的网格划分。
(17)为了便于后面操作,将网格划分后的单元和节点编号进行适当的清理。
NumberingCtrls->
MergeItems,在MergeCoincidentorEquivalentlyDefinedItems窗口中选择All,清理所有重复的元素。
同样选择NumberingCtrls->
CompressNumber菜单,在CompressNumber中选择All,对节点和单元进行重新编号。
(18)下面建立集中质量单元,采取直接输入单元的方法建立。
Elements->
ElemAttributes,在ElementAttributes中设定单元类型编号为2MASS21,材料编号任意,实参数编号为1。
(19)在ANSYS主菜单中选择Preprocessor->
AutoNumbered->
ThruNodes,选择节点2,建立第一个集中质量。
(20)再次进入第18步ElementAttributes窗口,设定实参数(Realconstantsetnumber)为2。
(21)重复19步,选择节点8,14,20,建立其他的三个集中质量单元。
(22)到此完成所有建模工作,下面开始进行结构分析
(23)进入ANSYS主菜单中Solution功能模块,选择Solution->
DefineLoads->
Apply->
Structural->
Displacement->
OnNodes,选择节点1,设定约束所有的自由度。
(24)首先做一次静力分析,选择ANSYS主菜单Solution->
AnalysisType->
NewAnalysis,设定分析类型为Static
(25)选择ANSYS主菜单Solution->
Solve->
CurrentLS选项,进行一次静力分析
(26)分析完后,下面进行模态分析,在ANSYS主菜单中选择Solution->
NewAnalysis,选择分析类型为Model
(27)在ANSYS主菜单中选择Solution->
AnalysisOptions,输入模态分析方法为子空间法(Subspace),求解8阶模态,同时需要作模态扩展,扩展的模态为8阶,并计算单元应力和应变即改为yes,输入窗口如图。
(28)再次选择ANSYS主菜单Solution->
CurrentLS,计算当前问题
(29)这时,如果需要看结果,可以进入后处理模块,即ANSYS主菜单GeneralPostproc,可以看到计算的各阶频率和振型。
!
存储各模态频率
*DIM,FRE,,NMODE
*DO,I,1,NMODE
*GET,FRE(I),MODE,I,FREQ!
OBTAINMODEFREQENCYFORMODEI
*ENDDO
FINISH
ALLSEL
(30)完成结构自振分析后,下面就可以进行反应谱分析
(31)进入ANSYS主菜单Solution->
NewAnalysis,选择分析类型为Spectrum。
(32)在ANSYS主菜单Solution->
AnalysisOptions中,选择谱分析的类型为单点输入(Sing-ptresp)
(33)在ANSYS主菜单中选择Preprocessor->
Loads->
LoadStepOpts->
Time/Frequenc->
Damping,输入所有的阻尼为DAMPRATIO,如图所示
(34)下面需要定义地震的反应谱。
我国规范给定的是基于加速度的反应谱。
在ANSYS主菜单中选择Solution->
LoadStepOpts->
Spectrum->
SinglePoint->
Settings,设定反应谱类型为地震加速度,放大系数为1,输入方向为X方向(1,0,0)。
(35)接下来开始输入地震反应谱。
注意:
在进行反应谱前,若按7度多遇地震,即地震影响系数0.08第一组III场地TG=0.45,可以通过如下方式从国标给定的地震谱频率中的结果中算出相应的谱值,不过其频率需要提前定义数组。
过程如下:
利用前面模态分析得到的频率。
(因为这些频率的取值附近容易引起共振,破坏最大)
计算反应谱数据(依据规范GB50011-2001第5.1.5条)
地震影响系数
GRAV=9.81!
重力加速度G
TG=0.45!
特性周期
AMAX=0.08!
水平地震影响系数最大值
C=0.05!
阻尼比
*DIM,A,,NMODE
*DIM,T,,NMODE
T(I)=1.0/FRE(I)
R=0.9+(0.05-C)/(0.5+5.0*C)
P1=0.02+(0.05-C)/8
P2=1+(0.05-C)/(0.06+1.7*C)
*IF,T(I),GE,0.0,AND,T(I),LT,0.1,THEN
A(I)=(0.45+(10.0*P2-4.5)*T(I))*AMAX*GRAV
*ELSEIF,T(I),GE,0.1,AND,T(I),LE,TG
A(I)=P2*AMAX*GRAV
*ELSEIF,T(I),GT,TG,AND,T(I),LE,5*TG
A(I)=(TG/T(I))**R*P2*AMAX*GRAV
*ELSE
A(I)=(P2*0.2**R-P1*(T(I)-5*TG))*AMAX*GRAV
*ENDIF
这里输入的反应谱按7度多遇地震,取地震影响系数为0.08,第一组,III类场地,卓越周期Tg=0.45s。
选择ANSYS主菜单Solution->
FreqTable,输入频率反应谱Freq1~Freq12为0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667,2.222,10,100000,如图
(36)选择ANSYS主菜单Solution->
SpectrValues,输入对应的反应谱数值依次如下:
(37)最后选择ANSYS主菜单Solution->
ModeCombine,设定振型组合方式为SRSS法,如图所示
(38)选择ANSYS主菜单Solution->
CurrentLS,计算反应谱结果
(39)进入ANSYS主菜单后处理模块GeneralPostproc,在ANSYS窗口顶部菜单选择File->
Readinputfrom,选择文件后缀名为*.mcom的文件
(40)进入ANSYS主菜单GeneralPostproc->
PlotResults->
DeformedShapes,选择绘制变形后形状和结构形状,得到地震反应谱分析的结构变形如图
(41)最后我们来进行地震时程分析,进行地震时程分析以前,首先要有一个地震时程记录,本例子给定的地震时程记录总长20秒,记录点间隔0.02s,共有1001个记录点。
该地震记录存放在RECORD.TXT文件中。
(42)首先建立两个变量,在ANSYS窗口顶部菜单选择Parameters->
Scalarparameters,在窗口中输入NT=1001,即总共1001个记录点,DT=0.02,即记录点间隔0.02s
(43)然后还需要定义一个数组来存放地震记录,选择ANSYS窗口顶部菜单Parameters->
Arrayparameters->
Define/Edit,在弹出的ArrayParameters窗口中点击Add按钮,在AddNewArrayParameter窗口中输入数组的名称为AC,Parametertype类型选择默认的Array,再I,J,K,的行数设为NT个,其余默认。
如图所示。
(44)下面从数据文件中读入地震时程记录。
选择ANSYS窗口顶部菜单Parameters->
Readfromfile,选择类型,Array,点击OK在弹出窗口中依次输入将数据读入的数组名称AC,可以用Browse选择目标数据文件RECORD.TXT,Label中选择IJK项,其中ROWS行的数改为NT,其余保持默认,并在最后要给出数据FormatSurroundedBy()读入的格式(F8.3)既可。
并且可以参考Fortran相应的文件输入输出要求。
(45)完成数据读入后就可以进行地震反应分析。
由于本次地震分析一共要进行1001步,超过了ANSYS默认的最大1000步限制,因此首先需要进行调整。
在ANSYS的命令输入窗口内,输入以下命令:
Finish
/CONFIG,NRES,20000
更改后可以进入,List-〉Status-〉configuration,中查看。
(46)进入ANSYS主菜单Solution->
NewAnalysis,指定分析类型为瞬态分析Transient,Solutionmethod可以选择Full。
(47)进入ANSYS主菜单Solution->
Sol'
nControls,在SolutionControls窗口里面选择Transient页面,设定瑞雷阻尼的数值。
由于阻尼的机理十分复杂,因此不同问题可能各不相同,这里仅介绍一种阻尼的取法。
输入质量阻尼系数为2*DAMPRATIO*FREQ1*2*3.1415926,刚度阻尼系数为2*DAMPRATIO/(FREQ1*2*3.1415926),FREQ1为结构的第一阶自振频率。
(48)下面将用一组循环控制语句来定义作用在结构上的加速度并进行计算。
读者可以将下面这部分语句用任意编辑器(比如windows的记事本)输入好以后复制ANSYS的输入栏中。
*DO,I,1,1001!
对变量I循环1001
ACEL,AC(I),0,0!
对结构施加X方向加速度
TIME,I*0.02!
计算时间步长为0.02s
OUTRES,ALL,ALL!
输出所有结果
SOLVE!
求解
*ENDDO!
循环结束
(49)下面介绍适用ANSYS的时程后处理器TimeHistPostpro来处理计算结果,在ANSYS主菜单上选择TimeHistPostpro,这时会弹出时程变量窗口TimeHistoryVariables。
点击工具栏上第一个绿色加号按钮,添加时程变量。
在AddTime-HistoryVariable窗口中选择要添加的变量为节点的X方向位移。
(50)选择第20号节点,回到时程变量窗口TimeHistoryVariables,点击工具栏上第三个按钮,即绘制出节点的位移时程曲线如图
一个四层弹簧-质点模型的地震ANSYS分析
ExampleofseismicanalysisinANSYS
包括频率分析,谱分析和时程分析
Withmodel&
frequencyanalysis,spectrumanalysisandtransientanalysis
FINI
/CLEAR
/UNITS,SI
/PREP7
DAMPRATIO=0.02!
振型阻尼系数都为0.02
ET,1,BEAM188
KEYOPT,1,7,1
KEYOPT,1,8,1
ET,2,MASS21,,,4
MP,EX,1,210E9
MP,NUXY,1,0.3
MP,DENS,1,7800
MPDATA,DAMP,1,,DAMPRATIO
SECTYPE,3,BEAM,I
SECDATA,0.2,0.2,0.5,0.01,0.01,0.008
R,1,1.6E2,1.6E2
R,2,1.2E2,1.2E2
L,1,2
L,2,3
L,3,4
L,4,5
LSEL,,,,1,4
LATT,1,0,1,,6,,3
ALLS
LESIZE,ALL,,,3,,,,,1
LMESH,ALL
NUMMRG,ALL
NUMCMP,ALL
TYPE,2
REAL,1
E,2
REAL,2
E,8
E,14
E,20
*SET,NT,1001
*SET,DT,0.02
*DIM,AC,,NT
*VREAD,AC,RECORD,txt,,IJK,NT,1,1,,
(F8.3)
/SOLU
D,1,ALL,
ANTYPE,0
SOLVE
/SOLU!
模态分析
ANTYPE,2
MODOPT,SUBSP,8
MXPAND,8,,,1
得到自振频率1
*GET,FREQ1,MODE,1,FREQ
/ESHAPE,1.0
/SOLU
ANTYPE,SPECTR!
谱分析
SPOPT,SPRS
MDAMP,1,DAMPRATIO,,,,,,
SED,1
SVTYPE,2,1
FREQ,0.167,0.25,0.333,0.44444,0.5,0.667,1,1.25,1.667
SV,,0.154350625,0.191590625,0.210210625,0.224175625,0.250716714,0.329514922,0.484352764,0.598723486,0.786897371
FREQ,2.222,10,100000
SV,,1.034212766,1.034212766,0.3528
按7度多遇地震,地震影响系数0.08第一组III场地TG=0.45
SRSS,0.01,DISP
/POST1
/INP,,MCOM
NSORT,U,SUM
*GET,DMAX_SPE,SORT,,MAX!
反应谱最大变形
开始地震计算
/CONFIG,NRES,20000
/config,fsplit,750这样的话在生成结果文件的时候,单个结果文件最大控制在3G
ANTYPE,TRANS
TRNOPT,FULL
ALPHAD,2*DAMPRATIO*FREQ1*2*3.1415926
BETAD,2*DAMPRATIO/(FREQ1*2*3.1415926)
*DO,I,1,1001
ACEL,AC(I),0,0
TIME,I*0.02
OUTRES,ALL,ALL
(注意:
若是输出东西格式太大,所以可以仔细利用输出项:
kbc,0
nsubst,2,10,1
设置子步数
outres,all,none !
输出控制
outres,nsol,last
outres,v,last
outres,a,last
outres,nload,last,
outres,strs,last
)
SOLVE
*ENDDO
FINISH
/POST26
FILE
NSOL,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
-0.0095
-0.012
-0.0142
-0.0128
-0.011
-0.0085
-0.0131
-0.0176
-0.0194
-0.0162
-0.0144
-0.0082
-0.0042
-0.0066
-0.019
-0.0196
0.003
0.0141
-0.0049
-0.0203
-0.026
-0.0325
-0.0306
-0.0172
-0.0197
-0.0163
-0.0164
-0.0067
0.0025
0.015
0.0236
0.0252
0.0336
0.0463
0.0492
0.0419
0.0359
0.0271
0.0235
0.0339
0.0412
0.053
0.0639
0.0732
0.0652
0.0599
0.04
0.0063
-0.0515
-0.0787
-0.0603
-0.0484
-0.025
-0.0059
0.0134
0.0308
0.0499
0.071
0.0995
0.1219
0.1529
0.1449
0.1155
0.0935
0.0892
0.0926
0.0839
0.0901
0.0993
0.1209
0.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.0017
0.036
0.0785
0.1164
0.1598
0.196
0.2412
0.2729
0.3036
0.32
0.3417
0.2821
0.2324
-0.1198
-0
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 地震 分析