This paper introduces a novel formulation, called Hybrid Discrete-Finite Element (HybriDFEM) method, for modeling one-directional continuous and discontinuous planar beam-like members, including nonlinear geometric and material effects. In this method, the structure is modeled as a series of distinct rigid blocks, connected to each other through contact pairs distributed along the interfaces. Each of those contact pairs are composed of two nonlinear multidirectional springs in series, which can represent either the deformation of the blocks themselves, or the deformation of their interface. Unlike the Applied Element Method, in which contact pairs are composed of one single spring, the current approach allows capturing phenomena such as sectional deformations or relative deformations between two blocks composed of different materials. This method shares similarities with the Discrete Element Methods in its ability to model contact interfaces between rigid or deformable units, but does not require a numerical time-domain integration scheme. More importantly, its formulation resembles that of the classical Finite Elements Method, allowing one to easily couple the latter with HybriDFEM. Following the presentation of its formulation, the method is benchmarked against analytical solutions selected from the literature, ranging from the linear-elastic response of a cantilever beam to the buckling and rocking response of continuous flexible columns, and rigid block stackings. One final example showcases the coupling of a HybriDFEM element with a linear beam finite element.

This paper presents a new discrete parametrization method for simultaneous topology and material optimization of composite laminate structures, referred to as Hyperbolic Function Parametrization (HFP). The novelty of HFP is the way the candidate materials are parametrized in the optimization problem. In HFP, a filtering technique based on hyperbolic functions is used, such that only one design variable is used for any given number of material candidates. Compared to state-of-the-art methods such Discrete Material and Topology Optimization (DMTO) and Shape Function with Penalization (SFP), HFP has much fewer optimization variables and constraints but introduces additional non-linearity in the optimization problems. A comparative analysis of HFP, DMTO and SFP are performed based on the problem of maximizing the stiffness of composite plates under a total volume constraint and multiple manufacturing constraints using various loads, boundary conditions and input parameters. The comparison shows that all three methods are highly sensitive to the choice of input parameters for the optimization problem, although the performance of HFP is overall more consistent. HFP method performs similarly to DMTO and SFP in terms of the designs obtained and computational cost. However, HFP obtains similar or better objective function values compared to the DMTO and SFP methods. © 2021 The Author(s)

We propose a new deterministic robust design optimization method for composite laminate structures under worst-case material uncertainty. The method is based on a simultaneous parametrization of topology and material and combines a design problem and a material uncertainty problem into a single min–max optimization problem which provides an efficient approach to handle variation of material properties in stiffness driven design optimization problems. An analysis is performed using a design problem based on a failure criterion formulation to evaluate the ability of the proposed method to generate robust composite designs. The design problem is solved using various loads, boundary conditions and manufacturing constraints. The designs generated with the proposed method have improved objective responses compared to the worst-case response of designs generated with nominal material properties and are less sensitive to the variation of material properties. The analysis indicates that the proposed method can be efficiently applied in a robust structural optimization framework. © 2023 The Author(s)

We propose a method to analyze effects of material uncertainty in composite laminate structures optimized using a simultaneous topology and material optimization approach. The method is based on computing worst-case values for the material properties and provides an efficient way of handling variation in material properties of composites for stiffness driven optimization problems. An analysis is performed to evaluate the impact of material uncertainty on designs from two design problems: Maximization of stiffness and minimization of a failure criteria index, respectively. The design problems are solved using different loads, boundary conditions and manufacturing constraints. The analysis indicates that the influence of material uncertainty is dependent on the type of optimization problem. For compliance problems the impact on the objective value is proportional to the changes of the constitutive properties and the effect of material uncertainty is consistent and predictable for the generated designs. The strength-based problem shows that material uncertainty has a significant impact on the response, and the effects of material uncertainty is not consistent and changes for different design requirements. In addition, the results show an increase of up to 25% of the maximum failure index when considering the worst-case deviation of the constitutive properties from their nominal values. © 2022 The Author(s)

We consider nonlinear shape effects appearing in the lumped electromechanical model of a bimorph piezoelectric bridge structure due to the interaction between the electromechanical constitutive model and the geometry of the structure. At finite proof-mass displacement and electrode voltage, the shape of the beams is no longer given by Euler-Bernoulli theory which implies that shape effects enter in both the electrical and mechanical domains and in the coupling between them. Accounting for such effects is important for the accurate modelling of, e.g., piezoelectrical energy harvesters and actuators in the regime of large deflections and voltages. We present a general method, based on a variational approach minimizing the Gibbs enthalpy of the system, for computing corrections to the nominal shape function and the associated corrections to the lumped model. The lowest order correction is derived explicitly and is shown to produce significant improvements in model accuracy, both in terms of the Gibbs enthalpy and the shape function itself, over a large range of displacements and voltages. Furthermore, we validate the theoretical model using large deflection finite element simulations of the bridge structure and conclude that the lowest order correction substantially improve the model, obtaining a level of accuracy expected to be sufficient for most applications. Finally, we derive the equations of motion for the lowest order corrected model and show how the coupling between the electromechanical properties and the geometry of the bridge structure introduces nonlinear interaction terms. © 2021, The Author(s).

We introduce numerical algebraic geometry methods for computing lower bounds on the reach, local feature size, and weak feature size of the real part of an equidimensional and smooth algebraic variety using the variety’s defining polynomials as input. For the weak feature size, we also show that nonquadratic complete intersections generically have finitely many geometric bottlenecks, and we describe how to compute the weak feature size directly rather than a lower bound in this case. In all other cases, we describe additional computations that can be used to determine feature size values rather than lower bounds. 

Motivated by the substantial interest in various fillers to enhance the barrier properties of polymeric films, especially graphene derivatives, we perform a computational screening of obstructed diffusion to explore the design parameter space of nanoplatelet-filled composites synthesized in silico. As a model for the nanoplatelets, we use circular and elliptical nonoverlapping and impermeable flat disks, and diffusion is stochastically simulated using a random-walk model, from which the effective diffusivity is calculated. On the basis of ∼1000 generated structures and diffusion simulations, we systematically investigate the impact of different nanoplatelet characteristics such as orientation, layering, size, polydispersity, shape, and amount. We conclude that the orientation, size, and amount of nanoplatelets are the most important parameters and show that using nanoplatelets oriented perpendicular to the diffusion direction, under reasonable assumptions, with approximately 0.2% (w/w) graphene, we can reach 90% reduction and, with approximately 1% (w/w) graphene, we can reach 99% reduction in diffusivity, purely because of geometrical effects, in a defect-free matrix with perfect compatibility. Additionally, our results suggest that the existing analytical models have some difficulty with extremely large aspect ratio (extremely flat) nanoplatelets, which calls for further development.

This paper presents ISOPE’s 2020 comparative study on the interaction between focused waves and a fixed cylinder. The paper discusses the qualitative and quantitative comparisons between 20 different numerical solvers from various universities across the world for a fixed cylinder. The moving cylinder cases are reported in a companion paper as part B (Agarwal, Saincher, et al., 2021). The numerical solvers presented in this paper are the recent state of the art in the field, mostly developed in-house by various academic institutes. The majority of the participants used hybrid modeling (i.e., a combination of potential flow and Navier–Stokes solvers). The qualitative comparisons based on the wave probe and pressure probe time histories and spectral components between laminar, turbulent, and potential flow solvers are presented in this paper. Furthermore, the quantitative error analyses based on the overall relative error in peak and phase shifts in the wave probe and pressure probe of all the 20 different solvers are reported. The quantitative errors with respect to different spectral component energy levels (i.e., in primary, sub-, and superharmonic regions) capturing capability are reported. Thus, the paper discusses the maximum, minimum, and median relative errors present in recent solvers as regards application to industrial problems rather than attempting to find the best solver. Furthermore, recommendations are drawn based on the analysis. 

The spectral/hp element method combines the geometric flexibility of the classical h-type finite element technique with the desirable numerical properties of spectral methods, employing high-degree piecewise polynomial basis functions on coarse finite element-type meshes. The spatial approximation is based upon orthogonal polynomials, such as Legendre or Chebychev polynomials, modified to accommodate a C0 - continuous expansion. Computationally and theoretically, by increasing the polynomial order p , high-precision solutions and fast convergence can be obtained and, in particular, under certain regularity assumptions an exponential reduction in approximation error between numerical and exact solutions can be achieved. This method has now been applied in many simulation studies of both fundamental and practical engineering flows. This paper briefly describes the formulation of the spectral/hp element method and provides an overview of its application to computational fluid dynamics. In particular, it focuses on the use of the spectral/hp element method in transitional flows and ocean engineering. Finally, some of the major challenges to be overcome in order´to use the spectral/hp element method in more complex science and engineering applications are discussed

As one of the innovative technologies of intelligent transportation systems (ITS), Connected and Autonomous Vehicles (CAVs) have been deployed gradually. Given that there will be a long transition period before reaching a fully CAVs environment, it is crucial to assess the potential impacts of CAVs on mixed traffic flow. Considering platoon formation process, this study develops a platoon cooperation strategy based on “catch-up” mechanism, and then analyzes the impact on fundamental diagram, traffic oscillation, and traffic safety within mixed traffic. Simulation results show that with an increasing market penetration rate (MPR) of CAVs, road capacity shows an increasing trend. Compared with base scenario, a clear increase in road capacity is also observed under platoon scenario. With an increasing MPR, traffic oscillation is shown to reduce largely. Furthermore, the proposed platoon strategy could dampen frequent shockwaves and shorten the propagation range of waves. Regarding traffic safety, multiple surrogate safety measures (SSMs) are used to evaluate the traffic risk: including Criticality Index Function (CIF), Potential Index for Collision with Urgent Deceleration (PICUD), and Deceleration Rate to Avoid a Crash (DRAC). With increasing MPR, collision risk identified by CIF and DRAC shows an increase tendency, while that identified by PICUD has no apparent trend. Furthermore, the platoon strategy is shown to increase the severity of traffic conflicts significantly. Overall, this study provides novel insights into CAVs deployment through the analysis of platoon strategy.