厚壁圆筒有限元分析报告
有限元与CAE分析报告
专业：
班级：
学号：
姓名：
指导教师：
实习时间：
年月日
平面问题的厚壁圆筒问题
一、问题提出
如图所示为一厚壁圆筒，其内半径为r1=50mm,外半径为r2=100mm,作用在内孔上的压力p=10 Mpa,无轴向压力，轴向长度很长可视为无穷，要求对其进行结构静力分析，并计算厚壁圆筒径向应力和切向应力沿半径r方向的分布。弹性模量E=200 Gpa,泊松比μ=0.3。
图1 厚壁圆筒
二、建模步骤
1 定义工作文件名
依次单击Utility Menu>File>Change Jobname，在文本框中输入：1245523229，在“New Log and error files”处选中“yes”，单击“OK”。
2 定义工作标题
依次单击Utility Menu>File>Change Title ，在文本框中输入：1245523229，单击“OK”。依次单击Plot>Replot,
3 定义单元类型
1)依次单击Main Menu>Prefrences，选中“Structural”，单击“OK”。
2)依次单击Main Menu>Preprocessor>Element type>Add/Edit/Delete，出现对话框，单击“Add”，出现一个“Library of Element Type”对话框，。在“Library of Element Type”左面的列表栏中选择“Structural Solid”，在右面的列表栏中选择“Quard 4node 182”，单击“OK”。
2) 单击对话框中的“Options”，在弹出的单元属性对话框中，选择K3关键字element behavior为“Plane strain”,再单击“Close”，完成单元的设置。
3 定义实常数
依次单击Main Menu>Preprocessor>RealConstants>Add/Edit/Delete,出现对话框define element constant sets，点击“set1”，单击“Add”；出现对话框choose element type，点击“Type1 Plane182”，出现Real Constant Set Number 1 for Plane 182对话框，在“Thickness”中输入1，单击“OK”，再单击“Close”，完成厚度实常数的设置。
4 设置材料属性
依次单击Main Menu>Preprocessor>MaterialProps>Material Modes，出现“Define Material ModelBehavior”对话框，在“Material Model Available”下面的对话框中，双击打开“Structural>Linear>Elastic>Isotropic”,出现对话框，输入弹性模量EX＝2e5，泊松比PRXY＝0.3，单击“OK”，单击“Material>Exit”。
5 建立几何模型
1) 创建矩形面
拾取菜单Main Menu>Preprocessor>Molding>Create>Areas>Rectangle>By Dimensions。弹出如下图所示的对话框，在“X1,X2”文本框中输入0,40，在“Y1,Y2”文本框中输入0,40，单击OK。
2) 创建圆形面
拾取菜单Main Menu>Preprocessor>Molding>Create>Areas>Circle>Solid Circle。弹出如下图所示的对话框，在“WPX”、“WPY”、“Radius”文本框中输入0,0,10，单击OK。
3) 布尔运算，形成带孔方板
拾取菜单Main Menu>Preprocessor>Molding>Operate>Booleans>Subtract>Areas。弹出拾取窗口，拾取矩形面，单击OK，再次弹出拾取窗口，拾取圆形面，单击OK。
4)合并相同几何特征，压缩编号
拾取菜单Main Menu>Preprocessor>Numbering Controls，选择Merge Items，在弹出菜单中，[NUMMRG]区域“lab”下拉列表中选择all，单击OK；拾取菜单Main Menu>Preprocessor>Compress Numbers，在弹出菜单中，[NUMCMP]区域“lab”下拉列表中选择all，单击OK。
6 划分单元
1）拾取菜单Main Menu>Mesh>Mesh Tool。弹出对话框，点击size controls部分中中lines的set，点击图形窗口中的圆弧边，弹出如下图所示的对话框，在NDIV 选项中输入10，单击OK。按照同样的方法选择其他各边，在NDIV选项中输入5。
2）在Mesh tool中，将Mesh部分中的shape选择为“Quard”、“Mapped”，在下拉菜单中，选择“Pick corners”，在弹出菜单中首先选中整个面，点击OK，然后在弹出菜单中按照如下图的顺序选择四个关键点，再次点击OK。划分出所示的规则映射网格。
7 施加载荷和边界条件
1) 拾取菜单Utility Menu>Select>Entities 。弹出select entities 对话框，选择实体类型为Node ，选择创建选择集的方法为“By Location ”，选中“X coordinates ”，在“Min ，Max ”文本框中输入0，选中“From Full ”，单击OK 按钮。
2) 拾取菜单Main Menu>Solution>Define Loads>Apply>Structural>Displacement>On nodes，弹出拾取窗口，拾取“Pick All”，弹出Apply U,ROT on nodes对话框，在列表中选择“UX”，单击OK。
3) 按照同样方式，拾取菜单Utility Menu>Select>Entities。选择实体类型为Node，选择创建选择集的方法为“By Location”，选中“Y coordinates”，在“Min，Max”文本框中输入0，选中“From Full”，单击OK按钮。
4) 拾取菜单Main Menu>Solution>Define Loads>Apply>Structural>Displacement>On nodes，弹出拾取窗口，拾取“Pick All”，在列表中选择“UY”，单击OK。
5) 拾取菜单Utility Menu>Select>Everything。
6) 拾取菜单Main Menu>Solution>Define Loads>Apply>Structural>Pressure>On lines，弹出拾取窗口，选择2号线，弹出如下图所示的对话框，在“Value”文本框中输入-10，点击OK按钮。
8 求解模型
依次点击Main Menu>Solution>Solve>Current LS进行求解，当提示Solution is Done说明求解完成。
三、仿真结果分析
依次单击Main Menu>GeneralPostproc>Plot Results>Contour Plot>Nodal Solu，弹出对话框，在“Item to be Contoured”栏中选择“Stress”，在其右面的栏中，选择“Sy”后，单击“OK”，这时在图形窗口处显示出平面的应力云图，如图所示。
尝试将圆弧边划分为20、40、80、160份，其他直边对应划分为10、20、40、80份，查看并记录y方向的应力，绘制出图标，得出随着网格加密，应力集中系数趋于稳定的结论。
四、学习《有限元与CAE技术》的心得体会
签名：
日期
有限元分析实验报告
武汉理工大学 学生实验报告书 实验课程名称机械中的有限单元分析 开课学院机电工程学院 指导老师姓名 学生姓名 学生专业班级机电研 1502班 2015—2016 学年第2学期
实验一方形截面悬臂梁的弯曲的应力与变形分析 钢制方形悬臂梁左端固联在墙壁，另一端悬空。工作时对梁右端施加垂直向下的30KN的载荷与60kN的载荷，分析两种集中力作用下该悬臂梁的应力与应变，其中梁的尺寸为10mmX10mmX100mm的方形梁。 1.1方形截面悬臂梁模型建立 建模环境：DesignModeler 15.0。 定义计算类型：选择为结构分析。 定义材料属性：弹性模量为2.1Gpa，泊松比为0.3。 建立悬臂式连接环模型。 （1）绘制方形截面草图：在DesignModeler中定义XY平面为视图平面，并正视改平面，点击sketching下的矩形图标，在视图中绘制10mmX10mm的矩形。（2）拉伸：沿着Z方向将上一步得到的矩阵拉伸100mm，即可得到梁的三维模型，建模完毕，模型如下图1.1所示。 图1.1 方形截面梁模型 1.2 定义单元类型： 选用6面体20节点186号结构单元。 网格划分：通过选定边界和整体结构，在边界单元划分数量不变的情况下，通过分别改变节点数和载荷大小，对同一结构进行分析，划分网格如下图1.2所示：
图1.2 网格划分 1.21 定义边界条件并求解 本次实验中，讲梁的左端固定，将载荷施加在右端，施以垂直向下的集中力，集中力的大小为30kN观察变形情况，再将力改为50kN，观察变形情况，给出应力应变云图，并分析。 （1）给左端施加固定约束； （2）给悬臂梁右端施加垂直向下的集中力； 1.22定义边界条件如图1.3所示： 图1.3 定义边界条件 1.23 应力分布如下图1.4所示： 定义完边界条件之后进行求解。
有限元分析实验报告
学生学号1049721501301实验课成绩 武汉理工大学 学生实验报告书 实验课程名称机械中的有限单元分析机电工程学院开课学院 指导老师姓名
学生姓名 学生专业班级机电研1502班 学年第学期2016—20152 实验一方形截面悬臂梁的弯曲的应力与变形分析 钢制方形悬臂梁左端固联在墙壁，另一端悬空。工作时对梁右端施加垂直 向下的30KN的载荷与60kN的载荷，分析两种集中力作用下该悬臂梁的应力与应变，其中梁的尺寸为10mmX10mmX100mm的方形梁。 方形截面悬臂梁模型建立1.1 建模环境：DesignModeler15.0。 定义计算类型：选择为结构分析。 定义材料属性：弹性模量为 2.1Gpa，泊松比为0.3。 建立悬臂式连接环模型。 （1）绘制方形截面草图：在DesignModeler中定义XY平面为视图平面，并正 视改平面，点击sketching下的矩形图标，在视图中绘制10mmX10mm的矩形。 （2）拉伸：沿着Z方向将上一步得到的矩阵拉伸100mm，即可得到梁的三维模型，建模完毕，模型如下图 1.1所示。
图1.1方形截面梁模型 ：定义单元类型1.2 选用6面体20节点186号结构单元。 网格划分：通过选定边界和整体结构，在边界单元划分数量不变的情况下，通过分别改变节点数和载荷大小，对同一结构进行分析，划分网格如下图 1.2
所示： 图1.2网格划分 1.21定义边界条件并求解 本次实验中，讲梁的左端固定，将载荷施加在右端，施以垂直向下的集中 力，集中力的大小为30kN观察变形情况，再将力改为50kN，观察变形情况，给出应力应变云图，并分析。 （1）给左端施加固定约束； （2）给悬臂梁右端施加垂直向下的集中力； 1.22定义边界条件如图1.3所示：
第三节 厚壁圆筒应力分析
第三节厚壁圆筒应力分析 3.3厚壁圆筒应力分析 3.3.1弹性应力 3.3.2弹塑性应力 3.3.3屈服压力和爆破压力 3.3.4提高屈服承载能力的措施 3.3.1弹性应力 i i c c o o 本小节的目的：求弹性区和塑性区里的应力 假设：a.理想弹塑性材料 b.圆筒体只取远离边缘区 图2-2
1、塑性区应力 平衡方程： r r d r dr θσσσ-= （2-26） M i s e s 屈服失效判据 ：r s θσσ-= （2-40） 联立积分，得 ln r s r A σ= + （2-41） :i r i r R p σ==-内壁边界条件，求出A 后带回上式得 ln r s i i r p R σ= - （2-42） 将（2-42）带入（2-40）得 1ln s i i r p R θσ?? = +- ??? （2-43） 12ln 2 r z i i r p R θ σσσ? ? += = +-?? （2-44） 将:c r c r R p σ==-代入（2-42）得 ln c c s i i R p p R =+ （2-45） 结论: ①(,//)i i s f R r p σσ= ②,(ln ) r r f r r θθσσσ=↑↑，， ③1()2z r const θσσσ=+≠（区别: 弹区1 ()2 z r const θσσσ=+=） 弹性区内壁处于屈服状态: ()( )Kc=Ro/Rc c c r s r R r R θσσ==-= 由表2-1拉美公式得出 ：22 c p = （2-46） 与2-45联立导出弹性区与塑性区交界面的p i 与R c 的关系 2202ln )c c i i R R p R R =-+ （2-47）
有限元分析报告
《有限元基础理论》报告 学院： 班级： 姓名： 学号： 任课老师： 二〇一一年十二月
题目一：三维托架实体受力分析 题目：1、三维托架实体受力分析：托架顶面承受50psi的均匀分布载荷。托架通过有孔的表面固定在墙上，托架是钢制的，弹性模量E=29×106psi，泊松比v=0.3.试通过ANSYS输出其变形图及其托架的von Mises应力分布。 题目1的分析：先进行建模，此建模的难点在对V3的构建（既图中的红色部分）。要想构建V3，首先应将A15做出来，然后执行Main Menu>Preprocessor>Modeling>Operate>Booleans>Add>V olumes命令，将所有的实体合并为一个整体。建模后，就对模型进行网格的划分，实行Main Menu>Preprocessor>Meshing>MeshTool,先对网格尺寸进行编辑，选0.1，然后点Meshing，Pick all进行网格划分，所得结果如图1.1。划分网格后，就可以对模型施加约束并进行加载求解了。施加约束时要注意，由于三维托架只是通过两个孔进行固定，故施加约束应该只是针对两孔的内表面，执行Main Menu>Solution>Define Loads>Apply>Structrual>Displacement>Symmetry B.C>On Areas命令，然后拾取两孔的内表面，单击OK就行了。施加约束后，就可以对实体进行加载求解了，载荷是施加在三维托架的最顶上的表面的，加载后求解运算，托架的变形图如图1.2。
图1.1、托架网格图 图1.2输出的是原型托架和施加载荷后托架变形图的对比，虚线部分即为托架的原型，从图1.2可看出，由于载荷的作用，托架上面板明显变形了，变形最严重的就是红色部分，这是因为其离托板就远，没有任何物体与其分担载荷，故其较容易变形甚至折断。这是我们 在应用托架的时候应当注意的。
有限元上机实验报告
有限元上机实验报告结构数值分析与程序设计 上机实验 院系: 土木工程与力学学院专业: 土木工程 班级: 姓名: 学号: 指导教师: 1、调试教材P26-30程序FEM1。 1.1、输入数据文件为: 6,4,12,6,1.0E0,0.0,1.0,0.0,1 3,1,2 5,2,4 3,2,5 6,3,5 0.0,2.0 0.0,1.0 1.0,1.0 0.0,0.0 1.0,0.0 2.0,0.0
1,3,7,8,10,12 1.2、输出数据文件为: NN NE ND NFIX E ANU T GM NTYPE 6 4 12 60.1000E+01 0.000 1.0000.0000E+00 1 NODE X-LOAD Y-LOAD 1 0.000000E+00 -0.100000E+01 2 0.000000E+00 0.000000E+00 3 0.000000E+00 0.000000E+00 4 0.000000E+00 0.000000E+00 5 0.000000E+00 0.000000E+00 6 0.000000E+00 0.000000E+00 NODE X-DISP Y-DISP 1 -0.879121E-15 -0.325275E+01 2 0.879121E-16 -0.125275E+01 3 -0.879121E-01 -0.373626E+00 4 0.117216E-1 5 -0.835165E-15 5 0.175824E+00 -0.293040E-15 6 0.175824E+00 0.263736E-15 ELEMENT X-STR Y-STR XY-STR 1 -0.879121E-01 -0.200000E+01 0.439560E+00 2 0.175824E+00 -0.125275E+01 0.256410E-15 3 -0.879121E-01 -0.373626E+00 0.307692E+00 4 0.000000E+00 -0.373626E+00 -0.131868E+00 2、修改FEM1,计算P31例2-2。
有限元分析报告样本
《有限元分析》报告基本要求： 1. 以个人为单位完成有限元分析计算，并将计算结果上交；（不允许出现相同的分析模型，如相 同两人均为不及格） 2. 以个人为单位撰写计算分析报告； 3. 按下列模板格式完成分析报告； 4. 计算结果要求提交电子版，报告要求提交电子版和纸质版。（以上文字在报告中可删除） 《有限元分析》报告 一、问题描述 （要求：应结合图对问题进行详细描述，同时应清楚阐述所研究问题的受力状况和约束情况。图应清楚、明晰，且有必要的尺寸数据。） 一个平面刚架右端固定，在左端施加一个y 方向的-3000N 的力P1，中间施加一个Y 方向的-1000N 的力P2，试以静力来分析，求解各接点的位移。已知组成刚架的各梁除梁长外，其余的几何特性相同。 横截面积：A=0.0072 m2 横截高度：H=0.42m 惯性矩：I=0.0021028m4ｘ 弹性模量： E=2.06ｘ10n/ m2/ 泊松比：u=0.3 二、数学模型 （要求：针对问题描述给出相应的数学模型，应包含示意图，示意图中应有必要的尺寸数据；如进行了简化等处理，此处还应给出文字说明。） （此图仅为例题）
三、有限元建模（具体步骤以自己实际分析过程为主，需截图操作过程） 用ANSYS 分析平面刚架 1．设定分析模块 选择菜单路径：MainMenu—preference 弹出“PRreferences for GUI Filtering”对话框，如图示，在对话框中选取：Structural”,单击[OK]按钮，完成选择。 2．选择单元类型并定义单元的实常数 （1）新建单元类型并定 （2）定义单元的实常数在”Real Constants for BEAM3”对话框的AREA中输入“0。0072”在IZZ 中输入“0。0002108”，在HEIGHT中输入“0.42”。其他的3个常数不定义。单击[OK]按 钮，完成选择 3．定义材料属性 在”Define Material Model Behavier”对话框的”Material Models Available”中，依次双击“Structural→Linear→Elastic→Isotropic”如图
机械零件有限元分析——实验报告
中南林业科技大学机械零件有限元分析 实验报告 专业：机械设计制造及其自动化 年级： 2013级 班级：机械一班 姓名：杨政 学号：20131461 I
一、实验目的 通过实验了解和掌握机械零件有限元分析的基本步骤；掌握在ANSYS 系统环境下，有限元模型的几何建模、单元属性的设置、有限元网格的划分、约束与载荷的施加、问题的求解、后处理及各种察看分析结果的方法。体会有限元分析方法的强大功能及其在机械设计领域中的作用。 二、实验内容 实验内容分为两个部分：一个是受内压作用的球体的有限元建模与分析，可从中学习如何处理轴对称问题的有限元求解；第二个是轴承座的实体建模、网格划分、加载、求解及后处理的综合练习，可以较全面地锻炼利用有限元分析软件对机械零件进行分析的能力。
实验一、受内压作用的球体的有限元建模与分析 对一承受均匀内压的空心球体进行线性静力学分析，球体承受的内压为 1.0×108Pa ，空 心球体的内径为 0.3m ，外径为 0.5m ，空心球体材料的属性：弹性模量 2.1×1011，泊松比 0.3。 承受内压：1.0×108 Pa 受均匀内压的球体计算分析模型（截面图） 1、进入 ANSYS →change the working directory into yours →input jobname: Sphere 2、选择单元类型 ANSYS Main Menu : Preprocessor →Element Type →Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK (back to Element Types window)→ Options… →select K3: Axisymmetric →OK →Close (the Element Type window) 3、定义材料参数 ANSYS Main Menu : Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3→ OK 4、生成几何模型生成特征点 ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标：input ：1(0.3，0)，2(0.5，0)，3(0，0.5)，4(0，0.3)→OK 生成球体截面 ANSYS 命令菜单栏: Work Plane>Change Active CS to>Global Spherical ANSYS Main Menu: Preprocessor →Modeling →Create →Lines →In ActiveCoord → 依次连接 1,2,3,4 点生成 4 条线→OK Preprocessor →Modeling →Create →Areas →Arbitrary →By Lines →依次拾取四条线→OK ANSYS 命令菜单栏: Work Plane>Change Active CS to>Global Cartesian 5、网格划分 ANSYS Main Menu : Preprocessor →Meshing →Mesh Tool →(Size Controls) lines: Set
有限元分析报告
有限元仿真分析实验 一、实验目的 通过刚性球与薄板的碰撞仿真实验，学习有限元方法的基本思想与建模仿真的实现过程，并以此实践相关有限元软件的使用方法。本实验使用HyperMesh 软件进行建模、网格划分和建立约束及载荷条件，然后使用LS-DYNA软件进行求解计算和结果后处理，计算出钢球与金属板相撞时的运动和受力情况，并对结果进行可视化。 二、实验软件 HyperMesh、LS-DYNA 三、实验基本原理 本实验模拟刚性球撞击薄板的运动和受力情况。仿真分析主要可分为数据前处理、求解计算和结果后处理三个过程。前处理阶段任务包括：建立分析结构的几何模型，划分网格、建立计算模型，确定并施加边界条件。 四、实验步骤 1、按照点－线－面的顺序创建球和板的几何模型 （1）建立球的模型：在坐标(0,0,0)建立临时节点，以临时节点为圆心，画半径为5mm的球体。 （2）建立板的模型：在tool-translate面板下node选择临时节点，选择Y-axis,magnitude输入5.5，然后点击translate+，return;再在2D－planes-square 面板上选择Y-axis,B选择上一步移下来的那个节点，surface only ,size=30。 2、画网格
（1）画球的网格：以球模型为当前part，在2D－atuomesh面板下，surfs 选择前面建好的球面，element size设为0.5mm，mesh type选择quads，选择elems to current comp,first order，interactive。 （2）画板的网格：做法和设置同上。 3、对球和板赋材料和截面属性 （1）给球赋材料属性：在materials面板内选择20号刚体，设置Rho为2.000e-08,E为200000,NU为0.30。 （2）给球赋截面属性：属性选择SectShll,thickness设置为0.1，QR设为0。 （3）给板赋材料属性：材料选择MATL1，其他参数：Rho为2.000e-08,E 为100000,Nu为0.30，选择Do Not Export。 （4）给板赋截面属性：截面选择SectShll，thickness设为0.2。其他参数:SHRE为8.333-01，QR为0，T1为0.2。 （5）给板设置沙漏控制：在Properties-Create面板下Card image选择HourGlass，IHQ为4，QM为0.100。更新平板。 4、加载边界条件 （1）将板上最外面的四行节点分别建成4个set。 （2）建立一个load collector。 （3）Analysis－constraints面板中，设置SIZE为1，nodes通过by sets 选择set_1、set_2、set_3、set_4，然后点击creat即可，边界条件加载完毕。 5、建立载荷条件（给球一个3mm的位移） （1）建立一个plot: post－xy plots-plots-creat plot,然后点击return;
有限元实验报告模板
有限元实验报告 T1013-5 20100130508 蔡孟迪
ANSYS有限元上机报告（一） 班级:T1013-5 学号:20100130508 姓名:蔡孟迪 上机题目： 图示折板上端固定，右侧受力F=1000N，该力均匀分布在边缘各节点上；板厚t=2mm 材料选用低碳钢，弹性模量E=210Gpa，μ=0.33. 一、有限元分析的目的： 1．利用ANSYS构造实体模型； 2．根据结构的特点及所受载荷的情况，确定所用单元类型；正确剖分网格并施加外界条件；3．绘制结构的应力和变形图，给出最大应力和变形的位置及大小；并确定折板角点Ａ处的应力和位移； 4．研究网格密度对Ａ处角点应力的影响； 5．若在Ａ处可用过渡圆角，研究Ａ处圆角半径对Ａ处角点应力的影响。 二、有限元模型的特点： 1．结构类型 本结构属于平面应力类型 2.单位制选择 本作业选择N（牛），mm（毫米），MPa（兆帕）。 3.建模方法 采用自左向右的实体建模方法。 4.定义单元属性及类型 1）材料属性：弹性模量：EX=2.10E5MPa, 泊松比：PRXY=0.33 2）单元类型：在Preferences选Structural，Preprocessor>ElemmentType>Add/Edit/Delete中定义单元类型为：Quad4 node 182，K3设置为：平面薄板问题（Plane strs w/thk） 3）实常数：薄板的厚度THK=2mm 5.划分网格 在MeshTool下选set,然后设置SIZE Element edge length的值，再用Mesh进行网格划分。6.加载和约束过程：在薄板的最上端施加X、Y方向的固定铰链，在薄板的最右端施加1000N 的均匀布置的载荷。
ANSYS实体建模有限元分析-课程设计报告
南京理工大学 课程设计说明书(论文) 作者:学号： 学院(系):理学院 专业:工程力学 题目:ANSYS实体建模有限元分析 指导者： (姓名) (专业技术职务) 评阅者： (姓名) (专业技术职务) 20 年月日
练习题一 要求： 照图利用ANSYS软件建立实体模型和有限元离散模型，说明所用单元种类、单元总数和节点数。 操作步骤： 拟采用自底向上建模方式建模。 1.定义工作文件名和工作标题 1)选择Utility Menu>File>Change Jobname命令，出现Change Jobname对话框，在[/FILNAM ] Enter new jobname文本框中输入工作文件名learning1，单击OK按钮关闭该对话框。 2)选择Utility Menu>File>Change Title命令，出现Change Title对话框，在[/TITLE] Enter new title文本框中输入08dp,单击OK按钮关闭该对话框。 2.定义单元类型 1)选择Main Menu>Preprocessor>Element Type>Add/Edit/Delete命令，出现Element Types对话框，单击Add按钮,出现 Library of Element Types 对话框。在Library of Element Types 列表框中选择 Structural Solid, Tet 10node 92，在Element type reference number文本框中输入1，单击OK按钮关闭该对话框。 2)单击Element Types对话框上的Close按钮，关闭该对话框。 3.创建几何模型 1)选择Utility Menu>P1otCtrls>Style>Colors>Reverse Video命令，设置显示颜色。 2)选择Utility Menu>P1otCtrls>View Settings>Viewing Direction命令，出现Viewing Direction对话框，在XV,YV,ZV Coords of view point文本框中分别输入1, 1, 1，其余选项采用默认设置，单击OK按钮关闭该对话框。 3)建立支座底块 选择Main Menu>Preprocessor> Modeling>Create>volumes>Block>By Demensios 命令，出现Create Block by Demensios对话框，在X1,X2 X-coor dinates文本框
厚壁圆筒有限元分析报告
有限元与CAE分析报告 专业： 班级： 学号： 姓名： 指导教师： 实习时间： 年月日
平面问题的厚壁圆筒问题 一、问题提出 如图所示为一厚壁圆筒，其内半径为r1=50mm,外半径为r2=100mm,作用在内孔上的压力p=10 Mpa,无轴向压力，轴向长度很长可视为无穷，要求对其进行结构静力分析，并计算厚壁圆筒径向应力和切向应力沿半径r方向的分布。弹性模量E=200 Gpa,泊松比μ=0.3。 图1 厚壁圆筒 二、建模步骤 1 定义工作文件名 依次单击Utility Menu>File>Change Jobname，在文本框中输入：1245523229，在“New Log and error files”处选中“yes”，单击“OK”。 2 定义工作标题 依次单击Utility Menu>File>Change Title ，在文本框中输入：1245523229，单击“OK”。依次单击Plot>Replot, 3 定义单元类型 1)依次单击Main Menu>Prefrences，选中“Structural”，单击“OK”。
2)依次单击Main Menu>Preprocessor>Element type>Add/Edit/Delete，出现对话框，单击“Add”，出现一个“Library of Element Type”对话框，。在“Library of Element Type”左面的列表栏中选择“Structural Solid”，在右面的列表栏中选择“Quard 4node 182”，单击“OK”。 2) 单击对话框中的“Options”，在弹出的单元属性对话框中，选择K3关键字element behavior为“Plane strain”,再单击“Close”，完成单元的设置。
有限元实验报告
一、实验目的 通过上机对有限元法的基本原理和方法有一个更加直观、深入的理解；通过对本实验所用软件平台Ansys 的初步涉及，为将来在设计和研究中利用该类大型通用CAD/CAE 软件进行工程分析奠定初步基础。 二、实验设备 机械工程软件工具包Ansys 三、实验内容及要求 1) 简支梁如图3.1.1所示，截面为矩形，高度h=200mm ，长度L=1000mm ，厚 度t=10mm 。上边承受均布载荷，集度q=1N/mm2，材料的E=206GPa ，μ=0.29。平面应力模型。 X 方向正应力的弹性力学理论解如下： 图3.1.1 ①在Ansys 软件中用有限元法探索整个梁上x σ，y σ的分布规律。 ②计算下边中点正应力x σ的最大值；对单元网格逐步加密，把x σ的计算值与理论解对比，考察有限元解的收敛性。 ③针对上述力学模型，对比三节点三角形平面单元和4节点四边形平面等参元的求解精度。 2) 一个正方形板，边长L = 1000mm ，中心有一小孔，半径R = 100mm ，左右边 受均布拉伸载荷，面力集度q = 25MPa ，如图 3.2.1所示。材料是 206E GPa =，0.3μ=，为平面应力模型。当边长L 为无限大时，x = 0截面上理论解为： ) 534()4 (6222 23-+-=h y h y q y x L h q x σ
)32(2|44 220r R r R q x x ++==σ 其中R 为圆孔半径，r 为截面上一点距圆心的距离。x = 0截面上孔边（R r =）应力q x 3=σ。所以理论应力集中系数为3.0。 图3.2.1 用四边形单元分析x = 0截面上应力的分布规律和最大值，计算孔边应力集中系数，并与理论解对比。利用对称性条件，取板的四分之一进行有限元建模。 3) 如图3.3.1所示，一个外径为0.5m ，内径为0.2m ，高度为0.4m 的圆筒，圆 筒的外壁施加100MPa 的压强，圆筒的内部约束全部的自由度，材料参数是密度。 使用平面单元，依照轴对称的原理建模分析。 q
厚壁圆筒应力分析
厚壁圆筒应力分析 1、概述 K>1.2的壳体成为厚壁圆筒。厚壁容器承压的应力特点有（此处不考虑热应力）：一、不能忽略径向应力，应做三向应力分析；二、厚壁容器的应力在厚度方向不是均匀分布，而是应力梯度。所以，在求解的时候需要联立几何方程、物理方程、平衡方程才能确定厚壁各点的应力大小。 2、解析解 一、内压为i p ，外压为0p 的厚壁圆筒，需要求出径向应力r σ、周向应力θσ和轴向应力z σ，其中轴向应力z σ不随半径r 变化。 （1）几何方程 如图所示，取内半径r ，增量为dr 的一段区域两条弧边的径向位移为ω和ωωd +，其应变的表达式为： r rd rd d r dr d dr d r ω θθθωεω ωωωεθ= -+== -+= ))(（周向应力：径向应力：（1） θσ对r 求导，得： ()θθσσωωωω ωσ-=??? ??-=-='??? ??=r r r dr d r r r dr d r dr d 112 （2） （2）物理方程 根据胡克定理表示为
[]z E σσμσεθθ+-= r (1 （3） 两式相减，消去z σ得： []θθσσμεε-+= r E ） （1-r []z r E σσμσεθ+-=(1r （4） 将（4）代入（2）得： []） z r E dr d σσμσεθθ+-=(1 （5） 对（3）的θε求导得，z σ看做常数： ?? ? ??-=dr r d dr d E dr d σμσεθθ1 （6） 联立（5）、（6）得： []θθθσσμσμσ-) 1-r r dr d dr d += （ （7） （3）平衡方程 如图所示，沿径向和垂直径向建立坐标 系，把θσ向x 轴和y 轴分解，得： ?? ? ??=-+2sin 2θθd p p p r dr r （8） 其中 ()θσσd dr r d p r r dr r ++=+)（ （9） θσrd p r r = 由于θd 很小，2 2sin θ θd d ≈??? ??，略去二阶微量r r d d σ，得 dr d r r r σσσθ=- （10） 联立（7）（10）得 0322=+dr d dr d r r r σσ （11）
有限元分析报告
创新实习报告 题目名称基于Solidworks simulation的潜孔冲击器前接头有限元分析学院（系）机械工程学院 专业班级材料成型及控制工程0801班 学生姓名（10） 指导教师杨雄教授 日期2012.2.27 至2012.3.23
基于Solidworks simulation的潜孔冲击器前接头有限元分析 目录 1.有限元分析软件简介 (2) 2.潜孔冲击器前接头实物及断口相片 (5) 3.潜孔冲击器前接头的基本属性，工作情况，受力情况的分析 (6) 4.利用三维画图软件建模 (7) 5. 利用solidworkd sinulation对零件进行有限元分析 (14) 5.1 分析原理及步骤…………………………………………………………… 5.2 算例属性…………………………………………………………………… 5.3 单位………………………………………………………………………… 5.4 材料属性…………………………………………………………………… 5.5 载荷和约束………………………………………………………………… 5.6 载荷………………………………………………………………………… 5.7 接触………………………………………………………………………… 5.8 网格信息…………………………………………………………………… 5.9 反作用力，自由实体力，自由体力矩…………………………………… 5.10 算例结果………………………………………………………………… 6.分析结论 (15) 6.1失效分析…………………………………………………………………… 6.2提出优化方案………………………………………………………………… 6.3对优化方案进行有限元分析………………………………………………… 6.4分析比较并得出结论………………………………………………………… 7.小结 (18) 8.参考文献 (18)
ansys实验报告
有限元上机实验报告 姓名柏小娜 学号0901510401
实验一 一 已知条件 简支梁如图所示，截面为矩形，高度h=200mm ，长度L=1000mm ，厚度t=10mm 。上边承受均布载荷，集度q=1N/mm 2，材料的E=206GPa ，μ=0.29。平面应力模型。 X 方向正应力的弹性力学理论解如下： )534()4 (6222 23-+-=h y h y q y x L h q x σ 二 实验目的和要求 （1）在Ansys 软件中用有限元法探索整个梁上x σ，y σ的分布规律。 （2）计算下边中点正应力x σ的最大值；对单元网格逐步加密，把x σ的计算值与理论解对比，考察有限元解的收敛性。 （3）针对上述力学模型，对比三节点三角形平面单元和4节点四边形平面等参元的求解精度。 三 实验过程概述 （1） 定义文件名 （2） 根据要求建立模型：建立长度为1m ，外径为0.2m ，平行四边行区域 （3） 设置单元类型、属性及厚度，选择材料属性： （4） 离散几何模型，进行网格划分 （5） 施加位移约束 （6） 施加载荷 （7） 提交计算求解及后处理 （8） 分析结果 四 实验内容分析 （1）根据计算得到应力云图，分析本简支梁模型应力分布情况和规律。主要考察x σ和y σ，并分析有限元解与理论解的差异。 由图1看出沿X 方向的应力呈带状分布，大小由中间向上下底面递增，上下底面应力方向相反。由图2看出应力大小是由两侧向中间递增的，得到X 方向
上最大应力就在下部中点，为0.1868 MPa 。根据理论公式求的的最大应力值为0.1895MPa 。由结果可知，有限元解与理论值非常接近。由图3看出Y 的方向应力基本相等，应力主要分布在两侧节点处。 图 1 以矩形单元为有限元模型时计算得出的X 方向应力云图 图 2 以矩形单元为有限元模型时计算得出的底线上各点x 方向应力图 （2）对照理论解，对最大应力点的x σ应力收敛过程进行分析。列出各次计算 应力及其误差的表格，绘制误差－计算次数曲线，并进行分析说明。 答：在下边中点位置最大应力理论值为： MPa h y h y q y x L h q x 1895.0)5 34()4(622223=-+-=σ
有限元分析报告
班级：土木1204 学号：19 姓名：廖枭冰
班级：土木1204 学号：23 姓名：梅雨辰
混凝土上承式空腹式拱桥研究 一引言 本文通过SAP2000软件，对混凝土上承式空腹式拱桥在上部车辆荷载作用下，各个部位的内力和应力的分布进行分析，对强度和刚度进行校核，提出存在的问题，最后进行改进。 工程实例图 模型三维图
二 模型尺寸及构件截面 该拱桥总跨度L=80m,高H=20m,宽度10m ，分为五个构件 1拱肋：一段圆弧线，水平投影长度80m,采用箱型截面，高1.6m,宽2m ，翼缘厚度0.22m,腹板厚度0.15m 2主梁：长80m,采用箱型截面，高6m,宽2m,翼缘厚度1.1m,腹板厚度0.55m 3立柱：拱桥与主梁的之间的竖向构件，采用矩形截面，长宽均为1.2m ，分别在桥的每隔10m 布置1根 4横系梁：拱肋之间的横向构件，采用矩形截面，高0.6m,宽0.4m 5桥面：长80m,宽10m,厚度为0.6m,保护层厚度30mm 三 材料定义 所有构件均采用C50混凝土，配置钢筋，抗压强度,50cu k f MPa =,弹性模量43.4510E MPa =? 四 计算模型假设与简化 ⑴由于拱肋，主梁，立柱，横系梁长度远大于宽度及高度，将其定义为杆件单元。 ⑵由于桥面的厚度远小于其长度和宽度，将其定义为平面厚壳单元。 ⑶圆弧拱肋采用在圆弧线上取点，用折线杆件进行逼近。 ⑷由于拱肋伸入桥台或桥墩，位移和转角均被束缚，两端采用固定端
约束，形成无铰拱模型。 ⑸由于主梁支撑在刚度较其大的多的桥台或桥墩上，又考虑到主梁长度方向的热胀冷缩，将其一端定义为固定铰支座，另一端定义为辊轴支座。 ⑹由于工程实际多采用混凝土现浇工艺，所有构件的连接处视为刚接 ⑺由于拱顶与主梁之间的混凝土的厚度较小，可忽略这部分混凝土，让拱顶与主梁直接接触。 ⑻由于桥面的重量较其它杆件大得多，故只考虑桥面的重量。 ⑼计算车辆对桥面的荷载时，不考虑车辆的具体尺寸，将其定义为均布荷载加在桥面上。 五模型受力分析 在桥面上施加规范规定的2 kN m的公路一级荷载，来模拟车辆对 10.5/ 桥的压力。 六结果展示（分析与校核） 1 强度分析 桥面单元 桥 面 弯 矩
有限元编程大作业报告
本科生实验报告书 四节点等参单元有限元分析的FORTRAN程序
目录 1．问题概述 (1) 2．四节点四边形等参单元介绍 (1) 3．单元应力磨平方法介绍 (4) 4．程序流程设计 (6) 程序设计概述 流程图 5．程序结构及程序说明 (8) 6．程序应用及算例分析 (9) 算例概述 算例ANSYS求解 算例程序数值解 算例分析 7. 总结 (15)
1. 问题概述 等参单元是有限元方法中使用最广泛的单元类型。等参单元的位移模式和坐标变换均采用相同的形函数，这种坐标变换叫做等参变换。通过等参变换可以将自然（局部）坐标中几何形状规则的单元转换成总体（笛卡尔）坐标中形状扭曲的单元，因而使得单元有较好的适应性。 本问题首先对平面四节点四边形等参单元的形函数、应力矩阵和等效节点力矩阵、应力磨平公式等的推导和计算求解。并通过设计FORTRAN 求解程序进行编程求解，最后给出算例（受集中荷载的悬臂梁）并进行求解，将解与ANSYS 的解进行比较。在这个过程中，采用了高斯三点积分和高斯两点积分，这种积分方法的求解效率较高而且精度也较好。在问题的最后，尝试去分析引起数值解误差的原因，并分析四节点等参单元的若干特性。 2. 四节点四边形等参单元介绍 边长为2的正方形单元（如下图所示），在其形心处安置一个局部坐标。单元角 结点i 的坐标（,）分别为 ,因此单元四条边界的方程可用简单公式 和 逐一给出。 图2-1 母单元 1 2 3 4 0
图2-2 四边形单元 形函数的表达式： 位移函数： ∑==41 i N u i i u ∑==4 1 i N v i i v 坐标变换式： ∑==41 i x N x i i ，∑==4 1 i y N y i i 单元应变矩阵 {}[]{}[]{}e e B B B B B x v y u y v x u δδε4321==??? ? ??? ?????????????+??????= 式中{ }[]T T T T T 4 321e δδδδδ=——单元节点的位移列阵；
有限元上机实验报告
有限元法基础及应用 上机报告 南京理工大学 2015年12月 上机实验一
1 实验题目 设计一个采用减缩积分线性四边形等参元的有限元模型，通过数值试验来研究网格密度、位移约束条件与总刚度矩阵奇异性、沙漏扩展、求解精度的关系，并验证采用减缩积分时保证总刚度矩阵非奇异的必要条件。总结出你的研究结论，撰写实验报告。 2 实验目的 通过实验来研究减缩积分方案中网格密度和位移约束条件对总体刚度矩阵奇异性和求解精度的影响，以此加深对有限元减缩积分的理解，和对减缩积分中保证总体刚度矩阵非奇异性的认识。 3建模概述 先保持位移约束条件不变，研究网格密度对总体刚度矩阵奇异性和求解精度的影响，并验证采用减缩积分时保证总刚度矩阵非奇异的必要条件。如下图1所示，建立一个简支和链杆的约束条件，然后不断增加网格密度，通过ABAQUS 来计算位移和应力的变化规律。 个独立关系式）节点（两个自由度）
4 计算结果分析讨论与结论 1）1*1单元四边形减缩积分实验 载荷布种/单元 应力云图 2）2*1单元四边形减缩积分实验 载荷单元
应力云图3）4*4单元四边形减缩积分实验 载荷布种单元 应力云图
结果分析 5 实验体会与小结 单元刚度矩阵的特征： (1)对称性 (2)奇异性 (3)主元恒正 K相同 (4)平面图形相似、弹性矩阵D、厚度t相同的单元，e K的分块子矩阵按结点号排列，每一子矩阵代表一个结点，占两行两 (5)e 列，其位置与结点位置对应。 整体刚度矩阵的特征： (1)对称性 (2)奇异性 (3)主元恒正 (4)稀疏性 (5)非零元素呈带状分布。 [K]的物理意义是任意给定结构的结点位移所得到的结构结点力总体上满足力和力矩的平衡。为消除[K]的奇异性，需要引入边界条件，至少需给出能限制刚体位移的约束条件。 对于一个给定形式的单元，如果采用精确积分，则插值函数中所有项次在|J|=常数的条件下能被精确积分，并能保证刚度矩阵的非奇异性。如果采用减缩积分，因为插值函数中只有完全多项式的项次能被精确积分，因此需要进行刚度矩阵非奇异必要条件的检查。
有限元分析报告格式
有限元分析及应用研究报告宋体三号字 ----------题目自定×××××宋体四号字 作者姓名宋体五号字 扬州大学机械工程学院××专业江苏，扬州 225009宋体五号字 摘要：宋体五号字 摘要的写法：可以是陈述式的，也可以是信息式的，或者二者兼而有之。 作为一般的学术论文，通常采用信息式的摘要，其内容主要包括：研究课题的目的、研 究方法、所获结果及结论； 关键词：宋体五号字 [1] [2] [3] 1、问题的引出（引言、前言）正文宋体小四号字，标题加粗，多倍行距1.2 内容 应含有：是一个什么问题？为什么要做有限元分析？目的和作用？等等？在工程上的重要 性？要解决什么问题？（研究目标、研究内容和拟解决的关键问题） 2、理论研究内容研 究方法：语句（二次开发）；程序流程图；…… 3、有限元模型确定有限元计算模型（所 分析问题的数学建模） (1)取对称结构（1/4结构） (2)引人支承条件 (3)载荷移置 (4)单 元的选取、单元、节点的个数 4、结果及讨论（数据处理、图表、对比结果、误差分析） 5、 结论 6、参考文献要规范（按正规杂志要求）参考文献××× 备注： (1) 可以模拟书上的列子。 (2) 时间2015年底完成。 (3) 交打印稿给班长，收齐送s403.篇二：有限元分析报告格式 标准有限元分析报告格式(请勿转载) ***************************************************************************** 如 果你的计算结果想通过en/iso/suv等认让的话，就要按下面的格式写 **************************************************************************** 3. specification for the layout of the safety proof: a. index of content b. page numbering c. cross references, if used f. the description of the classification system’s information has to be precisely. 4. content of the safety proof, issued by the author (in this case by us): a. short introduction: i. task which has to be solved ii. should be illustrated by a sketch iii. short description of the methods which are used to solve the problem (theories, calculation methods...) iv. followed regulations v. terms, formula symbols, units b. programme information: i. name ii. version and version names