We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
A numerical method is proposed for a class of one-dimensional stochastic control problems with unbounded state space. This method solves an infinite-dimensional linear program, equivalent to the original formulation based on a stochastic differential equation, using a finite element approximation. The discretization scheme itself and the necessary assumptions are discussed, and a convergence argument for the method is presented. Its performance is illustrated by examples featuring long-term average and infinite horizon discounted costs, and additional optimization constraints.
Large truss structures have many potential applications in space, such as antennas, telescopes and space solar power plants. In this scenario, a natural concern is the susceptibility of these lightweight structures to be damaged during their operational life, due to impacts, transient thermal states and fatigue phenomena. The inclusion of active elements, equipped with sensor/actuator systems capable of modulating their shape and strength, makes it possible to transform the truss into a smart structure capable of remedying the damage, once it is detected. In this paper, a procedure is described that is capable of restoring at least the basic functionality of a composite truss for space applications, starting with the observation that damage has occurred, regardless of its specific location. The system eigenstructure is used as a benchmark for damage detection, as well as a target characteristic for the subsequent restoration activity. The observer/Kalman filter identification algorithm (OKID), in cascade with the eigensystem realization algorithm (ERA), is adopted to reconstruct, from sensor recordings, the dynamic response of the truss in terms of system state-space representation and eigen-characteristics. Finally, a static output feedback control is developed to recover the low-frequency dynamic behaviour of the truss. The entire procedure is tested using finite element analysis. All activities are coordinated in an innovative procedure that, within a unique Python language code, automatically generates finite element (FE) models, launches finite element analysis (FEA), extracts output data, implements OKID-ERA, processes the control law and applies it to the final FE simulation.
Until quite recently, discussions on “polycrystals” have been rather concentrated on or confined to how to realistically evaluate the averaged (macroscopic) stress-strain response, focusing on, e.g., relaxed constraint even with FEM simulations. This chapter discusses new perspectives related to Scale C and the attendant theory and modeling for polycrystalline materials including nanocrystals based on the field theory (they mostly are the latest achievements). Emphasis here is placed on the collective effects brought about by a large number of composing grains on the meso- and macroscopic deformation behavior of polycrystals, in the context of hierarchy of polycrystalline plasticity. For this purpose, a series of systematically designed finite element simulations have been conducted.
Typical inhomogeneities to be evolved in this scale level are the deformation-induced structures, normally yielding lamellar or band-like morphologies accompanied by relatively large “misorientation” across the bands or the walls. Since the “misorientation” is introduced so as for the grain of interest to accommodate the imposed geometrical constraint from its surroundings, these substructures are roughly categorized in “geometrically-necessary” types of bands (GNBs), in contrast to the dislocation cells in Scale A (being “mechanically-necessary”), which mediates the other two scales, therefore, may be expressed as “absorber.” Also, the chapter discusses the inhomogeneity evolutions in Scale B based on FE-based simulations, which incorporates the incompatibility-tensor field model in its constitutive framework. Starting from showing preliminary simulation results, some advanced outcomes are presented, including modeling of metallurgical microstructures (e.g., martensite lath block and packet) as a further extension.
The finite element method (FEM) is widely used to simulate a variety of physics phenomena. Approaches that integrate FEM with neural networks (NNs) are typically leveraged as an alternative to conducting expensive FEM simulations in order to reduce the computational cost without significantly sacrificing accuracy. However, these methods can produce biased predictions that deviate from those obtained with FEM, since these hybrid FEM-NN approaches rely on approximations trained using physically relevant quantities. In this work, an uncertainty estimation framework is introduced that leverages ensembles of Bayesian neural networks to produce diverse sets of predictions using a hybrid FEM-NN approach that approximates internal forces on a deforming solid body. The uncertainty estimator developed herein reliably infers upper bounds of bias/variance in the predictions for a wide range of interpolation and extrapolation cases using a three-element FEM-NN model of a bar undergoing plastic deformation. This proposed framework offers a powerful tool for assessing the reliability of physics-based surrogate models by establishing uncertainty estimates for predictions spanning a wide range of possible load cases.
This chapter provides a very brief summary of the types of heterogeneous materials considered in this monograph: fiber-reinforced composites, particulate composites, nanocomposites, porous composites, and so on. A succinct summary is given of analytical homogenization methods to determine the overall properties of particulate composites based on the upper and lower bounds of Hashin and Shtrikman; the Eshelby ellipsoidal inclusion theory and the Self-consistent Method of Eshelby; and the Mori-Tanaka Method and some other semi-analytical methods. Numerical methods such as the finite element method, the boundary element method, XFEM, and so on to model a representative volume element (RVE) of a heterogeneous material are reviewed, and thus the motivation for the Computational Grains method discussed in the rest of this book is presented.
This chapter discusses the role of a representative volume element (RVE) in the computational homogenization of heterogeneous materials. The use of the finite element method in modeling an RVE is discussed. The role of using the Hill-Mandel boundary conditions, and the use of periodic displacement and aperiodic traction boundary conditions on an RVE are discussed. The advantages of using the present Computational Grains method in modeling an RVE, not only to determine the macro properties of a heterogeneous material but also to determine the detailed interfacial stress states which are damage precursors at the micro level are discussed.
This is the first book that systemically introduces the theory and implementation of Computational Grains for micromechanical modeling of heterogeneous materials. This book covers the specifically designed mathematics embedded in Computational Grains, and the entire process of microstructure construction, tessellation, CG simulation and homogenization. The Computational Grains discussed in this book consider elastic, non-elastic, and multi-physics solids. Materials damage development are also preliminarily discussed, with CGs considering matrix-inclusion debonding as well as embedded microcracks. Presenting the theory step-by-step and with detailed examples and MATLAB codes, the material is accessible and practical for readers. This will be ideal for graduate students and researchers in mechanical and aerospace engineering and applied mechanics.
In Chapter 4, firstly a few basic terms (object and configuration, stress, strain, and constitutive relation between stress tensor and strain tensor), three coordinate systems (shape coordinate, lattice coordinate, and laboratory coordinate), deformation gradient as well as fundamental equations in continuum mechanics are briefly recalled for the sake of understanding fundamental equations of the crystal plasticity finite element method (CPFEM). A few advantages of CPFEM (including its abilities to analyze multiparticle problems and solve crystal mechanics problems with complex boundary conditions) are highlighted. Then, representative mechanical constitutive laws of crystal plasticity including dislocation-based constitutive models and constitutive models for displacive transformation are briefly described, followed by a short introduction to the finite element method (FEM), several FEM software packages (including Adina, ABAQUS, Deform, and ANSYS) and a procedure for CPFEM simulation. Finally, a case study of plastic deformation-induced surface roughening in Al polycrystals is demonstrated to show important features of crystal plasticity finite element method in materials design.
Chapter 7 briefly introduces steels, including classification, production processes, microstructure, and properties as well as computational tools for design of steels. Two case studies for S53 and AISI H13 steels are demonstrated. For S53 steel, high strength and good corrosion resistance are needed. For that purpose, plots of thermodynamic driving forces for precipitates were established, guaranteeing the accurate precipitation of M2C strengthener in steels. In addition, a martensite model is developed, designing maximal strengthening effect and appropriate martensite start temperature to maintain an alloy with lath martensite as the matrix. The corrosion resistance was designed by analyzing thermodynamic effects to maximize Cr partitioning in spinel oxide and enhance the grain boundary cohesion. In the case of AISI H13 steel, precipitations of carbides were simulated. Then simulated microstructure was coupled with structure–property models to predict the stress–strain curve and creep properties. Subsequently, those simulated properties were coupled with FEM to predict the relaxation of internal stresses and deformation behavior at the macroscopic scale during tempering of AISI H13
Predicting the onset of shear localization is among the most challenging problems in machining. This phenomenon affects the process outputs, such as machining forces, surface quality, and machined part tolerances. To predict this phenomenon, analytical, experimental, and numerical methods (especially finite element analysis) are widely used. However, the limitations of each method hinder their industrial applications, demanding a reliable and time-saving approach to predict shear localization onset. Additionally, since this phenomenon largely depends on the type and parameters of the constitutive material model, any change in these parameters requires a new set of simulations, which puts further restrictions on the application of finite element modeling. This study aims to overcome the computational efficiency of the finite element method to predict the onset of shear localization when machining Ti6Al4V using machine learning methods. The obtained results demonstrate that the FCM (fuzzy c-means) clustering ANFIS (adaptive network-based fuzzy inference system) has given better results in both training and testing when it is compared to the ANN (artificial neural network) architecture with an R2 of 0.9981. Regarding this, the FCM-ANFIS is a good candidate to calculate the critical cutting speed. To the best of the authors’ knowledge, this is the first study in the literature that uses a machine learning tool to predict shear localization.
Losses induced by tip clearance limit decisive improvements in the system efficiency and aerodynamic operational stability of aero-engine axial compressors. The tendency towards even lower blade heights to compensate for higher fluid densities aggravates their influence. Generally, it is emphasised that the tip clearance should be minimised but remain large enough to prevent collisions between the blade tip and the casing throughout the entire mission. The present work concentrates on the development of a preliminary aero-engine axial compressor casing design methodology involving meta-modelling techniques. Previous research work at the Institute for Turbomachinery and Flight Propulsion resulted in a Two-Dimensional (2D) axisymmetric finite element model for a generic multi-stage high-pressure axial compressor casing. Subsequent sensitivity studies led to the identification of significant parameters that are important for fine-tuning the tip clearance via specific flange design. This work is devoted to an exploration of the potential of surrogate modelling in preliminary compressor casing design with respect to rapid tip clearance assessments and its corresponding precision in comparison with finite element results. Reputed as data-driven mathematical approximation models and conceived for inexpensive numerical simulation result reproduction, surrogate models show even greater capacity when linked with extensive design space exploration and optimisation algorithms.
Compared with high-fidelity finite element simulations, the reductions obtained in computational time when using surrogate models amount to 99.9%. Validated via statistical methods and dependent on the size of the training database, the precision of surrogate models can reach down to the range of manufacturing tolerances. Subsequent inclusion of such surrogate models in a parametric optimisation process for tip clearance minimisation rapidly returned adaptions of the geometric design variables.
This innovative approach to teaching the finite element method blends theoretical, textbook-based learning with practical application using online and video resources. This hybrid teaching package features computational software such as MATLAB®, and tutorials presenting software applications such as PTC Creo Parametric, ANSYS APDL, ANSYS Workbench and SolidWorks, complete with detailed annotations and instructions so students can confidently develop hands-on experience. Suitable for senior undergraduate and graduate level classes, students will transition seamlessly between mathematical models and practical commercial software problems, empowering them to advance from basic differential equations to industry-standard modelling and analysis. Complete with over 120 end-of chapter problems and over 200 illustrations, this accessible reference will equip students with the tools they need to succeed in the workplace.
The finite element method is a powerful technique that can be used to transform any continuous body into a set of governing equations with a finite number of unknowns called degrees of freedom (DOF). In this chapter, we will introduce the fundamentals of the finite element method using a system of linear springs and a slender linear elastic body undergoing axial deformation as examples. These simple problems are chosen to describe the essential features of the finite element method which are common to analysis of more complicated structural systems such as 3D bodies.
To reveal the thermal shock resistance of double-layer thermal barrier coatings (TBCs), two types of TBCs were prepared via atmospheric plasma spraying, i.e., Gd2Zr2O7/yttria-stabilized zirconia (GZ/YSZ) TBCs and La2Zr2O7 (LZ)/YSZ TBCs, respectively. Subsequently, thermal cycling tests of the two TBCs were conducted at 1100 °C and their thermal shock resistance and failure mechanism were comparatively investigated through experiments and the finite element method. The results showed that the thermal shock failure of the two TBCs occurred inside the top ceramic coating. However, the GZ/YSZ TBCs had longer thermal cycling life. It was the mechanical properties of the top ceramic coating, and the thermal stresses arising from the thermal mismatch between the top ceramic coating and the substrate that determined the thermal cycling life of the two TBCs together. Compared with the LZ layer in the LZ/YSZ TBCs, the GZ layer in the GZ/YSZ TBCs had smaller elastic modulus, larger fracture toughness, and smaller thermal stresses, which led to the higher crack propagation resistance and less spallation tendency of the GZ/YSZ TBCs. Therefore, the GZ/YSZ TBCs exhibited superior thermal shock resistance to the LZ/YSZ TBCs.
The expressions for the elastic strain energy and its volumetric and deviatoric parts are derived for three-dimensional states of stress and strain. Betti's reciprocal theorem of linear elasticity is formulated, which yields the Maxwell coefficients, frequently used in structural mechanics. Castigliano's theorem is formulated and applied to axially loaded rods and trusses, twisted bars, and bent beams and frames. The principle of virtual work and the variational principle of linear elasticity are introduced. The differential equation of the deformed shape of the bent beam is derived from the consideration of the principle of virtual work. The approximate Rayleigh–Ritz method is introduced and applied to selected problems of structural mechanics. An introduction to the finite element method in the analysis of beam bending, torsion, and axial loading is then presented. The corresponding stiffness matrices and load vectors are derived for each element and are assembled into the global stiffness matrix and load vector of the entire structure.
Precise bone cut is fundamental in total knee arthroplasty. However, notching of anterior femoral is not uncommon in clinical practice. Reviewing the article, notching and its complication may reach up to 30% and 2.5%, and there is scanty study of notching on the femoral strength. We therefore conduct the finite element analysis to elucidate the effect of notching on femoral mechanical strength. The computerized tomography images were used as the basis to develop the knee model, which was assumed mainly to consist of cortical and cancellous bones. For the implant joint, Zimmer data was considered partly as the basis to develop the model. This study investigated the femoral improper cut effect on the surgery with a static standing condition. The results show that the anterior femoral cut should be undercut 2 mm to overcut 1 mm during the surgery, in order to prevent bone materials from yielding. The exposure of the cancellous bone may cause bone materials to yield when the femur overcut was 2 mm; the cancellous bone may load too much and result in a fracture when the undercut was 3 mm. The effect of undercut, which was rarely discussed, was particularly addressed in our study. Precise femoral cut is crucial for the longevity of total knee arthroplasty.