A low cost numerical tool for the calculation of mode I stress intensity factors K in embedded and surface irregular cracks is presented in this paper. The proposed tool is an extension of the O-integral algorithm due...
详细信息
A low cost numerical tool for the calculation of mode I stress intensity factors K in embedded and surface irregular cracks is presented in this paper. The proposed tool is an extension of the O-integral algorithm due to Oore and Burns for the assessment of embedded plane cracks using the weight function methodology. The performance of the O-integral is assessed first by comparing its K results to exact solutions for embedded elliptical and rectangular cracks. From the analysis of this data it is found that the error in the K results systematically depends on the crack aspect ratio and the local crack front curvature. Based on this evidence a corrective function is derived in order to remediate the limitations of the O-integral. Solutions due to Newman and Raju are used to account for the effects of free surfaces and finite thickness. The accuracy of the proposed procedure is assessed by solving a number of examples and by comparing the obtained results to those available in the literature. (c) 2006 Elsevier Ltd. All rights reserved.
This paper presents a new numerical spectral domain algorithm devoted to blocking unsuccessful automatic reclosing onto permanent faults and fault distance calculation. Are voltage amplitude and fault distance are cal...
详细信息
This paper presents a new numerical spectral domain algorithm devoted to blocking unsuccessful automatic reclosing onto permanent faults and fault distance calculation. Are voltage amplitude and fault distance are calculated from the fundamental and third harmonies of the terminal voltages and currents phasors. The electric arc is modeled with its voltage waveshape defined numerically on the basis of a number of arc voltage records obtained in the high-voltage laboratory. From the calculated arc voltage amplitude, it can be concluded whether the fault is a transient arcing fault or a permanent arcless fault. If the fault is permanent, automatic reclosure should be blocked. The algorithm can be applied for adaptive autoreclosure, distance protection, and fault location. The results of algorithm testing through computer simulation and real field record are given.
A continuum model is presented for the motion of a domain wall in a plane 90 degrees-domain configuration subjected to an isolated extrinsic charge near the surface of a ferroelectric single crystal. Local pinning is ...
详细信息
A continuum model is presented for the motion of a domain wall in a plane 90 degrees-domain configuration subjected to an isolated extrinsic charge near the surface of a ferroelectric single crystal. Local pinning is postulated for the kinetic law. Before the appearance of the extrinsic charge, all polarization surface charges are taken to be neutralized by environmental charges. The domain wall motion after the appearance of the extrinsic charge is assumed to proceed sufficiently fast without any significant conductive currents on the surface or in the interior of the crystal such that new surface and interface polarization charges remain unscreened and contribute to the ferroelectric anisotropy energy. A non-admissible divergence of the electric field and consequently of the local thermodynamic driving force and of the domain wall velocity appears in the model if the domain wall charged by interface polarization charges intersects the crystal surface charged by surface polarization charges under an arbitrary angle. The physically possible domain wall angle is identified using the condition of a non-divergent driving force. The ferroelectric anisotropy energy and an intrinsic surface energy of the domain wall, however, do not provide stability of the domain wall trajectory against an unlimited increase of its curvature at the surface. The problem has been solved conceptually by proper account of the domain wall bending energy. numerical and dimensional analysis explain also why domain walls driven by extrinsic charges remain almost straight in soft ferroelectrics. (c) 2005 Elsevier Ltd. All rights reserved.
The generalized Mittag-Leffler function E-alpha,E-beta (z) has been studied for arbitrary complex argument z is an element of C and parameters alpha is an element of R+ and beta is an element of R. This function plays...
详细信息
The generalized Mittag-Leffler function E-alpha,E-beta (z) has been studied for arbitrary complex argument z is an element of C and parameters alpha is an element of R+ and beta is an element of R. This function plays a fundamental role in the theory of fractional differential equations and numerous applications in physics. The Mittag-Leffler function interpolates smoothly between exponential and algebraic functional behaviour. A numerical algorithm for its evaluation has been developed. The algorithm is based on integral representations and exponential asymptotics. Results of extensive numerical calculations for E-alpha,E-beta (z) in the complex z -plane are reported here. We find that all complex zeros emerge from the point z =1 for small alpha. They diverge towards -infinity+(2 k -1)pi i for alpha -> 1(-) and towards -infinity+2 k pi i for alpha -> 1(+) ( k is an element of Z). All the complex zeros collapse pairwise onto the negative real axis for alpha -> 2. We introduce and study also the inverse generalized Mittag-Leffler function L-alpha,L-beta (z) defined as the solution of the equation L-alpha,L-beta (E-alpha,E-beta (z )) = z . We determine its principal branch numerically.
This Note presents a time-integration strategy for computing the evolution of structures embedding shape-memory alloys in a thermomechanical setting. A variational formulation is associated with the scheme proposed, w...
详细信息
This Note presents a time-integration strategy for computing the evolution of structures embedding shape-memory alloys in a thermomechanical setting. A variational formulation is associated with the scheme proposed, which allows one to study the existence and unicity of solutions depending on the material model considered. A numerical example is presented to illustrate the method and discuss the influence of the thermomechanical coupling.
The aim of this paper is the analysis of the bias contained in star catalogues, these errors can be obtained by means of kinematics or dynamical methods. In this paper a general method, suitable when the sample is not...
详细信息
The aim of this paper is the analysis of the bias contained in star catalogues, these errors can be obtained by means of kinematics or dynamical methods. In this paper a general method, suitable when the sample is not homogeneous is presented.
A new algorithm is proposed for computing the intersection of two plane curves given in rational parametric form. It relies on the Ehrlich-Aberth iteration complemented with some computational tools like the propertie...
详细信息
A new algorithm is proposed for computing the intersection of two plane curves given in rational parametric form. It relies on the Ehrlich-Aberth iteration complemented with some computational tools like the properties of Sylvester and Bezout matrices, a stopping criterion based on the concept of pseudo-zero, an inclusion result and the choice of initial approximations based on the Newton polygon. The algorithm is implemented as a Fortran 95 module. From the numerical experiments performed with a wide set of test problems it shows a better robustness and stability with respect to the Manocha-Demmel approach based on eigenvalue computation. In fact, the algorithm provides better approximations in terms of the relative error and performs successfully in many critical cases where the eigenvalue computation fails.
Two algorithms for robust pole assignment by state feedback, proposed by Kautsky, Nichols and Van Dooren (1985) and by Tits and Yang (1996) are briefly reviewed. MATLAB code implementations of these algorithms, place ...
详细信息
ISBN:
(纸本)9780780397972
Two algorithms for robust pole assignment by state feedback, proposed by Kautsky, Nichols and Van Dooren (1985) and by Tits and Yang (1996) are briefly reviewed. MATLAB code implementations of these algorithms, place (from the MATLAB Control System Toolbox) and robpole (from SLICOT), are then numerically compared on randomly generated test data sets, as well as on examples from two benchmark collections, in terms of the robustness (insensitivity of poles to variations in plant parameters) of the closed-loop systems they produce. The functions place and robpole are also compared with each other, as well as with the (non robust) pole assignment code pass (from SLICOT) in terms of CPU time and accuracy of the pole assignment.
Different extrapolation strategies are presented for the construction of starting values in implicit Newmark time-stepping schemes used for the solution of the electromagnetic wave equation. The implicit character of ...
详细信息
Different extrapolation strategies are presented for the construction of starting values in implicit Newmark time-stepping schemes used for the solution of the electromagnetic wave equation. The implicit character of this schemes implies the solution of a system of linear equations in each time step. The choice of a good start value derived from the extrapolation of the equation system is shown to improve the convergence of the iterative conjugate gradient solver which reduces the total simulation time.
A health monitoring system based on analytical redundancy is developed for satellites on elliptical orbits. Analytical redundancy, which reduces the need for hardware redundancy, uses the modeled dynamic relationship ...
详细信息
A health monitoring system based on analytical redundancy is developed for satellites on elliptical orbits. Analytical redundancy, which reduces the need for hardware redundancy, uses the modeled dynamic relationship between system inputs and measured system outputs to form a residual process that is used for detecting and identifying faults. First, the dynamics of the satellite including orbital mechanics and attitude dynamics is modeled as a periodic system. Then, periodic fault detection filters are designed to detect and identify the satellite's actuator and sensor faults. In addition, parity equations are constructed using the algebraic redundant relationship among actuators and sensors. Furthermore, a residual processor is designed to generate the probability of each fault by using a sequential probability test. Finally, the health monitoring system, consisting of periodic fault detection filters, parity equations, and residual processor, is evaluated in the simulation in the presence of disturbances and uncertainty.
暂无评论