有限元平面钢架.docx

第一章问题重述一、题目内容：图示平面钢架结构图1.1 题目内容二、题目要求：（1）采用平面梁单元进行有限元法手工求解，要求写出完整的求解步骤，包括：a）离散化：单元编号、节点编号；b）单元分析：单元刚度矩阵，单元节点等效载荷向量；c）单元组装：总体刚度矩阵，总体位移向量，总体节点等效载荷；d）边界条件的引入及总体刚度方程的求解；e）B点的位移，A、C处支撑反力，并绘制该结构的弯矩图、剪力图和轴力图。（2）编制通用平面钢架分析有限元Matlab程序，并计算盖提，与手工结果进行比较；（3）利用Ansys求解，表格列出B点的位移，A、C处支反力，绘制弯矩图、剪力图和轴力图，并与手算和Matlab程序计算结果比较。（4）攥写报告，利用A4纸打印；（5）心得体会，并简要说明各成员主要负责完成的工作。第二章原理说明一、平面刚架的解题思路（步骤）：（1）?结构离散化：将每个杆件作为一个单元，杆件与杆件之间的联结点为节点，并进行编码。（2）单元刚度矩阵的建立：具体内容有：建立局部坐标系、单元刚度矩阵的导出、坐标转换、单元刚度矩阵的坐标变换、单元刚度矩阵的特性。（3）求总体刚度矩阵：总体刚度矩阵的分块形式、总体刚度矩阵的集成、总体刚度矩阵的特性。（4）位移的求解：支撑条件的引入、非节点载荷的处理。（5）求应变、应力二、单元刚度矩阵的建立（1）局部坐标系与整体坐标系：由于各杆件的取向不同，所以对各杆单元必须先建立各自的局部坐标系。以节点1(e)为原点，由1(e)到2(e)的方向x轴正向的右手坐标系就叫做该单元的局部坐标系，用英文小写字母x和y表示。如图2-1所示。图2-1 一维单元（钢架杆单元）为便于进行整体分析，要求整个结构的所有杆件都采用同一个坐标系，如图所示。该坐标系就称为整体坐标系，用英文大写字母X和Y表示：在局部坐标系中，杆单元每个节点都有三个位移分量，如图(a)所示。每个节点都有三个力分量，如图（b）所示图中——轴向位移、轴向力——横向位移、横向力——角位移、弯矩图中所示各力的方向为正。每个单元在两端共有六个节点位移和六个节点力分量（2）单元刚度矩阵的导出：对于刚架杆单元，由于位移和力向量都不再是二维而是六维，所以单元刚度矩阵也不再是2×2阶而是6×6阶，其基本形式为：求单元刚度矩阵的核心是求各单元刚度系数。下面求单元刚度矩阵。的第一列对应于只有节点1(e)在轴向方向产生单位位移时的情况，如下图所示。这时除了外，两节点的其余位移分量都为零。要产生这个变形并能保持该杆单元的平衡，只需要在两个节点处施加两个轴向力他们分别为刚度系数。很明显，第一列的其它元素均为零。的第二列对应于只有节点1(e)在横向方向产生单位位移时的情况，如下图所示。这时除了外，两节点的其余位移分量都为零。它相当于自由端绕度为1、转角为0的一根悬臂梁。要产生这个变形并能保持该梁的平衡，在节点1(e)处要施加横向力在节点2(e)处要施加横向力和弯矩，他们分别为刚度系数、。同样。的第三列对应于只有节点1(e)在产生单位角移时的情况，如下图所示。这时除了外，两节点的其余位移分量都为零。它相当于一根悬臂梁，自由端有铰链支撑，在节点1(e)处产生单位转角。要产生这个变形并能保持该梁的平衡，在节点1(e)处要施加横向力：在节点2(e)处要施加横向力和弯矩：他们分别为刚度系数、。同样。中的后三列元素也可以用同样的方法导出。局部坐标系中杆单元(e)的节点力和节点位移之间的关系式为：单元刚度矩阵是一个6×6阶的矩阵（3）坐标变换：前面推导时采用的是局部坐标系，在研究刚架或桁架结构时，由于各杆的取向不同，存在若干个局部坐标系，这就给整体分析带来了麻烦。因此必须通过坐标变幻，将在局部坐标系中获得的单元刚度矩阵转换成整体坐标系中的单元刚度矩阵。假设一个杆单元两个节点的局部码1(e)和2(e)所对应的总码分别为I和j，则可写出一个坐标转换的矩阵方程（过程略）：式中——单元（e）的坐标转换矩阵，它是一个正交矩阵，即其逆阵等于其转置阵所以上式又可写为——局部坐标系x、y和整体坐标系X、Y两者之间相差的角度。规定由X轴到x轴以逆时针方向为正。——单元（e）在局部坐标系中按节点局部码标记的节点位移向量——单元（e）在整体坐标系中按节点局部码标记的节点位移向量对于单元节点力向量也存在上述关系，即式中分别为与相对应的节点力向量（4）单元刚度矩阵的坐标变换整体坐标系中单元节点力和节点位移之间的关系式式中——整体坐标系中单元(e)的刚度矩阵。因为所以。三、总体刚度矩阵前面首先推导出了局部坐标系中的单元刚度矩阵，然后通过坐标转换获得了整体坐标系中的单元刚度矩阵。在得到所有单元在整体坐标系中的单元刚度矩阵的基础上，可以建立总体刚度矩阵。总体刚度矩阵的集成由整体坐标系中的单元刚度矩阵的子矩阵集成总体刚度矩阵的步骤如下