发明名称 一种钢筋混凝土梁截面弯矩曲率关系识别的方法
摘要 本发明公开了一种钢筋混凝土梁截面弯矩曲率关系识别的方法,该方法利用受弯构件实测的宏观位移来识别其截面弯矩-曲率关系,即利用各级荷载下的荷载-位移关系,结合共轭梁方法,通过最小二乘法识别受弯构件的截面弯矩与曲率的关系,继而有益于分析非线性受弯构件力学性能。由于采用本发明所述的钢筋混凝土梁截面弯矩曲率关系识别的方法,具有以下有益效果:通过准确的识别受弯构件的截面弯矩曲率关系,能够根据工程结构形式准确地挑选相应材料,对构件进行合理设计。
申请公布号 CN102914473A 申请公布日期 2013.02.06
申请号 CN201210375018.0 申请日期 2012.09.27
申请人 大连民族学院 发明人 张青霞;覃丽坤
分类号 G01N3/20(2006.01)I 主分类号 G01N3/20(2006.01)I
代理机构 大连东方专利代理有限责任公司 21212 代理人 姜玉蓉;李洪福
主权项 1.一种钢筋混凝土梁截面弯矩曲率关系识别的方法,其特征在于,包括以下步骤:S1、将需要测量截面弯矩曲率关系的梁l的两端铰接在刚性支座上,选取梁l截面底部的A、B、C三点,在A点处设置位移计D1,在B点处设置位移计D2,在C点处设置位移计D3,采用二集中力三分点加载方式对梁l逐级加荷载,左荷载距离梁l左端的距离为l<sub>b</sub>,右荷载距离梁l右端的距离为l<sub>b</sub>;S2、对梁l分级加荷载P,求取梁l的截面弯矩M,其中,初始荷载为零,第i级荷载级对应的荷载为P<sub>i</sub>,在荷载P<sub>i</sub>的作用下,第j点的实测位移,即实测挠度为<img file="FDA00002223585900011.GIF" wi="73" he="54" />(j=1,2,3),梁l的截面弯矩为公式一<img file="FDA00002223585900012.GIF" wi="841" he="255" />x是所关注的截面距离梁左段的距离,通过位移计D1,位移计D2和位移计D3分别测得前i级荷载作用下的三点的位移,获取前i级荷载-位移关系曲线;S3、根据前i级荷载-位移关系曲线识别出的截面弯矩M小于P<sub>i</sub>l<sub>b</sub>时对应的截面弯矩M与曲率<img file="FDA00002223585900013.GIF" wi="143" he="71" />的关系,k级荷载P<sub>k</sub>作用下,梁l跨中截面弯矩为<img file="FDA00002223585900014.GIF" wi="252" he="67" />对应截面曲率为公式二<img file="FDA00002223585900015.GIF" wi="543" he="84" />(k=1,2,...,i);S4、当两个相邻荷载级对应的截面弯矩<img file="FDA00002223585900016.GIF" wi="86" he="67" />和<img file="FDA00002223585900017.GIF" wi="100" he="67" />间的弯矩-曲率关系为线性时,根据已测的前i级弯矩-曲率关系,以及测量的第i+1级荷载的位移,识别第i+1级荷载作用下对应的各截面曲率;S5、对梁l加第i+1级荷载P<sub>i+1</sub>时,根据公式一求取梁l的截面弯矩M<sub>i+1</sub>(x),梁l的跨中截面弯矩为<img file="FDA00002223585900018.GIF" wi="142" he="73" />当截面弯矩M<sub>i+1</sub>(x)小于等于<img file="FDA00002223585900019.GIF" wi="96" he="73" />时,对计算得到的弯矩曲率关系<img file="FDA000022235859000110.GIF" wi="191" he="73" />(k=1,2,...,i)进行线性插值,获得截面曲率;截面的弯矩大于<img file="FDA000022235859000111.GIF" wi="86" he="67" />时,设定梁l弯矩在<img file="FDA000022235859000112.GIF" wi="296" he="95" />之间时截面曲率为线性变化,斜率为k<sub>i+1</sub>,得到第i+1级荷载P<sub>i+1</sub>下对应的截面曲率为公式三<img file="FDA000022235859000113.GIF" wi="1574" he="205" />给定任意值k<sub>i+1</sub>,根据公式一和公式三,确定在第i+1级荷载P<sub>i+1</sub>作用下梁l的截面弯矩和截面曲率,然后根据共轭梁法求得梁l的挠度δ(x),为公式四<img file="FDA00002223585900021.GIF" wi="700" he="234" />z为变量,表示所在截面距离梁左段支座的距离;R表示共轭梁的支座反力;在第i+1级荷载P<sub>i+1</sub>作用下,梁l截面底部三点的挠度记为δ<sub>1,i+1</sub>(k<sub>i+1</sub>),δ<sub>2,i+1</sub>(k<sub>i+1</sub>),δ<sub>3,i+1</sub>(k<sub>i+1</sub>),δ<sub>j,i+1</sub>(k<sub>i+1</sub>)(j=1,2,3)为关于k<sub>i+1</sub>的一次函数,将该函数记为公式五δ<sub>j,i+1</sub>(k<sub>i+1</sub>)=a<sub>j1</sub>k<sub>i+1</sub>+a<sub>j2</sub>(j=1,2,3);分别将<img file="FDA00002223585900022.GIF" wi="96" he="100" />赋值为0,将<img file="FDA00002223585900023.GIF" wi="79" he="77" />赋值为1,将<img file="FDA00002223585900024.GIF" wi="113" he="95" />和<img file="FDA00002223585900025.GIF" wi="111" he="92" />带入公式三,求取梁l的截面曲率,将k<sub>i+1</sub>带入公式四,求取公式五中的系数a<sub>j1</sub>和a<sub>j2</sub>为公式六<img file="FDA00002223585900026.GIF" wi="451" he="208" />(j=1,2,3)求取在第i+1级荷载Pi+1作用下梁l在第j点处的挠度,公式七<img file="FDA00002223585900027.GIF" wi="1080" he="104" />(j=1,2,3);S6、通过最小化目标函数公式八,使δ<sub>j,i+1</sub>(k<sub>i+1</sub>)与实测挠度<img file="FDA00002223585900028.GIF" wi="87" he="59" />最接近,求取k<sub>i+1</sub>,公式八<img file="FDA00002223585900029.GIF" wi="939" he="117" />(j=1,2,3),<img file="FDA000022235859000210.GIF" wi="391" he="57" />为梁l在第i+1级荷载P<sub>i+1</sub>作用下的A,B,C三点的实测位移;S7、利用最小二乘法公式九k<sub>i+1</sub>=(A<sup>T</sup>A)A<sup>T</sup>b,<img file="FDA000022235859000211.GIF" wi="451" he="313" /><img file="FDA000022235859000212.GIF" wi="448" he="313" />求取系数k<sub>i+1</sub>;S8、根据k<sub>i+1</sub>,以及第i+1级荷载P<sub>i+1</sub>作用下对应的截面曲率,根据公式十求取对应的跨中截面曲率,公式十为<img file="FDA000022235859000213.GIF" wi="576" he="81" />S9、循环上述方法,求取各荷载级的荷载下对应的截面曲率,获取梁l的弯矩曲率关系获得简支梁截面的弯矩曲率关系<img file="FDA000022235859000214.GIF" wi="267" he="89" />(i=1,2,...,n) 。
地址 116600 辽宁省大连市开发区辽河西路18号