ANSYS地震反应谱SRSS分析

ANSYS地震反应谱SRSS分析
我的具体操作采用“一次X方向，一次Y方向，他们要求是独立互不干扰的”应该这样完成才对：1、点击Solution-〉Load Step Opts-〉Spectrum-〉Settings，在对话框中将SEDX，SEDY，SEDZ设置为1，0，0后；点击Solution-〉Solve-〉Current LS；读入生成*.mcom文件完成X方向的反应谱内力计算。2、点击Solution-〉Load Step Opts-〉Spectrum-〉Settings，在对话框中将SEDX，SEDY，SEDZ设置为0，1，0后；点击Solution-〉Solve-〉Current LS；读入生成*.mcom文件完成Y方向的反应谱内力计算。"做后：他把Y方向的同样写入到同一个*.mcom中了，而且也没有晴空数据库，形成如下的*.mcom
所以这样操作后的*.mcom包含两个方向的作用。

所以为了制作一次模态分析就可以计算2个方向的地震的单独作用，我的办法是“把最初的模态分析的结果copy，等到第一次地震分析结束后再copy回来进行第二次分析”。

不知道大侠们有何高见？

下面就是我在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,MODE,I,FREQ ! OBTAIN MODE FREQENCY FOR MODE I
*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=1.0/fre
*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,ge,0.0,and,t,lt,0.1,then
a=(0.45+(10.0*p2-4.5)*t)*amax*grav
*elseif,t,ge,0.1,and,t,le,tg
a=p2*amax*grav
*elseif,t,gt,tg,and,t,le,5*tg
a=(tg/t)**r*p2*amax*grav
*else
a=(p2*0.2**r-p1*(t-5*tg))*amax*grav
*endif
*enddo
!
! X-方向谱分析 Spectrum analysis along Global X-axis direction

/SOLU
ANTYPE,SPECTR ! Spectrum analysis
SPOPT,SPRS ! Single point spectrum
SED,1,, ! Global X-axis as spectrum direction
SVTYP,2 ! Seismic acceleration response spectrum

! Frequency points and Spectrum values for SV vs. freq. table
FREQ,fre(1),fre(2),fre(3),fre(4),fre(5),fre,fre(7),fre,fre(9)
FREQ,fre(10)
SV,,a(1),a(2),a(3),a(4),a(5),a,a(7),a,a(9)
SV,,a(10)
FINISH

!/SOLU
!ANTYPE,MODAL ! Mode-frequency analysis
!EXPASS,ON
!MXPAND,nmode,,,YES,0.0 ! Expand nmode shapes, calculate element stresses
!SOLVE
!FINISH

/SOLU
ANTYPE,SPECTR
SRSS,0.0,DISP ! Square Root of Sum of Squares Mode combination
! with signif=0.0 and displacement solution requested
SOLVE
FINISH

/POST1
SET,LIST
/INPUT,,mcom
!***************EARTHQUAKE X******************
ALLSEL,ALL
FINISH

! Y-方向谱分析 Spectrum analysis along Global X-axis direction
!!**********************************************
!/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 ! Spectrum analysis
SPOPT,SPRS ! Single point spectrum
SED,,1, ! Global Y-axis as spectrum direction
SVTYP,2 ! Seismic acceleration response spectrum
FREQ

! Frequency points and Spectrum values for SV vs. freq. table
FREQ,fre(1),fre(2),fre(3),fre(4),fre(5),fre,fre(7),fre,fre(9)
FREQ,fre(10)
SV,,a(1),a(2),a(3),a(4),a(5),a,a(7),a,a(9)
SV,,a(10)
SOLVE
FINISH

!/SOLU
!ANTYPE,MODAL ! Mode-frequency analysis
!EXPASS,ON
!MXPAND,nmode,,,YES,0.0 ! Expand nmode shapes, calculate element stresses
!SOLVE
!FINISH

/SOLU
ANTYPE,SPECTR
SRSS,0.0,DISP ! Square Root of Sum of Squares Mode combination
! with signif=0.0 and displacement solution requested
SOLVE
FINISH

/POST1
SET,LIST
/INPUT,,mcom
!***************EARTHQUAKE Y******************
ALLSEL,ALL
FINISH
这里在进行X方向的反应谱分析以后，进行Y方向的分析，可是他生成的*.mcom文件如下：
/COM,ANSYS RELEASE 8.0 UP20030930 09:28:42 07/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,ANSYS RELEASE 8.0 UP20030930 09:28:42 07/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,ANSYS RELEASE 8.0 UP20030930 08:46:23 07/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

楼主您好，能在桥头堡看到您写的这篇资料感觉很有帮助，但是到了最后您输入不同方向的反应谱产生的结果文件，您也贴出来了，感觉您也比较疑惑，按照规范的要求，尤其是跨度比较大的结构，DX,DY,DDZ三个方向的地震引起的振动都是需要考虑的，如何在ansys的计算结果中体现出来呢？谢谢您！

这个有点太贵了吧，大概说明一下有什么内容？

看了你的帖子，既然是和别人交流，为什么要收那么贵？

刚来，没有钱买啊，
只要是自己需要的东东，花点钱也没什么吗，最好不要太贵啊

楼主,虽然是好东西,但是太贵了,买不起

这个有点太贵了吧，大概说明一下有什么内容？

楼主应该要回答别人的问题吧，既然步能免费提供，至少也应给个目录什么的是吧？

这么贵，我可买不起啊

哈哈，感觉这是我见过的最另类的论坛

什么东西这么贵啊

我是穷人啊，先收藏吧＞

好贵呀,楼主
买不起
赚点钱再说吧

:(太贵了吧，都是穷人

是太贵了点的吧  太TM狠了

很无聊，我不买，买不起气压