第6章 在patran界面中建立dytran模型及后處理的操作過程實(shí)例6

1、定義用于儲(chǔ)存歐拉單元的壓力、密度及耦合能的時(shí)間歷程文件。首先創(chuàng)建一個(gè)叫做“two_euler_elem”的組將兩個(gè)分別處于兩個(gè)撞擊中心點(diǎn)附近的歐拉單元包含進(jìn)來。Group/Create/New Group Name=alegrids/Entity Selection=拾取兩個(gè)分別處于兩個(gè)撞擊中心點(diǎn)附近的歐拉單元/ApplyAnalyze/Input Deck/Translate/Output Request/Result Name=

2、pressure/File Type=Time History/Result Type=Element Output/Time for Outptu=0 thru end by 1.0E-6/Add/Select Groups for Output=“two_euler_elem/Result Types=density,specific internal energy,pressure/Apply/Ok/Apply,數(shù)據(jù)文件產(chǎn)生后，

3、用文本編輯器加以修改，將PATRAN界面中不支持的卡片包含進(jìn)來。... $$ -------- Contact : master-slave-surface$$CONTACT 17 SURF SURF 1 2 +A012850$+A012850 V4 BOTH

4、 +A012851$+A012851 +A012852$+A012852 ON ALE,1,2,1...$$----------- S

5、ETS ------------$Setname: alegrids with Set ID: 6 contains Nodes$SET1 6 579 THRU 8670ALEGRID,1,0.,1E20,SPECIAL,,,,,++,579,THRU,8670$...,提交作業(yè)。,創(chuàng)建一個(gè)新的PATRAN數(shù)據(jù)庫(kù)名為“post_bird”。讀入模型及結(jié)果的檔案文件。同時(shí)讀入時(shí)間歷

6、程文件。步驟如下：Analysis/Read Archive File/Model And Results /Translate/Select Archive Files/“MBIRD_BIRD_0”/Add/ “MBIRD_PLATE_0”/Add/Apply/ApplyAnalysis/Read Time History/Result/Translate/Select Time History Files/“MBIRD_PRE

7、SSURE_0.THS”/Add/Apply/Apply,通過創(chuàng)建關(guān)于應(yīng)力、變形及鳥體的等值面工具對(duì)結(jié)果進(jìn)行動(dòng)畫顯示。步驟如下：Insight/Create/IsoSurface/Result Selection/選中所有情況/Update Results/Isosurface Results=FMAT (Material Fraction)/Animation Attributes/點(diǎn)開Enable Animation/Global

8、 Variable=Load Case Index/OK/Result Options/Select Default Load Cases=拾取最后一步/OK/OApply,創(chuàng)建應(yīng)力云圖。步驟如下：Create/Fringe/Result Selection/Fringe Results=拾取 EFFST-OUT/Animation Attributes/Global Variable=Load Case Index/OK/OK/Ap

9、ply,,創(chuàng)建變形工具。步驟如下：Create/Deformation/Result Selection/Deformation Result=Dispacement/Animation Attributes/Global Variable=LOADCASE INDEX/OK/Deformation Attributes/Scale=True/OK/OK/Apply,進(jìn)行動(dòng)畫顯示。步驟如下：Insight Control （工具欄）

10、/Animation Control/Setup/（可以改變動(dòng)畫圖幅數(shù)）/Animate,畫出壓力的時(shí)間歷程。步驟如下：XY PLOT/Create/XYWindow/Enter XYWindow Name = w1 /ApplyPost/Curve/Post Curves =選中PRESSURE_ELEMENT_PRESSURE_MBIRD.curve6/Apply,6.7 圓柱形鈑金件沖壓成形分析

11、 有關(guān)工藝數(shù)據(jù)如下：,板料的有關(guān)材料參數(shù)如下：,在PATRAN中創(chuàng)建凸模的幾何。步驟如下：Group/Create/New Group Name=upper_die/ApplyGeometry/Create/Point/XYZ/Point Coordinate List=[0 0 0]/Apply/ Point Coordinate List=[37 0 0]/Apply Point Coordinate List

12、=[50 13 0]/Apply Point Coordinate List=[50 50 0]/Apply,Geometry/Create/Curve/Point/2Point/Starting Point List=Point 1/Ending Point List=Point 2/Apply/2D Arc2Point/Radius/Construction Plane List=Coord0.3/Arc Radius=13/

13、Starting Point List=Point 2/Ending Point List=Point 3/Apply/Starting Point List=Point 3/Ending Point List=Point 4/Apply,Geometry/Create/Surface/Revolve/Axis=Coord0.2/Total Angle=90/Curve List=Curve 1/Apply/Curve List=Su

14、rface 1.2/Apply/Curve List= Surface 2.2 /Apply/Curve List= Surface 3.2 /ApplyAngle=360/Curve List=Curve 2 3/Apply,在PATRAN中創(chuàng)建凹模的幾何。步驟如下：Group/Create/New Group Name=lower_die/ApplyGeometry/Create/Point/XYZ/ Point Coordi

15、nate List=[51.25 -60 0]/Apply/Point Coordinate List=[51.25 -6 0]/Apply/Point Coordinate List=[56.25 -1 0]/Apply/Point Coordinate List=[105.05 -1 0]/Apply/Point Coordinate List=[56.25 -6 0]/ApplyGeometry/Create/Curv

16、e/Point/2Point/Starting Point List=Point 5/Ending Point List=Point 6/Apply/Starting Point List=Point 7/Ending Point List=Point 8/Apply/2D ArcAngles/Radius=5/Start Angle=90/End Angle=180/Center Point List=Point 9/Apply,

17、Geometry/Create/Surface/Revolve/Axis=Coord0.2/TotalAngle=360/Curve List=Curve 4 6 5/Apply,用同樣的方法創(chuàng)建壓邊圈及板料的幾何。步驟如下：Group/Create/New Group Name=lower_die/ApplyGeometry/Create/Point/XYZ/ Point Coordinate List=[56.25 –0.5

18、 0]/Apply/Point Coordinate List=[105.05 –0.5 0]/Apply/Curve/Point/2Point/Starting Point List=Point 10/Ending Point List=Point 11/ApplyGeometry/Create/Surface/Revolve/Axis=Coord0.2/TotalAngle=360/Curve List=Curve 7/Ap

19、ply,Group/Create/New Group Name=lower_die/ApplyGeometry/Create/Point/XYZ/ Point Coordinate List=[0 –0.75 0]/Apply/Point Coordinate List=[100 –0.75 0]/Apply/Curve/Point/2Point/Starting Point List=Point 12/Ending Point

20、 List=Point 13/Apply/Surface/Revolve/Axis=Coord 0.2/Total Angle=90/Curve List=Curve 8/Apply/Curve List=Surface 11.2/Apply/Curve List= Surface 12.2 /Apply/Curve List= Surface 13.2 /Apply,在凸模幾何面上創(chuàng)建網(wǎng)格種子點(diǎn)。步驟如下：Group/Post/選

21、中“upper_die”/ApplyFinite Elements/Create/Mesh Seed/Uniform/Number of Elements=20/Curve List=Surface 1.2 2.2 3.2 4.2 1.3 2.3 3.3 4.3/Apply/Number of Elements=10/Curve List=Curve 2 3/Apply /Number of Elements=80/Curve L

22、ist=Surface 5.3 6.3/Apply,劃分有限元網(wǎng)格：Finite Elements/Create/Mesh/Surface/Quad/ISOMesh/Quad4/Curve List=Surface 1:6/Apply,消除重復(fù)結(jié)點(diǎn)：Finite Element/Equivalence/All/Tolerance Cube/Equivalencing Tolerance=0.005/Apply,檢查板元的法線方向:

23、 Finite Elements/Verify/Element/Normals/Draw Normal Vectors/Test Control=Display Only/Apply,在凹模幾何面上創(chuàng)建網(wǎng)格種子點(diǎn)。步驟如下：Group/Post/選中“l(fā)ower_die”/ApplyFinite Elements/Create/Mesh Seed/Uniform/Number of Elements=5/Curve L

24、ist=Curve 4:6/Apply /Number of Elements=80/Curve List=Surface 7.1 8.1 9.3/Apply,劃分有限元網(wǎng)格：Finite Elements/Create/Mesh/Surface/Quad/ISOMesh/Quad4/Curve List=Surface 7:9/Apply,消除重復(fù)結(jié)點(diǎn)：Finite Element/Equivalence/All/Tolera

25、nce Cube/Equivalencing Tolerance=0.005/Apply,檢查板元的法線方向: Finite Elements/Verify/Element/Normals/Color Code Element/Test Control=Display Only/Apply從圖中可以看出，網(wǎng)格中單元法線方向不一致。,進(jìn)一步檢查板元的法線方向: Finite Elements/Verify/

26、Element/Normals/Draw Normal Vectors/Test Control=Display Only/Apply從圖中可以看出，倒角部分單元法線指向外。,調(diào)整板元的法線方向: Finite Elements/Verify/Element/Normals/Draw Normal Vectors/Test Control=Reverse Elements/ Guiding Element=Elemen

27、t 3620/Apply,再次檢查板元的法線方向: Finite Elements/Verify/Element/Normals/Color Code Element/Test Control=Display Only/Apply從圖中可以看出，網(wǎng)格中單元法線方向一致。,在壓邊圈幾何面上創(chuàng)建網(wǎng)格種子點(diǎn)。步驟如下：Group/Post/選中“binder”/ApplyFinite Elements/Create/Me

28、sh Seed/Uniform/Number of Elements=5/Curve List=Curve 7/Apply /Number of Elements=80/Curve List=Surface 10.3/Apply,劃分有限元網(wǎng)格：Finite Elements/Create/Mesh/Surface/Quad/ISOMesh/Quad4/Curve List=Surface 10/Apply,消除重復(fù)結(jié)點(diǎn)：Fin

29、ite Element/Equivalence/All/Tolerance Cube/Equivalencing Tolerance=0.005/Apply,檢查板元的法線方向: Finite Elements/Verify/Element/Normals/Draw Normal Vectors/Test Control=Display Only/Apply從圖中可以看出，網(wǎng)格中單元法線方向一致，指向上。,在板料幾何面

30、上創(chuàng)建網(wǎng)格種子點(diǎn)。步驟如下：Group/Post/選中“plate”/ApplyFinite Elements/Create/Mesh Seed/Uniform/Number of Elements=40/Curve List=Surface 11:14.4 11:14.3/Apply,劃分有限元網(wǎng)格：Finite Elements/Create/Mesh/Surface/Quad/ISOMesh/Quad4/Curve Li

31、st=Surface 11:14/Apply,消除重復(fù)結(jié)點(diǎn)：Finite Element/Equivalence/All/Tolerance Cube/Equivalencing Tolerance=0.005/Apply,用MATRIG材料本構(gòu)模型定義凸模的材料特性:Materials/Create/Isotropic/Manual Input/Material name=rigid_upper/Input Properties/

32、Constitutive Model=Rigid(MATRIG)/ElementType=Shell/Density=7.8E-6/Elastic Modulus=2.1E8/Poisson Ratio=0.3/Ok/Apply,定義凸模的物理特性:Properties/Create/2D/Shell/Property Set Name=upper/Input properties/Material Name=rigid_upper

33、 /Shell Thickness=1./OK/Select Members = 拾取組“upper_die”中的所有幾何面/Apply,用MATRIG材料本構(gòu)模型定義凹模的材料特性:Materials/Create/Isotropic/Manual Input/Material name=rigid_upper/Input Properties/Constitutive Model=Rigid(MATRIG)/ElementTyp

34、e=Shell/Density=7.8E-6/Elastic Modulus=2.1E8/Poisson Ratio=0.3/Ok/Apply,定義凹模的物理特性:Properties/Create/2D/Shell/Property Set Name=upper/Input properties/Material Name=rigid_lower /Shell Thickness=1./OK/Select Members = 拾取組

35、“l(fā)ower_die”中的所有幾何面/Apply,用MATRIG材料本構(gòu)模型定義壓邊圈的材料特性:Materials/Create/Isotropic/Manual Input/Material name=rigid_upper/Input Properties/Constitutive Model=Rigid(MATRIG)/ElementType=Shell/Density=7.8E-6/Elastic Modulus=2.1E

36、8/Poisson Ratio=0.3/Initial Vel. Of CG(Vy)=-1000./Ok/Apply,定義壓邊圈的物理特性:Properties/Create/2D/Shell/Property Set Name=binder/Input properties/Material Name=rigid_binder/Shell Thickness=1./OK/Select Members = 拾取組“binder”中的所

37、有幾何面/Apply,用SHEETMAT材料本構(gòu)模型定義板料的材料特性:Materials/Create/3d Orthotropic/Manual Input/Material name=metal/Input Properties/Constitutive Model=Sheetmaterial(SHEEETMAT)/Elasticity Type=Isotropic/Yield Criterion Type=Normal Ani

38、sotropic/Hardening Rule Type=Normal Anisotropic /Density=8.413E-4/Elastic Modulus=1.E8/Poisson Ratio=0.33/X-Coord. of rolling direction=1./Y-Coord. of rolling direction=0./Z-Coord. of rolling direction=0./a=0./b=895000./

39、c=0.000294/N=0.412/R0=0.85/Apply,定義板料的物理特性:Properties/Create/2D/Shell/Property Set Name=plate/Input properties/Material Name=metal/Shell Thickness=0.7/OK/Select Members = 拾取組“plate”中的所有幾何面/Apply,下面定義工件剛性體的約束、強(qiáng)迫運(yùn)動(dòng)以及載荷。為

40、了定義凸模的速度-時(shí)間歷程，用“場(chǎng)”定義一個(gè)時(shí)間函數(shù)：Fields/Field Name=upper_die_velocity/Time/Input Data/(0 0);(1.E-3 0);(2.E-3 4.E3);(11.E-3 4.E3);(12.E-3 0); (1 0)/OK/Apply,凹模的固定：Loads/BCs/Create/Rigid Body Object/New Set Name=rigi

41、d_lower/Input Data /Select Rigid Material=rigid_lower /選中UX UY UZ RX RY RZ/Select Application Region/點(diǎn)擊鼠標(biāo)任選板上一結(jié)點(diǎn)/OK/Apply,定義凸模在三個(gè)轉(zhuǎn)動(dòng)方向的約束以及在X、Z方向的平動(dòng)的約束，在-Y方向的強(qiáng)迫運(yùn)動(dòng)：Loads/BCs/Create/Rigid Body Object/New Set Name=rigid_upp

42、er/Input Data/Select Rigid Material=rigid_upper/選中UX UZ RX RY RZ/Enforced Transl. Vel. Vector=/Time Dependence=在“Time Dependent Fields”框中選中“upper_die_velocity”/Select Application Region/點(diǎn)擊鼠標(biāo)任選板上一結(jié)點(diǎn)/OK/Apply,定義壓邊圈在三個(gè)轉(zhuǎn)動(dòng)方向的

43、約束以及在X、Z方向的平動(dòng)的約束，在-Y方向的壓邊力：Loads/BCs/Create/Rigid Body Object/New Set Name=binder/Input Data/Select Rigid Material=binder/選中UX UZ RX RY RZ/Enforced Transl. Vel. Vector=/Force Vector=/OK/Select Application Region/點(diǎn)擊鼠標(biāo)任選板

44、上一結(jié)點(diǎn)/OK/Apply,定義板料與凸模之間的主從接觸（凸模為主，板料為從):Loads/BCs/Create/Contact/Element Uniform/Option=Master-Slave Surface/New Set Name=plate_upper/Input Data/Contact Deactivation Time=12.E-3/Select Application Region/Type=Master/Ele

45、ment Type=2D/Surface=Bottom/Geometry Filter=FEM/Select Entities=點(diǎn)擊鼠標(biāo)拖動(dòng)光標(biāo)將組“upper_die”中所有板元框入/Add/OK 然后將類型改為從屬型 /Type=Slave/Element Type=2D/Surface=Bottom/Geometry Filter=FEM/Select Entities=點(diǎn)擊鼠標(biāo)拖動(dòng)光標(biāo)將組“plate”中所有板元框入/Add/O

46、K/Apply,定義板料與凹模之間的主從接觸（凹模為主，板料為從):Loads/BCs/Create/Contact/Element Uniform/Option=Master-Slave Surface/New Set Name=plate_lower/Input Data/Contact Deactivation Time=12.E-3/Select Application Region/Type=Master/Element T

47、ype=2D/Surface=top/Geometry Filter=FEM/Select Entities=點(diǎn)擊鼠標(biāo)拖動(dòng)光標(biāo)將組“l(fā)ower_die”中所有板元框入/Add/OK 然后將類型改為從屬型 /Type=Slave/Element Type=2D/Surface=top/Geometry Filter=FEM/Select Entities=點(diǎn)擊鼠標(biāo)拖動(dòng)光標(biāo)將組“plate”中所有板元框入/Add/OK/Apply,定義板料

48、與壓邊圈之間的主從接觸： Loads/BCs/Create/Contact/Element Uniform/Option=Master-Slave Surface/New Set Name=plate_binder/Input Data/Contact Deactivation Time=12.E-3/Select Application Region/Type=Master/Element Type=2D/Surface=Botto

49、m/Geometry Filter=FEM/Select Entities=點(diǎn)擊鼠標(biāo)拖動(dòng)光標(biāo)將組“binder”中所有板元框入/Add/OK 然后將類型改為從屬型 /Type=Slave/Element Type=2D/Surface=Bottom/Geometry Filter=FEM/Select Entities=點(diǎn)擊鼠標(biāo)拖動(dòng)光標(biāo)將組“plate”中所有板元框入/Add/OK/Apply,定義有關(guān)程序運(yùn)行的控制參數(shù)：Analys

50、is/Analyze/Input Deck/Translate/Execution Controls/Execution Control Parameters/End Time=0.012/Time Step Size at Start=1.0E-7/OK/OK,定義用于儲(chǔ)存有關(guān)板料的輸出數(shù)據(jù)的檔案文件。步驟如下：Analyze/Input Deck/Translate/Output Requests/Result Name=Plat

51、e/File Type=Archive/Result Type=Element Output/Times for Output=0 thru end by 0.5E-3/Add/Select Groups for Output=Plate/Entity Type=Sublayer Variable/選中所有層 (Mid/Outer/Inner)/Result Types=effective plastic strain & ef

52、fective stress/Apply,定義用于儲(chǔ)存有關(guān)工件的輸出數(shù)據(jù)的檔案文件。步驟如下：Analyze/Input Deck/Translate/Output Requests/Result Name=Tools/File Type=Archive/Result Type=Element Output/Times for Output=0 thru end by 0.5E-3/Add/Select Groups for Outp

53、ut=Tools/Result Types=User-specified/Apply,定義再啟動(dòng)文件。步驟如下：Analyze/Input Deck/Translate/Output Requests/Result Name=rst/File Type=Restart File/Steps for Output=end/Add/OK產(chǎn)生數(shù)據(jù)文件：/Apply,運(yùn)行作業(yè)。,創(chuàng)建一個(gè)新的PATRAN數(shù)據(jù)庫(kù)名為“cyl_cup_post

54、”。讀入模型及結(jié)果的檔案文件。步驟如下：Analysis/Read Archive File/Model And Results /Translate/Select Archive Files/“CYL_CUP_PLATE_0”/Add/ “CYL_CUP_TOOLS_0”/Add/Apply/Apply,用Quick Plot顯示最后一步的變形。首先創(chuàng)建一個(gè)組存放板料上的所有單元：Group/Create/New Group

55、Name=plate/Entity Selection=Elem 1:4800(板料網(wǎng)格單元號(hào)）/Apply張貼“plate”組：/Post/選中“plate”組/Apply畫出變形：Results/Create/Quick Plot/選擇最后一步/Deformation Results=Displacement/Apply,用Quick Plot顯示最后一步的等效塑性應(yīng)變：Results/Create/Quick Plot/