Search

In this work, a known Eddy Current (EC) model is adapted to characterize subsurface defects in carbon fiber-reinforced polymer (CFRP) plates intended for the civil aerospace industry. The considered defects include delaminations, microcracks, porosity, fiber breakage, and the simultaneous presence of these defects. Each defect is modeled as an additive variation in the material’s electrical conductivity tensor, allowing for a detailed mathematical representation of the defect’s influence on the CFRP’s electromagnetic behavior. The additivity of the variations in the conductivity tensor is justified by the assumption that the defects are not visible to the naked eye, implying that the material does not require non-destructive testing. The adapted EC model admits a unique and stable solution by verifying that all analytical steps are satisfied. To reconstruct 2D maps of the magnetic flux density amplitude, a FEM formulation is adopted, based on the energy functional because it ensures a stable and consistent numerical formulation given its coercivity. Moreover, the numerical approach allows precise and reliable numerical solutions, enhancing the capability to detect and quantify defects. The numerical results show that the obtained 2D maps are entirely superimposable on those highlighting the distribution of mechanical stress states known in the literature, offering a clear advantage in terms of detection costs. This approach provides an effective and economical solution for the non-destructive inspection of CFRP, ensuring accurate and timely defect diagnosis for maintaining structural integrity.

An innovative non-homogeneous dynamic model is presented for the recovery of temperature during the industrial laser welding process of Al-Si 5% alloy plates. It considers that, metallurgically, during welding, the alloy melts with the presence of solid/liquid phases until total melt is achieved, and afterwards it resolidifies with the reverse process. Further, a polynomial substitute thermal capacity of the alloy is chosen based on experimental evidence so that the volumetric solid-state fraction is identifiable. Moreover, to the usual radiative/convective boundary conditions, the contribution due to the positioning of the plates on the workbench is considered (endowing the model with Cauchy–Stefan–Boltzmann boundary conditions). Having verified the well-posedness of the problem, a Galerkin-FEM approach is implemented to recover the temperature maps, obtained by modeling the laser heat sources with formulations depending on the laser sliding speed. The results achieved show good adherence to the experimental evidence, opening up interesting future scenarios for technology transfer.

In this paper, starting from a well-known nonlinear hyperbolic integro-differential model of the fourth order describing the dynamic behavior of an electrostatic MEMS with a parallel plate, the authors propose an upgrade of it by formulating an additive term due to the effects produced by the fringing field and satisfying the Pelesko–Driscoll theory, which, as is well known, has strong experimental confirmation. Exploiting the theory of hyperbolic equations in Hilbert spaces, and also utilizing Campanato’s Near Operator Theory (and subsequent applications), results of existence and regularity of the solution are proved and discussed particularly usefully in anticipation of the development of numerical approaches for recovering the profile of the deformable plate for a wide range of applications.

It is now clear that cooperation between academia and industries is crucial for social, cultural, technological and economic progress and innovation [...]

In this paper, we prove the existence and uniqueness of solutions for a nonlocal, fourth-order integro-differential equation that models electrostatic MEMS with parallel metallic plates by exploiting a well-known implicit function theorem on the topological space framework. As the diameter of the domain is fairly small (similar to the length of the device wafer, which is comparable to the distance between the plates), the fringing field phenomenon can arise. Therefore, based on the Pelesko–Driscoll theory, a term for the fringing field has been considered. The nonlocal model obtained admits solutions, making these devices attractive for industrial applications whose intended uses require reduced external voltages.

An important problem in membrane micro-electric-mechanical-system (MEMS) modeling is the fringing-field phenomenon, of which the main effect consists of force-line deformation of electrostatic field E near the edges of the plates, producing the anomalous deformation of the membrane when external voltage V is applied. In the framework of a 2D circular membrane MEMS, representing the fringing-field effect depending on |∇u|2 with the u profile of the membrane, and since strong E produces strong deformation of the membrane, we consider |E| proportional to the mean curvature of the membrane, obtaining a new nonlinear second-order differential model without explicit singularities. In this paper, the main purpose was the analytical study of this model, obtaining an algebraic condition ensuring the existence of at least one solution for it that depends on both the electromechanical properties of the material constituting the membrane and the positive parameter δ that weighs the terms |∇u|2. However, even if the the study of the model did not ensure the uniqueness of the solution, it made it possible to achieve the goal of finding a stable equilibrium position. Moreover, a range of admissible values of V were obtained in order, on the one hand, to win the mechanical inertia of the membrane and, on the other hand, to ensure that the membrane did not touch the upper disk of the device. Lastly, some optimal control conditions based on the variation of potential energy are presented and discussed.

When a straight cylindrical conductor of finite length is electrostatically charged, its electrostatic potential ϕ depends on the electrostatic charge qe, as expressed by the equation L(qe)=ϕ, where L is an integral operator. Method of moments (MoM) is an excellent candidate for solving L(qe)=ϕ numerically. In fact, considering qe as a piece-wise constant over the length of the conductor, it can be expressed as a finite series of weighted basis functions, qe=∑n=1Nαnfn (with weights αn and N, number of the subsections of the conductor) defined in the L domain so that ϕ becomes a finite sum of integrals from which, considering testing functions suitably combined with the basis functions, one obtains an algebraic system Lmnαn=gm with dense matrix, equivalent to L(qe)=ϕ. Once solved, the linear algebraic system gets αn and therefore qe is obtainable so that the electrostatic capacitance C=qe/V, where V is the external electrical tension applied, can give the corresponding electrostatic capacitance. In this paper, a comparison was made among some Krylov subspace method-based procedures to solve Lmnαn=gm. These methods have, as a basic idea, the projection of a problem related to a matrix A∈Rn×n, having a number of non-null elements of the order of n, in a subspace of lower order. This reduces the computational complexity of the algorithms for solving linear algebraic systems in which the matrix is dense. Five cases were identified to determine Lmn according to the type of basis-testing functions pair used. In particular: (1) pulse function as the basis function and delta function as the testing function; (2) pulse function as the basis function as well as testing function; (3) triangular function as the basis function and delta function as the testing function; (4) triangular function as the basis function and pulse function as the testing function; (5) triangular function as the basis function with the Galerkin Procedure. Therefore, five Lmn and five pair qe and C were computed. For each case, for the resolution of Lmnαn=gm obtained, GMRES, CGS, and BicGStab algorithms (based on Krylov subspaces approach) were implemented in the MatLab® Toolbox to evaluate qe and C as N increases, highlighting asymptotical behaviors of the procedures. Then, a particular value for N is obtained, exploiting both the conditioning number of Lmn and considerations on C, to avoid instability phenomena. The performances of the exploited procedures have been evaluated in terms of convergence speed and CPU-times as the length/diameter and N increase. The results show the superiority of BcGStab, compared to the other procedures used, since even if the number of iterations increases significantly, the CPU-time decreases (more than 50%) when the asymptotic behavior of all the procedures is in place. This superiority is much more evident when the CPU-time of BicGStab is compared with that achieved by exploiting Gauss elimination and Gauss–Seidel approaches.

In this paper, a stable numerical approach for recovering the membrane profile of a 2D Micro-Electric-Mechanical-Systems (MEMS) is presented. Starting from a well-known 2D nonlinear second-order differential model for electrostatic circular membrane MEMS, where the amplitude of the electrostatic field is considered proportional to the mean curvature of the membrane, a collocation procedure, based on the three-stage Lobatto formula, is derived. The convergence is studied, thus obtaining the parameters operative ranges determining the areas of applicability of the device under analysis.

The evaluation of the cutoff wavenumbers of metallic waveguides can be related to the numerical resolution of a suitable nonlinear eigenproblem defined on the domain C described by the contour of its transverse cross section. In this work, we show that the symmetries of C can be exploited to obtain a block diagonal matrix representation of the nonlinear eigenproblem, which enables a remarkable reduction in the computational effort involved in its resolution.