! 7.10  矩形截面薄壁懸臂梁的自由振動
! 本程序來源于邢靜忠編著的《ANSYS應用實例與分析》，科學出版社，2006年
!
FINISH
/CLEAR, NOSTART	! 清除數據, 不加載初始化文件
! （1）定義工程選項、單元類型、材料參數和實常數
/FILNAME, EX7.10	! 定義工程文件名稱
/PREP7	! 加載前處理模塊
/TITLE, EX7.10, TRANSVERSE VIBRATIONS OF A SHEAR BEAM.
ANTYPE, MODAL	! 指定分析類型為模態分析
ET, 1, SHELL28	! 使用考慮剪切和扭曲效應的殼單元SHELL28
R, 1, 0.1  	! 定義殼體厚度為0.1m
MP, EX, 1, 200E9	! 定義第1類材料的彈性模量EX=200GPa
MP, NUXY, 1, 0.27	! 定義第1類材料的泊松比NUXY=0.27
MP, DENS, 1, 7860.0	! 定義第1類材料的密度7860kg/m^3
! （2）定義節點位置和單元連接
/VIEW, 1, 1, 1, 1	! 改變窗口1的視點為等軸測位置
N, 1  	! 在(0, 0, 0)定義結點1
N, 2, 10 	! 在(10m, 0, 0)定義結點2
N, 3, 10, 10	! 在(10m, 10m, 0)定義結點3
N, 4, , 10	! 在(0, 10m, 0)定義結點4
NGEN, 9, 4, 1, 4, 1, , , 3.75	! 循環9次, 每次節點號增量4, 用前面定義的1號到4號節點模式, 增量為1, Z坐標增量3.75m, 自動生成其余節點
E, 1, 2, 6, 5	! 過節點1, 2, 6, 5定義單元
EGEN, 3, 1, -1	! 循環3次, 按照前面單元模式, 節點號增加1, 生成3個單元（包括前面一個單元）
E, 4, 1, 5, 8	! 過節點4, 1, 5, 8定義單元
EGEN, 8, 4, -4	! 按循環8次, 每次節點號增加4, 按照前面定義的4個單元模式, 生成其余單元
! （3）設置位移約束和耦合自由度
D, ALL, UZ 	! 限制所有節點的沿著Z方向的自由度UZ
D, 1, UX, , , 4, 1, UY	! 限制1號到4號節點的沿著X和Y方向的平移自由度UX和UY
NSEL, S, LOC, Z, 3.75	! 選擇位于Z=3.75米位置上的所有節點
CP, 1, UX, ALL	! 將選擇集中的節點在X方向的自由度UX耦合為一個自由度, 定義為集合1
CP, 2, UY, ALL	! 將選擇集中的節點在Y方向的自由度UY耦合為一個自由度, 定義為集合2
NSEL, S, LOC, Z, 7.5	! 選擇位于Z=7.5米位置上的所有節點
CP, 3, UX, ALL	! 將選擇集中的節點在X方向的自由度UX耦合為一個自由度, 定義為集合3
CP, 4, UY, ALL	! 將選擇集中的節點在Y方向的自由度UY耦合為一個自由度, 定義為集合4
NSEL, S, LOC, Z, 11.25	! 選擇位于Z=11.25米位置上的所有節點
CP, 5, UX, ALL	! 將選擇集中的節點在X方向的自由度UX耦合為一個自由度, 定義為集合5
CP, 6, UY, ALL	! 將選擇集中的節點在Y方向的自由度UY耦合為一個自由度, 定義為集合6
NSEL, S, LOC, Z, 15	! 選擇位于Z=15米位置上的所有節點
CP, 7, UX, ALL	! 將選擇集中的節點在X方向的自由度UX耦合為一個自由度, 定義為集合7
CP, 8, UY, ALL	! 將選擇集中的節點在Y方向的自由度UY耦合為一個自由度, 定義為集合8
NSEL, S, LOC, Z, 18.75	! 選擇位于Z=18.75米位置上的所有節點
CP, 9, UX, ALL	! 將選擇集中的節點在X方向的自由度UX耦合為一個自由度, 定義為集合9
CP, 10, UY, ALL	! 將選擇集中的節點在Y方向的自由度UY耦合為一個自由度, 定義為集合10
NSEL, S, LOC, Z, 22.5	! 選擇位于Z=22.5米位置上的所有節點
CP, 11, UX, ALL	! 將選擇集中的節點在X方向的自由度UX耦合為一個自由度, 定義為集合11
CP, 12, UY, ALL	! 將選擇集中的節點在Y方向的自由度UY耦合為一個自由度, 定義為集合12
NSEL, S, LOC, Z, 26.25	! 選擇位于Z=26.25米位置上的所有節點
CP, 13, UX, ALL	! 將選擇集中的節點在X方向的自由度UX耦合為一個自由度, 定義為集合13
CP, 14, UY, ALL	! 將選擇集中的節點在Y方向的自由度UY耦合為一個自由度, 定義為集合14
NSEL, S, LOC, Z, 30	! 選擇位于Z=30米位置上的所有節點
CP, 15, UX, ALL	! 將選擇集中的節點在X方向的自由度UX耦合為一個自由度, 定義為集合15
CP, 16, UY, ALL	! 將選擇集中的節點在Y方向的自由度UY耦合為一個自由度, 定義為集合16
NSEL, ALL	! 選擇所有結點
FINISH	! 退出前處理模塊PREP7
! （4）指定求解選項, 并求解
/SOLU	! 加載求解模塊
MXPAND, 4	! 展開前面兩階重復的4個模態
MODOPT, REDUC	! 用降階方法計算特征值
TOTAL, 8	! 設置自動主自由度生成到8
SOLVE	! 開始求解
! （5）提取頻率計算結果
*GET, F1, MODE, 1, FREQ	! 獲取第1階頻率, 存入變量F1
*GET, F2, MODE, 3, FREQ	! 獲取第3階頻率, 存入變量F2
*STATUS, PARM	! 列表顯示所有內存變量的值
FINISH    	! 退出求解模塊SOLU
! （6）在一般后處理模塊中, 顯示模態形狀, 列表顯示內存變量的值
/POST1	! 進入一般后處理模塊
SET, 1, 1	! 讀入第1個荷載步的計算結果(第1階振型)
/VUP, , Z	! 設置Z坐標軸向上
/VIEW, , , 1	! 從（0, 0, 1）位置向（0, 0, 0）方向觀察模型
PLDISP, 1	! 圖形顯示第1階模態
SET, 1, 3	! 讀入第3個荷載步的計算結果(第3階振型)
PLDISP, 1	! 圖形顯示第2階模態
FINISH	! 退出后處理模塊