有限元分析与应用_第7讲、有限元方法的一般步骤精选.ppt

有限元方法的一般步骤 结构离散化 选择位移函数 单元刚度局阵的建立 单元刚度局阵的投放形成全结构总体刚度矩阵 结构平衡方程求解 计算单元的各种结果(应力、应变、位移、温度等) 一、结构离散化 有限自由度系统代替具有无限自由度数目的系统,离散基本是靠工程经验,在选择单元的形状、尺寸、数量和排列时必须谨慎，以便尽可能精确地模拟原物体，而不增加求解的计算工作量。 结构离散应该考虑的问题： 1、单元的选择 通常,应根据物理问题本身来选择单元的类型.例如,如果问题属于分析在给定的一组载荷条件下的桁架结构,那么,为了理想化而选用的单元类型显然是如图所示的杆单元或线单元 对下图所示的短梁作应力分析时,可以用三维体单元进行有限元理想化 但是,在某些情况下，用作理想化的单元类型可能不明显。此时，必须谨慎地选择单元的类型.在分析如图所示的薄壁壳体问题时,可以用几种类型的单元把壳体理想化,此时,需要的自由度数目,预期的精度,容易推导所需的方程以及实际结构可模拟的准确程度将决定理想化所用单元类型的选择。 2、单元的尺寸 　　 由于单元的尺寸直接影响解的收敛性，因而必须小心地加以选择，单元尺寸越小，最后的解就越精确。但应当记住，采用小尺寸的单元也就意味着需要更长的计算时间.有时，对同一物体可能要使用不同尺寸的单元。 另一个影响有限元解的与单元尺寸有关的性质是单元的纵横比.纵横比描述了单元在单元集合中的形状,对一个二维单元来说,纵横比取为单元的最长尺寸与最短尺寸之比.纵横比几乎等于1的单元往往能得出最好的结果. 4、单元的数量 　　 为了理想化而选择的单元数量,与所要求的精度,单元尺寸以及所涉及的自由度数目有关.虽然增加单元数量通常意味着有更精确的结果,但对某一个给定的问题来说,存在着某个单元数,超过了这个数目,再也不会在有效数上增加精度,图中曲线示出了这种特点.此外,由于采用大量单元包含了大量的自由度,因此,在可用的计算机内存中可能没有能力来存放由此产生的矩阵. 5、利用物体的物理条件进行简化 　　 对称面上的边界条件可以按以下规则确定: ⑴在不同的对称面上,将位移分量区分为对称分量和反对称分量,例如在ox面上,u是对称分量, ν是反对称分量;而在oy面上, ν是对称分量,而u是反对称分量. ⑵将载荷也按不同的对称面,分别区分为对称分量和反对称分量. ⑶对于同一个对称面,如载荷是对称的,则位移的反对称分量为0;如载荷是反对称的,则位移的对称分量为0. 上述规则表示:一个结构若外载及结构相对于某轴(面)对称,则对称面上自由节点相对于对称轴(面)作为反对称移动的位移为零.反之,若结构对称,载荷反对称,则相对于对称轴(面)作为对称移动的位移分量为零.利用这个性质:可以确定对称轴(面)上结点的支持情况,只需要计算结构的一部分. 6、总刚度矩阵的性质与结点编号 用有限元法分析结构时，通常结点的自由度数目巨大,从而使方程组有大型的特点，[K]往往有数百阶甚至几千阶，而且[K]中的每一行势必含有大量的零元素，叫做稀疏矩阵。总刚度矩阵中非零元素集中在主对角线两侧，又使总刚度矩阵呈带状。考虑到刚度矩阵的对称性，因而定义包括主对角元素在内的一侧非零元素最大延伸长度，叫刚度矩阵的半带宽。 二、选择位移函数 * 孔洞附近(此处存在应力集中)的单元尺寸应远小于远离孔洞的单元尺寸.通常,每当预料会有陡的应力梯度时,在这些区域总是要采用更密的网格. 　　 3、结点的设置 　 如果物体在几何形状、材料性质和外部条件(如载荷、温度)方面没有突然变化,则可把物体分为相等的小部分,从而可使结点的间距均匀.另一方面,如果问题存在有任何间断,则显然应当在这些间断处设置结点,见图 工程实际中,很多结构具有对称性.如能恰当地加以利用,可以使结构的有限元计算模型以及相应的计算规模得到减缩,从而使数据准备工作和计算工作量大幅度地降低.如果物体的形状以及外加条件都是对称的,则在进行有限元理想化时就可以只考虑物体的一半.但是在求解过程中,必须加入对称条件,如图所示一具有中心圆孔的矩形板,由于板不可能在对称线A-A有水平位移,所以在求解时必须在对称线上引入u=0的条件. 如图(a)所示的桁架,全结构四个结点,五个杆单元.但杆元1—3的轴线是对称线,而且载荷p也是对称的.对称线上的结点为1、3,利用上述性质,在对称线上(即结点1、3) 为反对称位移,它必等于零.这也相当于在 方向上,在1,3结点处加了刚性支座,于是计算简化为图(b)所示结构 一般说来，计算最大半带宽得公式为 其中，B是最大半带宽，D是从全部