发明名称 一种基于混合模型的混凝土强度预测方法
摘要 本发明一种基于混合模型的混凝土强度预测方法,先在现场按照标准混凝土强度检测方法对不同混凝土配比进行强度实验,获得多组“水泥x<sub>1</sub>、高炉矿渣粉x<sub>2</sub>、粉煤灰x<sub>3</sub>、水x<sub>4</sub>、减水剂x<sub>5</sub>、粗集料x<sub>6</sub>和细集料x<sub>7</sub>、养护龄期x<sub>8</sub>混泥土组分配比信息‑混凝土强度y”的学样本,并对混合模型中的极限学机、人工神经网路和支持向量机进行训练,通过以最小化相对误差为优化目标确定最佳的极限学机、人工神经网路和支持向量机,在此基础上,基于自适应权重的决策函数可根据三种建模方法的预测值,确定最佳的混凝土预测强度y,用于判断该混凝土组分输入信息是否达到工程设计要求;本发明对三种建模方法的优缺点扬长避短,综合预测效果更好,从而提高对不同实际工况的适应性,即鲁棒性,这对于混凝土快速配合比设计和质量控制具有重要意义。
申请公布号 CN104991051B 申请公布日期 2017.02.22
申请号 CN201510372029.7 申请日期 2015.06.30
申请人 华侨大学 发明人 赖雄鸣;黄河;赖琴芳;王成;张勇;缑锦;言兰
分类号 G01N33/38(2006.01)I 主分类号 G01N33/38(2006.01)I
代理机构 泉州市文华专利代理有限公司 35205 代理人 陈智海
主权项 一种基于混合模型的混凝土强度预测方法,其特征在于包括如下步骤:步骤1、在现场按照标准混凝土强度检测方法对n组不同水泥x<sub>1</sub>、高炉矿渣粉x<sub>2</sub>、粉煤灰x<sub>3</sub>、水x<sub>4</sub>、减水剂x<sub>5</sub>、粗集料x<sub>6</sub>和细集料x<sub>7</sub>、养护龄期x<sub>8</sub>混泥土组分配比信息下进行强度实验,获得混凝土的强度y,其实验样本数据如下:<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mo>{</mo><mo>&lt;</mo><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>3</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>4</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>5</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>6</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>7</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>8</mn></msub><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>,</mo><mi>y</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>&gt;</mo><mo>;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>&lt;</mo><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>3</mn></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>4</mn></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>5</mn></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>6</mn></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>7</mn></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>8</mn></msub><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>,</mo><mi>y</mi><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>&gt;</mo><mo>;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mn>...</mn><mo>;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>&lt;</mo><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>5</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>6</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>7</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>8</mn></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>,</mo><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>&gt;</mo><mo>;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mn>...</mn><mo>;</mo></mrow></mtd></mtr><mtr><mtd><mrow><mo>&lt;</mo><msub><mi>x</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>3</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>4</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>5</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>6</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><msub><mi>x</mi><mn>7</mn></msub><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mi>x</mi><mn>8</mn><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>,</mo><mi>y</mi><mrow><mo>(</mo><mi>n</mi><mo>)</mo></mrow><mo>&gt;</mo><mo>}</mo></mrow></mtd></mtr></mtable><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>n</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001058817850000011.GIF" wi="1654" he="493" /></maths>令x(i)=[x<sub>1</sub>(i),x<sub>2</sub>(i),x<sub>3</sub>(i),x<sub>4</sub>(i),x<sub>5</sub>(i),x<sub>6</sub>(i),x<sub>7</sub>(i),x<sub>8</sub>(i)]表示式(1)中任意一组混凝土配比成分实验数据;步骤2、利用步骤1获得的n组实验样本数据,以式(2)表示的最小化相对误差为优化目标,分别对混合模型中的极限学习机、人工神经网路和支持向量机进行训练,优化目标如下:<maths num="0002"><math><![CDATA[<mrow><mi>m</mi><mi>i</mi><mi>n</mi><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>&lsqb;</mo><mfrac><mrow><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>y</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow><mrow><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mrow></mfrac><mo>&rsqb;</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001058817850000012.GIF" wi="926" he="134" /></maths>式(2)中,y<sub>k</sub>(i),k=1~3,分别为极限学习机、人工神经网路和支持向量机在x(i)的拟合值;在训练过程中,利用常规的优化算法即可获得各个模型中的最佳参数,训练完成后,可建立x(i)→y(i)的3个映射关系,即3个混凝土强度预测模型,分别表示为f<sub>k</sub>(x)=f<sub>k</sub>(x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>,x<sub>4</sub>,x<sub>5</sub>,x<sub>6</sub>,x<sub>7</sub>,x<sub>8</sub>),k=1~3,这里x=[x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>,x<sub>4</sub>,x<sub>5</sub>,x<sub>6</sub>,x<sub>7</sub>,x<sub>8</sub>],同时,对于式(1)中任意一组实验数据{x(i),y(i)},可以得到对应的极限学习机、人工神经网路和支持向量机拟合值,分别为y<sub>k</sub>(i)=f<sub>k</sub>(x(i)),k=1~3;步骤3、将新的不同混凝土配比成分,<img file="FDA0001058817850000021.GIF" wi="621" he="87" />带入上述三个预测模型f<sub>k</sub>(x)=f<sub>k</sub>(x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>,x<sub>4</sub>,x<sub>5</sub>,x<sub>6</sub>,x<sub>7</sub>,x<sub>8</sub>),k=1~3中,得到对应的预测结果,即:<img file="FDA0001058817850000022.GIF" wi="69" he="63" />k=1~3;步骤4、基于自适应决策函数,计算出对应的混凝土强度预测值:<img file="FDA0001058817850000023.GIF" wi="964" he="132" />其中α<sub>k</sub>为响应权重,k=1~3,<img file="FDA0001058817850000024.GIF" wi="310" he="126" />这里为△<sub>k</sub>为拟合误差|y(i)‑y<sub>k</sub>(i)|在x<sup>*</sup>位置附近区域的总误差,i=1~n,具体计算方法如下:<maths num="0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>&Delta;</mi><mi>k</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mo>|</mo><mi>y</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>y</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>|</mo><mo>*</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>d</mi><mo>(</mo><mi>i</mi><mo>)</mo><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><mi>d</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mo>|</mo><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>x</mi><mo>*</mo></msup><mo>|</mo><msup><mo>|</mo><mn>2</mn></msup><mo>=</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>p</mi><mo>=</mo><mn>1</mn></mrow><mn>8</mn></munderover><mo>|</mo><msub><mi>x</mi><mi>p</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>x</mi><mi>p</mi><mo>*</mo></msubsup><msup><mo>|</mo><mn>2</mn></msup></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001058817850000025.GIF" wi="990" he="278" /></maths>d(i)表示任意一组混凝土配比成分实验数据x(i)与新的不同混凝土配比成分x<sup>*</sup>在笛卡尔空间坐标系内直线距离的平方,△<sub>k</sub>的表达式意味着在在笛卡尔空间坐标系内实验数据点x(i)距离点x<sup>*</sup>越远,即d越大,则各个点的拟合误差|y(i)‑y<sub>k</sub>(i)|对x<sup>*</sup>位置附近区域总误差△<sub>k</sub>的贡献越小,换言之,与点x<sup>*</sup>越近的实验数据点x(i),其拟合误差|y(i)‑y<sub>k</sub>(i)|对x<sup>*</sup>位置附近区域总误差△<sub>k</sub>的贡献越大;步骤5、利用在现场按照标准混凝土强度检测方法实测得到的混凝土强度值y(i)对预测得到的混凝土强度值y′(i)进行检验,以验证其预测精度。
地址 362000 福建省泉州市丰泽区城东华侨大学