动力学分析几个入门实例

lbt18872

适合初学的人研究之用，非常经典！

问题一：悬臂梁受重力作用发生大变形，求其固有频率。
  
基本过程：
1、建模
2、静力分析  
   NLGEOM,ON
   STRES,ON     
3、求静力解
4、开始新的求解：modal
   STRES,ON
   UPCOORD,1,ON  修正坐标  
   SOLVE...      
5、扩展模态解
6、察看结果   
  

1. /PREP7   
   
2. ET,1,BEAM189                          !使用beam189梁单元
   
3. MPTEMP,,,,,,,,   
   
4. MPTEMP,1,0   
   
5. MPDATA,EX,1,,210e9   
   
6. MPDATA,PRXY,1,,0.3   
   
7. MPDATA,DENS,1,,7850  
   
8. SECTYPE,   1, BEAM, RECT, secA, 0     !定义梁截面secA
   
9. SECOFFSET, CENT  
   
10. SECDATA,0.005,0.01,0,0,0,0,0,0,0,0   
    
11. K, ,,,,                               !建模与分网
    
12. K, ,2,,,
    
13. K, ,2,1,,   
    
14. LSTR,       1,       2   
    
15. LATT,1, ,1, ,       3, ,1   
    
16. LESIZE,1, , ,20, , , , ,1  
    
17. LMESH,       1   
    
18. FINISH   
    
19.   
    
20. /SOL                                  !静力大变形求解
    
21. ANTYPE,0
    
22. NLGEOM,1                              
    
23. PSTRES,ON                             !计及预应力效果
    
24. DK,1, , , ,0,ALL, , , , , ,   
    
25. ACEL,0,9.8,0,                         !只考虑重力作用
    
26. TIME,1   
    
27. AUTOTS,1
    
28. NSUBST,20, , ,1  
    
29. KBC,0   
    
30. SOLVE   
    
31. FINISH   
    
32.   
    
33. /SOLUTION   
    
34. ANTYPE,2                              !进行模态求解
    
35. MSAVE,0  
    
36. MODOPT,LANB,10   
    
37. MXPAND,10, , ,0                       !取前十阶模态
    
38. PSTRES,1                              !打开预应力效应
    
39. MODOPT,LANB,10,0,0, ,OFF
    
40. UPCOORD,1,ON                          !修正坐标以得到正确的应力
    
41. PSOLVE,TRIANG                         !三角化矩阵
    
42. PSOLVE,EIGLANB                        !提取特征值和特征向量
    
43. FINISH   
    
44. /SOLU   
    
45. EXPASS,1                              !扩展模态解
    
46. PSOLVE,EIGEXP   
    
47. FINISH   
    
48.   
    
49. /POST1   
    
50. SET,LIST                              !观察结果
    
51. FINISH  

复制代码

问题二：循环对称结构模态分析
这是ANSYS HELP里的例子，但那个命令流似乎有些问题，下面是整理过的命令流。
  
基本过程：
1、建模
2、define cyclic
3、定义约束
4、求模态解
5、展开并察看结果
  

1. r1=5                                 !建模
   
2. r2=10
   
3. d1=2
   
4. nsect=24
   
5. alpha_deg=360/nsect
   
6. alpha_rad=2*acos(-1)/nsect
   
7.   
   
8. /prep7
   
9. csys,1
   
10. k,1,0,0,0
    
11. k,2,0,0,d1
    
12. k,3,r1,0,0
    
13. k,4,r1,0,d1
    
14. l,3,4
    
15. arotat,1,,,,,,1,2,alpha_deg/2
    
16. k,7,r2,0,0
    
17. k,8,r2,0,d1
    
18. l,7,8
    
19. arotat,5,,,,,,1,2,alpha_deg/2
    
20. arotat,2,,,,,,1,2,alpha_deg/2
    
21. arotat,6,,,,,,1,2,alpha_deg/2
    
22. a,5,6,10,9
    
23. mshkey,1
    
24. et,1,181
    
25. r,1,0.20
    
26. r,2,0.1
    
27. mp,ex,1,10e6
    
28. mp,prxy,1,0.3
    
29. mp,dens,1,1e-4
    
30. esize,0.5
    
31. asel,,,,1,4
    
32. aatt,,1
    
33. asel,,,,5
    
34. aatt,,2
    
35. allsel
    
36.   
    
37. CYCLIC, , , ,'CYCLIC'                !定义基本扇区
    
38. AMESH,all                            !分网
    
39. FINISH   
    
40.   
    
41. /SOL                                 !求模态解
    
42. ANTYPE,2
    
43. MODOPT,LANB,10   
    
44. MXPAND,10, , ,0  
    
45. PSTRES,0
    
46. MODOPT,LANB,10,0,0, ,OFF
    
47. DK,5, , , ,0,UZ, , , , , ,   
    
48. SOLVE   
    
49. FINISH  
    
50.   
    
51. /POST1   
    
52. SET,LIST
    
53. SET, , ,1, , , ,4,   
    
54. /CYCEXPAND, ,ON  
    
55. PLNS,U,SUM                           !观察扩展模态
    
56. FINISH

复制代码

问题三：三角平台受谐波载荷作用,求结构响应
谐波载荷为作用在平台上面一点的集中力，幅值为10N，频率范围5Hz~70Hz
  
基本过程：
1、建模
2. 求模态解
3、用模态叠加法作谐分析
4. 观察结果
  

1. /PREP7                          !建模
   
2. ET,1,BEAM189
   
3. ET,2,SHELL93
   
4. R,1,0.01, , , , , ,  
   
5. MP,EX,1,210e9   
   
6. MP,PRXY,1,0.3   
   
7. MP,DENS,1,7850   
   
8. SECTYPE,   1, BEAM, RECT, secA, 0   
   
9. SECOFFSET, CENT  
   
10. SECDATA,0.005,0.008,0,0,0,0,0,0,0,0  
    
11. K, ,-0.5,,,  
    
12. K, ,0.5,,,   
    
13. K, ,,,1,
    
14. K, ,,1,1,   
    
15. K, ,-0.5,1,,
    
16. K, ,0.5,1,,  
    
17. K, ,,,0.5,   
    
18. A,4,5,6  
    
19. LSTR,       1,       5   
    
20. LSTR,       3,       4   
    
21. LSTR,       2,       6   
    
22. LSEL,S,LINE,,4,6
    
23. LATT,1,1,1, ,       7, ,1   
    
24. LSEL, , , ,ALL   
    
25. LESIZE,ALL, , ,10, , , , ,1  
    
26. LMESH,4,6   
    
27. TYPE,   2   
    
28. MSHAPE,0,2D  
    
29. MSHKEY,1
    
30. AMESH,1  
    
31. FINISH   
    
32.   
    
33. /SOL                     !为了使用模态叠加法谐分析
    
34. ANTYPE,2
    
35. MODOPT,LANB,20           !先取结构前20阶模态  
    
36. MXPAND,20, , ,0  
    
37. MODOPT,LANB,20,0,0, ,OFF
    
38. KSEL,S,KP,,1,3  
    
39. DK,ALL, , , ,0,ALL, , , , , ,   
    
40. SOLVE   
    
41. FINISH  
    
42.   
    
43. !/POST1                  !最好事先看一下模态结果
    
44. !SET,LIST                !以了解模态频率范围
    
45. !FINISH   
    
46.   
    
47. /SOL                     !模态叠加法谐响应分析
    
48. ANTYPE,3
    
49. HROPT,MSUP   
    
50. HROUT,ON
    
51. HROPT,MSUP,20, ,
    
52. HROUT,ON,OFF,0   
    
53. F,177,FY,-10,            !施加10N的力在节点177的负Y方向
    
54. HARFRQ,5,70,             !载荷的强制频率范围
    
55. NSUBST,200,  
    
56. KBC,1   
    
57. ALPHAD,5,   
    
58. SOLVE   
    
59. FINISH   
    
60.   
    
61. /POST26                 !察看位移响应
    
62. FILE,,rfrq
    
63. NUMVAR,20  
    
64. NSOL,2,139,U,Y,  
    
65. PLVAR,2
    
66. FINISH

复制代码

问题四：三角平台受一地震谱激励,求应力分布和支反力
  
基本过程：
1、建模
2、求模态解
3、求谱解
4、扩展模态
5、模态合并
6、观察结果
  

1. /PREP7             !建模
   
2. ET,1,BEAM189
   
3. ET,2,SHELL93
   
4. R,1,0.01, , , , , ,  
   
5. MP,EX,1,210e9   
   
6. MP,PRXY,1,0.3   
   
7. MP,DENS,1,7850  
   
8. SECTYPE,   1, BEAM, RECT, secA, 0   
   
9. SECOFFSET, CENT  
   
10. SECDATA,0.005,0.008,0,0,0,0,0,0,0,0  
    
11.    
    
12. K, ,-0.5,,,  
    
13. K, ,0.5,,,   
    
14. K, ,,,1,
    
15. K, ,,1,1,   
    
16. K, ,-0.5,1,,
    
17. K, ,0.5,1,,  
    
18. K, ,,,0.5,   
    
19. A,4,5,6  
    
20. LSTR,       1,       5   
    
21. LSTR,       3,       4   
    
22. LSTR,       2,       6   
    
23. LSEL,S,LINE,,4,6  
    
24. LATT,1,1,1, ,       7, ,1   
    
25. LSEL, , , ,ALL
    
26. LESIZE,ALL, , ,10, , , , ,1  
    
27. LMESH,4,6   
    
28. TYPE,   2   
    
29. MSHAPE,0,2D  
    
30. MSHKEY,1
    
31. AMESH,1  
    
32. FINISH   
    
33.   
    
34. /SOL              !取前十阶模态
    
35. ANTYPE,2
    
36. MODOPT,LANB,10   
    
37. KSEL,S,KP,,1,3  
    
38. DK,ALL, , , ,0,ALL, , , , , ,   
    
39. SOLVE   
    
40. FINISH   
    
41.   
    
42. /SOL                 !谱分析
    
43. ANTYPE,8
    
44. SPOPT,SPRS,10,1
    
45. SVTYP,3
    
46. SED,0,1,0,         !给出激励方向
    
47. FREQ,0.25,1.34,6.73,14.6,28.9,0,0,0,0    !激励谱
    
48. SV,0,0.00073,0.00016,0.00034,0.00034,0.00052,   
    
49. SOLVE   
    
50. FINISH   
    
51.   
    
52. /SOL                  !扩展模态
    
53. ANTYPE,2
    
54. EXPASS,1
    
55. MXPAND,10,0,0,1,0.001,   
    
56. SOLVE   
    
57. FINISH   
    
58.   
    
59. /SOL                 !模态合并
    
60. ANTYPE,8
    
61. SRSS,0.001,DISP  
    
62. SOLVE   
    
63.   
    
64. /POST1            !观察结果
    
65. SET,LIST  
    
66. /INPUT,,mcom
    
67. PRRSOL,F
    
68. FINISH

复制代码

问题五：三角平台受时程载荷作用,求应力分布和变形过程
  
基本过程：
1、建模
2. 施加随时间变化载荷，定义载荷步
3、求解
4. 在POST1和POST26中观察结果
  

1. /PREP7                             !建模
   
2. ET,1,BEAM189
   
3. ET,2,SHELL93  
   
4. R,1,0.01, , , , , ,  
   
5. MP,EX,1,210e9   
   
6. MP,PRXY,1,0.3   
   
7. MP,DENS,1,7850   
   
8. SECTYPE,   1, BEAM, RECT, secA, 0   
   
9. SECOFFSET, CENT  
   
10. SECDATA,0.005,0.008,0,0,0,0,0,0,0,0  
    
11. K, ,-0.5,,,  
    
12. K, ,0.5,,,   
    
13. K, ,,,1,
    
14. K, ,,1,1,   
    
15. K, ,-0.5,1,,
    
16. K, ,0.5,1,,  
    
17. K, ,,,0.5,   
    
18. A,4,5,6  
    
19. LSTR,       1,       5   
    
20. LSTR,       3,       4   
    
21. LSTR,       2,       6   
    
22. LSEL,S,LINE,,4,6
    
23. LATT,1,1,1, ,       7, ,1   
    
24. LSEL, , , ,ALL   
    
25. LESIZE,ALL, , ,10, , , , ,1  
    
26. LMESH,4,6   
    
27. TYPE,   2   
    
28. MSHAPE,0,2D  
    
29. MSHKEY,1
    
30. AMESH,1  
    
31. FINISH   
    
32.   
    
33. /SOL                             
    
34. ANTYPE,4
    
35. TRNOPT,FULL                     !完全法瞬态分析
    
36. KSEL,S,KP,,1,3
    
37. DK,ALL, , , ,0,ALL, , , , , ,   
    
38. OUTRES,ALL,1
    
39. ALPHAD,5,                       !α阻尼
    
40.   
    
41. TIME,2                          !定义载荷曲线
    
42. AUTOTS,1
    
43. NSUBST,50, , ,1      
    
44. KBC,0   
    
45. SFA,ALL,1,PRES,500  
    
46. LSWRITE,1,                       
    
47.    
    
48. TIME,3   
    
49. LSWRITE,2,   
    
50.   
    
51. TIME,4
    
52. SFA,ALL,1,PRES,150  
    
53. KBC,1   
    
54. LSWRITE,3,   
    
55.   
    
56. TIME,5
    
57. SFA,ALL,1,PRES,
    
58. LSWRITE,4,   
    
59.   
    
60. LSSOLVE,1,4,1,                  !求解
    
61. FINISH   
    
62.    
    
63. /POST26                         !观察变形随时间的变化
    
64. NSOL,2,177,U,Y,uyy   
    
65. PLVAR,2
    
66. FINISH

复制代码

[ 本帖最后由 yejet 于 2006-9-1 09:05 编辑 ]

dxtemp

例子都很简单，不知道售价为什么这么贵？难道是作者自己原创的？

lbt18872

不好意思哦，本来也就是不怎么难，这个价钱是随便定的，不知道高了，给个地址我给大家发过去吧！

多情清秋

原帖由 lbt18872 于 2006-8-31 23:46 发表
不好意思哦，本来也就是不怎么难，这个价钱是随便定的，不知道高了，给个地址我给大家发过去吧！

可以编辑修改价格，论坛不提倡出售帖子

vehicle

谢谢楼主，不错的例题！

sfx

还不错

navigator

谢谢了

wschaoren

不错，谢谢分享

hyl2323

ding,谢谢

dazhi9956

谢谢分享阿
不错的阿

enqi01

seems good

hello2001

谢谢，好东西啊！

三千里

顶一个，赚点钱

paopao

现在初学ANSYS,只会图上选,还不会编程呢,先顶上.