963 resultados para Adaptive Finite Element Methods
Resumo:
This work presents the application of a scalar finite element formulation for Ex (TE-like) modes in anisotropic planar and channel waveguides with diagonal permittivity tensor, diffused in both transversal directions. This extended formulation considers explicitly both the variations of the refractive index and their spatial derivates inside of each finite element. Dispersion curves for Ex modes in planar and channel waveguides are shown, and the results compared with solutions obtained by other formulations.
Resumo:
The main feature of partition of unity methods such as the generalized or extended finite element method is their ability of utilizing a priori knowledge about the solution of a problem in the form of enrichment functions. However, analytical derivation of enrichment functions with good approximation properties is mostly limited to two-dimensional linear problems. This paper presents a procedure to numerically generate proper enrichment functions for three-dimensional problems with confined plasticity where plastic evolution is gradual. This procedure involves the solution of boundary value problems around local regions exhibiting nonlinear behavior and the enrichment of the global solution space with the local solutions through the partition of unity method framework. This approach can produce accurate nonlinear solutions with a reduced computational cost compared to standard finite element methods since computationally intensive nonlinear iterations can be performed on coarse global meshes after the creation of enrichment functions properly describing localized nonlinear behavior. Several three-dimensional nonlinear problems based on the rate-independent J (2) plasticity theory with isotropic hardening are solved using the proposed procedure to demonstrate its robustness, accuracy and computational efficiency.
Resumo:
We consider a model eigenvalue problem (EVP) in 1D, with periodic or semi–periodic boundary conditions (BCs). The discretization of this type of EVP by consistent mass finite element methods (FEMs) leads to the generalized matrix EVP Kc = λ M c, where K and M are real, symmetric matrices, with a certain (skew–)circulant structure. In this paper we fix our attention to the use of a quadratic FE–mesh. Explicit expressions for the eigenvalues of the resulting algebraic EVP are established. This leads to an explicit form for the approximation error in terms of the mesh parameter, which confirms the theoretical error estimates, obtained in [2].
Resumo:
Table of Contents
1 | Introduction | 1 |
1.1 | What is an Adiabatic Shear Band? | 1 |
1.2 | The Importance of Adiabatic Shear Bands | 6 |
1.3 | Where Adiabatic Shear Bands Occur | 10 |
1.4 | Historical Aspects of Shear Bands | 11 |
1.5 | Adiabatic Shear Bands and Fracture Maps | 14 |
1.6 | Scope of the Book | 20 |
2 | Characteristic Aspects of Adiabatic Shear Bands | 24 |
2.1 | General Features | 24 |
2.2 | Deformed Bands | 27 |
2.3 | Transformed Bands | 28 |
2.4 | Variables Relevant to Adiabatic Shear Banding | 35 |
2.5 | Adiabatic Shear Bands in Non-Metals | 44 |
3 | Fracture and Damage Related to Adiabatic Shear Bands | 54 |
3.1 | Adiabatic Shear Band Induced Fracture | 54 |
3.2 | Microscopic Damage in Adiabatic Shear Bands | 57 |
3.3 | Metallurgical Implications | 69 |
3.4 | Effects of Stress State | 73 |
4 | Testing Methods | 76 |
4.1 | General Requirements and Remarks | 76 |
4.2 | Dynamic Torsion Tests | 80 |
4.3 | Dynamic Compression Tests | 91 |
4.4 | Contained Cylinder Tests | 95 |
4.5 | Transient Measurements | 98 |
5 | Constitutive Equations | 104 |
5.1 | Effect of Strain Rate on Stress-Strain Behaviour | 104 |
5.2 | Strain-Rate History Effects | 110 |
5.3 | Effect of Temperature on Stress-Strain Behaviour | 114 |
5.4 | Constitutive Equations for Non-Metals | 124 |
6 | Occurrence of Adiabatic Shear Bands | 125 |
6.1 | Empirical Criteria | 125 |
6.2 | One-Dimensional Equations and Linear Instability Analysis | 134 |
6.3 | Localization Analysis | 140 |
6.4 | Experimental Verification | 146 |
7 | Formation and Evolution of Shear Bands | 155 |
7.1 | Post-Instability Phenomena | 156 |
7.2 | Scaling and Approximations | 162 |
7.3 | Wave Trapping and Viscous Dissipation | 167 |
7.4 | The Intermediate Stage and the Formation of Adiabatic Shear Bands | 171 |
7.5 | Late Stage Behaviour and Post-Mortem Morphology | 179 |
7.6 | Adiabatic Shear Bands in Multi-Dimensional Stress States | 187 |
8 | Numerical Studies of Adiabatic Shear Bands | 194 |
8.1 | Objects, Problems and Techniques Involved in Numerical Simulations | 194 |
8.2 | One-Dimensional Simulation of Adiabatic Shear Banding | 199 |
8.3 | Simulation with Adaptive Finite Element Methods | 213 |
8.4 | Adiabatic Shear Bands in the Plane Strain Stress State | 218 |
9 | Selected Topics in Impact Dynamics | 229 |
9.1 | Planar Impact | 230 |
9.2 | Fragmentation | 237 |
9.3 | Penetration | 244 |
9.4 | Erosion | 255 |
9.5 | Ignition of Explosives | 261 |
9.6 | Explosive Welding | 268 |
10 | Selected Topics in Metalworking | 273 |
10.1 | Classification of Processes | 273 |
10.2 | Upsetting | 276 |
10.3 | Metalcutting | 286 |
10.4 | Blanking | 293 |
Appendices | 297 | |
A | Quick Reference | 298 |
B | Specific Heat and Thermal Conductivity | 301 |
C | Thermal Softening and Related Temperature Dependence | 312 |
D | Materials Showing Adiabatic Shear Bands | 335 |
E | Specification of Selected Materials Showing Adiabatic Shear Bands | 341 |
F | Conversion Factors | 357 |
References | 358 | |
Author Index | 369 | |
Subject Index | 375 |
Resumo:
[EN]Longest edge (nested) algorithms for triangulation refinement in two dimensions are able to produce hierarchies of quality and nested irregular triangulations as needed both for adaptive finite element methods and for multigrid methods. They can be formulated in terms of the longest edge propagation path (Lepp) and terminal edge concepts, to refine the target triangles and some related neighbors. We discuss a parallel multithread algorithm, where every thread is in charge of refining a triangle t and its associated Lepp neighbors. The thread manages a changing Lepp(t) (ordered set of increasing triangles) both to find a last longest (terminal) edge and to refine the pair of triangles sharing this edge...
Resumo:
This work sets forth a `hybrid' discretization scheme utilizing bivariate simplex splines as kernels in a polynomial reproducing scheme constructed over a conventional Finite Element Method (FEM)-like domain discretization based on Delaunay triangulation. Careful construction of the simplex spline knotset ensures the success of the polynomial reproduction procedure at all points in the domain of interest, a significant advancement over its precursor, the DMS-FEM. The shape functions in the proposed method inherit the global continuity (Cp-1) and local supports of the simplex splines of degree p. In the proposed scheme, the triangles comprising the domain discretization also serve as background cells for numerical integration which here are near-aligned to the supports of the shape functions (and their intersections), thus considerably ameliorating an oft-cited source of inaccuracy in the numerical integration of mesh-free (MF) schemes. Numerical experiments show the proposed method requires lower order quadrature rules for accurate evaluation of integrals in the Galerkin weak form. Numerical demonstrations of optimal convergence rates for a few test cases are given and the method is also implemented to compute crack-tip fields in a gradient-enhanced elasticity model.
Resumo:
A residual based a posteriori error estimator is derived for a quadratic finite element method (FEM) for the elliptic obstacle problem. The error estimator involves various residuals consisting of the data of the problem, discrete solution and a Lagrange multiplier related to the obstacle constraint. The choice of the discrete Lagrange multiplier yields an error estimator that is comparable with the error estimator in the case of linear FEM. Further, an a priori error estimate is derived to show that the discrete Lagrange multiplier converges at the same rate as that of the discrete solution of the obstacle problem. The numerical experiments of adaptive FEM show optimal order convergence. This demonstrates that the quadratic FEM for obstacle problem exhibits optimal performance.
Resumo:
We have successfully extended our implicit hybrid finite element/volume (FE/FV) solver to flows involving two immiscible fluids. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centered FV method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. This updating strategy can be rigorously proven to be able to eliminate the unphysical pressure boundary layer and is crucial for the correct temporal convergence rate. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centers and the auxiliary variable at vertices. The fluid interface is captured by solving an advection equation for the volume fraction of one of the fluids. The same matrix-free FV method, as the one used for momentum equations, is used to solve the advection equation. We will focus on the interface sharpening strategy to minimize the smearing of the interface over time. We have developed and implemented a global mass conservation algorithm that enforces the conservation of the mass for each fluid.
Resumo:
A novel three-dimensional finite volume (FV) procedure is described in detail for the analysis of geometrically nonlinear problems. The FV procedure is compared with the conventional finite element (FE) Galerkin approach. FV can be considered to be a particular case of the weighted residual method with a unit weighting function, where in the FE Galerkin method we use the shape function as weighting function. A Fortran code has been developed based on the finite volume cell vertex formulation. The formulation is tested on a number of geometrically nonlinear problems. In comparison with FE, the results reveal that FV can reach the FE results in a higher mesh density.
Resumo:
Stability charts for soil slopes, first produced in the first half of the twentieth century, continue to be used extensively as design tools, and draw the attention of many investigators. This paper uses finite-element upper and lower bound limit analysis to assess the short-term stability of slopes in which the slopematerial and subgrade foundation material have two distinctly different undrained strengths. The stability charts are proposed, and the exact theoretical solutions are bracketed to within 4.2% or better. In addition, results from the limit-equilibrium method (LEM) have been used for comparison. Differences of up to 20% were found between the numerical limit analysis and LEM solutions. It also shown that the LEM sometimes leads to errors, although it is widely used in practice for slope stability assessments.
Resumo:
The authors describe a literature revision on assessing stresses in buccomaxillary prostheses photoelasticity, finite element technique, and extensometry. They describe the techniques and the importance for use of each method in buccomaxillary prostheses with implants and the need of accomplishing more studies in this scarce literary area.