基于malab语言的高性能混凝土配合比的优化设计上.docxVIP

基于malab语言的高性能混凝土配合比的优化设计上 江景平同志在《基于nb语言的高性能混凝土配合比优化设计》（以下简称“迄今为止”）一文中的关注，作者欢迎了江景平同志讨论“上”（以下简称“江景平同志介绍”）。“前文”采用Matlab编程实现高性能混凝土配合比的优化。Matlab是美国MathWorks公司研制开发的软件产品。它是当今世界上使用最为广泛的科学计算软件之一。Matlab具有强大、丰富的内置函数以及高度灵活的可编程性, 所以特别适合数据处理以及结果的图形化显示。另外, Matlab包含多种根据应用领域划分的工具箱, 这些工具箱极大地扩展了Matlab在这些领域的实用功能, 其中应用最为广泛的工具箱包括:数字信息处理、图像处理等工具箱等。现在, 无论是学术界还是在工业界, Matlab逐渐成为大多数科研人员的首选计算语言。而Matlab本身也逐渐成为一种工业界通用的计算语言。一些工业产品的研究与开发已开始完全使用Matlab, 比如日本丰田汽车公司的引擎控制单元 (ECU) 系统的研发等。在MathWorks公司的网站可以看到更多关于Matlab的工业应用的实例。由于篇幅所限, 现仅就讨论中涉及的几个问题答复如下。 1 关于临界水胶比的确定 众所周知, 新标准是在2001年以后实施的。“前文”没有注意文献发表的时间, 以新的水泥强度值代入旧的规范公式中导致结果出错。在此向江京平同志表示致谢。但在此也想对江京平同志对于水胶比和胶水比的概念提出如下意见。“原文”中对于42.5等级的水泥进行强度约束条件计算:采用“前文”的方法算出的C50、C60等级的混凝土临界胶水比为0.277、0.233, 经验算临界胶水比应该为3.610、4.295, 而0.277、0.233应该是C50、C60等级混凝土的临界水胶比;采用江京平同志的方法中算出的C50、C60等级的混凝土临界胶水比为0.356、0.298, 经验算临界胶水比应该是2.804、3.359, 而0.356、0.298应该是C50、C60等级的混凝土临界水胶比。 2 约束条件 经检查, Matlab源代码矩阵A第4行[0, 0, 0, 0.6, 0.4, 0]应该修改为[0, 0, 0, 0.6, -0.4, 0], 约束条件170≤x2≤210应该为130≤x2≤210, 这两处均为笔误, 向江京平同志表示致谢, 并向广大读者致歉。 3 多目标函数与fmincon函数的应用区别 作者在“前文”中采用序列二次规划法求解约束非线性问题, 它是目前求解约束非线性问题最有效的方法。该方法在计算机快捷运算速度的支持下, 用一系列的二次规划问题, 已逐渐在工程中获得应用, 因此选用Fgoalattain函数。虽然Fgoalattain是处理多目标的函数, 也没有明确说明不可以优化单目标问题;而fmincon函数是解决多变量约束问题的一种优化方法, 也是可以应用的。两种方法在计算结果上差别不大。用多目标函数处理单目标问题, 其实质是转化为单目标函数求解单目标问题。我们在“前文”中之所以采用多目标函数, 也是因为通过它, 我们可以比较便捷地对混凝土配合比优化问题从单目标拓展到多目标。 4 化学混凝土配合比优化方法的发展 随着社会的发展, 科学研究都需在不断革新, 这是必然规律。“前文”发表于2005年, 成稿于2004年3月, 限于当时采用数学规划模型解决混凝土配合比成本的优化问题还是基本处于初探阶段, 国内外都有研究者采用不同的数学规划模型对该问题求解, 作者也在这方面做了一些工作。对于江京平同志提出的优化方法表示赞同, 能让高性能混凝土配合比设计的优化方法发展得更加完美, 是我们共同的期待。再次感谢江京平同志的讨论, 并期望今后有机会与江京平同志再作交流。