有限元实验内容_带孔板受拉分析实例
实验一ANSYS软件环境
一、实验目的：
熟悉ANSYS软件菜单、窗口等环境、软件分析功能及解题步骤。
二、实验设备：
微机（P4配置），ANSYS软件（教学版）。
三、实验内容：
ANSYS软件功能、菜单、窗口及解题步骤介绍。
四、实验步骤：
1、ANSYS界面介绍：
ANSYS软件功能非常强大，应用范围很广，并具有友好的图形用户界面（GUI）和优秀和程序架构。基于Motif标注的GUI主要由主窗口和输出窗口组成。随着版本的不断升级，ANSYS界面不断改进，不同版本间的界面存在着较大差别。下面介绍ANSYS7.0的用户界面。
（1）主窗口
ANSYS7.0的主窗口主要由以下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个主要步骤：
①查看分析结果
②检验分析结果（验证结果是否正确）
ANSYS 典型实例分析
如图所示，使用ANSYS分析平面带孔平板，分析在均布载荷作用下板内的应力分布。
已知条件：F＝20N/mm，L＝200mm，b＝100mm，圆孔半径r＝20，圆心坐标为（100，
50），E＝200Gpa。板的左端固定。
图1-2 带孔平板模型
1．问题描述
难度级别：普通级别。
所需时间：一个小时或者更多（视ANSYS操作熟练程度而定）。
实例类型：ANSYS结构分析。
分析类型：线性静力分析。
单元类型：PLANE82
ANSYS功能示例：实体建模包括基本的建模操作，布尔运算和网格细化；施加均布载荷；显示变形后形状和应力等值线图、单元信息列表；基本的结果验证技巧。
ANSYS帮助文件：在ANSYS Structural Analysis Guide了解Structural Static Analysis 分析知识，在ANSYS Elements Reference部分了解Plane82单元的详细资料。
2．建立有限元模型
1)．建立工作目录并添加标题
以Interactive 方式进入ANSYS，File菜单中设置工作文件名为Plane、
标题为plane。
2)．创建实体模型
（1）创建矩形
通过定义原点、板宽和板高定义矩形，其操作如下：
GUI：PreProcessor > Modeling > Create > Areas > Rectangle > By 2
Corners
弹出Rectangle by 2 corners对话框（如图1-3所示），如图填写。WP X 和
WP Y表示左下角点坐标。
命令：BLC4,0,0,200,100
（2）生成圆面
首先在矩形面上生成圆，然后挖去生成圆孔。
生成圆面得操作如下：
GUI ：PreProcessor > Modeling > Create > Areas > Circle > Solid
Circle
弹出Solid Circular Area 对话框（如图1-4 所示），依图输入圆面几何参
数。
命令：CYL4,100,50,20
下面通过布尔“减”操作生成圆孔，其操作如下：
GUI ：Processor > Modeling > Operate > Booleans > Subtract > Areas
先选择矩形面为Base Area ，单击OK 按钮，然后选择圆，单击OK 按钮。
布尔操作完毕之后，实体模型为带孔平板。
3)．定义材料属性
材料属性是与几何模型无关的本构关系，如弹性模量、密度等。虽然
材料属性不是与单元直接相联系在一起，但是由于计算单元矩阵时需要材
料属性，ANSYS 为了用户分析过程中定义材料属性方便，对每个单元类型
进行了相应的分类。根据不同类型的应用，材料属性可以是线性或非线性的。与单元类型相似，材料也可以定义多个，系统自动根据材料定义的顺序编号。本问题只有一种材料，因此只需定义一种材料，而且只需定义弹性模量和泊松比，
其操作如下：
GUI ：PreProcessor > Material Props > Material models > Structural > Linear > Elastic >Isotropic 在弹出对话框中键入EX=200000（单位Mpa ），PRXY ＝0.3。
4)．划分网格
划分网格首先选择合理的单元类型，然后定义单元的实常数，最后根据分析问题的需要
划分网格。
（1）选择单元
对于任何分析，必须在单元类型库中选择一种或者多种合适的单元类型。单元类型决定了附加的自由度（位移，转角、温度等）。许多单元还需要设置一些单元选项，比如单元特性和假设。单元结果的打印输出选项等，对于本问题选择Plane82单元。选择单元得操作如下：
GUI ：PreProcessor Menu > Element Type > Add/Edit/Delete
选择Plane82，弹出单元类型对话
框（如图1-5 所示）。单击OK
按钮。
命令：ET ，1，plane82
（2）定义单元实常数
有限单元的几何特性，不能
仅用其节点的位置充分表达，这
时需要提供一些实常数来补充
几何参数。典型的实常数有壳单
元的厚度，梁单元的横截面参数，板单元的厚度等。这些单元类型所需要的实常数以实数数值的形式输入。本问题所用单元类型为带厚度平面应力分析，因此分析类型设定为Plane strs w/thk 类型，操作如下：
图1-3 生成矩形
图1-4 生成圆
图1-5 选择单元
GUI：PreProcessor Menu > Element Type > Add/Edit/Delete > Options
在K3项后面下拉菜单中选“Plane strs w/thk“。
命令：KEYOPT,1,3,3
单元厚度为20mm，定义单元厚度操作如下：
GUI：PreProcessor Menu > Real Constants > Add/Edit/Delete > Add
在弹出的对话框中THK后文本框中键入材料厚度值20。
命令：R，1，20
（3）设定网格尺寸
这里让ANSYS 知道需要划分多大网格。采用用户自定义网格尺寸参数，其操作如下：
GUI：PreProcessor > Meshing > Size Cntrls > Manual Size > Areas > All Areas
在弹出Element Size on All Selected Areas对话框，在SIZE栏键入20mm。
命令：AESIZE,ALL,20
（4）划分网格
让ANSYS 知道网格大小后，现在划分网格，操作如下：
GUI：Processor > Meshing > Mesh > Areas > Free > Pick All
命令：AMESH,ALL
（5）保存数据库
GUI：Utility Menu > File > Save as...
输入文件名为Mesh （表示分析进度：已完成网格划分）。
3．施加载荷并求解
在这里首先定义模型约束，然后施加载荷，最后求解，为后处理查看结果提供数据，具
体操作步骤如下所示：
1)．定义约束
由已知得，需要固定（Fix）板左边线，即需要约束线上节点所有自由度（All DOFs），其操作如下：
GUI：Solution > Define Loads > Apply > Structural > Displacement > On Lines
弹出Apply U，ROT on Lines 对话框。选择板左侧边线，在Lab2 栏选All DOF。单击Apply按钮。2)．施加载荷
在板右侧边施加均布载荷，载荷大小为20/20＝1Mpa，施加载荷操作如下：
GUI：Solution > Define Loads > Apply > Structural > Pressure > On Lines
弹出Apply PRES on Lines 对话框，键入载荷值-1（由于载荷方向离开板，为拉力所以为负值，反之为正），单击OK按钮。
3)．求解
GUI：Solution > Solve > Current LS
4．查看分析结果
下面查看分析结果，对于静力分析主要是模型位移及等效应力等值线图或者节点结果数据列表。
1)．显示模型变形图
其操作如下：
GUI：General Postproc > Plot Results > Deformed Shape
2)．显示位移等值线分布图
其操作如下：
GUI：General Postproc > Plot Results > Nodal Solution> DoF Solution > Displacement Vector sum
由图获知最大位移值。
3)．显示等效应力等值线图
其操作如下：
GUI：General Postproc > Plot Results > Contour Plot > Nodal Solu > Stress > Von Mises 由图发现最大应力出现在孔的上下顶点，与解析解吻合。
4)．显示变形动画
通过动画显示可以查看模型在载荷作用下的内力变化过程，以图形方式显示分析计算结果。其操作如下：
Utility Menu 中 PlotCtrls > Animate > Deformed Results
5)．列表显示位移结果数据
其操作如下：
GUI：General Postproc > List Results > Nodal Solution
查看左侧边上的节点位移是否为零（由于边上所有节点已被固定，所以任意节点的位移均应该为零，这也是验证的一个方面）。
6)．列表显示节点应力值
验证节点应力值是否合理，其操作如下：
GUI：General Postproc > List Results > Stress > Principals SPRIN
列表显示分析结果与参考数据相吻合，表明ANSYS分析结果可靠。
7)．列表显示反作用力值
在任何方向上，反作用力总是等于此方向上的载荷总和。通过显示反作用力可以检验分析结果是否合理，其操作如下：
GUI：General Postproc > List Results > Reaction Solution
专业ABAQUS有限元建模经验笔记
基于ABAQUS的有限元分析和应用 第一章绪论 1.有限元分析包括下列步骤： 2.为了将试验数据转换为输入文件，分析者必须清楚在程序中所应用的和由实验人员提供的材料数据的应力和应变的度量。 3.ABAQUS建模需注意以下内容： 4.对于许多包含过程仿真的大变形问题和破坏分析，选择合适的网格描述是非常重要的，需要认识网格畸变的影响，在选择网格时必须牢牢记住不同类型网格描述的优点。 第二章ABAQUS基础 1.一个分析模型至少要包含如下的信息：离散化的几何形体、单元截面属性、材料数据、载荷和边界条件、分析类型和输出要求。 ①离散化的几何形体：模型中所有的单元和节点的集合称为网格。 ②载荷和边界条件： 2.功能模块： （1）Assembly（装配）：一个ABAQUS模型只能包含一个装配件。 （2）Interaction（相互作用）：相互作用与分析步有关，这意味着用户必须规定相互作用是在哪些分析步中起作用。 （3）Load（载荷）：载荷和边界条件与分析步有关，这意味着用户指定载荷和边界条件是在哪些分析步中起作用。 （4）Job（作业）：多个模型和运算可以同时被提交并进行监控。 3.量纲系统 ABAQUS没有固定的量纲系统，所有的输入数据必须指定一致性的量纲系统，常用的一致性量纲系统如下：
4.建模要点 （1）创建部件：设定新部件的大致尺寸的原则必须是与最终模型的最大尺寸同一量级。（2）用户应当总是以一定的时间间隔保存模型数据（例如，在每次切换功能模块时）。（3）定义装配： 在模型视区左下角的三向坐标系标出了观察模型的方位。在视区中的第2个三向坐标系标出了坐标原点和整体坐标系的方向（X，Y和Z轴）。 （4）设置分析过程： （5）在模型上施加边界条件和荷载： 用户必须指定载荷和边界条件是在哪个或哪些分析步中起作用。 所有指定在初始步中的力学边界条件必须赋值为零，该条件是在ABAQUS/CAE中自动强加的。 在许多情况下，需要的约束方向并不一定与整体坐标方向对齐，此时用户可定义一个局部坐标系以施加边界条件。 在ABAQUS中，术语载荷通常代表从初始状态开始引起结构响应发生变化的各种因素，包括：集中力、压力、非零边界条件、体力、温度（与材料热膨胀同时定义）。
有限元分析实验报告
武汉理工大学 学生实验报告书 实验课程名称机械中的有限单元分析 开课学院机电工程学院 指导老师姓名 学生姓名 学生专业班级机电研 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所示： 定义完边界条件之后进行求解。
matlab有限元分析实例
MATLAB: MATLAB是美国MathWorks公司出品的商业数学软件，用于数据分析、无线通信、深度学习、图像处理与计算机视觉、信号处理、量化金融与风险管理、机器人，控制系统等领域。 MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室），软件主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JAVA的支持。 MATLAB有限元分析与应用:
《MATLAB有限元分析与应用》是2004年4月清华大学出版社出版的图书，作者是卡坦，译者是韩来彬。 内容简介： 《MATLAB有限元分析与应用》特别强调对MATLAB的交互应用，书中的每个示例都以交互的方式求解，使读者很容易就能把MATLAB用于有限分析和应用。另外，《MATLAB有限元分析与应用》还提供了大量免费资源。 《MATLAB有限元分析与应用》采用当今在工程和工程教育方面非常流行的数学软件MATLAB来进行有限元的分析和应用。《MATLAB有限元分析与应用》由简单到复杂，循序渐进地介绍了各种有限元及其分析与应用方法。书中提供了大量取自机械工程、土木工程、航空航天工程和材料科学的示例和习题，具有很高的工程应用价值。
abaqus有限元分析过程
一、有限单元法的基本原理 有限单元法（The Finite Element Method）简称有限元（FEM），它是利用电子计算机进行的一种数值分析方法。它在工程技术领域中的应用十分广泛，几乎所有的弹塑性结构静力学和动力学问题都可用它求得满意的数值结果。 有限元方法的基本思路是：化整为零，积零为整。即应用有限元法求解任意连续体时，应把连续的求解区域分割成有限个单元，并在每个单元上指定有限个结点，假设一个简单的函数（称插值函数）近似地表示其位移分布规律，再利用弹塑性理论中的变分原理或其他方法，建立单元结点的力和位移之间的力学特性关系，得到一组以结点位移为未知量的代数方程组，从而求解结点的位移分量. 进而利用插值函数确定单元集合体上的场函数。由位移求出应变, 由应变求出应力 二、ABAQUS有限元分析过程 有限元分析过程可以分为以下几个阶段 1.建模阶段: 建模阶段是根据结构实际形状和实际工况条件建立有限元分析的计算模型――有限元模型，从而为有限元数值计算提供必要的输入数据。有限元建模的中心任务是结构离散，即划分网格。但是还是要处理许多与之相关的工作：如结构形式处理、集合模型建立、单元特性定义、单元质量检查、编号顺序以及模型边界条件的定义等。
2.计算阶段:计算阶段的任务是完成有限元方法有关的数值计算。 由于这一步运算量非常大，所以这部分工作由有限元分析软件控制并在计算机上自动完成 3.后处理阶段: 它的任务是对计算输出的结果惊醒必要的处理， 并按一定方式显示或打印出来，以便对结构性能的好坏或设计的合理性进行评估，并作为相应的改进或优化，这是惊醒结构有限元分析的目的所在。 下列的功能模块在ABAQUS/CAE操作整个过程中常常见到，这个表简明地描述了建立模型过程中要调用的每个功能模块。 “Part（部件） 用户在Part模块里生成单个部件，可以直接在ABAQUS/CAE环境下用图形工具生成部件的几何形状，也可以从其它的图形软件输入部件。 Property（特性） 截面（Section）的定义包括了部件特性或部件区域类信息，如区域的相关材料定义和横截面形状信息。在Property模块中，用户生成截面和材料定义，并把它们赋于(Assign)部件。 Assembly（装配件） 所生成的部件存在于自己的坐标系里，独立于模型中的其它部件。用户可使用Assembly模块生成部件的副本（instance），并且在整体坐标里把各部件的副本相互定位，从而生成一个装配件。 一个ABAQUS模型只包含一个装配件。
带孔平板拉伸作业
带孔平板有限元分析 本文采用有限元法，对带圆孔的矩形平板进行了弹塑性受力分析，分析了圆孔处的应力集中现象，为其设计和应用提供了参考依据。 1. 研究问题概述 本文研究带圆孔矩形平板在轴对称拉力作用下的平面应力问题。平板开孔的应力问题是弹塑性力学平面中的一个经典的问题,也是实际工程中常见的问题。平板长200mm ，宽50mm ，厚8mm ，具体几何参数及受力见图1。 图1 平板几何参数及受力 2.弹性力学方法解答 由弹性力学知识知，在距圆孔圆心()r ρρ>处的径向正应力、环向正应力、切应力分别为： 222222 1c o s 211322p r p r r ρσψρρρ?????? =-+-- ? ????????? 22221cos 21322p r p r ?σψρρ????=+-+ ? ???? ? 2222sin 21132p r r ρψψρ ττψρρ???? ==--+ ?????? ? 沿着y 轴，90ψ=。，环向正应力为： 242413122r r p ?σρρ?? =++ ???
max 3q ?σ=由上表可知： ()max = 3K q ψ σ=故应力集中因子: 可见孔边最大应力比无孔时提高了3倍，应力集中系数k=3，如图2所示。 图2 孔边应力集中 3.有限元分析 3.1模型建立 图3 有限元模型 3.2边界条件和载荷 为避免在计算时平板产生移动引发计算问题，必须对试件的外部边界条件进行限定。对平板左侧进行铰接约束，示意图如下
图4 平板约束示意图 由于我们只关注孔附近的应力分布情况，根据圣维南原理，载荷的具体分布只影响载荷作用区附近的应力分布。故我们用均布力代替集中力施加在平板右侧的作用面上，其大小为225P MPa ，为负值。 图5 平板载荷示意图 3.3材料 平板的弹性模量为200GPa ，泊松比为0.3。其塑性的应力应变参数见下图 图6 塑性应力应变参数 3.4有限元网格划分 网格划分是非常重要的过程，它会对计算速度、精度、可靠性产生重要影响。网格划分主要包括两方面：尺寸、单元类型。
Matlab有限元分析操作基础
Matlab 有限元分析20140226 为了用Matlab 进行有限元分析，首先要学会Matlab 基本操作，还要学会使用Matlab 进行有限元分析的基本操作。 1. 复习：上节课分析了弹簧系统 x 推导了系统刚度矩阵 11221 21200k k k k k k k k -????-????--+??
2. Matlab有限元分析的基本操作 （1）单元划分（选择何种单元，分成多少个单元，标号）（2）构造单元刚度矩阵（列出…） （3）组装系统刚度矩阵（集成整体刚度矩阵） （4）引入边界条件（消除冗余方程） （5）解方程 （6）后处理（扩展计算）
3. Matlab有限元分析实战【实例1】
分析： 步骤一：单元划分
步骤二：构造单元刚度矩阵 >>k1=SpringElementStiffness(100) >>…?
步骤三：构造系统刚度矩阵 a) 分析SpringAssemble库函数function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵？ 今天的系统刚度矩阵是什么？ 因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ?? ?? -- ?? ？
abaqus 有限元分析(齿轮轴)
Abaqus分析报告 （齿轮轴） 名称：Abaqus齿轮轴 姓名： 班级： 学号： 指导教师：
一、简介 所分析齿轮轴来自一种齿轮泵，通过用abaqus软件对齿轮轴进行有限元分析和优化。齿轮轴装配结构图如图1，分析图1中较长的齿轮轴。 图1.齿轮轴装配结构图 二、模型建立与分析 通过part、property、Assembly、step、Load、Mesh、Job等步骤建立齿轮轴模型，并对其进行分析。 1.part 针对该齿轮轴，拟定使用可变型的3D实体单元，挤压成型方式。 2.材料属性 材料为钢材，弹性模量210Gpa，泊松比0.3。
3.截面属性 截面类型定义为solid，homogeneous。 4.组装 组装时选择dependent方式。 5.建立分析步 本例用通用分析中的静态通用分析（Static，General）。 6.施加边界条件与载荷 对于齿轮轴，因为采用静力学分析，考虑到前端盖、轴套约束，而且根据理论，对受力部分和轴径突变的部分进行重点分析。 边界条件：分别在三个轴径突变处采用固定约束，如图2。 载荷：在Abaqus中约束类型为pressure，载荷类型为均布载荷,分别施加到齿轮接触面和键槽面，根据实际平衡情况，两力所产生的绕轴线的力矩方向相反，大小按比例分配。 均布载荷比计算： 矩形键槽数据： 长度：8mm、宽度：5mm、高度：3mm、键槽所在轴半径：7mm 键槽压力面积：S1 = 8x3=24mm2 平均受力半径：R1=6.5mm 齿轮数据：= 齿轮分度圆半径：R2 =14.7mm、压力角：20°、 单个齿轮受力面积：S2 ≈72mm2 通过理论计算分析，S1xR1xP1=S2xR2xP2，其中，P1为键槽均布载荷
机械零件有限元分析——实验报告
中南林业科技大学机械零件有限元分析 实验报告 专业：机械设计制造及其自动化 年级： 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
支架的有限元分析ABAQUS
支架的线性静力学分析实例：建模和分析计算 在此实例中读者将学习ABAQUS/CAE的以下功能。 1) Sketch功能模块:导人CAD二维图形，绘制线段、圆弧和倒角，添加尺寸，修改平面图，输出平面图。 2) Part功能模块:通过拉伸来创建几何部件，通过切割和倒角未定义几何形状。 3) Property功能模块:定义材料和截面属性。 4) Mesh功能模块:布置种子，分割实体和面，选择单元形状、单元类型、网格划分 技术和算法，生成网格，检验网格质量，通过分割来定义承受载荷的面。 5) Assembly功能模块:创建非独立实体。 6) Step功能模块:创建分析步，设置时间增量步和场变量输出结果。 7) Interaction功能模块:定义分布榈合约束(distributing coupling constraint)。 8) Load功能模块:定义幅值，在不同的分析步中分别施加面载荷和随时间变化的集中力，定义边界条件。 9) Job功能模块:创建分析作业，设置分析作业的参数，提交和运行分析作业，监控运行状态。 10) Visualization功能模块:后处理的各种常用功能。 结构静力学分析（static analysis）是有限元法的基本应用领域，适用于求解惯性及阻尼对结构响应不显著的问题。主要用来分析由于稳态外载荷引起的位移，应力和应变等。本章的静力学分析实例按照ABAQUS工程分析的流程对支架进行线性静力学分析，通过实例基本掌握了分析的流程，同时了解接触的定义。 1.问题描述 所示的支架，一端牢固地焊接在一个大型结构上，支架的圆孔中穿过一个相对较软的杆件，圆孔和杆件用螺纹连接。材料的弹性模量E=2100000MPa,泊松比为0.3。
有限元实验报告模板
有限元实验报告 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 的均匀布置的载荷。
《有限元基础教程》_【MATLAB算例】3.3.7(2)__三梁平面框架结构的有限元分析(Beam2D2Node)
【MA TLAB 算例】3.3.7(2) 三梁平面框架结构的有限元分析 (Beam2D2Node) 如图3-19所示的框架结构，其顶端受均布力作用，结构中各个 截面的参数都为：113.010Pa E =?，746.510I m -=?，426.810A m -=?。试基 于MA TLAB 平台求解该结构的节点位移以及支反力。 图3-19 框架结构受一均布力作用 解答：对该问题进行有限元分析的过程如下。 （1） 结构的离散化与编号 将该结构离散为3个单元，节点位移及单元编号如图3-20所示， 有关节点和单元的信息见表3-5。 （a ） 节点位移及单元编号
（b）等效在节点上的外力 图3-20 单元划分、节点位移及节点上的外载 （2）各个单元的描述 首先在MA TLAB环境下，输入弹性模量E、横截面积A、惯性矩I和长度L，然后针对单元1，单元2和单元3，分别二次调用函数Beam2D2Node_ElementStiffness，就可以得到单元的刚度矩阵k1(6×6)和k2(6×6)，且单元2和单元3的刚度矩阵相同。 >> E=3E11; >> I=6.5E-7; >> A=6.8E-4; >> L1=1.44; >> L2=0.96; >> k1=Beam2D2Node_Stiffness(E,I,A,L1); >> k2=Beam2D2Node_Stiffness(E,I,A,L2); （3）建立整体刚度方程 将单元2和单元3的刚度矩阵转换成整体坐标下的形式。由于该结构共有4个节点，则总共的自由度数为12，因此，结构总的刚度矩阵为KK(12×12)，对KK清零，然后两次调用函数Beam2D2Node_Assemble进行刚度矩阵的组装。 >> T=[0,1,0,0,0,0;-1,0,0,0,0,0;0,0,1,0,0,0;0,0,0,0,1,0;0,0,0,-1,0,0;0,0,0,0,0,1] ; >> k3=T'*k2*T; >> KK=zeros(12,12); >> KK=Beam2D2Node_Assemble(KK,k1,1,2);
(完整版)有限元大作业matlab---课程设计例子
有限元大作业程序设计 学校：天津大学 院系：建筑工程与力学学院 专业：01级工程力学 姓名：刘秀 学号：\\\\\\\\\\\ 指导老师：
连续体平面问题的有限元程序分析 [题目]： 如图所示的正方形薄板四周受均匀载荷的作用，该结构在边界 上受正向分布压力， m kN p 1=，同时在沿对角线y 轴上受一对集中压 力，载荷为2KN ，若取板厚1=t ，泊松比0=v 。 [分析过程]： 由于连续平板的对称性，只需要取其在第一象限的四分之一部分参加分析，然后人为作出一些辅助线将平板“分割”成若干部分，再为每个部分选择分析单元。采用将此模型化分为4个全等的直角三角型单元。利用其对称性，四分之一部分的边界约束，载荷可等效如图所示。
[程序原理及实现]： 用FORTRAN程序的实现。由节点信息文件NODE.IN和单元信息文件ELEMENT.IN，经过计算分析后输出一个一般性的文件DATA.OUT。模型基本信息由文件为BASIC.IN生成。 该程序的特点如下： 问题类型：可用于计算弹性力学平面问题和平面应变问题 单元类型：采用常应变三角形单元 位移模式：用用线性位移模式 载荷类型：节点载荷，非节点载荷应先换算为等效节点载荷 材料性质：弹性体由单一的均匀材料组成 约束方式：为“0”位移固定约束，为保证无刚体位移，弹性体至少应有对三个自由度的独立约束 方程求解：针对半带宽刚度方程的Gauss消元法
输入文件：由手工生成节点信息文件NODE.IN，和单元信息文件ELEMENT.IN 结果文件：输出一般的结果文件DATA.OUT 程序的原理如框图：
带孔平板的线性静力分析
带孔平板的线性静力分析 本示例将对一个给定的带孔平板几何模型创建有限元模型、施加边界条件、进行有限元分析并在HyperView中观察受载平板的变形和应力结果。 本示例包括以下步骤： ?在HyperMesh中建立有限元模型 ?施加载荷和边界条件 ?求解 ?观察结果 1．在HyperMesh中建立有限元模型 （1）载入OptiStruct用户界面并打开模型文件 1）启动HyperMesh。 2）在User Profile对话框中选择OptiStruct，点击OK。 这就加载了OptiStruct用户界面，它包括OptiStruct模板、宏菜单等。简化了与OptiStruct 使用相关的HyperMesh功能。 User Profiles…可以从下拉式菜单中的Preferences中进入。 3）在工具条选择按钮。 弹出Open file…窗口。 4）选择plate_hole.hm文件，模型位于
Matlab有限元分析操作基础共11页
Matlab有限元分析20140226 为了用Matlab进行有限元分析，首先要学会Matlab基本操作，还要学会使用Matlab进行有限元分析的基本操作。 1. 复习：上节课分析了弹簧系统 x 推导了系统刚度矩阵
2. Matlab有限元分析的基本操作 （1）单元划分（选择何种单元，分成多少个单元，标号）（2）构造单元刚度矩阵（列出…） （3）组装系统刚度矩阵（集成整体刚度矩阵） （4）引入边界条件（消除冗余方程） （5）解方程 （6）后处理（扩展计算）
3. Matlab有限元分析实战【实例1】
分析： 步骤一：单元划分
>>k1=SpringElementStiffness(100)
a) 分析SpringAssemble库函数 function y = SpringAssemble(K,k,i,j) % This function assembles the element stiffness % matrix k of the spring with nodes i and j into the % global stiffness matrix K. % function returns the global stiffness matrix K % after the element stiffness matrix k is assembled. K(i,i) = K(i,i) + k(1,1); K(i,j) = K(i,j) + k(1,2); K(j,i) = K(j,i) + k(2,1); K(j,j) = K(j,j) + k(2,2); y = K; b) K是多大矩阵？ 今天的系统刚度矩阵是什么？ 因为 11 22 1212 k k k k k k k k - ?? ?? - ????--+ ?? 所以 1000100 0200200 100200300 - ?? ?? - ????-- ??？
有限元实验报告
一、实验目的 通过上机对有限元法的基本原理和方法有一个更加直观、深入的理解；通过对本实验所用软件平台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
带孔板的建模及有限元分析Word版
基于SolidWorks带孔板的建模及有限元分析 李军 摘要：利用SolidWorks对带孔矩形板进行虚拟建模，通过赋予板材材质、载荷后进行网格划分，进而进行有限元分析，得出其应力、应变和位移的分布图，并对结果进行分析研究对板材安全性的影响。 关键词：SolidWorks；带孔板；建模；有限元分析 0 SolidWorks简介 Solidworks是一款优秀的三维设计软件，具有十分强大的零件设计功能及装配模块，同时也拥有丰富的后置处理模块。由于其功能强大，新手上手快，应用领域广，所以成为了主流的三维造型软件。经过17年的发展，在全球已经拥有30多万的客户，最新版本为SolidWorks 2011版。在中国SolidWorks在计算机辅助设计、计算机辅助工程、计算机辅助制造、计算机辅助工艺、数据管理等方面为企业提供了强大的动力，使企业在管理、设计和制造方面有了很大的提升。 1 带孔板的模型建立 矩形板材的尺寸为300*180*10mm，孔位于中心，直径为50mm，模型如图1。 图1 带孔矩形板模型 2前置处理 2.1在Command Manager中点击SIMULATION选项，建立新算例，名称默认，确认。 2.2赋予板材材料属性 材料为AISI304，材料属性如表1
表1 材料的属性 模型参考属性零部件 名称:AISI 304 模型类型:线性弹性同向性 默认失败准则:最大von Mises 应力屈服强度: 2.06807e+008 N/m^2 张力强度: 5.17017e+008 N/m^2 弹性模量: 1.9e+011 N/m^2 泊松比:0.29 质量密度:8000 kg/m^3 抗剪模量:7.5e+010 N/m^2 热扩张系数: 1.8e-005 /Kelvin SolidBody 1(凸台-拉伸1)(aisi304带孔矩形钢板静力分析) 曲线数据:N/A 2.3网格生成 在SIMULATION选项中选择“运行”中的“生成网格”，使用默认网格划分。网格 信息如表2，网格信息细节如表3，网格划分后的模型如图2。 表2 网格信息 网格类型实体网格 所用网格器: 基于曲率的网格 雅可比点 4 点 最大单元大小7.44196 mm 最小单元大小7.44196 mm 网格品质高 表3 网格信息细节 节点总数23523 单元总数13612 最大高宽比例 3.9347 单元(%)，其高宽比例< 3 99.7 单元(%)，其高宽比例> 10 0 扭曲单元(雅可比)的% 0 完成网格的时间(时;分;秒): 00:00:03 计算机名: PC-201009062016
基于matlab的有限元法分析平面应力应变问题刘刚
姓名：刘刚学号：15 平面应力应变分析有限元法 Abstruct:本文通过对平面应力/应变问题的简要理论阐述，使读者对要分析的问题有大致的印象，然后结合两个实例，通过MATLAB软件的计算，将有限元分析平面应力/应变问题的过程形象的展示给读者，让人一目了然，快速了解有限元解决这类问题的方法和步骤！ 一.基本理论 有限元法的基本思路和基本原则以结构力学中的位移法为基础，把复杂的结构或连续体看成有限个单元的组合，各单元彼此在节点出连接而组成整体。把连续体分成有限个单元和节点，称为离散化。先对单元进行特性分析，然后根据节点处的平衡和协调条件建立方程，综合后做整体分析。这样一分一合，先离散再综合的过程，就是把复杂结构或连续体的计算问题转化简单单元分析与综合问题。因此，一般的有限揭发包括三个主要步骤：离散化单元分析整体分析。 二.用到的函数 1. LinearTriangleElementStiffness(E,NU,t,xi,yi,xj,yj,xm,ym,p) (K k I f) (k u) (k u A) (E NU t) 三.实例 例1.考虑如图所示的受均布载荷作用的薄平板结构。将平板离散化成两个线性三角元，假定E=200GPa，v=，t=0.025m,w=3000kN/m. 1.离散化 2.写出单元刚度矩阵
通过matlab 的LinearTriangleElementStiffness 函数，得到两个单元刚度矩阵1k 和2k ，每个矩阵都是6 6的。 >> E=210e6 E = >> k1=LinearTriangleElementStiffness(E,NU,t,0,0,,,0,,1) k1 = +006 * Columns 1 through 5 0 0 0 0 0 0 0 0 Column 6 >> NU= NU = >> t= t = >> k2=LinearTriangleElementStiffness(E,NU,t,0,0,,0,,,1)
ABAQUS有限元接触分析的基本概念
ABAQUS有限元接触分析的基本概念2009-11-24 00:06:28 作者：jiangnanxue 来源：智造网—助力中国制造业创新—https://www.sodocs.net/doc/e810275517.html, CAE(计算机辅助工程)是一门复杂的工程科学，涉及仿真技术、软件、产品设计和力学等众多领域。世界上几大CAE公司各自以其独到的技术占领着相应的市场。ABAQUS有限元分析软件拥有世界上最大的非线性力学用户群，是国际上公认的最先进的大型通用非线性有限元分析软件之一。它广泛应用于机械制造、石油化工、航空航天、汽车交通、土木工程、国防军工、水利水电、生物医学、电子工程、能源、地矿、造船以及日用家电等工业和科学研究领域。ABAQUS在技术、品质和可靠性等方面具有卓越的声誉，可以对工程中各种复杂的线性和非线性问题进行分析计算。 《ABAQUS有限元分析常见问题解答》以问答的形式，详细介绍了使用ABAQUS建模分析过程中的各种常见问题，并以实例的形式教给读者如何分析问题、查找错误原因和尝试解决办法，帮助读者提高解决问题的能力。 《ABAQUS有限元分析常见问题解答》一书由机械工业出版社出版。 16.1.1 点对面离散与面对面离散 【常见问题16-1】 在ABAQUS/Standard分析中定义接触时，可以选择点对面离散方法(node-to-surface-dis - cre-tization)和面对面离散方法(surface-to-surface discretization)，二者有何差别? 『解答』 在点对面离散方法中，从面(slave surface)上的每个节点与该节点在主面(master surface)上的投影点建立接触关系，每个接触条件都包含一个从面节点和它的投影点附近的一组主面节点。 使用点对面离散方法时，从面节点不会穿透(penetrate)主面，但是主面节点可以穿透从面。 面对面离散方法会为整个从面(而不是单个节点)建立接触条件，在接触分析过程中同时考虑主面和从面的形状变化。可能在某些节点上出现穿透现象，但是穿透的程度不会很严重。 在如图16-l和图16-2所示的实例中，比较了两种情况。
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=-+-=σ
有限元上机实验报告
有限元法基础及应用 上机报告 南京理工大学 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|=常数的条件下能被精确积分，并能保证刚度矩阵的非奇异性。如果采用减缩积分，因为插值函数中只有完全多项式的项次能被精确积分，因此需要进行刚度矩阵非奇异必要条件的检查。