CN115828762B - 受弯梁材料本构参数模型修正及判别方法 - Google Patents
受弯梁材料本构参数模型修正及判别方法 Download PDFInfo
- Publication number
- CN115828762B CN115828762B CN202211641562.5A CN202211641562A CN115828762B CN 115828762 B CN115828762 B CN 115828762B CN 202211641562 A CN202211641562 A CN 202211641562A CN 115828762 B CN115828762 B CN 115828762B
- Authority
- CN
- China
- Prior art keywords
- model
- section
- constitutive
- corrected
- parameters
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 239000000463 material Substances 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 85
- 230000004044 response Effects 0.000 claims abstract description 188
- 238000012937 correction Methods 0.000 claims abstract description 173
- 238000005452 bending Methods 0.000 claims abstract description 45
- 230000006835 compression Effects 0.000 claims abstract description 35
- 238000007906 compression Methods 0.000 claims abstract description 35
- 238000005457 optimization Methods 0.000 claims abstract description 21
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 6
- 238000004364 calculation method Methods 0.000 claims description 118
- 238000012360 testing method Methods 0.000 claims description 75
- 238000013461 design Methods 0.000 claims description 22
- 239000004567 concrete Substances 0.000 claims description 19
- 238000012986 modification Methods 0.000 claims description 12
- 230000004048 modification Effects 0.000 claims description 12
- 229910000831 Steel Inorganic materials 0.000 claims description 9
- 239000010959 steel Substances 0.000 claims description 9
- 230000007935 neutral effect Effects 0.000 claims description 6
- 239000011150 reinforced concrete Substances 0.000 claims description 4
- 238000012067 mathematical method Methods 0.000 claims description 2
- 238000013178 mathematical model Methods 0.000 claims description 2
- 239000011513 prestressed concrete Substances 0.000 claims description 2
- 238000003062 neural network model Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 description 53
- 238000004458 analytical method Methods 0.000 description 38
- 230000002787 reinforcement Effects 0.000 description 15
- 230000003068 static effect Effects 0.000 description 9
- 238000001514 detection method Methods 0.000 description 8
- 238000003825 pressing Methods 0.000 description 7
- 238000010206 sensitivity analysis Methods 0.000 description 7
- 238000004088 simulation Methods 0.000 description 6
- 230000003014 reinforcing effect Effects 0.000 description 5
- 230000001174 ascending effect Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 239000011374 ultra-high-performance concrete Substances 0.000 description 4
- 238000009795 derivation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000011241 protective layer Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000035945 sensitivity Effects 0.000 description 3
- 229910001294 Reinforcing steel Inorganic materials 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 2
- 238000013401 experimental design Methods 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 238000000540 analysis of variance Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000012850 discrimination method Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种受弯梁材料本构参数模型修正及判别方法。步骤如下：选择合适的材料本构方程；选取材料本构方程中的待修正参数，并确定合理的待修正参数取值区间；基于平截面假定，由截面应力平衡列轴力平衡方程，求解截面受压区高度；判断轴力平衡方程中是否存在线性相关项；由力矩平衡列弯矩平衡方程，计算截面抗弯刚度，截面受压区高度为常数时，判别式EIB≠0，能进行材料本构参数模型修正；截面受压区高度不为常数时，判断式与参数个数相关；采用理论模型或者响应面模型通过优化算法对待修正参数进行识别，完成模型修正。本发明提出一种受弯梁材料本构参数模型修正及判别方法，为非线性模型修正方法在实际工程中的应用提供了理论基础。
Description
技术领域
本发明属于结构模型修正领域，涉及梁结构本构模型参数修正技术，具体地说是涉及一种受弯梁材料本构参数模型修正及判别方法。
背景技术
模型修正是结构健康监测系统的重要组成部分，模型修正技术能有效减小模型与实际结构的误差，在土木工程领域中应用广泛，也是目前比较热门的研究方向。已有模型修正研究多以构件刚度参数为修正量，从材料本构模型的层面对本构参数进行修正的研究甚少。本发明开展混凝土梁桥静力响应数值模拟，建立本构关系参数与结构静力响应值(应变、刚度、挠度)的响应面，提出修正材料本构关系参数的方法，从本构层面进行有限元模型修正，识别本构参数，同时提出能否进行非线性本构模型参数修正的判据。
发明内容
针对受弯梁能否进行非线性模型修正，识别材料本构参数的问题，本发明提出一种受弯梁材料本构参数模型修正及判别方法。
本发明所述受弯梁材料本构参数模型修正及判别方法，步骤如下：
(1)根据工程经验，选择合适的材料本构方程；
(2)选取材料本构方程中的待修正参数，并确定合理的待修正参数取值区间，a)仅一个待修正参数时，可进行材料本构参数模型修正，转步骤(6)；b)两个及以上待修正参数时，继续下述步骤；
(3)基于平截面假定，假设截面曲率为由截面应力平衡，列轴力平衡方程，求解截面受压区高度表示关于自变量的函数；
(4)以待修正参数为变量，判断轴力平衡方程中是否存在线性相关项，若存在，则不能进行材料本构参数模型修正，转步骤(7)；若不存在，继续下述步骤；
(5)由力矩平衡，列弯矩平衡方程，求得弯矩M，计算截面抗弯刚度若对多个不同截面曲率k＝1,2,3...，a)截面受压区高度为常数时：若判别式EIB＝0，不能进行材料本构参数模型修正，转步骤(7)；判别式EIB≠0，能进行材料本构参数模型修正，继续下述步骤；
b)截面受压区高度不为常数时：两个待修正参数，判别式EIA＞0.4％，三个待修正参数，判别式EIR＞0.7％，四个待修正参数，判别式EIR＞1.5％，能进行材料本构参数模型修正，继续下述步骤；否则，不能进行材料本构参数模型修正，转步骤(7)；
(6)采用理论模型或者响应面模型，通过优化算法对待修正参数进行识别，完成材料本构参数的修正，模型修正时目标响应值的个数不小于待修正参数的个数；
(7)模型修正结束；
具体的，步骤(5)中，判别式分别为：
其中，带下标的EI表示对不同的值，通过理论模型计算得到的截面抗弯刚度，EI为响应值，n表示待修正参数的数量，下标k表示顺序号，k＝1,2,3...,n-1；
EIk：表示在选取的多组EI响应值中，按照曲率由小到大排序，曲率为时对应计算得到的截面抗弯刚度EIk；
EIk+1：表示在选取的多EI组响应值中，按照曲率由小到大排序，曲率为时对应计算得到的截面抗弯刚度EIk+1，与EIk相邻；
EImax：表示在选取的多EI组响应值中，按照曲率由小到大排序，最大的曲率值对应计算得到的响应值；
EImin：表示在选取的多组EI响应值中，按照曲率由小到大排序，最小的曲率值对应计算得到的响应值。
具体的，步骤(1)中材料本构方程可根据相应的设计规范选择，如混凝土、普通钢筋、预应力钢筋本构模型可参照《JTG 3362-2018公路钢筋混凝土及预应力混凝土桥涵设计规范》选用。
具体的，步骤(2)中材料本构方程中的待修正参数取值区间要包含待修正参数的实际值。
具体的，步骤(3)中，轴力平衡方程与梁截面形状、组成截面材料的种类数量、材料本构方程相关，以单一材料矩形截面梁为例：
∑N＝0；
其中，∑N为截面轴力之和，N表示轴力，B、H分别为矩形截面的宽度和高度，z为截面受压区高度，分别为材料的受压应力、受拉应力，由材料本构方程计算，为距离截面中性轴y位置的应变，为截面曲率，y为距离截面中性轴的位置；
假设截面曲率已知，可求得应变通过材料本构方程，可求得受压应力受拉应力轴力平衡方程中只有截面受压区高度z未知，故可求解截面受压区高度表示关于自变量的函数。
具体的，步骤(5)中，以单一材料矩形截面梁为例：
Mi＝∑M；
其中，∑M为截面弯矩之和，M表示弯矩，Mi为与截面曲率对应的总弯矩。
具体的，步骤(6)中，理论模型指采用数学方法建立的精确理论数学模型，响应面模型指由试验设计方法生成样本，拟合得到响应与待修正参数之间的响应面模型，是近似模型，与理论模型相比，有一定的误差，响应面模型可采用多项式模型、克里格模型、径向基函数、BP神经网络等模型。
具体的，步骤(6)中，模型修正的目标响应值取截面抗弯刚度EI、梁的挠度w或者梁的挠度倒数w-1等响应。
具体的，步骤(6)中，模型修正可采用单目标或多目标优化算法。
本发明基于理论模型和响应面模型，以截面抗弯刚度EI、梁的挠度w或者梁的挠度倒数w-1等响应为目标响应，进行非线性响应面模型修正，分析了多种材料本构方程，总结了受弯梁能否进行非线性模型修正的判别方法，为非线性模型修正方法在实际工程中的应用提供了理论基础。
附图说明
图1是本发明方法的流程框图。
图2是本发明方法的理论模型的参数修正步骤。
图3是本发明方法的矩形截面梁受弯分析。
图4是本发明方法的响应面模型的参数修正步骤。
图5是本发明方法的简支梁示意图(单位:mm)。
图6是本发明方法的梁截面示意图(单位:mm)。
图7是本发明方法的受弯状态下材料本构模型参数修正的判别方法示意图。
图8是本发明方法实施例一的结构设计布置示意图(单位:mm)。
图9是本发明方法实施例二的结构设计布置示意图(单位:mm)。
具体实施方式
下面结合附图及实施例对本发明做进一步说明，下面的描述涉及附图时，除非另有表示，不同附图的相同数字表示相同或相似的要素。
本发明所述受弯梁材料本构参数模型修正及判别方法，实现流程如图1，具体内容如下：
一、理论模型及响应面模型的材料本构参数修正方法
1)梁结构基于理论模型的材料本构参数修正方法(EI响应)
本节算例以UHPC本构(上升段)为例进行数值模拟分析，建立简支梁数值模型，以EI为响应采用基于理论模型的参数修正方法对本构参数进行修正。
数值算例采用简支梁(受弯结构)模型，通过数值算例对本构参数修正方法进行详细介绍，应用两种方法实现相关参数修正，并比较两种方法的修正效果。简支梁为不配筋的超高性能混凝土梁，矩形截面(B×H：200mm×400mm)，计算跨径为4100mm。
为分析方便算例采用本构模型的上升段(极限应变之前)作为分析本构模型，如下式：
式中式中fc、ft分别为混凝土的抗压、抗拉强度值；ε0为混凝土压应力刚达到fc时混凝土的应变值，即峰值应变；表示梁结构受弯时梁轴线变形后的曲率，yi表示距中性轴的距离。
通过上式可知：受压本构主要由三个参数决定，分别是fc、εc、εc0；受压本构主要也是由三个参数决定，分别是ft、εt、εt0；在本算例中将fc和ft选取为待修正参数。
①方法介绍及修正步骤
本节方法建立响应与本构参数之间的理论模型，利用数学中的最优化问题去解决工程应用中的参数修正问题，同时需要运用MATLAB软件辅助计算来实现。
以简支梁桥结构作为分析对象，该方法首先根据结构的材料本构关系选取待修正参数，然后根据平截面假定及结构的受力状态推导出截面的轴力平衡及弯矩平衡方程，以EI为目标响应值建立理论模型，然后根据结构设计已知条件及平衡方程运用MATLAB软件编程或非线性有限元分析软件对结构进行响应计算，最后根据理论模型运用MATLAB最优化问题求解办法，求解出待修正参数的修正值并进行误差分析。具体修正步骤如图2所示。
②理论分析
根据平截面假定以及本构模型，可得到本构上升段时截面的应力、应变分布图，如图3。
根据平衡条件，设截面受压区高度为z，距截面中性轴距离为y的位置，由静力关系可以写出两个平衡方程：
将本构方程代入式(2)计算可得轴力平衡方程：
将本构方程代入式(3)计算可得弯矩平衡方程：
通过式(4)可知受压区高度满足以下函数关系：
根据材料力学可知有以下关系，式中EI表示截面抗弯刚度：
那么通过式(4)、(5)和(7)可知EI需满足以下函数关系：
由以上推导，根据最优化问题(无约束最优化问题)，当EI0i为已知目标响应值时，EIi是根据以上推导得到的函数计算值，在参数给定的参数区间内求解参数fc和ft的修正值，可以建立一个关于EIi的数学理论模型，如式(9)，目标函数采用平方和最小，当目标函数最小时说明EIi与EI0i无限接近，而此时EIi是对应的参数值fc和ft通过理论计算得到，则此时的fc和ft就是拟合得到的参数的修正值，此拟合步骤要应用MATLAB实现相关计算。
式中：EI0i为已知目标响应值，EIi为计算响应值。
同时可以观察最后输出目标函数的数值检验模型的精度，当目标函数f越接近于0说明建立的理论模型的吻合度越高。
③试验数值计算
本试验设计两组理论试验，A组试验两参数的取值接近，模拟拉压本构参数接近的材料；B组试验两参数的取值相差比较大，模拟拉压本构参数相差大的材料(即模拟UHPC材料的拉压强度性能)。先求解两组EI响应，然后根据响应值结合理论模型求解参数的修正值：
(1)A组试验：以为已知量，峰值应变ε0＝0.002，给fc、ft一个包含理论值的区间为fc∈[10，40]、ft∈[10，30]，通过两平衡方程进行计算分别计算出响应值EI1、EI2和EI3。
(2)B组试验：以为已知量，峰值应变ε0＝0.002，给fc、ft一个包含理论值的区间为fc∈[120，270]、ft∈[10，30]，通过两平衡方程进行计算分别计算出响应值EI1、EI2和EI3。
A组中取一组理论值fc＝28MPa、ft＝21MPa，B组中取一组理论值fc＝150MPa、ft＝13MPa，应用MATLAB计算，截面受压区高度z的计算结果见表1，弯矩M的计算结果见表2，响应值EIi的计算结果见表3。
表1两组试验在不同值时计算得出的受压区高度z(单位：mm)
表2两组试验在不同值时计算得出的弯矩M(单位：kN·m)
表3两组试验在不同值时计算得出的EIi目标响应值(单位：N·mm2)
④参数灵敏度分析
灵敏度分析采用局部灵敏度分析，对某个参数变化增幅1％，比较响应值的变化，分析判断其显著性。两组试验中各参数的局部灵敏度计算结果如下表4和表5，表中响应为EIi(单位为：N·mm2)。
表4试验A组fc＝28MPa，ft＝21MPa参数的局部灵敏度分析
表5试验B组fc＝150MPa，ft＝13MPa参数的局部灵敏度分析
通过灵敏度分析，由表4、表5可知：在两组试验中，ft的灵敏度均大于fc的灵敏度，初步判断在应用该方法修正参数时，ft参数比fc参数对EI响应计算的影响更显著，而ft的修正效果比fc的修正效果较好。
⑤优化求解及效果
根据建立的理论模型，通过前面计算的响应值优化拟合计算得到fc、ft的修正解，比较修正解与理论解的误差。输入目标响应EIi(三个或任意两个)，计算结果见表6。
表6两组试验参数修正的计算结果表
两组试验理论模型修正参数时输出的目标函数皆符合f＝0；计算结果显示，经过理论模型对参数修正，参数修正结果非常好，参数相对误差均为0，修正值为精确解。综合说明，对梁结构的基于理论模型的本构材料参数修正方法是有效的，能够实现对材料参数进行精确修正，得到精确解。
2)梁结构基于响应面模型的材料本构参数修正方法(EI响应)
本节同样以UHPC本构(上升段)为例，建立简支梁数值模型，以EI作为响应值，采用基于响应面模型的参数修正方法对本构参数进行修正，详细介绍该方法并阐述该方法的有效性及修正效果。
①方法介绍及修正步骤
响应面法在模型修正中应用得相当广泛，本节该修正方法就是建立响应与本构参数之间的响应面模型，实现解决数值算例及工程应用中的参数修正问题。该方法首先根据设计结构的材料本构关系选取待修正参数，同时确定待修正参数的取值区间，然后选择试验设计方法进行样本设计，接着根据已有条件对结构计算及试验样本计算，再然后确定目标函数、以EI或w为目标响应值建立响应面模型，同时对模型进行评估，最后根据响应面模型优化求解出待修正参数的修正值并进行误差分析。结构计算需要运用MATLAB软件或有限元分析软件，求解修正值过程需要运用MATLAB软件和Design-Expert 12软件辅助计算来实现。具体修正步骤如图4所示。
②试验设计
根据本构模型，将fc和ft选取为待修正参数，同样设计两组试验，采用CCD实验设计方法进行试验设计，由两水平两因素、采用一个中心点，生成9组样本序列(通过Design-Expert12软件实现)：
(1)A组试验样本：理论值取fc＝28MPa、ft＝21MPa，给fc、ft一个包含理论值的区间为fc∈[10，40]、ft∈[10，30]；参数水平变量表见表7所示、CCD实验设计表见表8所示。
(2)B组试验样本，理论值取fc＝150MPa、ft＝13MPa，给fc、ft一个包含理论值的区间为fc∈[120，270]、ft∈[10，30]。参数水平变量表见表9所示、CCD实验设计表见表10所示。
表7A组参数水平变量表
表8A组CCD实验设计表
表9B组参数水平变量表
表10B组CCD实验设计表
③试验数值(EI响应)计算
根据所拟定简支梁模型的截面尺寸及跨径，以 为已知量，试验设计以EI作为响应值建立响应面模型做参数修正分析，对9组样本序列通过本构关系计算结构的各响应值，EIi响应值的计算可用MATLAB软件或有限元分析软件进行计算，设计两组理论试验的9组样本序列EIi响应计算结果如表11、表12。
表11A组试验样本值不同时的EIi响应值(N·mm2)
表12B组试验样本值不同时的EIi响应值(N·mm2)
④参数灵敏度分析
对试验样本的响应值做方差分析(应用Design-Expert 12软件)，通过软件得到F检测法各参数一次项、二次项和交叉项的F值和P值，进而判断各参数的灵敏度及显著性。
两组试验样本序列，EIi作为响应值时，各参数项的F、P值如表13、表14所示。
通过方差分析，由表13、表14可知：在A组试验中，fc与ft的显著性相当，初步判断在应用该方法修正参数fc、ft时，两参数修正效果相当；在B组试验中，ft的显著性大于fc的显著性，初步判断在应用该方法修正参数fc、ft时，ft的修正效果比fc的修正效果较好。
表13A组试验样本EIi为响应值各参数项的F、P值
表14B组试验样本EIi为响应值各参数项的F、P值
⑤刚度响应(EI)与参数的响应面拟合
响应面模型本节所应用的就是二次多项式(含交叉项)模型来拟合，形式如式(10)，式中y为响应特征值，xi(i＝1,2,…n)为待修正参数，α为待定系数，待定系数应用最小二乘法回归分析求得。
A组试验中fc、ft与EIi的二次多项式响应面方程如下：
B组试验中fc、ft与EIi的二次多项式响应面方程如下：
对拟合的响应面模型需要进行精度检测，精度检测用来评价模型与实际有限元模型的接近程度，进行精度检测时一般采用R-squared和Adjusted R-squared两个指标。R-squared(即R2)描述的变量之间相关程度的指标量；Adjusted R-squared(即Adj R2)是考虑了变量数目的影响根据样本数量和变量数量进行调整后的R2指标量。两个判定系数取值范围都为[0，1]，越接近1表明拟合的计算模型在试验样本空间内越能描述变量之间的相关程度，具体计算公式如下。
式中，yR(i)：表示模型(响应面模型)的某次计算值；y(i)：表示理论(有限元分析)的某次计算值；y：表示理论(有限元分析)计算结果的平均值；NN：表示试验样本空间的数量；n：表示变量(待修正参数)的数量。
A、B两组试验样本空间R2和Adj R2的计算结果见下表15响应面精度检测表：
表15A、B两组试验样本空间的响应面精度检测表
由上表可知，两组试验的R2和Adj R2的值都接近于1，说明拟合的响应面模型在试验样本空间内能代替有限元模型，能较准确的描述响应值和参数之间的相关关系。
⑥刚度响应(EI)优化求解及效果
根据建立的响应面模型，通过计算理论组的EI目标响应值，输入目标响应EIi(三个或任意两个)，应用Design-Expert 12软件优化拟合求解计算fc、ft的修正解，比较修正解与理论解的误差。
(1)A组试验样本：理论组取fc＝28MPa、ft＝21MPa，输入计算得到对应的EI响应值，通过优化拟合计算得到修正解，计算结果见表16。
表16A组试验：fc＝28MPa，ft＝21MPa参数修正的计算结果表
表17B组试验：fc＝150MPa，ft＝13MPa参数修正的计算结果表
(2)B组试验样本，理论组取fc＝150MPa、ft＝13MPa，输入计算得到对应的EI响应值，通过优化拟合计算得到修正解，计算结果见表17。
由表16、表17表计算结果显示，经过响应面模型修正法对参数拟合修正，两组试验参数修正结果都较好，相对误差控制在5％以下，其中A组相对误差控制在1％以下。综合说明，对梁结构的基于响应面模型的本构材料参数修正方法是有效的，能够实现对材料参数进行较精确修正，求出修正解。同时，在A组试验中，fc与ft两参数修正效果相当；在B组试验中，ft的修正效果比fc的修正效果较好，两组修正结果对比说明参数的修正效果与参数显著性有直接联系。
二、结构受弯状态下材料本构模型参数修正及判别方法
不同本构组算例采用简支梁(受弯状态)作为分析对象，尺寸拟定同上节，如图5所示，跨径为4100mm。截面都为矩形截面(200mm×400mm)，根据需要分为三种形式：
(1)结构中不配筋，截面如图6(1)所示；
(2)结构只在截面受拉侧配受力筋，采用钢筋类型为HRB335带肋钢筋、公称直径为d1，截面配筋面积为As，保护层厚度为t，截面如图6(2)所示；
(3)结构在截面受拉、受压双侧配受力筋，采用钢筋类型为HRB335带肋钢筋，受拉侧配筋公称直径为d1、截面配筋面积为As1、保护层厚度为t1，受压侧配筋公称直径为d2、截面配筋面积为As2、保护层厚度为t2，截面如图6(3)所示。
1)本构模型数值算例理论分析及响应计算
本构组一
本构组一假设由单种材料组成，其中拉压本构模型完全相同，均为非线性本构，选取两个待修正参数；假设的材料本构模型形式如式(19)，结构为受弯状态则式中
式中
根据平截面假定及平衡条件，设截面受压区高度为z，由前面推导的静力关系式(2)、(3)可得到的两个平衡关系如下：
联立式(20)、(21)进行化简并求解z、EI：
经计算，可得到z的三个解如式(23)，但有效解只有一个即H/2，由此可判断z与值的变化无关，而由式(22)可判断EI会随值得变化而改变。
采用图6(1)中拟定的不配筋截面，梁结构不配筋，选取fc、ft为待修正参数，设计一组理论值：fc＝28MPa、ft＝21MPa，峰值应变ε0＝0.002，对模型进行理论分析计算。
应用MATLAB对理论试验组进行计算分析，在不同值时计算得出有效受压区高度zi、弯矩Mi、EIi响应的计算结果见表18。
表18试验组在不同值时计算得出的zi、Mi、EIi
本构组二
本构组二假设由两种材料组成，两种材料拉压本构模型相同，均为线性本构，选取三个待修正参数；假设的材料本构模型形式如式(24)，结构为受弯状态则式中
式中
根据平截面假定及平衡条件，设截面受压区高度为z，由前面推导的静力关系式可得到的两个平衡关系如下：
对式(25)进行化简，式(26)已是最简形式：
联立式(26)、(27)进行化简并求解z、EI，经计算z为常数，EI线性相关：
采用图6(2)中拟定的拉侧配筋截面，在梁结构受拉侧配2根受力钢筋，带肋钢筋公称直径d＝22mm，截面配筋面积As＝760.265mm2，选取fc、ft、fs为待修正参数，设计一个理论值：fc＝30MPa、ft＝20MPa、fs＝25MPa，峰值应变ε0＝0.002，对模型进行理论分析计算。
应用MATLAB对理论试验组进行计算分析，在不同值时计算得出有效受压区高度zi、弯矩Mi、EIi响应的计算结果见表19。
表19试验组在不同值时计算得出的zi、Mi、EIi
本构组三
本构组三假设由单种材料组成，其中拉压本构模型形式上相同，表达形式中的参数不同，均为非线性本构，选取两个待修正参数；假设的材料本构模型形式如式(29)，结构为受弯状态则式中
式中
根据平截面假定及平衡条件，设截面受压区高度为z，由静力关系推导可得到的两个平衡关系如下：
联立式(30)、(31)进行化简并求解z、EI：
采用图6(1)中拟定的不配筋截面，梁结构不配筋，选取fc、ft为待修正参数，设计一个理论值：fc＝28MPa、ft＝21MPa，峰值应变ε0＝0.002，对模型进行理论分析计算。
应用MATLAB对理论试验组进行计算分析，在不同值时计算得出有效受压区高度zi、弯矩Mi、EIi响应。
本构组四
本构组四假设由单种材料组成，其中拉压本构模型不同，受压为非线性本构、受拉为线性本构，选取两个待修正参数；假设的材料本构模型形式如下式，结构为受弯状态则式中
式中
根据平截面假定及平衡条件，设截面受压区高度为z，由静力关系推导可得到的两个平衡关系同式(4)、(5)，联立进行化简并求解z、EI：
采用图6(1)中拟定的不配筋截面，梁结构不配筋，选取fc、ft为待修正参数，设计一个理论值：fc＝30MPa、ft＝21MPa，峰值应变ε0＝0.002，对模型进行理论分析计算。
应用MATLAB对理论试验组进行计算分析，在不同值时计算得出有效受压区高度zi、弯矩Mi、EIi响应。
本构组五
本构组五假设由两种材料组成，第一种材料拉压本构模型不同，受压为非线性本构、受拉为线性本构，第二种材料以截面受拉侧配筋体现为线性本构，选取三个待修正参数；假设的材料本构模型形式如下式，结构为受弯状态则式中
式中
根据平截面假定及平衡条件，设截面受压区高度为z，由静力关系推导可得到的两个平衡关系如下：
联立两式进行化简并求解z、EI：
采用图6(2)中拟定的拉侧配筋截面，在梁结构受拉侧配2根受力钢筋，带肋钢筋公称直径d＝22mm，截面配筋面积As＝760.265mm2，选取fc、ft、fs为待修正参数，设计一个理论值：fc＝30MPa、ft＝20MPa、fs＝25MPa，峰值应变ε0＝0.002，对模型进行理论分析计算。
应用MATLAB对理论试验组进行计算分析，在不同值时计算得出有效受压区高度zi、弯矩Mi、EIi响应。
本构组六
本构组六假设由两种材料组成，第一种材料拉压本构模型不同、但均为非线性本构，第二种材料以截面拉压两侧配筋体现、均为线性本构、但表达形式中参数不同，选取四个待修正参数；假设的材料本构模型形式如下式，结构为受弯状态则式中
式中
根据平截面假定及平衡条件，设截面受压区高度为z，由静力关系推导可得到的两个平衡关系如下：
联立两式进行化简并求解z、EI：
采用图6(3)中拟定的拉压双侧配筋截面，在梁结构受拉、受压双侧配筋，在受拉侧配2根受力钢筋，带肋钢筋公称直径d1＝22mm，截面配筋面积As1＝760.265mm2；在受压侧配2根受力钢筋，带肋钢筋公称直径d2＝10mm，截面配筋面积As2＝157.08mm2；选取fc、ft、Es1、Es2为待修正参数，设计一个理论值：fc＝26.64MPa、ft＝2.259MPa、Es1＝1.8×105MPa、Es2＝2×105MPa，峰值应变ε0＝0.002，对模型进行理论分析计算。
应用MATLAB对理论试验组进行计算分析，在不同值时计算得出有效受压区高度zi、弯矩Mi、EIi响应。
2)本构模型数值算例不同工况下的参数修正结果
基于理论模型对受弯梁材料本构模型参数进行修正，采用多种判别式对本构模型参数修正进行分析探讨。结合非线性度的概念，采用1-范数对响应进行分析计算得到非线性度判别式，如式(43)所示，式中ψ表示非线性度，N(EI)表示非线性系统的EI响应值，L(EI)表示线性系统的EI响应值。经计算分析发现，用非线性度判别式对修正效果进行判别时，判别效果不明显，无法准确分类判断。同时，还结合范德蒙行列式等提出多种判别式，但判别效果均不明显，因此均未采用。
最终提出三个对EI响应进行分析计算的判别公式，计算方法如式(44)～(46)，式(44)中EIB指标适用于受压区高度为定值情况，式(45)中EIA指标适用于受压区高度不为定值且待修正参数为两个时的情况，式(46)中EIR指标适用于受压区高度不为定值且待修正参数≥3时的情况：
其中，带下标的EI表示对不同的值，通过理论模型计算得到的截面抗弯刚度，EI为响应值，n表示待修正参数的数量，下标k表示顺序号，k＝1,2,3...,n-1。
EIk：表示在选取的多组EI响应值中，按照曲率由小到大排序，曲率为时对应计算得到的截面抗弯刚度EIk。
EIk+1：表示在选取的多组EI响应值中，按照曲率由小到大排序，曲率为时对应计算得到的截面抗弯刚度EIk+1，与EIk相邻。
EImax：表示在选取的多组EI响应值中，按照曲率由小到大排序，最大的曲率值对应计算得到的响应值。
EImin：表示在选取的多组EI响应值中，按照曲率由小到大排序，最小的曲率值对应计算得到的响应值。
本构组一
给fc、ft一个包含理论值的区间为fc∈[10，40]、ft∈[10，40]，对各组本构模型不同测点位置下的工况进行本构参数修正，每个工况分别取两个不同值对应的[EIk,EIk+1]作为目标响应值，并对修正效果进行判别，同时计算相应判别式的指标值，经大量分析计算，选取部分修正效果区别显著的修正组，如下表20所示。
表20部分修正组的效果表
注：表中●表示参数修正效果良好，○表示参数修正效果一般，×表示参数修正效果差。
本构组二
给fc、ft、fs一个包含理论值的区间为fc∈[10，40]、ft∈[10，40]、fs∈[10，40]，对各组本构模型不同测点位置下的工况进行本构参数修正，每个工况分别取三个不同值对应的响应值[EIk,EIk+1,EIk+2]作为目标响应值，并对修正效果进行判别，同时计算相应判别式的指标值，经大量分析计算，选取部分修正效果区别显著的修正组，如下表所示。
表21部分修正组的效果表
注：表中●表示参数修正效果良好，○表示参数修正效果一般，×表示参数修正效果差。
本构组三
给fc、ft一个包含理论值的区间为fc∈[10，40]、ft∈[10，40]，对各组本构模型不同测点位置下的工况进行本构参数修正，每个工况分别取两个不同值对应的响应值[EIk,EIk+1]作为目标响应值，并对修正效果进行判别，同时计算相应判别式的指标值，经大量分析计算，选取部分修正效果区别显著的修正组，如表22所示。
表22部分修正组的效果表
注：表中●表示参数修正效果良好，○表示参数修正效果一般，×表示参数修正效果差。
本构组四
给fc、ft一个包含理论值的区间为fc∈[10，40]、ft∈[10，40]，对各组本构模型不同测点位置下的工况进行本构参数修正，每个工况分别取两个不同值对应的响应值[EIk,EIk+1]作为目标响应值，并对修正效果进行判别，同时计算相应判别式的指标值，经大量分析计算，选取部分修正效果区别显著的修正组，如表23所示。
表23部分修正组的效果表
注：表中●表示参数修正效果良好，○表示参数修正效果一般，×表示参数修正效果差。
本构组五
给fc、ft、fs一个包含理论值的区间为fc∈[10，40]、ft∈[10，40]、fs∈[10，40]，对各组本构模型不同测点位置下的工况进行本构参数修正，每个工况分别取三个不同值对应的响应值[EIk,EIk+1,EIk+2]作为目标响应值，并对修正效果进行判别，同时计算相应判别式的指标值，经大量分析计算，选取部分修正效果区别显著的修正组，如表24所示。
表24部分修正组的效果表
注：表中●表示参数修正效果良好，○表示参数修正效果一般，×表示参数修正效果差。
本构组六
给fc、ft、Es1、Es2一个包含理论值的区间为fc∈[15，30]、ft∈[1.75，2.5]、Es1∈[160000，210000]、Es2∈[160000，210000]，对各组本构模型不同测点位置下的工况进行本构参数修正，每个工况分别取四个不同值对应的[EIk,EIk+1,EIk+2,EIk+3]作为目标响应值，并对修正效果进行判别，同时计算相应判别式的指标值，经大量分析计算，选取部分修正效果区别显著的修正组，如表25所示。
表25部分修正组的效果表
注：表中●表示参数修正效果良好，○表示参数修正效果一般，×表示参数修正效果差。
3)结构受弯状态下材料本构模型参数修正的判别方法
前文列举了多组表达形式不同的本构模型，经过理论推导及分析计算，对各组本构模型不同测点工况下的本构参数进行修正并对修正效果进行评判，分析发现：
通过轴力、弯矩平衡关系计算可知，本构组一和二受压区高度z恒为定值。本构组二无论何种工况修正效果都极差，受压区高度z为恒定值C，且轴力平衡方程中存在线性相关项，进一步EI不随值的变化而变化，因此该本构不能进行本构参数修正。本构组一只要满足EI不取相同的响应值各工况修正效果皆良好，受压区高度z恒等于H/2，不随值的变化而变化，但EI随值的变化而变化，采用判别式(44)进行分析，说明受压区高度z为恒值时，判别式满足EIB≠0，则该本构能够进行本构参数修正。该法则同样适用于本构二。
同样通过两平衡公式计算，其他本构组皆满足受压区高度z不为恒值，但修正效果参差不齐。经大量分析计算，对结构中各材料本构进行理论分析，若平衡方程不存在线性相关项时，且受压区高度z不为恒值时，若该本构可进行本构参数修正，且修正效果良好，则根据判别式(45)、(46)计算需满足以下条件：当待修正参数为两个时，满足判别式(阀值)EIA＞0.4％，如本构组三和本构组四；当待修正参数为三个时，满足判别式(阀值)EIR＞0.7％，如本构组五；当待修正参数为四个时，满足判别式(阀值)EIR＞1.5％，如本构组六。若确定本构能够进行参数修正时，解的特征是唯一解还是双解，需要由轴力、弯矩平衡方程及参数灵敏度进一步确定，如本构组三的修正解为双解。
根据以上分析对其他大量本构组进行理论分析计算并对参数进行修正，发现均能满足以上分析中给定的阀值，因此由以上的算例分析，针对结构受弯状态下材料本构模型能否进行参数修正提出的判别方法，具体流程如图7所示。
具体的，步骤6中，响应面模型还可采用克里格模型、径向基函数、BP神经网络等。
具体的，步骤6中，模型修正的目标响应值还可取梁的挠度w或者梁的挠度倒数w-1等响应。
具体的，步骤6中，模型修正可采用单目标或多目标优化算法。
实施例一：素混凝土梁本构模型参数修正
混凝土强度等级为C45，截面为300mm×450mm的矩形截面，跨径为5100mm，结构构造设计如图8所示。进行EI响应计算时，测点选择曲率为0.1×10-6mm-1、0.2×10-6mm-1、0.4×10-6mm-1时；进行挠度w响应计算时，加载方式为跨中集中力加载(P＝5kN)，测点位置为曲率对应的跨中位置。
本算例模拟素混凝土材料梁结构，先判别该材料本构能否进行参数修正，然后采用修正方法对强度参数(即本构模型参数)进行修正以及评判其修正效果。
步骤1：选择材料本构方程，在本算例中，混凝土受压本构采用Rüsch建议的非线性本构公式，混凝土受拉本构采用规范中建议的非线性本构公式，如式(47)。
式中
步骤2：选择fc和ft两个参数作为待修正参数，以EIi和wi作为响应值，混凝土峰值应变εc0＝0.002、εt0＝0.0001，给fc、ft一个包含理论组的区间为fc∈[15，50]MPa、ft∈[1.50，3.10]MPa。
步骤3：根据平截面假定，设截面受压区高度为z，列轴力平衡方程：
求解受压区高度：
步骤4：由式(48)可以判断待修正参数fc、ft线性无关。
步骤5：由力矩平衡，列弯矩平衡方程：
可求得截面抗弯刚度：
根据C45混凝土拉、压强度标准值为fc＝29.6MPa、ft＝2.51MPa，模拟理论组(模拟假设结构使用若干年后混凝土强度值下降10％的工况)拉、压强度取值为fc＝26.64MPa、ft＝2.259MPa。
应用MATLAB对理论试验组进行计算分析，在不同值时计算得出有效受压区高度zi、弯矩Mi、EIi响应的计算结果见表26。
表26试验组在不同值时计算得出的zi、Mi、EIi
对各组本构模型不同测点位置下的工况进行本构参数修正，每个工况分别取两个不同值对应的响应值[EIk,EIk+1]作为响应值，并对修正效果进行判别，同时计算相应判别式的指标值，经分析计算，选取部分修正效果区别显著的修正组，如下表27所示。
表27部分修正组的效果表
注：表中●表示参数修正效果良好，○表示参数修正效果一般，×表示参数修正效果差。
该本构组属于两参数型且受压区高度z≠C，根据判别方法需满足判别式EIA值大于0.4％，以上分析计算皆满足。而本算例中，对选定的 工况进行EIA值计算，由式(45)计算可得：
在该算例的条件下满足判别式EIA值大于0.4％的要求，所以在该工况下能对本构进行参数修正。
步骤6：模型修正
①理论模型优化
以为已知量，应用MATLAB或有限元分析软件进行结构计算，截面受压区高度z的计算结果见表28。
表28在不同值时计算得出的受压区高度z(单位:mm)
弯矩M的计算结果见表29，响应值EIi的计算结果见表30。
表29在不同值时计算得出的弯矩M(单位:kN·m)
表30在不同值时计算得出的EIi响应值(单位:N·mm2)
根据无约束最优化问题，根据建立的理论模型，通过目标响应值优化拟合计算得到fc、ft的修正解，比较修正解与理论解的误差。输入理论值对应的响应EIi，求解参数修正解，计算结果见表31，修正结果准确。
表31参数修正的计算结果表
②响应面方法修正(w响应)
试验设计及样本序列同前，根据拟定的结构，采用非线性分析程序对结构进行挠度w响应计算。以为已知量求出的Mi值映射到梁对应的测点位置为跨中0.128169L、0.25565L、0.5L处，计算得到三组挠度目标响应，响应计算结果如表32所示。
表32试验样本不同测点的wi响应值(mm)
试验中fc、ft与wi的二次多项式响应面方程形如式(53)，A、B、C、D、E表示各参数项对应的系数，三个响应面方程各参数项对应的系数值如表33所示：
对拟合得到的各响应面进行精度检测，如表34所示。
理论组取fc＝26.64MPa、ft＝2.259MPa，输入计算得到对应的w目标响应值，通过优化拟合计算得到修正解，修正计算结果见表35，计算可知参数ft的修正误差较大有7％。
表33挠度w响应面方程各参数项对应的系数值
表34参数fc、ft与响wi的响应面精度检测表
表35算例w响应参数修正计算结果表
由于直接采用w响应进行参数修正的误差较大，而采用EI响应进行参数修正的效果比较好，根据梁结构挠度w的计算方法式(54)可近似看出w与EI成反比，因此可以尝试着采用w-1作为响应进行参数修正。
在本算例采用挠度w-1响应进行修正计算，验证该猜想，选用跨中0.128169L、0.25565L、0.5L处的三组挠度的-1次方作为响应，响应计算结果如下表36所示：
表36试验样本不同测点的wi -1响应值(mm)
试验中fc、ft与wi -1的二次多项式响应面方程各参数项对应的系数值如表37所示：
表37挠度w-1响应面方程各参数项对应的系数值
对拟合得到的各响应面进行精度检测，如表38所示：
表38参数fc、ft与响wi -1的响应面精度检测表
表39算例w-1响应参数修正计算结果表
理论组取fc＝26.64MPa、ft＝2.259MPa，输入计算得到对应的w-1目标响应值，通过优化拟合计算得到修正解，修正计算结果见表39，通过修正计算可以观察到参数的修正误差都低于5％，参数修正效果较好，修正效果相比于直接用w响应有了较显著提高。
步骤7：模型修正结束。计算表明：采用EI目标响应值基于理论模型及响应面模型的参数修正方法能够对该材料本构模型进行参数修正且修正效果良好，基于理论模型的修正误差可控制在1％以内、基于响应面模型的修正误差可控制在2％以内；采用w、w-1目标响应值基于响应面模型的参数修正方法同样能够对该材料本构模型进行参数修正，但采用w-1响应比直接采用w响应的修正效果好，采用w响应的修正误差存在7％左右，采用w-1响应的修正误差可控制在3％左右。
实施例二：钢筋混凝土梁本构模型参数修正算例
结构拟定为钢筋混凝土简支梁，混凝土强度等级为C45，截面为300mm×450mm的矩形截面，在受拉侧配两根受力筋，采用钢筋类型为HRB335带肋钢筋、公称直径为28mm，配筋面积As＝1231.504mm2，配筋率为9.122％，跨径为5100mm，结构构造设计如图9所示。进行EI响应计算时，测点选择曲率为0.1×10-6mm-1、0.2×10-6mm-1、0.4×10-6mm-1。
步骤1：选择材料本构方程，在本算例中，混凝土本构方程的采用同算例一，钢筋本构方程采用上升段为线性的本构，如式(55)所示：
式中
步骤2：选择fc、ft和Es三个参数作为待修正参数，以EIi作为响应值对本构材料参数进行参数修正，混凝土峰值应变εc0＝0.002、εt0＝0.0001，给fc、ft和Es一个包含理论组的区间为fc∈[15，50]MPa、ft∈[1.50，3.10]MPa、Es∈[170000，220000]MPa。
步骤3：根据平截面假定，设截面受压区高度为z，列轴力平衡方程，并求解z：
步骤4：由式(56)可以判断待修正参数fc、ft和Es线性无关。
步骤5：由力矩平衡，列弯矩平衡方程，并求得截面抗弯刚度：
根据C45混凝土拉、压强度标准值为fc＝29.6MPa、ft＝2.51MPa，HRB335普通钢筋的弹性模量规范值为Es＝2.0×105MPa，模拟理论组(模拟假设结构使用若干年后混凝土强度值和钢筋的弹性模量下降10％的工况)fc、ft和Es取值分别为fc＝26.64MPa、ft＝2.259MPa、Es＝1.8×105MPa。
应用MATLAB对理论试验组进行计算分析，在不同值时计算得出有效受压区高度zi、弯矩Mi、EIi响应的计算结果见表40。
表40试验组在不同值时计算得出的zi、Mi、EIi
对各组本构模型不同测点位置下的工况进行本构参数修正，每个工况分别取三个不同值对应的响应值[EIk,EIk+1,EIk+2]作为目标响应值，并对修正效果进行判别，同时计算相应判别式的指标值，经大量分析计算，选取部分修正效果区别显著的修正组，如表41。
表41部分修正组的效果表
注：表中●表示参数修正效果良好，○表示参数修正效果一般，×表示参数修正效果差。
该本构组属于三参数型且受压区高度z≠C，根据判别方法需满足判别式EIR值大于0.7％，以上分析计算皆满足。而本算例中，对选定的 工况进行EIR值计算，由式(46)计算可得：
在该算例的条件下满足判别式EIR值大于0.7％的要求，所以在该工况下能对本构进行参数修正。
步骤6：模型修正
以为已知量，应用MATLAB或有限元分析软件进行结构计算，截面受压区高度z的计算结果见表42，弯矩M的计算结果见表43，目标响应值EIi的计算结果见表44。
表42在不同值时计算得出的受压区高度z(单位:mm)
表43在不同值时计算得出的弯矩M(单位:kN·m)
表44在不同值时计算得出的EIi响应值(单位:N·mm2)
根据无约束最优化问题，通过理论模型计算的目标响应值优化拟合计算得到fc、ft和Es的修正解，比较修正解与理论解的误差。输入理论值对应的目标响应EIi，求解参数修正解，计算结果见表45，修正结果准确。
表45参数修正的计算结果表
步骤7：模型修正结束。
以上所述仅为本发明的2个实施例，凡依本发明申请专利范围所做的均等变化与修饰，皆属于本发明的涵盖范围。
Claims (8)
1.一种受弯梁材料本构参数模型修正及判别方法，其特征在于包括如下步骤：
（1）根据工程经验，选择合适的材料本构方程；
（2）选取材料本构方程中的待修正参数，并确定合理的待修正参数取值区间，a) 仅一个待修正参数时，可进行材料本构参数模型修正，转步骤（6）；b) 两个及以上待修正参数时，继续下述步骤；
（3）基于平截面假定，假设截面曲率为，由截面应力平衡，列轴力平衡方程，求解截面受压区高度，表示关于自变量的函数；
（4）以待修正参数为变量，判断轴力平衡方程中是否存在线性相关项，若存在，则不能进行材料本构参数模型修正，转步骤（7）；若不存在，继续下述步骤；
（5）由力矩平衡，列弯矩平衡方程，求得弯矩，计算截面抗弯刚度，若对多个不同截面曲率，k=1,2,3...，a) 截面受压区高度为常数时：若判别式，不能进行材料本构参数模型修正，转步骤（7）；判别式，能进行材料本构参数模型修正，继续下述步骤；
b) 截面受压区高度不为常数时：两个待修正参数，判别式，三个待修正参数，判别式，四个待修正参数，判别式，能进行材料本构参数模型修正，继续下述步骤；否则，不能进行材料本构参数模型修正，转步骤（7）；
（6）采用理论模型或者响应面模型，通过优化算法对待修正参数进行识别，完成材料本构参数的修正，模型修正时目标响应值的个数不小于待修正参数的个数；
（7）模型修正结束；
具体的，步骤（5）中，判别式分别为：
；
；
；
其中，带下标的EI表示对不同的值，通过理论模型计算得到的截面抗弯刚度，，EI为响应值，n表示待修正参数的数量，下标k表示顺序号，k=1,2,3..., n-1；
EIk：表示在选取的多组EI响应值中，按照曲率由小到大排序，曲率为时对应计算得到的截面抗弯刚度EIk；
EIk+1：表示在选取的多组EI响应值中，按照曲率由小到大排序，曲率为时对应计算得到的截面抗弯刚度EIk+1，与EIk相邻；
EImax：表示在选取的多组EI响应值中，按照曲率由小到大排序，最大的曲率值对应计算得到的响应值；
EImin：表示在选取的多组EI响应值中，按照曲率由小到大排序，最小的曲率值对应计算得到的响应值。
2.根据权利要求1所述受弯梁材料本构参数模型修正及判别方法，其特征在于：步骤（1）中材料本构方程可根据相应的设计规范选择，混凝土、普通钢筋、预应力钢筋本构模型可参照《JTG 3362-2018公路钢筋混凝土及预应力混凝土桥涵设计规范》选用。
3.根据权利要求1所述受弯梁材料本构参数模型修正及判别方法，其特征在于：步骤（2）中材料本构方程中的待修正参数取值区间要包含待修正参数的实际值。
4.根据权利要求1所述受弯梁材料本构参数模型修正及判别方法，其特征在于：步骤（3）中，轴力平衡方程与梁截面形状、组成截面材料的种类数量、材料本构方程相关，单一材料矩形截面梁：
；
；
；
其中，为截面轴力之和，N表示轴力，B、H分别为矩形截面的宽度和高度，z为截面受压区高度，、分别为材料的受压应力、受拉应力，由材料本构方程计算，为距离截面中性轴y位置的应变，为截面曲率，y为距离截面中性轴的位置；
假设截面曲率已知，可求得应变，通过材料本构方程，可求得受压应力、受拉应力，轴力平衡方程中只有截面受压区高度z未知，故可求解截面受压区高度，表示关于自变量的函数。
5.根据权利要求1所述受弯梁材料本构参数模型修正及判别方法，其特征在于：步骤（5）中，单一材料矩形截面梁：
；
；
其中，为截面弯矩之和，M表示弯矩，Mi为与截面曲率对应的总弯矩。
6.根据权利要求1所述受弯梁材料本构参数模型修正及判别方法，其特征在于：步骤（6）中，理论模型指采用数学方法建立的精确理论数学模型，响应面模型指由试验设计方法生成样本，拟合得到响应与待修正参数之间的响应面模型，是近似模型，与理论模型相比，有一定的误差，响应面模型可采用多项式模型、克里格模型、径向基函数、BP神经网络模型。
7.根据权利要求1所述受弯梁材料本构参数模型修正及判别方法，其特征在于：步骤（6）中，模型修正的目标响应值取截面抗弯刚度EI、梁的挠度w或者梁的挠度倒数w-1响应。
8.根据权利要求1所述受弯梁材料本构参数模型修正及判别方法，其特征在于：步骤（6）中，模型修正可采用单目标或多目标优化算法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211641562.5A CN115828762B (zh) | 2022-12-20 | 2022-12-20 | 受弯梁材料本构参数模型修正及判别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211641562.5A CN115828762B (zh) | 2022-12-20 | 2022-12-20 | 受弯梁材料本构参数模型修正及判别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115828762A CN115828762A (zh) | 2023-03-21 |
| CN115828762B true CN115828762B (zh) | 2024-01-05 |
Family
ID=85517051
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211641562.5A Active CN115828762B (zh) | 2022-12-20 | 2022-12-20 | 受弯梁材料本构参数模型修正及判别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115828762B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117633986B (zh) * | 2023-12-05 | 2024-05-28 | 江苏中碳建筑产业研究院有限公司 | 基于功能梯度材料的复合叠合楼板正截面承载力计算方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008119090A1 (de) * | 2007-03-30 | 2008-10-09 | Trumpf Maschinen Austria Gmbh & Co. Kg. | Verfahren zur festlegung eines einstellparameterwerts einer biegepresse |
| CN102914473A (zh) * | 2012-09-27 | 2013-02-06 | 大连民族学院 | 一种钢筋混凝土梁截面弯矩曲率关系识别的方法 |
| CN106257541A (zh) * | 2016-07-14 | 2016-12-28 | 长安大学 | 一种桥梁有限元模型的修正方法 |
| CN110489924A (zh) * | 2019-08-31 | 2019-11-22 | 湘潭大学 | 基于响应面方法的混凝土梁非线性模型修正方法 |
| CN110929382A (zh) * | 2019-10-28 | 2020-03-27 | 安徽建筑大学 | 基于弹性约束梁应变影响线曲率的桥梁损伤检测方法 |
| CN113283143A (zh) * | 2021-06-09 | 2021-08-20 | 青岛理工大学 | 一种基于动静力数据的叠合梁有限元模型修正方法 |
| CN114547729A (zh) * | 2022-01-19 | 2022-05-27 | 江苏宿淮盐高速公路管理有限公司 | 一种钢-混凝土组合梁抗弯刚度定量识别方法 |
| CN115270239A (zh) * | 2022-06-16 | 2022-11-01 | 浙江工业大学 | 基于动力特性和智能算法响应面法的桥梁可靠性预测方法 |
| CN115270490A (zh) * | 2022-08-08 | 2022-11-01 | 昆明理工大学 | 一种分层充填体损伤动态本构模型构建方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070168081A1 (en) * | 2005-05-06 | 2007-07-19 | Sung-Ho Shin | Analytic integration of tolerances in designing precision interfaces for modular robotics |
-
2022
- 2022-12-20 CN CN202211641562.5A patent/CN115828762B/zh active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2008119090A1 (de) * | 2007-03-30 | 2008-10-09 | Trumpf Maschinen Austria Gmbh & Co. Kg. | Verfahren zur festlegung eines einstellparameterwerts einer biegepresse |
| CN102914473A (zh) * | 2012-09-27 | 2013-02-06 | 大连民族学院 | 一种钢筋混凝土梁截面弯矩曲率关系识别的方法 |
| CN106257541A (zh) * | 2016-07-14 | 2016-12-28 | 长安大学 | 一种桥梁有限元模型的修正方法 |
| CN110489924A (zh) * | 2019-08-31 | 2019-11-22 | 湘潭大学 | 基于响应面方法的混凝土梁非线性模型修正方法 |
| CN110929382A (zh) * | 2019-10-28 | 2020-03-27 | 安徽建筑大学 | 基于弹性约束梁应变影响线曲率的桥梁损伤检测方法 |
| CN113283143A (zh) * | 2021-06-09 | 2021-08-20 | 青岛理工大学 | 一种基于动静力数据的叠合梁有限元模型修正方法 |
| CN114547729A (zh) * | 2022-01-19 | 2022-05-27 | 江苏宿淮盐高速公路管理有限公司 | 一种钢-混凝土组合梁抗弯刚度定量识别方法 |
| CN115270239A (zh) * | 2022-06-16 | 2022-11-01 | 浙江工业大学 | 基于动力特性和智能算法响应面法的桥梁可靠性预测方法 |
| CN115270490A (zh) * | 2022-08-08 | 2022-11-01 | 昆明理工大学 | 一种分层充填体损伤动态本构模型构建方法 |
Non-Patent Citations (7)
| Title |
|---|
| Wei-Xin Ren,et al..Finite element model updating in structural dynamics by using the response surface method.《Engineering Structures》.2010,全文. * |
| Xiuming Yang,et al..Modal Strain Energy-Based Model Updating Method for Damage Identification on Beam-Like Structures.《Journal of Structural Engineering》.2020,全文. * |
| 二次积分迭代法的钢筋混凝土框架非线性全过程分析;吴春颖 等;《河海大学学报（自然科学版）》;全文 * |
| 基于响应面方法的混凝土梁非线性模型修正;简余;《中国优秀硕士学位论文全文数据库 工程科技II辑》;全文 * |
| 波形钢腹板PC组合箱梁桥有限元模型修正方法的试验研究;康志锐;《中国优秀硕士学位论文全文数据库 工程科技II辑》;全文 * |
| 混凝土结构材料非线性本构参数识别算法研究;辛景舟;周建庭;肖阳剑;李晓庆;苏欣;;材料导报(16);全文 * |
| 碳纤维布加固木梁抗弯性能全过程分析;张莉;于艳春;;山西建筑(16);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115828762A (zh) | 2023-03-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lin et al. | Advanced stress-strain model for FRP-confined concrete in square columns | |
| Naderpour et al. | Prediction of FRP-confined compressive strength of concrete using artificial neural networks | |
| Stramandinoli et al. | An efficient tension-stiffening model for nonlinear analysis of reinforced concrete members | |
| Ouladbrahim et al. | Experimental crack identification of API X70 steel pipeline using improved Artificial Neural Networks based on Whale Optimization Algorithm | |
| Feng et al. | Structural behaviour of cold-formed thin-walled short steel channel columns at elevated temperatures. Part 2: Design calculations and numerical analysis | |
| Shahbazian et al. | Application of the Direct Strength Method to local buckling resistance of thin-walled steel members with non-uniform elevated temperatures under axial compression | |
| Petryna et al. | Reliability of reinforced concrete structures under fatigue | |
| CN115828762B (zh) | 受弯梁材料本构参数模型修正及判别方法 | |
| Sayed et al. | Finite element modeling of the shear capacity of RC beams strengthened with FRP sheets by considering different failure modes | |
| Meyer et al. | Modeling of kinematic hardening at large biaxial deformations in pearlitic rail steel | |
| Pimentel et al. | Extended cracked membrane model for the analysis of RC panels | |
| CN110489924A (zh) | 基于响应面方法的混凝土梁非线性模型修正方法 | |
| Ism et al. | Flexural behavior of continuous RC beams strengthened with externally bonded CFRP sheets | |
| CN105243193A (zh) | 一种确定徐变试验棱柱体试件抗压强度折算系数的方法 | |
| Reguera et al. | Optimal design of composite thin-walled beams using simulated annealing | |
| Huang et al. | Validation of a numerical method for predicting shear deformation of reinforced concrete beams | |
| Rojíček et al. | Material model identification from set of experiments and validation by DIC | |
| Wang et al. | Analytical study on reinforced concrete beams strengthened for flexure and shear with composite plates | |
| Kang et al. | Full-range stress–strain relationship and fracture model for laser cladding additively manufactured 316L stainless steel sheets | |
| Miàs et al. | A simplified method to obtain time-dependent curvatures and deflections of concrete members reinforced with FRP bars | |
| Gribniak et al. | Improving efficiency of inverse constitutive analysis of reinforced concrete flexural members | |
| He et al. | A biparametric perturbation method for the Föppl–von Kármán equations of bimodular thin plates | |
| Flood et al. | Rapid analysis of externally reinforced concrete beams using neural networks | |
| Ferreira et al. | Shear strain influence in the service response of FRP reinforced concrete beams | |
| Singh et al. | Methodology for rapid estimation of deflections in two-way reinforced concrete slabs considering cracking |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |