Application of diflerential constraint method to exact solution of second-grade fluid

A differential constraint method is used to obtain analytical solutions of a second-grade fluid flow. By using the first-order differential constraint condition, exact solutions of Poiseuille flows, jet flows and Couette flows subjected to suction or blowing forces, and planar elongational flows are derived. In addition, two new classes of exact solutions for a second-grade fluid flow are found. The obtained exact solutions show that the non-Newtonian second-grade flow behavior depends not only on the material viscosity but also on the material elasticity. Finally, some boundary value problems are discussed.

Diagnosis and Resolution of Infeasibility in the Constraint Method for Solving Multi Objective Linear Programming Problems

In this paper we discuss about infeasibility diagnosis and infeasibility resolution, when the constraint method is used for solving multi objective linear programming problems. We propose an algorithm for resolution of infeasibility, which is a combination of interactive, weighting and constraint methods.Numerical examples are provided to illustrate the techniques developed.

Aerodynamic optimization of rotor airfoil based on multi-layer hierarchical constraint method

Rotor airfoil design is investigated in this paper. There are many difficulties for this highdimensional multi-objective problem when traditional multi-objective optimization methods are used. Therefore, a multi-layer hierarchical constraint method is proposed by coupling principal component analysis(PCA) dimensionality reduction and e-constraint method to translate the original high-dimensional problem into a bi-objective problem. This paper selects the main design objectives by conducting PCA to the preliminary solution of original problem with consideration of the priority of design objectives. According to the e-constraint method, the design model is established by treating the two top-ranking design goals as objective and others as variable constraints. A series of bi-objective Pareto curves will be obtained by changing the variable constraints, and the favorable solution can be obtained by analyzing Pareto curve spectrum. This method is applied to the rotor airfoil design and makes great improvement in aerodynamic performance. It is shown that the method is convenient and efficient, beyond which, it facilitates decision-making of the highdimensional multi-objective engineering problem.

Precipitation Retrieval from Himawari-8 Satellite Infrared Data Based on Dictionary Learning Method and Regular Term Constraint

In this paper,the application of an algorithm for precipitation retrieval based on Himawari-8 (H8) satellite infrared data is studied.Based on GPM precipitation data and H8 Infrared spectrum channel brightness temperature data,corresponding "precipitation field dictionary" and "channel brightness temperature dictionary" are formed.The retrieval of precipitation field based on brightness temperature data is studied through the classification rule of k-nearest neighbor domain (KNN) and regularization constraint.Firstly,the corresponding "dictionary" is constructed according to the training sample database of the matched GPM precipitation data and H8 brightness temperature data.Secondly,according to the fact that precipitation characteristics in small organizations in different storm environments are often repeated,KNN is used to identify the spectral brightness temperature signal of "precipitation" and "non-precipitation" based on "the dictionary".Finally,the precipitation field retrieval is carried out in the precipitation signal "subspace" based on the regular term constraint method.In the process of retrieval,the contribution rate of brightness temperature retrieval of different channels was determined by Bayesian model averaging (BMA) model.The preliminary experimental results based on the "quantitative" evaluation indexes show that the precipitation of H8 retrieval has a good correlation with the GPM truth value,with a small error and similar structure.

Partition method for impact dynamics of flexible multibody systems based on contact constraint

The impact dynamics of a flexible multibody system is investigated. By using a partition method, the system is divided into two parts, the local impact region and the region away from the impact. The two parts are connected by specific boundary conditions, and the system after partition is equivalent to the original system. According to the rigid-flexible coupling dynamic theory of multibody system, system＇s rigid-flexible coupling dynamic equations without impact are derived. A local impulse method for establishing the initial impact conditions is proposed. It satisfies the compatibility con- ditions for contact constraints and the actual physical situation of the impact process of flexible bodies. Based on the contact constraint method, system＇s impact dynamic equa- tions are derived in a differential-algebraic form. The contact/separation criterion and the algorithm are given. An impact dynamic simulation is given. The results show that system＇s dynamic behaviors including the energy, the deformations, the displacements, and the impact force during the impact process change dramatically. The impact makes great effects on the global dynamics of the system during and after impact.

Power System State Estimation Solution With Zero Injection Constraints Using Modified Newton Method and Fast Decoupled Method in Polar Coordinate

如何保证零注入节点的注入功率在状态估计结果中严格为0是电力系统状态估计研究中的重要问题。在直角坐标下,由于零注入约束为线性约束,可使用修正牛顿法来有效地解决这一问题。因此,借鉴直角坐标下修正牛顿法的思路,提出了极坐标下的修正牛顿法和修正快速解耦估计。这些方法的计算流程与传统的极坐标下的牛顿法和快速解耦估计非常相似,计算速度与大权重法相当,同时能够保证零注入约束严格满足。仿真结果验证了所得结论。

A LEVEL SET METHOD FOR STRUCTURAL TOPOLOGY OPTIMIZATION WITH MULTI-CONSTRAINTS AND MULTI-MATERIALS

Combining the vector level set model,the shape sensitivity analysis theory with the gradient projection technique,a level set method for topology optimization with multi-constraints and multi-materials is presented in this paper.The method implicitly describes structural material in- terfaces by the vector level set and achieves the optimal shape and topology through the continuous evolution of the material interfaces in the structure.In order to increase computational efficiency for a fast convergence,an appropriate nonlinear speed mapping is established in the tangential space of the active constraints.Meanwhile,in order to overcome the numerical instability of general topology opti- mization problems,the regularization with the mean curvature flow is utilized to maintain the interface smoothness during the optimization process.The numerical examples demonstrate that the approach possesses a good flexibility in handling topological changes and gives an interface representation in a high fidelity,compared with other methods based on explicit boundary variations in the literature.

Constraint Stabilization and One-Step Correction Method for Dynamical Systems

A computational method of constraint stabilization and correction is introduced. The method is based on the Baumgart＇s one-step method. Constraint conditions are addressed to stabilize and correct the solution. Two examples are given to illustrate the results of the method.

考虑湿度应力的深部软岩隧道大变形控制研究

为探明深部高地应力软岩隧道大变形产生机理,制定适应大变形控制措施,以月直山隧道为工程依托,首先明确求解围岩形变压力与松动压力的Kastner与Caquot公式,基于岩体弹塑性力学与连续介质理论建立围岩湿度应力解析解,采用收敛-约束法绘制出3种围岩应力作用下围岩与初期支护特征曲线,分析断面变形过程中围岩与支护结构相互作用规律及3种应力占比演化规律。分别以混凝土与型钢钢架作为二衬支护结构对月直山隧道围岩稳定性进行计算分析,明确考虑湿度应力与松动压力条件下隧道二衬最优支护时机与支护参数,以支护结构安全系数FS评判出最优支护方案并对隧道软岩大变形进行治理。结果表明:当围岩径向位移达到550 mm时,仅考虑形变压力Pi时围岩对支护结构的压力P仅为0.813 MPa,考虑湿度应力Pw与松动压力Pa时P增大为1.372 MPa,湿度应力与松动应力总占比达40.7%,仅考虑形变压力而设计的支护结构不满足围岩稳定性要求;根据“卸压支护”理念,确定以“位移释放峰值upeak=0.325 m”与“间距d=0.7 m”作为增设第二层钢架的最佳支护时机与支护参数,现场二次钢架设置24天后使围岩变形收敛于50.8 mm,围岩大变形得到控制,研究成果可为今后相关隧道工程设计与施工提供参考。

Novel Method to Handle Inequality Constraints for Nonlinear Programming

By redefining the multiplier associated with inequality constraint as a positive definite function of the originally-defined multiplier, say, u2_i, i=1, 2, ..., m, nonnegative constraints imposed on inequality constraints in Karush-Kuhn-Tucker necessary conditions are removed. For constructing the Lagrange neural network and Lagrange multiplier method, it is no longer necessary to convert inequality constraints into equality constraints by slack variables in order to reuse those results dedicated to equality constraints, and they can be similarly proved with minor modification. Utilizing this technique, a new type of Lagrange neural network and a new type of Lagrange multiplier method are devised, which both handle inequality constraints directly. Also, their stability and convergence are analyzed rigorously.

考虑个体客户满意度的多车程生鲜品配送路径优化

随着生鲜电商的兴起,生鲜品的道路运输量越来越大,其配送任务越来越复杂,提升客户满意度是生鲜品配送的重要目标组成。文章以一个具有配送自主权的生鲜品供应商为研究对象,以整个计划周期内生鲜品配送的车辆运输成本最小及客户满意度最大为双目标,且客户满意度从顾客的个体视角考虑,构建双目标多车程带时间窗及考虑订单打包时间的车辆路径优化模型。基于ε-约束方法,设计了两阶段变邻域搜索(VNS)和模拟退火(SA)相结合的元启发式算法,基于车辆构造初始路线,并针对性地提出了3点改进措施,包括通过数学性质推导下界、减少求解过程中不必要的搜索空间以及规避求解过程中被支配解的后优化处理。最后,以Solomon标准算例为基础构造算例对所设计的算法进行了计算实验,其求解结果表明:其设计的算法具有有效性,并通过敏感性分析得出相关管理启示。

薄壁杆件屈曲模态分析的力特征约束有限条法

整体、畸变、局部等基本变形模式的划分是条理化开展薄壁杆件工程计算的必要手段。薄壁杆件的基本变形模式传统上基于简化应力-应变关系,采用变形、应变等几何特征进行定义,不利于因应复杂化应力-应变关系。提出一套完全基于力特征和正交完备性原则的基本模式定义。与广义梁理论和约束有限条法相比,提出的整体、畸变、局部3种基本变形模式关于刚度严格正交,且完整覆盖薄壁杆件全变形域。基于该定义实现了针对薄壁杆件有限条模型的屈曲模态分解和识别。算例表明,提出的基本模式定义和屈曲模态分析方法具有开、闭口和折、曲线形截面的一致适应能力,剪应变和横向伸缩的影响在基本屈曲模式中得到合理表达,曲线形截面薄壁构件具有与多边形截面薄壁构件一致的整体、畸变和局部屈曲机理。

Levenberg-Marquardt Method for Mathematical Programs with Linearly Complementarity Constraints

In this paper, a new method for solving a mathematical programming problem with linearly complementarity constraints (MPLCC) is introduced, which applies the Levenberg-Marquardt (L-M) method to solve the B-stationary condition of original problem. Under the MPEC-LICQ, the proposed method is proved convergent to B-stationary point of MPLCC.

深部地层智能压井多解性分析与优化策略

开发深部地层油气资源时普遍存在地质条件复杂、钻井周期长和井筒压力控制困难等问题,采用智能压井方法结合多源实时信息反馈,可实现井筒内气液分布状态和压力变化规律的实时预测与更新,但不同修正系数组合可能得到相同的压力计算结果,导致模型存在多解性难题。为此,分析了不同历史时间节点解空间形态的演变规律,揭示了模型多解性的本质源于少量数据约束下模型训练的不完善性;并对应建立了基于实时信息序列的模型全局训练优化方法及动态随机种群训练优化方法,测试了其对于模型全局最优解的搜索能力及适用条件。测试结果表明,全局训练优化方法在压井初期能够实现精准调控,但计算耗时较长;而动态随机种群训练优化方法在压井初期与预期值略有差异,但计算耗时较少。根据可用计算资源情况选择合适的训练优化方法,可实现多源实时数据约束下模型关于井筒气液流动规律的深度学习。

金融能力、正规信贷约束与农户创业行为

基于金融能力和正规信贷约束视角,利用2017年和2019年CHFS微观数据,通过构建Probit模型实证分析金融能力对农户创业行为的影响,探讨正规信贷约束在金融能力和农户创业行为关系中发挥的中介效应。结果表明:金融能力对农户创业行为具有显著正向影响,且在考虑内生性的情况下,该结论仍然成立;中介机制检验表明,金融能力能够通过缓解正规信贷约束来促进农户创业行为;异质性分析表明,中、西部地区和三线及以下城市的农户金融能力对其创业行为的影响更为显著。据此,建议建立完善的金融能力培养体系,加大对欠发达地区的政策扶持,且提高正规信贷支持力度。

数字经济对实体企业金融化的影响研究——以云南省为例

数字经济是推动我国经济高质量发展的新动能,从实体企业金融化的角度出发,基于云南省2011——2019年沪深A股上市公司面板数据,运用主成分分析法构建城市数字经济水平综合指标,利用固定效应模型和中介效应模型,实证考察云南省地级市数字经济水平金融化效应及作用机制。研究发现数字经济发展对实体企业金融化具有正向影响。结果表明数字经济通过融资约束促进企业金融化,异质性分析发现数字经济发展对高新技术企业过度金融化表现出治理效应,对非高新技术企业金融化体现促进作用。

基于重加权L1的ATpV正则化叠前反演方法

地震叠前反演能够准确获取地下储层介质的各类参数,是油气的勘探与开发中重要技术之一。然而,地震反演是典型的病态问题,为了克服此问题,通常使用正则化约束目标函数,来减轻反演问题的病态性。但是正则化约束忽略了地层边界的振幅信息,使用重加权方法可以很好地克服这一问题,更好地恢复稀疏性。提出了一种基于重加权L1的ATpV正则化叠前三参数反演方法(ATpV-L1方法),首次将重加权L1方法与ATpV方法结合,并引入到叠前反演中。采用交替方向乘子算法(ADMM)建立反演框架,对目标函数进行分块优化,有效提高了收敛速度。首先,介绍ATpV-L1方法,建立了基于ATpV-L1的叠前反演目标函数;然后,应用理论模拟数据对比新方法和ATpV方法反演结果,验证了方法的效果;最后,使用实际数据进行实验分析,进一步验证了ATpV-L1方法的反演精度及可行性。实验结果表明,提出的ATpV-L1方法可以有效恢复反演结果的稀疏性,提高反演精度。

一种状态约束下的抗差自适应联邦滤波算法

充分利用先验约束信息可以提高多传感器组合导航的定位精度和可靠性。将状态约束下的卡尔曼滤波扩展到传统联邦滤波中,提出一种状态约束下的联邦滤波算法。当子传感器出现异常时,在状态约束下的联邦滤波基础上,采用Huber方法调整子滤波器观测噪声矩阵,同时在信息分配阶段引入自适应信息分配因子,实时调整子滤波器融合权重,得到一种状态约束下的抗差自适应联邦滤波算法,以进一步减少不准确的子滤波器估计对融合结果的影响。将该方法应用在捷联惯导、全球导航卫星系统和里程计的多传感器组合导航系统中。仿真实验表明,状态约束下的联邦滤波估计精度优于传统联邦滤波,状态约束下的抗差自适应联邦滤波能够进一步提高观测异常下的导航定位精度和可靠性。

考虑条件风险价值的交直流系统两阶段分布鲁棒低碳经济优化调度

考虑到海上风电出力的随机性以及日益突出的生态环境问题,以含柔性直流输电技术(voltagesource converter high voltage direct current,VSC-HVDC)的交直流系统为研究对象,提出了考虑条件风险价值(conditional valueatrisk,CVaR)的两阶段分布鲁棒低碳经济优化模型,构建了基于Kullback-Leibler(KL)散度的概率分布模糊集,同时利用条件风险价值量化了极端场景下的尾部风险,使得模型能够同时考虑概率分布不确定性以及处于最坏概率分布中极端场景下的尾部损失;此外,将阶梯型碳交易机制并入所提分布鲁棒模型中,通过合理利用柔性资源和储能装置,增强系统运行的灵活性,在兼顾运行风险的前提下,降低碳排放量的目标。再者,为了提高计算效率,在列和约束生成算法(column-and-constraint generation method,C&CG)和Multi-cut Benders分解算法的基础上提出了双循环分解算法。最后,在基于改进的IEEE RTS 79测试系统中验证了所提模型及算法的有效性。