有限元分析.docx

《有限元分析.docx》由会员分享，可在线阅读，更多相关《有限元分析.docx（12页珍藏版）》请在三一文库上搜索。

1、有限元上机分析报告学院：机械工程专业及班级：机械设计及其自动化08级7班姓名：王浩煜学号：20082798题目编号：21. 题目概况1.1 结构组成和基本数据结构：该结构为一个六根杆组成的桁架结构，其中四根杆组成了直径为800cm 的正方形，其他两根杆的两节点为四边形的四个角。材料：该六根杆截面面积均为100cm2,材料均为Q235,弹性模量为200GPa, 对于直径或厚度大于100mm的截面其强度设计值为190Mpa。载荷：结构的左上和左下角被铰接固定，限制了其在平面内x和y方向的位 移,右上角受到大小为 2000KN 的集中载荷。结构的整体状况如下图所示：2000KK800cm1.2 分析

2、任务该分析的任务是对该结构的静强度进行校核分析以验算该结构否满足强度要求。2. 模型建立2.1 物理模型简化及其分析由于该结构为桁架结构,故认为每根杆件只会沿着轴线进行拉压,而不会发生弯曲和扭转等变形。结构中每根杆为铰接连接，有集中载荷作用于最上方的杆 和最右方杆的铰接点。2.2 单元选择及其分析由于该结构的杆可以认为是只受拉压的杆件，故可以使用 LINK180 单元， 该单元是有着广泛工程应用的杆单元，它可以用来模拟桁架、缆索、连杆、弹簧 等等。这种三维杆单元是杆轴方向的拉压单元，每个节点具有三个自由度：沿节 点坐标系 X、Y、Z 方向的平动。就像铰接结构一样，不承受弯矩。输入的数据 有：两

3、个节点、横截面面积(AREA)、单位长度的质量(ADDMAS)及材料属性。 输出有：单元节点位移、节点的应力应变等等。由此可见， LINK180 单元适用 于该结构的分析。3.3 模型建立及网格划分(1)启动Ansys软件，选择PreferencesStructural，即将其他非结构菜单过滤掉。Peferences for GUI FilteringKEYW Preferences for GUI Filtering Individual discipline(s) to show in the GUI疋 Structural厂 Thermal r ANSYS FluidElectromag

4、netic:Note: If no individual disciplines arer FL0TRAN CFD厂 Magnetic-Nodal 厂 Magnetic-Edge I- High Frequency I- Electric selected they will all show.Discipline options h-Method(2)选择单元类型：选择 PreprocessorElement TypeAdd/Edit/Deletc Add,在出现的对话框中选择Link3d finit stn 180，即LINK180，点击“OK”(3)选择实常数：选择 Preprocess

5、orReal ConstantsAdd/Edit/DeleteAdd, 在出现的对话框中的Cross-sectional area中输入100,点击“OK”(4)材料定义： PreprocessorMaterial Props Material Models Structural LinearElasticIsotropic,出现以下窗口，在EX即弹性模量处输入2e7，点击 “ OK”。（5）建立线模型： PreprocessorModelingCreateKeypointsIn active CS， 由于坐标值在不输入情况下默认为0,故直接点击“Apply”则出现坐标为（0, 0）的 1 点

6、再在 X 项输入 800，则出现坐标为（80， 0）的 2 点，依次得到坐标 为（0, 80）和（80, 80）的3点和4点最后选择“OK”最后可在绘图区域看 到 4 个点。(6)划分网格：选择 PreprocessorMeshingMesh AttributesPicked Lines,选 择桁架结构的六根杆点击“OK”出现下图所示对话框，Real Constant Set Number 为1为要求的，故不需修改，直接点击“ OK”MAT Material numberNone defi门巳日TREAL Real constant set number1 -TYPE Element typ

7、e number1 UN KI 80JSECT Element sectionNone definedJPick Orientation Keypoint(s)I- NoA Line AttributesLATT Assign Attributes to Picked LinesApplyCancelHelp选择 PreprocessorMeshingMeshTool,在 Size Control 中的 global 中点击 Set进行单元设置，出现下图所示对话框，在No.of element divisions中填入1,即每个杆为一个单元。点击“0K”。选择PreprocessorMeshi

8、ngMeshTool Mesh,在弹出的对话框中，选择 “Pick All”，即将所有杆进行单元划分。此时选择ListElementsNodes+Attributes,即可列出单元，如下图：ELEM MAT TP REL ES SECNODES1112 113 114 115 116 11(7)施加载荷：1011 21012310134101141012410113施加约束：选择 SolutionDefined LoadsApplyStructuralDisplacement-On Nodes,选择节点1和2,在弹出的对话框中选择ALL DOF,点击“OK”施加力：选择 Solution-De

9、fined LoadsApplyStructural-Force/Moment On Nodes,输入数据，如下图所示：荷载施加完成后的图如下所示：3. 计算分析3.1 求解选择SolutionSolveCurrent LS进行求解，求解过程很快结束:八Note,*處U Solution is done 1CloseLJ选择 PlotCtrlsStyleSize And Shape，出现以下对话框，在 Display Of Element Shape Based On Real Constant Descriptions 后面的复选框中点击，贝“Off变为On” ,此时显示为实体模型。Size

10、 and Shape/SHRINK Shrink entities by/ESHAPE Display of elementshapes based on realconstant descriptionsSCALE Real constant multiplier/EFACET Facets/element edge/RATIO Di&tortion of GeometryWN Window numberRATOX X distortion ratioRATOY distortion ratio/CF0RMAT7I Component/Parameter Format N FIRST. N

11、LAST Widths/REPLOT Replot upon OK/Apply?Apply|o percentp* On01 facet/e，clgeiTCan匚已I再选择 General PostprocPlot ResultsDeformed Shape，在弹出的对话框中 选择Def+undef edge，点击“OK”，此时则可显示出桁架结构变化前后的位置, 如下图:DISPLACEMENTSIEP=1ANSYSAPR 21 2012 18:2:26TIME=1DMX =1.77042选择 General postprocList resultsNodal solution,在弹出的对di

12、spacement vector sum,列出的节点位移如下（单位为cm）：NODEUXUUZUSUM1 0.00000.00000.00000.00002 0.00000.00000.00000.000030.44621-1.70830.00001.76564 -0.35379-1.35450.00001.3999按照书上例题所计算的公式,将载荷、长度、弹性模量以及面积带入得到的 结果为：节点1和2的位移为0节点3的位移为(0.448,-1.712) cm,节点4 的位移为(-0.352,-1.352) cm,此结果与例题计算结果相符，说明了该计算是 正确的。选择 General Postp

13、rocPlot Results Countour Plot Nodal Solu，在弹出的 对话框中选择 Nodal SolutionDOF Solution Displacement vector sum，即出现 桁架杆件的总位移矢量图，最大位移出现在施加载荷的附近，为22.09cm,如下:3.2 应力和应变分析选择 General PostprocElement TableDefine Table,出现如下对话框，点击 “Add,出现如下对话框，在 User label for item 中输入 Axial stress-i，在 Item, Comp Results data item 中

14、选择 By sequence num，选择“LS,并在其后加 1,表 示节点i。同样的添加轴向应变项目，即在User label for item中输入Axial strain-i， 在 Item, Comp Results data item 中选择 By sequence num，选择 “LEPEL,并在其 后加1，表示节点i点击“OK”，如下图所示：（注:sequence num由帮助文件得到,LINK180单元的Sequence numbers如下图：）Table 180.2 LINK180 Item and Sequenee NumbersOutput Quantity NameET

15、ABLE and ESOL Command InputItemEIJSxxLS-12EPELxxLEPEL-12EFTOXXLEFTO-12EFTHXXLEPTH-12EPPLXXLEPPL-12EPCRXXLEPCR-12FORCESMISC1AREASMISC2TEMPLBFE-12选择 General PostprocElement TableList Elem Table，出现对话框:选择AXISSTRE和AXIALSTR，点击“OK”则列出了每个单元节点i的轴向应力和轴向应变，见下图：ELEMAKISSTREAXIALSTR10.00000.0000211155.0.55776E-B

16、33-8844.8-0.44224E-034-8844.8-0.44224E-B3512509.0.62543E-036-15776.-0.78879E-B3对于桁架杆单元，其应力在整个杆内的分布是完全相同的，即节点j和节点i 相同，对于应变也是如此。所以单元1的应力为0，单元2的应力为11155N/cm2, 即111.55Mpa,单元3为88.45Mpa，单元4的应力为88.45Mpa,单元5的应力 为125.09Mpa，单元6的应力为157.76Mpa。Q235的抗拉强度为190Mpa，故每 根杆的强度均满足要求。对单元 1 至 6，其轴向应变分别为 0， 0.558e-3， -0.442e-3， -0.442e-3， 0.625e-3， 0.789e-3。4.关于结构设计的改进意见及分析由计算结果可见，杆1没有受力，杆3和4的受力较小，从节约材料的角度 考虑，可以适当减小杆1、杆2以及杆3的截面尺寸，这不但不影响结构的可靠 程度，并尽可能的做到了等强度。