带肋筏板的无单元和有限元耦合方法.docVIP

带肋筏板的无单元和有限元耦合方法 张建辉，余晓雅 （河北大学机械与建筑工程学院，河北 保定 071002） 摘 要：本文将无单元法与传统的有限元法相耦合，建立了双参数地基上带肋筏板的一种耦合数值分析方法。其中，双参数地基上筏板由无单元法模拟，肋梁由梁的有限元方法模拟。本文方法可以方便的求得带肋筏板的挠度，肋梁内力与筏板内力，实例计算表明了其结果的正确性。 关键词：筏板 肋梁 带肋筏板 无单元法 引 言 筏板基础是高层建筑广泛采用的基础型式。为了加强筏板基础刚度，通常在筏板底部或上部添加肋梁，称之为带肋筏板基础。而有关带肋筏板基础与地基土相互作用的计算分析等问题是工程界倍受关注的问题。目前，关于带肋筏板的计算方法相对较少，通常有两种类型[1，2]，其一，将肋梁与筏板看成整体，建立刚度方程，再利用有限差分法或有限元等数值方法求解；其二，将带肋筏板分成筏板和肋梁两个部分：将肋梁视为交叉梁系，由有限元方法求解，筏板的分析是以地基上板的分析为基础 本文将带肋筏板分成筏板和肋梁两个部分，肋梁与筏板在节点上视为铰接，建立肋梁与筏板的耦合模型。采用筏板无单元方法[3]建立筏板的刚度矩阵，有限元法计算肋梁的刚度矩阵，根据虚功原理，建立了双参数弹性地基上带肋筏板的无单元和有限元耦合计算方法。并通过实例计算对这一方法进行验证，结果表明，该方法是合理可行的。 2. 肋梁有限元法 本文把带肋筏板的肋梁部分看作交叉的连续梁，如图1示。运用已经发展比较成熟完善的有限元法分析，通过最小势能原理，推导单元的刚度矩阵和等效节点荷载列阵以及结构的整体刚度矩阵[4]。 图1X轴的角位移θx、绕Y轴的角位移θy和沿Z轴方向的线位移wz；每个结点有三个结点力分量：绕X轴的集中力偶Mx、绕Y轴的集中力偶My和沿Z轴方向的集中力Pz。则交叉梁单元的每个端点有对应的三个杆端位移参数和三个杆端力参数，如图2、图3示。 图图 根据交叉梁系的虚功方程可导出其刚度方程， （1) 其中，为交叉梁结构刚度矩阵，是交叉梁结构的节点荷载列阵，是节点位移列阵， （2) 为整体坐标系下节点位移与局部坐标系下单元杆端位移间的变换矩阵，为局部坐标系下的交叉梁单元刚度矩阵。 3. 肋梁有限元法 无单元法（Element-free Method）采用移动最小二乘法（Moving Least-Square Method）所产生的光滑插值函数近似场函数，为得到离散的线性方程组仅需要对节点和边界条件进行描述，节点可以随机分布，并与积分网格无关。文献[5]首次将无单元法用于板弯曲问题中，成功地建立了弹性地基上筏板的无单元法。 在无单元法中，求解域场函数可被近似地表示为[6] (3) 其中为已知维完备多项式基函数向量；为维未知系数向量，可根据滑动最小二乘方法由节点的场函数值求得。 在板弯曲问题中，取为完备的二次多项式向量，MLS插值函数即可满足解收敛的完备性要求[双参数板]。设坐标面为板的中面，，取维多项式基函数向量。 根据滑动最小二乘（MLS）插值方法，可得到求解系数的线性方程组，即 (4) 其中，、分别为和矩阵，为节点位移参数向量，为节点影响域（的区域）上的权函数，是该节点处的场函数值，为其影响域包含点的节点数。 将式（4）代入式（3）得， (5) 其中为基于滑动最小二乘法的形函数， (6) 设双参数地基上边界自由的筏板承受面荷载，根据（3）式，其挠度函数可近似表示为 （7） 其中为形函数向量，为节点挠度向量。 根据双参数地基上筏板的虚功方程，得到筏板与地基共同作用的基本方程[3]： （8） 其中、分别为弹性板、地基的刚度矩阵，为节点荷载向量，即 ，， 其中，，其中，为表征地基土模型的两个弹性常数，由试验确定。 4. 筏板无单元与肋梁有限元耦合方法 4.1带肋筏板模型 带肋筏板模型见图4，以挠度作为耦合方程的基本位移未知量，将带肋筏板分为筏板和肋梁两部分，两者之间在节点处简化为铰接[1]（梁单元节点一定作为筏板节点），板与肋梁之间只传递节点力，不传递力矩 。 图4 带肋筏板模型 4.2耦合方程的建立与计算 设带肋筏板置于双参数地基上，边界自由，承受上部结构的荷载。根据带肋筏板的虚功方程，可建立带肋