偶然發現ANSYS命令流集錦,與大家共享:
Fini(退出四大模塊,回到BEGIN層)
/cle (清空內存,開始新的計算)
1. 定義參數、數組,并賦值.
2. /prep7(進入前處理)
定義幾何圖形:關鍵點、線、面、體
定義幾個所關心的節點,以備后處理時調用節點號。
設材料線彈性、非線性特性
設置單元類型及相應KEYOPT
設置實常數
設置網格劃分,劃分網格
根據需要耦合某些節點自由度
定義單元表
存盤
3./solu
加邊界條件
設置求解選項
定義載荷步
求解載荷步
4./post1(通用后處理)
5./post26 (時間歷程后處理)
6.PLOTCONTROL菜單命令
7.參數化設計語言
8.理論手冊
Fini(退出四大模塊,回到BEGIN層)
/cle (清空內存,開始新的計算)
1. 定義參數、數組,并賦值.
2. /prep7(進入前處理)
定義幾何圖形:關鍵點、線、面、體
定義幾個所關心的節點,以備后處理時調用節點號。
設材料線彈性、非線性特性
設置單元類型及相應KEYOPT
設置實常數
設置網格劃分,劃分網格
根據需要耦合某些節點自由度
定義單元表
存盤
3./solu
加邊界條件
設置求解選項
定義載荷步
求解載荷步
4./post1(通用后處理)
5./post26 (時間歷程后處理)
6.PLOTCONTROL菜單命令
7.參數化設計語言
8.理論手冊
2. 體映射網格劃分
(1)若將體劃分為六面體單元,必須滿足以下條件
a. 該體的外形為塊狀(六面體)、楔形或棱形(五面體)、四面體
b. 對邊必須劃分為相同的單元數,或分割符合過渡網格形式
c. 如果體是棱形或四面體,三角形面上的單元分割數必須是偶數
(2) 當需要減少圍成體的面數以進行映射網格劃分時,可以對面相加或連接。如果連接而有邊界線,線也必須連接在一起。
(3)體掃掠生成網格
步驟:
a. 確定體的拓撲是否能夠進行掃掠。側面不能有孔;體內不能有封閉腔;源面與目標面必須相對
b. 定義合適的單元類型
c. 確定掃掠操作中如何控制生成單元層的數目 lesize
d. 確定體的哪一個邊界面作為源面、目標面
e. 有選擇地對源面、目標面和邊界面劃分網格
3. 關于連接線和面的一些說明
連接僅是映射網格劃分的輔助工具
4. 用desize定義單元尺寸時單元劃分應遵守的級別
高:lesize
kesize
esize
desize
用smartzing定義單元尺寸時單元劃分應遵守的級別
高:lesize
kesize
smartsize
u LESIZE,NL1,Size, Angsiz,ndiv,space,kforc,layer1,layer2,kyndiv
為線指定網格尺寸
NL1: 線號,如果為all,則指定所有選中線的網格。
Size: 單元邊長,(程序據size計算分割份數,自動取整到下一個整數)?
Angsiz: 弧線時每單元跨過的度數?
Ndiv: 分割份數
Space: “+”: 最后尺寸比最先尺寸
“-“: 中間尺寸比兩端尺寸
free: 由其他項控制尺寸
kforc 0: 僅設置未定義的線,
1:設置所有選定線,
2:僅改設置份數少的,
3:僅改設置份數多的
kyndiv: 0,No,off 表示不可改變指定尺寸
1,yes,on 表示可改變
u ESIZE,size,ndiv 指定線的缺省劃分份數
(已直接定義的線,關鍵點網格劃分設置不受影響)
u desize, minl, minh,…… 控制缺省的單元尺寸
minl: n 每根線上低階單元數(缺省為3)
defa 缺省值
stat 列出當前設置
off 關閉缺省單元尺寸
minh: n 每根線上(高階)單元數(缺省為2)
u mshape, key, dimension 指定單元形狀
key: 0 四邊形(2D),六面體(3D)
1 三角形 (2D), 四面體(3D)
Dimension: 2D 二維
3D 三維
u smart,off 關閉智能網格
u mshkey, key 指定自由或映射網格方式
key: 0 自由網格劃分
1 映射網格劃分
2 如果可能的話使用映射,否則自由(即使自由smartsizing也不管用了)
u Amesh, nA1,nA2,ninc 劃分面單元網格
nA1,nA2,ninc 待劃分的面號,nA1如果是All,則對所有選中面劃分
u SECTYPE, ID, TYPE, SUBTYPE, NAME, REFINEKEY
定義一個截面號,并初步定義截面類型
ID: 截面號
TYPE: BEAM:定義此截面用于梁
SUBTYPE: RECT 矩形
CSOLID:圓形實心截面
CTUBE: 圓管
I: 工字形
HREC: 矩形空管
ASEC: 任意截面
MESH: 用戶定義的劃分網格
NAME: 8字符的截面名稱(字母和數字組成)
REFINEKEY: 網格細化程度:0~5(對于薄壁構件用此控制,對于實心截面用SECDATA控制)
u SECDATA, VAL1, VAL2, …….VAL10 描述梁截面
說明:對于SUBTYPE=MESH, 所需數據由SECWRITE產生,SECREAD讀入 #p#分頁標題#e#
u SECNUM,SECID 設定隨后梁單元劃分將要使用的截面編號
u LATT, MAT, REAL, TYPE, --, KB, KE, SECNUM
為準備劃分的線定義一系列特性
MAT: 材料號
REAL: 實常數號
TYPE: 線單元類型號
KB、KE: 待劃分線的定向關鍵點起始、終止號
SECNUM: 截面類型號
u SECPLOT,SECID,MESHKEY 畫梁截面的幾何形狀及網格劃分
SECID:由SECTYPE命令分配的截面編號
MESHKEY:0:不顯示網格劃分
1:顯示網格劃分
u /ESHAPE, SCALE 按看似固體化分的形式顯示線、面單元
SCALE: 0:簡單顯示線、面單元
1:使用實常數顯示單元形狀
u esurf, xnode, tlab, shape 在已存在的選中單元的自由表面覆蓋產生單元
xnode: 僅為產生surf151 或surf152單元時使用
tlab: 僅用來生成接觸元或目標元
top 產生單元且法線方向與所覆蓋的單元相同,僅對梁或殼有效,對實體單元無效
Bottom產生單元且法線方向與所覆蓋的單元相反,僅對梁或殼有效,對實體單元無效
Reverse 將已產生單元反向
Shape: 空 與所覆蓋單元形狀相同
Tri 產生三角形表面的目標元
注意:選中的單元是由所選節點決定的,而不是選單元,如同將壓力加在節點上而不是單元上
u Nummrg,label,toler, Gtoler,action,switch 合并相同位置的item
label: 要合并的項目
node: 節點, Elem,單元,kp: 關鍵點(也合并線,面及點)
mat: 材料,type: 單元類型,Real: 實常數
cp:耦合項,CE:約束項,CE: 約束方程,All:所有項
toler: 公差
Gtoler:實體公差
Action: sele 僅選擇不合并
空 合并
switch: 較低號還是較高號被保留(low, high)
注意:可以先選擇一部分項目,再執行合并。如果多次發生合并命令,一定要先合并節點,再合并關鍵點。合并節點后,實體荷載不能轉化到單元,此時可合并關鍵點解決問題。
u Lsel, type, item, comp, vmin, vmax, vinc, kswp 選擇線
type: s 從全部線中選一組線
r 從當前選中線中選一組線
a 再選一部線附加給當前選中組
au
none
u(unselect)
inve: 反向選擇
item: line 線號
loc 坐標
length 線長
comp: x,y,z
kswp: 0 只選線
1 選擇線及相關關鍵點、節點和單元
u Nsel, type, item, comp, vmin, vmax, vinc, kabs 選擇一組節點為下一步做準備
Type: S: 選擇一組新節點(缺省)
R: 在當前組中再選擇
A: 再選一組附加于當前組
U: 在當前組中不選一部分
All: 恢復為選中所有
None: 全不選
Inve: 反向選擇
Stat: 顯示當前選擇狀態
Item: loc: 坐標
node: 節點號
Comp: 分量
Vmin,vmax,vinc: ITEM范圍
Kabs: “0” 使用正負號
“1”僅用絕對值
u NSLL,type, nkey 選擇與所選線相聯系的節點
u nsla, type, nkey: 選擇與選中面相關的節點
type:s 選一套新節點
r 從已選節點中再選
a 附加一部分節點到已選節點
u 從已選節點中去除一部分
nkey: 0 僅選面內的節點
1 選所有和面相聯系的節點(如面內線,關鍵點處的節點)
u esel, type, item, comp, vmin, vmax, vinc, kabs 選擇一組單元
Type: S: 選擇一組單元(缺省)
R: 在當前組中再選一部分作為一組
A: 為當前組附加單元
U: 在當前組中不選一部分單元
All: 選所有單元
None: 全不選
Inve: 反向選擇當前組(?)
Stat: 顯示當前選擇狀態
Item: Elem: 單元號
Type: 單元類型號
Mat: 材料號
Real: 實常數號
Esys: 單元坐標系號
u ALLSEL, LABT, ENTITY 選中所有項目
LABT: ALL: 選所有項目及其低級項目
BELOW: 選指定項目的直接下屬及更低級項目
ENTITY: ALL: 所有項目(缺省)
VOLU:體 高級
AREA:面
LINE :線
KP:關鍵點
ELEM:單元
NODE:節點 低級
u Tshap,shape 定義接觸目標面為2D、3D的簡單圖形
Shape: line:直線
Arc:順時針弧
Tria:3點三角形
Quad:4點四邊形
…………. #p#分頁標題#e#
2.6 根據需要耦合某些節點自由度
u cp, nset, lab,,node1,node2,……node17
nset: 耦合組編號
lab: ux,uy,uz,rotx,roty,rotz
node1-node17: 待耦合的節點號。如果某一節點號為負,則此節點從該耦合組中刪去。如果node1=all,則所有選中節點加入該耦合組。
注意:1,不同自由度類型將生成不同編號
2,不可將同一自由度用于多套耦合組
u CPINTF, LAB, TOLER 將相鄰節點的指定自由度定義為耦合自由度
LAB:UX,UY,UZ,ROTX,ROTY,ROTZ,ALL
TOLER: 公差,缺省為0.0001
說明:先選中欲耦合節點,再執行此命令
2.7 定義單元表
說明:1,單元表僅對選中單元起作用,使用單元表之前務必選擇一種類型的單元
2,單元表各行為選中各單元,各列為每單元的不同數據
u ETABLE, LAB, ITEM, COMP 定義單元表,添加、刪除單元表某列
LAB:用戶指定的列名(REFL, STAT, ERAS 為預定名稱)
ITEM: 數據標志(查各單元可輸出項目)
COMP: 數據分量標志
2.8 存盤
u save, fname, ext,dir, slab 存盤
fname : 文件名(最多32個字符)缺省為工作名
ext: 擴展名(最多32個字符)缺省為db
dir: 目錄名(最多64個字符)缺省為當前
slab: “all” 存所有信息
“model” 存模型信息
“solv” 存模型信息和求解信息
3 /solu
u /solu 進入求解器
3.1 加邊界條件
u D, node, lab, value, value2, nend, ninc, lab2, lab3, ……lab6 定義節點位移約束
Node : 預加位移約束的節點號,如果為all,則所有選中節點全加約束,此時忽略nend和ninc.
Lab: ux,uy,uz,rotx,roty,rotz,all
Value,value2: 自由度的數值(缺省為0)
Nend, ninc: 節點范圍為:node-nend,編號間隔為ninc
Lab2-lab6: 將lab2-lab6以同樣數值施加給所選節點。
注意:在節點坐標系中討論
(1)若將體劃分為六面體單元,必須滿足以下條件
a. 該體的外形為塊狀(六面體)、楔形或棱形(五面體)、四面體
b. 對邊必須劃分為相同的單元數,或分割符合過渡網格形式
c. 如果體是棱形或四面體,三角形面上的單元分割數必須是偶數
(2) 當需要減少圍成體的面數以進行映射網格劃分時,可以對面相加或連接。如果連接而有邊界線,線也必須連接在一起。
(3)體掃掠生成網格
步驟:
a. 確定體的拓撲是否能夠進行掃掠。側面不能有孔;體內不能有封閉腔;源面與目標面必須相對
b. 定義合適的單元類型
c. 確定掃掠操作中如何控制生成單元層的數目 lesize
d. 確定體的哪一個邊界面作為源面、目標面
e. 有選擇地對源面、目標面和邊界面劃分網格
3. 關于連接線和面的一些說明
連接僅是映射網格劃分的輔助工具
4. 用desize定義單元尺寸時單元劃分應遵守的級別
高:lesize
kesize
esize
desize
用smartzing定義單元尺寸時單元劃分應遵守的級別
高:lesize
kesize
smartsize
u LESIZE,NL1,Size, Angsiz,ndiv,space,kforc,layer1,layer2,kyndiv
為線指定網格尺寸
NL1: 線號,如果為all,則指定所有選中線的網格。
Size: 單元邊長,(程序據size計算分割份數,自動取整到下一個整數)?
Angsiz: 弧線時每單元跨過的度數?
Ndiv: 分割份數
Space: “+”: 最后尺寸比最先尺寸
“-“: 中間尺寸比兩端尺寸
free: 由其他項控制尺寸
kforc 0: 僅設置未定義的線,
1:設置所有選定線,
2:僅改設置份數少的,
3:僅改設置份數多的
kyndiv: 0,No,off 表示不可改變指定尺寸
1,yes,on 表示可改變
u ESIZE,size,ndiv 指定線的缺省劃分份數
(已直接定義的線,關鍵點網格劃分設置不受影響)
u desize, minl, minh,…… 控制缺省的單元尺寸
minl: n 每根線上低階單元數(缺省為3)
defa 缺省值
stat 列出當前設置
off 關閉缺省單元尺寸
minh: n 每根線上(高階)單元數(缺省為2)
u mshape, key, dimension 指定單元形狀
key: 0 四邊形(2D),六面體(3D)
1 三角形 (2D), 四面體(3D)
Dimension: 2D 二維
3D 三維
u smart,off 關閉智能網格
u mshkey, key 指定自由或映射網格方式
key: 0 自由網格劃分
1 映射網格劃分
2 如果可能的話使用映射,否則自由(即使自由smartsizing也不管用了)
u Amesh, nA1,nA2,ninc 劃分面單元網格
nA1,nA2,ninc 待劃分的面號,nA1如果是All,則對所有選中面劃分
u SECTYPE, ID, TYPE, SUBTYPE, NAME, REFINEKEY
定義一個截面號,并初步定義截面類型
ID: 截面號
TYPE: BEAM:定義此截面用于梁
SUBTYPE: RECT 矩形
CSOLID:圓形實心截面
CTUBE: 圓管
I: 工字形
HREC: 矩形空管
ASEC: 任意截面
MESH: 用戶定義的劃分網格
NAME: 8字符的截面名稱(字母和數字組成)
REFINEKEY: 網格細化程度:0~5(對于薄壁構件用此控制,對于實心截面用SECDATA控制)
u SECDATA, VAL1, VAL2, …….VAL10 描述梁截面
說明:對于SUBTYPE=MESH, 所需數據由SECWRITE產生,SECREAD讀入 #p#分頁標題#e#
u SECNUM,SECID 設定隨后梁單元劃分將要使用的截面編號
u LATT, MAT, REAL, TYPE, --, KB, KE, SECNUM
為準備劃分的線定義一系列特性
MAT: 材料號
REAL: 實常數號
TYPE: 線單元類型號
KB、KE: 待劃分線的定向關鍵點起始、終止號
SECNUM: 截面類型號
u SECPLOT,SECID,MESHKEY 畫梁截面的幾何形狀及網格劃分
SECID:由SECTYPE命令分配的截面編號
MESHKEY:0:不顯示網格劃分
1:顯示網格劃分
u /ESHAPE, SCALE 按看似固體化分的形式顯示線、面單元
SCALE: 0:簡單顯示線、面單元
1:使用實常數顯示單元形狀
u esurf, xnode, tlab, shape 在已存在的選中單元的自由表面覆蓋產生單元
xnode: 僅為產生surf151 或surf152單元時使用
tlab: 僅用來生成接觸元或目標元
top 產生單元且法線方向與所覆蓋的單元相同,僅對梁或殼有效,對實體單元無效
Bottom產生單元且法線方向與所覆蓋的單元相反,僅對梁或殼有效,對實體單元無效
Reverse 將已產生單元反向
Shape: 空 與所覆蓋單元形狀相同
Tri 產生三角形表面的目標元
注意:選中的單元是由所選節點決定的,而不是選單元,如同將壓力加在節點上而不是單元上
u Nummrg,label,toler, Gtoler,action,switch 合并相同位置的item
label: 要合并的項目
node: 節點, Elem,單元,kp: 關鍵點(也合并線,面及點)
mat: 材料,type: 單元類型,Real: 實常數
cp:耦合項,CE:約束項,CE: 約束方程,All:所有項
toler: 公差
Gtoler:實體公差
Action: sele 僅選擇不合并
空 合并
switch: 較低號還是較高號被保留(low, high)
注意:可以先選擇一部分項目,再執行合并。如果多次發生合并命令,一定要先合并節點,再合并關鍵點。合并節點后,實體荷載不能轉化到單元,此時可合并關鍵點解決問題。
u Lsel, type, item, comp, vmin, vmax, vinc, kswp 選擇線
type: s 從全部線中選一組線
r 從當前選中線中選一組線
a 再選一部線附加給當前選中組
au
none
u(unselect)
inve: 反向選擇
item: line 線號
loc 坐標
length 線長
comp: x,y,z
kswp: 0 只選線
1 選擇線及相關關鍵點、節點和單元
u Nsel, type, item, comp, vmin, vmax, vinc, kabs 選擇一組節點為下一步做準備
Type: S: 選擇一組新節點(缺省)
R: 在當前組中再選擇
A: 再選一組附加于當前組
U: 在當前組中不選一部分
All: 恢復為選中所有
None: 全不選
Inve: 反向選擇
Stat: 顯示當前選擇狀態
Item: loc: 坐標
node: 節點號
Comp: 分量
Vmin,vmax,vinc: ITEM范圍
Kabs: “0” 使用正負號
“1”僅用絕對值
u NSLL,type, nkey 選擇與所選線相聯系的節點
u nsla, type, nkey: 選擇與選中面相關的節點
type:s 選一套新節點
r 從已選節點中再選
a 附加一部分節點到已選節點
u 從已選節點中去除一部分
nkey: 0 僅選面內的節點
1 選所有和面相聯系的節點(如面內線,關鍵點處的節點)
u esel, type, item, comp, vmin, vmax, vinc, kabs 選擇一組單元
Type: S: 選擇一組單元(缺省)
R: 在當前組中再選一部分作為一組
A: 為當前組附加單元
U: 在當前組中不選一部分單元
All: 選所有單元
None: 全不選
Inve: 反向選擇當前組(?)
Stat: 顯示當前選擇狀態
Item: Elem: 單元號
Type: 單元類型號
Mat: 材料號
Real: 實常數號
Esys: 單元坐標系號
u ALLSEL, LABT, ENTITY 選中所有項目
LABT: ALL: 選所有項目及其低級項目
BELOW: 選指定項目的直接下屬及更低級項目
ENTITY: ALL: 所有項目(缺省)
VOLU:體 高級
AREA:面
LINE :線
KP:關鍵點
ELEM:單元
NODE:節點 低級
u Tshap,shape 定義接觸目標面為2D、3D的簡單圖形
Shape: line:直線
Arc:順時針弧
Tria:3點三角形
Quad:4點四邊形
…………. #p#分頁標題#e#
2.6 根據需要耦合某些節點自由度
u cp, nset, lab,,node1,node2,……node17
nset: 耦合組編號
lab: ux,uy,uz,rotx,roty,rotz
node1-node17: 待耦合的節點號。如果某一節點號為負,則此節點從該耦合組中刪去。如果node1=all,則所有選中節點加入該耦合組。
注意:1,不同自由度類型將生成不同編號
2,不可將同一自由度用于多套耦合組
u CPINTF, LAB, TOLER 將相鄰節點的指定自由度定義為耦合自由度
LAB:UX,UY,UZ,ROTX,ROTY,ROTZ,ALL
TOLER: 公差,缺省為0.0001
說明:先選中欲耦合節點,再執行此命令
2.7 定義單元表
說明:1,單元表僅對選中單元起作用,使用單元表之前務必選擇一種類型的單元
2,單元表各行為選中各單元,各列為每單元的不同數據
u ETABLE, LAB, ITEM, COMP 定義單元表,添加、刪除單元表某列
LAB:用戶指定的列名(REFL, STAT, ERAS 為預定名稱)
ITEM: 數據標志(查各單元可輸出項目)
COMP: 數據分量標志
2.8 存盤
u save, fname, ext,dir, slab 存盤
fname : 文件名(最多32個字符)缺省為工作名
ext: 擴展名(最多32個字符)缺省為db
dir: 目錄名(最多64個字符)缺省為當前
slab: “all” 存所有信息
“model” 存模型信息
“solv” 存模型信息和求解信息
3 /solu
u /solu 進入求解器
3.1 加邊界條件
u D, node, lab, value, value2, nend, ninc, lab2, lab3, ……lab6 定義節點位移約束
Node : 預加位移約束的節點號,如果為all,則所有選中節點全加約束,此時忽略nend和ninc.
Lab: ux,uy,uz,rotx,roty,rotz,all
Value,value2: 自由度的數值(缺省為0)
Nend, ninc: 節點范圍為:node-nend,編號間隔為ninc
Lab2-lab6: 將lab2-lab6以同樣數值施加給所選節點。
注意:在節點坐標系中討論
3.2 設置求解選項
u antype, status, ldstep, substep, action
antype: static or 1 靜力分析
buckle or 2 屈曲分析
modal or 3 模態分析
trans or 4 瞬態分析
status: new 重新分析(缺省),以后各項將忽略
rest 再分析,僅對static,full transion 有效
ldstep: 指定從哪個荷載步開始繼續分析,缺省為最大的,runn數(指分析點的最后一步)
substep: 指定從哪個子步開始繼續分析。缺省為本目錄中,runn文件中最高的子步數
action, continue: 繼續分析指定的ldstep,substep
說明:繼續以前的分析(因某種原因中斷)有兩種類型
singleframe restart: 從停止點繼續
需要文件:jobname.db 必須在初始求解后馬上存盤
jobname.emat 單元矩陣
jobname.esav 或 .osav : 如果.esav壞了,將.osav改為.esav
results file: 不必要,但如果有,后繼分析的結果也將很好地附加到它后面
注意:如果初始分析生成了.rdb, .ldhi, 或rnnn 文件。必須刪除再做后繼分析
步驟: (1)進入anasys 以同樣工作名
(2)進入求解器,并恢復數據庫
(3)antype, rest
(4)指定附加的荷載
(5)指定是否使用現有的矩陣(jobname.trl)(缺省重新生成)
kuse: 1 用現有矩陣
(6)求解
multiframe restart:從以有結果的任一步繼續(用不著)
u pred,sskey, --,lskey….. 在非線性分析中是否打開預測器
sskey: off 不作預測(當有旋轉自由度時或使用solid65時缺省為off)
on 第一個子步后作預測(除非有旋轉自由度時或使用solid65時缺省為on)
-- : 未使用變量區
lskey: off 跨越荷載步時不作預測(缺省)
on 跨越荷載步時作預測(此時sskey必須同時on)
注意:此命令的缺省值假定solcontrol為on
u autots, key 是否使用自動時間步長
key
n: 當solcontrol為on時缺省為on
off: 當solcontrol為off時缺省為off
1: 由程序選擇(當solcontrol為on且不發生autots命令時在 .log文件中紀錄“1”
注意:當使用自動時間步長時,也會使用步長預測器和二分步長 #p#分頁標題#e#
u NROPT, option,--,adptky 指定牛頓拉夫遜法求解的選項
OPTION: AUT程序選擇
FULL:完全牛頓拉夫遜法
MODI:修正的牛頓拉夫遜法
INIT:使用初始剛陣
UNSYM:完全牛頓拉夫遜法,且允許非對稱剛陣
ADPTKY:ON: 使用自適應下降因子
OFF:不使用自適應下降因子
u NLGEOM,KEY
KEY: OFF:不包括幾何非線性(缺省)
ON:包括幾何非線性
u ncnv, kstop, dlim, itlim, etlim, cplim 終止分析選項
kstop: 0 如果求解不收斂,也不終止分析
1 如果求解不收斂,終止分析和程序(缺省)
2如果求解不收斂,終止分析,但不終止程序
dlim:最大位移限制,缺省為1.0e6
itlim: 累積迭代次數限制,缺省為無窮多
etlim:程序執行時間(秒)限制,缺省為無窮
cplim:cpu時間(秒)限制,缺省為無窮
u solcontrol ,key1, key2,key3,vtol 指定是否使用一些非線性求解缺省值
key1: on 激活一些優化缺省值(缺省)
CNVTOL Toler=0.5%Minref=0.01(對力和彎矩)
NEQIT 最大迭代次數根據模型設定在15~26之間
ARCLEN 如用弧長法則用較ansys5.3更先進的方法
PRED 除非有rotx,y,z或solid65,否則打開
LNSRCH 當有接觸時自動打開
CUTCONTROL Plslimit=15%, npoint=13
SSTIF 當NLGEOM,on時則打開
NROPT,adaptkey 關閉(除非:摩擦接觸存在;單元12,26,48,49,52存在;當塑性存在且有單元20,23,24,60存在)
AUTOS 由程序選擇
off 不使用這些缺省值
key2: on 檢查接觸狀態(此時key1為on)
此時時間步會以單元的接觸狀態(據keyopt(7)的假定)為基礎
當keyopt(2)=on 時,保證時間步足夠小
key3: 應力荷載剛化控制,盡量使用缺省值
空:缺省,對某些單元包括應力荷載剛化,對某些不包括(查)
nopl:對任何單元不包括應力剛化
incp:對某些單元包括應力荷載剛化(查)
vtol:
u outres, item, freq, cname 規定寫入數據庫的求解信息
item: all 所有求解項
basic 只寫nsol, rsol, nload, strs
nsol 節點自由度
rsol 節點作用荷載
nload 節點荷載和輸入的應變荷載(?)
strs 節點應力
freq: 如果為n,則每n步(包括最后一步)寫入一次
none: 則在此荷載步中不寫次項
all: 每一步都寫
last: 只寫最后一步(靜力或瞬態時為缺省)
3.3 定義載荷步
u nsubst, nsbstp, nsbmx, nsbmn, carry 指定此荷載步的子步數
nsbstp: 此荷載步的子步數
如果自動時間步長使用autots,則此數定義第一子步的長度;如果solcontrol打開,且3D面-面接觸單元使用,則缺省為1-20步;如果solcontrol打開,并無3D接觸單元,則缺省為1子步;如果solcontrol關閉,則缺省為以前指定值;如以前未指定,則缺省為1)
nsbmx, nsbmn:最多,最少子步數(如果自動時間步長打開)?
u time, time 指定荷載步結束時間
注意:第一步結束時間不可為“0”
u f, node, lab, value, value2, nend, ninc 在指定節點加集中荷載
node:節點號
lab: Fx,Fy,Fz,Mx,My,Mz
value: 力大小
value2: 力的第二個大小(如果有復數荷載)
nend,ninc:在從node到nend的節點(增量為ninc)上施加同樣的力
注意:(1)節點力在節點坐標系中定義,其正負與節點坐標軸正向一致
u sfa, area, lkey, lab, value, value2 在指定面上加荷載
area: n 面號
all 所有選中號
lkey: 如果是體的面,忽略此項
lab: pres
value: 壓力值
u SFBEAM, ELEM, LKEY, LAB, VALI, VALJ, VAL2I, VAL2J, IOFFST, JOFFST
對梁單元施加線荷載
ELEM: 單元號,可以為ALL,即選中單元
LKEY: 面載類型號,見單元介紹。對于BEAM188,1為豎向;2為橫向;3為切向
VALI,VALJ: I, J節點處壓力值
VAL2I,VAL2J: 暫時無用
IOFFST, JOFFST: 線載距離I, J 節點距離
u lswrite, lsnum 將荷載與荷載選項寫入荷載文件中
lsnum :荷載步文件名的后綴,即荷載步數
當 stat 列示當前步數
init 重設為“1”
缺省為當前步數加“1”
3.3.1 注意
1. 盡量加面載,不加集中力,以免奇異點 #p#分頁標題#e#
2. 面的切向荷載必須借助面單元
3.4 求解載荷步
u lssolve, lsmin, lsmax, lsinc 讀入并求解多個荷載步
lsmin, lsmax, lsinc :荷載步文件范圍
4 /post1(通用后處理)
u set, lstep, sbstep, fact, king, time, angle, nset 設定從結果文件讀入的數據
lstep :荷載步數
sbstep:子步數,缺省為最后一步
time: 時間點(如果弧長法則不用)
nset: data set number
u dscale, wn, dmult 顯示變形比例
wn: 窗口號(或all),缺省為1
dmult, 0或auto : 自動將最大變形圖畫為構件長的5%
u pldisp, kund 顯示變形的結構
kund: 0 僅顯示變形后的結構
1 顯示變形前和變形后的結構
2 顯示變形結構和未變形結構的邊緣
u *get, par, node, n, u, x(y,z) 獲得節點n的x(y,z)位移給參數par
等價于函數 ux(n),uy(n),uz(z)
node(x,y,z): 獲得(x,y,z)節點號
arnode(x,y,z):獲得和節點n相連的面
注意:此命令也可用于/solu模塊
u fsum, lab, item 對單元之節點力和力矩求和
lab: 空 在整體迪卡爾坐標系下求和
rsys 在當前激活的rsys坐標系下求和
item: 空 對所有選中單元(不包括接觸元)求和
cont: 僅對接觸節點求和
u PRSSOL, ITEM, COMP 打印BEAM188、BEAM189截面結果
說明:只有剛計算完還未退出ANSYS時可用,重新進入ANSYS時不可用
item comp 截面數據及分量標志
S COMP X,XZ,YZ應力分量
PRIN S1,S2,S3主應力SINT應力強度,SEQV等效應力
EPTO COMP 總應變
PRIN 總主應變,應變強度,等效應變
EPPL COMP 塑性應變分量
PRIN 主塑性應變,塑性應變強度,等效塑性應變
u plnsol, item, comp, kund, fact 畫節點結果為連續的輪廓線
item: 項目(見下表)
comp: 分量
kund: 0 不顯示未變形的結構
1 變形和未變形重疊
2 變形輪廓和未變形邊緣
fact: 對于接觸的2D顯示的比例系數,缺省為1
item comp discription
u x,y,z,sum 位移
rot x,y,z,sum 轉角
s x,y,z,xy,yz,xz 應力分量
1,2,3 主應力
Int,eqv 應力intensity,等效應力
epeo x,y,z,xy,yz,xz 總位移分量
1,2,3 主應變
Int,eqv 應變intensity,等效應變
epel x,y,z,xy,yz,xz 彈性應變分量
1,2,3 彈性主應變
Int,eqv 彈性intensity,彈性等效應變
eppl x,y,z,xy,yz,xz 塑性應變分量
u PRNSOL, item, comp 打印選中節點結果
item: 項目(見上表)
comp: 分量
u PRETAB, LAB1, LAB2, ……LAB9 沿線單元長度方向繪單元表數據
LABn : 空: 所有ETABLE命令指定的列名
列名: 任何ETABLE命令指定的列名
u PLLS, LABI, LABJ, FACT, KUND 沿線單元長度方向繪單元表數據
LABI:節點I的單元表列名
LABJ:節點J的單元表列名
FACT: 顯示比例,缺省為1
kund: 0 不顯示未變形的結構
1 變形和未變形重疊
2 變形輪廓和未變形邊緣
5 /post26 (時間歷程后處理)
u nsol, nvar, node, item, comp,name
在時間歷程后處理器中定義節點變量的序號
nvar:變量號(從2到nv(根據numvar定義))
node: 節點號
item comp
u x, y,z
rot x, y,z
u ESOL, NVAR, ELEM, NODE, ITEM, COMP, NAME 將結果存入變量
NVAR: 變量號,2以上
ELEM: 單元號
NODE: 該單元的節點號,決定存儲該單元的哪個量,如果空,則給出平均值
ITEM:
COMP:
NAME: 8字符的變量名, 缺省為ITEM加COMP
u rforce, nvar, node, item, comp, name 指定待存儲的節點力數據
nvar: 變量號
node: 節點號
item comp
F x, y.z
M x, y,z
name: 給此變量一個名稱,8個字符
u add, ir, ia,ib,ic,name,--,--,facta, factb, factc
將ia,ib,ic變量相加賦給ir變量
ir, ia,ib,ic:變量號
name: 變量的名稱
u /grid, key
key: “0” 或“off” 無網絡
“1”或“on” xy網絡
“2”或“x” 只有x線 #p#分頁標題#e#
“3”或“y” 只有y線
u xvar, n
n: “0”或“1” 將x軸作為時間軸
“n” 將x軸表示變量“n”
“-1” ?
u /axlab, axis, lab 定義軸線的標志
axis: “x”或“y”
lab: 標志,可長達30個字符
u plvar, nvar, nvar2, ……,nvar10 畫出要顯示的變量(作為縱坐標)
u prvar, nvar1, ……,nvar6 列出要顯示的變量
6 PLOTCONTROL菜單命令
u pbc, ilem, ……,key, min, max, abs 在顯示屏上顯示符號及數值
item: u 所加的位移約束
rot 所加的轉角約束
key: 0 不顯示符號
1 顯示符號
2 顯示符號及數值
u /SHOW, FNAME, EXT, VECT, NCPL 確定圖形顯示的設備及其他參數
FNAME: X11:屏幕
文件名:各圖形將生成一系列圖形文件
JPEG: 各圖形將生成一系列JPEG圖形文件
說明:沒必要用此命令,需要的圖形文件可計算后再輸出
u antype, status, ldstep, substep, action
antype: static or 1 靜力分析
buckle or 2 屈曲分析
modal or 3 模態分析
trans or 4 瞬態分析
status: new 重新分析(缺省),以后各項將忽略
rest 再分析,僅對static,full transion 有效
ldstep: 指定從哪個荷載步開始繼續分析,缺省為最大的,runn數(指分析點的最后一步)
substep: 指定從哪個子步開始繼續分析。缺省為本目錄中,runn文件中最高的子步數
action, continue: 繼續分析指定的ldstep,substep
說明:繼續以前的分析(因某種原因中斷)有兩種類型
singleframe restart: 從停止點繼續
需要文件:jobname.db 必須在初始求解后馬上存盤
jobname.emat 單元矩陣
jobname.esav 或 .osav : 如果.esav壞了,將.osav改為.esav
results file: 不必要,但如果有,后繼分析的結果也將很好地附加到它后面
注意:如果初始分析生成了.rdb, .ldhi, 或rnnn 文件。必須刪除再做后繼分析
步驟: (1)進入anasys 以同樣工作名
(2)進入求解器,并恢復數據庫
(3)antype, rest
(4)指定附加的荷載
(5)指定是否使用現有的矩陣(jobname.trl)(缺省重新生成)
kuse: 1 用現有矩陣
(6)求解
multiframe restart:從以有結果的任一步繼續(用不著)
u pred,sskey, --,lskey….. 在非線性分析中是否打開預測器
sskey: off 不作預測(當有旋轉自由度時或使用solid65時缺省為off)
on 第一個子步后作預測(除非有旋轉自由度時或使用solid65時缺省為on)
-- : 未使用變量區
lskey: off 跨越荷載步時不作預測(缺省)
on 跨越荷載步時作預測(此時sskey必須同時on)
注意:此命令的缺省值假定solcontrol為on
u autots, key 是否使用自動時間步長
key

off: 當solcontrol為off時缺省為off
1: 由程序選擇(當solcontrol為on且不發生autots命令時在 .log文件中紀錄“1”
注意:當使用自動時間步長時,也會使用步長預測器和二分步長 #p#分頁標題#e#
u NROPT, option,--,adptky 指定牛頓拉夫遜法求解的選項
OPTION: AUT程序選擇
FULL:完全牛頓拉夫遜法
MODI:修正的牛頓拉夫遜法
INIT:使用初始剛陣
UNSYM:完全牛頓拉夫遜法,且允許非對稱剛陣
ADPTKY:ON: 使用自適應下降因子
OFF:不使用自適應下降因子
u NLGEOM,KEY
KEY: OFF:不包括幾何非線性(缺省)
ON:包括幾何非線性
u ncnv, kstop, dlim, itlim, etlim, cplim 終止分析選項
kstop: 0 如果求解不收斂,也不終止分析
1 如果求解不收斂,終止分析和程序(缺省)
2如果求解不收斂,終止分析,但不終止程序
dlim:最大位移限制,缺省為1.0e6
itlim: 累積迭代次數限制,缺省為無窮多
etlim:程序執行時間(秒)限制,缺省為無窮
cplim:cpu時間(秒)限制,缺省為無窮
u solcontrol ,key1, key2,key3,vtol 指定是否使用一些非線性求解缺省值
key1: on 激活一些優化缺省值(缺省)
CNVTOL Toler=0.5%Minref=0.01(對力和彎矩)
NEQIT 最大迭代次數根據模型設定在15~26之間
ARCLEN 如用弧長法則用較ansys5.3更先進的方法
PRED 除非有rotx,y,z或solid65,否則打開
LNSRCH 當有接觸時自動打開
CUTCONTROL Plslimit=15%, npoint=13
SSTIF 當NLGEOM,on時則打開
NROPT,adaptkey 關閉(除非:摩擦接觸存在;單元12,26,48,49,52存在;當塑性存在且有單元20,23,24,60存在)
AUTOS 由程序選擇
off 不使用這些缺省值
key2: on 檢查接觸狀態(此時key1為on)
此時時間步會以單元的接觸狀態(據keyopt(7)的假定)為基礎
當keyopt(2)=on 時,保證時間步足夠小
key3: 應力荷載剛化控制,盡量使用缺省值
空:缺省,對某些單元包括應力荷載剛化,對某些不包括(查)
nopl:對任何單元不包括應力剛化
incp:對某些單元包括應力荷載剛化(查)
vtol:
u outres, item, freq, cname 規定寫入數據庫的求解信息
item: all 所有求解項
basic 只寫nsol, rsol, nload, strs
nsol 節點自由度
rsol 節點作用荷載
nload 節點荷載和輸入的應變荷載(?)
strs 節點應力
freq: 如果為n,則每n步(包括最后一步)寫入一次
none: 則在此荷載步中不寫次項
all: 每一步都寫
last: 只寫最后一步(靜力或瞬態時為缺省)
3.3 定義載荷步
u nsubst, nsbstp, nsbmx, nsbmn, carry 指定此荷載步的子步數
nsbstp: 此荷載步的子步數
如果自動時間步長使用autots,則此數定義第一子步的長度;如果solcontrol打開,且3D面-面接觸單元使用,則缺省為1-20步;如果solcontrol打開,并無3D接觸單元,則缺省為1子步;如果solcontrol關閉,則缺省為以前指定值;如以前未指定,則缺省為1)
nsbmx, nsbmn:最多,最少子步數(如果自動時間步長打開)?
u time, time 指定荷載步結束時間
注意:第一步結束時間不可為“0”
u f, node, lab, value, value2, nend, ninc 在指定節點加集中荷載
node:節點號
lab: Fx,Fy,Fz,Mx,My,Mz
value: 力大小
value2: 力的第二個大小(如果有復數荷載)
nend,ninc:在從node到nend的節點(增量為ninc)上施加同樣的力
注意:(1)節點力在節點坐標系中定義,其正負與節點坐標軸正向一致
u sfa, area, lkey, lab, value, value2 在指定面上加荷載
area: n 面號
all 所有選中號
lkey: 如果是體的面,忽略此項
lab: pres
value: 壓力值
u SFBEAM, ELEM, LKEY, LAB, VALI, VALJ, VAL2I, VAL2J, IOFFST, JOFFST
對梁單元施加線荷載
ELEM: 單元號,可以為ALL,即選中單元
LKEY: 面載類型號,見單元介紹。對于BEAM188,1為豎向;2為橫向;3為切向
VALI,VALJ: I, J節點處壓力值
VAL2I,VAL2J: 暫時無用
IOFFST, JOFFST: 線載距離I, J 節點距離
u lswrite, lsnum 將荷載與荷載選項寫入荷載文件中
lsnum :荷載步文件名的后綴,即荷載步數
當 stat 列示當前步數
init 重設為“1”
缺省為當前步數加“1”
3.3.1 注意
1. 盡量加面載,不加集中力,以免奇異點 #p#分頁標題#e#
2. 面的切向荷載必須借助面單元
3.4 求解載荷步
u lssolve, lsmin, lsmax, lsinc 讀入并求解多個荷載步
lsmin, lsmax, lsinc :荷載步文件范圍
4 /post1(通用后處理)
u set, lstep, sbstep, fact, king, time, angle, nset 設定從結果文件讀入的數據
lstep :荷載步數
sbstep:子步數,缺省為最后一步
time: 時間點(如果弧長法則不用)
nset: data set number
u dscale, wn, dmult 顯示變形比例
wn: 窗口號(或all),缺省為1
dmult, 0或auto : 自動將最大變形圖畫為構件長的5%
u pldisp, kund 顯示變形的結構
kund: 0 僅顯示變形后的結構
1 顯示變形前和變形后的結構
2 顯示變形結構和未變形結構的邊緣
u *get, par, node, n, u, x(y,z) 獲得節點n的x(y,z)位移給參數par
等價于函數 ux(n),uy(n),uz(z)
node(x,y,z): 獲得(x,y,z)節點號
arnode(x,y,z):獲得和節點n相連的面
注意:此命令也可用于/solu模塊
u fsum, lab, item 對單元之節點力和力矩求和
lab: 空 在整體迪卡爾坐標系下求和
rsys 在當前激活的rsys坐標系下求和
item: 空 對所有選中單元(不包括接觸元)求和
cont: 僅對接觸節點求和
u PRSSOL, ITEM, COMP 打印BEAM188、BEAM189截面結果
說明:只有剛計算完還未退出ANSYS時可用,重新進入ANSYS時不可用
item comp 截面數據及分量標志
S COMP X,XZ,YZ應力分量
PRIN S1,S2,S3主應力SINT應力強度,SEQV等效應力
EPTO COMP 總應變
PRIN 總主應變,應變強度,等效應變
EPPL COMP 塑性應變分量
PRIN 主塑性應變,塑性應變強度,等效塑性應變
u plnsol, item, comp, kund, fact 畫節點結果為連續的輪廓線
item: 項目(見下表)
comp: 分量
kund: 0 不顯示未變形的結構
1 變形和未變形重疊
2 變形輪廓和未變形邊緣
fact: 對于接觸的2D顯示的比例系數,缺省為1
item comp discription
u x,y,z,sum 位移
rot x,y,z,sum 轉角
s x,y,z,xy,yz,xz 應力分量
1,2,3 主應力
Int,eqv 應力intensity,等效應力
epeo x,y,z,xy,yz,xz 總位移分量
1,2,3 主應變
Int,eqv 應變intensity,等效應變
epel x,y,z,xy,yz,xz 彈性應變分量
1,2,3 彈性主應變
Int,eqv 彈性intensity,彈性等效應變
eppl x,y,z,xy,yz,xz 塑性應變分量
u PRNSOL, item, comp 打印選中節點結果
item: 項目(見上表)
comp: 分量
u PRETAB, LAB1, LAB2, ……LAB9 沿線單元長度方向繪單元表數據
LABn : 空: 所有ETABLE命令指定的列名
列名: 任何ETABLE命令指定的列名
u PLLS, LABI, LABJ, FACT, KUND 沿線單元長度方向繪單元表數據
LABI:節點I的單元表列名
LABJ:節點J的單元表列名
FACT: 顯示比例,缺省為1
kund: 0 不顯示未變形的結構
1 變形和未變形重疊
2 變形輪廓和未變形邊緣
5 /post26 (時間歷程后處理)
u nsol, nvar, node, item, comp,name
在時間歷程后處理器中定義節點變量的序號
nvar:變量號(從2到nv(根據numvar定義))
node: 節點號
item comp
u x, y,z
rot x, y,z
u ESOL, NVAR, ELEM, NODE, ITEM, COMP, NAME 將結果存入變量
NVAR: 變量號,2以上
ELEM: 單元號
NODE: 該單元的節點號,決定存儲該單元的哪個量,如果空,則給出平均值
ITEM:
COMP:
NAME: 8字符的變量名, 缺省為ITEM加COMP
u rforce, nvar, node, item, comp, name 指定待存儲的節點力數據
nvar: 變量號
node: 節點號
item comp
F x, y.z
M x, y,z
name: 給此變量一個名稱,8個字符
u add, ir, ia,ib,ic,name,--,--,facta, factb, factc
將ia,ib,ic變量相加賦給ir變量
ir, ia,ib,ic:變量號
name: 變量的名稱
u /grid, key
key: “0” 或“off” 無網絡
“1”或“on” xy網絡
“2”或“x” 只有x線 #p#分頁標題#e#
“3”或“y” 只有y線
u xvar, n
n: “0”或“1” 將x軸作為時間軸
“n” 將x軸表示變量“n”
“-1” ?
u /axlab, axis, lab 定義軸線的標志
axis: “x”或“y”
lab: 標志,可長達30個字符
u plvar, nvar, nvar2, ……,nvar10 畫出要顯示的變量(作為縱坐標)
u prvar, nvar1, ……,nvar6 列出要顯示的變量
6 PLOTCONTROL菜單命令
u pbc, ilem, ……,key, min, max, abs 在顯示屏上顯示符號及數值
item: u 所加的位移約束
rot 所加的轉角約束
key: 0 不顯示符號
1 顯示符號
2 顯示符號及數值
u /SHOW, FNAME, EXT, VECT, NCPL 確定圖形顯示的設備及其他參數
FNAME: X11:屏幕
文件名:各圖形將生成一系列圖形文件
JPEG: 各圖形將生成一系列JPEG圖形文件
說明:沒必要用此命令,需要的圖形文件可計算后再輸出
7 參數化設計語言
u *do, par, ival, fval, inc 定義一個do循環的開始
par: 循環控制變量
ival, fval, inc:起始值,終值,步長(正,負)
u *enddo 定義一個do循環的結束
u *if,val1, oper, val2, base: 條件語句
val1, val2: 待比較的值(也可是字符,用引號括起來)
oper: 邏輯操作(當實數比較時,誤差為1e-10)
eq, ne, lt, gt, le, ge, ablt, abgt
base: 當oper結果為邏輯真時的行為
lable: 用戶定義的行標志
stop: 將跳出anasys
exit: 跳出當前的do循環
cycle: 跳至當前do循環的末尾
then: 構成if-then-else結構
注意:不允許跳出、跳進一個do,if循環至label句 ?
8 理論手冊
1.方程組解法:(1)直接解法;(2)迭代解法
(1) 直接解法:a.稀疏矩陣法;b. 波前解法
a. 稀疏矩陣法:占內存大,但運算次數少;通過變換剛度矩陣的順序使得非零元素最少
b. 波前解法: 占內存小
波前是指在還沒有一個單元被解完的時候激活的方程數?
(2) 迭代解法:JCG法;PCG法;ICCG法
JCG法:可解實數、對稱、非對稱矩陣
PCG法:高效求解各種矩陣(包括病態),但僅解實、對稱矩陣
ICCG法:類似JCG,但更強
2. 應變密度,等效應變,應力密度,等效應力
(1)應變密度(strain intensity)
應變密度
是三個主應變
(2)等效應變
有效泊松比 :用戶由avprin 命令設定;0(如果不設定)
(3)應力密度(stress intensity)
應力密度
(4) 等效應力
等效應力
或
若
則有 (彈性狀態下)
u *do, par, ival, fval, inc 定義一個do循環的開始
par: 循環控制變量
ival, fval, inc:起始值,終值,步長(正,負)
u *enddo 定義一個do循環的結束
u *if,val1, oper, val2, base: 條件語句
val1, val2: 待比較的值(也可是字符,用引號括起來)
oper: 邏輯操作(當實數比較時,誤差為1e-10)
eq, ne, lt, gt, le, ge, ablt, abgt
base: 當oper結果為邏輯真時的行為
lable: 用戶定義的行標志
stop: 將跳出anasys
exit: 跳出當前的do循環
cycle: 跳至當前do循環的末尾
then: 構成if-then-else結構
注意:不允許跳出、跳進一個do,if循環至label句 ?
8 理論手冊
1.方程組解法:(1)直接解法;(2)迭代解法
(1) 直接解法:a.稀疏矩陣法;b. 波前解法
a. 稀疏矩陣法:占內存大,但運算次數少;通過變換剛度矩陣的順序使得非零元素最少
b. 波前解法: 占內存小
波前是指在還沒有一個單元被解完的時候激活的方程數?
(2) 迭代解法:JCG法;PCG法;ICCG法
JCG法:可解實數、對稱、非對稱矩陣
PCG法:高效求解各種矩陣(包括病態),但僅解實、對稱矩陣
ICCG法:類似JCG,但更強
2. 應變密度,等效應變,應力密度,等效應力
(1)應變密度(strain intensity)
應變密度
是三個主應變
(2)等效應變
有效泊松比 :用戶由avprin 命令設定;0(如果不設定)
(3)應力密度(stress intensity)
應力密度
(4) 等效應力
等效應力
或
若
則有 (彈性狀態下)
Fini(退出四大模塊,回到BEGIN層)
/cle (清空內存,開始新的計算)
1 定義參數、數組,并賦值.
u dim, par, type, imax, jmax, kmax, var1, vae2, var3 定義數組
par: 數組名
type: array 數組,如同fortran,下標最小號為1,可以多達三維(缺省)
char 字符串組(每個元素最多8個字符)
table
imax,jmax, kmax 各維的最大下標號
var1,var2,var3 各維變量名,缺省為row,column,plane(當type為table時)
2 /prep7(進入前處理)
2.1 定義幾何圖形:關鍵點、線、面、體
u csys,kcn
kcn , 0 迪卡爾zuobiaosi
1 柱坐標
2 球
4 工作平面
5 柱坐標系(以Y軸為軸心)
n 已定義的局部坐標系
u numstr, label, value 設置以下項目編號的開始
node
elem
kp
line
area
volu
注意:vclear, aclear, lclear, kclear 將自動設置節點、單元開始號為最高號,這時 如需要自定義起始號,重發numstr
u K, npt, x,y,z, 定義關鍵點
Npt:關鍵點號,如果賦0,則分配給最小號
u Kgen,itime,Np1,Np2,Ninc,Dx,Dy,Dz,kinc,noelem,imove
Itime:拷貝份數
Np1,Np2,Ninc:所選關鍵點
Dx,Dy,Dz:偏移坐標
Kinc:每份之間節點號增量
noelem: “0” 如果附有節點及單元,則一起拷貝。 #p#分頁標題#e#
“1”不拷貝節點和單元
imove: “0” 生成拷貝
“1”移動原關鍵點至新位置,并保持號碼,此時(itime,kinc,noelem)被忽略
注意:MAT,REAL,TYPE 將一起拷貝,不是當前的MAT,REAL,TYPE
u A, P1, P2, ……… P18 由關鍵點生成面
u AL, L1,L2, ……,L10 由線生成面
面的法向由L1按右手法則決定,如果L1為負號,則反向。(線需在某一平面內坐標值固定的面內)
u vsba, nv, na, sep0,keep1,keep2 用面分體
u vdele, nv1, nv2, ninc, kswp 刪除體
kswp: 0 只刪除體
1 刪除體及面、關鍵點(非公用)
u vgen, itime, nv1, nv2, ninc, dx, dy, dz, kinc, noelem, imove 移動或拷貝體
itime: 份數
nv1, nv2, ninc:拷貝對象編號
dx, dy, dz :位移增量
kinc: 對應關鍵點號增量
noelem,:0:同時拷貝節點及單元
1:不拷貝節點及單元
imove: 0:拷貝體
1:移動體
u cm, cname, entity 定義組元,將幾何元素分組形成組元
cname: 由字母數字組成的組元名
entity: 組元的類型(volu, area, line, kp, elem, node)
u cmgrp, aname, cname1, ……,cname8 將組元分組形成組元集合
aname: 組元集名稱
cname1……cname8: 已定義的組元或組元集名稱
u cmlist,name
u cmdele,name
u cmplot, label1
2.2 定義幾個所關心的節點,以備后處理時調用節點號。
u n,node,x,y,z,thxy, thyz, thzx 根據坐標定義節點號
如果已有此節點,則原節點被重新定義,一般為最大節點號。
2.3 設材料線彈性、非線性特性
u mp,lab, mat, co, c1,…….c4 定義材料號及特性
lab: 待定義的特性項目(ex,alpx,reft,prxy,nuxy,gxy,mu,dens)
ex: 彈性模量
nuxy: 小泊松比
alpx: 熱膨脹系數
reft: 參考溫度
reft: 參考溫度
prxy: 主泊松比
gxy: 剪切模量
mu: 摩擦系數
dens: 質量密度
mat: 材料編號(缺省為當前材料號)
c 材料特性值,或材料之特性,溫度曲線中的常數項
c1-c4: 材料的特性-溫度曲線中1次項,2次項,3次項,4次項的系數
u Tb, lab, mat, ntemp,npts,tbopt,eosopt 定義非線性材料特性表
Lab: 材料特性表之種類
Bkin: 雙線性隨動強化
Bis 雙線性等向強化
Mkin: 多線性隨動強化(最多5個點)
Mis 多線性等向強化(最多100個點)
Dp: dp模型
Mat: 材料號
Ntemp: 數據的溫度數
對于bkin: ntemp缺省為6
mis ntemp缺省為1,最多20
bis ntemp缺省為6,最多為6
dp: ntemp, npts, tbopt 全用不上
Npts: 對某一給定溫度數據的點數
u TBTEMP,temp,kmod 為材料表定義溫度值
temp: 溫度值
kmod: 缺省為定義一個新溫度值
如果是某一整數,則重新定義材料表中的溫度值
注意:此命令一發生,則后面的TBDATA和TBPT均指此溫度,應該按升序
若Kmod為crit, 且temp為空,則其后的tbdata數據為solid46,shell99,solid191中所述破壞準則
如果kmod為strain,且temp為空,則其后tbdata數據為mkin中特性。
u TBDATA, stloc, c1,c2,c3,c4,c5,c6
給當前數據表定義數據(配合tbtemp,及tb使用)
stloc: 所要輸入數據在數據表中的初始位置,缺省為上一次的位置加1
每重新發生一次tb或tbtemp命令上一次位置重設為1,
(發生tb后第一次用空閑此項,則c1賦給第一個常數)
u tbpt, oper, x,y 在應力-應變曲線上定義一個點
oper: defi 定義一個點
dele 刪除一個點
x,y:坐標
2.4 設置單元類型及相應KEYOPT
u ET, itype, ename, kop1……kop6, inopr 設定當前單元類型
Itype:單元號
Ename:單元名設置實常數
u Keyopt, itype, knum, value
itype: 已定義的單元類型號
knum: 單元的關鍵字號
value: 數值
注意:如果 ,則必須使用keyopt命令,否則也可在ET命令中輸入
2.5 設置網格劃分,劃分網格
2.5.1 映射網格劃分
1.面映射網格劃分
條件:a. 3或4條邊 #p#分頁標題#e#
b.面的對邊必須劃分為相同的單元或其劃分與一個過渡形網格的劃分相匹配
c. 該面如有3條邊,則劃分的單元不必須為偶數,并且各邊單元數相等
d. mahkey
e. mshpattern
* 如果多于四條邊,可將線合并成Lcomb
可用amap命令,先選面,再選4個關鍵點即可
* 指定面的對邊的分割數,以生成過渡映射四邊形網格,只適用于有四條邊的面?
相關文章
- 2021-09-08BIM技術叢書Revit軟件應用系列Autodesk Revit族詳解 [
- 2021-09-08全國專業技術人員計算機應用能力考試用書 AutoCAD2004
- 2021-09-08EXCEL在工作中的應用 制表、數據處理及宏應用PDF下載
- 2021-08-30從零開始AutoCAD 2014中文版機械制圖基礎培訓教程 [李
- 2021-08-30從零開始AutoCAD 2014中文版建筑制圖基礎培訓教程 [朱
- 2021-08-30電氣CAD實例教程AutoCAD 2010中文版 [左昉 等編著] 20
- 2021-08-30電影風暴2:Maya影像實拍與三維合成攻略PDF下載
- 2021-08-30高等院校藝術設計案例教程中文版AutoCAD 建筑設計案例
- 2021-08-29環境藝術制圖AutoCAD [徐幼光 編著] 2013年PDF下載
- 2021-08-29機械AutoCAD 項目教程 第3版 [繆希偉 主編] 2012年PDF