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.

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.

Fail-Safe Topology Optimization of Continuum Structures with Multiple Constraints Based on ICM Method

Traditional topology optimization methods may lead to a great reduction in the redundancy of the optimized structure due to unexpected material removal at the critical components.The local failure in critical components can instantly cause the overall failure in the structure.More and more scholars have taken the fail-safe design into consideration when conducting topology optimization.A lot of good designs have been obtained in their research,though limited regarding minimizing structural compliance(maximizing stiffness)with given amount of material.In terms of practical engineering applications considering fail-safe design,it is more meaningful to seek for the lightweight structure with enough stiffness to resist various component failures and/or to meet multiple design requirements,than the stiffest structure only.Thus,this paper presents a fail-safe topology optimization model for minimizing structural weight with respect to stress and displacement constraints.The optimization problem is solved by utilizing the independent continuous mapping(ICM)method combined with the dual sequence quadratic programming(DSQP)algorithm.Special treatments are applied to the constraints,including converting local stress constraints into a global structural strain energy constraint and expressing the displacement constraint explicitly with approximations.All of the constraints are nondimensionalized to avoid numerical instability caused by great differences in constraint magnitudes.The optimized results exhibit more complex topological configurations and higher redundancy to resist local failures than the traditional optimization designs.This paper also shows how to find the worst failure region,which can be a good reference for designers in engineering.

On the symplectic superposition method for free vibration of rectangular thin plates with mixed boundary constraints on an edge

A novel symplectic superposition method has been proposed and developed for plate and shell problems in recent years.The method has yielded many new analytic solutions due to its rigorousness.In this study,the first endeavor is made to further developed the symplectic superposition method for the free vibration of rectangular thin plates with mixed boundary constraints on an edge.The Hamiltonian system-based governing equation is first introduced such that the mathematical techniques in the symplectic space are applied.The solution procedure incorporates separation of variables,symplectic eigen solution and superposition.The analytic solution of an original problem is finally obtained by a set of equations via the equivalence to the superposition of some elaborated subproblems.The natural frequency and mode shape results for representative plates with both clamped and simply supported boundary constraints imposed on the same edge are reported for benchmark use.The present method can be extended to more challenging problems that cannot be solved by conventional analytic methods.

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.

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

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

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.

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.

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.

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.

深部软岩巷道围岩与锚喷U型钢支护结构相互作用研究

收敛-约束法(特性曲线法)是将理论解析、现场实测、工程经验相结合的一种地下工程结构设计方法,是目前分析围岩-支护相互作用关系及开展支护优化设计的常用方法。基于收敛-约束法的基本原理,总结给出了经典的支护特征与支护结构变形方程,理论计算获得了喷射混凝土、锚杆(索)、U型钢支架等支护结构的支护特征曲线,分析了支护结构的几何尺寸(直径、长度)、间排距、材料强度等参数对支护压力的影响特征;随着喷射混凝土厚度及强度等级的增加,喷射混凝土提供的支护刚度和支护压力逐渐增大;随着锚杆(索)直径、长度、杆体材料强度的增加及间排距的减小,锚杆的支护压力显著增加;U型钢支架的排距越小及材料强度越大,其提供的支护压力越大。采用FLAC^(3D)内嵌的莫尔库伦应变软化本构模型,建立了考虑岩石峰后应变软化与扩容特性的深部软岩巷道数值分析模型,计算获得了不同应力状态下巷道纵剖面变形曲线、围岩特征曲线,分析了锚喷、锚杆(索)喷、锚喷U型钢等3种联合支护技术对深部巷道围岩大变形控制的适用性,验证了锚喷U型钢联合支护技术应用于深部巷道支护工程的可行性。考虑岩石峰后应变软化与扩容特性的应变软化本构模型和经典的莫尔库伦本构模型的数值模拟结果相差较大,采用莫尔库伦本构模型的数值模拟结果保守,支护结构提供的支护压力无法满足深部巷道稳定性控制要求,支护后巷道围岩变形较大甚至会发生冒顶、片帮等安全事故。

基于调肝疏情法治疗肝气郁结型耳鸣

探讨基于调肝疏情法,以电针晕听区联合身心干预治疗肝气郁结型耳鸣的临床应用。肝气郁结型耳鸣源于情志不遂,肝气郁结,气机阻滞,阴阳升降失调而致耳鸣,治疗上应调肝疏情。电针晕听区可疏肝清窍、调和气血,身心干预治疗可调畅情志、平衡阴阳,二者联合共克耳鸣之难症。临床运用过程中,因其安全便捷,亦可根据患者不同伴随症状进行灵活加减,收效颇优。

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

为探明深部高地应力软岩隧道大变形产生机理,制定适应大变形控制措施,以月直山隧道为工程依托,首先明确求解围岩形变压力与松动压力的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,围岩大变形得到控制,研究成果可为今后相关隧道工程设计与施工提供参考。

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

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

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

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

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

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

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

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

对称约束的类间方差阈值方法

为了提高现有最大类间方差法(OTSU)的阈值化精度和适应性,提出了一种对称约束的类间方差阈值方法。该方法首先对输入图像使用Prewitt算子构建梯度幅值图像,并根据对称性原则提取对称采样区;然后,基于构建的对称约束类间方差目标函数最大化准则选取阈值,并判断在此阈值下对称采样区是否满足对称条件;当无法满足对称条件时,基于对称采样区对输入图像进行对称修正处理,并应用对称约束的类间方差目标函数对修正后的对称采样区选取阈值;最后,使用最终选取的阈值对输入图像阈值化。在28幅合成图像和70幅真实世界图像集上比较了提出的方法与OTSU法及4种OTSU的改进方法的阈值化性能。实验结果表明,提出方法的误分类率在合成图像和真实世界图像上分别为0.0106和0.016,相较于阈值化精度第2的方法在误分类方面分别降低了91.4%和86.1%。提出的方法虽然在计算效率方面不占有优势,但它对不同模态的测试图像具有更稳健的阈值化适应性和更高的阈值化精度。