有限元f复习1解说.doc

第  PAGE \* MERGEFORMAT 21 页 共  NUMPAGES \* MERGEFORMAT 21 页 一、简述两种平面问题及其基本方程。（10分） 弹性体在满足一定条件时，其变形和应力的分布规律可以用在某一平面内的变形和应力的分布规律来代替，这类问题称为平面问题。平面问题分为平面应力问题和平面应变问题。 1）.平面应力问题 设有张很薄的等厚薄板，只在板边上受到平行于板面并且不沿厚度变化的面力，体力也平行于板面且不沿厚度变化。设板的厚度为t ,在板面上： ， ， 。 由于平板很薄，外力不沿厚度变化，因此在整块板上有：，， 剩下平行于XY面的三个应力分量 未知 2）.平面应变问题 设有很长的柱体，支承情况不沿长度变化，在柱面上受到平行于横截面而且不沿长度变化的面力，体力也如此分布。 以柱体的任一截面为XY平面，任一纵线为Z轴，假定该柱体无限长，则任一截面都可以看作对称面，由对称性： ， ， 未知量为平行于Z平面的三个应力分量，物体在Z方向处于自平衡状态。 平面问题的基本方程为： 平衡方程 几何方程 物理方程 （弹性力学平面问题的物理方程由广义虎克定律得到） 平面应力问题的物理方程 平面应力问题有 平面应变问题的物理方程 平面应变问题有 在平面应力问题的物理方程中，将E替换为、替换为，可以得到平面应变问题的物理方程；在平面应变问题的物理方程中，将E替换为、替换为，可以得到平面应力问题的物理方程。 二、简述有限元中非线性的问题。（10分） 通常把结构非线性问题分为三大类：几何非线性和材料非线性和状态非线性。这主要包括三个方面：一、是在大位移问题中，尽管位移很大，结构的应变仍然不大，属于大位移小应变问题，材料的应力-应变关系仍是线性的，只是应变-位移关系是非线性的。物体经历大的刚体位移和转动，固连于物体坐标系中的应变分量仍假设为小量。二、是非线性效应由应变应力关系的非线性所引起，位移分量仍假设为小量，应力-应变关系是非线性的，即材料非线性问题；最一般的情况是位移、转动和应变都不再是小量，不但位移-应变是非线性的，而且应力-应变关系也是非线性的，即双重非线性问题。对于结构的几何非线性和材料非线性分析，可以归结为外力与内力的平衡方程，它是关于节点位移的非线性方程；非线性的稳态与瞬态温度场计算归结为热流平衡方程，它是关于节点温度的非线性方程；因此非线性分析的有限元计算最终归结为非线性方程求解。非线性分析简而言之就是：将系统的平衡方程式根据系统的非线性特性不断地进行修正，然后求平衡方程的增量解如果是几何非线性，则在新的一步增量求解之前，坐标系进行修正，然后去求解方程，并计算几何非线性对刚度阵和载荷阵的修正。若为材料非线性，则是将等效刚度阵和载荷阵不断地进行修正，然后进行求解。 三、阐述动力学有限元方程。（10分） 其中，[M]—结构的总质量矩阵； [C]—为阻尼矩阵； [K]—结构的总刚度矩阵； [u]—结构的位移向量； {R（t）}—强迫力列阵。 四、简述弹性薄板理论的基本假设和基本方程。（10分） 有限单元法第114-118页 五、就1D杆单元（30分） 节点位移（局部坐标下） 节点位移（整体坐标下） 写出和之间的关系 将该单元的位移场、应力场、应变场用整体坐标系下的节点位移q来表示。 推导出基于整体坐标下的刚度矩阵。 答：（1）如图所示： ∴ 其中 （2）在局部坐标下，设位移场模式（有两个节点）为： （，待定） 由边界位移知 解之，知：， ∴ 其中 ∴ 其中 在整体坐标系下有 （3）系统的势能为： ＝ 其中，K为刚度矩阵 六、如图所示为一对称均质薄板，两边分别为边长1m的等边三角形，厚度t=1cm,在对称中心处受垂直压力q=N,材料弹性模量为E,泊松比，不计自重，试用有限元法求其应力分量(40分)。 y 1m q 图一 图二 由于此结构长、宽远大于厚度，而载荷作用于板平面内，且沿板厚均匀分布，故可按平面应力问题处理，考虑到结构和载荷的对称性，可取结构的1/2来研究 划分网格 节点 坐标 1 2 3 x y 0 0 0 1 该1/2结构被离散为一个三角形 单元，节点编号，单元划分及 取坐标如图2所示，其各节点 表 1 的坐标值见表1。 求单元的刚度矩阵 计算单元的节点坐标差 及单元面积 单元１（i、j、m　1，2，3） 先计算用到的常数  所以刚度矩阵为: 引入约束条件，修改刚