An F-bar aided edge-based smoothed finite element methods with four-node tetrahedral elements (F-barES-FEM-T4) in large deformation explicit dynamics for nearly incompressible materials is proposed. Because F-barES-FEM-T4 is based on on the standard displacement-based variational formulations, it can be easily applied to dynamic analysis by adding the term of inertia to the static equilibrium equation. Some example analyses of explicit dynamics for simple and complex shapes reveal that the proposed method can suppress volumetric locking and pressure oscillation in dynamic explicit cases as well as in static implicit cases. Meanwhile, these analyses also reveal that the proposed method causes energy divergence in long-term analyses. As the speed of energy divergence can be suppressed by increasing the smoothing level, the proposed method with sufficient smoothing presents good results in relatively short-term analyses.

View full abstract

For large-scale numerical simulations on supercomputers, data transfer and storage present significant efficiency and productivity issues. Therefore, the jointed hierarchical precision compression number-data format (JHPCN-DF) technique was proposed for efficient visualization and analysis of plasma particle-in-cell simulation data. It is also available for lossless and lossy compression with user-defined errors. We implement a lossy compression method of JHPCN-DF in finite element code and evaluate the compression effectiveness and compression data accuracy in linear static and dynamic structural analyses. Our technique achieves the required accuracy, even for dynamic problems, and provides a significant increase in compression performance for variable datasets.

View full abstract

In order to educate teenager internet literacy on social network service, we have developed a Problem Solving Environment to evaluate the literacy-level of their messages on twitter for their teachers and them. We propose a method the system provides effective recognition for their risks. And we adapt the Naive Bayes classifier to evaluation for tweets on Twitter based on pattern-based classifier. In this result, the classification accuracy for word patterns increases from 39.6-57.6% to 68.0-79.9% using Naive Bayes classifier on a set of 3000 training data sets, and users obtain internet literacy skills base on this system.

View full abstract

To enable the prediction of the mechanical degradation of solid oxide fuel cell (SOFC) during operation, we have developed an analysis system for the electrochemical and mechanical coupling phenomena by incorporating general-purpose finite element analysis (FEA) software along with its pre- and post-processing functions and conducted validation analyses to examine its quantitative performance in reproducing the measurement data. The dependencies of material properties on both temperature and oxygen’s chemical potentials are empirically obtained and stored in separate spreadsheets and read into the system so that the actual operation environment can be realized. The thermal and reduction strains calculated in the transient analysis for the potentials are also stored in spreadsheets beforehand and used in the stress analyses to make quantitative evaluation of creep deformations. Through a numerical example, the mathematical model implemented into the developed system has been initially validated in a quantitative way.

View full abstract

A new nonlinear film resistance model and a new nonlinear film growth model in electrodeposition (ED) coating simulation are presented. The new resistance model considers the electric resistance of turbid paint on cathode before the film deposition. The new growth model considers history dependency in the film growth rate varied with current density on cathode when the initial film deposits. The parameters of the new models are identified with a series of simple one-plate ED test data. Some finite element analyses show that the results with the new models are more accurate compared to those with the conventional models.

View full abstract

Manufacturers must provide products that are robust under use by various types of customers. In this study, the quality engineering was used to optimize the shape of the crank-case considering the influence of the error. The oil behavior in the crank-case was calculated by the MPS method. The optimization was conducted using the orthogonal array L18. The design parameters were related to the shape of the crank-case, and the error was the viscosity of the oil and the revolving speed of the crank-shaft. Furthermore, the crank-case was prototyped and the oil temperature in crank-case was measured to validate the calculation result.

View full abstract

Accurate numerical integration technique for a discontinuous integrand by the Delaunay decomposition has been proposed to calculate interactive terms in finite element mesh superposition method (FEMS). In this paper, the high accuracy numerical integration technique for high order local elements with curved surface is proposed. High order elements are divided into segments only with plane boundaries using some additional vertices to deal with hollow domains adequately. The previous proposed technique that is based on grouping of vertices to obtain a tetrahedralization is applied to all segments. In the first example, the performance of the proposed technique is estimated in 2-D analysis. In the second example, the proposed technique is applied to 3-D zooming analysis. In the 3rd example, smooth microstructures is add in VOXEL model using FEMS with the proposed technique.

View full abstract

The quasi-three-dimensional XFEM is applied to damage propagation analyses of CFRP (Carbon Fiber Reinforced Plastics) laminate. Six-node triangular interface element and six-node pentahedral continuum element enriched with Heaviside function are employed to model delamination and matrix crack, respectively. Bi-linear type cohesive zone model is introduced between delamination and/or matrix crack and then implicit static or explicit dynamic method is utilized to solve system equations considering materially nonlinearity. Code verification was performed through DCB, ENF and TCT test specimen analyses and damage propagation analyses of NHT (No Hole Tension) and OHT (Open Hole Tension) test specimens were validated by comparison with experiment results. In addition, computation conditions for explicit dynamic analysis including mass scaling and energy balance, Zig-zag CZM, and threshold parameter of crack shape definition were examined. It was shown that implicit static method using Zig-zag CZM is the most efficient for damage propagation analysis using the proposed method.

View full abstract

There is a trade-off problem between the distance and the repeatability in golf. Our goal is to obtain the pareto solution by formulating the problem as a multi-objective optimization, and to develop the optimal golf shaft for the player. Club speed and body burden were used as indexes for the distance and the repeatability. Finite element method and multibody dynamics were used for the calculation of the club speed and the body burden. Response surface method was applied to get a pareto solution. As an analytical result, it was found out that the optimized shaft increased club speed by 2.2 m/s and decreased the body burden by around 7% compared to the control shaft. The optimized shaft was manufactured and conducted some experiments for validation. As average data of 10 shots, the distance increased by 8 yds and repeatability was improved by 16 yds with the optimal shaft.

View full abstract

The purpose of this study is to develop a snow model emphasizing the splashing phenomenon and the impact force acting on the train surface to ensure safe driving. Snow is assumed to be a Bingham fluid due to its peculiar shear thinning. A Bingham numerical model is proposed by using the Moving Particle Semi-Implicit (MPS) method and tested for collision between a train and snow. As well, verification of this special non-Newtonian fluid is carried out and good agreement is obtained by comparing the MPS numerical result with a reference solution. Furthermore, a collision simulation between train and snow is implemented in which a natural splashing phenomenon occurs by introducing the present Bingham snow model. Finally, the time history of pressure in the foremost area of the train head is analyzed, which can provide a reference for safe train driving.

View full abstract

This study presents a numerical method of snow avalanches. The method is based on a stabilized finite element method in SUPG/PSPG formulations with unstructublack mesh. The flow property of the snow is represented by a Bingham type fluid model; the shear strength is described using a Coulomb’s yield criterion. In order to verify the method, simulations of a model experiment were performed, and simulated run-out distance and impact force were compablack with the experimental results. The method was finally applied to an actual snow avalanche that took place at the foot of the Mount Gassan in 2011. Simulated results reproduced the actual snow avalanche accurately, and it was found the method can represent flow behaviors of snow avalanches even in the complex calculation conditions, such as those containing complex geometry, retaining walls and vegetation.

View full abstract

Recently, numerous researches on loudspeaker vibration have been conducted. However, little attention has been to the damping effect, which is of great importance concerning air viscosity in the vibration analysis of cone loudspeaker. Consequently, without consideration of the air viscosity, the accuracy in the estimation of the sound pressure frequency characteristic is deteriorated. To solve this problem, this research considers the viscosity of air in the narrow space around voice coil, and captures the air surrounding the voice coil. In this way, the high sound pressure frequency characteristics with high accuracy could be obtained. Furthermore, by considering the viscosity of the air, the vibration characteristics of the paper cone and the edge of the vibration system part can be reproduced.

View full abstract

The mass flow rate and the spray characteristics of an automotive fuel injector are substantially influenced by cavitation. It is therefore very useful to be able to perform the fuel flow analysis in nozzles using a numerical analysis method. To our knowledge, there is only one cavitation model until now that is proposed for Moving Particle Semi-implicit (MPS) method. This model, however, is relatively complicated and unstable due to using particles with a variable particle diameter. This paper describes a more simplified and stable cavitation model: the calculation algorithm for the growth of bubbles and the contraction and collapse of bubbles is not used. The validity of the proposed model is examined by the fuel flow analysis in a Schmitt’s nozzle. The simulation results for a mass flow rate are in good agreement with the experiment results. The images showing the flow state are generally consistent with them.

View full abstract

Ultimate strength evaluation for cracked plates is examined. A rectangular plate including a through-crack under uni-axial compression load is assumed. In analyzing the cracked plate, shell finite element (FE) modeling and shell-solid FE modeling are respectively adopted. Contact condition is included in the shell-solid modeling to introduce the crack-closure effect. Numerical examples are presented for analyzing accuracy in the FE modeling and investigating buckling and post-buckling behaviors of the cracked plates by changing position and length of the through-crack.

View full abstract

In this paper, we present an interpolating a ground surface method for a Digital Elevation Model (DEM) data based on the Heaviside B-Spline interpolation. To interpolate the DEM data, we modify a Heaviside function used for Heaviside B-Spline interpolation. By calculating the Heaviside function using the elevation value contained in the DEM data, it is possible to obtain the Spline function and an interpolating surface without explicitly giving boundary of sea water bodies, inner waters and void regions. Therefore, our proposed Heaviside function enables interpolating a ground surface including sea water bodies, inland waters or void regions in a DEM data by using a least square method. We show the interpolating results for several DEM data which include elevation discontinuous.

View full abstract

Domain-decomposition method (DDM) is an efficient parallelization scheme for the fintie element method. To optimize the performance of the DDM subdomain local solver for multi-core CPU architecture, we propose iterative solver-based approaches, which enable a DDM-based ”on-cache” iterative solver. A subdomain local solver using SSOR pre-conditioned CG with Eisenstat’s trick, which is OpenMP parallelized, is not only more efficient in memory usage but also comparable or even faster in computational time than the existing direct solver-based approach on multi-core PC clusters.

View full abstract

In this paper, a computationally implementable and feasible two-mode asymptotic bifurcation theory is proposed for structural buckling analysis. Using this theory, snap-through, asymmetric bifurcation, unstable and stable symmetric bifurcation might be reliably diagnosed by solely solving two simultaneous polygonal equations which can be derived by exact first and second order directional derivatives of singular tangent stiffness matrices with respect to nodal degree-of-freedoms along critical and non-critical modes. The error-free computation of the derivatives is accomplished by using hyper-dual numbers. Graphical tool such as Matlab is useful to visualize the existence of the solutions of the two resulting simultaneous polygonal equations and number of possible real roots can well predict post-buckling behavior. Numerical examples of nonlinear finite element stability analysis verify the robustness of the proposed two-mode asymptotic bifurcation theory.

View full abstract