ANSYS地震反应谱SRSS分析.docx

ANSYS地震反应谱SRSS分析.docx

《ANSYS地震反应谱SRSS分析.docx》由会员分享，可在线阅读，更多相关《ANSYS地震反应谱SRSS分析.docx（13页珍藏版）》请在冰豆网上搜索。

ANSYS地震反应谱SRSS分析

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,00,