计算固体力学第三章_1.ppt

进一步考察 上式可写为: 说明 1: 上式中{F}是将体积力和分布载荷堆积在适当节点上得出的等价总体节点载荷,包括集中节点载荷. 有关体积力和面积力的处理，在下节给出。 说明2: 上面的推导是对于总体坐标中的一般方向,因此其中的所有单元矩阵可以适用于所有单元. 不需要从局部坐标向总体方程转换. 如局部坐标轴和整体坐标轴不平行时,必须利用转轴变换方程来建立单元刚度矩阵,单元节点力和位移矩阵。参看弹性力学. 节点位移解出后,可以得到应变和应力. 为了保证解答的收敛性，要求位移模式必须满足以下三个条件： 同时满足3个条件的，称为完备的协调单元。 满足1和2，但是不满足3的，称为完备非协调单元。 完备非协调单元的缺点是不能保证计算结果的收敛性，但是，有时候这类单元的收敛速度反而比协调元快。 已经证明，对于完备的协调单元，起刚度系数的数值比精确的要大。这样，在给定载荷之下，计算模型的辨析功能比实际结构的要小。因此，当单元网格分割的越来越细时，位移的近似解由下方收敛于精确解。 选择多项式位移模式的阶次时，要考虑完备性和协调性的要求。还要考虑位移模式应该与局部坐标系的方位无关，这种性质称为几何各向同性。 例题: 计算图中的单元的刚度矩阵? 单位为英寸. 假定平面应力状态.假定单元节点位移已确定,确定单元应力?具体参数如下: (*) 利用方程 参看本节步骤 2, 3,可以得到 代入方程 (*) 同理,利用下面方程计算应力 * * （５）主轴坐标系内的梁单元的刚度矩阵 主轴坐标系内的梁单元的刚度矩阵K可以通过组装前面　（１），（２），（３），（４）步中的刚度，按照排列组成 12*12 的完整梁单元刚度矩阵． （6）节点坐标系内的力-位移关系式 （7）基准坐标系内的力-位移关系式 一般情况下,梁的形心和节点在空间不重合,他们之间的差距称为”偏心”。因此，要根据虚功原理，把主轴坐标系内导出的梁元关于形心的力-位移模式转换成关于节点的力-位移模式。（相当于移轴） 为了得到整个结构的刚度特性，必须要对单元进行组合，这使得所有的节点位移和他们对应的节点力系将参考一个公共的坐标系，这个坐标系就是基准坐标系。要建立各个单元节点坐标系内的力-位移关系式与基准坐标系内的力-位移关系式之间的转换关系。（相当于移轴加转轴） 例：如图，有中间支点的悬臂梁，端点受载荷 P 作用，假定梁的 EI 是常数，长度为2L ，假定中点有滚轴支撑，右端为固支． 确定节点位移，转动和反作用力，并作剪力和弯矩图． 按图建立离散单元，２个单元，并建立总体坐标轴，得到总体刚度矩阵 及 将已知的位移和转角回代入控制方程， 可以确定总体坐标节点力。 作业1： 参照上面方法，对于图P4.3所示的梁，确定铰支点A的转动以及载荷P作用处的转动和位移．确定反作用力，画出剪力和弯矩图，令EI沿梁的长度是常数． ３. 矩阵位移法 自学 作业2： ４. 平面应力和平面应变概念 （1） 平面应力 平面应力定义为一种应力状态，在这种应力状态下，假定垂直于该平面的法向应力和剪应力为零。即应力主要发生在平面内。例如很多薄板，如： （2） 平面应变 平面应变定义为一种应变状态，在这种应变状态下，假定垂直于该平面的线应变和剪切应变为零。 对于有固定横截面的长物体（如 z 方向），如果仅仅作用有(x,y方向)的载荷，并且沿长度方向不变化，就可以假定为平面应变状态。如： 5. 平面三角形单元 为了能用节点位移表示应变和应力，首先必须假设一个位移模式，也就是假定位移分量为坐标的某种简单函数，这些函数在节点上的数值，应当等于节点上的位移分量的数值。这些函数习惯上称为“形状函数”，或“插值函数”。 在协调模型中，多用多项式作为插值函数。 考察平面三角形单元。 考虑一个薄板受拉伸作用，如图 对薄板采用限元法分析: 考虑从右图中取出基本的三角形单元，每个单元有节点 i, j, 和m。每个节点有两个自由度。这里采用逆时针方向记录节点号。 选择每一个单元的位移函数为线形函数 式中u(x,y)和v(x,y)描述单元内部任何一点 (xi, yi)的位移。 线形函数保证满足协调关系． 有指定端点的线形函数只有一条路径（直线）通过两个端点．因此，线形函数可以保证沿边界和在相邻单元共享的节点处的位移是相等的．利用方程（１），可以得到包容函数　u 和 v 的广义位移函数为： （１） 回84 将节点坐标代入方程（２），可以得到 （２） （３） 将方程（3）的前，后三个方程分别写成矩阵形式： 代回方程 （１）得到 （4） A 为三角形面积 其中 （5） 可以将式（４） 写成 记为 将方程　　　　　　　　代入上式， 得到以节点位移表示的应变：　　　　　　 其中 包含平面应变和平面应力 其中 其中{f} 是广义位移函数, {X}是单位体积重量或重量密度矩阵