Finite Elements in Analysis and Design最新文献

Acoustic metamaterial structures have received extensive attention for sound and vibration engineering applications from the scientific community in recent years. However, the real-life application of conventional acoustic metamaterial structures is frequently limited by fixed frequency bands and increased structural thicknesses in low-frequency noise reduction. In this study, we introduce an origami-based acoustic metamaterial structure that consists of a Miura-ori foldcore, along with a perforated and an unperforated panel. The proposed Miura-ori foldcore sandwich acoustic metastructure (MOF-SAM) exhibits adjustable low-frequency sound absorption capacities due to the foldability of the origami foldcore. Moreover, we employ numerical methods to investigate the sound absorption properties of the MOF-SAM, quantified by the sound absorption coefficient. The results indicate that the structure has a single absorption peak which is superior to that of acoustic structures composed of conventional honeycomb cores. The dissipation of acoustic energy is due to the structural vibrations of the metastructure and the losses in the folding process of the origami foldcore. The numerical results of this study show that the proposed sound absorption mechanism enables tunable low-frequency sound absorption. The geometric design and periodicity of the origami unit fragments offer multiple distinct absorption peaks and thus tunable acoustic performance. These findings of this study are expected to inspire novel designs for next-generation acoustic devices.

The Hilbert spaces $H\left(\mathrm{curl}\right)$ and $H\left(\mathrm{div}\right)$ are employed in various variational problems formulated in the context of the de Rham complex in order to guarantee well-posedness. Seeing as the well-posedness follows automatically from the continuous setting to the discrete setting in the presence of commuting interpolants as per Fortin’s criterion, the construction of conforming subspaces becomes a crucial step in the formulation of stable numerical schemes. This work aims to introduce a novel, simple method of directly constructing semi-continuous vectorial base functions on the reference element via template vectors associated with the geometric polytopes of the element and an underlying $H^1$-conforming polynomial subspace. The base functions are then mapped from the reference element to the element in the physical domain via consistent Piola transformations. The method is defined in such a way, that the underlying $H^1$-conforming subspace can be chosen independently, thus allowing for constructions of arbitrary polynomial order. We prove a linearly independent construction of Nédélec elements of the first and second type, Brezzi–Douglas–Marini elements, and Raviart–Thomas elements on triangulations and tetrahedralizations. The application of the method is demonstrated with two examples in the relaxed micromorphic model.

This paper presents results and convergence study of the Global–Local Iterative Coupling through the implementation in the commercial software Abaqus making use of the co-simulation engine. A hierarchical modeling and simulation approach is often required to alleviate modeling burdens. Particular focus has been devoted here on convergence acceleration and performance optimization. Two applications in statics with nonlinear material behavior and geometrically nonlinear formulation are considered here: first a holed curved plate under traction with elastic–plastic material, then a pre-stressed bolted joint connecting two plates between each other and subjected to traction load. Three different convergence acceleration techniques are compared in terms of convergence performance and accuracy. An inexact solver strategy is proposed to improve computing time performance. The results show promising results for the coupling technology and constitute a step forward in the availability of non-intrusive multi-scale modeling capabilities for complex structures and assemblies.

This study investigates the numerical simulation of fracture behaviour in quasi-brittle materials like magnesia spinel refractories using the Gradient-Enhanced Damage (GED) model. It focuses on the complex modelling of these materials non-linear responses and compares conventional and variant GED models through a wedge splitting test. The results demonstrate that all GED models show a good fit to experimental data. However, the conventional GED model falls short in accurately depicting the fracture process zone. In contrast, the localizing GED model more accurately represents the fracture process zone, limiting spurious damage distribution, but requires finer meshing, elevating computational demands. The stress-based variant reduces spurious damage but is less effective comparatively. The study also assesses the role of heterogeneous strength distribution in replicating realistic crack patterns as observed in experiments.

In the present study, an elastoplastic phase-field model of quasi-static fracture in ductile materials is proposed in the variational framework for J₂ plasticity with isotropic hardening, which is suitable to describe the quasi-static behavior of metals as investigated in the performed experiments. These contributions include: (1) the free energy functions for coupling elastic response, plastic yielding and damage evolution are established. (2) The new elastic and plastic energy degradation functions are constructed to quantitatively describe the relationship between energy release and phase-field evolution of elastoplastic materials. (3) Damage evolution and plastic yielding criteria are derived. (4) From a numerical point of view, we derive the governing equations and the corresponding weak forms and the overall solution procedure for the phase-field model is given via the use of a return-mapping algorithm. This phase-field model was validated by a series of tensile experiments on Inconel 718 nickel-based super-alloys standard specimens. In order to compare the simulation results with the experimental results more comprehensively, the digital image correlation (DIC) technique is applied to experimentally investigate the specimen deformation information. In addition, to verify the potential of the model to capture complex cracks, we performed Nooru-Mohamed tests. The numerical simulation results are in good agreements with the results of previous experimental work.

This article presents the application of cubic Statistical Volume Elements (SVEs) to homogenize the elasticity tensor of epoxy matrix chopped glass fiber composites using displacement boundary conditions. A virtual microstructure reconstruction algorithm is used to reconstruct three large domains of the composites with different fiber orientation distributions. A non-iterative parallel meshing algorithm, named CISAMR, is then implemented to generate high-fidelity finite element models and simulate the linear elastic response of 1536 SVEs extracted from these domains. While the fiber orientations imply transversely isotropic elasticity stiffness matrices, for the SVE sizes considered, the composite is not quite transversely isotropic. We propose two indices of transverse isotropy to (1) determine the orientation at which a given property most closely matches the transversely isotropic assumption for an SVE, (2) quantify the corresponding transversely isotropic discrepancy, and (3) state the extent of transverse isotropy by measuring the difference between transverse and average normal quantities. The former can be applied to any orientation-dependent quantity such as strength, whereas the latter only applies to the elasticity tensor. We demonstrate the superiority of the latter for elastic properties and use the former to show that a proposed initiation fracture strength is farther away from its transversely isotropic limit compared to the directional elasticity normal stiffness.

This study examines the phenomenon of intrinsic nature in wave mitigation, specifically focusing on the concept of destructive interference (DI). When waves interact, they can exhibit either destructive interference or constructive interference depending on the phase difference. In the case of mechanical waves propagating through a mechanical structure, their characteristics such as wave speed, wavelength, and wave attenuation are influenced by the properties of the structure. As waves travel within the structure, the resonance phenomenon of the mechanical structure induces a phase shift of approximately 180 degrees in the wave. Consequently, what was initially destructive interference can transit into constructive interference, and vice versa. To address this challenge and systematically enhance the performance of mechanical structures employing destructive interference, a topology optimization scheme is applied. The concept of the present optimization scheme and the advantages are highlighted for several in-plane vibrations.

This article introduces a mixed domain decomposition method (DDM) designed to meet the requirements of advanced numerical optimization in electrical machines. The primary objective is to adapt the multiscale LATIN method, primarily used for mechanical studies, to the magnetostatic context. The proposed method offers an effective iterative scheme that relies on a mixed formulation of the equations on the domain interfaces, considering both primal and dual fields. A parallel computation strategy for the algorithm is implemented. Computational experiments performed demonstrate the scalability and efficiency of the algorithm.

This research continues the research of Deng et al. (2022) [1], using Discrete Element Method (DEM) coupled with Finite Element Analysis to solve shotgun exterior ballistics. The simulation examples in this research are using an Italian-made 24 gm #9½ birdshot with 433 pellets fired from 30” long, 12-gauge cylinder and full choke barrels. The simulations of shotgun exterior ballistics of this research included pellet swarm velocity and the pellet dispersion at different distances until 50 yards away from the muzzle. The ballistics simulation of the pellet swarm is completed from interior to exterior consecutively after the shotshell is fired inside the chamber so all ballistics performances can be calculated at one time. Three forces were applied to the pellets for exterior ballistics simulation: the contact force between pellets, the aerodynamic separation force between pellets, and the drag force. Because of the complexity of the aerodynamic forces exerted on pellets, this research used an equivalent aerodynamic force to simulate this complex phenomenon. Two birdshot models with different pellet formations were created; the first one was simulated to calibrate the separation scale factor defined in aerodynamic separation force, and the second one was used for validation and sensitivity of the model. The simulation results show that for #9½ birdshots fired by cylinder barrel, the average Effective Shot Dispersion (ESD) of pellet dispersion of both birdshots inside the 30” diameter of the target circle at 40 yards from the muzzle is 398.53, which is remarkably close to 396.98 of the experiment result. The simulation of the average pellets' target hit rate is 77.14% (inside the 30” diameter of the target circle), which is also remarkably close to the experiment hit rate of 77.57%. The same birdshot fired from a full choke barrel shows that ESD and hit rate rose to 406.34 and 83.14%, respectively. These results demonstrate the effectiveness of using the discrete element method in conjunction with the proposed equivalent aerodynamic force to predict the shotgun's interior and exterior ballistics.

Topology optimization applied to fluid–structure interaction problems is challenging because the physical phenomenon in real engineering applications is usually transient and strongly coupled. This leads to costly solutions for the forward and adjoint problems, the computational bottleneck of the topology optimization method. Thus, this paper proposes a topology optimization problem formulated in the steady state with post-processing and verification in the transient state. The objective is to design a stiff structure with lower effects of vibrations induced by the transient fluid vortices. For that, the compliance minimization problem is solved subject to a natural frequency constraint (without any volume constraint). The TOBS-GT (Topology Optimization of Binary Structures with geometry trimming) method is used to solve the problem. To observe the vortex-shedding around the structure, a transient simulation is performed considering an incompressible fluid flow under a laminar regime and the structure subject to large displacements. For topology optimization, the fluid flow is at a steady state and the structure is modeled considering small displacements, i.e., a one-way coupled analysis. The finite element method is used to solve the governing equations and obtain the direct/adjoint sensitivities for the compliance and natural frequency functions. In this approach, the natural frequency of the structure is shifted away from the fluid flow vortex-shedding frequency, avoiding resonance. Numerical examples show that the proposed method can be effectively applied to design 2D structures in FSI problems with lower effects of Flow-Induced Vibration, attenuating the levels of displacement at the analyzed points of the structure.