ANSYS地震反应谱SRSS分析.docx

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 ! 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(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-方向谱分析 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(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-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(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-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 可是在X后作Y他不清空数据库，需要进行两次模态分析，这很耗时间对于大型结构，请大家讨论讨论如何处理呢？Re:讨论：ANSYS地震反应谱SRSS分析本人是学土木工程的，平时主要用Patran+Nastran对结构做线性分析，偶尔使用Ansys对结构做地震反应谱分析，但对Ansys的命令流不熟悉。依据本人的使用经验看，“一次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方向的反应谱内力计算。需要注意的是，在对话框中将SEDX，SEDY，SEDZ为1，1，0表示地震动的方向为与X轴成45度的方向，这显然不是所希望的Re:讨论：ANSYS地震反应谱SRSS分析我也知道是这样，可是他形成的*.mcom好像没有不是这个意思啊。按你说的，“一次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/COM,ANSYS RELEASE8.0 UP20030930 09:28:42 07/23/2005/COM, truss.mcom ！X方向LCOPER,ZEROLCDEFI,1, 1, 1LCFACT,1,0.263825E-17LCASE,1LCOPER,SQUARELCDEFI,1, 1, 2LCFACT,1, 8.55778 LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 3LCFACT,1, -0.188669E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 4LCFACT,1, -0.871099E-15LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 5LCFACT,1, -0.757013 LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 6LCFACT,1,0.967307E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 7LCFACT,1,0.533141E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 8LCFACT,1, -0.203699 LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 9LCFACT,1,0.445795E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 10LCFACT,1, -0.387808E-13LCOPER,ADD,1,MULT,1LCOPER,SQRT/COM,ANSYS RELEASE8.0 UP20030930 09:28:42 07/23/2005/COM, truss.mcom ！Y方向的LCOPER,SQUARE ！这里没有用zero清空LCDEFI,1, 1, 1LCFACT,1, 50.7528 LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 2LCFACT,1,0.887017E-14LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 3LCFACT,1,0.612824E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 4LCFACT,1,-1.96484 LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 5LCFACT,1, -0.331613E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 6LCFACT,1,0.330459E-13LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 7LCFACT,1,0.366569 LCOPER,ADD,1,MULT,1LCDEFI,1, 1, 8LCFACT,1, -0.976991E-13LCOPER,ADD,1,M