杆冲击实验

1新建文件

首先点击File>new，新建如图所示的文件。其中Folder是文件的保存目录，通过Browse可以更改保存的目录。点击Folderlist可以设置添加和去除常用的保存目录（带加号的是添加常用的目录，减号是去除选中目录）。文件名为talor，head和description可以不写（有前面没有红色感叹号的可以忽略的）。

对称性采用的是2D模型（Axial表示轴对称，Planar表示的是平面对称）。Axial实际是将三维的旋转后的图形简化为2D模型。而Planar是用于解决二维平面应变问题。对于杆，采用轴对称模型。

2定义材料

选中Materials>load选中的材料为STEEL4340（装甲钢）和TANTALUM（钽）

选完材料后，可以通过review查看材料的相关参数。Modify修改相应的参数。Delete删除所选的材料（只是在本例子中删除加载的材料，不会永久删除库中的材料），通过copy选项可以很方便的将材料的状态方程失效模式等复制进另一材料，比如新建的材料，再通过修改少量与原来不同的参数即可定义新的材料，不过不要乱修改材料的参数，这些库中的材料的参数都是经过试验验证的。

3定义初始条件

选中Init.Cond.>new，定义泰勒杆的x方向速度100mm/ms，选中IncludeMaterial是指只要是这个材料都是100mm/ms，对于材料相同速度不同时就不行。同时，轴向没有旋转速度，若有可以在Radialvelocity设置旋转速度。

同样，设定后可以通过modify进行修改。

4设置边界条件

点击Boundaries>new,可以定义靶板最右边的x方向的速度为0（即右端固

定），可以通过modify修改。

5模型及网格

点击part，输入部件的名字，采用lagrange算法，点击next，进入模型建立坐标系。

建立一个taylor杆的part泰勒杆长100，宽10（采用对称所以dy=5）这是网格划分，表示I方向化5格，J方向100格。对泰勒杆的填充钽材料（勾选FillwithinitialConditionset，表示将初始定义的速度赋予泰勒杆，可以发现x-velocity变为了100新建靶板part注意靶板x起始点，靶板长20宽100（对称性）靶板的网格靶板固定不动没有初始速度条件不用勾选Fillwith...（选择4340钢材料)这是建立好的模型下面给靶板加边界条件。对于靶板的最右边我们希望是靶板最右方x方向没有速度的。我们选择图示中靶板，点击Boundary，我们希望的是最右边的J方向上是固定的。

为了显示出此边界，我们点击plots勾选Boundaries,可以看到边界条件在图中显示了。

通过勾选Grid可以查看网格状况，勾选Vector查看速度矢量线。Plot这个菜单主要是负责显示的。

6定义接触

通常来说，不同拉格朗日part之间要定义接触。依次点击Interaction>Largrange/Largrange>ExternalGap>calculate>check，然后一直点击确定。因为本例子采用的是largrange和largrange接触。外部间隙，采用自动计算并检验（计算的值一般是网格大小的十分之一）。

结果我们发现是这样的，这是因为在前面建模的时候我们没有考虑间隙，此时我们可以回到建模时。点击part，选中靶板，点击Zoning，点击Transformation>translate>Selectedpart选中靶板，在xdistance中填写0.11表示向x正方向移动0.11mm。使之大于或等于caculate计算出的间隙（CurrentIJKRange是指移动整个坐标，我们需要的是将杆与靶板之间制造出缝隙）。通过放大（tips：按住ctrl+鼠标左键框选需要放大的区域【和Ansys一样】）可以看出间隙了。再次回到Interaction，再点击check。

将靶板与杆产生大于check的间隙放大后的间隙重新check此时表示的是间隙符合要求。点击确定即可。

7设置求解步骤

点击controls。设置循环为100000（如果计算发现循环过小，可以再加），设置的求解时间是1ms（这是指作用过程是1ms），其他不变。如果有时候模型很间单，网格也很好了，但是就是经常出现能量错误是可以通过修改能量循环为较大值（比如1000000）。

点击output（输出选项），如图，表示的是保存的AD文件是从0ms开始，每隔0.01ms保存一次，直到1ms。

8运行及结果

点击run，（中途可以stop），最后结果如图。通过plot>contour点击那个大于号，选择需要显示的内容（图示的是显示MIS.STRESS）

这时我们发现自己可能需要某个点随时间应力速度或者其他的东西的值。这时我们可以点击Plots>Cycles。选择回到最初的0循环。通过定义高斯点来查看节点信息。

点击parts>Gauges>add。勾选InteractiveSelection。按住alt+鼠标右键选择你想要的点。图中给出仅供参考。如果觉得图标太大，影响捕捉。可以进入plots。勾选Gauges，点击右端的大于号。将Markscale改为0。

再次点击run。一直点确定，即覆盖原来求解数据。（推荐把原来计算产生的文件删了，只留下t=0时的文件）。建议在模型建立的一开始要考虑求解的东西，

而不是后来再重新修改。

9History历史曲线

运算结束后，点击History，点击Gaugepoint，点击Singlevariableplot,弹出SingleVariablePlot对话框，按住ctrl+鼠标右键选中所有的gauge点，yvar选择PRESSURE,XVar选择TIME点击勾后即可出现所有高斯点压力随时间变化的图像。

当然我们可以将不同使用MultipleVariablePlots将不同y变量值绘制在图中，如需要将gauge3y变量改为compress（这个单位不一样是没有意义的）。

我们比如要比较gauge1x方向速度和y方向速度随时间变化的情况。可以单击MultipleVariablePlots.通过Deleteplot选项删除所有的plot（可能剩最后一个删除不了，可以通过modify改为你需要的gauge点和相应的变量）然后addplot添加你需要的变量。

当我们需要将这些数据导出（比如导入到origin或者excel中），可以点击writetxtoutput，弹出如下的对话框。点击确定。我们打开计算保存目录（如果不知道自己保存路径，可以通过autodyn中file>open点击后你会发现弹出保存目录。用写字板打开Ident0-talor01.uhs即可得到数据。

可以将uhs后缀改为txt，新建一个excel文件，打开后点击数据>自文本。

选中重命名后的txt文件，点击确定，按照图示进行操作（因为uhs文件

中用逗号作为分割符号，所以要勾选上逗号，这样就能分列显示数据）

最后得到的数据如图，然后可以进行一些相应的修改和做一些数理分析。

有可能你觉得这些数据过多，那么在最开始还未计算的时候可以通过output，点击History展开菜单，可以设置按循环数，将Increment设置为50，那么再次运行过后的数据减少为原来的1/50。或者你选中Times，用作用时间来进行History数据的保存。

或者你运行结束的时候，可以通过History>选中Reduce如图，表示的是每隔一个循环移除1个循环的数据，输入3表示隔一个循环移除2个循环数据。（事实上，这个功能经常会出错，建议在一开始就设置好保存History的数据间隔。如果间隔过大，那么数据可能太少，有可能出现没有数据的现象）。

10视角

最后可以通过Plots，勾选Rotate，此时可以查看90°时的模型，这也是说明轴对称2D模型事实上是3D模型的一种转化。

通过History,Partsummaries，选中靶板，yvar下拉菜单选中TotMass，可

以查看靶板的质量为1.2299e6mg。而通过Materials>review。可以查看4340钢密度为7.83g/cm3。那么总质量M=3.14*502*20*7.83=1229310mg。这和仿真的计算是一样的。