当前位置：文档视界 › 计算结构力学有限元方法_三维结构和轴对称

计算结构力学有限元方法_三维结构和轴对称

结构力学思考题答案

1、结构的动力特性一般指什么？答：结构的动力特性是指：频率（周期）、振型和阻尼。动力特性是结构固有的，这是因为它们是由体系的基本参数（质量、刚度）所确定的、表征结构动力响应特性的量。动力特性不同，在振动中的响应特点亦不同。 2、什么是阻尼、阻尼力，产生阻尼的原因一般有哪些？什么是等效粘滞阻尼？答：振动过程的能量耗散称为阻尼。产生阻尼的原因主要有：材料的内摩擦、构件间接触面的摩擦、介质的阻力等等。当然，也包括结构中安装的各种阻尼器、耗能器。阻尼力是根据所假设的阻尼理论作用于质量上用于代替能量耗散的一种假想力。粘滞阻尼理论假定阻尼力与质量的速度成比例。粘滞阻尼理论的优点是便于求解，但其缺点是与往往实际不符，为扬长避短，按能量等效原则将实际的阻尼耗能换算成粘滞阻尼理论的相关参数，这种阻尼假设称为等效粘滞阻尼。 3、采用集中质量法、广义位移法（坐标法）和有限元法都可使无限自由度体系简化为有限自由度体系，它们采用的手法有何不同？答：集中质量法：将结构的分布质量按一定规则集中到结构的某个或某些位置上，认为其他地方没有质量。质量集中后，结构杆件仍具有可变形性质，称为“无重杆”。广义坐标法：在数学中常采用级数展开法求解微分方程，在结构动力分析中，也可采用相同的方法求解，这就是广义坐标法的理论依据。所假设的形状曲线数目代表在这个理想化形式中所考虑的自由度个数。考虑了质点间均匀分布质量的影响（形状函数），一般来说，对于一个给定自由度数目的动力分析，用理想化的形状函数法比用集中质量法更为精确。有限元法：有限元法可以看成是广义坐标法的一种特殊的应用。一般的广义坐标中，广义坐标是形函数的幅值，有时没有明确的物理意义，并且在广义坐标中，形状函数是针对整个结构定义的。而有限元法则采用具有明确物理意义的参数作为广义坐标，且形函数是定义在分片区域的。在有限元分析中，形函数被称为插值函数。综上所述，有限元法综合了集中质量法和广义坐标法的特点： (l) 与广义坐标法相似，有限元法采用了形函数的概念。但不同于广义坐标法在整体结构上插值（即定义形函数），而是采用了分片的插值，因此形函数的表达式（形状）可以相对简单。 (2) 与集中质量法相比，有限元法中的广义坐标也采用了真实的物理量，具有直接、直观的优点，这与集中质量法相同。 4、直接动力平衡法中常用的有哪些具体方法？它们所建立的方程各代表什么条件？答：常用方法有两种：刚度法和柔度法。刚度法方程代表的是体系在满足变形协调条件下所应满足的动平衡条件；而柔度法方程则代表体系在满足动平衡条件下所应满足的变形协调条件。 5、刚度法与柔度法所建立的体系运动方程间有何联系？各在什么情况下使用方便？答：刚度法与柔度法建立的运动方程在所反映的各量值之间的关系上是完全一致的。由于刚度矩阵与柔度矩阵互逆，刚度法建立的运动方程可转化为柔度法建立的方程。一般来，对于单自由度体系，求[δ]和求[k]的难易程度是相同的，因为它们互为倒数，都可以用同一方法求得，不同的是一个已知力求位移，一个已知位移求力。对于多自由度体系，若是静定结构，一般情况下求柔度系数容易些，但对于超静定结构就要根据具体情况而定。若仅从建立运动方程来看，当刚度系数容易求时用刚度法，柔度系数容易求时用柔度法。 6、计重力与不计重力所得到的运动方程是一样的吗？答：如果计与不计重力时都相对于无位移的位置来建立运动方程，则两者是不一样的。但如果计重力时相对静力平衡位置来建立运动方程，不计重力仍相对于无位移位置来建立，

有限元ansys静力分析的一个小例子

有限元学院：机电学院专业：姓名：学号：

一、问题描述如图所示的平面，板厚为0.01m，左端固定，右端作用50kg的均布载荷，对其进行静力分析。弹性模量为210GPa，泊松比为0.25. 二、分析步骤 1．启动ansys，进入ansys界面。 2．定义工作文件名进入ANSYS/Multiphsics的的程序界面后，单击Utility Menu菜单下File中Change Jobname的按钮，会弹出Change Jobname对话框，输入gangban为工作文件名，点击ok。 3．定义分析标题选择菜单File-Change Title在弹出的对话框中，输入Plane Model作为分析标题，单击ok。 4．重新显示选择菜单Plot-Replot单击该按钮后，所命令的分析标题工作文件名出现在ANSYS 中。 5．选择分析类型在弹出的对话框中，选择分析类型，由于此例属于结构分析，选择菜单Main Menu：Preferences，故选择Structural这一项，单击ok。 6．定义单元类型选择菜单Main Menu-Preprocessor-Element Type-Add/Edit/Delete单击弹出对话框中的Add按钮，弹出单元库对话框，在材料的单元库中选Plane82单元。即在左侧的窗口中选取Solid单元，在右侧选择8节点的82单元。然后单击ok。 7．选择分析类型定义完单元类型后，Element Type对话框中的Option按钮被激活，单击后弹出一个对话框，在Elenment behavior中选择Plane strs w/ thk，在Extra Element output 中，选择Nodal stress，单击close，关闭单元类型对话框。 8．定义实常数选择菜单Main Menu-Preprocessor-Real Constants Add/Edit/Delete执行该命令后，在弹出Real Constants对话框中单击Add按钮，确认单元无误后，单击ok，弹出Real Constants Set Number 1，for Plane 82对话框，在thickness后面输入板的厚度0.01单击ok，单击close。 9．定义力学参数选择菜单Main Menu-Preprocessor-Material Props-Material Model 在弹出的对

三维有限元法计算过程

三维有限元法计算过程三维有限元法的计算过程： 1）网格单元剖分； 2）线性插值； 3）单元分析； 4）总体刚度矩阵合成； 5）求解线性方程组等部分组成。一、偏微分方程对应泛函的极值问题矿井稳恒电流场分布示意图主要任务是分析在给定边界条件下，求解稳定电流场的Laplace 方程或Poisson方程的数值解，即三维椭圆型微分方程的边值问题：

) ()((0)(0 )()()(000z z y y x x I F u n u n u F z u z y u y x u x Lu w D ---=???? ?????=+??=??=????+????+????≡ΓΓ+Γδδδγσσσ 上述微分方程边值问题等价于下面泛函的极小值问题： dS U dxdydz fU z U y U x U U J w D ?????Γ+Γ+ΓΩ +-??+??+??=222221 }])()()[(2{][γσσ 二、网格剖分 ∞1 ρi i h ρ......... ... 1、网格单元的类型图2-5 网格单元类型 2、网格单元剖分原则及其步长选择因此，网格内的单元剖分应按以下剖分原则 1）、各单元节点（顶点）只能与相邻单元节点（顶点）重合，而

不能成为其它单元内点； 2）、如果求解区域对称，那么单元剖分也应该对称； 3）、在场变化剧烈的区域网格剖分单元要密一些，在场变化平缓的区域单元密度应小。 4）、网格单元体的大小变化应逐步过渡。根据上述剖分原则，以x 、y 、z 坐标轴原点o 为中心，分别向x 、y 、z 方向的两侧作对称变步长剖分，距o 越远，步长应越大。常用的变步长方法有： c i x x i i )1(1+=?-?+ c x x i i =??+/1(i ≠0) c x x i i =?-?+1 1 1(i ≠0) 以上各式中c 为常数，1+?i x 、i x ?为同一坐标轴上相邻步长值。以x 方向为例，可知，x 正方向与负方向对称，只相差一负号。若令00=?x ，只要给出距原点最近节点的坐标1x ?，由上式即可求出其它相应的步长i x ?。同理可求得y 、 z 方向上的变步长i y ?、i z ?。 3、网格剖分方法图2-6 平面内节点编号示意图

弹性力学及有限元基础复习权威版(最新)

《弹性力学及有限元基础》复习思考题 ★1．对弹性体所做的基本假设? 答：连续性假设；均匀性假设；各向同性假设；弹性假设；小变形假设； ★2．用D'Alember 原理由平衡方程推导运动微分方程? 答：微元体的平衡微分方程的表达式为： 31 112111 2332 122221 23 132333 31 23000f x x x f x x x f x x x σσσσσσσσσ????+++=?????????+++=? ????????+++=? ???? 根据D'Alember 原理，将运动物体看成是静止的，将惯性力22()u t ρ?-?当作体力加到微元体上，由上式可以直接写出弹性动力学问题的运动微分方程： 23111211 12123232 12222221 2321323333321 23()()() u f x x x t u f x x x t u f x x x t σσσρσσσρσσσρ?????+++=????????????+++=? ???????????+++=?????? ☆3．什么是应力张量? 我们说一点的应力状态是什么涵义? 答：应力张量是一点应力状态的完整描述，它有面元方向和分解方向两个方向性，共有九个分量，由于存在对称性，其独立分量只有六个。应力张量是与坐标选择无关的不变量，但其分量与坐标有关，当已知某坐标系中的九个分量时，其他坐标系中的分量均可由应力转换公式确定。一点的应力状态是一个具有双重方向性的物理量，其中第一个是面元的方向，用其法矢量ν表示，第二个是作用在该面元上的应力矢量方向，一般用其三个分量来表示。 4．在引出 Cauchy 应力公式时，我们假设四面体处于平衡状态，如不处在平衡状态则如何? 答：如果不处在平衡状态，Cauchy 应力公式仍然满足，关系式的成立与是否平衡无关。 5．在什么情况下剪应力互等定律不成立? 答：无论在变形体的内部或者表面上，若存在体力偶时，剪应力互等定律不成立。 6．任意斜截面上的正应变和剪应变的意义是什么? 答：应变张量的三个对角分量x ε、y ε、z ε称为正应变，分别等于坐标轴方向三个线元的单位伸长率，伸长为正，缩短为负。应变张量的三个非对角分量xy ε、yz ε、zx ε称为剪应变，分别等于变形前沿该分量下标所示两坐标方向的、相互正交的线元在变形后的夹角减小量之半。 7．刚性位移，刚性转动，刚体位移，刚体转动有何区别? 答：（1）刚性位移：物体内任意两点间无相对位移；（2）刚性转动：应变张量为0，转动张量不为0；（3）刚体位移：运动分为变形运动和刚体运动，每点都发生相同的位移就叫作刚体位移；（4）刚体转动：用刚性

结构力学

张永发结构力学（二）·平时作业， 1. 请简述结构力学在土木工程领域的作用。答：结构力学是固体力学的一个分支，它主要研究工程结构受力和传力的规律，以及如何进行结构优化的学科，它是土木工程专业和机械类专业学生必修的学科。结构力学研究的内容包括结构的组成规则，结构在各种效应（外力，温度效应，施工误差及支座变形等）作用下的响应，包括内力（轴力，剪力，弯矩，扭矩）的计算，位移（线位移，角位移）计算，以及结构在动力荷载作用下的动力响应（自振周期，振型）的计算等。结构力学通常有三种分析的方法：能量法，力法，位移法，由位移法衍生出的矩阵位移法后来发展出有限元法，成为利用计算机进行结构计算的理论基础。结构力学是一门古老的学科，又是一门迅速发展的学科。新型工程材料和新型工程结构的大量出现，向结构力学提供了新的研究内容并提出新的要求。计算机的发展，又为结构力学提供了有力的计算工具。另一方面，结构力学对数学及其他学科的发展也起了推动作用。有限元法这一数学方法的出现和发展就和结构力学的研究有密切关系。在固体力学领域中，材料力学给结构力学提供了必要的基本知识，弹性力学和塑性力学是结构力学的理论基础。另外，结构力学与流体力学相结合形成边缘学科——结构流体弹性力学。评定结构的优劣，从力学角度看，主要是结构的强度和刚度。工程结构设计既要保证结构有足够的强度，又要保证它有足够的刚度。强度不够，结构容易破坏；刚度不够，结构容易皱损，或出现较大的振动，或产生较大的变形。皱损能够导致结构的变形破坏，振动能够缩短结构的使用寿命，皱损、振动、变形都会影响结构的使用性能，例如，降低机床的加工精度或减低控制系统的效率等。观察自然界中的天然结构，如植物的根、茎和叶，动物的骨骼，蛋类的外壳，可以发现它们的强度和刚度不仅与材料有关，而且和它们的造型有密切的关系。很多工程结构是受到天然结构的启发而创制出来的。人们在结构力学研究的基础上，不断创造出新的结构造型。加劲结构(见加劲板壳)、夹层结构(见夹层板壳)等都是强度和刚度比较高的结构。结构设计不仅要考虑结构的强度和刚度，还要做到用料省、重量轻。减轻重量对某些工程尤为重要，如减轻飞机的重量就可以使飞机航程远、上升快、速度大、能耗低。 2. 请简述对称性对力学分析中的影响，并举例说明答：在工程实际问题中，有很多结构都具有对称性。我们对这些结构进行受力分析的时候，常常将结构简化为杆系模型，而结构力学研究的就是结构的杆系模型，因此对称性在结构力学中有着广泛的应用。特别是在求解超静定结构问题中，无论力法还是位移法，都是繁杂的。但对于对称结构，利用结构的对称性，可使结构内力计算大为简化。现在本文章就对称性在结构力学中的应用做一简单的总结。

111ANSYS进行有限元静力学分析

经典理论一、设计大纲概述 1、设计目的（1）熟悉有限元分析的基本原理和基本方法；（2）掌握有限元软件ANSYS的基本操作；（3）对有限元分析结果进行正确评价。 2、设计原理利用ANSYS进行有限元静力学分析。 3、设计仪器设备 1）安装windows 2000以上版本的微机； 2）ANSYS 8.0以上版本软件。 4、实验内容与步骤 1）熟悉ANSYS的界面和分析步骤； 2）掌握ANSYS前处理方法，包括平面建模、单元设置、网格划分和约束设置； 3）掌握ANSYS求解和后处理的一般方法； 4）实际应用ANSYS软件对平板结构进行有限元分析。二、题目：如图试样期尺寸为100mm*5mm*5mm,下端固定，上端受拉力10000N作用。已知该试样材料的应力-应变曲线如图所示。计算试样的位移分布。

三、分析步骤：分析：从应力-应变关系可以看出该材料的屈服极限是225MPa 左右，弹性部分曲线的斜率为常数75GPa。之后材料进入塑性变形阶段，应力-应变关系为非线性的。估计本题应力10000/（0.05*.005）=400MPa,因此材料屈服进入塑性，必须考虑材料非线性影响。（1）建立关键点。单击菜单Main Menu>Preprocessor>Modeling>Create>Keypoints>In ActiveCS，建立两个关键点（0，0，0）和（0，100， 0）。（2）建立直线。单击菜单Main Menu>Preprocessor>Modeling>Create>Lines>Staight Line,在关键点1、2之间建立直线。（3）定义单元类型。单击菜单Main Menu>Preprocessor>ElementType>Add/Edit/Delete, 定义单元Structural>Link>2D spar1(LINK1) （4）定义单元常数。单击菜单Main Menu>Preprocessor>RealConstants>Add/Edit/Delete,

弹性力学及有限元法学习总结

弹性力学及有限元法学习总结摘要：本文就弹性力学的研究对象与方法，弹性力学的基本假设，研究方法，有限元法的基本思想，数学基础，有限元分析的基本步骤进行阐述。正文：弹性力学是固体力学的一个分支学科，是研究固体材料在外部作用下（外部作用一般包括：荷载、温度变化以及固体边界约束改变），弹性变形及应力状态的一门学科。弹性力学的研究对象：材料力学--研究杆件（如梁、柱和轴）材料力学的拉压、弯曲、剪切、扭转和组合变形等问题。结构力学--在材料力学基础上研究杆系结构结构力学（如桁架、刚架等）。弹性力学--研究各种形状的弹性体，如杆弹性力学件、平面体、空间体、板壳、薄壁结构等问题。弹性力学研究方法：在研究方法上，弹力和材力也有区别：弹力研究方法：在区域V内严格考虑静力学、几何学和物理学三方面条件，建立三套方程; 三套方程在边界s 上考虑受力或约束条件，建立边界条件并在边界条件下求解上边界条件; 边界条件述方程，得出较精确的解答。弹性力学的基本假设： 1）连续性，假定物体是连续的。连续性因此，各物理量可用连续函数表示。 2）均匀性与各向同性假设假定固体材料是均匀的，并且在各个方向上物理特性相同，也即材料的物理性质在空间分布上是均匀的（或不变的）3）小变形假设假定固体材料在受到外部作用（荷载、温度等）后的位移（或变形）与物体的尺寸相比是很微小的，在研究物体受力后的平衡状态时，物体尺寸及位置的改变可忽略不计，物体位移及形变的二次项可略去不计，由此得到的弹性力学微分方程将是线性的。 4）完全弹性假设假设固体材料是完全弹性的。 5）无初始应力假设假定外部作用（荷载、温度等）之前，物体处于无应力状态，由弹性力学所求得的应力仅仅是由外部作用（荷载、温度等）所引起的。有限元法的基本思想：有限元是一种结构分析的方法，先把所有系统分解为他们的元件或单元，这些元件的行为已经被充分的了解，再把元件重新组装成原来的系统。及将连续的求解区域离散为一组由有限个单元组成并按一定方式相互连接在一起的单元组

三维有限元方法-为一种新型的研究方法,是利用数学的形式概括事件的各种条件和性能

三维有限元方法-为一种新型的研究方法，是利用数学的形式概括事件的各种条件和性能三维有限元方法-为一种新型的研究方法，是利用数学的形式概括事件的各种条件和性能，并进行重复分析计算的研究方法。有限元数值模型分析技术将现代数学、力学的基础理论与有限元分析技术、计算机图形学和优化技术相结合，具有丰富、完善的单元库、材料模型库和求解器，可利用数值模拟技术高效求解各类结构动力、静力和线性、非线性问题。将其应用于骨科领域，可以更好的进行各种骨科生物力学分析，对各种生物力学强度进行数值模拟分析，较精确地掌握各点的受力情况，了解内部应力应变的分布规律，获得应力应变分布图等，从而更好的指导临床治疗。学术术语来源—— 锁骨中段骨折修复：重建钢板前置与上置的生物力学差异文章亮点： 1 文章结果显示，不论怎样的载荷条件，骨折断端均会存在一定的应力。而且，前置位和上置位不同内固定方式对骨折端愈合的影响方面不存在明显差别，但在骨折断端应力和内固定应力方面，前置位均显著大于上置位。即提示，较之上置位，前置位固定具有更明显的应力集中效应。 2 临床对锁骨中段骨折进行修复的过程中，利用不同重建钢板位置进行内固定修复会产生不同的生物力学情况。其中，较之重建钢板上置内固定，前置内固定修复效果更佳，是一种较为可靠的治疗方法。 3 文章仅对不同重建钢板位置的内固定效果进行了分析研究，并未考虑到不同骨折类型力学特性以及不同钢板类型等因素的影响。并在研究过程中假设螺钉为圆形杆，因此最终的研究结果可能存在导致内固定装置最大等效应力下降的情况，计算精确度存在一定的误差。另外，文章中对所使用的各种生物材料的力学特性均进行了假设，与客观情况存在较大的差异。因此，文章还存在一定的不足之处，还需要在今后的研究中不断予以完善，以提高研究结果的准确性和可信性，更好的为临床治疗提供可参考的依据。关键词：

(完整版)结构力学问答题总结

概念题 1.1 结构动力计算与静力计算的主要区别是什么？答：主要区别表现在：(1) 在动力分析中要计入惯性力，静力分析中无惯性力；(2) 在动力分析中，结构的内力、位移等是时间的函数，静力分析中则是不随时间变化的量；(3) 动力分析方法常与荷载类型有关，而静力分析方法一般与荷载类型无关。 1.2 什么是动力自由度，确定体系动力自由度的目的是什么？答：确定体系在振动过程中任一时刻体系全部质量位置或变形形态所需要的独立参数的个数，称为体系的动力自由度（质点处的基本位移未知量）。确定动力自由度的目的是：(1) 根据自由度的数目确定所需建立的方程个数（运动方程数=自由度数），自由度不同所用的分析方法也不同；(2) 因为结构的动力响应（动力内力和动位移）与结构的动力特性有密切关系，而动力特性又与质量的可能位置有关。 1.3 结构动力自由度与体系几何分析中的自由度有何区别？答：二者的区别是：几何组成分析中的自由度是确定刚体系位置所需独立参数的数目，分析的目的是要确定体系能否发生刚体运动。结构动力分析自由度是确定结构上各质量位置所需的独立参数数目，分析的目的是要确定结构振动形状。 1.4 结构的动力特性一般指什么？答：结构的动力特性是指：频率（周期）、振型和阻尼。动力特性是结构固有的，这是因为它们是由体系的基本参数（质量、刚度）

所确定的、表征结构动力响应特性的量。动力特性不同，在振动中的响应特点亦不同。 1.5 什么是阻尼、阻尼力，产生阻尼的原因一般有哪些？什么是等效粘滞阻尼？答：振动过程的能量耗散称为阻尼。产生阻尼的原因主要有：材料的内摩擦、构件间接触面的摩擦、介质的阻力等等。当然，也包括结构中安装的各种阻尼器、耗能器。阻尼力是根据所假设的阻尼理论作用于质量上用于代替能量耗散的一种假想力。粘滞阻尼理论假定阻尼力与质量的速度成比例。粘滞阻尼理论的优点是便于求解，但其缺点是与往往实际不符，为扬长避短，按能量等效原则将实际的阻尼耗能换算成粘滞阻尼理论的相关参数，这种阻尼假设称为等效粘滞阻尼。 1.6 采用集中质量法、广义位移法（坐标法）和有限元法都可使无限自由度体系简化为有限自由度体系，它们采用的手法有何不同？答：集中质量法：将结构的分布质量按一定规则集中到结构的某个或某些位置上，认为其他地方没有质量。质量集中后，结构杆件仍具有可变形性质，称为“无重杆”。广义坐标法：在数学中常采用级数展开法求解微分方程，在结构动力分析中，也可采用相同的方法求解，这就是广义坐标法的理论依据。所假设的形状曲线数目代表在这个理想化形式中所考虑的自由度个数。考虑了质点间均匀分布质量的影响（形状函数），一般来说，

ANSYS进行有限元静力学分析

一、设计大纲概述 1、设计目的（1）熟悉有限元分析的基本原理和基本方法；（2）掌握有限元软件ANSYS的基本操作；（3）对有限元分析结果进行正确评价。 2、设计原理利用ANSYS进行有限元静力学分析。 3、设计仪器设备 1）安装windows 2000以上版本的微机； 2）ANSYS 8.0以上版本软件。 4、实验内容与步骤 1）熟悉ANSYS的界面和分析步骤； 2）掌握ANSYS前处理方法，包括平面建模、单元设置、网格划分和约束设置； 3）掌握ANSYS求解和后处理的一般方法； 4）实际应用ANSYS软件对平板结构进行有限元分析。二、题目：如图试样期尺寸为100mm*5mm*5mm,下端固定，上端受拉力10000N作用。已知该试样材料的应力-应变曲线如图所示。计算试样的位移分布。

有限元计算原理与方法..

1.有限元计算原理与方法有限元是将一个连续体结构离散成有限个单元体，这些单元体在节点处相互铰结，把荷载简化到节点上，计算在外荷载作用下各节点的位移，进而计算各单元的应力和应变。用离散体的解答近似代替原连续体解答，当单元划分得足够密时，它与真实解是接近的。 1.1. 有限元分析的基本理论有限元单元法的基本过程如下： 1.1.1.连续体的离散化首先从几何上将分析的工程结构对象离散化为一系列有限个单元组成，相邻单元之间利用单元的节点相互连接而成为一个整体。单元可采用各种类型，对于三维有限元分析，可采用四面体单元、五西体单元和六面体单元等。在Plaxis 3D Foundation 程序中，土体和桩体主要采用包含6个高斯点的15节点二次楔形体单元，该单元由水平面为6 节点的三角形单元和竖直面为四边形8节点组成的，其局部坐标下的节点和应力点分布见图3.1，图3.1 15节点楔形体单元节点和应力点分布界面单元采用包含9个高斯点的 8个成对节点四边形单元。在可能出现应力集中或应力梯度较大的地方，应适当将单元划分得密集些；

若连续体只在有限个点上被约束，则应把约束点也取为节点：若有面约束，则应把面约束简化到节点上去，以便对单元组合体施加位移边界条件，进行约束处理；若连续介质体受有集中力和分布荷载，除把集中力作用点取为节点外，应把分布荷载等效地移置到有关节点上去。最后，还应建立一个适合所有单元的总体坐标系。由此看来，有限单元法中的结构已不是原有的物体或结构物，而是同样材料的由众多单元以一定方式连接成的离散物体。因此，用有限元法计算获得的结果只是近似的，单元划分越细且又合理，计算结果精度就越高。与位移不同，应力和应变是在Gauss 积分点(或应力点)而不是在节点上计算的，而桩的内力则可通过对桩截面进行积分褥到。 1.1. 2. 单元位移插值函数的选取在有限元法中，将连续体划分成许多单元，取每个单元的若干节点的位移作为未知量，即{}[u ,v ,w ,...]e T i i i δ=，单元体内任一点的位移为{}[,,]T f u v w =。引入位移函数N (x,y,z )表示场变量在单元内的分布形态和变化规律，以便用场变量在节点上的值来描述单元内任一点的场变量。因此在单元内建立的位移模式为： {}[]{}e f N δ= （3-1）其中：12315[][,,......]N IN IN IN IN =，I 为单位矩阵。按等参元的特性，局部坐标(,,)ξηζ到整体坐标,,x y z （）的坐标转换也采用与位移插值类似的表达式。经过坐标变化后子单元与母单元(局部坐标下的规则单元)之间建立一种映射关系。不管内部单元或边界附近的单元均可选择相同的位移函数，则为它们建立单元特性矩阵的方法是相同的。因此，对于15节点楔形体单元体内各点位移在整体坐标系,,x y z （）下一般取：

弹性力学及有限元试题

弹性力学及有限元试题 (一) 问答题（20分） 1、什么是圣维南原理？举例说明怎样把它应用于工程问题的简化中。 2、什么叫做一点的应力状态？如何表示一点的应力状态（要求具体说明或表达）。 3、何谓逆解法和半逆解法？它们的理论依据是什么？ 4、什么是平面应力问题？什么是平面应变问题？分别写出弹性力学平面应力问题和平面应变问题的物理方程。 5、要保证有限元方法解答的收敛性，位移模式必须满足那些条件？ (二) (10分) 1.利用坐标变换从直角坐标的平衡方程推导极坐标下平衡方程（无体力）。 2.利用坐标变换从直角坐标下几何方程推导极坐标下几何方程。（三）已知，其他应力分量为零，求位移场。（10分）（四）设有矩形截面的悬臂粱，在自由端受有集中荷载F；体力可以不

计。试根据材料力学公式，写出弯应力σx和切应力τxy的表达式，并取挤压应力σy=0，然后证明，这些表达式满足平衡微分方程和相容方程，再说明，这些表达式是否就表示正确的解答（10分）。（五）设半平面体在直边界上受有集中力偶，单位宽度上力偶矩为M，试求应力分量（10分）。提示：单位厚度上的力偶矩M的量纲是LMT-2，应力只能是M/ρ2的形式，所以可假设应力函数由：Φ=Φ(φ). (六) 铅直平面内的正方形薄板，边长为2a，四边固定，图5—18，只受重力的作用。设μ=0，试取位移分量的表达式为用瑞利—里茨法求解（15分）。

（七）试按图示网格求解结点位移，取t =1m，μ= 0（15分）。（八）用刚度集成法求下图所示结构的整体刚度矩阵K。（10分）要求：单元刚度矩阵元素用e k形式表示；单元刚度矩阵用e K形式表 ij 示，其中e为单元号。

有限元计算结构力学大作业

SHANGHAI JIAO TONG UNIVERSITY 平面应力问题解的Matlab实现姓名: heiya168 学号: 帆哥班级：指导老师：

目录 1绪论 (1) 2平面问题的四节点四边形单元 (2) 2.1单元的构造 (2) 2.2等参变换 (5) 2.3边界条件的处理——置“1”法 (8) 3有限元分析流程 (10) 3.1程序原理和流程 (10) 3.2使用的函数 (11) 3.3文件管理 (11) 3.4数据文件格式 (11) 4算例——开方孔的矩形板拉伸分析 (13) 4.1问题的具体参数与载荷 (13) 4.2Matlab程序计算 (13) 4.3ANSYS建模计算 (15) 4.4误差分析 (17) 5总结 (18) 参考文献 (18) 附录 (19)

1绪论有限元方法(finite element method)，是求取复杂微分方程近似解的一种非常有效的工具，是现代数字化科技的一种重要基础性原理。将它用于在科学研究中，可成为探究物质客观规律的先进手段。将它应用于工程技术中，可成为工程设计和分析的可靠工具。弹性体在载荷作用下，其基本方程可写成以下的三类方程和两种边界条件。平衡方程——应力与外载荷的关系；几何方程——应变位移关系；物理方程——应力应变关系；力的边界条件；几何边界条件。应用最小位能原理，并利用上述关系，最终建立由刚度方程，节点位移和等效节点载荷所构成的求解方程。带入边界条件求解方程，就可以得出弹性力学问题的一般性解答。本次大作业基于有限元方法的基本原理，使用Matlab这一平台，针对平面应力问题，采用四节点四边形单元编写了求解单元节点位移的程序。主要内容包括：1）介绍有限元的基本原理；2）编程基本思路及流程介绍；3）程序原理及说明； 4）具体算例这四个部分。

CATIA有限元分析计算实例

CATIA有限元分析计算实例 11.1例题1 受扭矩作用的圆筒 11.1－1划分四面体网格的计算（1）进入【零部件设计】工作台启动CATIA软件。单击【开始】→【机械设计】→【零部件设计】选项，如图11－1所示，进入【零部件设计】工作台。图11－1单击【开始】→【机械设计】→【零部件设计】选项单击后弹出【新建零部件】对话框，如图11-2所示。在对话框内输入新的零件名称，在本例题中，使用默认的零件名称【Part1】。点击对话框内的【确定】按钮，关闭对话框，进入【零部件设计】工作台。（2）进入【草图绘制器】工作台在左边的模型树中单击选中【xy平面】, 如图11-3所示。单击【草图编辑器】工具栏内的【草图】按钮，如图11-4所示。这时进入【草图绘制器】工作台。图11－2【新建零部件】对话框图11－3单击选中【xy平面】（3）绘制两个同心圆草图点击【轮廓】工具栏内的【圆】按钮，如图11-5所示。在原点点击一点，作为圆草图的圆心位置，然后移动鼠标，绘制一个圆。用同样分方法再绘制一个同心圆，如图11-6所示。图11－4【草图编辑器】工具栏图11－5【轮廓】工具栏下面标注圆的尺寸。点击【约束】工具栏内的【约束】按钮，如图11-7所示。点击选择圆，就标注出圆的直径尺寸。用同样分方法标注另外一个圆的直径，如图11-8所示。

图11－6两个同心圆草图图11－7【约束】工具栏双击一个尺寸线，弹出【约束定义】对话框，如图11－9所示。在【直径】数值栏内输入100mm，点击对话框内的【确定】按钮，关闭对话框，同时圆的直径尺寸被修改为100mm。用同样的方法修改第二个圆的直径尺寸为50mm。修改尺寸后的圆如图11-10所示。图11－8标注直径尺寸的圆草图图11－9【约束定义】对话框（4）离开【草图绘制器】工作台点击【工作台】工具栏内的【退出工作台】按钮，如图11-11所示。退出【草图绘制器】工作台，进入【零部件设计】工作台。图11－10修改直径尺寸后的圆图11－11【工作台】工具栏（5）拉伸创建圆筒点击【基于草图的特征】工具栏内的【凸台】按钮，如图11-12所示。弹出【凸台定义】对话框，如图11-13所示。在【第一限制】选项组内的【长度】数值栏内输入50mm，点击对话框内的【确定】按钮，生成一个圆筒体，如图11-14所示。在左边的模型树上出现【填充器.1】元素。

有限元分析过程

有限元分析过程有限元分析过程可以分为以下三个阶段：1.建模阶段: 建模阶段是根据结构实际形状和实际工况条件建立有限元分析的计算模型——有限元模型，从而为有限元数值计算提供必要的输入数据。有限元建模的中心任务是结构离散，即划分网格。但是还是要处理许多与之相关的工作：如结构形式处理、集合模型建立、单元特性定义、单元质量检查、编号顺序以及模型边界条件的定义等。2.计算阶段: 计算阶段的任务是完成有限元方法有关的数值计算。由于这一步运算量非常大，所以这部分工作由有限元分析软件控制并在计算机上自动完成。3.后处理阶段: 它的任务是对计算输出的结果惊醒必要的处理，并按一定方式显示或打印出来，以便对结构性能的好坏或设计的合理性进行评估，并作为相应的改进或优化，这是惊醒结构有限元分析的目的所在。注意：在上述三个阶段中，建立有限元模型是整个有限分析过程的关键。首先，有限元模型为计算提供所以原始数据，这些输入数据的误差将直接决定计算结果的精度；其次，有限元模型的形式将对计算过程产生很大的影响，合理的模型既能保证计算结构的精度，又不致使计算量太大和对计算机存储容量的要求太高；再次，由于结构形状和工况条件的复杂性，要建立一个符合实际的有限元模型并非易事，它要考虑的综合因素很多，对分析人员提出了较高的要求；最后，建模所花费的时间在整个分析过程中占有相当大的比重，约占整个分析时间的70%，因此，把主要精力放在模型的建立上以及提高建模速度是缩短整个分析周期的关键。原始数据的计算模型，模型中一般包括以下三类数据：1.节点数据: 包括每个节点的编号、坐标值等；2.单元数据：a.单元编号和组成单元的节点编号；b.单元材料特性，如弹性模量、泊松比、密度等；c.单元物理特征值，如弹簧单元的刚度系数、单元厚度、曲率半径等；d.一维单元的截面特征值，如截面面积、惯性矩等；e.相关几何数据3.边界条件数据：a.位移约束数据；b.载荷条件数据；c.热边界条件数据；d.其他边界数据.建立有限元模型的一般过程：1.分析问题定义在进行有限元分析之前，首先应对结果的形状、尺寸、工况条件等进行仔细分析，只有正确掌握了分析结构的具体特征才能建立合理的几何模型。总的来说，要定义一个有限元分析问题时，应明确以下几点：a.结构类型； b.分析类型； c.分析内容； d.计算精度要求； e.模型规模； f.计算数据的大致规律2.几何模型建立几何模型是从结构实际形状中抽象出来的，并不是完全照搬结构的实际形状，而是需要根据结构的具体特征对结构进行必要的简化、变化和处理，以适应有限元分析的特点。3.单元类型选择划分网格前首先要确定采用哪种类型的单元，包括单元的形状和阶次。单元类型选择应根据结构的类型、形状特征、应力和变形特点、精度要求和硬件条件等因素综合进行考虑。4.单元特性定义有限元单元中的每一个单元除了表现出一定的外部形状外，还应具备一组计算所需的内部特征参数，这些参数用来定义结构材料的性能、描述单元本身的物理特征和其他辅助几何特征等.5.网格划分网格划分是建立有限元模型的中心工作，模型的合理性很大程度上可以通过所划分的网格形式反映出来。目前广泛采用自动或半自动网格划分方法，如在Ansys中采用的SmartSize网格划分方法就是自动划分方法。6.模型检查和处理一般来说，用自动或半自动网格划分方法划分出来的网格模型还不能立即应用于分析。由于结构和网格生成过程的复杂性，划分出来的网格或多或少存在一些问题，如网格形状较差，单元和节点编号顺序不合理等，这些都将影响有限元计算的计算精度和计算时间。7.边界条件定义在对结构进行网格划分后称为离散模型，它还不是有限元模型，只有在网格模型上定义了所需要的各类边界条件后，网格模型才能成为完整的有限元模型。计算机几何建模方法㈠.几何模型的形式1.线框模型：用组成结构的棱边表示结构形状和大小的模型称为线框模型，或线架模型。它是使用最早的几何模型，其特点是数据量少、数据结构简单、算法处理方便，模型输入可以通过定义线段端点坐标来实现。但是这种模型有很大的局限性，它的几何描述能力差，只能提供一个框架，对几何形状的理解很容易产生多义性，也不能计算结构的重量、

弹性力学与有限元法分析及实例讲解

弹性力学与有限元法分析弹性力学是固体力学的一个重要分支，是研究弹性固体在受外力作用、温度改变、边界约束或其他外界因素作用下而发生的应力、形变和位移状态的科学。有限单元法是力学、数学、物理学、计算方法、计算机技术等多种学科综合发展和结合的产物，是随着计算机技术的广泛应用而迅速发展起来的一种数值分析方法。有限元法的基本思想就是化整为零，分散分析，再集零为整。即用结构力学方法求解弹性力学问题，实质是将复杂的连续体划分为有限多个简单的单元体，单元体之间仅仅通过结点相连，实现化无限自由度问题为有限稀有度问题，将连续场函数的（偏）微分方程的求解问题转化为有限个参数的代数方程组的求解问题。有限元方法经过近半个世纪的发展，目前已经成为各种工程问题特别是结构分析问题的标准分析方法，而有限元软件也已成为现代结构设计中不可缺少的工具。有限元软件是有限元理论通向实际工程应用的桥梁，它的应用极大地提高了力学学科解决自然科学和工程实际问题的能力，进一步促进了有限元方法的发展。ANSYS软件是融结构、流体、电场、磁场、声场分析于一体的大型通用有限元分析软件，广泛用于机械制造、石油化工、航空航天、汽车交通、土木工程、造船、水利等一般工业及科学研究。 ANSYS软件的组成：（一）前处理模块该模块为用户提供了一个强大的实体建模及网格划分工具，可以方便的构造有限元模型，软件提高了100种以上的单元类型，用来模拟工程中的各种结构和材料。包括： 1.实体建模：参数化建模，布尔运算及体素库，拖拉、旋转、拷贝、蒙皮、倒角等。 2.自动网格划分，自动进行单元形态、求解精度检查及修正。 3.在集合模型上加载：点加载、分布载荷、体载荷、函数载荷。 4.可扩展的标准梁截面形状库。（二）分析计算模块该模块包括结构分析（可进行线性分析、非线性分析和高度非线性分析）、流体动力学分析、电磁场分析、声场分析、压电分析以及多物理场的耦合分析，可模拟多种物理介质的相互作用，具有灵敏度分析及优化分析能力。（三）后处理模块将计算结果以彩色等值线、梯度、矢量、粒子流、立体切片、透明及半透明等图形方式显示出来，也可以用图表、曲线形式显示或输出。由于现在只是对ANSYS工程软件有初步的了解和掌握，所以本次作业仅以（1）结构静力学分析为例，运用ANSYS软件对汽车连杆进行受力分析；（2）

悬臂梁—有限元ABAQUS线性静力学分析实例

线性静力学分析实例——以悬臂梁为例线性静力学问题是简单且常见的有限元分析类型，不涉及任何非线性（材料非线性、几何非线性、接触等），也不考虑惯性及时间相关的材料属性。在ABAQUS中，该类问题通常采用静态通用（Static，General）分析步或静态线性摄动（Static，Linear perturbation）分析步进行分析。线性静力学问题很容易求解，往往用户更关系的是计算效率和求解效率，希望在获得较高精度的前提下尽量缩短计算时间，特别是大型模型。这主要取决于网格的划分，包括种子的设置、网格控制和单元类型的选取。在一般的分析中，应尽量选用精度和效率都较高的二次四边形/六面体单元，在主要的分析部位设置较密的种子；若主要分析部位的网格没有大的扭曲，使用非协调单元（如CPS4I、C3D8I）的性价比很高。对于复杂模型，可以采用分割模型的方法划分二次四边形/六面体单元；有时分割过程过于繁琐，用户可以采用精度较高的二次三角形/四面体单元进行网格划分。悬臂梁的线性静力学分析 1.1 问题的描述一悬臂梁左端受固定约束，右端自由，结构尺寸如图1-1所示，求梁受载后的Mises应力、位移分布。 ν 材料性质：弹性模量3 = E=，泊松比3.0 2e 均布载荷：F=103N 图1-1 悬臂梁受均布载荷图 1.2 启动ABAQUS 启动ABAQUS有两种方法，用户可以任选一种。（1）在Windows操作系统中单击“开始”--“程序”--ABAQUS 6.10 --

ABAQUS/CAE。（2）在操作系统的DOS窗口中输入命令：abaqus cae。启动ABAQUS/CAE后，在出现的Start Section（开始任务）对话框中选择Create Model Database。 1.3 创建部件在ABAQUS/CAE顶部的环境栏中，可以看到模块列表：Module：Part，这表示当前处在Part（部件）模块，在这个模块中可以定义模型各部分的几何形体。可以参照下面步骤创建悬臂梁的几何模型。（1）创建部件。对于如图1-1所示的悬臂梁模型，可以先画出梁结构的二维截面（矩形），再通过拉伸得到。单击左侧工具区中的（Create Part）按钮，或者在主菜单里面选择Part--Create，弹出如图1-2所示的Create Part对话框。图1-2 Create Part对话框在Name（部件名称）后面输入Beam，Modeling Space（模型所在空间）设

有限元分析的基本步骤

一个典型的ANSYS分析过程可分为以下6个步骤： 1定义参数 2创建几何模型 3划分网格 4加载数据 5求解 6结果分析 1定义参数 1.1指定工程名和分析标题启动ANSYS软件，选择File→Change Jobname命令选择File→Change Title菜单命令 1.2定义单位 (2) 设置计算类型 ANSYS Main Menu: Preference→Material Props →Material Models →Structural →OK (3) 定义分析类型 ANSYS Main Menu: Preprocessor →Loads →Analysis Type →New Analysis→STATIC →OK 1.3定义单元类型选择Main Menu→Preprocessor→Element Type→Add/Edit/Delete命令单击[Options]按钮，在[Element behavior]下拉列表中选择[Plane strs w/thk]选项，单击确定 1.4定义单元常数在ANSYS程序主界面中选择Main Menu→Preprocessor→Real Constants→Add/Edit/Delete命令单击[Add]按钮，进行下一个[Choose Element Type]对话框 1.5定义材料参数在ANSYS程序主界面，选择Main Menu→Preprocessor→Material Props→Material Models命令 (1)选择对话框右侧Structural→Linear→Elastic→Isotropic命令，并单击[Isotropic]选项，接着弹出如下所示[Linear Isotropic Properties for Material Number 1]对话框。在[EX]文本框中输入弹性模量“200000”，在[PRXY]文本框中输入泊松比“0.3”，单击OK 2创建几何模型在ANSYS程序主界面，选择Main Menu→Preprocessor→Modeling→Creat→Areas→Rectangle →By 2Corners命令选择Main Menu→Preprocessor→Modeling→Creat→Areas→Circle→Solid Circle命令 3网格划分(之前一定要进行材料的定义和分配) 选择Main Menu→Preprocessor→Modeling→Operate→Booleans→Subtract→Arears Circle命令选择Main Menu→Preprocessor→Meshing→Mesh→Areas→Free命令，弹出实体选择对话框，单击[Pick All]按钮，得到如下所示网格 4加载数据 (1)选择Main Menu→Preprocessor→Loads→Define Loads→Apply→Structural→Displacement→On Lines命令，出现如下所示对话框，选择约束[ALL DOF]选项，并设置[Displacement value]为0，单击OK。