patran&nastran問(wèn)題集錦

1、關(guān)于軟件安裝的問(wèn)題去掉了
2、請(qǐng)問(wèn)在PATRAN中輸出圖片能將黑色背景去掉
在display 下的color plat....下面調(diào)整。把上面得黑條,變成白得,點(diǎn)擊apply就行了
方法二:
1. 用文字編輯器開啟c:MSCpatran2003r2sharewaremscunsupported
utilitiesextra_filesbv_p3toolbar_ntgui.def
2. 選取所有的文字并復(fù)制。
3. 用文字編輯器開啟 c:MSCpatran2003r2p3toolbar.def ,在最後面的位置貼上剛剛復(fù)制的文字, 最後存檔離開
4. 復(fù)制c:MSCpatran2003r2sharewaremscunsupportedutilities icons*.bmp 到c:MSCpatran2003r2icons工具欄出現(xiàn)三個(gè)圖標(biāo),背景顏色輕松改變。
3、計(jì)算完畢后,只想顯示應(yīng)力超過(guò)某個(gè)值的單元,而其它單元不想顯示,如何設(shè)置?
tool-list-creat,方法選attribute,設(shè)f>你要求的應(yīng)力,apply以后選中在list里面的即為你要求的.再用plot/erase不顯示你不要的單元。
4、一個(gè)四邊形平板,一端的邊上節(jié)點(diǎn)6個(gè)自由度全約束住,另外一端上幾個(gè)節(jié)點(diǎn)上施加z方向強(qiáng)迫位移<, , 1E-5>,沒有別的條件。計(jì)算完以后看F06文件,看看那些節(jié)點(diǎn)的位移是否加上了!用的是loads/BCs中的creat-diaplacement,我很奇怪的是:我試了幾次這個(gè)強(qiáng)迫位移值,如0.1,0.01,0.001,0.0001,f06文件中顯示正確,節(jié)點(diǎn)位移值確實(shí)就是輸入值!但是這個(gè)值在變小時(shí)如1e-5,1e-6,F(xiàn)06文件顯示結(jié)果為0?。。「杏X好像是nastran的識(shí)別問(wèn)題,把10的-5次方一下的數(shù)默認(rèn)為0!
問(wèn)題出在translation parameter里面的一個(gè)參數(shù)numerical,幫助文件里面說(shuō)它用于比較兩個(gè)數(shù)是否相等,其默認(rèn)為1e-4.writing才是判斷一個(gè)數(shù)是否近似為零,默認(rèn)為1e-21.但實(shí)際上當(dāng)你給出的強(qiáng)迫位移量小于numerical時(shí),它就認(rèn)為近似為0,在bdf文件中就給忽略掉了.你修改numerical為1e-5,你上面說(shuō)的1e-5就可以算了。
5、自重怎么加到模型上去?
自重是在load/BCS里加的create->inertial load->element 在input data->load/BC set scale Factor [輸入加速度的值一般取9.8] Trans Accel(a1 a2 a3)<0 -1 0> (力是沿Y軸向下 )后就ok了
tools下面有個(gè)mass properties是計(jì)算模型質(zhì)量和慣量的,不知對(duì)你有沒有幫助
/solu
/output,mass,txt
psolve,elform
/output
finish
將在當(dāng)前目錄下生成mass.txt文件里面會(huì)么都有。
6、共振時(shí)的應(yīng)力
模態(tài)分析得到的位移是模態(tài)位移,沒有輸入載荷時(shí)的結(jié)果其值沒有意義,同一模態(tài)下模態(tài)坐標(biāo)位移之間的比值才有意義,因此,通過(guò)NORMALMODES分析的結(jié)果不表示真實(shí)的位移以及應(yīng)力等。
需要計(jì)算FREQENCY RESPONSE或者 TRANSIENT RESPONSE才能獲得真實(shí)應(yīng)力。
7、體單元的節(jié)點(diǎn)沒有轉(zhuǎn)動(dòng)自由度,只需考慮三個(gè)平移自由度。
8、我在用PATRAN做計(jì)算時(shí),分析出的結(jié)果文件很大,往往提示磁盤空間不足,結(jié)果提示需要980M空間。
你試試不用 patran直接遞交分析,進(jìn)入nastran 用手工遞交,可輸入scr=yes 之后進(jìn)行分析。因?yàn)橹苯舆f交會(huì)產(chǎn)生個(gè)dball文件 很大的 ,沒什么大用,用手遞交后,用scr=yes 就可以不用產(chǎn)生這個(gè)文件了, 另外所產(chǎn)生的dball文件可以delete它,它占的空間太大了
在patran下不進(jìn)行分析,不用full run 用analysic deck 之后運(yùn)行 nastran.ex桌面沒有 在msc 文件夾下找到,用它開打你剛才生成的文件就行 在最底一行輸入scr=yes 就行,之后和在patran下直接調(diào)用nastran的界面一樣。
patran向nastran遞交運(yùn)算時(shí)的參數(shù)可以在nastxxxx.rcf文件中設(shè)定,這個(gè)文件可以用文本編輯軟件修改,加一句scr=yes就行了,運(yùn)算式通用的參數(shù)都可以在這里面定制。
9、做模態(tài)分析需要密度。
10、如何畫橢圓
http://www.caesky.com/dispbbs.asp?BoardID=25&ID=1722
1) Create a circle, then use Transform/Curve/Scale to stretch it in one direction.
The easiest way to create an ellipse is with the Utilities menu of Patran. Utilities--Geometry--Create curve by function
Exemple : If you want to create an ellipse with a=110 and b=90 your equation will be: Y = SQRT(90**2*(1-(`X**2/110**2)))
2) 畫圓,將圓旋轉(zhuǎn)一個(gè)角度
具體旋轉(zhuǎn)角度:設(shè)a為長(zhǎng)軸長(zhǎng); b為短軸長(zhǎng)
則滿足關(guān)系 cos(angle)=(b/a)
旋轉(zhuǎn)該角度后,將其投影到園形原來(lái)的所在平面即得所需橢圓
3) 利用局部坐標(biāo)系
創(chuàng)建局部坐標(biāo)系,選cylindrical 1,設(shè)置長(zhǎng)短軸即可
4) 用 curve ?一??,然後在 move/scale ?不同比例,??curve 就可以.
5) 據(jù)說(shuō)通過(guò)creat-curve-conic可以畫橢圓
11、MPC
MPC是一個(gè)很好的工具
用來(lái)擬和多點(diǎn)自由度的
比如
要對(duì)一個(gè)截面上的多個(gè)點(diǎn)施加位移或者力
特別是對(duì)整個(gè)截面施加彎矩
就可以使用MPC
把多個(gè)點(diǎn)的位移
力
彎矩施加到一個(gè)點(diǎn)上就可以了
特別是對(duì)彎矩
你知道,對(duì)實(shí)體模型的橫截面施加一個(gè)彎矩
如果沒有MPC的話
會(huì)有多么地困難
12、模態(tài)分析的時(shí)候是不能夠分析應(yīng)力的,分析應(yīng)力要單獨(dú)用別的求解類型再分析。
14、注意轉(zhuǎn)速的單位是 n 轉(zhuǎn)/s,不是弧度
一些小技巧
1、在 Patran里 如何Move 一組Points 的位置 , 而不改變這組 Points 的 ID 編號(hào)?
Group/Transform/Translate的功能,
這樣不但編號(hào)不會(huì)變, 連property跟邊界條件都會(huì)保留
2、Patran如何執(zhí)行多次Undo?
所有Patran的操作步驟, 都記錄在最新的一個(gè)patran.ses.xx中,
如果需要多次undo, 可以?除最后不需要的步驟指令行,
再利用 File -> Session -> Play 的方式, 執(zhí)行改過(guò)的patran.ses.xx ,
這樣可以無(wú)限制的undo。
3、Patran中如何定義桿件之間的鉸接?
用的是rod元素,不需要定義鉸接,因?yàn)閞od元間本身就是以鉸接形式連接。
如果用bar或是beam,需在properties里的Pinned DOFs @ Node 1/2定義鉸接。
4、Patran 如何把不小心Equivalence的node分開
用 Utility/FEM-Elements/Separate Elements
在equivalence時(shí), 可以?選項(xiàng)切換為"List",
只選擇特定某些節(jié)點(diǎn)作equivalence, 可以避免不小心的失鍺。
5、Patran如何?Tri3單元轉(zhuǎn)換為Tri6單元
在Patran -> Element -> Modify/Element/Edit , ?Type選項(xiàng)打勾,
在Shape中選Tri, New Shape 選Tri6, 最后選取想要改變的Tri3單元,
6、Patran 如何定義材料庫(kù)
Patran除了可以直接讀取MSC.Mvision的材料材料庫(kù)外,
還可利用執(zhí)行Session File的方式,直接使用以前已經(jīng)定義過(guò)的材料。
編輯patran.ses.xx,將定義材料的PCL指令剪下,
粘貼到另外一個(gè)文件中(如mat.ses)。
之后便可直接由Patran的File/ Session/ Play來(lái)加入此一材料的定義。
也可以直接加入Patran菜單的中:
把剛剛定義的session file 復(fù)制到
C:MSCpatran2001r3sharewaremscunsupportedutilitiesdata_filesbv_material_data ,
之后就會(huì)在Utilities/ Material/ Material Session File Library 中出現(xiàn)剛才的材料名稱。.
7、Patran 的完整信息輸出:
在執(zhí)行Patran的時(shí)候出現(xiàn)齊怪的錯(cuò)誤信息時(shí),
可以先把Patran關(guān)閉, 接著啟動(dòng)DOS窗口,
在DOS下直接輸入patran -stdout ,執(zhí)行
再重新啟動(dòng)Patran, 會(huì)多了一個(gè)信息窗口。
16、超單元
超單元的應(yīng)用跟模態(tài)分析沒有關(guān)系,它只是有限元中的一種分支方法,用于解決運(yùn)算量過(guò)大的問(wèn)題的,現(xiàn)在基本上已經(jīng)不用了!
17、約束
我Patran/Nastran在做項(xiàng)目時(shí),對(duì)結(jié)構(gòu)應(yīng)力進(jìn)行分析,發(fā)現(xiàn)在結(jié)構(gòu)被約束的部位往往產(chǎn)生比較大的應(yīng)力,甚至是最大的應(yīng)力值就在這些部位出現(xiàn),好像這種結(jié)果并不太符合實(shí)際情況。請(qǐng)問(wèn)約束該如何加,才能消除這種影響?
這種現(xiàn)象很正常的,只要你的約束反映了實(shí)際情況,約束部位的過(guò)大應(yīng)力你可根據(jù)圣維南原理視而不見。
圣維南原理看看力學(xué)方面的書!如果一定要把約束區(qū)域分析準(zhǔn)確的話,你的約束一定要盡可能地與真實(shí)情況相符,航空經(jīng)常用的方法是加彈簧元,彈簧元的剛度是用的經(jīng)驗(yàn)值或者試驗(yàn)值!
18、MSC.Flds 2004
計(jì)算飛行載荷-flightload and dynamics,主要用來(lái)計(jì)算氣動(dòng)彈性問(wèn)題。
19、如何找到頻譜分析的響應(yīng)最大點(diǎn)?
在Tools--List---Create就可以將你想要的點(diǎn)的結(jié)果導(dǎo)入到一個(gè)組中,然后可以從結(jié)果文件中調(diào)出再在origin中編輯,或是直接顯示。
20、自由模態(tài)
所謂自由模態(tài),就是無(wú)約束模態(tài)
前6階都是剛體模態(tài)
21、能否設(shè)置計(jì)算結(jié)束后不輸出f06文件?
可以設(shè)置的,Solution Parameters里面Maximum Printed Lines =
就是f06輸出的最大行數(shù)
22、一個(gè)利用DMAP提取質(zhì)量、剛度矩陣的方法
對(duì)于一個(gè)結(jié)構(gòu)靜力的線性分析或非線性分析,總剛度矩陣是每次求解過(guò)程中所必須的。如果想再一次分析完成之后在結(jié)果中查看結(jié)構(gòu)的總剛度矩陣和總質(zhì)量矩陣,可以在nastran中加入下列代碼:
PARAM,EXTOUT,DMIGPCH
在patran中直接輸入,會(huì)在結(jié)果中得到一個(gè).pch文件,里面包含了所需的總剛度矩陣和總質(zhì)量矩陣。
若要在求解過(guò)程中輸出總剛度矩陣或總質(zhì)量矩陣,就需要用DMAP語(yǔ)言。下面給出一個(gè)簡(jiǎn)單的利用DMAP語(yǔ)言提取單元?jiǎng)偠染仃嚕傎|(zhì)量矩陣和總質(zhì)量矩陣的例子:
ID MSC cxh77 $
DIAG 8
TIME 5
$ BEGIN (開始)
SOL 100
MALTER 'MALTER:USERDMAP'
GP1 GEOM1,GEOM2,,,/GPL,EQEXIN,GPDT,CSTM,BGPDT,SIL,/S,N,LUSET/0/0 $
GP2 GEOM2,EQEXIN,EPT,,/ECT, $
PLTSET PCDB,BGPDT,ECT/PLTX,PLTPAR,GPSET,ELSET/S,N,NSLLS/S,N,JPLOT $
COND P1S,JPLOT $
LABEL P1S $
GP3 GEOM3,bgpdt,GEOM2,,,,,/SLT,ETT/0/0/0 $
TA1, ,ECT,EPT,BGPDT,SIL,ETT,CSTM,,,/EST,,GEI,GPECT,,/LUSET/-1/0/1/0/0 $
EMG EST,CSTM,MPT,DIT,GEOM2,,,,,,BGPDT,,,,/
KELM,KDICT,MELM,MDICT,BELM,BDICT/
S,N,NOKGG $
$
$ KELM為所需要的單元?jiǎng)偠染仃?br>
EMA GPECT,KDICT,KELM,BGPDT,SIL,CSTM,,/KGG,/ $
$
$ 組集單元?jiǎng)偠染仃?br>$
$ KGG為總剛度矩陣
EMA GPECT,MDICT,MELM,BGPDT,SIL,CSTM,,/MGG,/-1/V,Y,WTMASS $
$
$ MGG為總質(zhì)量矩陣
matprn kelm// $ print element stifiness matricies (打印單元?jiǎng)偠染仃?
matprn kgg// $ print global stifiness matrix (打印總剛度矩陣)
matprn mgg// $ print global mass matrix (打印總質(zhì)量矩陣)
ENDALTER
LINK USERDMAP
$ 下面為一段具體的結(jié)構(gòu)模型
CEND
TITLE=
SUBTITLE=
SUBCASE 1
LOAD = 1
SPC = 1
DISP = ALL
STRESS = ALL
BEGIN BULK
$
GRID,1,,0.,0.,0.
GRID,2,,5.,0.,0.
GRID,3,,10.,0.,0.
CROD,1,1,1,2
CROD,2,1,2,3
PROD,1,1,.2
MAT1,1,1.+7,,.32,2.7
FORCE 1,2,,1000.,1.,0.,0.
SPC1,1,123456,1
ENDDATA
將之保存為matrix.bdf直接用nastran分析調(diào)用,在生成的.f06文件中可以查看輸出矩陣:
一些片斷:
0 MATRIX KELM (GINO NAME 101 ) IS A DB PREC 2 COLUMN X 78 ROW RECTANG MATRIX.
0COLUMN 1 ROWS 1 THRU 58 --------------------------------------------------
ROW
1) 4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00
ROW POSITIONS 11 THRU 50 NOT PRINTED - ALL ARE NULL.
51) 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 4.0000D+05
0COLUMN 2 ROWS 1 THRU 58 --------------------------------------------------
ROW
1) 4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00
ROW POSITIONS 11 THRU 50 NOT PRINTED - ALL ARE NULL.
51) 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 4.0000D+05
0THE NUMBER OF NON-ZERO TERMS IN THE DENSEST COLUMN = 3
0THE DENSITY OF THIS MATRIX IS 3.85 PERCENT.
1 DECEMBER 6, 2004 MSC.NASTRAN 9/23/03 PAGE 8
0
0 MATRIX KGG (GINO NAME 101 ) IS A DB PREC 18 COLUMN X 18 ROW SYMMETRC MATRIX.
0COLUMN 1 ROWS 1 THRU 7 --------------------------------------------------
ROW
1) 4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 -4.0000D+05
0COLUMNS 2 THRU 6 ARE NULL.
0COLUMN 7 ROWS 1 THRU 13 --------------------------------------------------
ROW
1) -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 8.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00
11) 0.0000D+00 0.0000D+00 -4.0000D+05
0COLUMNS 8 THRU 12 ARE NULL.
0COLUMN 13 ROWS 7 THRU 13 --------------------------------------------------
ROW
7) -4.0000D+05 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 0.0000D+00 4.0000D+05
0COLUMNS 14 THRU 18 ARE NULL.
0THE NUMBER OF NON-ZERO TERMS IN THE DENSEST COLUMN = 3
0THE DENSITY OF THIS MATRIX IS 2.16 PERCENT.
1 DECEMBER 6, 2004 MSC.NASTRAN 9/23/03 PAGE 9
0
23、MSC.Nastran可形成如下輸出文件
Job_name.DBALL 包含數(shù)據(jù)庫(kù)運(yùn)行的永久數(shù)據(jù)
Job_name.MASTER 數(shù)據(jù)庫(kù)運(yùn)行的總詞典
Job_name.F04 包含數(shù)據(jù)庫(kù)文件信息和模塊執(zhí)行摘要
Job_name.F06 包含MSC.Nastran的分析結(jié)果,為文本文件。
Job_name.LOG 包含系統(tǒng)信息和系統(tǒng)錯(cuò)誤信息
Job_name.OP2 MSC.Nastran的分析結(jié)果,為二進(jìn)制文件
Job_name.XDB MSC.Nastran的分析結(jié)果,為二進(jìn)制文件
24、GRID POINT ID DEGREE OF FREEDOM MATRIX/FACTOR DIAGONAL RATIO MATRIX DIAGONAL
102 T1 -2.78284E+14 3.28730E+01
102 T2 -9.16132E+13 3.28730E+01
113 R3 -8.35181E+08 2.72727E-03
USER FATAL MESSAGE 9050 (SEKRRS)
RUN TERMINATED DUE TO EXCESSIVE PIVOT RATIOS IN MATRIX KLL
USER ACTION: CONSTRAIN MECHANISMS WITH SPCI OR SUPORTI ENTRIES OR SPECIFY PARAM,BAILOUT,-1 TO
在bdf文件中,找到solution parameter字段,加入“PARAM,BAILOUT,-1”再提交計(jì)算!
不過(guò)這種方法用于調(diào)試局部模型可以,調(diào)試總體模型是不能用的。我們工作的規(guī)定文件中明確說(shuō)了,計(jì)算總體模型不能夠用bailout卡,也不能打開automatic constraint選項(xiàng)!
你的加載板面的1,2,6自由度沒有剛度(即加載板沒有約束剛體位移)
你可以任意選擇一點(diǎn)加
spc1,1,126,**
同時(shí)在BEGIN BULK卡下一行加卡片
PARAM,K6ROT,100.0
即可
25、彈簧
彈簧的位置由兩個(gè)節(jié)點(diǎn)定義,彈簧是只能承受指定自由度方向剛度的單元,所以在定義彈簧的時(shí)候,還必須指定彈簧的自由度方向(UX,UY等)比如彈簧只能拉壓是在x方向,則定義DOF為UX。上述所說(shuō)的自由度依賴于你定義節(jié)點(diǎn)時(shí)用的坐標(biāo)系,也就是節(jié)點(diǎn)若定義在總體坐標(biāo)下,則DOF也是總體系下的方向。
上述所說(shuō)的自由度依賴于你定義節(jié)點(diǎn)時(shí)用的坐標(biāo)系,也就是節(jié)點(diǎn)若定義在總體坐標(biāo)下,則DOF也是總體系下的方向。這個(gè)節(jié)點(diǎn)的坐標(biāo)系是可以改的,可以用element->modify->node->edit之下就有改變坐標(biāo)系的命令。
對(duì)于一維的彈簧單元,兩端頭的自由度選擇,我認(rèn)為是為它選的可以釋放的(即可自由的)的自由度。可能描述的不是很清楚,你去IDEAS9.0的幫助文件里面看看,關(guān)于彈簧單元的解釋比較詳細(xì)。
在定義彈簧單元的property時(shí),有兩個(gè)框框dof at node 1:UX,dof at node 2:UX:節(jié)點(diǎn)一的自由度UX(約束X方向的平動(dòng)),節(jié)點(diǎn)2的自由度UX約束;除了x方向約束(受剛度影響)外,其他自由
CURVE是有方向性的,,每條CURVE生成的時(shí)候方向是固定的,這個(gè)可以從GEOMETRY里面進(jìn)行確認(rèn)。以CURVE生成網(wǎng)格的時(shí)候,沿著方向依次是NODE1和NODE2。如果直接生成一維單元先選的是NODE1后選的是NODE2。
模擬一個(gè)y方向可被壓縮彈簧:在用SPRING屬性的時(shí)候兩個(gè)自由都選則UY,然后在邊界條件里將你需要的節(jié)點(diǎn)全約束就可以了。STRING 是NASTRAN的使用的一個(gè)字符串。在這里指SPRING的方向。
零維彈簧指的是GROUND SPRING即地彈簧比如在結(jié)構(gòu)的某個(gè)部位需要施加一個(gè)對(duì)地的彈性約束,就可以直接使用這種單元屬性。當(dāng)然用一維的同樣可以模擬,用零維的比較方便而已。
定義彈簧元:property的1D單元中有一個(gè)項(xiàng)是彈簧,只要再輸入這個(gè)彈簧單元的方向和剛度值即可。注意彈簧單元的性質(zhì)只取決于作用方向和剛度系數(shù)這兩個(gè)參數(shù),而于單元的兩個(gè)節(jié)點(diǎn)位置無(wú)關(guān)。