有限元分析报告.docx

有限元分析报告.docx

《有限元分析报告.docx》由会员分享，可在线阅读，更多相关《有限元分析报告.docx（9页珍藏版）》请在冰豆网上搜索。

有限元分析报告

有限元仿真分析实验

一、实验目的

通过刚性球与薄板的碰撞仿真实验，学习有限元方法的基本思想与建模仿真的实现过程，并以此实践相关有限元软件的使用方法。

本实验使用HypeiMesh软件进行建模、网格划分和建立约束及载荷条件，然后使用LS-DYNA软件进行求解计算和结果后处理，计算出钢球与金属板相撞时的运动和受力情况，并对结果进行可视化。

二、实验软件

HypeiMesh>LS-DYNA

三、实验基本原理

本实验模拟刚性球撞击薄板的运动和受力情况。

仿真分析主要可分为数据前处理、求解计算和结果后处理三个过程。

前处理阶段任务包括：

建立分析结构的几何模型，划分网格、建立计算模型，确定并施加边界条件。

4.实验步骤

1、按照点一线一面的顺序创建球和板的几何模型

（1）建立球的模型：

在坐标（0,0,0）建立临时节点，以临时节点为圆心，画半径为5111111的球体。

（2）建立板的模型：

在tool-translate面板下node选择临时节点，选择Y-axis,magnitude输入5.5,然后点击translate+,letuin;再在2D—planes-squaie面板上选择Y-axis,B选择上一步移下來的那个节点，surfaceonly,size=30。

2、画网格

（1）画球的网格：

以球模型为当前pan,在2D—atiiomesh面板下，surfs选择前面建好的球面，elementsize设为0.5mm,meshtype选择quads,选择elemstocunentcomp.firstolder,interactive。

（2）画板的网格：

做法和设置同上。

3、对球和板赋材料和截面属性

（1）给球赋材料属性：

在materials面板内选择20号刚体，设置Rlio为

2.000e-08,E为200000,NU为0.30。

（2）给球赋截面属性：

属性选择SectShll,thickness设置为0.1,QR设为0。

（3）给板赋材料属性：

材料选择MATL1,其他参数：

Rlio为2.000e-08,E为100000'Nu为0.30,选择DoNotExport。

（4）给板赋截面属性：

截面选择SectShlbtluckiiess设为0.2。

其他参数:

SHRE为8.333-01,QR为0,T1为0.2。

（5）给板设置沙漏控制：

在Piopemes-Create面板下Cardmiage选择HouiOlass,IHQ为4,QM为0.100o更新平板。

4、加载边界条件

（1）将板上最外面的四行节点分别建成4个set。

（2）建立一个loadcollectoro

（3）Analysis—constraints面板中，设置SIZE为1,nodes通过bysets选择set_l、set_2、set_3、set_4,然后点击cieat即可，边界条件加载完毕。

5、建立载荷条件（给球一个3mm的位移）

（1）建立一个plot:

post—xyplots-plots-creatplot,然后点击return;

（2）在post-xyplots-editcuives面板中输入X{0,0,0.0001},Y{0,3}。

（3）给刚性球一个3mm的沿y正方向的位移：

cardimage设为PiciRgd,DOF为2,VAD为2,LCID为1,SF为1,option选择Rigido

6、接触处理

（1）做两个用于接触的segment：

在Analysis-set_segment面板中，Cardmiage选择setSegment,elems选择球这个pan。

重复操作对平板创建segmento要同时保证球的setsegment的方向朝外，plane的setsegment方向朝上。

（2）建立接触：

1对称接触：

在Aiialysis-mterfaces面板中,tvpe选择SuifaceToSuface,Cardmiage选择SingleSurfaceomaster的contactsurfs选择球的setgment;slave的contactsuifs选择平板的setgmento点击edit后，设置FS为10,FD为10,SFS为100,SFM为100,Automatic一项勾选OneWay。

2非对称接触：

在Analysis-interfaces面板中,type选择SmgleSurface,Cardmiage选择SingleSurfaceomaster的contactsurfs选择球的setgment;slave的contactsuifs选择平板的setgmento点击edit后，设置FS为10,FD为10,SFS为100,SFM为100,Automatic一项勾选smootho

7、定义控制卡片

在Aiialysis-contiolcards面板中，

（1）选择Contiol_Enegv,将hgen设置为2；

（2）选择ContioLTernuiiation,将ENDTIM设为0.0001s;

（3）选择Contiol_Time_step,将DTINIT设为l*10-6s,将TSSFAC设置为0.6;

（4）选择DATABASE_BINARY_D3PLOT,将DT设置为5*10-6;

（5）选择DATABASE_OPTION,将MATSUM设置为1*10-6,将RCFORC设置为l*10-6o

8、删除临时节点

在geom中选择tempnodes,点击node,选择all,然后点击clearo

9、节点重新排号。

tool中选择renumber,然后选择all,最后点击lenumbei。

10、将文件导出成KEY文件。

11、生成的KEY文件导入LS-DYNA中，并运行。

12、打开后处理程序Ls-Piepost,点击file-open-ls-dvnabinaiyplot,选择计算得到的d3plot结果。

13、点击下方运行按钮，观看碰撞动画。

选择右边的ASCII按钮，选择icforc后，点击左边的Load按钮,再选择S1或者Ma,然后选择最下方的Resultantforce,最后点击Plot,观看冲击力随时间的变化曲线。

五、实验结果

1、对称接触情况

对称接触定义下，冲击过程中薄板的弹塑性变形过程如图1至图5所示。

问c遵逊EF

■母Ellylfepu!

■鼻电知&r>l>ltv

tnt

Time=9.9604M06

93t*:

BKIB►仆►«

state-1；

ior

•nnwt*

•^nf.dJpkX

▼庐-金・禹右

m«电血&臨>怜

▲er，ae

fistfUndw«f

•^nf.dJpkX

6.4986e£05

MJ

X1Ur22tic1叶：

6M&5QX5t>t*：

14

sQF

30

Prwtou&SUtA

gWe-1；:

FastfUn3nM

9.9894e4）05

bjcri

，Tr：

•-殳曲

5tateV

iur

图5

冲击过程中刚性球与薄板层接触力时间历程如图6所示。

厲PlotWindow-

Time（E・03）

SRS

TltteScaleAttr

Alter

Print

Save

Load

Oper

Hide

Close

Quit

倉少B歼-ni£j丁

PE蚪門jj曲八UQ■翳妙厲为斥色丄

700e商圈^一jE*3一・n砂#•裁®l®le®lESa"

79

1010

Sts

阖囹▲）{■▼■▼-

d主wow.kKf2逬bxEX^Vti>上l£

。

※汪二匝細Z座吕盟H4怒尝•期stg坚ftHsaH旳fe•宀汉很屋舉極Ex#

S匝

«•S2

♦0—5一

feKM

・|,?

0的曙ll-tll®吩倉I0P醴9

-■wygiOCRM一0

图9

图10

图12

六、实验总结

根据对称接触定义与非对称接触定义的对比实验计算过程和实验结果可以看出，对称接触定义下刚性球撞击薄板的冲击现象更加明显，接触力更大，计算所花费的时间也更长，这与理论上双接触面与单接触面的结论是相同的。