有限元上机报告一
T1343-2 周坤诚 20130420216
作业三要求
1、属于力学的那类问题？
答：属于薄板弯曲问题。
2、单位制；单元类型；单元描述；实常数；材料参数
答：长度单位为m，力的单位为N，应力的单位为Pa。
选择shell壳单元中shell63单元，为4节点四边形单元。实常数是板的厚度为0.1m。Shell63 弹性壳。具有弯矩和薄膜特性。可承受与平面同方向及法线方向的荷载。每个节点6个自由度：x,y,z方向和绕x,y,z轴方向。有应力强化和大变形能力。提供用于大变形分析的连续性相切矩阵。材料为低碳钢，弹性模量为210GPa。
3、划分网格的方案，施加载荷和边界条件
答：自由网格剖分。四边加载固定约束，中间加载均布载荷。
4、计算结果的体现：MISIS应力和位移，需要知道最大值以及位置。
5、1/4模型的计算结果。
6、总体结论感想。
在实验指导书中的步骤的指导下，完成了上机操作的内容，虽然不是很理解，但是那种自己分析出来的结果很是满意，通过上机操作，我对有限元有了一定的认识，原本一个问题，用材料力学或理论力学计算好长时间才得到的结果，到现在只需要简单的几步操作就可以的出来了，感觉十分神奇，方便。并且在自己完成报告的过程中对有限元分析更加透彻，感觉有限元分析作用真的很大，总之，通过在理论课上的学习及上机操作再加上老师布置的报告这一学习过程，让我受益匪浅，更加完善了我对这一课程及软件的理解。
有限元上机实验指导书
弹性力学及有限元实验上机指导书 土木工程教学部 20１5年6月 一、ANSYS软件安装及有限元建模方法演示１、实验目的 掌握Ansys商用有限元软件的安装及了解有限元建模方法２、实验任务 １）Ansys商用软件8.1安装过程详解。
２）采用有限元直接建模法创建杆系模型演示。 3）采用有限元间接建模法创实体模型演示。 模型一：平面桁架如下图所示，长度单位为m，求支座反力和各杆内力。弹性模量2.1e+11，泊松比0.3，杆件截面面积均为0.01m2。 1/6 1/3 1/2 图1 平面桁架模型 模型二：正方形带孔平板，边长1m，小孔直径0.1m，板厚0.05m。弹性模量2.1e+11，泊松比0.3。上下边受均匀压力1000N。 图2 带孔正方形平板 3、实验方法 实验方法同课堂操作演示。 命令流(模型一) /PREP7 ET,1,LINK1 R,1,0.01, , MP,EX,1,,2.1e+11
MP,PRXY,1,,0.3 n,1,0,0 n,2,1,0 n,3,1,0.5 n,4,2,0 n,5,2,0.833 n,6,3,0 n,7,3,1 n,8,4,0 n,9,4,0.833 n,10,5,0 n,11,5,0.5 n,12,6,0 e,1,2 e,1,3 e,2,3 e,2,4 e,3,5 e,4,5 e,3,4 e,4,6 e,5,7 e,6,7 e,4,7 e,6,8 e,7,9 e,8,9 e,7,8 e,8,10
e,9,11 e,10,11 e,8,11 e,10,12 e,11,12 save 命令流(模型二) 二、利用ANSYS创建杆系或实体结构有限元模型 １、实验目的 掌握有限元建模的基本方法并能创建简单的杆系结构和实体结构有限元模型。 ２、实验任务 （1）采用直接建模法创建上节平面桁架结构模型。 （2）采用间接建模法创建上节带孔平板实体模型。 3、实验方法 实验方法同操作演示。 三、有限元求解及结果后处理演示 １、实验目的 掌握基本参数的设置、荷载施加方法及后处理操作。 ２、实验任务 （1）读入数据文件（命令流）的形式生成杆系结构有限元模型 （2）实常数、材料参数、求解参数设置演示 （3）位移约束、集中荷载施加方法演示 （4）计算求解与后处理操作演示。 3、实验方法 实验方法同课堂操作演示 附：后处理GUI操作及命令流操作 A、通过后处理提取节点计算结果 三种后处理操作： 1、plot results
有限元分析实验报告
武汉理工大学 学生实验报告书 实验课程名称机械中的有限单元分析 开课学院机电工程学院 指导老师姓名 学生姓名 学生专业班级机电研 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所示： 定义完边界条件之后进行求解。
有限元分析程序设计
结构有限元分析程序设计 绪论 §0.1 开设“有限元程序设计”课程的意义和目的 §0.2 课程特点 §0.3 课程安排 §0.4 课程要求 §0.5 基本方法复习 $0.1 意义和目的 1.有限元数值分析技术本身要求工程设计研究人员掌握 1). 有限元数值分析技术的完善标志着现代计算力学的真正成熟和实用化,已在各种 力学中得到了广泛的应用。比如：，已杨为工程结构分析中最得以收敛的技术手段，现代功用大致有： a). 现代结构论证。对结构设计从内力，位移等方面进行优劣评定，从而进 行结构优化设计。 b)可取代部份实验，局部实验+有限元分析，是现代工程设计研究方法的一大 特点。 c)结构的各种功能分析（疲劳断裂，可靠性分析等）都以有限元分析工具作为 核心的计算工具。 2). 有限元数值分析本身包括着理论+技术实现（本身功用所绝定的） 有限元数值分析本身包括着泛函理论+分片插值函数+程序设计 2. 有限元分析的技术实现（近十佘年的事）更依赖于计算机程序设计 有限元分析的技术取得的巨大的成就，从某种意义上说，得益于计算机硬件技术的发展和程序设计技术的发展，这两者的依赖性在当代表现得更加突出。（如可视化技术） 3.从学习的角度，不仅要学习理论，而且要从程序设计设计角度对这些理论的技术实现有 一个深入的了解，应当致力于掌握这些技术实现能力，从而开发它，发展它。（理论本身还有待于进一步完美相应的程序设计必须去开发） 4.程序设计不仅是实现有限元数值分析的工具和桥梁，而且在以下诸方面也有意义： 1). 精通基本概念，深化理论认识； 2). 锻炼实际工程分析，实际动手的能力； 3). 获得以后工作中必备的工具。（作业+老师给元素库） 目的：通过讲述有限元程序设计的技术与技巧，便能达到自编自读的能力。 §0.2 课程特点 总描述：理论+算法+数据结构（程序设计的意义） 理论：有限元算法，构造，步骤，解的等外性，收敛性，稳定性，误差分析 算法；指求解过程的技术方法，含两方面的含义；a. 有限元数值分析算法，b, 与数据结构有关的算法（总刚稀疏存贮，提取，节点优化编号等） 数据结构：指各向量矩阵存贮管理与实现，辅助管理结构（指针，数据记录等） 具体特点： 理论性强：能量泛函理论+有限元构造算法+数据结构构造算法 内容繁杂：理论方法+技术方法+技术技巧 技巧性强：排序，管理结构（指针生成，整型运算等）
有限元上机实验报告
有限元上机实验报告结构数值分析与程序设计 上机实验 院系: 土木工程与力学学院专业: 土木工程 班级: 姓名: 学号: 指导教师: 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。
大作业报告参考2有限元学习心得
有限元学习心得 吴清鸽车辆工程 50110802411 短短八周的有限元课已经结束。关于有限元，我一直停留在一个很模糊的概念。我知道这是一个各个领域都必须涉及的点，只要有关于CAE分析的，几乎都要涉及有限元。总体来说，这是一门非常重要又有点难度的课程。 有限元方法(finite element method) 或有限元分析(finite element analysis)，是 求取复杂微分方程近似解的一种非常有效的工具，是现代数字化科技的一种重要 基础性原理。将它用于在科学研究中，可成为探究物质客观规律的先进手段。将 它应用于工程技术中，可成为工程设计和分析的可靠工具。本课程教学基本内容 有固体力学和结构力学简介；有限元法基础；桁架、梁、刚架、二维固体、板和 壳、三维固体的有限元法；建模技术；热传导问题的有限元分析；PATRAN软件 的使用. 通过有限元分析课程学习使我了解和掌握了一些有限元知识： 1．简要了解二维和三维固体以及桁架、梁和板结构的三组基本力学方程，即表示位移-应变关系的几何方程，表示应力-应变关系的本构方程和表示内力-外力关系的平衡方程。 2．了解利用能量法形成有限元离散系统方程的基本原理，即哈密尔顿原理。掌握有限元分 析的基本方法及步骤，包括域的离散、位移插值、构造形函数、单元有限元方程 的建立、坐标变换、整体有限元方程的组装、整体有限元方程的求解技术。 3．具体深入的了解并掌握桁架结构、梁结构、刚架结构、二维固体、板和壳结构、三维固体的有限元法分析技术，包括他们具体的形函数构造，应变矩阵，局部坐标系和整体坐标系中的单元矩阵。各种结构的实例研究。 4．了解并掌握建立高质量建模所涉及的各种关键技术。包括单元类型的选择，单元畸形的限制，不同阶数单元混用时网格的协调性问题，对称性的应用（平面对称、轴对称、旋转对称、重复对称），由多点约束方程形成刚域及应用（模拟偏移、不同自由度单元的连接、网格协调性的施加）等，以及多点约束方程的求解。以PATRAN有限元通用软件为例了解一般商业有限元软件的组成及结构。掌握PATRAN软件的基本使用。利用PATRAN软件上机实践完成两个上机练习：刚架结构有限元分析和三维固体有限元分析。 课程的具体学习内容： 内容： 1.三节点三角形单元：单元分析、总刚度矩阵组装、引入约束条件修正总刚度 矩阵、载荷移置、方程求解； 2.四边形单元分析、四节点四面体单元分析、八节点六面体单元分析；
中南大学ANSYS上机实验报告
ANSYS上机实验报告 小组成员：郝梦迪、赵云、刘俊 一、实验目的和要求 本课程上机练习的目的是培养学生利用有限单元法的商业软件进行数值计算分析，重点是了解和熟悉ANSYS的操作界面和步骤，初步掌握利用ANSYS建立有限元模型，学习ANSYS分析实际工程问题的方法，并进行简单点后处理分析，识别和判断有限元分析结果的可靠性和准确性。 二、实验设备和软件 台式计算机，ANSYS10.0软件 三、基本步骤 1）建立实际工程问题的计算模型。实际的工程问题往往很复杂，需要采用适当的模型在计算精度和计算规模之间取得平衡。常用的建模方法包括：利用几何、载荷的对称性简化模型，建立等效模型。 2）选择适当的分析单元，确定材料参数。侧重考虑一下几个方面：是否多物理耦合问题，是否存在大变形，是否需要网格重划分。 3）前处理（Preprocessing）。前处理的主要工作内容如下：建立几何模型（Geometric Modeling），单元划分（Meshing）与网格控制，给定约束（Constraint）和载荷（Load）。在多数有限元软件中，不能指定参数的物理单位。用户在建模时，要确定力、长度、质量及派生量的物理单位。在建立有限元模型时，最好使用统一的物理单位，这样做不容易弄错计算结果的物理单位。建议选用kg，N，m，sec；常采用kg，N，mm，sec。 4）求解（Solution）。选择求解方法，设定相应的计算参数，如计算步长、迭代次数等。 5）后处理（Postprocessing）。后处理的目的在于确定计算模型是否合理、计算结果是否合理、提取计算结果。可视化方法（等值线、等值面、色块图）显
有限元作业第二次作业
土木工程专业 有限元第二次作业 姓名： 班级： 学号： 指导教师： 二〇一五年6月12日
习 题：平面应力问题的八节点等参元，已给定8个节点 的坐标。试查资料并论述： 1、单元中位移函数u （ξ，η），v （ξ，η）和单元节点位 移{δe }的关系式； 2、[ B ]矩阵的计算步骤和计算式； 3、单元刚度矩阵[ k e ]的一般计算方法和计算步骤； 4、论述相邻单元间公共边界上位移的连续性； 5、如果给定母单元中点A ， （ξ，η），怎样求实际单元中与 A ， 相对应的点A （x ，y ）；反之，如果给定实际单元中的点A （x ，y ），怎样求其在母单元中对应点A ， （ξ，η）？ 6、如果已经求解得到单元8个节点的位移值{δe }怎样求单 元中某一点B （x ，y ）的应力？ 实际单元 1 2 6 7 Y 1 2 43 67 8η= 1η=﹣1 母单元 ξ= 1 ξ=﹣1
解： 1、此题分两步进行： 单元位移场的表达： 如图1所示,在任意四边形的每边中间设一附加节点，则单元边界就变成二次曲线的了。如果直接在整体坐标系(),x y 下，像八节点矩形元那样，构造双二次多项式的位移插值函数，则因曲边四边形单元边界是二次曲线，故边界上的位移是()x y 或的五次多项式， 它不能由曲边上三个节点的位移分量唯一地决定，从而不能保证相邻两个单元在公共边上位移的协调条件，所以在整体坐标系(),x y 下构造完全协调的位移插值 函数是很困难的，利用坐标变换，可将曲边四边形单元变换成基本单元，如图2所示的在自然坐标(),ξη下具有边长为2的八节点正方形单元，自然坐标系(),ξη是外节点坐标值为±1的局部坐标系。在自然坐标系的单元上构造 协调的位移插值函数，其形状函数是较普通的，取位移分量为,ξη的双二次多项式, 即： 2222 12345678222 2910111213141516u a a a a a a a a v a a a a a a a a ξηξξηηξηξηξηξξηηξηξη?=+++++++??=+++++++?? （1-1） 利用8 个节点的16 个位移分量可唯一确定16 个待定常数1216,,a a a …,，若代入8个节点的局部坐标值，得： 图1：在总坐标系中具有二 次曲边的四边形单元 图2：在自然坐标系中的 曲边四边形的基本单元
机械零件有限元分析——实验报告
中南林业科技大学机械零件有限元分析 实验报告 专业：机械设计制造及其自动化 年级： 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
有限元 第二次作业
2-2 图示悬臂板,属于平面应力问题,其网格图及单元、节点编号见图2-1,E=2、1×1011,u=0、28,演算其单刚阵到总刚阵得组集过程,并用MATLAB 软件计算总刚阵。 图2-1 答:根据图2-1所示列出单元节点列表: i j k 1 3 5 4 2 2 5 3 3 2 6 5 4 1 6 2 (1)计算单元刚度阵 单元1得刚度矩阵: ,; 单元2得刚度矩阵:,; 单元3得刚度矩阵:,; 单元4得刚度矩阵:,; 总刚度矩阵: []??????? ?? ? ????? ?????++++++++++++=4 6,636,635 ,642 ,632,641 ,636 ,535 ,525,515,514,523 ,513,53 2,522,515,414,413,425 ,315,314,323 ,313,322 ,346,236,235 ,225,223 ,242 ,232,222,241,246,142 ,141 ,10 00000 00 000k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k K 节点 单元
Matlab 程序语言得编写: function Idex global gNode gElement gMaterial gNode=[0、0 0、01 0、5 0、01 1、0 0、01 1、0 0、0 0、5 0、0 0、0 0、0] %gNode 同样就是一个矩阵,每一行表示一个结点,第1 列就是结点得x 坐标,第2 列就是结点得y坐标 gElement=[3 4 5 2 3 5 2 5 6 1 2 6 ]; %gElement 就是一个矩阵,每一行表示一个单元,第1 行就是单元得第1 个结点号,第2 行就是单元得第2个结点号。 Return function k=StiffnessMatrix(ie) %计算单元刚度矩阵函数 global gNode gElement k=zeros(6,6); %6x6单元刚阵 E=2、1*10^11; %材料特性 u=0、28 ; %材料特性 t=0、01; %材料特性 xi=gNode(gElement(ie,1),1);
有限元实验报告模板
有限元实验报告 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 的均匀布置的载荷。
有限元08上机作业
调试书本26到30页程序 开列数组维数： DIMENSION LOC(4,3),CX(6),CY(6),IFIX(6),F(12), 1GK(12,12),STRES(4,3),BAK(4,3,6) 结点集中力输入： DO 10 I=1,ND 10 F(I)=0.0 F(2)=-1.0 数据文件输入： 6,4,12,6,1.0E0,0.0,1.0,0.0,1 3,1,2 5,2,4 2,5,3 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 数据文件输出： NN NE ND NFIX E ANU T GM NTYPE 6 4 12 6 0.1000E+01 0.000 1.000 0.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
有限元课程实训结课上机参考例题
有限元课程实训结课上机参考例题 题 1：图1所示薄板左边固定，右边受均布压力P=100Kn/m 作用，板厚度为0.3cm ；试采用如下方案，对其进行有限元分析，并对结果进行比较。 (1)三节点三角形单元；（2个和200个单元） (2)四节点矩形单元； （1个和50个单元） (3)八节点等参单元。 （1个和20个单元） 图 1 题 2：图2所示为一带圆孔的单位厚度（1M ）的正方形平板，在x 方向作用均布压力0.25Mpa ，试用三节点常应变单元和六节点三角形单元对平板进行有限元分析，并对以下几种计算方案的计算结果进行比较： （1） 分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算； （2）分别采用不同数量的三节点常应变单元计算； 注：在y 轴上，孔边应力的精确解为：MPa x 75.0-=σ，在x 轴上，孔边应力的精确解为：MPa y 25.0=σ 图 2
题 3：图3所示为带方孔（边长为80mm）的悬臂梁，其上受部分均布载荷（p=10Kn/m）作用，试采用一种平面单元，对图示两种结构进行有限元分析，并就方孔的布置进行分析比较，如将方孔设计为圆孔，结果有何变化？（板厚为1mm，材料为钢） 图 3 题 4：图4所示为一隧道断面，其内受均布水压力q，外受土壤均布压力p；试采用不同单元计算断面内的位移及应力，并分别分析q=0或p=0时的位移和应力分布情况。(材料为钢，隧道几何尺寸和压力大小自行确定) 图 4
题 5：图5所示一简化直齿轮轮齿截面，高h=60mm，齿根宽b=60mm，齿顶宽c=25mm，齿顶作用力P=10Kn；试采用不同单元分析轮齿上位移及应力分布，并只指出最大应力位置。 图 5 题 6：图6所示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： (1)分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算； (2)分别采用不同数量的三节点常应变单元计算； (3)当选常应变三角单元时，分别采用不同划分方案计算。 图 6
有限元实验报告
一、实验目的 通过上机对有限元法的基本原理和方法有一个更加直观、深入的理解；通过对本实验所用软件平台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
ansys上机作业
实验一坝体的有限元建模及应力应变分析 一、实验目的： 1、掌握ANSYS软件基本的几何形体构造方法、网格划分方法、边界条件施加方法及各 种载荷施加方法。 2、熟悉有限元建模、求解及结果分析步骤和方法。 3、能利用ANSYS软件对结构进行有限元分析。 二、实验设备： 微机，ANSYS软件。 三、实验内容： 计算分析模型如图所示，分析坝体的应力、应变。 四、实验步骤： 1 进入ANSYS 程序→ANSYS →change the working directory into yours →input Initial jobname: dam 2设置计算类型 ANSYS Main Menu: Preferences →select Structural → OK 3选择单元类型
ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 42 →OK (back to Element Types window)→Options… →select K3: Plane Strain →OK→Close (the Element Type window) 4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear→Elastic→Isotropic→input EX:2.1e11, PRXY:0.3→ OK 5生成几何模型 生成特征点 ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标：input:1(0,0),2(1,0),3(1,5),4(0.45,5)→OK 生成坝体截面 ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS→依次连接四个特征点，1(0,0),2(1,0),3(1,5),4(0.45,5) →OK 生成坝体截面如图一 图一 6 网格划分 ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条横边:OK→input NDIV: 15 →Apply→依次拾取两条纵边:OK →input NDIV: 40 →OK →(back to the mesh tool window)Mesh: Areas, Shape: Quad, Mapped →
有限元上机题
注：题中E表示材料弹性模量，μ表示泊松比，ρ表示密度。 一静力结构分析 1 如图1所示为普通订书钉，E=2.1×105MPa，μ=0.3，横截面尺寸为宽B=0.64mm，高H=0.402mm。当订书钉被压入纸张时，约需要120N的载荷，载荷均匀地分布在订书钉上部。以下面两种情况进行有限元分析。（单位：mm）（1）订入时A、B点为铰支条件； （2）订入时A、B点为固定约束。 图1 载荷和尺寸情况 2、小型铁路桥由横截面积为3250mm2的钢制杆件组装而成。一辆火车停在桥上，其载荷施加在桥梁两侧的桁架上，单侧的桁架如图2所示，等效载荷为F1，F2，E=2.1×105MPa，μ=0.3，ρ=7.8×103kg/m3。试计算位置R处由于载荷作用而沿水平方向移动的距离以及支反力，同时，分析各个节点的位移和非单元应力。 图2 铁路桥单侧桁架及载荷情况 3如图3所示，模型参数为：E=3.0×1010Pa，A1=30m，A2=10m，B=80m，t=20m，p=2200Pa。 有关风载的确定，按照海洋井架行业标准，有以下方法： 风压（Pa）=0.6115×风速（m/s）×高度系数×形状系数对于一般的海洋井架及建筑物，高度在30m左右，高度系数取为1.1，形状系数取为1.25，风速取为47.8m/s。换算出来后得到的风压为2200Pa。
图3 高层建筑物及其风载荷 4对于含裂纹体的结构及材料，若按照线弹性力学分析，会在裂纹的尖端处产生应力的奇异性，这时需要计算裂纹尖端处的应力强度因子（对于Ⅰ型裂纹，有K1=σ（πa）1/2），并以应力强度因子作为准则来对材料的裂纹是否扩展进行判断。图4所示为一块矩形平板，其边缘存在长为a的裂纹，板的两端承受拉应作用。利用结构上下的对成性，取矩形的一半建立有限元模型，完成看一力σ 下工作： （1）球裂纹的张角θ（在施加载荷前=0，θ=0） （2）沿直线AO，画出y方向应力σy沿x变化的曲线图。假设 σy= K1/(2πx)1/2利用回归方法估计K1。将计算结果与计算无限大平板的修正结果进行比较，其中基于无限大平板的修正关系式为K1=1.2σ0（πa）1/2）（3）在裂缝尖端处，进行网格的细化，重新求解（2）中的问题。 模型中的相关参数为 E=2.1×105MPa，μ=0.3，L=400mm，a=9.5mm，b=95mm，σ0=450MPa。 图4 承受拉升载荷的带裂纹平面结构
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=-+-=σ
ansys 上机实验指导书
实验一ANSYS软件环境 一、实验目的： 熟悉ANSYS软件菜单、窗口等环境、软件分析功能及解题步骤。 二、实验设备： 微机，ANSYS软件。 三、实验内容： ANSYS软件功能、菜单、窗口及解题步骤介绍。 四、实验步骤： 1、ANSYS界面介绍： ANSYS软件功能非常强大，应用范围很广，并具有友好的图形用户界面（GUI）和优秀和程序架构。基于Motif标注的GUI主要由主窗口和输出窗口组成。随着版本的不断升级，ANSYS界面不断改进，不同版本间的界面存在着较大差别。下面介绍ANSYS的用户界面。 （1）主窗口 图1-1 ANSYS主窗口 ANSYS的主窗口主要由以下5个部分组成。 ①Utility菜单 这些菜单主要通过ANSYS的相关功能组件起作用，比如文件控制、参数选择、图像参数控制及参数输入等。 ②Input Line(Input Window命令输入窗口) 命令输入窗口（也称为命令栏）用于显示程序的提示信息并允许用户直接输入命令，简化分析过程。
③工具栏（Toolbar） 工具栏主要由按钮组成，这些按钮都是ANSYS中的常用命令。用户可以根据工作类型定义自己的工具栏以提高分析效率。 ④主菜单（Main Menu） 主菜单包括了ANSYS最主要的功能，分为前处理器（Preprocessor）、求解器（Solution）、通用后处理器（General Postprocessor）、设计优化器（Design Optimizer）。展开主菜单可以看到非常多的树状建模命令，这也是ANSYS7.0版本和以前版本的一个显著差别。虽然菜单的外观改变了，但是菜单结构没有变化，这对ANSYS用户平滑升级非常有利。 ⑤图形窗口（Graphic Windows） 图形窗口用于显示分析过程的图形，实现图形的选取。在这里可以看到实体建模各个过程的图形并可查看随后分析的结果。 （2）输出窗口（Output Windows） 输出窗口用于显示程序的文本信息，即以简单表格形式显示过程数据等信息。通常，输出窗口被主窗口遮盖，当然，如果需要随时可以将输出窗口拖到前面。 注意： 应该在ANSYS分析的各个步骤中随时查看输出窗口中的信息，检验分析过程是否正确，以便及时调整。 通过GUI可以方便地交互式访问程序的各种功能、命令、用户手册和参考材料，一步步地完成整个分析，很好地体现出ANSYS的易用性。同时，ANSYS软件也提供了完整的在线说明和帮助文件，以协助有经验的用户进行高级应用。 在用户界面中，ANSYS软件提供了4种通用的命令输入方式：菜单、对话框、工具栏和直接输入命令。 2、ANSYS分析过程： 一个典型的ANSYS分析过程包含3个主要步骤，每个主要步骤及其子步骤如下： （1）建立有限元模型 在ANSYS中建立有限元模型的过程大致可分为以下3个主要步骤： ①建立或导入几何模型 ②定义材料属性 ③划分网格建立有限元模型 （2）施加载荷并求解 在ANSYS中施加载荷及求解的过程大致可以分为以下3个主要步骤： ①定义约束 ②施加载荷 ③设置分析选项并求解 （3）查看分析结果 在ANSYS中查看分析结果的过程大致可以分为以下2个主要步骤： ①查看分析结果 ②检验分析结果（验证结果是否正确）
有限元上机实验报告
有限元法基础及应用 上机报告 南京理工大学 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|=常数的条件下能被精确积分，并能保证刚度矩阵的非奇异性。如果采用减缩积分，因为插值函数中只有完全多项式的项次能被精确积分，因此需要进行刚度矩阵非奇异必要条件的检查。
有限元第二次作业
2-2图示悬臂板，属于平面应力问题，其网格图及单元、节点编号见图2-1，E=×1011，u=，演算其单刚阵到总刚阵的组集过程，并用MATLAB软件计算总刚阵。 图2-1 答：根据图2-1所示列出单元节点列表： i j k 1354 2253 3265 4162 （1）计算单元刚度阵 单元1的刚度矩阵：[] ? ? ? ? ? ? ? ? ? ? = 1 5,5 1 4,5 1 3,5 1 5,4 1 4,4 1 3,4 1 5,3 1 4,3 1 3,3 1 k k k k k k k k k k，[] ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = 1 4,5 1 4,5 1 3,5 1 5,4 1 4,4 1 3,4 1 5,3 1 4,3 1 3,3 1 k k k k k k k k k k； 单元2的刚度矩阵：[] ? ? ? ? ? ? ? ? ? ? = 2 5,5 2 3,5 2 2,5 2 5,3 2 3,3 2 2,3 2 5,2 2 3,2 2 2,2 2 k k k k k k k k k k，[] ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? = 2 4,5 2 3,5 2 2,5 2 5,3 2 3,3 2 2,3 2 5,2 2 3,2 2 2,2 2 k k k k k k k k k k； 节点 单元
单元3的刚度矩阵：[] ??? ? ????? ?=36 ,635 ,632 ,636,535,532,536,235,232,23 k k k k k k k k k k ，[] ????????? ? ????? ?? ???=36,635 ,632,636,535,532 ,536,23 5 .23 2,2300000000 00 0000000000000000k k k k k k k k k k ； 单元4的刚度矩阵：[] ???? ? ???? ?=46,642 ,641,646,242,241 ,246,142,141 ,14 k k k k k k k k k k ，[] ?? ? ??? ? ?? ? ??????????=46,641 ,646,242,241.246,142 ,141,140000 000000000000000000000000k k k k k k k k k ； 总刚度矩阵：[][][][][][]4 3 2 1 4 1 k k k k k K e e +++=∑== []??????? ?? ? ????? ?????++++++++++++=4 6,636,635 ,642 ,632,641 ,636 ,535 ,525,515,514,523 ,513,53 2,522,515,414,413,425 ,315,314,323 ,313,322 ,346,236,235 ,225,223 ,242 ,232,222,241,246,142 ,141 ,10 00000 00 000k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k K Matlab 程序语言的编写： function Idex global gNode gElement gMaterial gNode=[