温馨提示！由于新浪微博认证机制调整，您的新浪微博帐号绑定已过期，请重新绑定！&&|&&
LOFTER精选
网易考拉推荐
用微信&&“扫一扫”
将文章分享到朋友圈。
用易信&&“扫一扫”
将文章分享到朋友圈。
本文转载自SAP2000
1.上部结构输入：约束结构基础节点，采用ACEL方式对所有质量点施加加速度时程；2.支座输入：释放地震波方向的约束，只针对基础节点施加加速度，同时为了避免结构出现整体的刚体位移导致数值计算不能进行，需要对基础节点施加“大质量‘，也就是常说的”大质量法“。事实上，两种方法在原理上都是一致的，即一致输入。只不过方法一固定基础节点，对上部结构施加加速度。可以这样理解，地震动本来是基础震动时程，现在约束了基础，只要对上部结构施加一个反方向的加速度时程，即可认为等同于基础输入，只不过得到的节点位移均为相对位移而非绝对位移。但是应力是准确的；方法二，真实的模拟了地震动的作用机理，得到的结构位移为绝对位移；但是由于”大质量“的引入，影响到了结构的阻尼比计算（若采用锐利阻尼的话），这样对结构到底能产生多大的影响还需要进一步研究。总体来讲，对一致输入，两种办法并没有什么对错之分，都是可以的。
阅读(135)|
用微信&&“扫一扫”
将文章分享到朋友圈。
用易信&&“扫一扫”
将文章分享到朋友圈。
历史上的今天
loftPermalink:'',
id:'fks_',
blogTitle:'【转载】地震时程分析中加速度的输入方式——以ansys为例',
blogAbstract:'d地震时程的输入常用的两类办法：上部结构输入和支座输入。（这里暂时不考虑行波效应等，均以一致输入为例）这两种办法反应到ansys中，实现有两种主流的做法：1.上部结构输入：约束结构基础节点，采用ACEL方式对所有质量点施加加速度时程；2.支座输入：释放地震波方向的约束，只针对基础节点施加加速度，同时为了避免结构出现整体的刚体位移导致数值计算不能进行，需要对基础节点施加“大质量‘，也就是常说的”大质量法“。事实上，两种方法在原理上都是一致的，即一致输入。只不过方法一固定基础节点，对上部结构施加加速度。可以这样理解，地震动本来是基础震动时程，现在约束了基础，只要对上部结构施加一个反方向的加速度时程，即可认为等同于基础输入，只不过得到的节点位移均为相对位移而非绝对位移。但是应力是准确的；',
blogTag:'',
blogUrl:'blog/static/',
isPublished:1,
istop:false,
modifyTime:1,
publishTime:0,
permalink:'blog/static/',
commentCount:0,
mainCommentCount:0,
recommendCount:1,
bsrk:-100,
publisherId:0,
recomBlogHome:false,
currentRecomBlog:false,
attachmentsFileIds:[],
groupInfo:{},
friendstatus:'none',
followstatus:'unFollow',
pubSucc:'',
visitorProvince:'',
visitorCity:'',
visitorNewUser:false,
postAddInfo:{},
mset:'000',
remindgoodnightblog:false,
isBlackVisitor:false,
isShowYodaoAd:false,
hostIntro:'',
hmcon:'0',
selfRecomBlogCount:'0',
lofter_single:''
{list a as x}
{if x.moveFrom=='wap'}
{elseif x.moveFrom=='iphone'}
{elseif x.moveFrom=='android'}
{elseif x.moveFrom=='mobile'}
${a.selfIntro|escape}{if great260}${suplement}{/if}
{list a as x}
推荐过这篇日志的人：
{list a as x}
{if !!b&&b.length>0}
他们还推荐了：
{list b as y}
转载记录：
{list d as x}
{list a as x}
{list a as x}
{list a as x}
{list a as x}
{if x_index>4}{break}{/if}
${fn2(x.publishTime,'yyyy-MM-dd HH:mm:ss')}
{list a as x}
{if !!(blogDetail.preBlogPermalink)}
{if !!(blogDetail.nextBlogPermalink)}
{list a as x}
{if defined('newslist')&&newslist.length>0}
{list newslist as x}
{if x_index>7}{break}{/if}
{list a as x}
{var first_option =}
{list x.voteDetailList as voteToOption}
{if voteToOption==1}
{if first_option==false},{/if}&&“${b[voteToOption_index]}”&&
{if (x.role!="-1") },“我是${c[x.role]}”&&{/if}
&&&&&&&&${fn1(x.voteTime)}
{if x.userName==''}{/if}
网易公司版权所有&&
{list x.l as y}
{if defined('wl')}
{list wl as x}{/list}新手来~地震时程反应分析实例~自己写的~【ansys吧】_百度贴吧
&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&签到排名：今日本吧第个签到，本吧因你更精彩，明天继续来努力！
本吧签到人数：0成为超级会员，使用一键签到本月漏签0次！成为超级会员，赠送8张补签卡连续签到：天&&累计签到：天超级会员单次开通12个月以上，赠送连续签到卡3张
关注：59,192贴子：
新手来~地震时程反应分析实例~自己写的~收藏
近日自己写了针对一个小框架的ansys计算实例，内容包含地震时程反应分析，地震反应谱分析，欢迎新手进来学习交流，高手勿喷！问题描述：一平面框架结构，设弹性模量为3.0e10Mpa，截面尺寸为0.2x0.3m，密度为2500kg/m^3；7度设防。采用Elcentro波对其进行地震时程反应分析；按照《建筑结构抗震规范》对其进行反应谱分析。这里先看第一问。按照规范，七度地区时程最大值为875px/s2,故先对波进行峰值调整，调整后的波形如下，最大值为35cm/s^2，时间间隔0.02s，这里选取前50s，也即是共有2500个记录数据调整后的波形图如下：
作为ANSYS中国区专业资质合作伙伴,20多年服务支持经验,我们提供结构/流体/电磁等一体化解决方案及工程咨询服务.咨询热线:
2、建模finish/clear/filname,earquake/prep7et,1,beam3tt=0.02N=2500
!数据记录间隔以及数据个数mp,ex,1,3e10mp,prxy,1,0.2mp,dens,1,2500r,1,0.6,0.k,1k,2,,4k,3,4,4k,4,4l,1,2l,2,3l,3,4esize,0.5lmesh,alldk,1,alldk,4,allfinish
3、读取地震加速度数据并运行分析/solualphad,3betad,0.001
!定义阻尼系数antype,transtimint,off
!首先做结构在重力荷载作用下的响应，关闭时间积分acel,,-9.8time,1e-6solve timint,on
!打开时间效应/input,output,txt*do,i,1,N,1acel,acc(i)time,tt*isolve*enddo
!地震时程计算
4、结果读取/post26nsol,2,14,u,x
!节点位移响应deriv,3,2,1
!节点速度响应deriv,4,3,1
!节点加速度响应xvar,1plvar,2
plvar,3plvar,4
位移响应曲线
后续会更新其他关于这个小框架的计算内容。自己顶一个，高手勿喷~~~
我也在做框架的时程分析，有时间可以交流一下,qq
不错，就是利用求到的 方式求加速度，求出来的值不准确。
楼主，按你这命令流运行不了啊
楼主你的地震加速度的时程数据哪里来的啊
楼主，如何对波进行峰值调整。我做的是刚架的地震响应分析，时程分析时用到天津波、EI centro波，但是我不会对波的峰值进行调整。好像调整和地震烈度和场地有关，求大神帮助！
欢迎加入达索系统举办的“体验时代的制造业”大会.-8日在上海W酒店,与智能制造领域专家讨论智能制造与工业4.0.更有CATIA,DELMIA等智能制造软件现场演示!
说好的更新呢
楼主 你的地震波是施加在整个框架上吗 不应该施加在框架底部吗
楼主，请教下，这个命令流最后结果读取要是需要y方向的加速度，需要修改哪些数值啊？
偶遇我大水哥
原来这是水哥开的贴呀，真是眼拙呀，。没看出来。
我想问一下你这个曲线怎么直接在ansys里面看啊？
楼主您好，最近也在做时程振动分析，但我对经典ANSYS不熟悉想询问您能否以ANSYSWORKBENCH做，又是以哪个模块做(transient structure)??
楼主，请教个问题，我正在做地铁车站的地震响应，计算过程中出现了警告：Shell element 72 references section 1 which is not a shell section.
will be associated with the real type ( 1 ).
但是我的壳单元181并没有定义截面啊，用得实常数定义的，定义的截面是梁单元188的，怎么会出现这个，搞不懂。我的命令流如下（可以看看问题在哪里吗）：/PREP7ET,1,SOLID65
模拟土体ET,2,SHELL181
模拟车站ET,3,BEAM188
模拟车站中柱R,1,0.027
车站顶板厚度R,2,0.014
车站中楼板厚R,3,0.030
车站底板厚R,4,0.020
车站侧壁厚SECTYPE,1,BEAM,RECT,,0
定义梁单元188的截面SECOFFSET,CENTSECDATA,0.026,0.026,0,0,0,0,0,0,0,0,0,0MP,EX,1,3E10MP,NUXY,1,0.2MP,DENS,1,2500MP,EX,2,15.3E6MP,NUXY,2,0.36MP,DENS,2,1910TB,DP,2,,,TBDATA,2,61E3,19.5,AA=0.656BB=0.442CC=0.014DD=0.020EE=0.026FF=0.266GG=0.027HH=0.030LL=2.261MM=0.3325KK=2.5NN=1.2K,1,(AA/2-DD/2),(BB/2-GG/2)K,2,-(AA/2-DD/2),(BB/2-GG/2)K,3,-(AA/2-DD/2),-(BB/2-GG/2)K,4,(AA/2-DD/2),-(BB/2-GG/2)K,5,(AA/2-DD/2)K,6,-(AA/2-DD/2)L,1,2L,2,3L,3,4L,4,1L,6,5K,7,AA,,,K,8,AA,,-LLL,7,8ADRAG,1,2,3,4,,,6ADRAG,5,,,,,,6LDELE,6,,,1AGLUE,1,2,3,4,5NUMCMP,ALLK,13,0,0,-MMK,14,0,(BB/2-GG/2),-MML,13,14LGEN,7,21,,,,,-FFLGEN,2,21,27,1,,-(BB/2-GG/2)K,1801,(AA/2-DD/2)K,1802,KK/2K,1803,KK/2,NN/2K,1804,(AA/2-DD/2),NN/2K,1805,,NN/2K,1806,,(BB/2-GG/2)K,1807,,(BB/2-GG/2),-LLK,1808,(AA/2-DD/2),(BB/2-GG/2)L,L,L,L,L,L,L,L,L,AL,35,38,37,36AL,39,40,41,42ASEL,S,,,8,9,1CM,AA2,AREAVDRAG,AA2,,,,,,43VSEL,S,,,1,2,1CM,VV1,VOLUVSYMM,X,VV1VSEL,S,,,1,4,1CM,VV2,VOLUVSYMM,Y,VV2ALLSEL,ALLAOVLAP,ALLVGLUE,ALLNUMMRG,KPNUMCMP,ALLLSEL,S,,,21,34,1CM,LL1,LINELATT,1,,3,,,,1LESIZE,LL1,,,5LMESH,ALLLSEL,ALLLSEL,S,,,6,7,1LSEL,A,,,42CM,LL2,LINELESIZE,LL2,0.0665LSEL,ALLLSEL,S,,,1,2,1LSEL,A,,,40,52,12CM,LL3,LINELESIZE,LL3,,,5ASEL,ALLASEL,S,,,2,14,12CM,AAA1,AREAAATT,1,1,2AMESH,AAA1LSEL,ALLLSEL,S,,,12,13,1CM,LL4,LINELESIZE,LL4,0.0665LSEL,ALLLSEL,S,,,3,11,8CM,LL5,LINELESIZE,LL5,,,10ASEL,ALLASEL,S,,,3AATT,1,2,2AMESH,3LSEL,ALLLSEL,S,,,5,9,4LSEL,A,,,74,76,2CM,LL6,LINELESIZE,LL6,,,5LSEL,ALLLSEL,S,,,8,10,2LSEL,A,,,75CM,LL7,LINELESIZE,LL7,0.0665ASEL,ALLASEL,S,,,24,30,6CM,AAA2,AREAAATT,1,3,2AMESH,AAA2LSEL,ALLLSEL,S,,,14,20,2LSEL,A,,,15,19,2LSEL,A,,,4CM,LL8,LINELESIZE,LL8,,,5ASEL,ALLASEL,S,,,4,7,1CM,AAA3,AREAAATT,1,4,2AMESH,AAA3ALLSEL,ALLSMRT,1VSEL,S,,,1,8,1CM,VV3,VOLUVATT,2,,1VSWEEP,VV3NUMMRG,KPNUMCMP,ALLFINISH/SOLUNSEL,ALLACEL,0,9.8,0,ASEL,S,,,21,22ASEL,A,,,28,35,7DA,ALL,UYASEL,S,,,10,18,8ASEL,A,,,31,34,3DA,ALL,UXASEL,S,,,1,8,7ASEL,A,,,16,26,10ASEL,A,,,27,33,6ASEL,A,,,36,37,1DA,ALL,UZFINISHALLSEL*DIM,SHUIPING,ARRAY,2,1000*CREATE,ANSUITMP*VREAD,SHUIPING(1,1),'SHUIPING','TXT','E:\ansys-kc\t1',1000,(E9.3,E11.3)*END/INPUT,ANSUITMP/SOLUANTYPE,4TRNOPT,FULLLUMPM,0*DO,T,1,50,1TIME,0.02*TKBC,1NSUB,1ACEL,SHUIPING(2,T)ALPHAD,0.069BETAD,0.035SOLVE*ENDDOFINISH
登录百度帐号推荐应用ansys中如何输入地震波?跪求一详细的命令流，附加读入地震波数据文件的格式。。谢谢了
全部答案（共1个回答）
试试把源安装文件的路径都改为英文，不要有中文出现。
众所周知，ansys有静态，谐波，瞬态三磁场种场求解方式。Bertotti提出过常系数三相式铁耗计算模型：
1：加载计算之后，可以看看单元节点间的应力差，如果这个差值占到这个单元最大应力的10%以上的时候，证明这个单元的尺寸比较大，应该在这个单元所在区域进行细分2：一...
地震波分为纵波和横波，纵波可以通过固体流体，横波只通过固体，两者随着介质的刚性和密度的变化而变速，根据地震波的传播数据，可以制成地球内部的震波传播速度曲线图，可...
答: 如何选择好的建筑防水材料？
答: 上网找
答: “鸟巢”外形结构主要由巨大的门式钢架组成，共有24根桁架柱，现已完成20根桁架柱整柱及2根下柱吊装。国家体育场建筑顶面呈鞍形，长轴为332.3米，短轴为296....
答: 建造师必须是通过国家资格监考人员进行查看
大家还关注
确定举报此问题
举报原因(必选)：
广告或垃圾信息
激进时政或意识形态话题
不雅词句或人身攻击
侵犯他人隐私
其它违法和不良信息
报告，这不是个问题
报告原因(必选)：
这不是个问题
这个问题分类似乎错了
这个不是我熟悉的地区在ANSYS里施加地震惯性力的方法
在ANSYS里做地震分析时，需要对结构施加地震惯性荷载，地震惯性力是通过加速度的方式输入进结构的，然后与结构的质量一起形成动力计算时的惯性荷载，下面说一下在ANSYS里施加地震惯性力的方法。
首先，将三个方向的地震加速度放到一个文本文件里，如accexyz.txt，在这个数据文件里共放三列数据，每列为一个方向的地震加速度值，这里仅给出数据文件中前几行的数据：
.......................
然后，再建一个文本文件用来存放三个方向的地震加速度时间点，如time.txt，在这个数据文件里仅一列数据，对应于加速度数据文件里每一行的时间点，这里给出数据文件中前几行数据：
.......................
编写如下的命令流文件，并命名为acce.inp
*dim,ACCEXYZ,TABLE,2000,3&&&&&&&&&&&&&&&&&&&&
*vread,ACCEXYZ(1,1),accexyz,txt,,JIK,3,2000&&
!02行(3e16.6)&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
*vread,ACCEXYZ(1,0),time,txt&&&&&&&&&&&&&&&&&
(e16.6)&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
ACCEXYZ(0,1)=1&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
ACCEXYZ(0,2)=2&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
!07行，同上
ACCEXYZ(0,3)=3&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&
!08行，同上
ANTYPE,trans&
btime=0.01&&&&&&&&&&&&&&
!定义计算起始时间
etime=15.00&&&&&&&&&&&&&&!定义计算结束时间
dtime=0.01&&&&&&&&&&&&&&
!定义计算时间步长
*DO,itime,btime,etime,dtime
time,itime
NSUBST,1, , ,1&
acel,ACCEXYZ(itime,1),ACCEXYZ(itime,2),ACCEXYZ(itime,3)&&&&
!施加三个方向的地震加速度
最后，在命令窗口里输入/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&&
2.9938”这样的格式，此行的数据格式也就修改为“(3f10.4)”这样的格式。
04行：从数据文件time.txt里读时间值到数组ACCEXYZ的第零列
06行：将数组ACCEXYZ的第零行赋值，如果不对行插值的话也可以不赋值
2、ANSYS中施加地震动荷载的讨论
已投稿到：
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
下载作业帮安装包
扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
ansys如何对位移时程取两次导数,来求出加速度曲线
扫二维码下载作业帮
拍照搜题，秒出答案，一键查看所有搜题记录
ANSYS里面好像没法求这个导数的,用matlab或者mathematics应该可以做的
为您推荐：
其他类似问题
扫描下载二维码}