第三章杆梁结构的有限元原理讲述.ppt

* * （c）剪力图 （b）弯矩图 （d）轴力图 （a）位移曲线 * 例 4 * (1)单元划分，建立局部坐标系和整体坐标系，其中局部坐标下的x方向如箭头所示，y方向向下，z轴依据右手法则 （2) 形成局部坐标下的单刚 * 将以上数据代入单刚矩阵 且有 * * * （3） 单元位移分量：依据先处理法：即先引入支承条件，根据支承条件仅对未知的自由结点位移分量进行编号，得到位移分量中不包含已知约束的位移分量 （4） 整体刚度矩阵：依据位移分量的对应情况，可以得到整体刚度矩阵，叠加同位置的单刚（对应零位移分量去掉） * 依据位移分量的对应情况，可以得到坐标下的结点力（对应零位移分量去掉） * * * * * 简 写成： 　　　　 这里 为结构的整体刚度矩阵，有： 　　　　 注意，在建立方程的过程中，我们假设所有结点都有位移。因此整个结构在外力作用下，除了发生弹性变形外，还可能发生刚体平动位移，这样各结点位移不能唯一确定。这说明整体刚度矩阵为一奇异矩阵，不能求逆矩阵，即根据整体刚度方程可得到无穷多个解。 * 实际上，在图所示刚架中，结点1和结点4均为固定端，其三个位移分量均为0，即有： 　　　　 这样，将上述支承条件引入到方程中，对整体刚度方程进行修改，可得： 　　　　 * 对上述方程进行化简，可以得到两组方程： 　　　　 和 　　　　 这样，利用第1式可以求得结点位移 和 ，再根据第2式可以求得支座反力 和 。 * 2　先（前）处理法 　　所谓先处理法，就是先引入支承条件，根据支承条件仅对未知的自由结点位移分量编号，得到的位移矩阵中不包含已知的约束位移分量，即可以直接得到方程求解自由结点位移分量。 * 　　如图所示的平面刚架结构ABCD，由于在A处和D处均为固定端，其位移为0，故位移编码均为0，在C处为铰接，故BC杆在C端的角位移与DC杆在C端的角位移不相同，因此在C处编两个结点3和4，但结点3和4的横向位移和竖向位移相同，故采用相同的编号，各结点位移编码如图所示。 图　先处理法位移编码示意图 * 3.4　等效结点荷载与边界条件的处理 非结点等效荷载和边界条件的处理是有限元分析中必须考虑的的两个重要方面。由于只考虑结点荷载，因此必须将作用于单元上的非结点荷载转换到节点上。 有限元的结点荷载来自两部分 （1）作用于结点处的集中力、集中力偶，前面多次提到的 直接叠加到结点上即可（整体坐标下） （2）非结点荷载首先须在局部坐标下等效到结点荷载 ， 然后再转到整体坐标系下的 * 1　非结点荷载的处理 　在前面的分析中，我们已经介绍了求等效结点荷载的方法，如3-7式、3-15式、3-22式分别可用来求不同情况下的等效结点荷载。 　　 此外，可以这样来考虑： 　　第一步，在局部坐标系下求单元的固端力 。对于某个单元，我们假定单元的两端均固定，然后根据静力平衡求得固定端的反力。 　　第二步，根据单元固端力求单元的等效结点荷载 。根据局部坐标系与整体坐标系单元杆端力的变换式，固端内力在两种坐标系下的变换形式可以写成： 　　　　 因此，整体坐标系下的等效结点荷载矩阵可以由下式计算： 　　　　 （3-7） （3-15） （3-22） （拉压杆） 纯弯杆 扭转杆 * 平面刚架单元固端力（遵行结构力学） * * 2　边界条件的处理 　 （1）　铰结点 　　在杆系结构中，除了刚性结点外，通常会遇到一些杆件通过铰结点与其它杆件联结，如下图所示杆件系统，有4根杆件汇交于D点，其中BD杆在D端通过铰支座与其它杆件铰接，其余3根杆为刚性接触。对于这样的铰结点，具有如下的性质： 　　①　铰结点上各杆具有相同的线位移，但截面的转角位移不相同； 　　②　结点上具有铰接杆端不承受弯矩作用。如下图所示结构中，BD杆在D端的杆端弯矩为0，只有CD、ED、GD杆在结点D上与外弯矩保持平衡。 * 　　对于这样的结点，我们在对其进行单元划分时，通常考虑在D处设置2个结点。按照先处理法，对图示结构进行位移编码，如图2b所示。 图　铰结点的处理示意图 ANSYS通过结点耦合实现 * 2　弹性支承点 　　在实际工程中，有时会遇到弹性支承的情况（如图），这时一般将弹性支座看作是在结构约束点沿约束方向的一个弹簧，弹簧的刚度系数为 ， 在数值上等于使弹簧支座沿约束方向产生单位位移时所需施加的力。 ANSYS引入弹簧单元即可 * 具体做法可以归结为：先解除弹性支承点约束，在i处给一个结点号，形成总刚度矩阵，然后在总刚度矩阵中将第i行的主元素加上弹性支承的刚度系数，此时第行变为