计算力学第4课(2015).pptVIP

* 在等参元的坐标变换中，局部坐标(自然坐标)中的正方形或立方体通常称为母单元，而整体坐标(直角坐标)内曲边图形则称为子单元。 * 用单元的直边边界去拟合物体的曲线边界时，总不可避免地降低了计算模型逼近实际结构的精度。因此构造这些具有曲边的高精度单元就成为必要。它可使结构在离散时采用较少的单元而达到较高计算精度，以边长为2的8节点正方形单元为例来说明这种映射过程。 * 用 * 采用试凑法，设，它可以保证N1在253边的节点2、5、3处的值为零，在46线上的节点4、6处的值为零。根据N1在1节点处为1的条件，可得 * * * * * 式中，，是容易求得的，代入上式即可求得[J]，求逆，然后可求得，，即可求得应变矩阵[B]。 * * * * 对协调元(包括应用成功的非协调元)，节点位移无须区分，也不能区分是属于哪个单元的，必须注意。 * 对协调元(包括应用成功的非协调元)，节点位移无须区分，也不能区分是属于哪个单元的，必须注意。 * 对于任一个节点，可能受有两种力：单元给予节点的反作用力（即节点力反号）和节点荷载，前者属于内力，而后者属于外力，每个节点在这两种力的作用下处于平衡。 * 考虑整体能量泛函极值，同样可以得到整体刚度方程。 * 上述两种考虑，实质上是同一个整体平衡问题采用两种不同的形式去反映的。前者，用各个节点的平衡去推演整体平衡；后者，用各个单元的能量去组集成整体能量，并根据能量极值原理导出整体平衡。 * * 的计算问题，由于不能列出x，y表示?，?的显式，故 也无法求得 的算式，因此无法求得 。 由复合函数求导法则，改用求逆的方法 于是形函数对x，y的导数应为 平面问题的有限元法 等参元应力矩阵 平面问题的有限元法 等参元的刚度矩阵 等参元的刚度矩阵仍用式 计算，只不过式中的微面积dxdy应换成 所包围的四边微元面积。 图中，在Oxy坐标系内的任一微矢量 ，因为x，y是?，?的函数，所以 故 平面问题的有限元法 故 平面问题的有限元法 整体刚度方程 节点平衡 单元贡献方程 整体节点位移为 整体节点荷载为 平面问题的有限元法 这样对每个单元写出相应的平衡方程，例如，对①号单元 称为单元贡献刚度方程 平面问题的有限元法 节点平衡方程 以节点3为例 同理由其它四个节点，则整体节点荷载与节点位移之间的关系式为 式中 [K]为整体刚度矩阵 平面问题的有限元法 能量极值 单元在节点力作用下的势能为 将各单元能量叠加，可得整体能量 平面问题的有限元法 令 则 整体刚度矩阵的特性 对称性 正定性 稀疏性 平面问题的有限元法 整体节点载荷列阵 集中力：把集中力的作用点取为节点，不需作任何特殊处理，就可直接把集中力加入到节点荷载中去。 体积力 ： 设单元体内单位体积上作用的体积力为 将移置到节点上后的单元节点荷载记为 根据静力等效原则，给单元一个虚位移 表面力 ：与体积力类似 平面问题的有限元法 整体荷载列阵 平面问题的有限元法 小结 公共信箱：heu.tulixue@163.com 密码：tlx000 * 两端铰接的任一直杆单元，长为l，横截面积为A，置于x轴上。杆两端分别受有轴向节点力Xei，Xej作用，相应的杆端节点位移为ui、uj。 将位移函数用节点位移乘上一个已知插值函数之积的和来表示。 * * 等截面梁单元，是分析平面刚架系统的基本单元，由于载荷处于同一平面内，故单元处于轴向受拉(或受压)和平面弯曲的组合变形状态。取x轴与梁的轴线重合，而y轴和z轴为梁的截面主惯性轴。设I为截面惯性矩，A为截面积，l为梁长度。 * * 一般说来，梁在横向载荷作用下所产生的变形的影响是微不足道的，但是当截面高度大于长度的1/5时，剪切变形的影响必须考虑。贝努利梁。(铁摩辛柯最著名的工作是在梁的振动问题中计入了旋转惯性与剪力,这种模型后来被称为“铁摩辛柯梁”。) * 三角形单元的优点是它的“适应性”，