我想看一级黄色大片_久久亚洲国产精品一区二区_久久精品免视看国产明星_91久久青青青国产免费

您的位置:網站首頁 > Ansys教程

ansys命令流----前后處理和求解常用命令之求解與后處理

時間:2010-11-14 11:34:33 來源:未知

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 以同樣工作名
#p#分頁標題#e#(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:on: 當solcontrol為on時缺省為on
off: 當solcontrol為off時缺省為off
1: 由程序選擇(當solcontrol為on且不發生autots命令時在 .log文件中紀錄“1”
注意:當使用自動時間步長時,也會使用步長預測器和二分步長
u NROPT, option,--,adptky 指定牛頓拉夫遜法求解的選項
OPTION: AUTO:程序選擇
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:程序執行時間(秒)限制,缺省為無窮
#p#分頁標題#e#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 在指定節點加集中荷載
#p#分頁標題#e#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. 盡量加面載,不加集中力,以免奇異點
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 顯示變形前和變形后的結構
#p#分頁標題#e#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 不顯示未變形的結構
#p#分頁標題#e#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線
“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 不顯示符號
#p#分頁標題#e#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) 等效應力
等效應力或若則有 (彈性狀態下)
#p#分頁標題#e#