大体积混凝土浇筑过程中的温度场应力场分析第四章.docVIP

4. 大体积混凝土应力场模型的建立及其有限元离散 4.1热应力和热弹性力学基本方程 所谓热应力，是指由于温度变化在弹性体内所产生的应力。热弹性力学的基本方程包括平衡方程、几何方程、和物理方程，平衡方程和几何方程同等温情况一样。设不计体力，平衡方程为 （4-1-1） 几何方程为 （4-1-2） 物理方程和等温情况有所不同。因为在变温情况下，弹性体的应变分量应由两部分叠加而成：其一，是由于自由膨胀引起的应力分量；其二，在热膨胀时由于弹性体内各部分之间的相互约束所引起的应变分量，它们和热应力之间服从胡克定律。因此，变温情况下的物理方程为 （4-1-3） 也可用应变分量来表示应力分量，即 （4-1-4） 如果将式（4-2）代入，则应力分量又可通过位移分量表示，有 （4-1-5） 式中： 在物体的边界处，还必须满足面力为零的应力边界条件 （4-1-6） 和位移边界条件 （4-1-7） 将式（4-1-5）代入式（4-1-1）中，得到变温情况下一位移表示的平衡微分方程： （4-1-8） 如果将式（4-1-5）代入应力边界条件（4-1-6），得到变温情况下以位移表示的应力边界条件 （4-1-9） 在一定的位移边界条件下，弹性体内由于变温而引起的位移，等于等温情况下受假想体力 和假想的法向面力 作用时的位移。在应力边界条件（4-1-9）和位移边界条件（4-1-7）下求得了方程（4-1-8）的解后，利用式（4-1-5）可求得应力分量。 4.2大体积混凝土施工期温度应力的计算 大体积混凝土施工期的温度应力，包括外约束即地基约束引起的外约束应力和混凝土内部对外部引起的内约束应力，目前施工单位多使用下式来计算温度应力: （4-2-1） 式中:——混凝土的温度（包括收缩）应力 ——混凝土龄期t的弹性模量; ——混凝土的线膨胀系数，取; ——混凝土的最大综合温度差; ——考虑徐变影响的松弛系数; ——混凝土的外约束系数: ——混凝土的泊松比。 实际上这只是计算外约束应力的公式。在浇筑后的早期，由于混凝土的内外温差较大，内约束应力往往要大于外约束应力[23]。某高层建筑板厚为2.5m的筏式基础的温度应力验算表明:浇筑混凝土后的6—9天，外约束应力为0.15—0.25，而内约束应力为0.336—0.611，内约束应力远大于外约束应力，因而仅通过公式（4-1-1）验算温度应力是不安全的。还应考虑内约束应力，计算内约束应力的公式为: （4-2-2） 式中:——混凝土块体的平均温度; ——混凝土的表面温度。 用（4-2-1）式计算外约束应力时，只取截面中均匀降温差。即取截面中部水化热温升作为全截面的均匀温度，而事实上任何大体积混凝土温度分布都是不均匀的。（4-2-2）式计算内约束应力时，假设混凝土不受外界约束，上下散热条件相同，且假定温度沿混凝土的厚度为二次抛物线分布。其计算假设显然与实际条件及温度分布情况有出入[24]。 按照（4-2-1）、（4-2-2）式把大体积混凝土的温度应力分为外约束应力与内约束应力分别计算，虽然具有简化计算、易于运用等特点，但由于公式推导过程做了如上的假设，其计算结果与实际值会有差别，并且它不能反映出大体积混凝土各处的温度应力分布。由于大体积混凝土的温度应力是由分布在整个混凝土体积的温度场所致，且由于混凝土的弹性徐变性质等因素常引起应力“异常”现象。譬如，按热弹性应力理论计算，应该出现压应力的地方，有时却出现了拉应力;反之，在拉应力区却出现了压应力[25]。 综上所述，用公式（4-2-1）、（4-2-2）计算温度应力，但对实际工程仅实现“了解温度应力，及时采取有效措施”是不够的。 4.3大体积混凝土应力场的有限元分析 4.3.1连续介质的离散化 我们求解的对象是二维连续介质，它本来具有无限个自由度，如图4-2所示，把它划分成有限个三角形单元，全部单元都在它们的角点上互相连接起来，这些点称为结点，以结点位移作为未知量，在每一个结点建立两个平衡方程，联立解之，即可求出各结点的位移，由结点位移可求出各单元的应力。在单元足够小时，这些单元的应力与原结构的应力很接近。可见有限单元法的核心是把本来具有无限个自由度的连续介质离散成有限个自由度的结构，就可以在计算机上求解[26]。 图 4-2 混凝土平面单元划分 图 4-2 混凝土平面单元划分 4.3.2位移函数 典型的三角形单元e，三个结点i，j，m，按逆时针方向排列。每