查看: 1121|回复: 3|关注: 0
Matlab 地震波的编写 希望高手指点
<h1 style="color:# 麦片财富积分
新手, 积分 5, 距离下一级还需 45 积分
如附件所示蓝色区域所示的2个公式中，i属于x,J属于y,n为快照时间属于t，差分法化公式已经化好，我现在只知道要3重循环好像，但要根据那2个公式画出波场快照图，心有余而力不足，希望各位帮忙下，万分感谢。
15:14 上传
点击文件名下载附件
1.54 MB, 下载次数: 11629
地震波2个公式和图像
<h1 style="color:# 麦片财富积分
帮顶:lol :lol :lol
<h1 style="color:# 麦片财富积分
求解答啊！！！
<h1 style="color:# 麦片财富积分
寻求帮助啊:Q
站长推荐 /2
利用MATLAB进行传感器数据分析
MATLAB中文论坛是全球最大的 MATLAB & Simulink 中文社区。用户免费注册会员后，即可下载代码，讨论问题，请教资深用户及结识书籍作者。立即注册加入我们吧！
MATLAB官方社交平台
MATLAB中文论坛微社区建好模型以后
outr,all,all
anty,trans
*do,n,1,444
acel,,aa(n)
lssolve,1,444
接着后处理
关于如何定义及给
数组赋值参见
UtilityMenu——&Parameters——&ArrayParameters——&Define/Edit
UtilityMenu——&Parameters——&ArrayParameters——&ReadFrom
、注意：若输入的格式不正确，会导致
、弹出菜单的含义为（例如要输入一个
矩阵的第一个元素
Fname:1.txt/*
Ncol:5/*aa
矩阵的列数
列，指定从数据文件的每行读
个数据赋予矩
赋值完毕以后，点
UtilityMenu——&Parameters——&ArrayParameters
——&Define/Edit
已被赋值为
(window.slotbydup=window.slotbydup || []).push({
id: '4540225',
container: s,
size: '910,250',
display: 'inlay-fix'
在此可输入您对该资料的评论~
(window.slotbydup=window.slotbydup || []).push({
id: '4540180',
container: s,
size: '250,200',
display: 'inlay-fix'
资料评价：
所需积分：1君，已阅读到文档的结尾了呢~~
ANSYS地震作用分析应用实例解析&#x2d;&#x2d;楼房地震波瞬态分析
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
ANSYS地震作用分析应用实例解析--楼房地震波瞬态分析
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer--144.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口(30) 完成结构自振分析后，下面就可以进行反应谱分析
(31) 进入ANSYS主菜单Solution-&Analysis Type-&New Analysis，选择分析类型为Spectrum。
(32) 在ANSYS主菜单Solution-&Analysis Type-&Analysis Options中，选择谱分析的类型为单点输入(Sing-pt resp)
（33) 在ANSYS主菜单中选择Preprocessor-&Loads-&Load Step Opts-& Time/Frequenc-& Damping，输入所有的阻尼为DAMPRATIO，如图所示 （
34) 下面需要定义地震的反应谱。我国规范给定的是基于加速度的反应谱。在ANSYS主菜单中选择Solution-& Load Step Opts-&Spectrum-&Single Point-&Settings，设定反应谱类型为地震加速度，放大系数为1，输入方向为X方向（1，0，0）。
(35) 接下来开始输入地震反应谱。这里输入的反应谱按7度多遇地震，取地震影响系数为0.08，第一组，III类场地，卓越周期Tg=0.45s。值得注意的是，我国规范给的反应谱横坐标是周期，ANSYS定义的反应谱横坐标是频率，应该注意上述区别。选择ANSYS主菜单Solution-& Load Step Opts-&Spectrum-&Single Point-&Freq Table，输入频率反应谱Freq1~Freq12为0.167, 0.25, 0.333, 0.4, 0.667, 1, 1.25, 1.667, 2.222, 10, 100000, 如图
(36) 选择ANSYS主菜单Solution-& Load Step Opts-&Spectrum-&Single Point-&Spectr Values，输入对应的反应谱数值依次如下：
(37) 最后选择ANSYS主菜单Solution-& Load Step Opts-&Spectrum-&Single Point-&Mode Combine，设定振型组合方式为SRSS法，如图所示
(38) 选择ANSYS主菜单Solution-&Solve-&Current LS，计算反应谱结果
(39) 进入ANSYS主菜单后处理模块General Postproc，在ANSYS窗口顶部菜单选择File-&Read input from，选择文件后缀名为*.mcom的文件
(40) 进入ANSYS主菜单General Postproc-&Plot Results-&Deformed Shapes，选择绘制变形后形状和结构形状，得到地震反应谱分析的结构变形如图
(41) 最后我们来进行地震时程分析，进行地震时程分析以前，首先要有一个地震时程记录，本例子给定的地震时程记录总长20秒，记录点间隔0.02s，共有1001个记录点。该地震记录存放在RECORD.TXT文件中。 (42) 首先建立两个变量，在ANSYS窗口顶部菜单选择Parameters-&Scalar parameters，在窗口中输入NT=1001，即总共1001个记录点，DT=0.02，即记录点间隔0.02s
(43) 然后还需要定义一个数组来存放地震记录，选择ANSYS窗口顶部菜单Parameters-&Array parameters-&Define/Edit，在弹出的Array Parameters窗口中点击Add按钮，在Add New Array Parameter窗口中输入数组的名称为AC，行数为NT个，如图所示。
(44) 下面从数据文件中读入地震时程记录。选择ANSYS窗口顶部菜单Parameters-&Array parameters-&Read from file，在弹出窗口中依次输入将数据读入的数组名称AC，可以用Browse选择目标数据文件，最后要给出数据读入的格式（F8.3），可以参考Fortran相应的文件输入输出要求。
(45) 完成数据读入后就可以进行地震反应分析。由于本次地震分析一共要进行1001步，超过了ANSYS默认的最大1000步限制，因此首先需要进行调整。在ANSYS的命令输入窗口内，输入以下命令： Finish /CONFIG, NRES, 20000
(46) 进入ANSYS主菜单Solution-&Analysis Type-&New Analysis，指定分析类型为瞬态分析Transient，Solution method 可以选择Full。
(47) 进入ANSYS主菜单Solution-&Analysis Type-&Sol&#39;n Controls，在Solution Controls窗口里面选择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的时程后处理器TimeHist Postpro来处理计算结果，在ANSYS主菜单上选择TimeHist Postpro，这时会弹出时程变量窗口Time History Variables。点击工具栏上第一个绿色加号按钮，添加时程变量。在Add Time-History Variable窗口中选择要添加的变量为节点的X方向位移。
(50) 选择第20号节点，回到时程变量窗口Time History Variables，点击工具栏上第三个按钮，即绘制出节点的位移时程曲线如图
*SET,NT,500
假如有500个地震波数据点
*SET,DT,0.02
每个地震波数据点之间的时间间隔是0.02
*dim,AC1,,NT,1
!定义X方向加速波的数组
!*dim,AC2,,NT,1
!定义Y、Z方向加速波的数组 如果有的话
!*dim,AC3,,NT,1
*vread,AC1(1),elcentro_EW,txt,,JIK,1,NT !读入X方向地震波文件(elcentro_EW.txt)到数组
!*vread,AC2(1),elcentro_NS,txt,,JIK,1,NT !读入Y、Z方向地震波文件、如果有的话
!*vread,AC3(1),elcentro_EW,txt,,JIK,1,NT
*DO,I,1,NT
ACEL,AC1(I),AC2(I),AC3(I)
!************读入地震波数据************
*DIM,ELC,ARRAY,1,150,0, , ,
*CREATE,ansuitmp
*VREAD,ELC(1,1),&#39;ELC&#39;,&#39;txt&#39;,&#39; &#39;,150, , , , , ,
(e9.3,e11.3)
/INPUT,ansuitmp
!****************求解****************
!eqslv,pcg,1e-4
!指定分析类型为瞬态动力学分析
TRNOPT,full
!瞬态动力学分析采用模态叠加法 !**************地震载荷**************
*DO,T,1,150,1
!循环读入地震数据并求解
TIME,0.001*T
!设置时间步
!指定载荷为递增载荷
!设定子步数为1
DAMPRATIO=0.05
!结构的阻尼比
FREQ1=1.0721*2*pi
!角频率变圆频率
freq2=1.0927*2*pi
!通过模态求解得出频率
ALPHAD,2*DAMPRATIO*FREQ1*freq2/(freq1+freq2)
!质量阻尼 BETAD,2*DAMPRATIO/(FREQ1+freq2)
ACEL,ELC(1,T)/347,,
!设定x，y方向加速度 ALLSEL
!选中所有元素
建立acelx、acely两个地震波文本（txt），文本中第一列为时间，第二列为加速度，如果两列之间为两空格，则加速度为负，一个空格则加速度为正，190行，19秒，间隔0.1记录一个值
麻烦你看看错在哪，谢谢
*dim,tjx,array,190,2
*dim,tjy,array,190,2
*create,ansuitmp
*vread,tjx(1,1),ACELX,txt,,jik,2,190
(E9.3,E11.3)
/input,ansuitmp
*create,ansuitmp
*vread,tjy(1,1),ACELY,txt,,jik,2,190
(e9.3,e11.3)
/input,ansuitmp
!对于地震波的输入，可以把荷载记录做成文件，利用apdl的读取功能读入倒数据库中。下面的例子是自己编的一个小文件。修改一下可以更简洁。有用到的朋友自己作一下把。
/config,nres,1000
*dim,aceX,TABLE,3000,1
*dim,aceY,TABLE,3000,1
*dim,aceZ,TABLE,3000,1
*vread,aceX(1,1),acex,txt,,1
*vread,aceX(1,0),ACETT,,,1
ACEX(0,1)=1
*vread,aceY(1,1),acey,txt,,1
*vread,aceY(1,0),ACETT,,,1
ACEY(0,1)=1
*vread,aceZ(1,1),acez,txt,,1
*vread,aceZ(1,0),ACETT,,,1
ACEZ(0,1)=1
！地震波时程记录分成了3个文件，每个文件是一列。分别记录x,y,z方向的加速度。acett是时间记录。
这样就可以把加速度记录读取倒ansys数据库中作为数组。
也可以把加速度记录做成一个文件，这样程序就简单多了。大家可以试看看修改一下。资源分类：
简介：在ansys 中如何 施加 地震波_法律资料_人文社科_专业资料。三向输入简化后的单向...U4 }&#39; T3 ~ L4 d7 Z !对于地震波的输入,可以把荷载记录做成文件,利用 ......
简介：ANSYS地震波的输入_信息与通信_工程科技_专业资料 暂无评价|0人阅读|0次下载|举报文档 ANSYS地震波的输入_信息与通信_工程科技_专业资料。三向输入简化后的单向......
简介：对于地震波的输入,可以把荷载记录作成文件,利用 apdl 的读取功能读入数据 库中...这样就可以把加速度记录读取到 ansys 数据库中作为数组。 也可以把加速度记录作......
简介：使用ANSYS软件模拟地震荷载的方法_计算机软件及应用_IT/计算机_专业资料。使用 ANSYS...模拟地震荷载首先需要有地震的加速度数据,这里采用迁安波,迁安波为渤 海的地震......
简介：ANSYS地震波的输入_法律资料_人文社科_专业资料
暂无评价|0人阅读|0次下载|举报文档ANSYS地震波的输入_法律资料_人文社科_专业资料。三向输入简化后的单向输入 ......
简介：ansys载荷分析命令_互联网_IT/计算机_专业资料 暂无评价|0人阅读|0次下载|举报文档 ansys载荷分析命令_互联网_IT/计算机_专业资料。/filname,plane42 /title,......
简介：ansys 施加正弦载荷的例子 /PREP7 ET,1,SHELL63 R,1,0.1 MP,EX,1,2e11 MP,NUXY,1,0.3 MP,DENS,1,7800 BLC4,,,10,2 AATT,1,1,1,0 AESIZE,......
简介：一般分为以下两类:基于结构自振特性的地震反应谱分析和基于特定地震波的地震时程...(ANSYS) 土木工程中除了常见的静力分析以外, 动力分析, 特别是结构在地震荷载......
简介：Ansys 载荷步叠加验证———浮云 152 为验证 ANSYS 是否将载荷步一与载荷步二的力进行叠加计算的, 先分别加载对其进行验证分析。 对比验证步骤: 对比 1:载荷为......
简介：ansys 极限载荷实例时间: 19:41 来源:未知 作者:蚂蚁 点击: 527 次 极限分析认为在一次加载 (静载)情况下整个承载截面都达到屈服极限时结构所 能......
copyright@ 2013- Inc. All Rights Reserved 果果文库 版权所有 联系站长: ; 经营许可证编号:浙ICP备号