Search

The simulation of rubber-like materials is usually based on hyperelasticity. If strain-rate dependency has to be considered viscous dampers also have to be taken into account in the rheological model. A disadvantage of such a description is time-consuming parameter identification associated with the damping constants. In this paper, a tabulated formulation is presented which allows the fast generation of input data based on uniaxial static and dynamic tensile tests at different strain rates. Unloading, i.e. forming a hysteresis loop, can also be modeled easily based on a damage formulation. We show the theoretical background and algorithmic setup of our model that has been implemented in the explicit finite element program LS-DYNA [1-3]. Apart from purely numerical examples, the validation of a soft and a hard rubber under loading and subsequent unloading at different strain rates is shown. [ABSTRACT FROM AUTHOR]

This paper is concerned with mathematical and numerical analysis of the system of radiative integral transfer equations. The existence and uniqueness of solution to the integral system is proved by establishing the boundedness of the radiative integral operators and proving the invertibility of the operator matrix associated with the system. A collocation-boundary element method is developed to discretize the differential-integral system. For the non-convex geometries, an element-subdivision algorithm is developed to handle the computation of the integrals containing the visibility factor. An efficient iterative algorithm is proposed to solve the nonlinear discrete system and its convergence is also established. Numerical experiment results are also presented to verify the effectiveness and accuracy of the proposed method and algorithm. [ABSTRACT FROM AUTHOR]

This paper defines the self-covariance property for the solutions of transferable utility cooperative games (TU-games) as a weakening of their covariance. Self-covariant solutions are positively homogenous and satisfy a "restricted" translation covariance so that feasible shifts are only the solution vectors themselves and their multipliers. A description of all non-empty, single-valued, efficient, anonymous, weakly and self-covariant solutions in the class of twoplayer TU-games is given. As demonstrated below, among them there exist just three solutions admitting consistent extensions in the Davis-Maschler sense. They are the equal share solution, the standard solution, and the constrained egalitarian solution for superadditive two-player TUgames. For the third solution mentioned, characterizations of some consistent extensions to the class of all TU-games are given. [ABSTRACT FROM AUTHOR]

We consider the problem of tensor (i.e., multidimensional array) inpainting in this paper. By using higher-order singular value decomposition, we propose an iterative algorithm that performs soft thresholding on entries of the core tensor and then reconstructs via the directional orthogonal matrices. An inpainted tensor is obtained at the end of the iteration. Simulations conducted over color images, video frames, and MR images validate that the proposed algorithm is competitive with state-of-the-art completion algorithms. The evaluation is made in terms of quality metrics and visual comparison. [ABSTRACT FROM AUTHOR]

A typical interpolation strategy for line segments consists of a transition scheme, a Look-ahead ACC/DEC scheduling, and an interpolation algorithm. In these three parts, the main computation occurs in the first and second part. Some research work has been carried out to decrease the computation in the previous literatures, but these methods occupy a lot of computing resources for the optimization process during the calculation of transition curve parameters and feed rates. Consequently, the computational efficiency of interpolation strategy is greatly reduced. To deal with the issue, a real-time interpolation strategy is proposed in this paper. In the transition scheme, a Bézier curve is utilized to smooth the line segments. Based on the relationship among the approximation error, the approximation radius, and the transition curve, the curve can be directly generated when the approximation error is given. In the ACC/DEC scheduling, a 3-segment feed rate profile with jerk continuity is constructed. Meanwhile, a Look-ahead planning based on Backward Scanning and Forward Revision (BSFR) algorithm is utilized to eliminate redundant computation. Compared with Zhao’s and Shi’s strategy, the proposed strategy has the merits of C2 and G2 continuity for the tool path, jerk continuity for the tool movement, and distinguished real-time performance for interpolation. The experiments of 3D pentagram and 2D butterfly are carried out with different strategies and their results demonstrate that the interpolation efficiency can be greatly improved with the proposed strategy. [ABSTRACT FROM AUTHOR]

This paper proposes a stream selection algorithm to deal with the interference and increase the sum capacity in the context of heterogeneous networks where cells of different types coexist. Due to the different transmit powers between the macro and small cells, interference levels are also different. Since the small cells are densely deployed, fully connected interference network between small cells is considered in this paper. The proposed algorithm performs the selection of a stream sequence among a predetermined set of sequences. Those selected sequences are the ones that mostly contribute to the sum rate when performing the exhaustive search. As a consequence, since the search space is reduced, the complexity is significantly decreased. These stream sequences form a regular structure where the first stream is associated to a pico user. Another distinguishing property of the proposed approach is that the stream sequences include at least one stream for each user. When selecting the streams, the channel matrices of the unselected streams are projected orthogonally to the virtual transmit and virtual receive channels of the selected stream in order to align the interference in the null space of these virtual channels. The performance evaluations are carried out by considering different scenarios with different numbers and locations of pico cells. It is shown that the proposed method can significantly reduce the computational complexity while achieving a very close performance to the exhaustive search. [ABSTRACT FROM AUTHOR]

In this paper, we introduce and analyze a general iterative algorithm for finding a common solution of a finite family of mixed equilibrium problems, a general system of generalized equilibria and a fixed point problem of nonexpansive mappings in a real Hilbert space. Under some appropriate conditions, we derive the strong convergence of the sequence generated by the proposed algorithm to a common solution, which also solves some optimization problem. The result presented in this paper improves and extends some corresponding ones in the earlier and recent literature. [ABSTRACT FROM AUTHOR]

In this paper, we introduce and analyze a multi-step hybrid extragradient algorithm by combining Korpelevich's extragradient method, the viscosity approximation method, the hybrid steepest-descent method, Mann's iteration method and the gradient-projection method (GPM) with regularization in the setting of infinite-dimensional Hilbert spaces. It is proven that, under appropriate assumptions, the proposed algorithm converges strongly to a solution of the convex minimization problem (CMP) with constraints of several problems: finitely many generalized mixed equilibrium problems (GMEPs), finitely many variational inclusions, and the fixed point problem of a strictly pseudocontractive mapping. In the meantime, we also prove the strong convergence of the proposed algorithm to the unique solution of a hierarchical variational inequality problem (over the fixed point set of a strictly pseudocontractive mapping) with constraints of finitely many GMEPs, finitely many variational inclusions and the CMP. The results presented in this paper improve and extend the corresponding results announced by many others. [ABSTRACT FROM AUTHOR]

In this paper, we construct a second order algorithm based on the spectral deferred correction method for the time-dependent magnetohydrodynamics flows at a low magnetic Reynolds number. We present a complete theoretical analysis to prove that this algorithm is unconditionally stable, consistent and second order accuracy. Finally, two numerical examples are given to illustrate the convergence and effectiveness of our algorithm. [ABSTRACT FROM AUTHOR]

This paper proposes and evaluates a near admissible heuristic search strategy and its application to a kind of flexible manufacturing system (FMS) scheduling in a Petri net framework. Petri nets can concisely model the strict precedence constraint, multiple kinds of resources, and concurrent activities. To cope with the complexities for scheduling of FMS with alternative routings, this paper proposes an admissible heuristic function based on the execution of P-timed Petri nets and presents an improved dynamic weighting A* strategy using the proposed heuristic function. The search scheme does not need to predict the depth of solution in advance and the quality of the search result is also controllable. Some numerical experiments are carried out to demonstrate usefulness of the algorithm. [ABSTRACT FROM AUTHOR]

This paper discusses characteristics unique in the two-phase high-frequency-testing (HFT) environment in semiconductor manufacturing. We believe this paper is the first to define, formalize, and analyze the decision making problem associated with the two-phase HFT. Specifically, this paper defines the problem of minimizing the total HFT capacity usage by systematically finding the optimal number of preliminary HFT bin-1 chips subject to the main HFT, with the existence of the target service rate. We also propose a heuristic algorithm that exploits the special structure of the problem for efficiently obtaining a near-optimal solution. Finally, a numerical analysis and a case study have been conducted to gain more insights on the problem structure and the proposed algorithm. [ABSTRACT FROM AUTHOR]

In this paper, we propose two iterative methods for finding a common solution of a finite family of equilibrium problems for pseudomonotone bifunctions. The first is a parallel hybrid extragradientcutting algorithm which is extended from the previously known one for variational inequalities to equilibrium problems. The second is a new cyclic hybrid extragradient-cutting algorithm. In the cyclic algorithm, using the known techniques, we can perform and develop practical numerical experiments. [ABSTRACT FROM AUTHOR]

In this paper, a modified feasible semi-smooth asymptotically Newton method for nonlinear complementarity problems is proposed. The global convergence of the method and superlinear convergence are proved under some suitable assumptions. Numerical experiments are included to highlight the efficacy of the modified algorithm. [ABSTRACT FROM AUTHOR]

In this paper, we propose a new iterative algorithm to approximate fixed points of total asymptotically nonexpansive mappings in $\operatorname{CAT}(0)$ spaces. We also provide two examples to illustrate the convergence behavior of the proposed algorithm and numerically compare the convergence of the proposed iteration scheme with the existing schemes. [ABSTRACT FROM AUTHOR]

In this paper, we attempt to introduce an algorithm which considers not only time factor but also cost, risk, and quality criteria to determine the critical path under fuzzy environment. In this algorithm first, decision makers allocate time, cost, risk, and quality to each activity. We are lacking in data and information, so in the proposed algorithm, the ratings of each activity and the weight of each criterion are described by fuzzy numbers and linguistic variables, which can be expressed in triangular fuzzy number. Linguistic variables are applied to represent the intensity of preferences of one criterion over another. Then, we add up triangular fuzzy numbers to determine the final evaluation value of each criterion for paths. Next, we use fuzzy TOPSIS, a technique for order preferences by similarity to an ideal solution, a method proposed by the authors in another paper, to choose the best alternative. Finally, numerical example is solved to illustrate the procedure of proposed method at the end of this paper. [ABSTRACT FROM AUTHOR]

In any machining fixture, the workpiece elastic deformation caused during machining influences the dimensional and form errors of the workpiece. Placing each locator and clamp in an optimal place can minimize the elastic deformation of the workpiece, which in turn minimizes the dimensional and form errors of the workpiece. Design of fixture configuration (layout) is a procedure to establish the workpiece–fixture contact through optimal positioning of clamping and locating elements. In this paper, an ant colony algorithm (ACA) based discrete and continuous optimization methods are applied for optimizing the machining fixture layout so that the workpiece elastic deformation is minimized. The finite element method (FEM) is used for determining the dynamic response of the workpiece caused due to machining and clamping forces. The dynamic response of the workpiece is simulated for all ACA runs. This paper proves that the ACA-based continuous fixture layout optimization method exhibits the better results than that of ACA-based discrete fixture layout optimization method. [ABSTRACT FROM AUTHOR]

Modern constitutive models have the potential to improve the quality of engineering calculations involving non-linear anisotropic materials. The adoption of complex models in practice, however, depends on the availability of reliable and accurate solution methods for the stress point integration problem. This paper presents a modular implementation of explicit Runge–Kutta methods with error control, that is suitable for use, without change, with any rate-type constitutive model. The paper also shows how the complications caused by the algebraic constraint of conventional plasticity are resolved through a simple subloading modification. With this modification any rate-independent model can be implemented without difficulty, using the integration module as an accurate and robust standard procedure. The effectiveness and efficiency of the method are demonstrated through a comparative evaluation of second and fifth-order formulas, applied to a complex constitutive model for natural clay, full details of which are given. [ABSTRACT FROM AUTHOR]

A sub-cycling integration algorithm (or named multi-time-steps integration algorithm), which has been successfully applied to FEM dynamical analysis, was firstly presented by Belytschko et al. (Comput Methods Appl Mech Eng 17/18:259–275, 1979). However, the problem of how to apply this type of algorithm to flexible multi-body dynamics (FMD) problems still lacks investigation up to now. Similar to the region-partitioning method used in FEM, this paper presents a central-difference-based sub-cycling integral method by decomposing the variables of an FMD equation into several groups and adopting different integral step sizes to each group of the variables. Based on the condensed form of an FMD equation, a group of common update formulae and a sub-step update formula, which constitute the sub-cycling together, are established in the paper. Furthermore, an implementation flowchart of the sub-cycling is presented. Stability of the sub-cycling will be analyzed and numerical examples will be performed to verify availability and precision of the sub-cycling in part II of the paper. [ABSTRACT FROM AUTHOR]

Numerical stability is an important issue for any integral procedure. Since sub-cycling algorithm has been presented by Belytschko et al. (Comput Methods Appl Mech Eng 17/18: 259–275, 1979), various kinds of these integral procedures were developed in later 20 years and their stability were widely studied. However, on how to apply the sub-cycling to flexible multi-body dynamics (FMD) is still a lack of investigation up to now. A particular sub-cycling algorithm for the FMD based on the central difference method was introduced in detail in part I (Miao et al. in Comp Mech doi: 10.1007/s00466-007-0183-9) of this paper. Adopting an integral approximation operator method, stability of the presented algorithm is transformed to a generalized eigenvalue problem in the paper and is discussed by solving the problem later. Numerical examples are performed to verify the availability and efficiency of the algorithm further. [ABSTRACT FROM AUTHOR]

This paper explores various aspects of the image decomposition problem using modern variational techniques. We aim at splitting an original image f into two components u and ρ, where u holds the geometrical information and ρ holds the textural information. The focus of this paper is to study different energy terms and functional spaces that suit various types of textures. Our modeling uses the total-variation energy for extracting the structural part and one of four of the following norms for the textural part: L 2, G, L 1 and a new tunable norm, suggested here for the first time, based on Gabor functions. Apart from the broad perspective and our suggestions when each model should be used, the paper contains three specific novelties: first we show that the correlation graph between u and ρ may serve as an efficient tool to select the splitting parameter, second we propose a new fast algorithm to solve the TV − L 1 minimization problem, and third we introduce the theory and design tools for the TV-Gabor model. [ABSTRACT FROM AUTHOR]

In this paper a three-dimensional Voronoi cell finite element model is developed for analyzing heterogeneous materials containing a dispersion of ellipsoidal inclusions or voids in the matrix. The paper starts with a description of 3D tessellation of a domain with ellipsoidal heterogeneities, to yield a 3D mesh of Voronoi cells containing the heterogeneities. A surface based tessellation algorithm is developed to account for the shape and size of the ellipsoids in point based tessellation methods. The 3D Voronoi cell finite element model, using the assumed stress hybrid formulation, is developed for determining stresses and displacements in a linear elastic material domain. Special stress functions that introduce classical Lamé functions in ellipsoidal coordinates are implemented to enhance solution convergence. Numerical methods for implementation of algorithms and yielding stable solutions are discussed. Numerical examples are conducted with inclusions and voids to demonstrate the effectiveness of the model. [ABSTRACT FROM AUTHOR]

This paper proposes a numerical procedure that can predict the minimum duration of a rectangular, half-cycle sine wave, linearly increasing and triangular single-pulse excitation required to overturn the rigid block resting on a moving base. Since the linearization assumption in the derivation of an analytical solution, which has already been used by Housner and other researchers, cannot be valid for a block with various slenderness ratios and dimensions, a derivation of numerical solutions which would be valid for all types of block is necessary. The proposed numerical procedure takes into account the possibility of sliding and uplifting of the block from the base. The numerical algorithm for the proposed numerical procedure in a commercial code C++ is given in the 'Appendix.' [ABSTRACT FROM AUTHOR]

In this paper, we present a numerical method for solving reaction-diffusion equations on one dimensional branched structures. Through the use of a simple domain decomposition scheme, the many branches are decoupled so that the equations can be solved as a system of smaller problems that are tri-diagonal. This technique allows for locally adaptive time stepping, in which the time step used in each branch is determined by local activity. Though the method is presented in the specific context of electrical activity in neural systems, it is sufficiently general that it can be applied to other classes of reaction-diffusion problems and higher dimensions. Information in neurons, which can be effectively modeled as one-dimensional branched structures, is carried in the form of electrical impulses called action potentials. The model equations, based on the Hodgkin-Huxley cable equations, are a set of reaction equations coupled to a single diffusion process. Locally adaptive time stepping schemes are well suited to neural simulations due to the spatial localization of activity. The algorithm significantly reduces the computational cost compared to existing methods, especially for large scale simulations. [ABSTRACT FROM AUTHOR]

In this paper, we investigate the coupling of reduced models for the simulation of structures involving localized geometrical details. Herein, we use the Arlequin method, originally designed to deal with multimodel and multiscale analyses of mechanical problems, to mix reduced models built using the proper generalized decomposition. Instead of solving the global coupled problem in a monolithic way, the LATIN strategy is used to propose a decoupled algorithm. The numerical examples demonstrate the feasibility of the approach and in particular its potentiality in terms of flexibility. [ABSTRACT FROM AUTHOR]

In this paper a nonlinear backward parabolic problem in one dimensional space is considered. Using a suitable iterative algorithm, the problem is converted to a linear backward parabolic problem. For the corresponding problem, the backward finite differences method with suitable grid size is applied. It is shown that if the coefficients satisfy some special conditions, this algorithm not only is convergent, but also is conditionally stable. Moreover, it is proved that the estimated values converge to the exact solution of the problem. All these approaches examined in some numerical examples. corresponding theorems for the convergency and stability of the solution are studied. [ABSTRACT FROM AUTHOR]

This paper analyzes the routing problem in order to design automatic train operation in a station. The authors consider the routing problem for the receipt-departure of trains and their shunting movements. It is noted that, in the case of stations with a branched track infrastructure, the number of all admissible routes becomes large, leading to bulky computations in the search procedure. A certain method is proposed to partition stations into zones for reducing the amount of computations. And finally, the authors pose the graph cutting problem of a railway station network with due account of track infrastructure features and present its solution algorithm. [ABSTRACT FROM AUTHOR]

This paper presents a method to find all intersections between curved lines such as structural line elements and finite element meshes with intentions to generate smaller, non-compatible, line cells (e.g. bar elements) between crossings. The intersection finding algorithm works for two and three-dimensional meshes constituted by linear, quadratic or higher order elements. Using the proposed algorithm, meshes can then be automatically prepared for finite element analyses with techniques for embedding elements within others or analyses that require lines within solids. The application of the method is demonstrated by a number of numerical examples illustrating its capabilities in handling complex geometries, relative speed and convenience. [ABSTRACT FROM AUTHOR]

In the paper, first, we introduce a new hybrid projection algorithm and present its strong convergence theorem. Next, we analyze different hybrid algorithms in computing and conclude that our proposed algorithm has an advantage. Finally, the numerical experiments validate the efficiency and advantages of the new algorithm. [ABSTRACT FROM AUTHOR]

The split variational inclusion problem is an important problem, and it is a generalization of the split feasibility problem. In this paper, we present a descent-conjugate gradient algorithm for the split variational inclusion problems in Hilbert spaces. Next, a strong convergence theorem of the proposed algorithm is proved under suitable conditions. As an application, we give a new strong convergence theorem for the split feasibility problem in Hilbert spaces. Finally, we give numerical results for split variational inclusion problems to demonstrate the efficiency of the proposed algorithm. [ABSTRACT FROM AUTHOR]

Consider a data set of vector-valued observations that consists of noisy inliers, which are explained well by a low-dimensional subspace, along with some number of outliers. This work describes a convex optimization problem, called reaper, that can reliably fit a low-dimensional model to this type of data. This approach parameterizes linear subspaces using orthogonal projectors and uses a relaxation of the set of orthogonal projectors to reach the convex formulation. The paper provides an efficient algorithm for solving the reaper problem, and it documents numerical experiments that confirm that reaper can dependably find linear structure in synthetic and natural data. In addition, when the inliers lie near a low-dimensional subspace, there is a rigorous theory that describes when reaper can approximate this subspace. [ABSTRACT FROM AUTHOR]

Let $${\mathcal {A}}\subseteq {\mathbb {N}}^n$$ be a finite set, and $$K\subseteq {\mathbb {R}}^n$$ be a compact semialgebraic set. An $${\mathcal {A}}$$ -truncated multisequence ( $${\mathcal {A}}$$ -tms) is a vector $$y=(y_{\alpha })$$ indexed by elements in $${\mathcal {A}}$$ . The $${\mathcal {A}}$$ -truncated $$K$$ -moment problem ( $${\mathcal {A}}$$ -TKMP) concerns whether or not a given $${\mathcal {A}}$$ -tms $$y$$ admits a $$K$$ -measure $$\mu $$ , i.e., $$\mu $$ is a nonnegative Borel measure supported in $$K$$ such that $$y_\alpha = \int _K x^\alpha \mathtt {d}\mu $$ for all $$\alpha \in {\mathcal {A}}$$ . This paper proposes a numerical algorithm for solving $${\mathcal {A}}$$ -TKMPs. It aims at finding a flat extension of $$y$$ by solving a hierarchy of semidefinite relaxations $$\{(\mathtt {SDR})_k\}_{k=1}^\infty $$ for a moment optimization problem, whose objective $$R$$ is generated in a certain randomized way. If $$y$$ admits no $$K$$ -measures and $${\mathbb {R}}[x]_{{\mathcal {A}}}$$ is $$K$$ -full (there exists $$p \in {\mathbb {R}}[x]_{{\mathcal {A}}}$$ that is positive on $$K$$ ), then $$(\mathtt {SDR})_k$$ is infeasible for all $$k$$ big enough, which gives a certificate for the nonexistence of representing measures. If $$y$$ admits a $$K$$ -measure, then for almost all generated $$R$$ , this algorithm has the following properties: i) we can asymptotically get a flat extension of $$y$$ by solving the hierarchy $$\{(\mathtt {SDR})_k\}_{k=1}^\infty $$ ; ii) under a general condition that is almost sufficient and necessary, we can get a flat extension of $$y$$ by solving $$(\mathtt {SDR})_k$$ for some $$k$$ ; iii) the obtained flat extensions admit a $$r$$ -atomic $$K$$ -measure with $$r\le |{\mathcal {A}}|$$ . The decomposition problems for completely positive matrices and sums of even powers of real linear forms, and the standard truncated $$K$$ -moment problems, are special cases of $${\mathcal {A}}$$ -TKMPs. They can be solved numerically by this algorithm. [ABSTRACT FROM AUTHOR]

In this paper, a coupled hydro-mechanical formulation is developed for deformable porous media subjected to crack interfaces in the framework of extended finite element method. Governing equations of the porous medium consist of the momentum balance of the bulk together with the momentum balance and continuity equations of the fluid phase, known as [InlineEquation not available: see fulltext.] formulation. The discontinuity in fractured porous medium is modeled for both opening and closing modes that results in the fluid flow within the fracture, and/or contact behavior at the crack edges. The fluid flow through the fracture is assumed to be viscous and is modeled by employing the Darcy law in which the permeability of fracture is obtained using the cubic law. The contact condition in fractured porous medium is handled by taking the advantage from two different algorithms of LATIN method and penalty algorithm. The effect of contact on fluid phase is employed by considering no leak-off from/into the porous medium. The nonlinearity of coupled equations produced due to opening and closing modes is carried out using an iterative algorithm in the Newton-Raphson procedure. Finally, several numerical examples are solved to illustrate the performance of proposed X-FEM method for hydro-mechanical behavior of fractured porous media with opening and closing modes. [ABSTRACT FROM AUTHOR]

In this paper, we introduce two iterative algorithms for finding common solutions of a finite family of variational inclusions for maximal monotone and inverse-strongly monotone mappings with the constraints of two problems: a generalized mixed equilibrium problem and a common fixed point problem of an infinite family of nonexpansive mappings and an asymptotically strict pseudocontractive mapping in the intermediate sense in a real Hilbert space. We prove some strong and weak convergence theorems for the proposed iterative algorithms under suitable conditions. [ABSTRACT FROM AUTHOR]

In the present paper, we propose several kinds of adaptively relaxed iterative algorithms with a new step size for solving the split feasibility problem in real Hilbert spaces. The proposed algorithms never terminate, while the known algorithms existing in the literature may terminate. Several weak and strong convergence theorems of the proposed algorithms have been established. Some numerical experiments are also included to illustrate the effectiveness of the proposed algorithms. [ABSTRACT FROM AUTHOR]

The multiple-sets split equality problem (MSSEP) requires finding a point x ∈∩Ni=1 Ci, y ∈∩Mj=1 Qj, such that Ax = By, where N and M are positive integers, {C1, C2,..., CN} and {Q1,Q2,...,QM} are closed convex subsets of Hilbert spaces H1, H2, respectively, and A : H1→H3, B : H2→H3 are two bounded linear operators. When N = M = 1, the MSSEP is called the split equality problem (SEP). If let B = I, then the MSSEP and SEP reduce to the well-known multiple-sets split feasibility problem (MSSFP) and split feasibility problem (SFP), respectively. Recently, some authors proposed many algorithms to solve the SEP and MSSEP. However, to implement these algorithms, one has to find the projection on the closed convex sets, which is not possible except in simple cases. One of the purposes of this paper is to study the SEP and MSSEP for a family of quasi-nonexpansive multi-valued mappings in the framework of infinite-dimensional Hilbert spaces, and propose an algorithm to solve the SEP and MSSEP without the need to compute the projection on the closed convex sets. [ABSTRACT FROM AUTHOR]

Spatial outlier detection is an important research problem that has received much attentions in recent years. Most existing approaches are designed for numerical attributes, but are not applicable to categorical ones (e.g., binary, ordinal, and nominal) that are popular in many applications. The main challenges are the modeling of spatial categorical dependency as well as the computational efficiency. This paper presents the first outlier detection framework for spatial categorical data. Specifically, a new metric, named as Pair Correlation Ratio (PCR), is measured for each pair of category sets based on their co-occurrence frequencies at specific spatial distance ranges. The relevances among spatial objects are then calculated using PCR values with regard to their spatial distances. The outlierness for each object is defined as the inverse of the average relevance between an object and its spatial neighbors. Those objects with the highest outlier scores are returned as spatial categorical outliers. A set of algorithms are further designed for single-attribute and multi-attribute spatial categorical datasets. Extensive experimental evaluations on both simulated and real datasets demonstrated the effectiveness and efficiency of our proposed approaches. [ABSTRACT FROM AUTHOR]

In this paper, we present an improved one-step lookahead A* algorithm for the scheduling problem for parallel-mode cluster tools, where n wafer types are simultaneously processed on a cluster tool with m chambers. Specifically, we suggest a refined cost function, priority-based node selection, N-chamber cycle detection to identify deadlock situations, and a one-step lookahead method to reduce the number of unsafe states generated. Using numerical experiments, we show that the proposed algorithm leads to a significant reduction in the number of explored states and execution time necessary to find an optimal solution. [ABSTRACT FROM AUTHOR]

Compressed sensing ensures the accurate reconstruction of sparse signals from far fewer samples than required in the classical Shannon-Nyquist theorem. In this paper, a generalized hard thresholding pursuit (GHTP) algorithm is presented that can recover unknown vectors without the sparsity level information. We also analyze the convergence of the proposed algorithm. Numerical experiments are given for synthetic and real-world data to illustrate the validity and the good performance of the proposed algorithm. [ABSTRACT FROM AUTHOR]

In this paper we propose a numerical procedure for the prediction of capillary forces in polydisperse granular assemblies at a degree of moisture content that corresponds to the so-called pendular regime. The capillary force model is adopted within the Laplace-Young framework with a toroidal approximation of the liquid bridge geometry. Governing equations are first derived in a form that highlights the role of intrinsic parameters such as inter-particle separation distance, ratio of particle radii and liquid volume. A proper scaling of these equations is adopted so that the solution applies to any particle pair configuration. Numerical integration algorithms are provided in a way that facilitates implementations in macroscopic procedures for computer simulations. A qualitative evaluation is undertaken to highlight model predictions of the effects on capillary force of various intrinsic parameters that characterise the particle pair and liquid bridge. The model is validated against the experimental results provided by Willet et al. (Langmuir 16:9396-9405, ) for a wide range of liquid volumes and particle-pair polydispersity. [ABSTRACT FROM AUTHOR]

In this paper, we present a dynamic refinement algorithm for the smoothed particle Hydrodynamics (SPH) method. An SPH particle is refined by replacing it with smaller daughter particles, which positions are calculated by using a square pattern centered at the position of the refined particle. We determine both the optimal separation and the smoothing distance of the new particles such that the error produced by the refinement in the gradient of the kernel is small and possible numerical instabilities are reduced. We implemented the dynamic refinement procedure into two different models: one for free surface flows, and one for post-failure flow of non-cohesive soil. The results obtained for the test problems indicate that using the dynamic refinement procedure provides a good trade-off between the accuracy and the cost of the simulations. [ABSTRACT FROM AUTHOR]

In this paper we introduce a new smoothing function and show that it is coercive under suitable assumptions. Based on this new function, we propose a smoothing Newton method for solving the second-order cone complementarity problem (SOCCP). The proposed algorithm solves only one linear system of equations and performs only one line search at each iteration. It is shown that any accumulation point of the iteration sequence generated by the proposed algorithm is a solution to the SOCCP. Furthermore, we prove that the generated sequence is bounded if the solution set of the SOCCP is nonempty and bounded. Under the assumption of nonsingularity, we establish the local quadratic convergence of the algorithm without the strict complementarity condition. Numerical results indicate that the proposed algorithm is promising. [ABSTRACT FROM AUTHOR]

The second Zagreb index of a graph G is an adjacency-based topological index, which is defined as Σuv∊E(G)(d(u)d(v)), where uv is an edge of G, d(u) is the degree of vertex u in G. In this paper, we consider the second Zagreb index for bipartite graphs. Firstly, we present a new definition of ordered bipartite graphs, and then give a necessary condition for a bipartite graph to attain the maximal value of the second Zagreb index. We also present an algorithm for transforming a bipartite graph to an ordered bipartite graph, which can be done in O(n2+n21) time for a bipartite graph B with a partition ∣X∣ = n1 and ∣Y∣ = n2. [ABSTRACT FROM AUTHOR]

In this paper we present a new approach for solving elastoplastic problems as second order cone complementarity problems (SOCCPs). Specially, two classes of elastoplastic problems, i.e. the J plasticity problems with combined linear kinematic and isotropic hardening laws and the Drucker-Prager plasticity problems with associative or non-associative flow rules, are taken as the examples to illustrate the main idea of our new approach. In the new approach, firstly, the classical elastoplastic constitutive equations are equivalently reformulated as second order cone complementarity conditions. Secondly, by employing the finite element method and treating the nodal displacements and the plasticity multiplier vectors of Gaussian integration points as the unknown variables, we obtain a standard SOCCP formulation for the elastoplasticity analysis, which enables the using of general SOCCP solvers developed in the field of mathematical programming be directly available in the field of computational plasticity. Finally, a semi-smooth Newton algorithm is suggested to solve the obtained SOCCPs. Numerical results of several classical plasticity benchmark problems confirm the effectiveness and robustness of the SOCCP approach. [ABSTRACT FROM AUTHOR]

In this paper, an efficient algorithm is established for computing the maximum (minimum) angular separation ρ( ρ), the corresponding apparent position angles ( $\theta|_{\rho_{\rm max}}$, $\theta|_{\rho_{\rm min}}$) and the individual masses of visual binary systems. The algorithm uses Reed's formulae (1984) for the masses, and a technique of one-dimensional unconstrained minimization, together with the solution of Kepler's equation for $(\rho_{\rm max}, \theta|_{\rho_{\rm max}})$ and $(\rho_{\rm min}, \theta|_{\rho_{\rm min}})$. Iterative schemes of quadratic coverage up to any positive integer order are developed for the solution of Kepler's equation. A sample of 110 systems is selected from the Sixth Catalog of Orbits (Hartkopf et al.). Numerical studies are included and some important results are as follows: (1) there is no dependence between ρ and the spectral type and (2) a minor modification of Giannuzzi's () formula for the upper limits of ρ functions of spectral type of the primary. [ABSTRACT FROM AUTHOR]

Fast multipole boundary element methods (FMBEMs) are developed based on the couple of fast multipole algorithm and generalized minimal residual algorithm. The FMBEMs improve the efficiency of conventional BEMs, accelerate the computing, enlarge the solving scale, and it is applied in various engineering fields. The paper tried to do a brief review for the FMBEMs, and focus on the description of basic principles and applications in rolling engineering. The basic principles and main frameworks of two typical methods of FMBEMs (sphere harmonic function multipole BEM and Taylor series multipole BEM) are briefly described, and then the key numerical iterative and preconditioning techniques suitable for the FMBEMs are introduced. The typical numerical examples are presented, including the elasticity problems, the elastic contact problems and the elastoplasticity problems, etc. The validity and effectiveness of FMBEMs are effectively illustrated by engineering analysis examples. The numerical results suggest that the FMBEMs are suitable for the analysis and solution of large scale rolling engineering problems. The implementation process of numerical analysis can provide useful reference for the applications in other engineering fields. [ABSTRACT FROM AUTHOR]

This paper studies adaptive finite element methods (AFEMs), based on piecewise linear elements and newest vertex bisection, for solving second order elliptic equations with piecewise constant coefficients on a polygonal domain Ω⊂ℝ. The main contribution is to build algorithms that hold for a general right-hand side f∈ H( Ω). Prior work assumes almost exclusively that f∈ L( Ω). New data indicators based on local H norms are introduced, and then the AFEMs are based on a standard bulk chasing strategy (or Dörfler marking) combined with a procedure that adapts the mesh to reduce these new indicators. An analysis of our AFEM is given which establishes a contraction property and optimal convergence rates N with 0< s≤1/2. In contrast to previous work, it is shown that it is not necessary to assume a compatible decay s<1/2 of the data estimator, but rather that this is automatically guaranteed by the approximability assumptions on the solution by adaptive meshes, without further assumptions on f; the borderline case s=1/2 yields an additional factor log N. Computable surrogates for the data indicators are introduced and shown to also yield optimal convergence rates N with s≤1/2. [ABSTRACT FROM AUTHOR]

In this paper, we focus on real-life settings that require the development of new models of flowshop scheduling problems, where job processing times can increase with the number of processed jobs due to the aging effect and decrease by the allocation of additional resource. We analyse the makespan minimization flowshop problem with such model and also with the aging effect only. We prove that the considered problems and their special cases are still polynomially solvable under given conditions, and on their basis, we provide optimal polynomial time solution algorithms. [ABSTRACT FROM AUTHOR]

Introducing beads into thin-walled structures for enhancing bending stiffness is a common practice and has been used for many decades. Typically, forms and patterns of such beads are rather based on experience as well as practical results than on analytical or numerical investigations. Recently computational methods have been developed for designing or optimizing beads. The paper at hand contributes to these attempts. It deals with increasing the buckling resistance as well as the fundamental frequency of thin-walled structures by systematic bead design. The design strategy is based on the idea to disturb the buckling mode and the fundamental vibration mode, respectively, of a structure in an efficient way by laying beads along the direction of the occurring maximum principal curvature of the corresponding mode shape. Since the beads influence this mode shape, an incremental approach is required. Numerical investigations into plate and profile structures are performed by applying this new iterative design procedure in combination with the finite element method. [ABSTRACT FROM AUTHOR]

Measuring the bridge frequencies indirectly from a test vehicle is a potentially powerful technique for its advantage of efficiency and economy. However, the existence of road surface roughness may excite the vehicle-related frequencies to a higher level of amplitude than the bridge-related ones, making it difficult to identify the bridge frequencies. In order to overcome this drawback, this paper proposes the use of two connected vehicles to remove or reduce the blurring effect of road surface roughness, so as to ensure the successful identification of bridge frequencies. Underlaid by this idea, a Fourier-based algorithm for extracting the bridge frequencies is presented. Numerical studies are conducted for several different combinations of vehicle masses and stiffnesses, as well as for vehicle spacing. Through such a study, the feasibility of using two connected vehicles to extract the bridge frequencies is confirmed. [ABSTRACT FROM AUTHOR]

Often in the engineering practice, cutting tools bounded by primary peripheral surfaces of revolution are used because of their effectiveness. Among these, ring and tangential tools can be used for the generation of constant pitch cylindrical helical surfaces. In this paper, we present an algorithm for the profiling of these types of tools. The algorithm is based on the topological representation of the tool's primary peripheral surface. The main goal is to devise a methodology for the profiling of tools whose surfaces are reciprocally enveloping with cylindrical helical surfaces. We present a numerical example for the numerical determination of the axial section form for this type of tools. The application method for this algorithm was developed in the CATIA graphical design environment within which the procedure is developed as a vertical application. In addition, we present a solution for the shape correction of the tool's axial cross-section by considering the existence of singular points on the profile of the helical surface to be generated where multiple normals to the surface exist. [ABSTRACT FROM AUTHOR]