Tensor‐based Gauss–Jacobi numerical integration for high‐order mass and stiffness matrices

In this work, we choose the points and weights of the Gauss–Jacobi, Gauss–Radau–Jacobi and Gauss–Lobatto–Jacobi quadrature rules that optimize the number of operations for the mass and stiffness matrices of the high‐order finite element method. The procedure is particularly applied to the mass and stiffness matrices using the tensor‐based nodal and modal shape functions given in (Int. J. Numer. Meth. Engng 2007; 71 (5):529–563). For square and hexahedron elements, we show that it is possible to use tensor product of the 1D mass and stiffness matrices for the Poisson and elasticity problem. For the triangular and tetrahedron elements, an analogous analysis given in (Int. J. Numer. Meth. Engng 2005; 63 (2):1530–1558) was considered for the selection of the optimal points and weights for the stiffness matrix coefficients for triangles and mass and stiffness matrices for tetrahedra. Copyright © 2009 John Wiley & Sons, Ltd.     

Robust H 8 output feedback control for general nonlinear systems with structured uncertainty

Abstract

In this paper, the robust H ⁸ output feedback control problem for general nonlinear systems with L ₂‐norm‐bounded structured uncertainties is considered. Sufficient conditions for the solvability of robust performance synthesis problems are represented in terms of two Hamilton‐Jacobi inequalities with n independent variables. Based on these conditions, a state space characterization of a robust H ⁸ output feedback controller solving the considered problem is proposed. An example is provided for illustration.     

Partition of unity for localization in implicit a posteriori finite element error control for linear elasticity

The partition of unity for localization in adaptive finite element method (FEM) for elliptic partial differential equations has been proposed in Carstensen and Funken (SIAM J. Sci. Comput. 2000; 21 : 1465–1484) and is applied therein to the Laplace problem. A direct adaptation to linear elasticity in this paper yields a first estimator η_L based on patch‐oriented local‐weighted interface problems. The global Korn inequality with a constant C_Korn yields reliability for any finite element approximation u_h to the exact displacement u. In order to localize this inequality further and so to involve the global constant C_Korn directly in the local computations, we deduce a new error estimator µ_L. The latter estimator is based on local‐weighted interface problems with rigid body motions (RBM) as a kernel and so leads to effective estimates only if RBM are included in the local FE test functions. Therefore, the excluded first‐order FEM has to be enlarged by RBM, which leads to a partition of unit method (PUM) with RBM, called P₁+RBM or to second‐order FEMs, called P₂ FEM. For P₁+RBM and P₂ FEM (or even higher‐order schemes) one obtains the sharper reliability estimate . Efficiency holds in the strict sense of . The local‐weighted interface problems behind the implicit error estimators η_L and µ_L are usually not exactly solvable and are rather approximated by some FEM on a refined mesh and/or with a higher‐order FEM. The computable approximations are shown to be reliable in the sense of . The oscillations are known functions of the given data and higher‐order terms if the data are smooth for first‐order FEM. The mathematical proofs are based on weighted Korn inequalities and inverse estimates combined with standard arguments. The numerical experiments for uniform and adapted FEM on benchmarks such as an L‐shape problem, Cook's membrane, or a slit problem validate the theoretical estimates and also concern numerical bounds for C_Korn and the locking phenomena. Copyright © 2007 John Wiley & Sons, Ltd.     

Finite element analysis for a class of nonlinear variational inequalities

In this paper, we prove that the error estimates for the finite element approximation of a class of variational inequalities arising in elastostatics is of order h in the energy norm. In fact, our estimates improve all of the previously known error estimates for elliptic variational inequalities.     

On computational methods for variational inequalities

In this note, we focus on optimised mesh design for the Finite Element (FE) method for variational inequalities using global norm estimates for local error control. The strategies are based on the so called dual-weighted-residual (DWR) approach to a posteriori error control for FE-schemes (see, e.g., Rannacher et al. [19, 6, 2]), where error control for the primal problem is established by solving an auxiliary (dual) problem. In this context we blamed (cf. e.g., Rannacher and Suttmeier [18, 19]) global norm estimates being not that useful in applications. But having a closer look at the DWR-concept, one observes that in fact global (energy) error bounds can be employed to establish local error control. Our ideas and techniques are illustrated at the socalled obstacle problem.It turns out, that reliable and efficient energy error control is one main ingredient to establish useful a posteriori error bounds for local quantities. Therefore, in addition, we derive an unified approach to a posteriori error control in the energy norm for elliptic variational inequalities of first kind. Eventually, this framework is applied to Signorinis problem.     

Robust mixed H 2/H ∞ controller for uncertain neutral state‐input delayed systems via LMI optimization approach

Abstract

In this paper, the problem of designing a robust mixed H ₂/H _∞ controller for a class of uncertain neutral state‐input delays system is considered. Based on Lyapunov‐Krasovskii functional theory, a delay‐dependent criterion is derived for the existence of a desired mixed H ₂/H _∞ controller, which can be easily constructed by certain feasible linear matrix inequalities (LMIs). Furthermore, a convex optimization problem is formulated to solve for a robust mixed H ₂/H _∞ controller which achieves the minimization of an upper bound of the closed‐loop H ₂ perforance measure.     

An application of modified reductive perturbation method to long water waves

In this work, we extended the application of “the modified reductive perturbation method” to long water waves and obtained the governing equations as the KdV hierarchy. Seeking a localized travelling wave solutions to these evolution equations we determined the scale parameter c₁so as to remove the possible secularities that might occur. The present method is seen to be fairly simple as compared to the renormalization method [Kodama, Y., & Taniuti, T. (1977). Higher order approximation in reductive perturbation method 1. Weakly dispersive system. Journal of Physics Society of Japan, 45, 298–310] and the multiple scale expansion method [Kraenkel, R. A., Manna, M. A., & Pereira, J. G. (1995). The Korteweg–deVries hierarchy and long water waves. Journal of Mathematics Physics, 36, 307–320].     

A constrained non‐linear system approach for the solution of an extended limit analysis problem

A number of recent papers (see, e.g. (Int. J. Mech. Sci. 2007; 49 :454–465; Eur. J. Mech. A/Solids 2008; 27 :859–881; Eng. Struct. 2008; 30 :664–674; Int. J. Mech. Sci. 2009; 51 :179–191)) have shown that classical limit analysis can be extended to incorporate such important features as geometric non‐linearity, softening and various so‐called ductility constraints. The generic formulation takes the form of a challenging (nonconvex and nonsmooth) optimization problem referred to in the mathematical programming literature as a mathematical program with equilibrium constraints (MPEC). Similar to a classical limit analysis, the aim is to compute in a single step a bound (upper bound, in the case of the extended problem) to the maximum load. The solution algorithm so far proposed to solve the MPEC is to convert it into an iterative non‐linear programming problem and attempts to process this using a standard non‐linear optimizer. Motivated by the fact that no method is guaranteed to solve such MPECs and by the need to avoid the use of an optimization approach, which is unfamiliar to most practising engineers, we propose, in the present paper, a novel numerical scheme to solve the MPEC as a constrained non‐linear system of equations. We illustrate the application of this approach using the simple class of elastoplastic softening skeletal structures for which certain ductility conditions are prescribed. Copyright © 2009 John Wiley & Sons, Ltd.     

Re-entrant flow shop scheduling problem with time windows using hybrid genetic algorithm based on auto-tuning strategy

The re-entrant flow shop scheduling problem considering time windows constraint is one of the most important problems in hard-disc drive (HDD) manufacturing systems. In order to maximise the system throughput, the problem of minimising the makespan with zero loss is considered. In this paper, evolutionary techniques are proposed to solve the complex re-entrant scheduling problem with time windows constraint in manufacturing HDD devices with lot size. This problem can be formulated as a deterministic Fm?|?fmls, rcrc, temp?|?C_max problem. A hybrid genetic algorithm was used for constructing chromosomes by checking and repairing time window constraints, and improving chromosomes by a left-shift heuristic as a local search algorithm. An adaptive hybrid genetic algorithm was eventually developed to solve this problem by using fuzzy logic control in order to enhance the search ability of the genetic algorithm. Finally, numerical experiments were carried out to demonstrate the efficiency of the developed approaches.     

The importance of diagonal dominance in the iterative solution of equations generated from the boundary element method

The unsymmetric matrix equations generated from the boundary element method (BEM) can be solved iteratively, with convergence to the correct solution guaranteed, if the boundary element system of equations can be first transformed into an equivalent, diagonally dominant system. A transformation is presented which selectively annihilates terms in the coefficient matrix of the system Ax = b until an equivalent, row diagonally dominant system, if available, is obtained. The new, row diagonally dominant system is well suited for use with Jacobi and Gauss-Seidel point iterative equation solvers. The diagonal dominizing transformation presented in this work is not suitable for large systems of equations but is useful as a research tool for studying the importance of diagonal dominance in the iterative solution of equations generated from the BEM. A simple Laplacian problem is used to examine the structure of the BEM equations and to introduce the diagonal dominizing transformation. The importance of diagonal dominance is shown by comparing iterative convergence of positive-definite, symmetric positive-definite and diagonally dominant systems of BEM equations obtained from a plane strain elasticity problem.     

Optical solitons for cascaded system: Jacobi elliptic functions

Abstract

In this work, the coupled NLSE is studied that appears in a cascased system. Exact solutions of this coupled system acquire by using Jacobi elliptic functions. The presented system will be studied with two forms of nonlinearity. The forms are kerr law and power law. Dark, bright, singular solitons and Jacobi elliptic solutions will be acquired. We present constraint conditions for the existence of obtained solitons. Our obtained solitons will be controlled by changing the diffraction coefficients.     

An efficient numerical method for the solution of sliding contact problems

In this paper, an efficient numerical method to solve sliding contact problems is proposed. Explicit formulae for the Gauss–Jacobi numerical integration scheme appropriate for the singular integral equations of the second kind with Cauchy kernels are derived. The resulting quadrature formulae for the integrals are valid at nodal points determined from the zeroes of a Jacobi polynomial. Gaussian quadratures obtained in this manner involve fixed nodal points and are exact for polynomials of degree 2n ? 1, where n is the number of nodes. From this Gauss–Jacobi quadrature, the existing Gauss–Chebyshev quadrature formulas can be easily derived. Another apparent advantage of this method is its ability to capture correctly the singular or regular behaviour of the tractions at the edge of the region of contact. Also, this analysis shows that once if the total normal load and the friction coefficient are given, the external moment M and contact eccentricity e (for incomplete contact) in fully sliding contact are uniquely determined. Finally, numerical solutions are computed for two typical contact cases, including sliding Hertzian contact and a sliding contact between a flat punch with rounded corners pressed against the flat surface of a semi‐infinite elastic solid. These results provide a demonstration of the validity of the proposed method. Copyright © 2005 John Wiley & Sons, Ltd.     

Synthesis and Superior Optical‐Limiting Properties of Fluorene‐Thiophene‐Benzothiadazole Polymer‐Functionalized Graphene Sheets

A polymer based on fluorene, thiophene, and benzothiadazole as the donor–spacer–acceptor triad is covalently coupled to reduced graphene oxide (rGO) sheets via diazonium coupling with phenyl bromide, followed by Suzuki coupling. These polymer–graphene hybrids show good solubility in organic solvents, such as chloroform, tetrahydrofuran (THF), toluene, dichlorobenzene, and N,N‐dimethylformamide (DMF), and exhibit an excellent optical‐limiting effect with a 532‐nm laser beam. The optical‐limiting threshold energy values (0.93 J cm^?2 for G–polymer 1 and 1.12 J cm^?2 for G–polymer 2) of these G–polymer hybrids are better than that of carbon nanotubes (3.6 J cm^?2).     

A note on “an improved MCDM DEA model for technology selection”

This paper gives a modification on the improved technology selection DEA model proposed by Amin et al. (Amin, G.R., Toloo, M. and Sohrabi, B., An improved MCDM DEA model for technology selection. Int. J. Prod. Res., 2006, 44, 2681–2686). The paper shows the problem of using the existing model and then introduces a new modified one to obtain a single efficient DMU for the technology selection alternatives.     

Hierarchical high‐order conforming C1 bases for quadrangular and triangular finite elements

We present three new sets of C¹ hierarchical high‐order tensor‐product bases for conforming finite elements. The first basis is a high‐order extension of the Bogner–Fox–Schmit basis. The edge and face functions are constructed using a combination of cubic Hermite and Jacobi polynomials with C¹ global continuity on the common edges of elements. The second basis uses the tensor product of fifth‐order Hermite polynomials and high‐order functions and achieves global C¹ continuity for meshes of quadrilaterals and C² continuity on the element vertices. The third basis for triangles is also constructed using the tensor product of one‐dimensional functions defined in barycentric coordinates. It also has global C¹ continuity on edges and C² continuity on vertices. A patch test is applied to the three considered elements. Projection and plate problems with smooth fabricated solutions are solved, and the performance of the h‐ and p‐refinements are evaluated by comparing the approximation errors in the L₂‐ and energy norms. A plate with singularity is then studied, and h‐ and p‐refinements are analysed. Finally, a transient problem with implicit time integration is considered. The results show exponential convergence rates with increasing polynomial order for the triangular and quadrilateral meshes of non‐distorted and distorted elements. Copyright © 2016 John Wiley & Sons, Ltd.     

Cuprate superconductors: Universal properties and trends; evidence for Bose-Einstein condensation?

We explore the compatibility of empirical trends in various thermodynamic properties of cuprate superconductors with the Bose-Einstein condensation scenario. These trends include the relations between transition temperature, hole concentration and condensate density, the rise and the upper limit of the transition temperature, the dependence of pressure and isotope coefficients on transition temperature, as well as the observed critical behavior, which is reminiscent of three-dimensional systems with a scalar complex order parameter and short-range interactions. For this purpose we consider an interacting charged Bose gas. Due to the high polarizability of the cuprates, the Coulomb interaction is strongly screened. For this reason, the problem of calculating thermodynamic properties becomes essentially equivalent to that of the uncharged gas with short-range interactions. This problem, however, has not been solved either. Nevertheless, in the dilute limit the problem reduces to the ideal Bose gas treated by Schafroth, while in the dense regime condensation and superfluidity are suppressed because bosons of finite extension fill the available volume. This limiting behavior provides an interpolation scheme for the dependence of both transition temperature and zero temperature superfluid density on boson density. On this basis, and relating the hole concentration in the cuprates which corresponds to the boson density and the superfluid density to the square of the inverse London penetration depth, the compatibility of the empirical trends in the cuprates with the Bose gas behavior can be verified. Our analysis reveals remarkable agreement between these trends and the corresponding Bose gas behavior. There is even strong evidence of the most striking implication of this scenario, the dependence of the transition temperature on the zero-temperature superfluid density, which resembles the outline of a fly's wing. This evidence emerges from recentSR data for Tl₂Ba₂CuO₆₊ and kinetic inductance measurements for La_2–xSr_xCuO₄ films, revealing that the penetration depths of underdoped and overdoped samples atT _c do not differ significantly. In view of this we found considerable evidence of the nature of the superconducting transition in the cuprates, without invoking any specific pairing mechanism.The authors are grateful to J. G. Bednorz, D. Baeriswyl, H. Beck, J. I. Budnick, H. Keller, K. A. Müller, Ch. Niedermeier, and J. J. Rodriguez for valuable discussions.     

Study of the simple extension tear test sample for rubber with Configurational Mechanics

The simple extension tear test-piece also referred to as the trousers sample is widely used to study crack propagation in rubber. The corresponding energy release rate, called tearing energy for rubber materials, was first established by Rivlin and Thomas (J Polym Sci, 10:291–318, 1953); a second derivation was proposed later by Eshelby (In G.C. Sih, H. C. van Elst, and D. Broek, editors, Prospects of Fracture Mechanics, 69-84, Leyden, 1975). We show here that the derivation of this result can be advantageously revisited through the scope of Configurational Mechanics. Our approach is based on the rigorous definition of the configurations of the body and on the physical significance of the configurational stress tensor. More precisely, it is demonstrated that the change in energy due to crack growth, and then the tearing energy, is directly related to the components of the configurational stress tensor in the body.     

Fuzzy adaptive preventive maintenance in a manufacturing control system: a step towards self-maintenance

In this paper, a framework of integrating preventive maintenance (PM) and manufacturing control system is proposed. Fuzzy-logic control is used to enable an intelligent approach of integrating PM and a manufacturing control system. This will contribute to the novel development of an integrated and intelligent framework in those two fields that are sometimes difficult to achieve. This idea is based on combining work on an intelligent real-time controller for a failure-prone manufacturing system using a fuzzy-logic approach (Yuniarto, M.N. and Labib, A.W., Optimal control system of an unreliable machine using fuzzy logic control: from design to implementation. Int. J. Prod. Res. (in press a); Yuniarto, M.N. and Labib, A.W., Intelligent real time control of disturbances in manufacturing systems. Integr. Manuf. Syst.: Int. J. Manuf. Technol. Manage. (in press b) and the work on PM proposed by Labib et al. (Labib, A.W., Williams, G.B. and O’Connor, R.F., An intelligent maintenance model (system): an application of analytic hierarchy process and a fuzzy logic rule-based controller. J. Oper. Res. Soc., 1998, 49, 745–757)). The aim of the research is to control a failure-prone manufacturing system and at the same time propose which PM method is applicable to a specific failure-prone manufacturing system. The mean time to repair and mean time between failures of the system are used as integrator agents, by using them to couple the two areas to be integrated (i.e. a maintenance system and manufacturing system).     

On an SSOR matrix relationship and its consequences

A matrix relationship connecting the Jacobi and the Symmetric Successive Overrelaxation (SSOR) matrices associated with a k-cyclic consistently ordered matrix A is presented. Next the equivalence of the SSOR method and a certain monoparametric k-step one for the solution of the linear algebraic system Ax = b is established. The aforementioned equivalence can be exploited to derive regions of convergence, optimum parameters involved, etc. of the two iterative methods above. This is done by studying the simplest of the two methods that is the monoparametric k-step one. To show how the idea works the case k = 2 is very briefly discussed.