当前位置:研发设计首页 >> 软件世界 >> 数字化设计 >> solidThinking Inspire在端承桩布局优化设计中的应用
solidThinking Inspire在端承桩布局优化设计中的应用
2017-03-31 22:04:53  作者:李涛 王鹏飞 蔡坤   来源:互联网
  •   1 概述   端承桩作为深基础中桩基的一种类型,是指在竖向极限荷载作用下,桩顶荷载全部或主要由桩端阻力承受的桩。沉降量作为影响端承桩基础质量的关键因素之一,在其设计和施工中格外重要。理论和实测结 ...

  1 概述

  端承桩作为深基础中桩基的一种类型,是指在竖向极限荷载作用下,桩顶荷载全部或主要由桩端阻力承受的桩。沉降量作为影响端承桩基础质量的关键因素之一,在其设计和施工中格外重要。理论和实测结果均表明,端承桩基础的沉降主要是桩的弹性变形,基岩的压缩变形很小,所以需要通过控制桩径、桩体材料、桩的布局以及承台厚度等参数来减小承台沉降差和整体沉降。传统设计方法根据《建筑地基基础设计规范》《建筑桩基技术规范》等规范要求可以完成包括端承桩在内的各种桩基础参数的确定和校核,但是往往存在很大的优化潜能,造成材料的浪费和工程量的增加。

  而采用Altair solidThinking Inspire提供的解决方案,可以方便快捷地对端承桩群桩基础进行优化设计。solidThinking Inspire内置自动划分有限元网格功能,使用连续体拓扑优化技术,以OptiStruct为后台拓扑优化求解器,应用变密度法求解策略,对结构上的低效材料进行删减,可以得到更大程度发挥材料性能的结构拓扑方案,为工程节省工时和材料,带来直接的经济效益。

  OptiStruct使用的变密度法是目前应用最为广泛的拓扑优化方法之一,已经在多种商业有限元优化软件中得到实现和认可。变密度法以各向同性材料为基础,通过一定的插值模型人为地建立单元相对密度和材料杨氏模量之间的显式函数关系。其中常用的插值模型主要有:固体各向同性惩罚微结构模型(Solid Isotropic Microstructure with Penalization,简称SIMP)和材料属性的合理近似模型(Rational Approximation of Material Properties,简称RAMP)。通过在材料插值模型中引入惩罚因子,以结构的单元相对密度为设计变量,对中间密度进行惩罚,促使其向实体和空白材料转化,并运用数学规划法或优化准则法求解,把结构的拓扑优化问题转化为不同密度材料的最优分布问题。

  2 拓扑优化模型的建立

  通常由于端承桩的桩身会穿越软弱土层,桩端设置在密实砂层、碎石类土层、中等风化及微风化岩层等承载能力较强的持力层上,桩侧阻力相对于桩端阻力而言很小,可忽略不计。由此可知端承桩的荷载传递机理:荷载通过基础承台传递给桩体,桩体主要通过桩端的端承力将荷载传递给地基土体或岩层,基本上可以忽略桩间土的作用。所以,端承桩计算模型主要由承台、桩体和基岩三部分组成,如图1所示。

  

图1 端承桩计算模型示意图

  图1 端承桩计算模型示意图

  2.1模型概况

  为了更加清晰地介绍利用solidThinking Inspire进行端承桩布局优化的过程,本文参照实际工程案例,采用简化模型进行三维有限元分析,桩位平面布置如图2。基础承台面积为34.5m×21m,厚1m,作用均布荷载250kPa。采用圆形截面C30混凝土桩,桩径1.5m,桩长12m,桩中心距4.5m,承台边缘距边桩中心距离为1.5m。基岩厚8m,长宽均与承台相等。

  

图2 桩位平面布置图

  图2 桩位平面布置图

  2.2模型的建立

  为简化建模流程和节约计算成本,根据计算模型的对称性,采用半结构计算模型建模,断面处使用对称约束控制。solidThinking Inspire提供了方便简洁的建模方式,采用国际单位制(m,kg,N,Pa),以二维草绘图形为起点,通过拉伸、复制、粘贴、移动等操作,可以快速建立端承桩三维几何模型。使用默认的单元尺寸自动划分有限元网格,最小单元尺寸0.090449m,平均单元尺寸0.45225m,将几何模型转化为有限元计算模型。

  2.3材料与属性

  严格来说,混凝土和基岩是存在非线性性质的,但一方面在工作荷载下,桩体和基岩基本上处于弹性变形范围内,其非线性性质表现并不明显;另一方面,即使材料存在非线性,而用线弹性模型的分析并不改变由此得到的某些规律性的认识。因此,为简化计算,本文采用线弹性各向同性材料来模拟混凝土和基岩,结构中的材料参数设置见表1。

  表1 结构中的材料参数

  

表1 结构中的材料参数

  2.4约束和载荷

  由于计算模型根据其对称性采用半结构,并且基岩模型相对于整个岩层而言在结构面垂直方向也存在对称性,所以需在基岩的前、后、左、右、下部和承台断面处施加对称约束,控制结构面的法向位移。承台上表面承受均布荷载,施加大小为250kPa,方向竖直向下的均布压力。

  2.5优化设置

  设定所有桩体为设计空间,并施加竖直向下的单向拔模形状控制,防止桩体出现中间孔洞,如图3所示。以设计空间中桩体材料的伪密度为设计变量,不考虑频率约束和重力,使用默认厚度控制,以全部设计空间中桩体材料的体积剩余率为约束条件,以整体结构刚度最大为目标函数,基于初始设计方案运用拓扑优化技术进行桩位布局二次优化。

  

图3 设置单向拔模形状控制

  图3 设置单向拔模形状控制

  3 拓扑优化结果

  根据半结构模型的对称性,分别依次设定约束条件为全部设计空间体积的80%、70%、65%、60%和50%,运行优化,初始桩位布置(100%)和不同约束条件下的拓扑优化结果桩位布置如图4所示。

  

图4 不同体积分数下的桩位布置图

  图4 不同体积分数下的桩位布置图

  4 校核分析

  根据solidThinking Inspire几何模型,导出iges格式几何文件到HyperMesh中,设置单位卡片与solidThinking Inspire中相同(国际单位制),采用三角形单元多体映射方式划分体网格,单元尺寸0.5,节点总数50251,单元总数84764。使用MAT1材料卡片和PSOLID属性卡片分别创建混凝土和基岩材料,参数设置如表1。创建对称约束和压力荷载并添加到荷载步中,半结构有限元模型如图5所示。

  

图5 HyperMesh中的半结构有限元模型

  图5 HyperMesh中的半结构有限元模型

  根据图4中的桩位布置修改图5中的半结构有限元模型,使用RADIOSS求解器进行静力分析,并通过后处理模块HyperView查看不同体积分数下的承台的竖向位移云图和承台的von Mises应力云图,如图6、图7所示。

  

图6 不同体积分数下承台的竖向位移云图

  图6 不同体积分数下承台的竖向位移云图

  

图7 不同体积分数下承台的von Mises应力云图

  图7 不同体积分数下承台的von Mises应力云图

  5 结果讨论

  观察图6不同体积分数下承台的竖向位移云图可知,当体积分数为100%时,最小沉降出现在承台角部,最大沉降位于承台中部,且沉降较为均匀;当体积分数为80%时,最小沉降仍出现在承台角部,但面积有所扩大,最大沉降集中于承台边处失桩位置,承台中部沉降量居中,均匀性稍有降低;当体积分数为70%、65%和60%时,最小沉降位于承台中部有桩部位,由于角桩的丢失,最大沉降均集中在承台角部,且承台中部沉降均匀性越来越差;当体积分数为50%时,由于失桩过多,最小沉降位于承台对称轴附近有桩位置,最大沉降出现在承台中部,且出现漏斗形沉降,均匀性最差。观察图7不同体积分数下承台的von Mises应力云图可知,在桩顶与承台连接处会出现应力集中,在设计和施工中需要采取相应措施防止结构破坏。

  由RADIOSS校核分析结果可得不同体积分数下基础承台的最大沉降、最小沉降、沉降差和最大von Mise应力,具体数据见表2。并可根据表2中数据绘制最大沉降、最小沉降、沉降差和最大von Mises应力随体积分数约束的变化曲线,如图8、图9所示。

  表2 校核分析结果

  

表2 校核分析结果

  

图8 沉降量变化曲线

  图8 沉降量变化曲线

  

图9 最大von Mises应力变化曲线

  图9 最大von Mises应力变化曲线

  根据图8沉降量变化曲线可知,最小沉降关于体积分数并不敏感,主要表现为桩体和承台的弹性压缩变形;最大沉降和沉降差随体积分数的减小显著增加,并具有相似的变化趋势,最大沉降表现为承台的挠曲变形,当体积分数减少为80%时,最大沉降和沉降差分别为初始设计的2.4倍和3.4倍,当体积分数为70%时,最大沉降和沉降差分别激增为初始设计的5.2倍和7.5倍,此后随体积分数降低虽略有增减,但仍维持在较高水平。同样,根据图9最大von Mises应力变化曲线可知, 随着体积分数的降低,von Mises应力最大值不断增加,主要表现为桩体和承台连接处的应力集中,体积分数为80%时,为初始设计的1.8倍,体积分数为70%时,为初始设计的2.6,此后随体积分数的减小,最大von Mises应力出现震荡,并在 50%时出现激增,结构趋于破坏。

  综合分析以上数据, 对比不同体积分数的桩位布局对承台沉降和最大von Mises应力的影响,认为桩体体积分数为初始设计的80%时得到的桩位布置具有较好的性能指标。此后随着桩体体积分数的减小,最大沉降、沉降差和最大von Mises应力相对于初始设计均有较大的增加。

  6 结论

  随着建筑施工方法的改进和设计技术的提高,对设计方案经济合理性的要求也越来越高。本文在端承桩布局优化设计中,把桩体去留的离散优化问题转化为以桩体材料密度为设计变量的连续拓扑优化问题,通过使用solidThinking Inspire的拓扑优化功能对整个桩体进行优化,得到了不同体积约束下的桩位布置,并借助RADIOSS求解器进行静力分析校核,对比评估得到具有较好结构性能指标的桩位布置方案,为工程设计提供指导和参考。



版权所有:智造网 京ICP证100778号 京公网安备110102003025 虚假新闻举报电话:010-88379107