上傳人：7***

認(rèn)證信息

認(rèn)證類型：個(gè)人認(rèn)證

認(rèn)證主體：王**（實(shí)名認(rèn)證）

IP屬地：湖北

IP屬地：湖北 上傳時(shí)間：2022-01-11 格式：DOC 頁數(shù)：11 大小：22KB 積分：20 舉報(bào) 版權(quán)申訴

1、Ansys中單元?jiǎng)偠染仃?、總剛矩陣的提取一、單元?jiǎng)偠染仃嚨奶崛?DEBUG命令詳細(xì)說明:finish/clearPI=3.1415926w1=3w2=10w3=6w4=1.2r=.8t=0.08/PREP7!*ET,1,SHELL63R,1,tET,2,MASS21R,2,500,500,500,2000,2000,2000,!*UIMP,1,EX, , ,2e11UIMP,1,NUXY, , ,0.3,UIMP,1,DAMP, , ,0.2,UIMP,1,DENS, , ,7800,BLC4,0,0,w2,w1ESIZE,1.5,0,AMESH,allNSEL,S,LOC,X,0.0D,al

2、l, , , , , ,ALL, , , , ,allsel,allSFA,all,1,PRES,12FINISH/OUTPUT,cp,out, ! 將輸出信息送到cp.out文件/debug,-1,1 ! 指定輸出單元矩陣/SOLUSOLVEfinish/OUTPUT, TERM ! 將輸出信息送到output windows中這時(shí)用編輯器打開cp.out文件,可以看到按單元寫出的質(zhì)量、剛度等矩陣二、整體剛度矩陣的提取(有三種方法:用戶程序法、超單元法、HBMAT命令法1、用戶程序法:需要二次開發(fā)(略2、超單元法/soluantype,7 !substructuring分析類型seopt,m

3、atname,1 !設(shè)置文件名稱和剛度矩陣類型(剛度,質(zhì)量,阻尼等nsel,all !選擇所有節(jié)點(diǎn)m,all,all !定義所有節(jié)點(diǎn)自由度為主自由度solve !求解selist,matname,3 !列出整體剛度矩陣3、HBMAT命令法提取整體矩陣命令:HBMAT,fname,ext,-,form,matrx,rhs其中:Fname-輸出矩陣的路徑和文件名,缺省為當(dāng)前工作路徑和當(dāng)前工作文件名。ext-輸出矩陣文件的擴(kuò)展名,缺省為.matrix。form-定義輸出矩陣文件的格式,其值可取:=ASCII:ASCII碼格式;=BIN:二進(jìn)制格式。matrix-定義輸出矩陣的類型,其值可取:=STI

4、FF:輸出剛度矩陣?？捎糜趯懭肓?FULL文件的任何類型的分析。=MASS:輸出質(zhì)量矩陣。可用于特征值屈曲、子結(jié)構(gòu)分析、模態(tài)分析。=DAMP:輸出阻尼矩陣。僅用于有阻尼的模態(tài)分析。rhs-右邊項(xiàng)輸出控制(右邊項(xiàng)指用矩陣所表示方程的等號(hào)右端矢量,這里可為節(jié)點(diǎn)荷載向量,如rhs=YES則輸出,如rhs=NO則不輸出。模態(tài)分析時(shí),因僅LANB和QR法可生成完整的質(zhì)量矩陣,因此也僅采用這兩種方法時(shí)才可使用HBMAT命令得到質(zhì)量矩陣文件。Harwell-Boeing文件格式用HBMAT命令可輸出結(jié)構(gòu)剛度矩陣、質(zhì)量矩陣和阻尼矩陣,其文件記錄格式為大型稀疏矩陣的標(biāo)準(zhǔn)交換格式,采用索引存儲(chǔ)方法僅記錄矩陣的非零

5、元素。文件基本格式是前面有4或5行描述數(shù)據(jù),其后為單列矩陣元素值,說明如下:第1行:格式(A72,為文件頭的字符型解釋,如剛度矩陣或質(zhì)量矩陣等標(biāo)題。第2行:格式(5I14,分別表示該文件的總行數(shù)(不包括文件頭、矩陣列指針的總行數(shù)、矩陣行索引的總行數(shù)、矩陣元素?cái)?shù)值的總行數(shù)、右邊項(xiàng)總行數(shù)。第3行:格式(A3,11X,4I14,分別為矩陣類型、矩陣行數(shù)、矩陣列數(shù)、矩陣行索引數(shù)(對(duì)組裝后的矩陣,該值等于矩陣行索引數(shù)、單元元素?cái)?shù)(對(duì)組裝后的矩陣此值為0。第4行:格式(2A16,2A20,分別表示列指針格式、行索引格式、系數(shù)矩陣數(shù)值格式、右邊項(xiàng)數(shù)值格式。第5行:格式(A3,11X,2I14,A3各列分別表

6、示右邊項(xiàng)格式、應(yīng)用高斯起始矢量、應(yīng)用eXact求解矢量;兩個(gè)整數(shù)分別表示右邊項(xiàng)列數(shù)、行索引數(shù)。三個(gè)字符中的第1個(gè)字符可取:F-全部存貯(如節(jié)點(diǎn)荷載向量的全部元素、M-與系數(shù)矩陣相同方法。第6行后:矩陣元素值(單列。矩陣類型用3個(gè)字符表示,第1個(gè)字符可取:R-實(shí)數(shù)矩陣、C-復(fù)數(shù)矩陣、P-僅矩陣結(jié)構(gòu)(無元素?cái)?shù)值;第2個(gè)字符可取:S-對(duì)稱矩陣、U-不對(duì)稱矩陣、H-Hermitian 矩陣、Z-病態(tài)對(duì)稱矩陣;R-帶狀矩陣;第3個(gè)字符可取:A-組裝的矩陣、E-單元矩陣(未組裝。對(duì)稱矩陣只存儲(chǔ)下三角元素,如結(jié)構(gòu)剛度矩陣為對(duì)稱矩陣,Harwell-Boeing 格式則僅記錄下三角元素。命令流:/aux2fi

7、le,hbfile,fullhbmat,hbfile,txt,ascii,stiff,yesfinish!從hbfile.txt讀入數(shù)據(jù),并還原為滿矩陣存儲(chǔ)*dim,contline,5*vread,contline(1,hbfile,txt,5,1(5f14.0totcrd=contline(1ptrcrd=contline(2indcrd=contline(3valcrd=contline(4rhscrd=contline(5*vread,contline(1,hbfile,txt,4,2(a3,11x,4f14.0nrow=contline(2ncol=contline(3strline

8、=contline=*if,rhscrd,eq,0,thenls0=4*elsels0=5*endif*dim,pointr,ptrcrd*dim,rowind,indcrd*dim,values,valcrd*dim,rhsval,rhscrd*vread,pointr(1,hbfile,txt,ptrcrd,ls0(f14.0*vread,rowind(1,hbfile,txt,indcrd,ls0+ptrcrd(f14.0*vread,values(1,hbfile,txt,valcrd,ls0+ptrcrd+indcrd(d25.15*vread,rhsval(1,hbfile,txt

9、,rhscrd,ls0+ptrcrd+indcrd+valcrd (d25.15*dim,smatr,nrow,ncol*do,icol,1,ncolstacol=pointr(icolendcol=pointr(icol+1*do,irow,stacol,endcol-1truerow=rowind(irowsmatr(truerow,icol=values(irow*enddo*enddo*do,irow,1,nrow*do,icol,1,ncolsmatr(irow,icol=smatr(icol,irow*enddo*enddopointr=rowind=values=rhsval=i

10、col=irow=ls0=stacol=endcol=truerow=totcrd=ptrcrd=indcrd=valcrd=rhscrd=FINISH/cleFINISH/cle/PREP7ET,1,SOLID45mp,ex,1,1e7mp,nuxy,1,0.25KEYOPT,1,1,1KEYOPT,1,2,0BLC4,0,0,10,2 VOFFST,1,2FLST,5,12,4,ORDE,2FITEM,5,1FITEM,5,-12CM,_Y,LINELSEL, , , ,P51XCM,_Y1,LINECMSEL,_YLESIZE,_Y1,1, , , , , , ,1MSHAPE,0,3D

11、MSHKEY,1CM,_Y,VOLUVSEL, , , ,1CM,_Y1,VOLUCHKMSH,'VOLU' CMSEL,S,_YVMESH,_Y1CMDELE,_YCMDELE,_Y1CMDELE,_Y2FLST,2,9,1,ORDE,9FITEM,2,1FITEM,2,14FITEM,2,24FITEM,2,34 FITEM,2,47 FITEM,2,57 FITEM,2,68 FITEM,2,80 FITEM,2,90 !* /GO D,P51X, , , , , ,ALL, , , , , FLST,2,1,1,ORDE,1 FITEM,2,79 !* /GO F,P5

12、1X,FY,-1000 /solu /output,elemstiff,txt ! 進(jìn)入求解層，定義輸出單元?jiǎng)偠染仃囄募蛿U(kuò)展名 /debug,-1,1 ! 設(shè)置輸出單元?jiǎng)偠染仃嚭蛦卧奢d向量 solve /output finish ! 求解，設(shè)置輸出到終端 FINISH /OUTPUT,cp,out, ! 將輸出信息送到 cp.out 文件 /debug,-1,1 ! 指定輸出單元矩陣 /SOLU SOLVE finish /OUTPUT, TERM ! 將輸出信息送到 output windows 中 2-/solu antype,7 !substructuring 分析類型 seopt,matname,1 !設(shè)置文件名稱和剛度矩陣類型(剛度，質(zhì)量，阻尼等 nsel,all !選擇所有節(jié)點(diǎn) m,all,all !定義所有節(jié)點(diǎn)自由度為主自由度 solve !求解 selist,matname,3 !列出整體剛度矩陣 /OUTPUT,cp,out, ! 將輸出信息送到 cp.out 文件 /debug,-