ANSYS地震反应谱SRSS分析共24页.docx
《ANSYS地震反应谱SRSS分析共24页.docx》由会员分享，可在线阅读，更多相关《ANSYS地震反应谱SRSS分析共24页.docx（13页珍藏版）》请在冰豆网上搜索。
ANSYS地震反应谱SRSS分析共24页
ANSYS地震反应谱SRSS分析
我在ANSYS中作地震分解反应谱分析，一次X方向，一次Y方向，他们要求是独立互不干扰的，可是采用直进行一次模态分析的话，他生成的*.mcom文件好像是包含了前面的计算结果，命令流如下：
!
进入PREP7并建模
/PREP7
B=15!
基本尺寸
A1=1000!
第一个面积
A2=1000!
第二个面积
A3=1000!
第三个面积
ET,1,beam4!
二维杆单元
R,1,0.25,0.0052,0.0052,0.5,0.5!
以参数形式的实参
MP,EX,1,2.0E11!
杨氏模量
mp,PRXY,1,,0.3
mp,dens,1,7.8e3
N,1,-B,0,0!
定义结点
N,2,0,0,0
N,3,-B,0,b
N,4,0,0,b
N,5,-B,0,2*b
N,6,0,0,2*b
N,7,-B,0,3*b
N,8,0,0,3*b
E,1,3!
定义单元
E,2,4
E,3,5
E,4,6
E,3,4
E,5,6
e,5,7
e,6,8
e,7,8
D,1,ALL,0,,2
FINISH
!
!
进入求解器，定义载荷和求解
/SOLU
D,1,ALL,0,,2!
结点UX=UY=0
sfbeam,1,1,PRES,100000,
sfbeam,3,1,PRES,100000,
sfbeam,7,1,PRES,100000,
SOLVE
FINISH
allsel
NMODE=10
/SOL
!
*
ANTYPE,2
!
*
MSAVE,0
!
*
MODOPT,LANB,NMODE
EQSLV,SPAR
MXPAND,NMODE,,,1
LUMPM,0
PSTRES,0
!
*
MODOPT,LANB,NMODE,0,0,,OFF
SOLVE
*DIM,FRE,,NMODE
*DO,I,1,NMODE
*GET,FRE（I）,MODE,I,FREQ!
OBTAINMODEFREQENCYFORMODEI
*ENDDO
FINISH
!
地震影响系数
grav=9.81
tg=0.35
amax=0.08
c=0.05
!
*dim,a,,nmode
*dim,t,,nmode
*do,i,1,nmode
t（i）=1.0/fre（i）
*enddo
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）
*do,i,1,nmode
*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
*enddo
!
!
X-方向谱分析SpectrumanalysisalongGlobalX-axisdirection
/SOLU
ANTYPE,SPECTR!
Spectrumanalysis
SPOPT,SPRS!
Singlepointspectrum
SED,1,,!
GlobalX-axisasspectrumdirection
SVTYP,2!
Seismicaccelerationresponsespectrum
!
FrequencypointsandSpectrumvaluesforSVvs.freq.table
FREQ,fre
（1）,fre
（2）,fre（3）,fre（4）,fre（5）,fre（6）,fre（7）,fre（8）,fre（9）
FREQ,fre（10）
SV,,a
（1）,a
（2）,a（3）,a（4）,a（5）,a（6）,a（7）,a（8）,a（9）
SV,,a（10）
FINISH
!
/SOLU
!
ANTYPE,MODAL!
Mode-frequencyanalysis
!
EXPASS,ON
!
MXPAND,nmode,,,YES,0.0!
Expandnmodeshapes,calculateelementstresses
!
SOLVE
!
FINISH
/SOLU
ANTYPE,SPECTR
SRSS,0.0,DISP!
SquareRootofSumofSquaresModecombination
!
withsignif=0.0anddisplacementsolutionrequested
SOLVE
FINISH
/POST1
SET,LIST
/INPUT,,mcom
!
***************EARTHQUAKEX******************
ALLSEL,ALL
FINISH
!
Y-方向谱分析SpectrumanalysisalongGlobalX-axisdirection
!
!
**********************************************
!
/SOL
!
!
*
!
ANTYPE,2
!
!
*
!
MSAVE,0
!
!
*
!
MODOPT,LANB,NMODE
!
EQSLV,SPAR
!
MXPAND,NMODE,,,1
!
LUMPM,0
!
PSTRES,0
!
!
*
!
MODOPT,LANB,NMODE,0,0,,OFF
!
SOLVE
!
FINISH
!
!
**********************************************
/SOLU
LSCLEAR,LSOPT
ANTYPE,SPECTR!
Spectrumanalysis
SPOPT,SPRS!
Singlepointspectrum
SED,,1,!
GlobalY-axisasspectrumdirection
SVTYP,2!
Seismicaccelerationresponsespectrum
FREQ
!
FrequencypointsandSpectrumvaluesforSVvs.freq.table
FREQ,fre
（1）,fre
（2）,fre（3）,fre（4）,fre（5）,fre（6）,fre（7）,fre（8）,fre（9）
FREQ,fre（10）
SV,,a
（1）,a
（2）,a（3）,a（4）,a（5）,a（6）,a（7）,a（8）,a（9）
SV,,a（10）
SOLVE
FINISH
!
/SOLU
!
ANTYPE,MODAL!
Mode-frequencyanalysis
!
EXPASS,ON
!
MXPAND,nmode,,,YES,0.0!
Expandnmodeshapes,calculateelementstresses
!
SOLVE
!
FINISH
/SOLU
ANTYPE,SPECTR
SRSS,0.0,DISP!
SquareRootofSumofSquaresModecombination
!
withsignif=0.0anddisplacementsolutionrequested
SOLVE
FINISH
/POST1
SET,LIST
/INPUT,,mcom
!
***************EARTHQUAKEY******************
ALLSEL,ALL
FINISH
这里在进行X方向的反应谱分析以后，进行Y方向的分析，可是他生成的*.mcom文件如下：
/COM,ANSYSRELEASE8.0UP2003093009:
28:
4207/23/2005
/COM,truss.mcom
LCOPER,ZERO
LCDEFI,1,1,1
LCFACT,1,0.263825E-17
LCASE,1
LCOPER,SQUARE
LCDEFI,1,1,2
LCFACT,1,8.55778
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,3
LCFACT,1,-0.188669E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,4
LCFACT,1,-0.871099E-15
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,5
LCFACT,1,-0.757013
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,6
LCFACT,1,0.967307E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,7
LCFACT,1,0.533141E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,8
LCFACT,1,-0.203699
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,9
LCFACT,1,0.445795E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,10
LCFACT,1,-0.387808E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
/COM,ANSYSRELEASE8.0UP2003093009:
28:
4207/23/2005
/COM,truss.mcom
LCOPER,SQUARE!
注意这里没有清空数据库
LCDEFI,1,1,1
LCFACT,1,50.7528
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,2
LCFACT,1,0.887017E-14
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,3
LCFACT,1,0.612824E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,4
LCFACT,1,-1.96484
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,5
LCFACT,1,-0.331613E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,6
LCFACT,1,0.330459E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,7
LCFACT,1,0.366569
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,8
LCFACT,1,-0.976991E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,9
LCFACT,1,0.417313E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,10
LCFACT,1,0.401040E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
我感觉这样好像是X和Y两个方向地震的叠加，
可是如果在座Y方向的地震以前把注释掉的模态分析在做一下这样的Y方向的地震的*.mcom就是：
/COM,truss.mcom
LCOPER,ZERO!
注意这里清空数据库
LCDEFI,1,1,1
LCFACT,1,50.7528
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,2
LCFACT,1,0.887017E-14
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,3
LCFACT,1,0.612824E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,4
LCFACT,1,-1.96484
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,5
LCFACT,1,-0.331613E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,6
LCFACT,1,0.330459E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,7
LCFACT,1,0.366569
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,8
LCFACT,1,-0.976991E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,9
LCFACT,1,0.417313E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,10
LCFACT,1,0.401040E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
如果在X方向后不作Y方向的地震，他的*.mcom：
/COM,ANSYSRELEASE8.0UP2003093008:
46:
2307/23/2005
/COM,truss.mcom
LCOPER,ZERO
LCDEFI,1,1,1
LCFACT,1,0.263825E-17
LCASE,1
LCOPER,SQUARE
LCDEFI,1,1,2
LCFACT,1,8.55778
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,3
LCFACT,1,-0.188669E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,4
LCFACT,1,-0.871099E-15
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,5
LCFACT,1,-0.757013
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,6
LCFACT,1,0.967307E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,7
LCFACT,1,0.533141E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,8
LCFACT,1,-0.203699
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,9
LCFACT,1,0.445795E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1,1,10
LCFACT,1,-0.387808E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
可是在X后作Y他不清空数据库，需要进行两次模态分析，这很耗时间对于大型结构，请大家讨论讨论如何处理呢？
Re:
讨论：
ANSYS地震反应谱SRSS分析
本人是学土木工程的，平时主要用Patran+Nastran对结构做线性分析，偶尔使用Ansys对结构做地震反应谱分析，但对Ansys的命令流不熟悉。
依据本人的使用经验看，“一次X方向，一次Y方向，他们要求是独立互不干扰的”应该这样完成才对：
1、点击Solution-〉LoadStepOpts-〉Spectrum-〉Settings，在对话框中将SEDX，SEDY，SEDZ设置为1，0，0后；点击Solution-〉Solve-〉CurrentLS；读入生成*.mcom文件完成X方向的反应谱内力计算。
2、点击Solution-〉LoadStepOpts-〉Spectrum-〉Settings，在对话框中将SEDX，SEDY，SEDZ设置为0，1，0后；点击Solution-〉Solve-〉CurrentLS；读入生成*.mcom文件完成Y方向的反应谱内力计算。
需要注意的是，在对话框中将SEDX，SEDY，SEDZ为1，1，0表示地震动的方向为与X轴成45度的方向，这显然不是所希望的
Re:
讨论：
ANSYS地震反应谱SRSS分析
我也知道是这样，可是他形成的*.mcom好像没有不是这个意思啊。
按你说的"，“一次X方向，一次Y方向，他们要求是独立互不干扰的”应该这样完成才对：
1、点击Solution-〉LoadStepOpts-〉Spectrum-〉Settings，在对话框中将SEDX，SEDY，SEDZ设置为1，0，0后；点击Solution-〉Solve-〉CurrentLS；读入生成*.mcom文件完成X方向的反应谱内力计算。
2、点击Solution-〉LoadStepOpts-〉Spectrum-〉Settings，在对话框中将SEDX，SEDY，SEDZ设置为0，1，0后；点击Solution-〉Solve-〉CurrentLS；读入生成*.mcom文件完成Y方向的反应谱内力计算。
"做后：
他把Y方向的同样写入到同一个*.mcom中了，而且也没有晴空数据库，形成如下的*.mcom
/COM,ANSYSRELEASE 8.0 UP20030930 09:
28:
42 07/23/2005
/COM,truss.mcom ！
X方向
LCOPER,ZERO
LCDEFI,1, 1, 1
LCFACT,1, 0.263825E-17
LCASE,1
LCOPER,SQUARE
LCDEFI,1, 1, 2
LCFACT,1, 8.55778
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 3
LCFACT,1,-0.188669E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 4
LCFACT,1,-0.871099E-15
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 5
LCFACT,1,-0.757013
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 6
LCFACT,1, 0.967307E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 7
LCFACT,1, 0.533141E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 8
LCFACT,1,-0.203699
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 9
LCFACT,1, 0.445795E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 10
LCFACT,1,-0.387808E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
/COM,ANSYSRELEASE 8.0 UP20030930 09:
28:
42 07/23/2005
/COM,truss.mcom ！
Y方向的
LCOPER,SQUARE ！
这里没有用zero清空
LCDEFI,1, 1, 1
LCFACT,1, 50.7528
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 2
LCFACT,1, 0.887017E-14
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 3
LCFACT,1, 0.612824E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 4
LCFACT,1, -1.96484
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 5
LCFACT,1,-0.331613E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 6
LCFACT,1, 0.330459E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 7
LCFACT,1, 0.366569
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 8
LCFACT,1,-0.976991E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 9
LCFACT,1, 0.417313E-13
LCOPER,ADD,1,MULT,1
LCDEFI,1, 1, 10
LCFACT,1, 0.401040E-13
LCOPER,ADD,1,MULT,1
LCOPER,SQRT
Re:
讨论：
ANSYS地震反应谱SRSS分析
本版精华贴中好像有关于地震分析的东东，各位可以找找看，希望有用。
Re:
讨论：
ANSYS地震反应谱SRSS分析
我也发现那个问题，我是利用/copy命令去解决的，模态分析完成后
下面命令流：
（0_1Ydizh.txt与0_1Zdizh.txt为y、z方向的谱分析命令流。
）
/COPY,file,db,,00,db
/COPY,file,full,,00,full
/COPY,file,esav,,00,esav
/COPY,file,mode,,00,mode
/COPY,file,rst,,00,rst
!
/COPY,00,db,,file,db
/COPY,00,full,,file,full
/COPY,00,esav,,file,esav
/COPY,00,mode,,file,mode
/COPY,00,rst,,file,rst
RESUME,file,db,,0,0
!
/INPUT,0_1Ydizh,txt,,,0
!
LCWRITE,30,fY,L01,
/COPY,file,mcom,,fY,mcom
!
/CLEAR,START
/COPY,00,db,,file,db
/COPY,00,full,,file,full
/COPY
升级会员 