This paper presents a new dimension reduction strategy for medium and large-scale linear programming problems. The proposed method uses a subset of the original constraints and combines two algorithms: the weighted av...This paper presents a new dimension reduction strategy for medium and large-scale linear programming problems. The proposed method uses a subset of the original constraints and combines two algorithms: the weighted average and the cosine simplex algorithm. The first approach identifies binding constraints by using the weighted average of each constraint, whereas the second algorithm is based on the cosine similarity between the vector of the objective function and the constraints. These two approaches are complementary, and when used together, they locate the essential subset of initial constraints required for solving medium and large-scale linear programming problems. After reducing the dimension of the linear programming problem using the subset of the essential constraints, the solution method can be chosen from any suitable method for linear programming. The proposed approach was applied to a set of well-known benchmarks as well as more than 2000 random medium and large-scale linear programming problems. The results are promising, indicating that the new approach contributes to the reduction of both the size of the problems and the total number of iterations required. A tree-based classification model also confirmed the need for combining the two approaches. A detailed numerical example, the general numerical results, and the statistical analysis for the decision tree procedure are presented.展开更多
Based on the dynamical theory of multi-body systems with nonholonomic constraints and an algorithm for complementarity problems, a numerical method for the multi-body systems with two-dimensional Coulomb dry friction ...Based on the dynamical theory of multi-body systems with nonholonomic constraints and an algorithm for complementarity problems, a numerical method for the multi-body systems with two-dimensional Coulomb dry friction and nonholonomic constraints is presented. In particular, a wheeled multi-body system is considered. Here, the state transition of stick-slip between wheel and ground is transformed into a nonlinear complementarity problem (NCP). An iterative algorithm for solving the NCP is then presented using an event-driven method. Dynamical equations of the multi-body system with holonomic and nonholonomic constraints are given using Routh equations and a con- straint stabilization method. Finally, an example is used to test the proposed numerical method. The results show some dynamical behaviors of the wheeled multi-body system and its constraint stabilization effects.展开更多
This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zer...This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.展开更多
In this paper, a branch-and-bound method for solving multi-dimensional quadratic 0-1 knapsack problems was studied. The method was based on the Lagrangian relaxation and the surrogate constraint technique for finding ...In this paper, a branch-and-bound method for solving multi-dimensional quadratic 0-1 knapsack problems was studied. The method was based on the Lagrangian relaxation and the surrogate constraint technique for finding feasible solutions. The Lagrangian relaxations were solved with the maximum-flow algorithm and the Lagrangian bounds was determined with the outer approximation method. Computational results show the efficiency of the proposed method for multi-dimensional quadratic 0-1 knapsack problems.展开更多
Through detailed three-dimensional (3D) finite element (FE) calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are o...Through detailed three-dimensional (3D) finite element (FE) calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are obtained near the crack front with aspect ratios (a/c) of 0.2, 0.4, 0.5, 0.6, 0.8 and 1.0. Tz decreases from an approximate value of Poisson ratio v at the crack tip to zero with increasing normalized radial distances (r/a) in the normal plane of the crack front line, and increases gradually when the elliptical parameter angle φ changes from 0° to 90°at the same r/a. With a/c rising to 1.0, Tz is getting nearly independent of φ and is only related to r/a. Based on the present FE calculations for Tz, empirical formulas for Tz are obtained to describe the 3D distribution of Tz for embedded center-elliptical cracks using the least squares method in the range of 0.2 ≤ a/c ≤ 1.0. These Tz results together with the corresponding stress intensity factor K are well suitable for the analysis of the 3D embedded centerelliptical crack from field, and a two-parameter K-Tz principle is proposed.展开更多
Finding a suitable space is one of the most critical problems for dimensionality reduction. Each space corresponds to a distance metric defined on the sample attributes, and thus finding a suitable space can be conver...Finding a suitable space is one of the most critical problems for dimensionality reduction. Each space corresponds to a distance metric defined on the sample attributes, and thus finding a suitable space can be converted to develop an effective distance metric. Most existing dimensionality reduction methods use a fixed pre-specified distance metric. However, this easy treatment has some limitations in practice due to the fact the pre-specified metric is not going to warranty that the closest samples are the truly similar ones. In this paper, we present an adaptive metric learning method for dimensionality reduction, called AML. The adaptive metric learning model is developed by maximizing the difference of the distances between the data pairs in cannot-links and those in must-links. Different from many existing papers that use the traditional Euclidean distance, we use the more generalized l<sub>2,p</sub>-norm distance to reduce sensitivity to noise and outliers, which incorporates additional flexibility and adaptability due to the selection of appropriate p-values for different data sets. Moreover, considering traditional metric learning methods usually project samples into a linear subspace, which is overstrict. We extend the basic linear method to a more powerful nonlinear kernel case so that well capturing complex nonlinear relationship between data. To solve our objective, we have derived an efficient iterative algorithm. Extensive experiments for dimensionality reduction are provided to demonstrate the superiority of our method over state-of-the-art approaches.展开更多
The theoretical solutions are obtained for the three-dimensional(3-D)stress field in an infinite isotropic elastic plate with a through-the-thickness circular hole subjected to shear load at far field by using Kane an...The theoretical solutions are obtained for the three-dimensional(3-D)stress field in an infinite isotropic elastic plate with a through-the-thickness circular hole subjected to shear load at far field by using Kane and Mindlin′s assumption based on the stress function method.Based on the present solutions,the characteristics of 3-D stress field are analyzed and the emphasis is placed on the effects of the plate thickness and Poisson′s ratio on the deviation of the present 3-D in-plane stress from the related plane stress solutions,the stress concentration and the out-of-plane constraint.The present solutions show that the stress concentration factor reaches its peak value of about 8.9% which is higher than that of the plane stress solutions.As expected,the out-of-plane stress constraint factor can reach 1on the surface of the hole when the plate is a very thick one.展开更多
An explicit Bargmann symmetry constraint is computed and its associated binary nonlinearization of Lax pairs is carried out for the super NLS-MKdV hierarchy. Under the obtained symmetry constraint, the n-th flow of th...An explicit Bargmann symmetry constraint is computed and its associated binary nonlinearization of Lax pairs is carried out for the super NLS-MKdV hierarchy. Under the obtained symmetry constraint, the n-th flow of the super NLS-MKdV hierarchy is decomposed into two super finite-dimensional integrable Hamiltonian systems, defined over the super-symmetry manifold R4N|2N with the corresponding dynamical variables x and tn. The integrals of motion required for Liouville integrability are explicitly given.展开更多
文摘This paper presents a new dimension reduction strategy for medium and large-scale linear programming problems. The proposed method uses a subset of the original constraints and combines two algorithms: the weighted average and the cosine simplex algorithm. The first approach identifies binding constraints by using the weighted average of each constraint, whereas the second algorithm is based on the cosine similarity between the vector of the objective function and the constraints. These two approaches are complementary, and when used together, they locate the essential subset of initial constraints required for solving medium and large-scale linear programming problems. After reducing the dimension of the linear programming problem using the subset of the essential constraints, the solution method can be chosen from any suitable method for linear programming. The proposed approach was applied to a set of well-known benchmarks as well as more than 2000 random medium and large-scale linear programming problems. The results are promising, indicating that the new approach contributes to the reduction of both the size of the problems and the total number of iterations required. A tree-based classification model also confirmed the need for combining the two approaches. A detailed numerical example, the general numerical results, and the statistical analysis for the decision tree procedure are presented.
基金Project supported by the National Natural Science Foundation of China(Nos.11372018 and 11572018)
文摘Based on the dynamical theory of multi-body systems with nonholonomic constraints and an algorithm for complementarity problems, a numerical method for the multi-body systems with two-dimensional Coulomb dry friction and nonholonomic constraints is presented. In particular, a wheeled multi-body system is considered. Here, the state transition of stick-slip between wheel and ground is transformed into a nonlinear complementarity problem (NCP). An iterative algorithm for solving the NCP is then presented using an event-driven method. Dynamical equations of the multi-body system with holonomic and nonholonomic constraints are given using Routh equations and a con- straint stabilization method. Finally, an example is used to test the proposed numerical method. The results show some dynamical behaviors of the wheeled multi-body system and its constraint stabilization effects.
基金supported by the Development of airborne gravity gradiometer(No.2017YFC0601601)open subject of Key Laboratory of Petroleum Resources Research,Institute of Geology and Geophysics,Chinese Academy of Sciences(No.KLOR2018-8)
文摘This research proposes a novel three-dimensional gravity inversion based on sparse recovery in compress sensing. Zero norm is selected as the objective function, which is then iteratively solved by the approximate zero norm solution. The inversion approach mainly employs forward modeling; a depth weight function is introduced into the objective function of the zero norms. Sparse inversion results are obtained by the corresponding optimal mathematical method. To achieve the practical geophysical and geological significance of the results, penalty function is applied to constrain the density values. Results obtained by proposed provide clear boundary depth and density contrast distribution information. The method's accuracy, validity, and reliability are verified by comparing its results with those of synthetic models. To further explain its reliability, a practical gravity data is obtained for a region in Texas, USA is applied. Inversion results for this region are compared with those of previous studies, including a research of logging data in the same area. The depth of salt dome obtained by the inversion method is 4.2 km, which is in good agreement with the 4.4 km value from the logging data. From this, the practicality of the inversion method is also validated.
基金Project supported by the National Natural Science Foundation of China (Grant No.10571116)
文摘In this paper, a branch-and-bound method for solving multi-dimensional quadratic 0-1 knapsack problems was studied. The method was based on the Lagrangian relaxation and the surrogate constraint technique for finding feasible solutions. The Lagrangian relaxations were solved with the maximum-flow algorithm and the Lagrangian bounds was determined with the outer approximation method. Computational results show the efficiency of the proposed method for multi-dimensional quadratic 0-1 knapsack problems.
基金The project supported by the National Natural Science Foundation of China (50275073)
文摘Through detailed three-dimensional (3D) finite element (FE) calculations, the out-of-plane constraints Tz along embedded center-elliptical cracks in mode I elastic plates are studied. The distributions of Tz are obtained near the crack front with aspect ratios (a/c) of 0.2, 0.4, 0.5, 0.6, 0.8 and 1.0. Tz decreases from an approximate value of Poisson ratio v at the crack tip to zero with increasing normalized radial distances (r/a) in the normal plane of the crack front line, and increases gradually when the elliptical parameter angle φ changes from 0° to 90°at the same r/a. With a/c rising to 1.0, Tz is getting nearly independent of φ and is only related to r/a. Based on the present FE calculations for Tz, empirical formulas for Tz are obtained to describe the 3D distribution of Tz for embedded center-elliptical cracks using the least squares method in the range of 0.2 ≤ a/c ≤ 1.0. These Tz results together with the corresponding stress intensity factor K are well suitable for the analysis of the 3D embedded centerelliptical crack from field, and a two-parameter K-Tz principle is proposed.
文摘Finding a suitable space is one of the most critical problems for dimensionality reduction. Each space corresponds to a distance metric defined on the sample attributes, and thus finding a suitable space can be converted to develop an effective distance metric. Most existing dimensionality reduction methods use a fixed pre-specified distance metric. However, this easy treatment has some limitations in practice due to the fact the pre-specified metric is not going to warranty that the closest samples are the truly similar ones. In this paper, we present an adaptive metric learning method for dimensionality reduction, called AML. The adaptive metric learning model is developed by maximizing the difference of the distances between the data pairs in cannot-links and those in must-links. Different from many existing papers that use the traditional Euclidean distance, we use the more generalized l<sub>2,p</sub>-norm distance to reduce sensitivity to noise and outliers, which incorporates additional flexibility and adaptability due to the selection of appropriate p-values for different data sets. Moreover, considering traditional metric learning methods usually project samples into a linear subspace, which is overstrict. We extend the basic linear method to a more powerful nonlinear kernel case so that well capturing complex nonlinear relationship between data. To solve our objective, we have derived an efficient iterative algorithm. Extensive experiments for dimensionality reduction are provided to demonstrate the superiority of our method over state-of-the-art approaches.
基金Supported by the National Natural Science Foundation of China(11372269,10902057)
文摘The theoretical solutions are obtained for the three-dimensional(3-D)stress field in an infinite isotropic elastic plate with a through-the-thickness circular hole subjected to shear load at far field by using Kane and Mindlin′s assumption based on the stress function method.Based on the present solutions,the characteristics of 3-D stress field are analyzed and the emphasis is placed on the effects of the plate thickness and Poisson′s ratio on the deviation of the present 3-D in-plane stress from the related plane stress solutions,the stress concentration and the out-of-plane constraint.The present solutions show that the stress concentration factor reaches its peak value of about 8.9% which is higher than that of the plane stress solutions.As expected,the out-of-plane stress constraint factor can reach 1on the surface of the hole when the plate is a very thick one.
文摘An explicit Bargmann symmetry constraint is computed and its associated binary nonlinearization of Lax pairs is carried out for the super NLS-MKdV hierarchy. Under the obtained symmetry constraint, the n-th flow of the super NLS-MKdV hierarchy is decomposed into two super finite-dimensional integrable Hamiltonian systems, defined over the super-symmetry manifold R4N|2N with the corresponding dynamical variables x and tn. The integrals of motion required for Liouville integrability are explicitly given.