定义哑板元网格物理特性:
Properties/Create/2D/Dummy Shell/Property Set Name=Dum/Select Members=拾取所有板元/Add/Apply
定义钢的材料特性。采用多项式状态方程(EOSPOL)作为本构 模型来定义钢的材料特性: Materials/Create/Isotropic/Manual Input/Material Name=Steel/Input Properties/Constitutive Model=ElasPlas (DMAT)/Element Type=Lagrangian Solid/Yield Model=Johnson-Cook/Strain Rate=N/A / Failure Model=None/Spallation Model = Spallation Pressure/Input properties….etc./Apply/Apply
输出管内空气的有关物理参数，然后生成数据文件： Output Requests/Result Name=Tube/File Type=Archive/Result Type=Element Output/Steps for Output= 0 thru end by 5/Add/Select Groups for Output=Tube/Entity Type=Eulerian Elements/Result Types=x-vel,y-vel,z-vel,density of material,specific internal energy,pressure/Apply/OK/Apply
密度 = 7890 系数 A1 = 1.642E11 系数 A2 = 2.307E11 系数 A3 = 2.23E11 系数 B0 = 1.67 （下页续）
剪切模量 = 8.18E10 静屈服应力 = 350E6 强化参数 =275E6 强化指数 = 0.36 应变率参数 = 0.022
输入钢的材料特性（续）。
消除重复结点： Fenite Elements/Equivalence/Apply
创建一个名为“Tube”的组将刚刚生成的单元与结点包含进去(只 是体元与相关结点)。 Group/Create/New Group name=Tube/Entity Selection=All Solids & Grids/Apply
定义输出请求。 输出活塞上的应力、应变： Output Requests/Result Name=Piston/File Type=Archive/Result Type=Element Output/Steps for Output= 0 thru end by 5/Add/Select Groups for Output=piston/Entity Type=Lagrangian/Result Types=effective stress & effective plastic strain/Apply
在活塞表面构造哑板元网格用于定义耦合面的定义。 Finite Elements/Create/Mesh/ Surface/Shape=Quad/Mesher=IsoMesh/Topology=Quad4 /Surface List=拾取活塞几何体(Solid)的所有面/Apply 创建一个名为“Couple”的组把所有活塞上的哑元包含进去。
然后定义区域内的物理参数初始值： Load&Boundary Conditions/Create/Init. Cond. Euler/Initial Values/New Set Name=ticval/Input Data/Select Euler Material=air/Density=1.2887/Apply
提交作业。
创建一个名为“piston_post”的新数据库进行后处理。 读入结果文件： Analysis/Read Archive Files/Object=Model And Results /Translate/Select Archive Files/依次选中所有有关的Archive 文件/Add（重复）/Apply
将活塞与其表面的哑板元的重复结点消除:
Finite Elements/Equivalence/Object=Group/Method=Tolerance Cube/Select Groups = Piston&Couple/Equivalencing Toleraபைடு நூலகம்ce = 0.005/Apply
通过Quick Plot画出X方向的速度在第80步的云图: Results/Create/Quick Plot/Select Result Cases=Cycle 80/Select Fringe Results=XVELEL/Apply
通过Quick Plot画出压力在第80的云图: Results/Create/Quick Plot/Select Result Cases=Cycle 80/Select Fringe Results=Pressure/Apply
最后完成欧拉网格初始状态的完整定义： Load&Boundary Conditions/Create/Init. Cond. Euler/Region Definition/New Set Name=TICEUL1/Input Data/在“Existing PEULER1 Sets”框中点击“peuler1”/在“Existing Shapes Sets” 框中点击“tube”/在“Existing Initial Values Sets”框中点击 “ticva1”/Level Indicator=1/Add Row/OK/Apply
沿-X方向将这个面平移4个单位，生成另一个面： Geometry/Transform/Surface/Translate/Translation Vector=<-4 0 0>/Surface List=Surface 1:4/Apply
构造1个三维实体（Solid): Geometry/Create/Solid/Surface/Option=2 Surface/Starting Surface List=Surface 9/Ending Surface List=Surface 10/Apply
现在创建活塞的模型。 首先构造1个面: Geometry/Create/Surface/XYZ/Vector Coordinates List=<0 2 –2>/Origin Coordinates List=(-3.9 –1 1)/Apply
提示: 在一般耦合中初始状态耦合面应当与欧拉单元有相交的地方 。因此在该例中初始时刻模型中的管与作为耦合面的圆柱面 有一点重叠，以使耦合面被激活。
定义管端口处的流场边界条件：
Load&Boundary Conditions/Create/Flow/New Set Name=tube_end/Input Data/Select Material=Air/OK/Select Application Region/Geometry/Select Solid Faces=Solid1:4.6/Add/OK/Apply
检查活塞表面的哑板元的法线方向:
Finite Elements/Verify/Element/Normals/Draw Normal Vectors/Test Control=Display Only/Apply
发现方向不一致，调整哑板元的法线方向:
Finite Elements/Verify/Element/Normals/Draw Normal Vectors/Test Control=Reverse Elements/Guiding Element=Element 1234/Apply
张贴组“couple”。定义耦合面：
Load&Boundary Conditions/Create/Coupling/New Set Name=gcoup/Input Data/Sandard/Inside/Reverse Normal/OK/Select Application Region/FEM/Select Entities=选中所有哑板元/Add/OK/Apply
定义活塞初速度： Load/Boundary Conditions/Create/Initial Conditions/New Set Name=piston_velocity/Input Data/Trans Veloc=<172.75,0.,0.>/Select Application Region/Select Geometry Entities=Solid 5/Add/OK/Apply
温度指数 = 参考应变率 比热 = 452 融化温度 = 室温 = 0 分离压力 =
1.0 = 1.0 1538 -1.0E10
定义活塞物理特性: Properties/Create/3D/Lagrangian Solids/Property Set Name=Piston_Steel/Input Properties /Material Name = Steel/OK/Select Memebers = 拾取 “Piston”组内的所有体元 /Add/Apply
做压力云图的动画 … Results/Create/Fringe/Select Result Cases=选中所有时间步 的结果/Select Fringe Result=PRESSURE/Animate/打开Display Attribute 图标/Display=Free Edges/关闭 Max/Min Label/Show On Deformed/Apply/ 点击打开Animation 图标/Select Global Variable=Load case Index/Apply