Search results

This study presents a novel hp-version adaptive finite element method (FEM) to investigate the high-precision eigensolutions of the free vibration of moderately thick circular cylindrical shells, involving the issues of variable geometrical factors, such as the thickness, circumferential wave number, radius and length.

An hp-version adaptive finite element (FE) algorithm is proposed for determining the eigensolutions of the free vibration of moderately thick circular cylindrical shells via error homogenisation and higher-order interpolation. This algorithm first develops the established h-version mesh refinement method for detecting the non-uniform distributed optimised meshes, where the error estimation and element subdivision approaches based on the superconvergent patch recovery displacement method are introduced to obtain high-precision solutions. The errors in the vibration mode solutions in the global space domain are homogenised and approximately the same. Subsequently, on the refined meshes, the algorithm uses higher-order shape functions for the interpolation of trial displacement functions to reduce the errors quickly, until the solution meets a pre-specified error tolerance condition. In this algorithm, the non-uniform mesh generation and higher-order interpolation of shape functions are suitable for addressing the problem of complex frequencies and modes caused by variable structural geometries.

Numerical results are presented for moderately thick circular cylindrical shells with different geometrical factors (circumferential wave number, thickness-to-radius ratio, thickness-to-length ratio) to demonstrate the effectiveness, accuracy and reliability of the proposed method. The hp-version refinement uses fewer optimised meshes than h-version mesh refinement, and only one-step interpolation of the higher-order shape function yields the eigensolutions satisfying the accuracy requirement.

The proposed combination of methodologies provides a complete hp-version adaptive FEM for analysing the free vibration of moderately thick circular cylindrical shells. This algorithm can be extended to general eigenproblems and geometric forms of structures to solve for the frequency and mode quickly and efficiently.

The purpose of this paper is to introduce a general equation for eigensolution. Eigenvalues and eigenvectors of graphs have many applications in combinatorial optimization and structural mechanics. Some important applications of graph products consist of nodal ordering and graph partitioning for structuring the structural matrices and finite element subdomaining, respectively.

In the existing methods for the eigensolution of Laplacian matrices, members have been added to the model of a graph product such that for its Laplacian matrix an algebraic relation between blocks become possible. These methods are categorized as topological approaches. Here, using concepts of linear algebra a general algebraic method is developed.

A new algebraic method is introduced for calculating the eigenvalues of Laplacian matrices in graph products.

The present method provides a simple tool for calculating the eigenvalues of the Laplacian matrices without using the configurational model and merely by using the Laplacian matrices. The developed formula for calculating the eigenvalues contains approximate terms which can be managed by the analyst.

This paper presents an efficient methodology to calculate fuzzy eigenvalues and eigenvectors of finite element structures defined by imprecise parameters. The material and geometric parameters are then described by fuzzy numbers. The proposed methodology, based on α‐cut discretization of fuzzy numbers and Taylor's expansion, determines the extreme eigensolutions for each α‐cut. The study of a finite element model and the comparison of results with a combinatorial approach, based on Zadeh's extension principle, show the efficiency of this methodology.

The evaluation of linear dynamic response analysis of large structures by vector superposition requires, in its traditional formulation, the solution of a large and expensive eigenvalue problem. A method of solution based on a Ritz transformation to a reduced system of generalized coordinates using load dependent vectors generated from the spatial distribution of the dynamic loads is shown to maintain the high expected accuracy of modern computer analysis and significantly reduces the execution time over eigensolution procedures. New computational variants to generate load dependent vectors are presented and error norms are developed to control the convergence characteristics of load dependent Ritz solutions. Numerical applications on simple structural systems are used to show the relative efficiency of the proposed solution procedures.

This paper aims to focus on the accurate analysis of the fractional heat transfer in a two-dimensional (2D) rectangular monolayer tissue with three different kinds of lateral boundary conditions and the quantitative evaluation of the degree of thermal damage and burn depth.

A symplectic method is used to analytically solve the fractional heat transfer dual equation in the frequency domain (s-domain). Explicit expressions of the dual vector can be constructed by superposing the symplectic eigensolutions. The solution procedure is rigorously rational without any trial functions. And the accurate predictions of temperature and heat flux in the time domain (t-domain) are derived through numerical inverse Laplace transform.

Comparison study shows that the maximum relative error is less than 0.16%, which verifies the accuracy and effectiveness of the proposed method. The results indicate that the model and heat source parameters have a significant effect on temperature and thermal damage. The pulse duration (Δt) of the laser heat source can effectively control the time to reach the peak temperature and the peak slope of the thermal damage curve. The burn depth is closely correlated with exposure temperature and duration. And there exists the delayed effect of fractional order on burn depth.

A symplectic approach is presented for the thermal analysis of 2D fractional heat transfer. A unified time-fractional heat transfer model is proposed to describe the anomalous thermal behavior of biological tissue. New findings might provide guidance for temperature prediction and thermal damage assessment of biological tissues during hyperthermia.

This study aims to overcome the involved challenging issues and provide high-precision eigensolutions. General eigenproblems in the system of ordinary differential equations (ODEs) serve as mathematical models for vector Sturm-Liouville (SL) and free vibration problems. High-precision eigenvalue and eigenfunction solutions are crucial bases for the reliable dynamic analysis of structures. However, solutions that meet the error tolerances specified are difficult to obtain for issues such as coefficients of variable matrices, coincident and adjacent approximate eigenvalues, continuous orders of eigenpairs and varying boundary conditions.

This study presents an h-version adaptive finite element method based on the superconvergent patch recovery displacement method for eigenproblems in system of second-order ODEs. The high-order shape function interpolation technique is further introduced to acquire superconvergent solution of eigenfunction, and superconvergent solution of eigenvalue is obtained by computing the Rayleigh quotient. Superconvergent solution of eigenfunction is used to estimate the error of finite element solution in the energy norm. The mesh is then, subdivided to generate an improved mesh, based on the error.

Representative eigenproblems examples, containing typical vector SL and free vibration of beams problems involved the aforementioned challenging issues, are selected to evaluate the accuracy and reliability of the proposed method. Non-uniform refined meshes are established to suit eigenfunctions change, and numerical solutions satisfy the pre-specified error tolerance.

The proposed combination of methodologies described in the paper, leads to a powerful h-version mesh refinement algorithm for eigenproblems in system of second-order ODEs, that can be extended to other classes of applications in damage detection of multiple cracks in structures based on the high-precision eigensolutions.

The purpose of this paper is to present a numerically adaptive finite element (FE) method for accurate, efficient and reliable eigensolutions of regular second- and fourth-order Sturm–Liouville (SL) problems with variable coefficients.

After the conventional FE solution for an eigenpair (i.e. eigenvalue and eigenfunction) of a particular order has been obtained on a given mesh, a novel strategy is introduced, in which the FE solution of the eigenproblem is equivalently viewed as the FE solution of an associated linear problem. This strategy allows the element energy projection (EEP) technique for linear problems to calculate the super-convergent FE solutions for eigenfunctions anywhere on any element. These EEP super-convergent solutions are used to estimate the FE solution errors and to guide mesh refinements, until the accuracy matches user-preset error tolerance on both eigenvalues and eigenfunctions.

Numerical results for a number of representative and challenging SL problems are presented to demonstrate the effectiveness, efficiency, accuracy and reliability of the proposed method.

The method is limited to regular SL problems, but it can also solve some singular SL problems in an indirect way.

Comprehensive utilization of the EEP technique yields a simple, efficient and reliable adaptive FE procedure that finds sufficiently fine meshes for preset error tolerances on eigenvalues and eigenfunctions to be achieved, even on problems which proved troublesome to competing methods. The method can readily be extended to vector SL problems.

The purpose of this paper is to present an efficient method for dynamic analysis of structures utilizing a modal analysis with the main purpose of decreasing the computational complexity of the problem. In traditional methods, the solution of initial-value problems (IVPs) using numerical methods like finite difference method leads to step by step and time-consuming recursive solutions.

The present method is based on converting the IVP into boundary-value problems (BVPs) and utilizing the features of the latter problems in efficient solution of the former ones. Finite difference formulation of BVPs leads to matrices with repetitive tri-diagonal and block tri-diagonal patterns wherein the eigensolution and matrix inversion are obtained using graph products rules. To get advantage of these efficient solutions for IVPs like the dynamic analysis of single DOF systems, IVPs are converted to boundary-value ones using mathematical manipulations. The obtained formulation is then generalized to the multi DOF systems by utilizing modal analysis.

Applying the method to the modal analysis leads to a simple and efficient formulation. The laborious matrix inversion and eigensolution operations, of computational complexities of O(n2.373) and O(n3), respectively, are converted to a closed-form formulation with summation operations.

No limitation.

Swift analysis has become possible.

Suitability of solving IVPs and modal analysis using conversion and graph product rules is presented and applied to efficient seismic optimal analysis and preliminary design.

Spurious modes often appear in the computed spectrum when an electromagnetic eigenproblem is solved by the finite element method. Demonstrates that the inclusion condition, often claimed as the theoretical reason for the absence of (non‐zero frequency) spurious modes, is a sufficient but not necessary condition for that. Does this by proving that edge elements, which are spectrally correct, do not satisfy the inclusion condition. As intermediate steps towards this result, proves the equivalence of the inclusion condition to a less cryptic one and gives two more easily‐checked necessary conditions for the latter. Concludes that from this investigation, the inclusion condition seems too strong to be useful as a sufficient condition. Works out the present analysis in the framework of spectral approximation theory for non‐compact operators, which emerges as a basic tool for a deeper understanding of the whole question of spurious modes.

Considering typical self‐adjoint and non‐self‐adjoint problems that are governed by differential equations with predominant lower order derivatives, a comparison is presented of their finite element solutions by Ritz, Galerkin, least square (LSQ) and discrete least square (DLSQ) methods.