ANSYS入門教程加載求解及后處理技術B

1、ANSYS 入門教程 - 加載、求解及后處理技術B2011-01-09 14:20:38|  分類： ansys |  標簽： |字號大中小 訂閱 三、 施加面載荷 (續(xù))3.   在線上施加面荷載    命令：SFL, LINE, Lab, VALI, VALJ, VAL2I, VAL2J        LINE - 擬施加荷載的線號，也可為 ALL 或組件名。     &#

2、160;  Lab - 面荷載標識符，結構分析為 PRES。        VALI - 線始端關鍵點處的面荷載值，也可為表格型邊界條件的表格名。        VALJ - 線末端關鍵點處的面荷載值，也可為表格型邊界條件的表格名。如為空（缺?。┡c VALI 相等，否則采用輸入數據。        VAL2I,VAL2J - 為復數輸入時的虛部，而 VALI 和 VALJ 則為

3、實部。    該命令僅對 2D 面單元的邊界（線）、軸對稱單元本身、殼單元邊界（線）有效，對 3D 實體單元的線無效。    對于 2D 面單元，其輸入的面荷載值為“力/面積”；而對殼單元，其輸入的面荷載值為“力/長度”，這點需要特別注意。4.   在面上施加面荷載    命令：SFA, AREA, LKEY, Lab, VALUE, VALUE2        AREA - 擬施加面荷載的面號，也可為ALL或元件名

4、。        LKEY - 荷載施加的面號 (缺省為 1)。如果面為體單元的表面，則 LKEY 將被忽略；對殼單元 LKEY 可取 1 或 2，而其它值無效，                    單元幫助中有詳細說明。        Lab - 面荷載標

5、識符，結構分析為 PRES。        VALUE - 面荷載值，也可為表格名稱。        VALUE2 - 對結構分析無意義。    該命令對殼單元和 3D 體單元的面施加法向面荷載，對 2D  面單元無效。5. 面荷載梯度及其加載    定義面荷載梯度后，可在 SF、SFE、SFL 和 SFA 命令中使用。如前所述，SFE、SFL 及 SFBEAM 命令可以直接施加線性分布

6、荷載，采用 SFGRAD 命令定義荷載梯度后，使用 SF 和 SFA 命令施加線性分布荷載比較方便，如靜水壓力、圓柱體分布壓力等。    示例：        !  EX4.9  利用荷載梯度在直角坐標系下的施加方法        FINISH  $ /CLEAR  $ /PREP7        et,1,82

7、60; $ blc4,10,60                 ! 定義單元類型，創(chuàng)建面        esize,2  $ amesh,all                  

8、;   ! 定義單元尺寸，劃分網格        /PSF,PRES,NORM,2                    ! 設置荷載顯示方式        sfgrad,pres,y,0,-5    

9、0;                  ! 定義荷載梯度，SLZER=0，沿 Y 正方向遞減 5 單位/長度        nsel,s,loc,x,0                

10、60;              ! 選擇 X=0，且 Y=040 的節(jié)點群        nsel,r,loc,y,0,40        sf,all,pres,600             &

11、#160;               ! 對節(jié)點群施加面荷載，基值（Y=SLZER=0 處）為 600                ! 上述荷載施加結果：Y=0 處為 600，Y=40 處為 600+(40-0)×(-5)=400     

12、60;  sfgrad,pres,y,30,-20                  ! 再重新定義梯度荷載，SLZER=30，斜率為 -20        nsel,s,loc,x,10              &#

13、160;              ! 選擇 X=10 的節(jié)點群        sf,all,pres,0                         &#

14、160;       ! 對節(jié)點群施加面荷載，基值（Y=SLZER=30 處）為 0        allsel  $ eplot                           

15、60;   ! 如圖 4-1 所示 6.   表面效應單元施加面荷載    如前所述，施加具有 LKEY 參數的面荷載與單元類型相關，        對于 2D 面單元僅可在單元邊上或邊界上施加平行于單元面的荷載；        對于 3D 體單元，僅可施加單元面法向面荷載；        對于 3

16、D 殼單元，可施加單元面法向面荷載和在單元邊上或邊界上施加平行于單元面的荷載。    但有時所要施加的荷載不屬于上述情況，例如面的切向荷載或其它非法向面荷載等，此時可使用表面效應單元覆蓋所要施加荷載的表面，并用這些單元作為“管道”施加所需荷載。    如 2 D  面單元和 3 D  單元可分別使用 SURF153  單元和 SURF154 單元。    示例：        finish  $

17、/clear  $ /prep7        et,1,solid95  $ et,2,surf154         ! 定義 SOLID95 單元和表面效應單元 surf 154        blc4,10,10,40           &

18、#160;                ! 創(chuàng)建長方體        esize,5  $ vmesh,all                     ! 定

19、義單元尺寸，劃分網格        /psf,pres,tany,2                           ! 設置壓力顯式方式（單元坐標系 Y 切向）        nsel,s,lo

20、c,y,10                             ! 選擇 Y=10 的所有節(jié)點        type,2  $ esurf       

21、                      ! 設置單元類型 2, 生成表面效應單元        esel,s,type,2                &

22、#160;              ! 選擇單元類型為 2 的單元（表面效應單元）        sfe,all,3,pres,100                       

23、; ! 施加 LKEY=3 的面荷載（切向）        allsel                                       

24、     ! 可查看單元荷載（求解過程與 SURF154 無關）  四、  施加體載荷    在結構分析中，ANSYS 的體荷載只有溫度，其標識符為 TEMP。    常用體載荷命令如下表：    幾個主要的體荷載施加命令如下：        BF, NODE, Lab, VAL1      

25、60; BFE, ELEM, Lab, STLOC, VAL1, VAL2, VAL3, VAL4        BFK, KPOI, Lab, VAL1        BFL, LINE, Lab, VAL1        BFA, AREA, Lab, VAL1       

26、 BFV, VOLU, Lab, VAL1    其使用方法與面荷載施加命令類似，例如第 1 個參數均為圖素編號，也可為 ALL 或組件名；第 2 個參數 Lab = Ui、ROTi (i = x 或 y 或 z)、TEMP 或 FLUE；VAL1VAL4 為體荷載值，其中 VAL2VAL4 為單元不同位置上的體荷載值；STLOC 為 VAL1 指定一個對應的起始位置。 五、  施加慣性荷載    慣性荷載有加速度、角速度和角加速度。     慣

27、性荷載沒有刪除命令，要刪除慣性荷載，衷心將荷載值設為 0 即可；且慣性載荷為斜坡荷載。    ACEL、OMEGA 和 DOMEGA 命令分別用于施加在總體直角坐標系中的加速度、角速度和角加速度。需要注意的是 ACEL 命令施加的是加速度不是重力場，而重力加速度的方向與重力方向相反。因此要施加一個 -Y 方向的重力場，必須施加一個 Y 方向的加速度。    使用 CGOMGA 和 DCGOMG 命令定義一轉動物體的加速度和角加速度，但為相對于參考坐標系轉動時的物理量（該物體繞參考坐標系轉動）。CGLOC 命令用于指定參考坐標系相對于整

28、個笛卡爾坐標系的位置。    CMOMEGA 和 CMDOMGA 命令在單元元件上施加參考坐標系下的角速度和角加速度。    ANSYS 定義的三種類型轉動為：          整個結構繞總體直角坐標系轉動（OMEGA 和 DOMEGA 命令輸入）；          單元元件繞參考坐標系軸的轉動（CMOMEGA 和 CMDOMEGA 命令）；  

29、0;        整體直角坐標系繞加速度原點的轉動（ CGOMGA、DCGOMG 和 CGLOC命令）。    以上三種類型轉動中，可兩兩組合同時施加到結構上。    此處僅介紹 ACEL 命令及使用方法，命令如下：    命令：ACEL,  ACELX,  ACELY,  ACELZ        其中 ACELX, ACELY, ACELZ

30、分別為總體直角坐標系 X 軸、Y 軸和 Z 軸的結構線加速度值。 六、  施加耦合場荷載    在耦合場分析中，通常將包含一個分析中的結果施加在第二個分析中作為荷載，例如可將熱分析中計算得到的節(jié)點溫度，作為體積荷載施加到結構分析中，形成耦合場荷載。    施加耦合場荷載的命令為LDREAD命令，該命令是從一個結果文件讀出數據然后作為荷載施加到模型上。因此該命令不僅僅在施加耦合場荷載時使用，也可用于其它分析目的，例如可用于結構分析中讀入反作用力作為進一步分析的荷載等。    命令：LDR

31、EAD, Lab, LSTEP, SBSTEP, TIME, KIMG, Fname, Ext 七、  初應力荷載及施加    初應力（Initial Stress）可以指定為一種“荷載”進行施加，但僅在靜態(tài)分析和全瞬態(tài)分析中可以使用，可以用于線性分析或非線性分析。初應力荷載只能在第一個荷載步中施加。    ANSYS 中支持初應力荷載的單元類型有：PLANE2、PLANE42、PLANE82、PLANE182、PLANE183、SOLID45、SOLID92、SOLID95、SOLID185、SOLID

32、186、SOLID187、SHELL181、SHELL208、SHELL209、LINK180、BEAM188、BEAM189。    初應力荷載是單元坐標系下的值，如果單元坐標系與總體坐標系不同應謹慎。初應力荷載只能在求解層施加。初應力荷載的施加采用覆蓋方式，即多次施加時后面命令結果覆蓋前面的命令結果。    初應力荷載施加在被選擇的單元上，如果單元選擇集為空或不選擇某些單元，則不施加初應力荷載到這些單元上。    主要初應力載荷命令如下表： 1.   施加初始常應力荷載&

33、#160;   命令：ISTRESS,  Sx,  Sy,  Sz,  Sxy,  Syz,  Sxz,  MAT1,  MAT2,  MAT3,  MAT4,  MAT5,  MAT6,  MAT7,  MAT8,  MAT9,  MAT10        Sx,Sy,Sz,Sxy,Syz,Sxz - 初始的常應力值。  &#

34、160;     MAT1MAT10 - 初應力擬施加到的材料編號，如沒有指定，則施加到所有材料上。    該命令對所選擇的單元施加一組初始常應力值。2. 從文件施加初應力荷載    命令：ISFILE, Option, Fname, Ext, -, LOC, MAT1, MAT2, MAT3, MAT4, MAT5, MAT6, MAT7, MAT8, MAT9, MAT10        Option - 初應力荷載操作控

35、制參數，其值可?。?#160;           =READ（缺省）：從文件讀入初應力數據；            =LIST：列出已經讀入的初應力；            =DELE：刪除已經讀入的初應力。      

36、;  Fname - 當 Option=READ 時，F(xiàn)name 為一目錄和文件名。                      當 Option=LIST 或 DELE時，F(xiàn)name 為列表或刪除單元編號上的初應力。        Ext - 文件擴展名或層號，當 Fname 為空時，Ext 缺省為“IST”

37、。                    如 Option=LIST 或 DELE 則 Ext 為層殼單元的層號。        LOC - 總體位置標志，確定每個單元內初應力要施加的位置，其值可?。?#160;           =0（

38、缺?。涸趩卧|心上施加初應力；            =1：單元積分點上施加初應力；            =2：在單元指定位置上施加初應力。即由初應力文件確定將初應力荷載施加到什么位置，此時各個單元施加的位置可以不相同。            =3：常應力狀態(tài)。用初

39、應力文件中的第一個應力數據將所有單元初始化為一個常應力。        MAT1MAT10 - 初應力擬施加到的材料編號。    該命令對所選擇的單元施加初應力荷載，初應力的單元號與所選擇的單元號相對應。 3.   生成初應力文件    命令：ISWRITE, Switch    其中 Switch 參數控制初應力文件是否生成文件，其中可?。?#160;  

40、     ON：以工作文件及擴展名 IST 生成初應力文件，并寫入數據；        OFF：不生成初應力文件。    該命令僅在求解層有效，如果已有同名文件存在則覆蓋之。    該命令不支持 CDWRITE 命令。    用 ISWRITE 命令寫出的應力為單元積分點應力，對于非線性分析，寫入的應力數據為收斂后應力；對于線性分析，為求解完成后的

41、應力。因此其初應力文件標志區(qū)數據總為 eis,elemno,1，其中 elemno 為單元號，而 1 表示為積分點應力的位置標識。在用 ISFILE 命令讀入時，如果位置標志為 0，則采用各單元的第一個應力記錄；如果位置標志為 2，則采用初應力文件中的位置標志 （即1）；如果位置標志為 3，則采用應力文件的第一個應力數據。    示例：        !  EX4.12  初應力荷載      &

42、#160; finish  $ /clear  $ /filname,colu1  $ /prep7                 ! 定義工作文件名為 colu1        et,1,plane82  $ mp,ex,1,2e5  $ mp,nuxy,1,0.3 

43、0;        ! 定義單元類型和材料屬性        blc4,1,10  $ esize,2  $ amesh,all                         &

44、#160;   ! 創(chuàng)建面，定義網格尺寸，劃分網格        nsel,s,loc,y,0  $ d,all,uy  $ d,1,ux                            

45、60; ! 施加約束條件        nsel,s,loc,y,10  $ sf,all,pres,-10                                ! 施加

46、節(jié)點面荷載        allsel  $ finish        /solu  $ iswrite,on                         

47、0;                            ! 進入求解層，打開初應力文件生成開關        solve  $ finish        

48、                                                  

49、  ! 求解生成初應力文件（在當前工作目錄中）     Colu1.ist 文件格式及內容          ! 為說明問題，這里重新建模        finish  $ /clear  $ /filname,colu2  $ /prep7         

50、                ! 定義工作文件名為 colu2        et,1,plane82   $ mp,ex,1,2e5  $ mp,nuxy,1,0.3             &#

51、160;   ! 定義單元類型和材料屬性        blc4,1,10  $ esize,2  $ amesh,all                             

52、        ! 創(chuàng)建面，定義網格尺寸，劃分網格        nsel,s,loc,y,0  $ d,all,uy  $ d,1,ux  $ alls                      &

53、#160;    ! 施加約束條件        /solu                                   

54、                                              ! 進入求解層   

55、0;    loc=2                                             

56、;                                   ! 定義位置參數，改變此參數可得到不同的加載效果        isfile,read,colu1

57、,ist,loc                                                

58、0;    ! 從文件 colu1.ist 中讀入初應力并作為荷載施加        isfile,list                                

59、;                                             ! 查看施加的初應力荷載   

60、0;    solve                                             

61、;                                    !求解并可查看結果     初應力荷載不是施加“應力歷史”而是一種“荷載”。因此對于用 ISWRITE 命令生成的初應力文件，再用 IS

62、FILE 命令讀入后，當僅有初應力荷載時，其效果是模型中應力為零而位移與原荷載產生的位移反向。    要消除由于初應力荷載引起的位移且保持模型中應力不變，可將原荷載一并施加，此時模型中應力與原荷載產生的應力相同，但位移場為零（位移很小，可認為是零）。    一懸臂梁在端部受集中力作用，先生成初應力文件；然后再施加初應力荷載和集中力；計算后可得到荷載作用下的應力場但無位移場。    示例：        !  EX4.13 懸臂

63、梁荷載作用下的應力場但無位移場        finish  $ /clear  $ /filname,cant1  $ /prep7                          ! 定義工作文件名等    

64、60;   et,1,PLANE42  $ mp,ex,1,2e5  $ mp,nuxy,1,0.3                 ! 定義單元類型、材料特性等        blc4,10,1  $ esize,0.5  $ amesh,all     &

65、#160;                             ! 創(chuàng)建面、劃分網格等        nsel,s,loc,x,0  $ d,all,all  $ allsel    &

66、#160;                                    ! 施加約束等        f,2,fy,-10  $ finish 

67、60;                                                 

68、60;           ! 施加端部集中荷載等        /solu  $ iswrite,on  $ solve                       &

69、#160;                          ! 求解，生成初應力文件等        /post1  $ pldisp  $ plnsol,s,x        

70、                                       ! 查看求解結果等        /solu