In this paper a refined higher-order shear deformation theory for the instability finite element analysis of fibre reinforced shell like composite structures is developed. A higher order shear deformation theory allow...
详细信息
In this paper a refined higher-order shear deformation theory for the instability finite element analysis of fibre reinforced shell like composite structures is developed. A higher order shear deformation theory allows parabolic description of the transverse shear stresses and therefore the shear correction factors of the usual shear deformation theory are not required. The present formulation is based on a higher order shear theory in which in-plane displacements are expanded as cubic functions of the thickness coordinate. The conditions of zero transverse shear stresses on the top and bottom faces are satisfied. Laminate material is assumed to be linearly elastic, homogeneous and isotropic/orthotropic. The 4-node quadrilateral shell finite element with 8 degrees of freedom per node has been developed which eleviates most of the deficiencies associated with such types elements. The effects of the transverse shear deformation on the buckling loads are investigated. It is shown that the present theory predicts buckling loads more accurately then CTP or the first order shear deformation theory. To determine buckling loads here is used lanczos algorithm. The good agreement of the numerical and experimental results is indicative of a reliable presented shell element and procedure for practical instability analysis of the structures made of the fiber reinforced laminates.
An efficient algorithm is presented for the finite-element solution of electromagnetic eigenvalue problems associated with lossy and unbounded structures. The algorithm is based on the (E) over right arrow - (B) over ...
详细信息
An efficient algorithm is presented for the finite-element solution of electromagnetic eigenvalue problems associated with lossy and unbounded structures. The algorithm is based on the (E) over right arrow - (B) over right arrow formulation of the finite-element approximation of the electromagnetic equations. The special relationship between the vector bases used for the expansion of the electric field vector (E) over right arrow and the magnetic flux density vector (B) over right arrow is used to reduce the computational complexity of the proposed formulation. The occurrence of spurious de modes is avoided through the careful selection of divergence-free initial vectors in the Lanezos-Arnoldi-based iterative schemes used by the proposed algorithm. The resultant algorithm is only marginally more expensive than standard finite-element-based algorithms used for electromagnetic eigenvalue problems involving lossless structures. Numerical experiments from the application of the proposed algorithm to the eigenvalue analysis of both lossless and lossy cavities are used to demonstrate its accuracy, computational efficiency, and robustness.
The smoothed finite element methods (S-FEMs) have been found capable of producing softer models whose "stiffness" is closer to the true continuous model. This paper explores, for the first time, this unique ...
详细信息
The smoothed finite element methods (S-FEMs) have been found capable of producing softer models whose "stiffness" is closer to the true continuous model. This paper explores, for the first time, this unique feature of SFEMs to develop a complete formulism and procedure to study free vibration and forced vibration of solid structures, via (1) solving eigenvalue problems that produces vibration modes of a given structure;(2) using model superimposition techniques and the lanczos algorithm to obtain transient dynamic solution for structures subjected to arbitrary dynamics forces. The present S-FEM modeling takes the advantageous of the so-called softening effects achieved by establishing proper types of smoothing domains based on edges and faces of the mesh known as ES/FS-FEM, so as to obtain accurate modes for both free and forced vibration analysis. For easy automation in creating 3D solids, we use only the automatically generatable tetrahedral mesh, while to ensure excellent stress solution using the ES/FS-FEM models. A 3D code has been developed in the framework of SFEMs, and applied to solve a number of 3D solid structures. The results are compared with those from the commercial finite element analysis software ABAQUS in terms of accuracy and convergence.
Bounds for entries of matrix functions based on Gauss-type quadrature rules are applied to adjacency matrices associated with graphs. This technique allows to develop inexpensive and accurate upper and lower bounds fo...
详细信息
Bounds for entries of matrix functions based on Gauss-type quadrature rules are applied to adjacency matrices associated with graphs. This technique allows to develop inexpensive and accurate upper and lower bounds for certain quantities (Estrada index, subgraph centrality, communicability) that describe properties of networks. (C) 2010 Elsevier Inc. All rights reserved.
When using edge element basis functions for the solution of eigenmodes of the vector wave equation, "dc spurious modes" are introduced. The eigenvalues of these modes are zero and their corresponding eigenve...
详细信息
When using edge element basis functions for the solution of eigenmodes of the vector wave equation, "dc spurious modes" are introduced. The eigenvalues of these modes are zero and their corresponding eigenvectors are in the null space of the curl operator. These modes arise due to the irrotational vector space spanned by the edge element basis functions and lead to nonzero divergence of the electric flux. We introduce a novel method to eliminate the occurrence of such solutions using "divergence-free" constraint equations. The constraint equations are imposed efficiently by tree-cotree partitioning of the finite-element mesh and does not require any basis functions other than the edge elements. The constraint equations can be directly incorporated into any Krylov-subspace-based eigenvalue solver, such as the Lanezos/Arnoldi algorithm used widely for the solution of generalized sparse eigenvalue problems.
We construct a microscopic theory for the proton spin-lattice relaxation-rate 1/T-1 measurements around field-induced level crossings in a single crystal of the trivalent chromium ion wheel complex [Cr8F8(OOCtBu)(16)]...
详细信息
We construct a microscopic theory for the proton spin-lattice relaxation-rate 1/T-1 measurements around field-induced level crossings in a single crystal of the trivalent chromium ion wheel complex [Cr8F8(OOCtBu)(16)] at sufficiently low temperatures [E. Micotti et al., Phys. Rev. B 72 (2005) 020405(R)]. Exactly diagonalizing a well-equipped spin Hamiltonian for the individual clusters and giving further consideration to their possible interactions, we reveal the mechanism of 1/T-1] being single-peaked normally at the first level crossing but double-peaked intriguingly around the second level crossing. We wipe out the doubt about poor crystallization and find out a solution intramolecular alternating Dzyaloshinsky-Moriya interaction combined with intermolecular coupling of antiferromagnetic character, each of which is so weak as several tens of mK in magnitude. (C) 2015 Elsevier B.V. All rights reserved.
The full-wave analysis of the generalized microstrip line on an inhomogeneous anisotropic substrate is carried out by using the finite-difference method. The resulting sparse matrix equation is solved efficiently usin...
详细信息
The full-wave analysis of the generalized microstrip line on an inhomogeneous anisotropic substrate is carried out by using the finite-difference method. The resulting sparse matrix equation is solved efficiently using the bi-lanczos algorithm, The use of the inhomogeneous wave equation to formulate the problem makes it easy to analyze structures with multilayered substrates, The algorithm can analyze complicated structures with multiple conductors at arbitrary locations. A spatial interpolation scheme is used to evaluate the contribution from the off-diagonal terms in <(mu)over bar> and <(epsilon)over bar> The use of the bi-lanczos algorithm allows us to solve the problem at O(N-1.5) complexity. Storage requirements can be made to scale as O(N), This makes it possible to analyze large problems on a small computer. Very good agreement is seen between published results and results obtained using this technique.
The parallel computation of the pseudospectrum is presented. The Parallel Path following algorithm using Triangles (PPAT) is based on the Path following algorithm using Triangles (PAT). This algorithm offers total rel...
详细信息
The parallel computation of the pseudospectrum is presented. The Parallel Path following algorithm using Triangles (PPAT) is based on the Path following algorithm using Triangles (PAT). This algorithm offers total reliability and can handle singular points along the level curve without difficulty. Furthermore, PPAT offers a guarantee of termination even in the presence of round-off errors and makes use of the large granularity for parallelism in PAT. This results in large speedups and high efficiency. The PPAT is able to trace multiple level curves simultaneously and takes into account the symmetry of the pseudospectrum in the case of real matrices to reduce the total computational cost. (C) 2002 Elsevier Science B.V. All rights reserved.
Let A is an element of C-NxN be a nonsingular complex matrix and b and c be complex vectors of length N. The goal of this paper is to investigate approaches for efficient approximations of the bilinear form c* A(-1)b....
详细信息
Let A is an element of C-NxN be a nonsingular complex matrix and b and c be complex vectors of length N. The goal of this paper is to investigate approaches for efficient approximations of the bilinear form c* A(-1)b. Equivalently, we wish to approximate the scalar value c* x, where x solves the linear system Ax = b. Here the matrix A can be very large or its elements can be too costly to compute so that A is not explicitly available and it is used only in the form of the matrix-vector product. Therefore a direct method is not an option. For A Hermitian positive definite, b* A(-1)b can be efficiently approximated as a by-product of the conjugate-gradient iterations, which is mathematically equivalent to the matching moment approximations computed via the Gauss-Christoffel quadrature. In this paper we propose a new method using the biconjugate gradient iterations which is applicable to the general complex case. The proposed approach will be compared with existing ones using analytic arguments and numerical experiments.
Given a Hamiltonian matrix H = JS with S symmetric and positive definite, we analyze a symplectic lanczos algorithm to transform -H-2 in a symmetric and positive definite tridiagonal matrix of half size. By means of t...
详细信息
Given a Hamiltonian matrix H = JS with S symmetric and positive definite, we analyze a symplectic lanczos algorithm to transform -H-2 in a symmetric and positive definite tridiagonal matrix of half size. By means of two effective restarted procedures, this algorithm is then used to compute few extreme eigenvalues of H. Numerical examples are also reported to compare the presented techniques. (c) 2004 Elsevier B.V. All rights reserved.
暂无评论