复杂形态钢结构设计培训班

首页 结构设计-公众号 思考 | 有限元计算出来就是对的?如何保证有限元分析的正确性和合理性

思考 | 有限元计算出来就是对的?如何保证有限元分析的正确性和合理性

转自:结构设计-公众号


来源:Nikko结构空间

引言

最近遇上一个逆作法项目施工评审,支护结构为1m厚的地连墙。由于出土方量巨大,且工期紧张,施工单位提出把原来逐层开挖的方案调整为跳仓一次开挖至承台面。同时,为了最大限度减小对出土流线的影响,约11.5m高的地连墙不设置水平支撑,通过局部预留反压土堆减小地连墙的临空长度。

网络图片非本文描述项目

在查看施工方案和相应计算书后,小编给出的主要意见是“计算假定有误,支护结构验算偏不安全,建议局部增设临时水平支撑”。收到意见后,施工单位的Y工就来电沟通了,以下是对话摘录。

小编:你们计算假定不合理,得出结果有问题。

Y工:我们是采用了专业有限元软件PLAXIS进行分析的,结果是通过精细化有限元计算得出。

小编:哪模型本构参数是如何取值呢?这些数值是否适用于本项目岩土情况?

Y工:本构模型采用经典的莫尔库仑模型,参数取值后续补充给你。这是目前最常用的模型。

小编:项目靠近江边,地下常水位很高,这个模型里面是如何考虑呢?

Y工:(沉默30秒……)这个我再查查,好像没有考虑。

小编:有限元模型边界条件不对啊!地连墙顶怎么能按固接考虑呢?顶板压根受不了这弯矩?

Y工:(沉默30秒……)这个回去查查。

小编:计算模型中地连墙长度取值也有问题。留一个小土堆(3m高,6m长)作为反压措施,就起到约束作用?

Y工:(沉默中……)

与Y工沟通后,让我想起读研时导师的一句话:“有限元就是一个rubbish in,rubbish out的东西,没有经验证的有限元结果都必须谨慎判别”。持“有限元算出来就是对的”观点的工程师应该不少,但真是这样子吗?


什么是有限元?


要搞清楚“有限元算出来就是对的”这个问题,首先得知道什么是有限元。这里不涉及高深的基础理论(如变分原理、加权余量法等,主要小编也不太懂),仅分享小编的理解,供各位同行指正。
我们常说的“有限元”其实是一种数值计算方法,工作思路是将分析对象划分为有限个互不重叠单元,然后借助不同原理建立方程并求解。1956年,Clough等人将刚架分析的位移法应用到弹性力学的平面问题中,人们一般认为这是有限单元法的诞生时间。其实结构力学中的矩阵位移法就是有限单元法的雏形。
有限单元法最早应用于结构力学,逐渐发展到流体力学等领域。随着计算机的发展,一大批有限元软件涌现,其中部分软件开发商会以“操作简单”、“一键生成”作为卖点,以此降低软件使用门槛,扩大其市场占有份额。然而,对于我们这个涉及建筑安全的结构专业来说,某些情况下,这对工程师来说不一定是个好事情。
这里从王勖成老师《有限单元法》一书中摘录了几段文字供各位细读。

 有限单元法的其中一个特点是:“只要原问题的数学模型(基本变量、基本方程、求解域和边界条件等)是正确的,同时求解的有限元方程算法是稳定、可靠的,则随着单元数目增加,或者随着单元自由度数目增加及插值函数阶次提高,有限元解的近似程度将不断被改进。”

  应用有限元软件进行工程分析时,分析者需要具有必要的力学知识工程经验。同时,还要求分析者对于有限元的基本原理和离散方法,常用单元形式和求解方法的特点和应用条件有清晰的理解。
  分析者需对计算结果作出分析和评估,决定是否需要修改有限元模型和计算方案进行重新分析;是否需要修改力学数学模型;是否需要修改原设计方案。
 前面提及的第2、3点是成功应用有限元软件,特别是大型通用软件进行工程分析的前提条件。
由以上几段文字可知,如果结构工程师们不具备一定的理论知识、工程经验、判别能力,就使用有限元软件进行计算分析,这就很有可能得到无效的计算结果。


影响有限元结果的一些主要因素


接王勖成老师上述的话,结合小编一些使用和开发有限元的经验,建议结构工程师在使用有限元软件时,注意以下几个影响计算结果的主要因素。

● 本构模型及其参数

在有限元分析中,本构模型用于模拟材料受到外界激励时的反应(主要为应力-应变关系)。根据材料试验,人们建立本构数学模型,它们一般是由带参数的数学式来表示。大部分本构模型需使用者对其参数进行定义,有的参数要求在某个范围内取值,有的会建议经验值,而有些参数则仅有数学形式意义,没有物理意义

在弹性分析中,本构参数仅涉及一些可测量的物理量(弹模、泊松比、密度等),它们是可以准确定义。但当考虑材料非线性时,本构模型中部分不可测量或仅有数学形式意义的参数就无法准确定义了,且某些参数还会影响材料进入非线性状态后的反应。
当涉及到材料非线性行为时,材料本构参数少则四、五个,多则十来个,其中部分关键参数影响十分大。这么多参数需要使用者逐一确定其合理的取值,大大地增加了我们准确模拟材料响应的难度。如下图所示,同一个本构模型,仅某个关键参数取值不同时,应力-应变滞回曲线形状就完全改变。
Mander约束混凝土本构模型
(约束因子取值不同时应力-应变曲线有不同形态)

Giuffre-Menegotto-Pinto钢筋本构模型

(形状参数R取值不同时应力-应变曲线会有不同形态)

● 连接与接触

无论是连接定义,还是接触定义,它们都是反映节点之间、单元之间的相互作用。采用不同的连接方式(自由度的约束和释放)或接触方式(不同接触算法,如罚函数法、拉格朗日乘子法等),主要影响到应力或内力在分析对象中的大小分布。

例如混凝土与钢筋之间粘结滑移问题,采用不同粘结滑移关系定义,直接影响钢筋混凝土构件进入非线性状态后的滞回性能表现。又例如主次梁刚铰接定义,直接影响到主次梁的弯矩和扭矩的大小与分布。

ABAQUS中采用硬接触算法模拟两刚体相对运动

● 荷载与边界条件

荷载作用大小、分布形式对计算结果的影响是显然易见的,而边界条件的影响同样不能忽视。边界条件是指分析对象与外界之间的相互关系。在结构工程中,主要涉及结构支座是固定、铰支、滑动等

以本文开篇讲述的地连墙验算为例,若计算分析中考虑反压土堆的约束作用,地连墙临空长度将减少1/4。无论是承载力,还是变形、裂缝均能满足要求。但这个反压小土堆(高3m,长6m)对地连墙根本起不了有效的侧限作用。如何合理地、正确地定义分析对象的边界条件,也是影响有限元结果准确性的关键要素之一。

● 单元类型与网格划分

单元类型与网格划分主要影响有限元计算精度,单元类型的选取与网格划分的疏密是有一定联系的。在综合考虑分析对象几何属性、计算精度、计算成本等因素基础上,选取合理的单元类型和合适的单元大小是有限元分析的另一个关键点。

我们需要根据分析对象的物理属性,选用不同维度、不同形式的分析单元,如实体单元、壳单元、膜单元、梁单元等。但不管是什么单元,它们本质区别主要在于构造的单元插值函数不同。如果使用者不了解所使用的单元原理和属性,可能会产生一些不合理的结果。

如ABAQUS中的一阶完全积分单元,对于以弯曲变形为主的分析对象,采用这种单元会出现剪力自锁的数值问题,从而导致分析对象过刚,计算的变形偏小。此时,需选用高阶完全积分单元或者减缩积分单元来解决。

网格划分大小与单元的插值函数有密切关系。一般来说,由插值函数精度高低来确定单元划分疏密程度。如OpenSees中的基于刚度法或柔度法的纤维截面杆系单元。由于刚度法采用的是位移插值函数,只能近似描述变形沿杆长分布情况,因此需对杆件进行细分,以此获得较多单元,减少插值函数与实际变形之间的误差。而柔度法采用的是力插值函数,它能准确描述出力沿杆长分布情况,因此仅需小量单元就可满足计算精度要求。

某两端固定框架梁的竖向变形-荷载曲线

以上是小编经验之谈,实际影响有限元计算结果还有很多要素,诸如求解方法、收敛性等等。由此可见,要想通过有限元计算得到相对合理的结果并不容易,必须严谨把控各个环节的参数取值。但这还不够,如本构模型中部分参数仅有数学意义,其取值合理性无法校验。又如单元尺寸大小、网格疏密,需要进行敏感性分析,评估哪个尺度的计算结果趋于稳定。这时就需对有限元结果进行校验。


有限元计算结果校验


有限元计算结果的校验一般有两种方式:与理论计算值对比与试验结果对比。前者多出现在一些通用有限元软件的verification文档中,如ABAQUS、ANSYS。后者多用于进行参数化分析,校验基础数值模型的正确性。除了这两种主流验证方式外,还有就是通过盲测比赛检验各种商用、自编的有限元程序的可靠性。

● 软件自带的验证文件

像ABAQUS、ANSYS这类大型通用有限元软件一般都是提供一些相对简单构件的理论值与计算值的对比。然而实际分析对象往往是十分复杂的,尤其关于抗震性能研究。

这里特别介绍下SeismoStruct软件的验证文档。它不但提供了与理论值的验证结果,还提供了与试验结果的验证对比。它的verification文件内容涵括了混凝土结构、钢结构、桥梁、静力分析、动力分析等等。

官网SeismoStruct软件介绍


SeismoStruct自带的verification文档目录

● 盲测比赛

盲测比赛(blind prediction contest)源于日本防灾可研所下属的振动台试验机构E-Defense(专注于大比例尺振动台试验),大概从2005年开始每年都会举办,接受全球的科研团队、工程师、专业软件商、自编程序人员参加。它以先模拟分析、后进行实验的方式检验有限元软件的准确性。随后国内也开始参照该模式举办了多个盲测比赛,其中以清华大学和中国建筑学会抗震防灾分会联合组织得比较好。

E-Defense网页介绍


国内盲测比赛网页介绍


结语


前说了那么多,主要是希望不管是结构工程师们,还是在校学生们,在使用有限元软件时必须慎之又慎。它是个强大的计算工具,但并没有我们想象中那么容易操控。除了必备的基础理论知识,还需对计算结果合理性进行判断。有条件的话,还应对有限元模型进行校验。这才能最大限度保证我们由有限元分析、研究和设计出来的东西是合理的、正确的。


——-END——-

来源:Nikko结构空间如有侵权请联系我们。


转自:结构设计-公众号

公众号:MyStructure

“转自:结构设计-公众号”公众号,致力于推广建筑结构新理论、新技术、新工艺、新材料及大量工程实例,介绍建筑结构知识和行业资讯,搭建一个建筑结构行业交流的平台。

  

本文来自网络,不代表钢构人的立场,转载请注明出处。搜索工程类文章,就用钢构人网站。 https://www.ganggouren.com/2022/01/abef59302f/

钢结构地图

上一篇
下一篇

作者: ganggouren

为您推荐

发表回复

您的电子邮箱地址不会被公开。 必填项已用 * 标注

联系我们

联系我们

17717621528

在线咨询: QQ交谈

邮箱: 1356745727@qq.com

工作时间:周一至周五,9:00-17:30,节假日休息
关注微信
微信扫一扫关注我们

微信扫一扫关注我们

关注微博
返回顶部