一种桁架结构刚度和轴力的结构识别方法技术

本发明专利技术公开了一种基于MATLAB‑ABAQUS交互条件下关于桁架结构刚度和轴力同时识别以评估结构损伤及预测剩余承载力的技术，属于结构损伤识别领域，主要包括四个步骤：建立钢桁架损伤结构的有限元模型；修改inp文件并建立有限元模板文件；输入真实的待识别量并生成inp文件，计算并获得相应节点加速度和杆件应变；模拟估计的待识别量响应；使用目标函数不断优化待识别量直到满足收敛条件，输出最终结果。本发明专利技术的主要优势在于：(1)识别钢桁架某一杆件刚度的同时，也可识别同一杆件的轴力；(2)传感器位置没有特殊要求，需识别轴力的杆件不必安装传感器；(3)传感器的数量需求较少。因此，本发明专利技术的应用前景较为广阔。

A structural identification method of stiffness and axial force of truss structure

【技术实现步骤摘要】
一种桁架结构刚度和轴力的结构识别方法
本专利技术是一种基于MATLAB-ABAQUS交互条件下的关于桁架结构刚度和轴力同时识别以评估结构损伤及预测剩余承载力的技术，隶属于结构损伤识别领域。
技术介绍
桁架结构损伤识别一直是结构健康监测领域的前沿问题，对钢桁架桥梁结构或大型空间网架结构开展健康监测并评估极具必要性。此种结构单元的主要受力形式为轴向拉压，且杆件轴向力的大小是判断杆件是否正常工作的关键指标，一旦某一杆单元所受轴力超过其承载能力，此杆件有可能已经发生损伤。因此，正确识别杆件轴力为桁架结构的健康状况监测提供了至关重要的预警信息。目前桁架的结构，如大型钢桁架桥梁、大型体育场馆顶部网架等，特点是杆件多、节点多、单元尺寸大。从经济性考虑，结构健康监测无法为所有的杆件单元布设传感器；从使用寿命考虑，现有传感器大多远比建筑结构的使用寿命短。关键杆件上用于监测其响应的传感器发生损坏或因使用时间过长而造成较大测量精度误差且又难以更换时，转而利用安装在其他杆件上的传感器监测此处损伤情况是一个新途径。因此，用少量传感器监测整体结构损伤，或用正常工作的传感器监测结构关键部位的损伤，是桁架结构健康监测领域亟待解决的新问题之一。桁架结构相较于其他结构体系，其工作性能不仅依赖于其刚度，也与轴力相关。而当桁架结构受到损伤时，其刚度的变化必然会同时引起轴力的重分布，所以，两者的相互关系对结构单独的刚度识别造成了极大的影响。因此，同时识别桁架结构关键杆件的刚度和轴力的技术方案是解决这一问题的关键。
技术实现思路
<br>为建立桁架结构刚度和轴力同时识别过程，本部分将专利技术主要步骤陈述如下：步骤1：利用有限元分析计算软件ABAQUS建立待识别结构模型，提交计算后生成此模型的同名inp后缀文件。在不打开ABAQUSCAE界面的情况下可直接调用此inp文件进行计算，生成的结果可直接用于模型的后处理。调用inp文件进行计算在ABAQUS的后续处理过程中的实用性显著优于反复生成模型进行计算，且inp文件提供了一个MATLAB软件与ABAQUS软件共同工作的接口，为后续调用MATLAB实施结构损伤和轴力识别奠定了基础；步骤2：以可编辑形式打开inp文件，将文件中的待识别量以未知数形式表示。此待识别量可设置为模型中某杆的弹性模量，或是模型节点上的外荷载等，经修改后即可改变对应的几何或力学参数。随后，打开MATLAB建立一函数文件并写入步骤一中生成的inp文件，将待识别量作为函数文件中的自变量。在待识别量的搜索空间内随机取值，即可生成与其对应的新的inp文件，记此函数文件为有限元模板文件；步骤3：定义结构中某一杆件刚度和某一杆件(可异于前一杆件)所受轴力为待识别量，以期同时识别某一杆件的刚度和轴力。将待识别量实际值代入有限元模板文件，生成inp文件且提交至ABAQUS计算，从计算结果fil文件中提取桁架结构真实的应变和加速度响应；步骤4：将待识别量估计值代入有限元模板文件，此估计值可限于一定范围内，此范围可根据实际情况确定。在待识别量的搜索上下限内生成随机值，生成对应的inp文件，将此inp文件提交至ABAQUS中即可计算出对应桁架结构的应变和加速度响应；步骤5：设置一个目标函数用于评价步骤4中待识别量与真实值的误差，当步骤3和4生成的响应小于允许误差值时，可判断此时步骤4中桁架结构待识别量即实际值，将此待识别量的估计值作为最终结果输出。通过ABAQUS中的有限元计算和MATLAB中的优化计算交互循环即可实现刚度和轴力的同时识别，如图1所示，使用MATLAB中自带优化工具箱中的粒子群(PSO)算法实施优化计算，具体步骤如图2所示。步骤6：输出刚度和轴力识别结果，此结果即为识别结果。本专利技术的主要优势：(1)识别钢桁架的某一杆件刚度的同时，也可识别同一杆件的轴力，与传统方法分两次识别杆件的刚度和轴力相比效率更高；(2)传感器位置没有确切要求，要求识别的杆件上不必布设传感器。对某些不便安装传感器的杆件，或是在需要识别的杆件上安装的传感器因某种原因而损坏时，可借助结构其他杆件上安装的传感器识别此杆件的刚度和轴力；(3)传感器数量需求较少，用于识别刚度和轴力的加速度计和应变片均只布置了1个，即使本专利技术识别的结构相对简单，但其原理也适用于大型桁架结构或者空间网架结构的刚度与轴力识别，可以预计相应的传感器需求的数量也较少。因此，本专利技术的应用前景较为广阔。附图说明图1：本专利技术所提出的刚度和轴力同时识别的主要流程图。图2：粒子群优化(PSO)算法图。图3：二维桁架示意图。图4：识别结果图。具体实施方式下面结合附图及实例对本专利技术做进一步说明。步骤1：建立模型。桁架结构由10个节点和16根杆组成，包括上弦杆4根，下弦杆4根，腹杆4根，长度均为500mm，斜腹杆4根长度为707mm，截面均采用直径21.3mm，壁厚1.5mm的空心钢管。在ABAQUS中建立如图3所示的二维桁架模型，截面弹性模量和泊松比分别设置为210GPa和0.3，杆件使用beam模型。为考虑杆件的刚度削弱，可另建立一个截面属性，此截面的几何属性与原结构其他截面的几何属性相同，而刚度考虑为其他截面的0.5倍，以此表征杆件损伤。为结构划分网格时，由于施加荷载阶段ABAQUS软件无法为结构体系每一杆件都施加各自不同的轴力，因此，本专利技术考虑为每一杆件都划分两个等长的单元，用于在杆件中点处施加一个轴向的外荷载。该荷载虽与常规轴力有较大差别，但在桁架结构中，仅仅考虑小变形的前提下，施加在杆件中部的外荷载仅会对杆件截面上的轴力造成影响，而不会有其他额外的改变。因此，施加在杆件中部的外荷载即可等效为各杆件的轴力，实现通过改变施加在某一杆件的外荷载数值来等效地任意改变某一杆件所受的轴力；步骤2：修改inp文件并建立有限元模板文件。上一步骤在ABAQUS建立模型成功提交计算后，软件便自动生成对应模型的inp文件。修改inp文件的具体操作为：删除OUTPUTREQUESTS字段和EndStep字段之间的所有命令，并用*FILEFORMAT,ASCII*NODEFILEA*ELFILESE代替，其中A表示输出所选节点的加速度，SE表示输出所选杆件的截面应变。同时，将此inp文件每一行命令首尾分别添加fprintf(fileID,'和\n')；两个字段，即可得到有限元模板文件主体。然后，在inp文件中搜索需要替换的刚度和轴力，将两数值各乘以一个系数因子，把系数因子形成的一行两列向量作为有限元模板文件的自变量。使用系数因子而不直接使用待识别量的优点在于，第一，系数因子直接反映了识别的误差，比生成刚度和轴力更清晰直观；第二，利用系数因子可直接对比两个识别结果的精确度，便于权衡和改进目标函数。此文件即为生成所需刚度与轴力inp文件的有限元模板文件；步骤3：输入待识别量并生成inp文件，计算并获得真实的加速度和应变。将桁本文档来自技高网...

【技术保护点】
1.一种桁架结构刚度和应变的结构识别方法，其特征在于：/n步骤1：建立模型；/n在ABAQUS软件中建立二维桁架模型，并设置相应的属性；/n另建立一个杆件截面属性，此杆件截面属性的几何属性与原杆件其他截面的几何属性相同，而刚度小于为其他截面，以此表征杆件损伤；/n为杆件划分网格时，为每一杆件都划分两个等长的单元，用于在杆件中点处施加一个轴向的外荷载；施加在杆件中部的外荷载即可等效为各杆件的轴力，通过改变施加在某一杆件的外荷载数值来等效地任意改变某一杆件所受的轴力；/n步骤2：修改inp文件并制作有限元模板文件；/n步骤1在ABAQUS建立模型成功提交计算后，ABAQUS软件自动生成对应模型的inp文件；修改inp文件，将此修改后的inp文件每一行命令首尾分别添加fprintf(fileID,'和\n')；两个字段得到有限元模板文件主体；然后，在inp文件中搜索需要替换的刚度和轴力，将两数值各乘以一个系数因子，把系数因子形成的一行两列向量作为有限元模板文件的自变量；/n步骤3：输入关键量并生成inp文件，计算并获得真实的加速度和应变；/n将桁架某一杆件的真实[EI，N]输入有限元模板文件，EI和N分别表示杆件刚度和轴力，然后生成相应inp文件，与建模时生成的inp文件区别在于有限元模板文件生成的inp文件其输出设置经修改后，只能生成对应结果的二进制文件，其文件后缀为fil；MATLAB读取此fil文件，获得相应加速度和应变数据，然后将ABAQUS的结果文件fil中包含的ASCII格式的信息汇编为具有一行的字符串Rec，调用abaqus2matlab工具箱中的Rec103函数和Rec29函数从字符串Rec中提取相应节点的加速度和杆件的截面应变；/n步骤4：模拟估计的待识别量响应；/n将该搜索上下限即系数因子取值范围设置在一定范围内，随后在搜索区间中随机产生数值输入有限元模板文件，得某一刚度和轴力条件下加速度和应变；/n步骤5：设置一个目标函数，将步骤3和步骤4两次计算的数据输入目标函数，通过ABAQUS中的有限元计算和MATLAB中的优化计算交互循环即可计算估计刚度和轴力相应的误差，若计算过程中该值达到收敛条件即认为估计刚度和轴力即为实际刚度和轴力；/n步骤6：输出最终结果，此结果即为识别结果。/n...

【技术特征摘要】
1.一种桁架结构刚度和应变的结构识别方法，其特征在于：
步骤1：建立模型；
在ABAQUS软件中建立二维桁架模型，并设置相应的属性；
另建立一个杆件截面属性，此杆件截面属性的几何属性与原杆件其他截面的几何属性相同，而刚度小于为其他截面，以此表征杆件损伤；
为杆件划分网格时，为每一杆件都划分两个等长的单元，用于在杆件中点处施加一个轴向的外荷载；施加在杆件中部的外荷载即可等效为各杆件的轴力，通过改变施加在某一杆件的外荷载数值来等效地任意改变某一杆件所受的轴力；
步骤2：修改inp文件并制作有限元模板文件；
步骤1在ABAQUS建立模型成功提交计算后，ABAQUS软件自动生成对应模型的inp文件；修改inp文件，将此修改后的inp文件每一行命令首尾分别添加fprintf(fileID,'和\n')；两个字段得到有限元模板文件主体；然后，在inp文件中搜索需要替换的刚度和轴力，将两数值各乘以一个系数因子，把系数因子形成的一行两列向量作为有限元模板文件的自变量；
步骤3：输入关键量并生成inp文件，计算并获得真实的加速度和应变；
将桁架某一杆件的真实[EI，N]输入有限元模板文件，EI和N分别表示杆件刚度和轴力，然后生成相应inp文件，与建模时生成的inp文件区别在于有限元模板文件生成的inp文件其输出设置经修改后，只能生成对应结果的二进制文件，其文件后缀为fil；MATLAB读取此fil文件，获得相应加速度和应变数据，然后将ABAQUS的结果文件fil中包含的ASCII格式的信息汇编为具有一行的字符串Rec，调用abaqus2matlab工具箱中的Rec103函数和Rec29函数从字符串Rec中提取相应...

【专利技术属性】
技术研发人员：王小娟，陈峰，周宏元，倪萍禾，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：北京;11