[转帖]ansys分析的几个经典例子

AaronSpark

悬臂梁受重力作用发生大变形，求其固有频率。<BR>本过程： <BR>1、建模 <BR>2、静力分析 <BR>NLGEOM,ON <BR>PSTRES,ON <BR>3、求静力解 <BR>4、开始新的求解：modal <BR>PSTRES,ON <BR>UPCOORD,1,ON 修正坐标 <BR>PSOLVE... <BR>5、扩展模态解 <BR>6、察看结果 <BR><BR>/PREP7 <BR>ET,1,BEAM189 !使用beam189梁单元 <BR>MPTEMP,,,,,,,, <BR>MPTEMP,1,0 <BR>MPDATA,EX,1,,210e9 <BR>MPDATA,PRXY,1,,0.3 <BR>MPDATA,DENS,1,,7850 <BR>SECTYPE, 1, BEAM, RECT, secA, 0 !定义梁截面secA <BR>SECOFFSET, CENT <BR>SECDATA,0.005,0.01,0,0,0,0,0,0,0,0 <BR>K, ,,,, !建模与分网 <BR>K, ,2,,, <BR>K, ,2,1,, <BR>LSTR, 1, 2 <BR>LATT,1, ,1, , 3, ,1 <BR>LESIZE,1, , ,20, , , , ,1 <BR>LMESH, 1 <BR>FINISH <BR><BR>/SOL !静力大变形求解 <BR>ANTYPE,0 <BR>NLGEOM,1 <BR>PSTRES,ON !计及预应力效果 <BR>DK,1, , , ,0,ALL, , , , , , <BR>ACEL,0,9.8,0, !只考虑重力作用 <BR>TIME,1 <BR>AUTOTS,1 <BR>NSUBST,20, , ,1 <BR>KBC,0 <BR>SOLVE <BR>FINISH <BR><BR>/SOLUTION <BR>ANTYPE,2 !进行模态求解 <BR>MSAVE,0 <BR>MODOPT,LANB,10 <BR>MXPAND,10, , ,0 !取前十阶模态 <BR>PSTRES,1 !打开预应力效应 <BR>MODOPT,LANB,10,0,0, ,OFF <BR>UPCOORD,1,ON !修正坐标以得到正确的应力 <BR>PSOLVE,TRIANG !三角化矩阵 <BR>PSOLVE,EIGLANB !提取特征值和特征向量 <BR>FINISH <BR>/SOLU <BR>EXPASS,1 !扩展模态解 <BR>PSOLVE,EIGEXP <BR>FINISH <BR><BR>/POST1 <BR>SET,LIST !观察结果 <BR>FINISH<BR>

AaronSpark

问题二：循环对称结构模态分析 <BR><BR>这是ANSYS HELP里的例子，但那个命令流似乎有些问题，下面是整理过的命令流。 <BR><BR>基本过程： <BR>1、建模 <BR>2、define cyclic <BR>3、定义约束 <BR>4、求模态解 <BR>5、展开并察看结果 <BR><BR>r1=5 !建模 <BR>r2=10 <BR>d1=2 <BR>nsect=24 <BR>alpha_deg=360/nsect <BR>alpha_rad=2*acos(-1)/nsect <BR><BR>/prep7 <BR>csys,1 <BR>k,1,0,0,0 <BR>k,2,0,0,d1 <BR>k,3,r1,0,0 <BR>k,4,r1,0,d1 <BR>l,3,4 <BR>arotat,1,,,,,,1,2,alpha_deg/2 <BR>k,7,r2,0,0 <BR>k,8,r2,0,d1 <BR>l,7,8 <BR>arotat,5,,,,,,1,2,alpha_deg/2 <BR>arotat,2,,,,,,1,2,alpha_deg/2 <BR>arotat,6,,,,,,1,2,alpha_deg/2 <BR>a,5,6,10,9 <BR>mshkey,1 <BR>et,1,181 <BR>r,1,0.20 <BR>r,2,0.1 <BR>mp,ex,1,10e6 <BR>mp,prxy,1,0.3 <BR>mp,dens,1,1e-4 <BR>esize,0.5 <BR>asel,,,,1,4 <BR>aatt,,1 <BR>asel,,,,5 <BR>aatt,,2 <BR>allsel <BR><BR>CYCLIC, , , ,'CYCLIC' !定义基本扇区 <BR>AMESH,all !分网 <BR>FINISH <BR><BR>/SOL !求模态解 <BR>ANTYPE,2 <BR>MODOPT,LANB,10 <BR>MXPAND,10, , ,0 <BR>PSTRES,0 <BR>MODOPT,LANB,10,0,0, ,OFF <BR>DK,5, , , ,0,UZ, , , , , , <BR>SOLVE <BR>FINISH <BR><BR>/POST1 <BR>SET,LIST <BR>SET, , ,1, , , ,4, <BR>/CYCEXPAND, ,ON <BR>PLNS,U,SUM !观察扩展模态 <BR>FINISH<BR><BR>

AaronSpark

问题三：三角平台受谐波载荷作用,求结构响应 <BR>谐波载荷为作用在平台上面一点的集中力，幅值为10N，频率范围5Hz~70Hz <BR><BR>基本过程： <BR>1、建模 <BR>2. 求模态解 <BR>3、用模态叠加法作谐分析 <BR>4. 观察结果 <BR><BR>/PREP7 !建模 <BR>ET,1,BEAM189 <BR>ET,2,SHELL93 <BR>R,1,0.01, , , , , , <BR>MP,EX,1,210e9 <BR>MP,PRXY,1,0.3 <BR>MP,DENS,1,7850 <BR>SECTYPE, 1, BEAM, RECT, secA, 0 <BR>SECOFFSET, CENT <BR>SECDATA,0.005,0.008,0,0,0,0,0,0,0,0 <BR>K, ,-0.5,,, <BR>K, ,0.5,,, <BR>K, ,,,1, <BR>K, ,,1,1, <BR>K, ,-0.5,1,, <BR>K, ,0.5,1,, <BR>K, ,,,0.5, <BR>A,4,5,6 <BR>LSTR, 1, 5 <BR>LSTR, 3, 4 <BR>LSTR, 2, 6 <BR>LSEL,S,LINE,,4,6 <BR>LATT,1,1,1, , 7, ,1 <BR>LSEL, , , ,ALL <BR>LESIZE,ALL, , ,10, , , , ,1 <BR>LMESH,4,6 <BR>TYPE, 2 <BR>MSHAPE,0,2D <BR>MSHKEY,1 <BR>AMESH,1 <BR>FINISH <BR><BR>/SOL !为了使用模态叠加法谐分析 <BR>ANTYPE,2 <BR>MODOPT,LANB,20 !先取结构前20阶模态 <BR>MXPAND,20, , ,0 <BR>MODOPT,LANB,20,0,0, ,OFF <BR>KSEL,S,KP,,1,3 <BR>DK,ALL, , , ,0,ALL, , , , , , <BR>SOLVE <BR>FINISH <BR><BR>!/POST1 !最好事先看一下模态结果 <BR>!SET,LIST !以了解模态频率范围 <BR>!FINISH <BR><BR>/SOL !模态叠加法谐响应分析 <BR>ANTYPE,3 <BR>HROPT,MSUP <BR>HROUT,ON <BR>HROPT,MSUP,20, , <BR>HROUT,ON,OFF,0 <BR>F,177,FY,-10, !施加10N的力在节点177的负Y方向 <BR>HARFRQ,5,70, !载荷的强制频率范围 <BR>NSUBST,200, <BR>KBC,1 <BR>ALPHAD,5, <BR>SOLVE <BR>FINISH <BR><BR>/POST26 !察看位移响应 <BR>FILE,,rfrq <BR>NUMVAR,20 <BR>NSOL,2,139,U,Y, <BR>PLVAR,2 <BR>FINISH<BR>

AaronSpark

问题四：三角平台受一地震谱激励,求应力分布和支反力 <BR><BR><BR>基本过程： <BR>1、建模 <BR>2、求模态解 <BR>3、求谱解 <BR>4、扩展模态 <BR>5、模态合并 <BR>6、观察结果 <BR><BR>/PREP7 !建模 <BR>ET,1,BEAM189 <BR>ET,2,SHELL93 <BR>R,1,0.01, , , , , , <BR>MP,EX,1,210e9 <BR>MP,PRXY,1,0.3 <BR>MP,DENS,1,7850 <BR>SECTYPE, 1, BEAM, RECT, secA, 0 <BR>SECOFFSET, CENT <BR>SECDATA,0.005,0.008,0,0,0,0,0,0,0,0 <BR><BR>K, ,-0.5,,, <BR>K, ,0.5,,, <BR>K, ,,,1, <BR>K, ,,1,1, <BR>K, ,-0.5,1,, <BR>K, ,0.5,1,, <BR>K, ,,,0.5, <BR>A,4,5,6 <BR>LSTR, 1, 5 <BR>LSTR, 3, 4 <BR>LSTR, 2, 6 <BR>LSEL,S,LINE,,4,6 <BR>LATT,1,1,1, , 7, ,1 <BR>LSEL, , , ,ALL <BR>LESIZE,ALL, , ,10, , , , ,1 <BR>LMESH,4,6 <BR>TYPE, 2 <BR>MSHAPE,0,2D <BR>MSHKEY,1 <BR>AMESH,1 <BR>FINISH <BR><BR>/SOL !取前十阶模态 <BR>ANTYPE,2 <BR>MODOPT,LANB,10 <BR>KSEL,S,KP,,1,3 <BR>DK,ALL, , , ,0,ALL, , , , , , <BR>SOLVE <BR>FINISH <BR><BR>/SOL !谱分析 <BR>ANTYPE,8 <BR>SPOPT,SPRS,10,1 <BR>SVTYP,3 <BR>SED,0,1,0, !给出激励方向 <BR>FREQ,0.25,1.34,6.73,14.6,28.9,0,0,0,0 !激励谱 <BR>SV,0,0.00073,0.00016,0.00034,0.00034,0.00052, <BR>SOLVE <BR>FINISH <BR><BR>/SOL !扩展模态 <BR>ANTYPE,2 <BR>EXPASS,1 <BR>MXPAND,10,0,0,1,0.001, <BR>SOLVE <BR>FINISH <BR><BR>/SOL !模态合并 <BR>ANTYPE,8 <BR>SRSS,0.001,DISP <BR>SOLVE <BR><BR>/POST1 !观察结果 <BR>SET,LIST <BR>/INPUT,,mcom <BR>PRRSOL,F <BR>FINISH<BR>

AaronSpark

问题五：三角平台受时程载荷作用,求应力分布和变形过程<BR><BR><BR>基本过程： <BR>1、建模 <BR>2. 施加随时间变化载荷，定义载荷步 <BR>3、求解 <BR>4. 在POST1和POST26中观察结果 <BR><BR>/PREP7 !建模 <BR>ET,1,BEAM189 <BR>ET,2,SHELL93 <BR>R,1,0.01, , , , , , <BR>MP,EX,1,210e9 <BR>MP,PRXY,1,0.3 <BR>MP,DENS,1,7850 <BR>SECTYPE, 1, BEAM, RECT, secA, 0 <BR>SECOFFSET, CENT <BR>SECDATA,0.005,0.008,0,0,0,0,0,0,0,0 <BR>K, ,-0.5,,, <BR>K, ,0.5,,, <BR>K, ,,,1, <BR>K, ,,1,1, <BR>K, ,-0.5,1,, <BR>K, ,0.5,1,, <BR>K, ,,,0.5, <BR>A,4,5,6 <BR>LSTR, 1, 5 <BR>LSTR, 3, 4 <BR>LSTR, 2, 6 <BR>LSEL,S,LINE,,4,6 <BR>LATT,1,1,1, , 7, ,1 <BR>LSEL, , , ,ALL <BR>LESIZE,ALL, , ,10, , , , ,1 <BR>LMESH,4,6 <BR>TYPE, 2 <BR>MSHAPE,0,2D <BR>MSHKEY,1 <BR>AMESH,1 <BR>FINISH <BR><BR>/SOL <BR>ANTYPE,4 <BR>TRNOPT,FULL !完全法瞬态分析 <BR>KSEL,S,KP,,1,3 <BR>DK,ALL, , , ,0,ALL, , , , , , <BR>OUTRES,ALL,1 <BR>ALPHAD,5, !α阻尼 <BR><BR>TIME,2 !定义载荷曲线 <BR>AUTOTS,1 <BR>NSUBST,50, , ,1 <BR>KBC,0 <BR>SFA,ALL,1,PRES,500 <BR>LSWRITE,1, <BR><BR>TIME,3 <BR>LSWRITE,2, <BR><BR>TIME,4 <BR>SFA,ALL,1,PRES,150 <BR>KBC,1 <BR>LSWRITE,3, <BR><BR>TIME,5 <BR>SFA,ALL,1,PRES, <BR>LSWRITE,4, <BR><BR>LSSOLVE,1,4,1, !求解 <BR>FINISH <BR><BR>/POST26 !观察变形随时间的变化 <BR>NSOL,2,177,U,Y,uyy <BR>PLVAR,2 <BR>FINISH<BR>

qilong

hao

luckyriver

不错！

usefully

好东西，顶上去！

xiven2008

加分鼓励! 命令流的实例贴希望多发一些，这样就可以互相学习<BR><BR>[em01][em01]

漓江烟雨

UP~~~~~~~~~~~~

zkuihua

xinyuxf

谢谢分享！还没用过LANB方法呢，我试试

喂蜗牛

谢谢