谭继锦有限元法课件之八 5概述薄板单元

有限元基础知识归纳

有限元知识点归纳 1.、有限元解的特点、原因? 答:有限元解一般偏小,即位移解下限性 原因:单元原是连续体的一部分,具有无限多个自由度。在假定了单元的位移函数后,自由度限制为只有以节点位移表示的有限自由度,即位移函数对单元的变形进行了约束和限制,使单元的刚度较实际连续体加强了,因此,连续体的整体刚度随之增加,离散后的刚度较实际的刚度K为大,因此求得的位移近似解总体上将小于精确解。 2、形函数收敛准则(写出某种单元的形函数,并讨论收敛性)P49 (1)在节点i处N i=1,其它节点N i=0; (2)在单元之间,必须使由其定义的未知量连续; (3)应包含完全一次多项式; (4)应满足∑Ni=1 以上条件是使单元满足收敛条件所必须得。可以推证,由满足以上条件的形函数所建单元是完备协调的单元,所以一定是收敛的。 4、等参元的概念、特点、用时注意什么?(王勖成P131) 答:等参元—为了将局部坐标中几何形状规则的单元转换成总体(笛卡尔)坐标中的几何形状扭曲的单元,以满足对一般形状求解域进行离散化的需要,必须建立一个坐标变换。即: 为建立上述的变换,最方便的方法是将上式表示成插值函数的形式,即: 其中m是用以进行坐标变换的单元节点数,xi,yi,zi是这些结点在总体(笛卡尔)坐标内的坐标值,Ni’称为形状函数,实际上它也是局部坐标表示的插值函数。称前者为母单元,后者为子单元。 还可以看到坐标变换关系式和函数插值表示式:在形式上是相同的。如果坐标变换和函数插值采用相同的结点,并且采用相同的插值函数,即m=n,Ni’=Ni,则称这种变换为等参变换。 5、单元离散?P42 答:离散化既是将连续体用假想的线或面分割成有限个部分,各部分之间用有限个点相连。每个部分称为一个单元,连接点称为结点。对于平面问题,最简单、最常用的离散方式是将其分解成有限个三角形单元,单元之间在三角形顶点上相连。这种单元称为常应变三角形单元。常用的单元离散有三节点三角形单元、六节点三角形单元、四节点四边形单元、八节点四边形单元以及等参元。 6、数值积分,阶次选择的基本要求? 答:通常是选用高斯积分 积分阶次的选择—采用数值积分代替精确积分时,积分阶数的选取应适当,因为它直接影响计算精度,计算工作量。选择时主要从两方面考虑。一是要保证积分的精度,不损失收敛性;二是要避免引起结构总刚度矩阵的奇异性,导致计算的失败。

有限元知识点汇总

有限元知识点汇总 第一章 1、何为有限元法?其基本思想是什么? 》有限元法是一种基于变分法而发展起来的求解微分方程的数值计算方法。 》基本思想:化整为零,化零为整 2、为什么说有限元法是近似的方法,体现在哪里? 》有限元法的基本思想是几何离散和分片插值; 》用离散单元的组合来逼近原始结构,体现了几何上的近似;用近似函数逼近未知量在单元内的真实解,体现了数学上的近似;利用与问题的等效的变分原理建立有限元基本方程,又体现了明确的物理背景。 3、单元、节点的概念? 》单元:把参数单元划分成网格,这些网格就称为单元。 》节点:网格间相互连接的点称为节点。 4、有限元法分析过程可归纳为几个步骤? 》3大步骤;——结构离散化;——单元分析;——整体分析。 5、有限元方法分几种?本课程讲授的是哪一种? 》有限元方法分3种;——位移法、力法、混合法。 》本课程讲授的:位移法 6、弹性力学的基本变量是什么?何为几何方程、物理方程及虚功方程?弹性矩阵的特点?》弹性力学的基本变量是——{外力、应力、应变、位移} 》几何方程——{描述弹性体应变分量与位移分量之间关系的方程} 》物理方程——{描述应力分量与应变分量之间的关系} 》虚功方程——{描述内力和外力的关系的方程} 》弹性矩阵特点——{ } 7、何为平面应力问题和平面应变问题? 》平面应力问题——{满足(1)几何条件——所研究的是一根很薄的等厚度薄板,即一个方向上的几何尺寸远远小于其余两个面上的几何尺寸;(2)载荷条件——作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布,而在两板面上无外力作用} 》平面应变问题——{满足(1)几何条件——所研究的是长柱体,即长度方向的尺寸远远大于横截面的尺寸,且横截面沿长度方向不变;(2)载荷条件——作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力} 第二章 7、形函数的特点? 》1形函数Ni再节点i处等于1,在其他节点上的值等于0,对于Nj、Nm也有同样的性质。》2在单元内任一点的各形函数之和等于1,即Ni+Nj+Nm=1 8、单元刚度矩阵的性质? 》1 K^e中每个元素都有明确的物理意义,每个元素都是一个刚度系数,他是单位节点位移分量所引起的节点力分量 》2 k^e是对称矩阵,具有对称性。 》3 K^e的每一行或每一列元素之和为零,是奇异矩阵

有限元期末复习题资料

1、弹性力学与材料力学主要不同在于:研究方法。 2、利用Ansys 进行结构分析时,结果文件是什么文件:jobname.rst文件。 3、在Ansys单元库中,Plane42属于结构实体单元。 4、在一个分析中可能有多个材料特性组,Ansys通过独特的( C )来识别每个材料的特性组。 A. 特性 B. 说明 C. 参考号 D. 方法 5、载荷包括所有边界条件以及外部或者内部的作用效应,下列不属于Ansys载荷的是( D )。 A. DOF约束 B. 力 C. 体载荷 D. 应力 【解析】:应力是结果,不是条件。 6、( B )什么要求面或者体有规则的形状,即必须满足一定的准则。 A. 自由网格 B. 映射网格 C. Sweep网格 D. 其他 7、什么样的载荷独立于有限元网格,即可以改变单元网格而不影响施加的载荷( C )。 A. 阶跃载荷 B. 有限元模型载荷 C. 实体模型载荷 D. 斜坡载荷 8、有限元法首先把求解出的解是( D ),单元应变和应力都可以由它来求得。 A. 节点坐标 B. 节点自由度 C. 节点载荷 D.节点位移 9、下列不属于Ansys产品当中求解联立方程的方法是( C )。 A. 稀疏矩阵直接解法 B. 直接解法 C. 变分法 D. 雅可比共轭梯度法 10、下列不属于/post1显示的图形类别的是( B )。 A. 等直线图 B. 灰度图 C. 形状变形图 D. 矢量图 11、对二维桁架进行强度校核时,选择的单元类型是( C )。 A. plane82 B. Beam3 C. Link2DSPrl D. Shell63 δ 12、δ为板的厚度,b为长度的最小值,当满足8/1 < -b < 1/100 5/1 / 1/80- 时,这样的板属于( B )。 A. 薄膜 B. 薄板 C. 厚板 D. 壳 13、下列哪个布尔运算的结果是由每个初始输入的图元的共同部分形成的新图元( A )

有限元分析及其应用思考题附答案2012

有限元分析及其应用-2010 思考题: 1、有限元法的基本思想是什么?有限元法的基本步骤有那些?其中“离散”的含义是什 么?是如何将无限自由度问题转化为有限自由度问题的? 答:基本思想:几何离散和分片插值。 基本步骤:结构离散、单元分析和整体分析。 离散的含义:用假想的线或面将连续物体分割成由有限个单元组成的集合,且单元之间仅在节点处连接,单元之间的作用仅由节点传递。当单元趋近无限小,节点无限多,则这种离散结构将趋近于实际的连续结构。 2、有限元法与经典的差分法、里兹法有何区别? 区别:差分法:均匀离散求解域,差分代替微分,要求规则边界,几何形状复杂精度较低; 里兹法:根据描述问题的微分方程和相应的定解构造等价的泛函表达式,求得近似解; 有限元:基于变分法,采用分片近似进而逼近总体的求解微分方程的数值计算方法。 3、一根单位长度重量为q的悬挂直杆,上端固定,下端受垂直向下的外力P,试 1)建立其受拉伸的微分方程及边界条件; 2)构造其泛函形式; 3)基于有限元基本思想和泛函求极值构造其有限元的计算格式(即最小势能原理)。4、以简单实例为对象,分别按虚功原理和变分原理导出有限元法的基本格式(单元刚度矩 阵)。 5、什么是节点力和节点载荷?两者有何区别? 答:节点力:单元与单元之间通过节点相互作用 节点载荷:作用于节点上的外载 6、单元刚度矩阵和整体刚度矩阵各有何特点?其中每个矩阵元素的物理意义是什么(按自 由度和节点解释)? 答:单元刚度矩阵:对称性、奇异性、主对角线恒为正 整体刚度矩阵:对称性、奇异性、主对角线恒为正、稀疏性、带状性。 Kij,表示j节点产生单位位移、其他节点位移为零时作用i节点的力,节点力等于节点位移与单元刚度元素乘积之和。 7、单元的形函数具有什么特点?有哪些性质? 答:形函数的特点:Ni为x,y的坐标函数,与位移函数有相同的阶次。 形函数Ni在i节点的值为1,而在其他节点上的值为0; 单元内任一点的形函数之和恒等于1; 形函数的值在0~1间变化。 8、描述弹性体的基本变量是什么?基本方程有哪些组成? 答:基本变量:外力、应力、应变、位移 基本方程:平衡方程、几何方程、物理方程、几何条件 9、何谓应力、应变、位移的概念?应力与强度是什么关系? 答:应力:lim△Q/△A=S △A→0 应变:物体形状的改变 位移:弹性体内质点位置的变化 10、问题的微分方程提法、等效积分提法和泛函变分提法之间有何关系?何谓“强形 式”?何谓“弱形式”,两者有何区别?建立弱形式的关键步骤是什么?

华南理工大学有限元考试试题

华南理工大学广州学院有限单元法期末试题大纲 一、选择题: 1 在加权余量法中,若简单地利用近似解的试探函数序列作为权函数,这类方 法称为________________。 (A)配点法(B)子域法(C)伽辽金法 2 等参变换是指单元坐标变换和函数插值采用______的结点和______的插值 函数。 (A)不相同,不相同(B)相同,相同(C)相同,不相同(D)不相同,相同3 有限元位移模式中,广义坐标的个数应与___________相等。 (A)单元结点个数(B)单元结点自由度数(C)场变量个数 4 采用位移元计算得到应力近似解与精确解相比较,一般___________。(A)近似解总小于精确解(B)近似解总大于精确解(C)近似解在精确解上下震荡(D)没有规律 5 如果出现在泛函中场函数的最高阶导数是m阶,单元的完备性是指试探函数 必须至少是______完全多项式。 (A)m-1次(B)m次(C)2m-1次 6 与高斯消去法相比,高斯约当消去法将系数矩阵化成了_________形式,因 此,不用进行回代计算。 (A)上三角矩阵(B)下三角矩阵(C)对角矩阵 7 对称荷载在对称面上引起的________________分量为零。 (A)对称应力(B)反对称应力(C)对称位移(D)反对称位移 8 对分析物体划分好单元后,__________会对刚度矩阵的半带宽产生影响。(A)单元编号(B)单元组集次序(C)结点编号 9 n个积分点的高斯积分的精度可达到______阶。 (A)n-1 (B)n(C)2n-1 (D)2n 10 引入位移边界条件是为了消除有限元整体刚度矩阵K的__________。(A)对称性(B)稀疏性(C)奇异性 C B B C B C D C C C 二、填空题:(课本···黑色字体)····仿题 1、有限元网格划分的过程中应注意:网格数目、网格疏密、单元阶次、网格质量

有限元法基本原理与应用

有限元法基本原理与应用 班级机械2081 姓名方志平 指导老师钟相强 摘要:有限元法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。 关键词:有限元法;变分原理;加权余量法;函数。 Abstract:Finite element method is based on the variational principle and the weighted residual method, the basic idea is to solve the computational domain is divided into a finite number of non-overlapping units, each unit, select some appropriate function for solving the interpolation node points as , the differential variables rewritten or its derivative by the variable value of the selected node interpolation functions consisting of linear expressions, by means of variational principle or weighted residual method, the discrete differential equations to solve. Different forms of weight functions and interpolation functions, it constitutes a different finite element method. Keywords:Finite element method; variational principle; weighted residual method; function。 引言 有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中,常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同,有限元方法也分为多种计算格式。从权函数的选择来说,有配置法、矩量法、最小二乘法和伽辽金法,从计算单元网格的形状来划分,有三角形网格、四边形网格和多边形网格,从插值函数的精度来划分,又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数,伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数;最小二乘法是令权函数等于余量本身,而内积的极小值则为对代求系数的平方误差最小;在配置法中,先在计

ANSYS 有限元分析基本流程

第一章实体建模 第一节基本知识 建模在ANSYS系统中包括广义与狭义两层含义,广义模型包括实体模型和在载荷与边界条件下的有限元模型,狭义则仅仅指建立的实体模型与有限元模型。建模的最终目的是获得正确的有限元网格模型,保证网格具有合理的单元形状,单元大小密度分布合理,以便施加边界条件和载荷,保证变形后仍具有合理的单元形状,场量分布描述清晰等。 一、实体造型简介 1.建立实体模型的两种途径 ①利用ANSYS自带的实体建模功能创建实体建模: ②利用ANSYS与其他软件接口导入其他二维或三维软件所建立的实体模型。 2.实体建模的三种方式 (1)自底向上的实体建模 由建立最低图元对象的点到最高图元对象的体,即先定义实体各顶点的关键点,再通过关键点连成线,然后由线组合成面,最后由面组合成体。 (2)自顶向下的实体建模 直接建立最高图元对象,其对应的较低图元面、线和关键点同时被创建。 (3)混合法自底向上和自顶向下的实体建模 可根据个人习惯采用混合法建模,但应该考虑要获得什么样的有限元模型,即在网格划分时采用自由网格划分或映射网格划分。自由网格划分时,实体模型的建立比较1e单,只要所有的面或体能接合成一体就可以:映射网格划分时,平面结构一定要四边形或三边形的面相接而成。 二、ANSYS的坐标系 ANSYS为用户提供了以下几种坐标系,每种都有其特定的用途。 ①全局坐标系与局部坐标系:用于定位几何对象(如节点、关键点等)的空间位置。 ②显示坐标系:定义了列出或显示几何对象的系统。 ③节点坐标系:定义每个节点的自由度方向和节点结果数据的方向。 ④单元坐标系:确定材料特性主轴和单元结果数据的方向。 1.全局坐标系 全局坐标系和局部坐标系是用来定位几何体。在默认状态下,建模操作时使用的坐标系是全局坐标系即笛卡尔坐标系。总体坐标系是一个绝对的参考系。ANSYS提供了4种全局坐标系:笛卡尔坐标系、柱坐标系、球坐标系、Y-柱坐标系。4种全局坐标系有相同的原点,且遵循右手定则,它们的坐标系识别号分别为:0是笛卡尔坐标系(cartesian),1是柱坐标系 (Cyliadrical),2是球坐标系(Spherical),5是Y-柱坐标系(Y-aylindrical),如图2-1所示。

有限元复习题答案

1、何为有限元法?其基本思想是什么? 有限元法是一种基于变分法而发展起来的求解微分方程的数值计算方法,该方法以计算机为手段,采用分片近似,进而逼近整体的研究思想求解物理问题。 基本思想是化整为零集零为整。 2、为什么说有限元法是近似的方法,体现在哪里? 有两点:用离散单元的组合体来逼近原始结构,体现了几何上的近似;而用近似函数逼近未知变量在单元内的真实解,体现了数学上的近似。 3、单元、节点的概念? 节点:表达实际结构几何对象之间相互连接方式的概念 单元:网格划分中的每一个小部分称为单元,网格间相互联结点称为节点 4、有限元法分析过程可归纳为几个步骤? 结构离散化、单元分析、整体分析 5、有限元方法分几种?本课程讲授的是哪一种? 位移法、力法、混合法本课程讲授位移法 6、弹性力学的基本变量是什么?何为几何方程、物理方程及虚功方程?弹性矩阵的特点? 弹性力学变量:外力、应力、应变和位移。 描述弹性体应变分量与位移分量之间的方程称为几何方程;物理方程描述应力分量与应变分量之间的关系;弹性体上外力在虚位移发生过程中所做的虚功与储存在弹性体内的需应变能相等。 弹性矩阵由材料的弹性模量和泊松比确定,与坐标位置无关。 7、何为平面应力问题和平面应变问题? 平面应力问题:在结构上满足a几何条件:研究对象是等厚度薄板。b载荷条件:作用于薄板上的载荷平行于板平面且沿厚度方向均匀分布,而在两板面无外力作用。 平面应变问题:满足a几何条件:长柱体,即长度方向的尺寸远远大于横截面的尺寸,且横截面沿长度方向不变。b载荷条件:作用于长柱体结构上的载荷平行于横截面且沿纵向方向均匀分布,两端面不受力两条件的弹性力学问题。 1、何为结构的离散化?离散化的目的?何为有限元模型? ①离散化:把连续的结构看成由有限个单元组成的集合体。②目的:建立有限元计算模型③通常把由节点,单元及相应的节点载荷和节点约束构成的模型称为有限元模型2、结构离散化时,划分单元数目的多少以及疏密分布,将直接影响到什么?确定单元数量的原则?通常如何设置节点?

有限元分析基本理论问答 基础理论知识

1. 诉述有限元法的定义 答:有限元法是近似求解一般连续场问题的数值方法 2. 有限元法的基本思想是什么 答:首先,将表示结构的连续离散为若干个子域,单元之间通过其边界上的节点连接成组合体。其次,用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。 3. 有限元法的分类和基本步骤有哪些 答:分类:位移法、力法、混合法;步骤:结构的离散化,单元分析,单元集成,引入约束条件,求解线性方程组,得出节点位移。 4. 有限元法有哪些优缺点 答:优点:有限元法可以模拟各种几何形状复杂的结构,得出其近似解;通过计算机程序,可以广泛地应用于各种场合;可以从其他CAD软件中导入建好的模型;数学处理比较方便,对复杂形状的结构也能适用;有限元法和优化设计方法相结合,以便发挥各自的优点。 缺点:有限元计算,尤其是复杂问题的分析计算,所耗费的计算时间、内存和磁盘空间等计算资源是相当惊人的。对无限求解域问题没有较好的处理办法。尽管现有的有限元软件多数使用了网络自适应技术,但在具体应用时,采用什么类型的单元、多大的网络密度等都要完全依赖适用者的经验。 5. ?梁单元和平面钢架结构单元的自由度由什么确定 答:每个节点上有几个节点位移分量,就称每个节点有几个自由度 6. ?简述单元刚度矩阵的性质和矩阵元素的物理意义 答:单元刚度矩阵是描述单元节点力和节点位移之间关系的矩阵 单元刚度矩阵中元素aml的物理意义为单元第L个节点位移分量等于1,其他节点位移分量等于0时,对应的第m个节点力分量。 7. 有限元法基本方程中的每一项的意义是什么 答:整个结构的节点载荷列阵(外载荷、约束力),整个结构的节点位移列阵,结构的整体刚度矩阵,又称总刚度矩阵。 8. 位移边界条件和载荷边界条件的意义是什么 答:由于刚度矩阵的线性相关性不能得到解,从而引入边界条件。 9. ?简述整体刚度矩阵的性质和特点 答:对称性;奇异性;稀疏性;对角线上的元素恒为正。 11. 简述整体坐标的概念 答:单元刚度矩阵的坐标变换式把平面刚架的所有单元在局部坐标系X’Y’Z’下的单元刚度矩阵变换到一个统一的坐标系xOy下,这个统一的坐标系xOy称为整体坐标系。 13. 简述平面钢架问题有限元法的基本过程 答:力学模型的确定,结构的离散化,计算载荷的等效节点力,计算各单元的刚度矩阵,组集整体刚度矩阵,施加边界约束条件,求解降价的有限元基本方程,求解单元应力,计算结果的输出。 14. 弹性力学的基本假设是什么。 答:连续性假定,弹性假定,均匀性和各向同性假定,小变形假定,无初应力假定。 15.弹性力学和材料力学相比,其研究方法和对象有什么不同。 答:研究对象:材料力学主要研究杆件,如柱体、梁和轴,在拉压、剪切、弯曲和扭转等作用下的应力、形变和位移。弹性力学研究各种形状的弹性体,除杆件外,还研究平面体、空间体,板和壳等。因此,弹性力学的研究对象要广泛得多。研究方法:弹性力学和材料力学

有限元复习重点

●有限元起源于20世纪50年代中期航空工程中飞机结构的矩阵分析。 ●有限元基本思想:在力学模型上将一个原来连续的物体离散成为有限个具有一定大小的单元,这些单元仅在有限个节点上相连接,并在节点上引进等效力以代替实际作用于单元上的外力。对于每个单元,根据分块近似的思想,选择一种简单的函数来表示单元内位移的分布规律,并按弹性理论中的能量原理(或用变分原理)建立单元节点力和节点位移之间的关系。最后,把所有单元的这种关系式集合起来,就得到一组以节点位移为未知量的代数方程组,解这些方程组就可以求出物体上有限个离散节点上的位移。 “一分一合”,化整为零,集零为整,把复杂的结构看成由有限个单元组成的整体。 ●单元、节点、边界:采用8节点四边形等参数单元把受力体划分成网格,这些网格称为单元;网格间互相连接的点称为节点;网格与网格的交界线称为边界。节点数和单元数目是有限的。 ●有限元法的优点:(1)理论基础简明,物理概念清晰,且可在不同的水平上建立起对该法的理解。(2) 具有灵活性和适用性,应用范围极为广泛。(3) 该法在具体推导运算中,广泛采用了矩阵方法,便于实现程序设计的自动化。 ●有限单元法分为三类:位移法(以节点位移为基本未知量)、力法(以节点力为基本未知量)和混合法(一部分以节点位移,另一部分以节点力作为基本未知量)。 ●有限元法分析计算的基本步骤可归纳如以下五点。1.结构的离散化(将某个机械结构划分为由各种单元组成的计算模型)在平面问题用三角形、矩形或任意四边形单元。在空间问题用四面体、长方体或任意六面体单元2.单元分析①选择位移模式(位移模式是表示单元内任意点的位移随位置变化的函数式,由于所采用的函数是一种近似的试函数,一般不能精确地反映单元中真实的位移分布)位移模式或位移函数:i n i i a y φ∑=②建立单元刚度方程e e e F k =δ,e 为单元编号;e δ为单元的节点位移向量;e F 为单元的节 点力向量 ;e k 为单元刚度矩阵.③计算等效节点力:用等效的节点力来代替所有作用在单元上的力。3.整体分析:整体的有限元方程F K =δ。K 为整体结构的刚度矩阵;δ为整体节点位移向量;F 为整体载荷向量。4.求解方程,得出节点位移5.由节点位移计算单元的应变与应力 ●有限元中得一个基本近似性是几何近似性 ●有限元中的变量:应力、应变、变形。基本方程有:平衡方程、物理方程、几何方程。边界条件:力边界、位移边界。 ●弹性力学的任务是分析弹性体在受外力作用并处于平衡状态下产生的应 力、应变和位移状态及其相互关系等。 ●外力:体力(分布在物体体积内的力---重力、惯性力、电磁力)、面力(分布在物体表面上的力---流体压力、接触力、风力) ●应力:物体受外力的作用,或由于温度有所改变,其内部将发生内力。

有限元法及其在工程中的应用

机械与汽车学院 曹国强 主要内容: 1、有限元法的基本思想。 2、结构力学模型的简化和结构离散化。 3、有限元法的实施过程。 一、有限元法的基本思想 有限元法是随着计算机的发展而发展起来的一种有效的数值方法。其基本思想是:将连续的结构分割成数目有限的小单元体(称为单元),这些小单元体彼此之间只在数目有限的指定点(称为节点)上相互连接。用这些小单元体组成的集合体来代替原来的连续结构。再把每个小单元体上实际作用的外载荷按弹性力学中的虚功等效原理分配到单元的节点上,构成等效节点力,并按结构实际约束情况决定受约束节点的约束。这一过程称为结构的离散化。其次,对每个小单元体选择一个简单的函数来近似地表示其位移分量的分布规律,并按弹性力学中的变分原理建立起单元节点力和节点位移之间的关系(单元刚度方程),最后,把全部单元的节点力和节点位移之间的关系组集起来,就得到了一组以结构节点位移为未知量的代数方程组(总体刚度方程),同时考虑结构的约束情况,消去那些结构节点位移为零的方程,再由最后的代数方程组就可求得结构上有限个离散节点的各位移分量。求得了结构上各节点的位移分量之后,即可按单元的几何方程和物理方程求得各单元的应变和应力分量。 有限元法的实质就是把具有无限个自由度的连续体,理想化为有限个自由度的单元的集合体,使问题简化为适合于数值解法的结构型问题。 经典解法(解析法)与有限元法的区别 解析法 { } 建立一个描述连续体性质的偏微分方程组 有限元解法 连续体 数目增加到∞ 大小趋于0 微元 有限元 离散化 (单元分析)集合 总体分析 求得近似解

二、结构力学模型的简化和结构离散化 (一)结构力学模型的简化 用有限元法研究实际工程结构问题时,首先要从工程实际问题中抽象出力学模型,即要对实际问题的边界条件、约束条件和外载荷进行简化,这种简化应尽可能地反映实际情况,不至于使简化后的解答与实际差别过大,但也不要带来计算上的过分复杂,在力学模型的简化过程中,必须判断实际结构的问题类型,是二维问题还是三维问题。如果是平面问题,是平面应力问题,还是平面应变问题。同时还要搞清楚结构是否对称,外载荷大小和作用位置,结构的几何尺寸和力学参数(弹性模量E、波松比μ等)。 (二)结构的离散化 将已经简化好的结构力学模型划分成只在一些节点连续的有限个单元,把每个单元看成是一个连续的小单元体,各单元之间只在一些点上互相联结,这些点称作节点,每个单元体称为一个单元。用只在节点处连接的单元的集合体代替原来的连续结构,把外载荷按虚功等效原理移置到有关受载的节点上,构成节点载荷,把连续结构进行这样分割的过程称为结构的离散化。现举例说明。 设一平面薄板,中间有一个园孔,其左端固定,右端受面力载荷q,试对其进行有限元分割和力学模型简化。

有限元方法理论及其应用

1 课程论文:弹性力学有限元位移法原理(30分) 撰写一篇论文,对有限元位移法的原理作一般性概括和论述。要求论文论及但不限于下列内容:1)弹性力学有限元位移法的基本思想和数学、力学基础;2)有限元法求解的原理和过程,推导计算列式;对基本概念和矩阵符号进行解释和讨论;3)等参单元的概念、原理和应用。 1.1 对一维杆单元有限元形式的理解 我对此提出了几点疑问: 1)为什么边界条件u1=0,就要划去刚度矩阵[K]中对应的行列再解方程? 2)为什么刚度矩阵[K]会奇异? 3)为什么平衡方程本身是矛盾的,而加上边界条件u1=0之后就能解出一 个唯一的近似解? 4)为什么刚度矩阵[K]是对称的? 下面我谈谈自己的理解:节点平衡方程是在u1不定的前提下,假设单元内位移都是线性变化推导出来的,由此u1相当于一个不确定的定值约束,再加上中间两个节点的连续性要求,系统实际上只有三个独立的自由度(广义坐标)。 对于第一个问题,其实刚度矩阵[K]中的元素不是一成不变的,相反它是伴随边界条件动态变化的。当u1=0时由刚度矩阵的推导过程可以知道,刚度矩阵的第一行和第一列都会变为0,所以此时第一行和第一列对于求解方程是没有作用的。 对于第二个问题,由于系统自由度(广义坐标)只有三个,而我们的方程却列出

了四个,显然

这四个方程不可能线性无关,所以刚度矩阵奇异。 对于第三个问题,首先我们应该明确方程区别于等式,虽然左右两边都是用“=”连接,但是方程只在特殊条件下取得定解。由于平衡方程是在没有约束的条件下推导出来的,显然它不可能满足等式要求。宏观上看,系统在没有外部约束,而又施加有外力,显然系统会产生加速度而绝不会平衡。所以平衡方程本身是矛盾的。而加上边界条件之后,不但满足了平衡的前提,还改变了矩阵的结构和性质,所以有解。但是,由于我们提前假设了位移线性变化,相当于人为对单元施加了额外约束,让位移按照我们假设的规律变化,所以得到的解是过刚的近似解。但对于方程本身而言是精确解。 对于第四个问题,其力学的作用机理类似于作用力与反作用力,由于刚度矩阵不表征方向,所以其大小是相等的。 1.2 有限元法的思想 有限元法是求解连续介质力学问题的数值方法,更一般意义是一种分析结构问题和连续场数学物理问题的数值方法。 有限元法的基本思想是离散化和分片插值。 即把连续的几何机构离散成有限个单元,并在每一个单元中设定有限个节点,从而将连续体看作仅在节点处相连接的一组单元的集合体,同时选定场函数的节点值作为基本未知量并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律,再建立用于求解节点未知量的有限元方程组,从而将一个连续域中的无限自由度问题转化为离散域中的有限自由度问题。 求解得到节点值后就可以通过设定的插值函数确定单元上以至个集合体上的场函数。对每个单元,选取适当的插值函数,使得该函数在子域内部、在子域分界面上以及子域与外界面上都满足一定的条件。单元组合体在已知外载荷作用下处于平衡状态时,列出一系列以节点、位移为未知量的线性方程组,利用计算机解出节点位移后,再用弹性力学的有关公式,计算出各单元的应力、应变,当各单元小到一定程度,那么它就代表连续体各处的真实情况。

有限元法基础试题

有限元法基础试题(A ) 一、填空题(5×2分) 1.1单元刚度矩阵e T k B DBd Ω = Ω? 中,矩阵B 为__________,矩阵D 为___________。 1.2边界条件通常有两类。通常发生在位置完全固定不能转动的情况为_______边界,具体指定有限的非零值位移的情况,如支撑的下沉,称为_______边界。 1.3内部微元体上外力总虚功: ()(),,,,e x x xy y bx xy x y y by d W F u F v dxdy δστδτσδ??=+++++??+(),,,,x x y y xy y x u v u u dxdy σδσδτδδ??+++??的表达式中,第一项为____________________的虚功,第二项为____________________的虚功。 1.4弹簧单元的位移函数1N +2N =_________。 1.5 ij k 数学表达式:令j d =_____,k d =_____,k j ≠,则力i ij F k =。 二、判断题(5×2分) 2.1位移函数的假设合理与否将直接影响到有限元分析的计算精度、效率和可靠性。( ) 2.2变形体虚功原理适用于一切结构(一维杆系、二维板、三位块体)、适用于任何力学行为的材料(线性和非线性),是变形体力学的普遍原理。 ( ) 2.3变形体虚功原理要求力系平衡,要求虚位移协调,是在“平衡、协调”前提下功的恒等关系。 ( ) 2.4常应变三角单元中变形矩阵是x 或y 的函数。 ( ) 2.5 对称单元中变形矩阵是x 或y 的函数。 ( ) 三、简答题(26分) 3.1列举有限元法的优点。(8分) 3.2写出有限单元法的分析过程。(8分) 3.3列出3种普通的有限元单元类型。(6分) 3.4简要阐述变形体虚位移原理。(4分) 四、计算题(54分) 4.1对于下图所示的弹簧组合,单元①的弹簧常数为10000N/m ,单元②的弹簧常数为20000N/m ,单元③的弹簧常数为10000N/m ,确定各节点位移、反力以及单元②的单元力。(10分) 4.2对于如图所示的杆组装,弹性模量E 为10GPa ,杆单元长L 均为2m ,横截面面积A 均为2×10-4m 2,弹簧常数为2000kN/m ,所受荷载如图。采用直接刚度法确定节点位移、作用力和单元②的应力。(10分)

有限元法课后习题复习资料

1、有限元是近似求解一般连续场问题的数值方法 2、有限元法将连续的求解域离散为若干个子域,得到有限个单元,单元和单元之间用节点连接 3、直梁在外力的作用下,横截面的内力有剪力和弯矩两个. 4、平面刚架结构在外力的作用下,横截面上的内力有轴力、剪力、弯矩 . 5、进行直梁有限元分析,平面刚架单元上每个节点的节点位移为挠度和转角 6、平面刚架有限元分析,节点位移有轴向位移、横向位移、转角。 7、在弹性和小变形下,节点力和节点位移关系是线性关系。 8、弹性力学问题的方程个数有15个,未知量个数有15个。 9、弹性力学平面问题方程个数有8,未知数8个。 10、几何方程是研究应变和位移之间关系的方程 11、物理方程是描述应力和应变关系的方程 12、平衡方程反映了应力和体力之间关系的 13、把经过物体内任意一点各个截面上的应力状况叫做一点的应力状态 14、9形函数在单元上节点上的值,具有本点为_1_.它点为零的性质,并且在三角形单元的任一节点上,三个行函数之和为_1_ 15、形函数是_三角形_单元内部坐标的_线性_函数,他反映了单元的_位移_状态

16、在进行节点编号时,同一单元的相邻节点的号码差尽量小 . 17、三角形单元的位移模式为_线性位移模式 18、矩形单元的位移模式为双线性位移模式_ 19、在选择多项式位移模式的阶次时,要求_所选的位移模式应该与局部坐标系的方位无关的性质为几何_各向同性 20、单元刚度矩阵描述了_节点力_和_节点位移之间的关系 21、矩形单元边界上位移是连续变化的 1. 诉述有限元法的定义 答:有限元法是近似求解一般连续场问题的数值方法 2. 有限元法的基本思想是什么 答:首先,将表示结构的连续离散为若干个子域,单元之间通过其边界上的节点连接成组合体。其次,用每个单元内所假设的近似函数分片地表示求解域内待求的未知厂变量。 3. 有限元法的分类和基本步骤有哪些 答:分类:位移法、力法、混合法;步骤:结构的离散化,单元分析,单元集成,引入约束条件,求解线性方程组,得出节点位移。 4. 有限元法有哪些优缺点 答:优点:有限元法可以模拟各种几何形状复杂的结构,得出其近似解;通过计算机程序,可以广泛地应用于各种场合;可以从其他软件中导入建好的模型;数学处理比较方便,对复杂形状的

结构分析及有限元分析基础知识

第一章结构分析及有限元分析基础知识 注:摘自《NX知识工程应用技术——CAD/CAE篇》 洪如瑾编译 清华大学出版社 [目标] 本章将简述结构分析及有限元分析的基础知识,为学习与应用结构分析做好准备,包括: ※ 结构与结构分析定义 ※ 结构的线性静态分析 ※ 材料行为与故障 ※ 有限元分析的基本概念 ※ 有限元模型 1.1结构分析基础知识 1.1.1结构基本概念 1.结构定义 结构可以定义为一个正承受作用的载荷处于平衡中的系统。平衡条件意味着结构是不移动的。一个自由的支架不是一个结构,它未被连接到任一物体上并无载荷作用与它。仅当它附着到外部世界,并且有作用力、压力或力矩时,支架成为一个结构。 例如横跨江面的大桥就是一个普通的结构,一个支架通过它的支撑连接到地面上,桥的重量是在结构上的一种载荷(力)。当汽车通过桥时,附加的力作用于桥的不同位置。 一个好的结构必须满足以下标准: (1) 当预期的载荷作用时,结构必须不出现故障。这个似乎是显而易见的,并意味着结构必须是“强度足够的”。故障意味着结构破裂、分离、弯曲,以及支撑作用载荷失败。 注意:考虑到意外的载荷,通常在设计中提供安全余量。余量常常利用安全因素来描述。例如,如果在结构上期待载荷是10 000磅,规定安全因素是2.0,则结构将设计成能经受住20 000磅载荷。 (2) 当载荷作用时,结构必须不产生过分变形。这意味着结构必须“刚度足够”。 变形可接受的极限(弯曲度、挠度、拉伸等)取决于特定情况。例如,在通常住宅中的地板由足够的吊带支撑,以防止当人在地板岸上行走时有“柔软”的感觉。 (3) 在它的服务生命周期,结构的行为应不会恶化。这意味着结构必须“足够耐用”,必须考虑环境影响和“磨损与破裂”。如果一座桥假定维持50年,则桥的设计必须提供整个50年寿命的结构完整性与充分的安全余量。2.结构分析 结构分析是用于决定一个结构是否将正确完成任务的工程分析过程。结构将在某些方式中进行模拟和求解描述它的行为的数学方程。分析可以人工方法或用计算机方法来完成。 结构分析的结果(答案)用于评估性能,摘要如下: (1)“强度足够吗?”:应力必须是在一可接受的范围内。 (2)“刚度足够吗?”:位移必须是在一可接受的范围内。 (3)“耐用度足够?”:对一个长的疲劳周期应力必须足够低。

有限元法理论及应用参考答案分析

有限元法理论及应用大作业 1、试简要阐述有限元理论分析的基本步骤主要有哪些? 答:有限元分析的主要步骤主要有: (1)结构的离散化,即单元的划分; (2)单元分析,包括选择位移模式、根据几何方程建立应变与位移的关系、根据虚功原理建立节点力与节点位移的关系,最后得到单元刚度方程; (3)等效节点载荷计算; (4)整体分析,建立整体刚度方程; (5)引入约束,求解整体平衡方程。 2、有限元网格划分的基本原则是什么?指出图示网格划分中不合理的地方。 题2图 答:一般选用三角形或四边形单元,在满足一定精度情况,尽可能少一些单元。 有限元划分网格的基本原则: 1.拓扑正确性原则。即单元间是靠单元顶点、或单元边、或单元面连接 2.几何保持原则。即网络划分后,单元的集合为原结构近似 3.特性一致原则。即材料相同,厚度相同 4.单元形状优良原则。单元边、角相差尽可能小 5.密度可控原则。即在保证一定精度的前提下,网格尽可能的稀疏一些。(a)(b)中节点没有有效的连接,且(b)中单元边差相差很大。 (c)中没有考虑对称性,单元边差很大。 3、分别指出图示平面结构划分为什么单元?有多少个节点?多少个自由度?

题3图 答:(a )划分为杆单元, 8个节点,12个自由度。 (b )划分为平面梁单元,8个节点,15个自由度。 (c )平面四节点四边形单元,8个节点,13个自由度。 (d )平面三角形单元,29个节点,38个自由度。 4、什么是等参数单元?。 答:如果坐标变换和位移插值采用相同的节点,并且单元的形状变换函数与位移插值的形函数一样,则称这种变换为等参变换,这样的单元称为等参单元。 5、在平面三节点三角形单元中,能否选取如下的位移模式,为什么? (1). ?????++=++=2 65432 21),(),(y x y x v y x y x u αααααα (2). ?????++=++=2 65242 3221),(),(y xy x y x v y xy x y x u αααααα 答:(1)不能,因为位移函数要满足几何各向同性,即单元的位移分布不应与人为选取的 坐标方位有关,即位移函数中的坐标x,y 应该是能够互换的。所以位移多项式应按巴斯卡三角形来选择。 (2)不能,位移函数应该包括常数项和一次项。

有限元分析复习资料打印版

有限元复习资料 1.简述有限单元法的应用范围 答:①工程地质现象机制的研究;②工程区岩体应力边界条件或区域构造力的反馈;③工程岩土体位移场和应力场的模拟;④岩土体稳定性模拟 2.简述有限元单元法的基本原理 答:有限元单元法是随着电子计算机的发展而迅速发展起来的一种现代计算方法。它是50年代首先在连续体力学领域----飞机结构静,动态特性分析中应用的一种由此奥的数分析方法,随后很快广泛的应用于求解热传导。电磁场、流体力学等连续性问题。有限元分析计算的思路和做法可归纳如下: ①物体离散化 将整个工程结构离散为由各个单元组成的计算模型,这一步称作单元剖分。离散散后单元与单元之间利用单元的节点相互连接起来;单元节点的设置、性质、树木等应是问题的性质,描述变形形态的需要和计算进度而定(一般情况但愿划分月息则描述变形情况月精确,及月接近实际变形,但计算两越大)。所以有限元中分析的结构已不是原有的物体或结构物,而是同新材料的由众多单元以一定方式连接成的离散物体。这样,用有限元分析计算所获得的结果只是近似的。如果划分单元数目非常多而又合理,则所获得的结果就与实际情况相符合。 ②单元特性分析 A.选择位移模式 在有限单元法中,选择节点位移为基本未知量称为位移法;选择节点力作为基本未知量时称为力法;取一部分节点力和一部分节点位移作为基本未知量时称为混合法。位移法易于实现计算机自动化,所以,在有限单元法中位移法应用范围最广。当采用位移法时,物体或结构离散化之后,就可把单元总的一些物理量如位移,应变和应力等由节点位移来表示。这时可以对单元中位移的分布采用一些能逼近原原函数的近似函数予以描述。通常,有限元法我们就将位移作为坐标变量的简单函数。这种函数称为位移模式或位移函数,如y=a其中a 是待定系数,y是与坐标有关的某种函数。 B.分析但愿的力学性质 根据单元的材料性质、形状、尺寸、节点数目、位置及其含义等,找出单元节点力和节点位移的关系式,折中单元分析中的关键一部。此时需要应用弹性力学中的几何方程和物理方程来来建立力和位移的方程式,从而导出单元刚度矩阵,这是有限元法的基本步骤之一。C.计算等效节点力 物体离散化后,假定力是通过节点从一个单元传递到另一个单元。但是,对于实际的连续题,力是从单元的公共边传递到另一个单元中去的。因而,这种作用在单元辩解上的表面力、体积力和集中力都需要等效的移动节点上去,也就是用等效的节点力来代替所有作用在单元上的力。 ③单元组集 利用结构力的平衡条件和边界条件把各个单元按原来的结构重新连接起来,形成整体的有限元方程 ④求解未知节点位移 解有限元方程式得出位移。这里,可以根据方程的具体特点来选择合适的计算方法。通过上述分析,可以看出,有限单元法的基本思想是“一分一合”,分是为了进行单元分

有限元单元法复习资料

1.1有限单元法中“离散”的含义是什么?有限单元法是如何将具 有无限自由度的连续介质的问题转变为有限自由度问题的?位移有限单元法的标准化程式是怎样的?(1)离散的含义即将结构离散化,即用假想的线或面将连续体分割成数目有限的单元,并在其上设定有限个节点;用这些单元组成的单元集合体代替原来的连续体,而场函数的节点值将成为问题的基本未知量。(2)给每个单元选择合适的位移函数或称位移模式来近似地表示单元内位移分布规律,即通过插值以单元节点位移表示单元内任意点的位移。因节点位移个数是有限的,故无限自由度问题被转变成了有限自由度问题。(3)有限元法的标准化程式:结构或区域离散,单元分析,整体分析,数值求解。 1.2单元刚度矩阵和整体刚度矩阵各有哪些性质?各自的物理意义是什么?两者有何区别? 单元刚度矩阵的性质:对称性、奇异性(单元刚度矩阵的行列式为零)。整体刚度矩阵的性质:对称性、奇异性、稀疏性。 单元刚度矩阵Kij物理意义Kij即单元节点位移向量中第j个自由度发生单位位移而其他位移分量为零时,在第i个自由度方向引起的节点力。 整体刚度矩阵K中每一列元素的物理意义是:要迫使结构的某节点位移自由度发生单位位移,而其他节点位移都保持为零的变形状态,在所有个节点上需要施加的节点荷载。 2.1 为了使计算结果能够收敛于精确解,位移函数需要满足什么条件?为什么? 满足完备性和协调性。 原因:完备性包括两个条件:即刚体位移条件与常应变条件。首先,位移函数必须包含单元的刚体位移。结构中的单元不仅产生与该单元本身变形相应的位移,还可能因其他单元变形而通过节点位移产生单元刚体位移。为了正确反映单元的实际位移形态,位移函数必须具有反映刚体位移的能力。 其次,由于单元位移函数采用多项式,故在单元内部协调条件总能满足,要求反映在相邻单元之间。实质上来说,要求相邻单元间协调是为了保证单元交界面上应变有限。 3.1构造单元形函数有那些基本原则?试采用构造单元几何方法,构造T10单元的形函数,并对其收敛性进行讨论。 答:形函数是定义于单元内坐标的连续函数。通常单元位移函数采用多项式,其中的待定常数由节点位移参数确定,因此其个数应与单元节点自由度数相等。根据实体结构的几何方程,单元的应变是位移的一次导数。为了反映单元刚体位移和常应变即满足完备性要求,位移函数中必须包含常数项和一次项,即完全一次多项式。 3.3 何谓面积坐标?其特点是什么? 答:三角形单元中,任一点P(x,y)与其3个角点相连形成3个子三角形,其位置可以用下述称为面积坐标的三个比值来确定: L1=A1/A L2=A2/A L3=A3/A 其中A1,A2,A3分别为P23,P31,P12的面积。 各三角形面积为:Ai=1/2* =(ai+bi+ci)/2 由于A1+A2+A3=A,所以有L1+L2+L3=1,Li=(ai+bi+ci)/(2A) 特点:①T3单元的形函数Ni就是面积坐标Li ②面积坐标与三角形在整体坐标系中的位置无关,故称为局部坐标。 ③三个节点的面积坐标分别为节点1(1,0,0),节点2(0,1,0),节点3(0,0,1),形心的面积坐标(1/3,1/3 ,1/3)。④单元边界方程为Li=0 (i=1,2,3); ⑤在平行于2,3边的一条直线上,所有点都要相同的面积坐标。⑥面积坐标与直角坐标互为线性关系。 体积坐标:P点与四面体四个面围成的四个子四面体的体积与原来四面体体积的比值。即 剪切闭锁现象:当梁的高度与梁的长度之比t/l趋于零时,这种单元将出现这种现象,算得的挠度趋于零。 为克服剪切闭锁,使C0型单元适用于各种高度的梁。采用减缩积分方案与假设减应变法。 零能模式:对应于某种非刚体位移模式,减缩积分时高斯点上的应变正好等于零,此时的应变能当然也为零,这种非刚体位移模式称为零能模式。采用减缩积分时会发生零能模式。 5.1、等参单元:将整体坐标系中xy中形状中较复杂的真实单元变换成局部坐标系xy中规则的标准单元,然后在标准单元中构造形函数。由于坐标变换式与单元位移函数中用了相同的形函数N i(ξ,η),故称这种变换为等参变换,相应的单元称为等参单元。 2、等参单元的优越性:①有些工程较复杂,用直边单元离散这些结构需要大量的单元才能得到较好的近似,而曲边的等参单元可方便地离散复杂结构。②如在单元内多取些节点,单元便具有较多的位移自由度,从而就能够插值表示较复杂的单元内部位移场,这样就提高了单元本身的精度。③等参单元刚度矩阵、荷载矩阵的计算是在规则单元域内进行的因此不管被积函数多么复杂,都可以方便地采用标准化数值积分。 3、数值积分的阶次:对于N点积分,当被积函数为m次多项式且m<=2N-1时,可得精确积分值。反之,对于m次多项式的被积函数,精确积分要求的积分点数N>=(m+1)/2。 6.1、工程梁和剪切梁的基本假设?有两种梁弯曲理论①工程梁理论基本假设:平截面假设与横向纤维无挤压假设。前者认为梁横截面变形后仍为平面,且垂直于变形后的中性轴。该假设意味着横向剪切应变γxy =0,后者认为梁的横向纤维无挤压,即εy=0。②剪切梁理论基本假设:横向纤维元无挤压与另一假设认为法平面变形后仍为平面,但不再垂直于变形后的中性轴。 6.2. 剪切梁怎么考虑剪切影响:在结构单元分析中,可在工程梁单元的基础上考虑剪切变形的影响,也可通过挠度与转角各自独立插值直接构造剪切梁单元。 6.3对于杆系结构单元,为什么要在局部坐标系内建立单元刚度矩阵?为什么还要坐标变换?(1)在局部坐标系内可以更方便的建立单 元刚度矩阵。(2)在整体分析中,对所有单元都应采用同一个坐标系即整体坐标系X Y,否则围绕同一节点的不同单元对节点施加的节点力不能直接相加。因此,在进行整体分析之前,还需要进行坐标转换工作,把局部坐标系中得出的单元刚度方程转换成整体坐标系中的单元刚度方程,从而得出整体坐标系中的单元刚度矩阵。 7.1. 薄板弯曲理论基本假定:第一条:板厚方向的挤压变形可忽略不计,即εz=0,。这项假设类似于梁的横向纤维间无挤压假设。第二条:在板弯曲变形中,中面法线保持为直线且仍为弹性曲面(挠度曲面)的法线。第三条:薄板中面只发生弯曲变形,没有面内的伸缩变形,即中面水平位移(u)z=0=(v) z=0=0. 7.2. 厚板理论基本假设:板的中面法线变形后基本保持为直线,但因横向变形的缘故,该直线不在垂直于变形后的中面。因此,法线绕坐标轴的转角θx、θy不再是挠度的导数,而是独立变量;中面内的线位移和板厚方向的挤压变形也可忽略。 7.3. 薄板、厚板基本假定的不同:薄板:板弯曲变形中,中面法线保持为直线且仍为弹性曲面法线。厚板:板中面法线变形后仍基本保持为直线,但该直线不再垂直于变形后的中面。 7.4. DKT单元:离散Kirchhoff理论的基本思想是在若干离散点上满足Kirchhoff直法线假设。基于这种理论构造薄板单元时,w,θx,,θy 也各自独立插值;然后在若干离散点上引入直法线假设。这样构造的单元叫做DKT单元 8.1. 薄壳单元基本假设:薄壳理论假设:薄壳发生微笑变形时,忽略沿壳体厚度方向的挤压变形;且认为直法线假设成立,即变形后中面法线保持为直线且仍为中面的法线;壳体变形时中面不但发生弯曲,而且面内也将产生面内伸缩变形;折板假设;非耦合假设。 薄壳与薄板理论的假设的异同点:相同点:直法线假设和法向(板厚度方向)的纤维无挤压假设均成立。不同点:薄板中面只发生弯曲变形,没有面内的伸缩变形,即中面水平位移为零,而壳体变形时中面不但发生弯曲,而且也将产生面内伸缩变形。 厚壳分析的假设:变形前后的中曲面法线变形后仍基本保持为直线,但因横向剪切变形的缘故,该直线不再垂直于变形后的中曲面,此外,壳体厚度方向的挤压变形可以忽略。 与厚板理论的假设的 相同点:中面法线变形后仍基本保持为直线,但因横向剪切变形的缘故,该直线不在垂直于变形后的中面。厚度方向的挤压变形忽略不计。不同点:厚板理论的假设中,中面内的线位移可以忽略,而厚壳理论的假设中,中面内的位移不可忽略,并且厚壳的位移场可用中面位移表示。 8.2. 平板型单元:组成的折板系统去代替原来的壳体,由平面应力状态与平板弯曲应力状态加以组合而得壳体的应力单元。 分析这种单元时所提出的假设:理论假设:薄壳发生微笑变形时,忽略沿壳体厚度方向的挤压变形,且认为直法线假设成立,即变形后中面法线保持为直线且仍为中面的法线。,折板假设,非耦合假设。 应用平板型壳单元可能会出现的问题,如何解决:1.单元共面问题,解法:引入唯一边界条件可解方程Ka=P 。2.虚拟旋转刚度,解法:在特殊节点上给以任意的虚拟刚度系数。Kθzθzθzi=0,经坐标变换,整体坐标系中该节点平衡方程将满足唯一解条件。赋予Kθzθz任何值。3.新型平面膜元:在平面膜元角点上增加旋转自由度θz,使其有对应的刚度。 8.3. 面内变形与弯曲变形之间非耦合的假设是针对什么提出来的?试说明单元组装时,面内效应与弯曲效应的耦合将会出现。 答:面内变形与弯曲变形之间非耦合的假设是针对局部坐标系下的单元提出的。 9.1. 减少自由度的措施有哪些?各自基本概念如何? 答:1.恰当利用结构对称性。基本思想:利用结构的对称性,取结构一部分建立有限元模型。根据荷载对称性,分析对称面上的位移状态,以确定对称面上节点的位移边界条件。2.采用子结构技术。基本思想:在大型复杂结构的有限元分析中,可将原结构分成若干区域,每个区域作为一个子结构,这些子结构在其公共边界上互相连接起来。 2. 为什么说位移法中应力解的精度低于位移解? 答:在位移有限单元法中,沿单元边界是连续的,而位移的导数通常不连续,因此,在单元边界上应力是不连续的;基本未知量是位移,而单元应变和应力是由位移求导得到的,因此应力精度低于位移精度。 3. 在无法获得精确解的条件下,如何进行误差估计? 答:有限元解法的误差估计有:残值法,后处理法。后处理法:由于无法获得精确解,一般以修匀后的改进值σ*作为“精确解”进行误差估计,通过与精确值误差范数对比,这样做非常有效。

相关文档
最新文档