Analysis of elastic-viscoplastic creep model based on variable-order differential operator

Based on the variable-order differential operator (VODO), a variable-order elastic-viscoplastic (VO-EVP) creep model of geomaterials is established. Firstly, the variable-order strain rate model is established according to the newly defined Riesz-Caputo (RC) VODO, and the creep equation of the VO-EVP model is solved by Laplace transform. Then, the order function expression suitable for the creep of geomaterials is proposed by analyzing the existing order functions. A series of numerical analyses are carried out on the VO-EVP model and the order function. Finally, combined with the results of triaxial creep tests of Hong Kong marine deposits and Sackville clay, the effectiveness of the model under static load test and the sensitivity of the order function are analyzed. The results show that the VO-EVP model is significantly better than the conventional EVP model. Compared with the constant fractional EVP (CF-EVP) model, the VO-EVP model considering only the internal variables in time does not have obvious advantages in predicting the creep behavior of soil under static load, nor does it show strong time memory and parameter sensitivity. However, the VO-EVP model also has important significance in studying creep behavior.     

饱和粘土非线性渗流规律与径向固结

以多孔介质等效渗流概念得出粘土平均孔隙微尺度数量级范围在0.01 μm至0.1μm,与测试结果一致．实验结果表明饱和粘土微尺度孔隙渗流为非线性流．理论推导表明固液界面作用与渗透率平方根或孔隙半径成反比,固液界面相互作用是导致饱和粘土非线性渗流的重要原因．提出了精确描述饱和粘土微尺度孔隙非线性渗流基本规律的数学模型,其参数量纲明确,物理意义清楚．建立了饱和粘土非线性渗流新定律,统一描述了从较低到较高水力梯度全过程渗流曲线特征,达西定律是其特例．基于新定律,建立了饱和粘土非线性渗流定流量径向固结数学模型．以粘性边界层思想与稳态依序替换法及积分方程法,导出了粘土非线性渗流平均质量守恒方程及活动边界运动方程,给出了饱和粘土非线性渗流超孔隙水压力分布公式与平均固结度计算公式,获得了粘土层压力分布规律和平均固结度随时间变化规律．结果表明:饱和粘土非线性渗流使活动边界运动速度减小．研究结果为粘土地质工程与岩土工程应用提供了新的科学依据．达西渗流径向固结计算是新的非线性渗流固结计算的特例．     

A novel method for non-probabilistic convex modelling based on data from practical engineering

In this paper, a novel method for non-probabilistic convex modelling with the bounds to precisely encircle all the data of uncertain parameters extracted from practical engineering is developed. The method is based on the traditional statistical method and the correlation analysis technique. Mean values and correlation coefficients of uncertain parameters are first calculated by utilizing the information of all the given data. Then, a simple yet effective optimization procedure is first introduced in the mathematical modelling process for uncertain parameters to obtain their precise bounds. This procedure works by optimizing the area of the convex model, at the same time, covering all the given data. Thus, the effective mathematical expression of the convex models are finally formulated. To test the prediction capability and generalization ability of the proposed convex modelling method, evaluation criteria, i.e. volume ratio, standard volume ratio, and prediction accuracy are established. The performance of the proposed method is systematically studied and compared with other existing competitive methods through test standards. The results demonstrate the effectiveness and efficiency of the present method.     

考虑订单合并和货物转运的多式联运路径优化研究

为了获得运输的规模经济效应,本文研究了一种考虑订单合并和货物转运的零担多式联运路径优化问题。首先,以总运输成本为目标函数,以网络中的运输工具容量、可以提供的运输工具最大数量、运输工具服务的关闭时间以及订单时间窗为约束,构建混合整数规划模型,在模型中允许多个订单进行合并运输并考虑运输过程中的转运成本。其次,由于多式联运路径优化问题是典型的NP-hard问题,为了快速求解该模型,开发了一种可以快速为该问题提供近似最优解和下界的列生成启发式算法。最后,生成并测试了大量算例,结果表明所开发的列生成启发式算法可以在较短的时间内提供高质量的近似最优解。文章所构建的模型和开发的列生成启发式算法可以为零担自营多式联运物流企业提供高效的决策支持。     

低维混沌时序非线性动力系统的预测方法及其应用研究

主要研究由低维混沌时序所确定的非线性动力系统的预测方法及其应用。在国外学者研究工作的基础上,应用一种非线性混沌模型在相空间内对时序进行重构工作,先通过改进的最小二乘方法来估计模型的参数,满足一定精度后,再采用最优化方法来估计模型的参数,并用所求得的混沌时序模型在其相空间内对时序的未来值进行预测。给出了非常有代表性的实例对文中模型和算法进行验证。结果发现采用该算法能较准确地求得模型的参数,在相空间中对混沌时序进行预测,将传统方法中的外推变成了相空间中的内插,及选取最佳的模型阶数等工作都能增加预测的准确程度,且混沌时序不可能进行长期的预测。     

Mathematical modelling and study of the consolidation of an elastic saturated soil with an incompressible fluid by FEM

The main aim of this paper is to validate and to solve a model for consolidation of an elastic saturated soil with incompressible fluid. Firstly, we prove the existence and uniqueness of the solution of the variational problem corresponding to an initial and boundary value problem (IBVP): a special case of the Biot’s ‘consolidation of clay’ model (where the applied forces depend on time). Secondly, we prove the stability of the method as well as the estimation of the error by using semi-discretization in time. Finally, we then solved this one by the finite element method (FEM) employing repeated fixed point techniques in order to obtain the results for displacement and pore water pressure. The pore fluid is considered incompressible. The results of the numerical experiments are compared with analytical solutions and, in cases where such solutions do not exist, with experimental data.     

An Improved Spectral Conjugate Gradient Algorithm for Nonconvex Unconstrained Optimization Problems

In this paper, an improved spectral conjugate gradient algorithm is developed for solving nonconvex unconstrained optimization problems. Different from the existent methods, the spectral and conjugate parameters are chosen such that the obtained search direction is always sufficiently descent as well as being close to the quasi-Newton direction. With these suitable choices, the additional assumption in the method proposed by Andrei on the boundedness of the spectral parameter is removed. Under some mild conditions, global convergence is established. Numerical experiments are employed to demonstrate the efficiency of the algorithm for solving large-scale benchmark test problems, particularly in comparison with the existent state-of-the-art algorithms available in the literature.     

Empirical mode decomposition-based time grid refinement optimization approach for optimal control problems

This paper presents a novel and efficient grid refinement approach for choosing discretization level of control vector parameterization (CVP) method. Different from the uniform time grid CVP method, the proposed method adaptively refines the time grid knots by applying the empirical mode decomposition at each iteration. An attractive property of the proposed method is that fewer parameters are needed to achieve better optimization results, which means that a high-quality solution can be obtained with lower computation cost. The related derivation shows the advantage of simplicity for implementation. Compared with the existing approaches, for which all the control variables are parameterized on the same time grid, the proposed method is more flexible. A standard optimal control problem is tested as an illustration to demonstrate the validity of the proposed method, results indicate better performance index and lower computation cost of the proposed method.     

聚氯乙烯接近玻璃化转变温度时的蠕变行为研究

利用动态力学分析仪,详细研究了小应力作用下聚氯乙烯（PVC）接近玻璃化转变温度（63℃）时的短期蠕变性能．根据Struik-Kohlrausch (SK)方程及Struik移位方法,分析了物理老化对PVC瞬时蠕变柔量的影响,并提出了基于多松弛机制的新模型定量描述初始蠕变柔量、松弛时间和形状因子．结果表明:SK方程虽可表征不同老化时间下PVC的蠕变性能,但通过Struik移位方法却无法获得理想的主曲线,且初始蠕变柔量、松弛时间、形状因子和蠕变率均受物理老化显著影响,并在双对数坐标下随老化时间呈明显的非线性变化;新模型不仅揭示了物理老化控制参量的演化规律,而且极好地预测了PVC在物理老化进程中的高温蠕变行为．     

Creep of Heat-Resistant Composites of an Oxide-Fiber/Ni-Matrix Family

A creep model of a composite with a creeping matrix and initially continuous elastic brittle fibers is developed. The model accounts for the fiber fragmentation in the stage of unsteady creep of the composite, which ends with a steady-state creep, where a minimum possible average length of the fiber is achieved. The model makes it possible to analyze the creep rate of the composite in relation to such parameters of its structure as the statistic characteristics of the fiber strength, the creep characteristics of the matrix, and the strength of the fiber-matrix interface, the latter being of fundamental importance. A comparison between the calculation results and the experimental ones obtained on composites with a Ni-matrix and monocrystalline and eutectic oxide fibers as well as on sapphire fiber/TiAl-matrix composites shows that the model is applicable to the computer simulation of the creep behavior of heat-resistant composites and to the optimization of the structure of such composites. By combining the experimental data with calculation results, it is possible to evaluate the heat resistance of composites and the potential of oxide-fiber/Ni-matrix composites. The composite specimens obtained and tested to date reveal their high creep resistance up to a temperature of 1150°C. The maximum operating temperature of the composites can be considerably raised by strengthening the fiber-matrix interface.     

基于交互式模糊规划方法的可持续闭环供应链网络规划

在市场需求、设施开设成本和产品回收率不确定的条件下,采用一种交互式可能性规划方法,研究由多个工厂、分销点、市场和废旧点构成的可持续闭环供应链网络设计问题。基于可持续闭环供应链网络结构,构建以企业运营成本和环境伤害最小、社会效益最大为目标的混合整数规划模型。同时,引入改进Epsilon约束方法将多目标优化问题转化为单目标优化问题,在此基础上提出一种两阶段可能性规划方法,基于TH模糊方法对不确定性参数进行处理。最后,通过数值实例,验证本文所建可持续闭环供应链网络模型的有效性,并对悲观-乐观值、不确定参数最低可接受水平β、可调参数γ进行敏感性分析;通过与其他模糊方法对比表明,采用TH模糊方法能得到稳定的最优解。     

Distributionally robust parameter identification of a time-delay dynamical system with stochastic measurements

In this paper, we consider a parameter identification problem involving a time-delay dynamical system, in which the measured data are stochastic variable. However, the probability distribution of this stochastic variable is not available and the only information we have is its first moment. This problem is formulated as a distributionally robust parameter identification problem governed by a time-delay dynamical system. Using duality theory of linear optimization in a probability space, the distributionally robust parameter identification problem, which is a bi-level optimization problem, is transformed into a single-level optimization problem with a semi-infinite constraint. By applying problem transformation and smoothing techniques, the semi-infinite constraint is approximated by a smooth constraint and the convergence of the smooth approximation method is established. Then, the gradients of the cost and constraint functions with respect to time-delay and parameters are derived. On this basis, a gradient-based optimization method for solving the transformed problem is developed. Finally, we present an example, arising in practical fermentation process, to illustrate the applicability of the proposed method.     

Modeling of nonlinear viscoelastic contact problems with large deformations

Contact problems are one of the most important engineering problems. These problems become much more tedious when one of the contacting bodies behaves nonlinear viscoelasticity and large deformations. This paper presents an incremental-iterative finite element model for the analysis of two dimensional quasistatic frictionless contact problems. Nonlinear viscoelastic behavior and large deformations are considered. The Schapery’s single-integral creep model with stress-dependent properties is used for nonlinear viscoelasticity. The constitutive equations are transformed into an incremental form resulting in a recursive relationship. Thereby, the need to store the entire strain histories is eliminated, except that from the previous time increment. The updated Lagrangian formulation is used to model the material and geometrical nonlinearities. Also, the Lagrange multiplier method is adopted to enforce the contact constraints. The converged solution is obtained using the Newton–Raphson iterative technique. The developed model has been verified with the previously published works and found a good agreement with them. To demonstrate the efficient capability of the developed computational model, three contact problems with different nature are analyzed.     

On the Dynamics of Bundled Conductors in Overhead Transmission Lines

The Energy Balance Principle (EBP) is well established for estimating the vibration levels of wind induced oscillations of single overhead transmission lines. The mathematical model, wherein a conductor is treated as a continuous system, results in a transcendental eigenvalue problem (EVP), which gives numerical difficulties in the case of bundled conductors. In this paper, different approaches for solving transcendental EVP and their relative merits are discussed. A new method named continuous spectrum approach provides a good engineering solution. Results from different approaches are compared. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Application of robust optimization to automated test assembly

In automated test assembly (ATA), 0-1 linear programming (0-1 LP) methods are applied to select questions (items) from an item bank to assemble an optimal test. The objective in this 0-1 LP optimization problem is to assemble a test that measures, in as precise a way as possible, the ability of candidates. Item response theory (IRT) is commonly applied to model the relationship between the responses of candidates and their ability level. Parameters that describe the characteristics of each item, such as difficulty level and the extent to which an item differentiates between more and less able test takers (discrimination) are estimated in the application of the IRT model. Unfortunately, since all parameters in IRT models have to be estimated, they do have a level of uncertainty to them. Some of the other parameters in the test assembly model, such as average response times, have been estimated with uncertainty as well. General 0-1 LP methods do not take this uncertainty into account, and overestimate the predicted level of measurement precision. In this paper, alternative robust optimization methods are applied. It is demonstrated how the Bertsimas and Sim method can be applied to take this uncertainty into account in ATA. The impact of applying this method is illustrated in two numerical examples. Implications are discussed, and some directions for future research are presented.     

Experimental Evaluation of Heuristic Optimization Algorithms: A Tutorial

Heuristic optimization algorithms seek good feasible solutions to optimization problems in circumstances where the complexities of the problem or the limited time available for solution do not allow exact solution. Although worst case and probabilistic analysis of algorithms have produced insight on some classic models, most of the heuristics developed for large optimization problem must be evaluated empirically—by applying procedures to a collection of specific instances and comparing the observed solution quality and computational burden.This paper focuses on the methodological issues that must be confronted by researchers undertaking such experimental evaluations of heuristics, including experimental design, sources of test instances, measures of algorithmic performance, analysis of results, and presentation in papers and talks. The questions are difficult, and there are no clear right answers. We seek only to highlight the main issues, present alternative ways of addressing them under different circumstances, and caution about pitfalls to avoid.     

目标规划在建立水泥最佳配料方法中的应用

本文在用正交设计理论得出水泥最佳配料条件下，根据此最佳配料条件参数和原燃料成分的化验数据，应用线性规划理论与方法建立数学模型，利用计算机算出水泥最佳配料方案和技术上要求的一些重要参数。本课题的完成，改变了国内水泥配料沿用的经验公式方法，建立了一种崭新的水泥配料计算方法。     

镍基定向结晶合金蠕变损伤的细观模型

针对镍基定向结晶合金的蠕变变形和损伤,本文建立了有限元和自洽两种晶体滑移细观模型。在这两种模型中同时考虑了晶粒滑移损伤和晶界损伤。通过一组镍基单晶合金和一组特殊结晶取向的定向结晶合金标定了模型参数。数值分析结果表明,这两种模型都可以比较准确地描述定向结晶合金的蠕变损伤变形和失效规律,尤其是定向结晶合金蠕变性能的结晶取向相关性。     

饱和土体自重固结问题的相似解

采用了比以往更为普遍的土体物理力学性质假设,利用Hopf-Cole变换方法求得了厚层土体在自重应力作用下的非线性固结问题的完整解析解答．通过试验数据,将该解答与传统的大应变线性固结理论解答和基于实验数据的有限元数值解答相比较,结果表明,该解答能够更好地描述土体的实际固结过程,而由线性化固结理论所得的解答对固结过程中的沉降量和固结度的估计偏小．