Search results

This study implements the fourth-order phase field method (PFM) for modeling fracture in brittle materials. The weak form of the fourth-order PFM requires C¹ basis functions for the crack evolution scalar field in a finite element framework. To address this, non-Sibsonian type shape functions that are nonpolynomial types based on distance measures, are used in the context of natural neighbor shape functions. The capability and efficiency of this method are studied for modeling cracks.

The weak form of the fourth-order PFM is derived from two governing equations for finite element modeling. C⁰ non-Sibsonian shape functions are derived using distance measures on a generalized quad element. Then these shape functions are degree elevated with Bernstein-Bezier (BB) patch to get higher-order continuity (C¹) in the shape function. The quad element is divided into several background triangular elements to apply the Gauss-quadrature rule for numerical integration. Both fourth-order and second-order PFMs are implemented in a finite element framework. The efficiency of the interpolation function is studied in terms of convergence and accuracy for capturing crack topology in the fourth-order PFM.

It is observed that fourth-order PFM has higher accuracy and convergence than second-order PFM using non-Sibsonian type interpolants. The former predicts higher failure loads and failure displacements compared to the second-order model due to the addition of higher-order terms in the energy equation. The fracture pattern is realistic when only the tensile part of the strain energy is taken for fracture evolution. The fracture pattern is also observed in the compressive region when both tensile and compressive energy for crack evolution are taken into account, which is unrealistic. Length scale has a certain specific effect on the failure load of the specimen.

Fourth-order PFM is implemented using C¹ non-Sibsonian type of shape functions. The derivation and implementation are carried out for both the second-order and fourth-order PFM. The length scale effect on both models is shown. The better accuracy and convergence rate of the fourth-order PFM over second-order PFM are studied using the current approach. The critical difference between the isotropic phase field and the hybrid phase field approach is also presented to showcase the importance of strain energy decomposition in PFM.

The main objective of this study is to develop a numerical model based on Isogeometric Analysis to study the dynamic behavior of multi-directional functionally graded plates with variable thickness.

A numerical study was conducted on the dynamic behavior of multi-directional functionally graded plates. Rectangular and circular plates with variable thickness are taken into investigation. The third-order shear deformation plate theory of Reddy is used to describe the displacement field, while the equation of motion is developed based on the Hamilton's principle. Isogeometric Analysis approach is employed as a discretization tool to develop the system equation, where NURBS basis functions are used. The famous Newmark method is used to solve time-dependent problems.

The results obtained from this study indicated that the thickness gradation has a more considerable effect than in-plane variation of materials in MFGM plates. Additionally, the influence of the damping factor is observed to affect the vibration amplitude of the plate. The results obtained from this study could be used for future investigations, where the viscous elasticity and other dynamic factors are considered.

Although there have been a number of studies in the literature devoted to analyzing the linear static bending and free vibration of FGM and MFGM plates with variable thickness, the study on dynamic response of FGM and MFGM plate is still limited. Therefore, this study is dedicated to the investigation of the dynamic behavior of multi-directional functionally graded plates.

This research presents a design method for designing greenhouse structures based on topology optimization. Moreover, the structural design of a gothic greenhouse is proposed in which its structural strength has been improved by using this proposed method. In this method, the design of the structure is done mathematically; therefore, in the design process, more attention can be focused on the constraint space and boundary conditions. It was also shown how the static reliability and fatigue coefficients will change as a result of the design of the greenhouse structure with this method. Another purpose of this study is to find the weakest part of the greenhouse structure against lateral winds and other general loads on the greenhouse structure.

In the proposed method, the outer surface and the allowable volume as a constraint domain were considered. The desired loads can be located on the constraint domain. The topology optimization was used to minimize the mass and structural compliance as the objective function. The obtained volume was modified for simplifying the construction. The changes in the shape of the greenhouse structure were investigated by choosing three different penalty numbers for the topology optimization algorithm. The final design of the proposed structure was performed based on the total simultaneous critical loads on the structure. The results of the proposed method were compared in the order of different volume fractions. This showed that the volume fraction approach can significantly reduce the weight of the structure while maintaining its strength and stability.

Topology optimization results showed different strut and chords composition because of the changes in maximum mass limit and volume fraction. The results showed that the fatigue was more hazardous, and it decreased the strength of structure nearly three times more than a static analysis. Further, it was noticed that how the penalty numbers can affect topology optimization results. An optimal design based on topology optimization results was presented to improve the proposed greenhouse design against destruction and demolition. Furthermore, this study shows the most sensitive part of the greenhouse against the standard loads of wind, snow, and crop.

The obtained designs were compared with a conventional arch greenhouse, and then the structural performances were shown based on standard loads. The results showed that in designing the proposed structure, the optimized changes increased the structure strength against the standard loads compared to a simple arch greenhouse. Moreover, the stress safety factor and fatigue safety factor because of different designs of this structure were also compared with each other.

This study presents finite element (FE) and generalized regression neural network (GRNN) approaches for modeling multiple crack growth problems and predicting crack-growth directions under the influence of multiple crack parameters.

To determine the crack-growth direction in aluminum specimens, multiple crack parameters representing some degree of crack propagation complexity, including crack length, inclination angle, offset and distance, were examined. FE method models were developed for multiple crack growth simulations. To capture the complex relationships among multiple crack-growth variables, GRNN models were developed as nonlinear regression models. Six input variables and one output variable comprising 65 training and 20 test datasets were established.

The FE model could conveniently simulate the crack-growth directions. However, several multiple crack parameters could affect the simulation accuracy. The GRNN offers a reliable method for modeling the growth of multiple cracks. Using 76% of the total dataset, the NN model attained an R2 value of 0.985.

The models are presented for static multiple crack growth problems. No material anisotropy is observed.

In practical crack-growth analyses, the NN approach provides significant benefits and savings.

The proposed GRNN model is simple to develop and accurate. Its performance was superior to that of other NN models. This model is also suitable for modeling multiple crack growths with arbitrary geometries. The proposed GRNN model demonstrates its prediction capability with a simpler learning process, thus producing efficient multiple crack growth predictions and assessments.

The purpose of this study is to extract the response of the simultaneous low-velocity impact of multiple impactors on a porous functionally graded (FG) aluminum plate.

To design a porous FG structure, a series of functions are applied using the porosity coefficient, and mechanical properties including Young’s modulus, shear modulus and the density of the porous structure are presented as a function of the axis placed in the direction of the plate thickness. The first-order shear deformation theory of the plate is used. To simulate the contact process between each impactor and the plate, a nonlinear Hertz contact force is considered for that impactor independently.

ABAQUS finite element software is used for the verification process of the theorical equations. The effects of porous function type, radius and initial velocity of impactor are investigated for the simultaneous impact of five impactors on porous FG aluminum plate with a simply supported boundary condition. Histories of contact force and displacement of the impactor placed in the center of the beam are analyzed in detail with the changes of the mentioned parameters.

Due to the advantages of porous aluminum plate such as high energy absorption and low weight, such structures may be subjected to the simultaneous impact of multiple impactors, which is studied in this research.

This study aims to propose and numerically assess different ways of discretising a very weak formulation of the Poisson problem.

We use integration by parts twice to shift smoothness requirements to the test functions, thereby allowing low-regularity data and solutions.

Various conforming discretisations are presented and tested, with numerical results indicating good accuracy and stability in different types of problems.

This is one of the first articles to propose and test concrete discretisations for very weak variational formulations in primal form. The numerical results, which include a problem based on real MRI data, indicate the potential of very weak finite element methods for tackling problems with low regularity.

Bayesian cubature (BC) has emerged to be one of most competitive approach for estimating the multi-dimensional integral especially when the integrand is expensive to evaluate, and alternative acquisition functions, such as the Posterior Variance Contribution (PVC) function, have been developed for adaptive experiment design of the integration points. However, those sequential design strategies also prevent BC from being implemented in a parallel scheme. Therefore, this paper aims at developing a parallelized adaptive BC method to further improve the computational efficiency.

By theoretically examining the multimodal behavior of the PVC function, it is concluded that the multiple local maxima all have important contribution to the integration accuracy as can be selected as design points, providing a practical way for parallelization of the adaptive BC. Inspired by the above finding, four multimodal optimization algorithms, including one newly developed in this work, are then introduced for finding multiple local maxima of the PVC function in one run, and further for parallel implementation of the adaptive BC.

The superiority of the parallel schemes and the performance of the four multimodal optimization algorithms are then demonstrated and compared with the k-means clustering method by using two numerical benchmarks and two engineering examples.

Multimodal behavior of acquisition function for BC is comprehensively investigated. All the local maxima of the acquisition function contribute to adaptive BC accuracy. Parallelization of adaptive BC is realized with four multimodal optimization methods.

The purpose of this paper is to investigate the free vibration response of a laminated honeycomb sandwich panels (LHSP) for aerospace applications. Higher order shear deformation theory (HSDT) was simplified for the dynamic analysis of LHSP. Furthermore, the effects of honeycomb parameters on the value of natural frequency (NF) of vibration were explored.

This paper applies HSDT to the analysis of composite LHSP to derive four vibration differential equations of motion and solve it to find the NF of vibration. Two analytical models (Nayak and Meunier models) were selected from literature for comparison of the NF of vibration. In addition, a numerical model was built by using ABAQUS and the results were compared. Furthermore, parametric studies were conducted to explore the effect of honeycomb parameters on the value of the NF of vibration.

The present model is successful in simplifying HSDT for the analysis of LHSP. The first five natural frequencies of vibration were calculated analytically and numerically. In the parametric study, increasing core height or young’s modulus or changing laminate layup will increase the value of NF of vibration. Furthermore, increasing plate constraint (using clamped edge boundary condition) will increase the value of NF of vibrations.

The current analysis is suitable for all-composite symmetric LHSP. However, for isotropic or non-symmetric materials, minor modifications might be adopted.

The application of simplified HSDT to the analysis of LHSP is one of the important values of this research. The other is the successful and complete dynamic analysis of all-composite LHSP.

This article aims to present a direct solution to handle linear constraints in finite element (FE) analysis without penalties or the Lagrange multipliers introduced.

First, the system of linear equations corresponding to the linear constraints is solved for the leading variables in terms of the free variables and the constants. Then, the reduced system of equilibrium equations with respect to the free variables is derived from the finite-dimensional virtual work equation. Finally, the algorithm is designed.

The proposed procedure is promising in three typical cases: (1) to enforce displacement constraints in any direction; (2) to implement local refinements by allowing hanging nodes from element subdivision and (3) to treat non-matching grids of distinct parts of the problem domain. The procedure is general and suitable for 3D non-linear analyses.

The algorithm is fitted only to the Galerkin-based numerical methods.

The proposed procedure does not need Lagrange multipliers or penalties. The tangential stiffness matrix of the reduced system of equilibrium equations reserves positive definiteness and symmetry. Besides, many contemporary Galerkin-based numerical methods need to tackle the enforcement of the essential conditions, whose weak forms reduce to linear constraints. As a result, the proposed procedure is quite promising.

Accurate and automatic clothing pattern making is very important in personalized clothing customization and virtual fitting room applications. Clothing pattern generating as well as virtual clothing simulation is an attractive research issue both in clothing industry and computer graphics.

Physics-based method is an effective way to model dynamic process and generate realistic clothing animation. Due to conceptual simplicity and computational speed, mass-spring model is frequently used to simulate deformable and soft objects follow the natural physical rules. We present a physics-based clothing pattern generating framework by using scanned human body model. After giving a scanned human body model, first, we extract feature points, planes and curves on the 3D model by geometric analysis, and then, we construct a remeshed surface which has been formatted to connected quad meshes. Second, for each clothing piece in 3D, we construct a mass-spring model with same topological structures, and conduct a typical time integration algorithm to the mass-spring model. Finally, we get the convergent clothing pieces in 2D of all clothing parts, and we reconnected parts which are adjacent on 3D model to generate the basic clothing pattern.

The results show that the presented method is a feasible way for clothing pattern generating by use of scanned human body model.

The main contribution of this work is twofold: one is the geometric algorithm to scanned human body model, which is specially conducted for clothing pattern design to extract feature points, planes and curves. This is the crucial base for suit clothing pattern generating. Another is the physics-based pattern generating algorithm which flattens the 3D shape to 2D shape of cloth pattern pieces.