The majority of topology optimization methods for porous infill designs is based on the assumption of deterministic loads. However, in practice, quantities such as positions, weights, and directions ...of applied loads may change accidentally. Deterministic load-based designs might deliver poor structural performance under loading uncertainties. Such uncertain factors need to be taken into account in topological optimization to seek robust results. This paper presents a novel robust concurrent topology optimization method for the design of uniform/non-uniform porous infills under the accidental change of loads. A combination of moving morphable bars (MMBs) and loading uncertainties is proposed to directly model multiscale structures and seek robust designs. The macro- and microscopic structures can be simultaneously optimized through the minimization of the weighted sum of the expected compliance and standard deviation. The geometries of adaptive geometric components (AGCs) are straightforwardly optimized. The AGCs consist of two classes of geometric components: macroscopic bars describing the overall structure and microscopic bars describing the material microstructures. Automatic mesh-refinement is utilized to enhance computing efficiency. Numerical examples demonstrate that robust porous design can be obtained with only one global volume constraint while the material continuity of neighboring unit cells and the structural porosity can be maintained without additional constraints. The robust designs yield a more robust structural performance along with a smaller standard deviation compared with deterministic porous designs under loading uncertainties.
This paper proposes an adaptive atomistic- continuum numerical method for quasi-static crack growth. The phantom node method is used to model the crack in the continuum region and a molecular statics ...model is used near the crack tip. To ensure self-consistency in the bulk, a virtual atom cluster is used to model the material of the coarse scale. The coupling between the coarse scale and fine scale is realized through ghost atoms. The ghost atom positions are interpolated from the coarse scale solution and enforced as boundary conditions on the fine scale. The fine scale region is adaptively enlarged as the crack propagates and the region behind the crack tip is adaptively coarsened. An energy criterion is used to detect the crack tip location. The triangular lattice in the fine scale region corresponds to the lattice structure of the (111) plane of an FCC crystal. The Lennard-Jones potential is used to model the atom–atom interactions. The method is implemented in two dimensions. The results are compared to pure atomistic simulations; they show excellent agreement.
A thorough understanding of brain metabolism is essential to tackle neurodegenerative diseases. Astrocytes are glial cells which play an important metabolic role by supplying neurons with energy. In ...addition, astrocytes provide scaffolding and homeostatic functions to neighboring neurons and contribute to the blood–brain barrier. Recent investigations indicate that the complex morphology of astrocytes impacts upon their function and in particular the efficiency with which these cells metabolize nutrients and provide neurons with energy, but a systematic understanding is still elusive. Modelling and simulation represent an effective framework to address this challenge and to deepen our understanding of brain energy metabolism. This requires solving a set of metabolic partial differential equations on complex domains and remains a challenge. In this paper, we propose, test and verify a simple numerical method to solve a simplified model of metabolic pathways in astrocytes. The method can deal with arbitrarily complex cell morphologies and enables the rapid and simple modification of the model equations by users also without a deep knowledge in the numerical methods involved. The results obtained with the new method (
CutFEM
) are as accurate as the finite element method (FEM) whilst
CutFEM
disentangles the cell morphology from its discretisation, enabling us to deal with arbitrarily complex morphologies in two and three dimensions.
In this paper we develop a framework for solving inverse deformation problems using the FEniCS Project finite-element software. We validate our approach with experimental imaging data acquired from a ...soft silicone beam under gravity. In contrast with inverse iterative algorithms that require multiple solutions of a standard elasticity problem, the proposed method can compute the undeformed configuration by solving only one modified elasticity problem. This modified problem has a complexity comparable to the standard one. The framework is implemented within an open-source pipeline enabling the direct and inverse deformation simulation directly from imaging data. We use the high-level unified form language (UFL) of the FEniCS Project to express the finite-element model in variational form and to automatically derive the consistent Jacobian. Consequently, the design of the pipeline is flexible: for example, it allows the modification of the constitutive models by changing a single line of code. We include a complete working example showing the inverse deformation of a beam deformed by gravity as supplementary material.
With ageing populations, understanding life course factors that raise the risk of depression in old age may help anticipate needs and reduce healthcare costs in the long run. We estimate the risk of ...depression in old age by combining adult life course trajectories and childhood conditions in supervised machine learning algorithms. Using data from the Survey of Health, Ageing and Retirement in Europe (SHARE), we implement and compare the performance of six alternative machine learning algorithms. We analyse the performance of the algorithms using different life-course data configurations. While we obtain similar predictive abilities between algorithms, we achieve the highest predictive performance when employing semi-structured representations of life courses using sequence data. We use the Shapley Additive Explanations method to extract the most decisive predictive patterns. Age, health, childhood conditions, and low education predict most depression risk later in life, but we identify new predictive patterns in indicators of life course instability and low utilization of dental care services.
•We predict old age depression risk by combining retrospective life course data and machine learning models.•Models are estimated on SHARElife biographical data.•The best-performing models are those that include high-dimensional, semi-structured data.•Depression in women is more foreseeable than in men.•Our method uncovers under-studied predictors such as emotional-status changes and low lifetime utilization of dental care services.
A smoothed finite element method for shell analysis Nguyen-Thanh, N.; Rabczuk, Timon; Nguyen-Xuan, H. ...
Computer methods in applied mechanics and engineering,
12/2008, Letnik:
198, Številka:
2
Journal Article
Recenzirano
A four-node quadrilateral shell element with smoothed membrane-bending based on Mindlin–Reissner theory is proposed. The element is a combination of a plate bending and membrane element. It is based ...on mixed interpolation where the bending and membrane stiffness matrices are calculated on the boundaries of the smoothing cells while the shear terms are approximated by independent interpolation functions in natural coordinates. The proposed element is robust, computationally inexpensive and free of locking. Since the integration is done on the element boundaries for the bending and membrane terms, the element is more accurate than the MITC4 element for distorted meshes. This will be demonstrated for several numerical examples.
We provide a primer to numerical methods based on Taylor series expansions such as generalized finite difference methods and collocation methods. We provide a detailed benchmarking strategy for these ...methods as well as all data files including input files, boundary conditions, point distribution and solution fields, so as to facilitate future benchmarking of new methods. We review traditional methods and recent ones which appeared in the last decade. We aim to help newcomers to the field understand the main characteristics of these methods and to provide sufficient information to both simplify implementation and benchmarking of new methods. Some of the examples are chosen within a subset of problems where collocation is traditionally known to perform sub-par, namely when the solution sought is non-smooth, i.e. contains discontinuities, singularities or sharp gradients. For such problems and other simpler ones with smooth solutions, we study in depth the influence of the weight function, correction function, and the number of nodes in a given support. We also propose new stabilization approaches to improve the accuracy of the numerical methods. In particular, we experiment with the use of a Voronoi diagram for weight computation, collocation method stabilization approaches, and support node selection for problems with singular solutions. With an appropriate selection of the above-mentioned parameters, the resulting collocation methods are compared to the moving least-squares method (and variations thereof), the radial basis function finite difference method and the finite element method. Extensive tests involving two and three dimensional problems indicate that the methods perform well in terms of efficiency (accuracy versus computational time), even for non-smooth solutions.