有限元分析和应用_第7讲、有限元方法一般步骤.pptVIP

有限元方法的一般步骤;一、结构离散化 有限自由度系统代替具有无限自由度数目的系统,离散基本是靠工程经验,在选择单元的形状、尺寸、数量和排列时必须谨慎，以便尽可能精确地模拟原物体，而不增加求解的计算工作量。 结构离散应该考虑的问题： 1、单元的选择 通常,应根据物理问题本身来选择单元的类型.例如,如果问题属于分析在给定的一组载荷条件下的桁架结构,那么,为了理想化而选用的单元类型显然是如图所示的杆单元或线单元;对下图所示的短梁作应力分析时,可以用三维体单元进行有限元理想化; 但是,在某些情况下，用作理想化的单元类型可能不明显。此时，必须谨慎地选择单元的类型.在分析如图所示的薄壁壳体问题时,可以用几种类型的单元把壳体理想化,此时,需要的自由度数目,预期的精度,容易推导所需的方程以及实际结构可模拟的准确程度将决定理想化所用单元类型的选择。; 2、单元的尺寸 　　 由于单元的尺寸直接影响解的收敛性，因而必须小心地加以选择，单元尺寸越小，最后的解就越精确。但应当记住，采用小尺寸的单元也就意味着需要更长的计算时间.有时，对同一物体可能要使用不同尺寸的单元。; 另一个影响有限元解的与单元尺寸有关的性质是单元的纵横比.纵横比描述了单元在单元集合中的形状,对一个二维单元来说,纵横比取为单元的最长尺寸与最短尺寸之比.纵横比几乎等于1的单元往往能得出最好的结果. ; 4、单元的数量 　　 为了理想化而选择的单元数量,与所要求的精度,单元尺寸以及所涉及的自由度数目有关.虽然增加单元数量通常意味着有更精确的结果,但对某一个给定的问题来说,存在着某个单元数,超过了这个数目,再也不会在有效数上增加精度,图中曲线示出了这种特点.此外,由于采用大量单元包含了大量的自由度,因此,在可用的计算机内存中可能没有能力来存放由此产生的矩阵.;5、利用物体的物理条件进行简化 　　;对称面上的边界条件可以按以下规则确定: ⑴在不同的对称面上,将位移分量区分为对称分量和反对称分量,例如在ox面上,u是对称分量, ν是反对称分量;而在oy面上, ν是对称分量,而u是反对称分量. ⑵将载荷也按不同的对称面,分别区分为对称分量和反对称分量. ⑶对于同一个对称面,如载荷是对称的,则位移的反对称分量为0;如载荷是反对称的,则位移的对称分量为0. 上述规则表示:一个结构若外载及结构相对于某轴(面)对称,则对称面上自由节点相对于对称轴(面)作为反对称移动的位移为零.反之,若结构对称,载荷反对称,则相对于对称轴(面)作为对称移动的位移分量为零.利用这个性质:可以确定对称轴(面)上结点的支持情况,只需要计算结构的一部分.;如图(a)所示的桁架,全结构四个结点,五个杆单元.但杆元1—3的轴线是对称线,而且载荷p也是对称的.对称线上的结点为1、3,利用上述性质,在对称线上(即结点1、3) 为反对称位移,它必等于零.这也相当于在 方向上,在1,3结点处加了刚性支座,于是计算简化为图(b)所示结构;6、总刚度矩阵的性质与结点编号;二、选择位移函数;位移函数的多项式形式;插值多项式阶次的选择 　　在选择插值函数多项式的阶次时,必须考虑到下列因素: 　　　　(1)多项式描述的位移形式与局部坐标系无关; 　　　　(2)ai的数目应等于单元结点自由度的数目. (3)插值多项式应当尽可能满足收敛性要求; 第(1)条要求即在不同局部坐标系中,位移函数(多项式)表达式不变,这种性质称为几何等向性,为获得几何等向性,多项式中应含有不违反对称性的那些项.如二维线性单元中,含有X和Y两项,而不仅仅是其中一项.或在二维高阶单元(三角形)情况下,如果由于种种原因忽略了X3(或X2Y)项,则为了保持模式的几何等向性,也不应包括Y3(或XY2)项.;第(2)条选择插值多项式阶次的要求,是使多项式中所含的总项数等于单元的结点自由度数目.如三结点三角形平面单元,插值多项式选为:;(1)位移函数中必须含有反映刚体运动的项数. 多项式形式的常数项即体现这一刚体位移.每个单元的位移一般总是包含两部分,一是由本单元形变引起,另一部分是与本单元形变无关的,即刚体位移,它是由其他单元发生形变而连带引起的,如悬臂梁自由端处本身形变小,位移主要是连带引起. (2)位移函数应反映单元的常应变,即位移函数的导数中必须有常数项存在. 　 当单元尺寸无限缩小时,单元应变将趋近于常量,因此单元位移函数中应包括常应变项平面应力和空间应力中,应变是位移的一阶导数,常应变即要求位移函数含有一次项　　 (3)位移函数必须保证在相邻单元的接触面上应变是有限的. 　 在有限单元法中,按位移(即按最小