Search results

The response of the Pine Flat dam–water–foundation rock system is studied by a new described approach (i.e. FE-(FE-TE)-FE). The initial part of study is focused on the time harmonic analysis. In this part, it is possible to compare the transfer functions against corresponding responses obtained by the FE-(FE-HE)-FE approach (referred to as exact method which employs a rigorous fluid hyper-element). Subsequently, the transient analysis is carried out. In that part, it is only possible to compare the results for low and high normalized reservoir length cases. Therefore, the sensitivity of results is controlled due to normalized reservoir length values.

In the present study, dynamic analysis of a typical concrete gravity dam–water–foundation rock system is formulated by the FE-(FE-TE)-FE approach. In this technique, dam and foundation rock are discretized by plane solid finite elements while, water domain near-field region is discretized by plane fluid finite elements. Moreover, the H-W (i.e. Hagstrom–Warburton) high-order condition is imposed at the reservoir truncation boundary. This task is formulated by employing a truncation element at that boundary. It is emphasized that reservoir far-field is excluded from the discretized model.

High orders of H-W condition, such as O5-5 considered herein, generate highly accurate responses for both possible excitations under both types of full reflective and absorptive reservoir bottom conditions. It is such that transfer functions are hardly distinguishable from corresponding exact responses obtained through the FE-(FE-HE)-FE approach in time harmonic analyses. This is controlled for both low and high normalized reservoir length cases (L/H = 1 and 3). Moreover, it can be claimed that transient analysis leads practically to exact results (in numerical sense) when one is employing high order H-W truncation element. In other words, the results are not sensitive to reservoir normalized length under these circumstances.

Dynamic analysis of concrete gravity dam–water–foundation rock systems is formulated by a new method. The salient aspect of the technique is that it utilizes H-W high-order condition at the truncation boundary. The method is discussed for all types of excitation and reservoir bottom conditions.

The vibration of the rails is a significant source of railway rolling noise, often forming the dominant component of noise in the important frequency region between 400 and 2000 Hz. The purpose of the paper is to investigate the influence of the ground profile and the presence of the train body on the sound radiation from the rail.

Two-dimensional boundary element calculations are used, in which the rail vibration is the source. The ground profile and various different shapes of train body are introduced in the model, and results are observed in terms of sound power and sound pressure. Comparisons are also made with vibro-acoustic measurements performed with and without a train present.

The sound radiated by the rail in the absence of the train body is strongly attenuated by shielding due to the ballast shoulder. When the train body is present, the sound from the vertical rail motion is reflected back down toward the track where it is partly absorbed by the ballast. Nevertheless, the sound pressure at the trackside is increased by typically 0–5 dB. For the lateral vibration of the rail, the effects are much smaller. Once the sound power is known, the sound pressure with the train present can be approximated reasonably well with simple line source directivities.

Numerical models used to predict the sound radiation from railway rails have generally neglected the influence of the ground profile and reflections from the underside of the train body on the sound power and directivity of the rail. These effects are studied in a systematic way including comparisons with measurements.

This study investigates the impact of three parameters such as: number of LED chips, pitch and LED power on the junction temperature of LEDs using a best heat sink configuration selected according to a lower temperature. This study provides valuable insights into how to design LED arrays with lower junction temperatures.

To determine the best configuration of a heat sink, a numerical study was conducted in Comsol Multiphysics on 10 different configurations. The configuration with the lowest junction temperature was selected for further analysis. The number of LED chips, pitch and LED power were then varied to determine the optimal configuration for this heat sink. A general equation for the average LED temperature as a function of these three factors was derived using Minitab software.

Among 10 configurations of the rectangular heat sink, we deduce that the best configuration corresponds to the first design having 1 mm of width, 0.5 mm of height and 45 mm of length. The average temperature for this design is 50.5 C. For the power of LED equal to 50 W–200 W, the average temperature of this LED drops when the number of LED chips reduces and the pitch size decreases. Indeed, the best array-LED corresponds to 64 LED chips and a pitch size of 0.5 mm. In addition, a generalization equation for average temperature is determined as a function of the number of LED chips, pitch and power of LED which are key factors for reducing the Junction temperature.

The study is original in its focus on three factors that have not been studied together in previous research. A numerical simulation method is used to investigate the impact of the three factors, which is more accurate and reliable than experimental methods. The study considers a wide range of values for the three factors, which allows for a more comprehensive understanding of their impact. It derives a general equation for the average temperature of the LED, which can be used to design LED arrays with desired junction temperatures.

This work compares the performance of the three boundary element techniques for solving Helmholtz problems: dual reciprocity, multiple reciprocity and direct interpolation. All techniques transform domain integrals into boundary integrals, despite using different principles to reach this purpose.

Comparisons here performed include the solution of eigenvalue and response by frequency scanning, analyzing many features that are not comprehensively discussed in the literature, as follows: the type of boundary conditions, suitable number of degrees of freedom, modal content, number of primitives in the multiple reciprocity method (MRM) and the requirement of internal interpolation points in techniques that use radial basis functions as dual reciprocity and direct interpolation.

Among the other aspects, this work can conclude that the solution of the eigenvalue and response problems confirmed the reasonable accuracy of the dual reciprocity boundary element method (DRBEM) only for the calculation of the first natural frequencies. Concerning the direct interpolation boundary element method (DIBEM), its interpolation characteristic allows more accessibility for solving more elaborate problems. Despite requiring a greater number of interpolating internal points, the DIBEM has presented higher-quality results for the eigenvalue and response problems. The MRM results were satisfactory in terms of accuracy just for the low range of frequencies; however, the neglected higher-order primitives impact the accuracy of the dynamic response as a whole.

There are safe alternatives for solving engineering stationary dynamic problems using the boundary element method (BEM), but there are no suitable comparisons between these different techniques. This paper presents the particularities and detailed comparisons approaching the accuracy of the three important BEM techniques, aiming at response and frequency evaluation, which are not found in the specialized literature.

This paper aims to investigate the performance of two novel numerical methods, the face-based smoothed finite element method (FS-FEM) and the edge-based smoothed finite element method (ES-FEM), which employ linear tetrahedral elements, for the purpose of strength assessment of a high-speed train hollow axle.

The calculation of stress for the wheelset, comprising an axle and two wheels, is facilitated through the application of the European axle strength design standard. This standard assists in the implementation of loading and boundary conditions and is exemplified by the typical CRH2 high-speed train wheelset. To evaluate the performance of these two methods, a hollow cylinder cantilever beam is first used as a benchmark to compare the present methods with other existing methods. Then, the strength analysis of a real wheelset model with a hollow axle is performed using different numerical methods.

The results of deflection and stress show that FS-FEM and ES-FEM offer higher accuracy and better convergence than FEM using linear tetrahedral elements. ES-FEM exhibits a superior performance to that of FS-FEM using linear tetrahedral elements, showing accuracy and convergence close to FEM using hexahedral elements.

This study channels the novel methods (FS-FEM and ES-FEM) in the static stress analysis of a railway wheelset. Based on the careful testing of FS-FEM and ES-FEM, both methods hold promise as more efficient tools for the strength analysis of complex railway structures.

Hexahedral meshing is one of the most important steps in performing an accurate simulation using the finite element analysis (FEA). However, the current hexahedral meshing method in the industry is a nonautomatic and inefficient method, i.e. manually decomposing the model into suitable blocks and obtaining the hexahedral mesh from these blocks by mapping or sweeping algorithms. The purpose of this paper is to propose an almost automatic decomposition algorithm based on the 3D frame field and model features to replace the traditional time-consuming and laborious manual decomposition method.

The proposed algorithm is based on the 3D frame field and features, where features are used to construct feature-cutting surfaces and the 3D frame field is used to construct singular-cutting surfaces. The feature-cutting surfaces constructed from concave features first reduce the complexity of the model and decompose it into some coarse blocks. Then, an improved 3D frame field algorithm is performed on these coarse blocks to extract the singular structure and construct singular-cutting surfaces to further decompose the coarse blocks. In most modeling examples, the proposed algorithm uses both types of cutting surfaces to decompose models fully automatically. In a few examples with special requirements for hexahedral meshes, the algorithm requires manual input of some user-defined cutting surfaces and constructs different singular-cutting surfaces to ensure the effectiveness of the decomposition.

Benefiting from the feature decomposition and the 3D frame field algorithm, the output blocks of the proposed algorithm have no inner singular structure and are suitable for the mapping or sweeping algorithm. The introduction of internal constraints makes 3D frame field generation more robust in this paper, and it can automatically correct some invalid 3–5 singular structures. In a few examples with special requirements, the proposed algorithm successfully generates valid blocks even though the singular structure of the model is modified by user-defined cutting surfaces.

The proposed algorithm takes the advantage of feature decomposition and the 3D frame field to generate suitable blocks for a mapping or sweeping algorithm, which saves a lot of simulation time and requires less experience. The user-defined cutting surfaces enable the creation of special hexahedral meshes, which was difficult with previous algorithms. An improved 3D frame field generation method is proposed to correct some invalid singular structures and improve the robustness of the previous methods.

This paper aims to focus on the finite element method in the frequency domain (FD-FEM) for the transient electric field in the non-sinusoidal steady state under the non-sinusoidal periodic voltage excitation.

Firstly, the boundary value problem of the transient electric field in the frequency domain is described, and the finite element equation of the FD-FEM is derived by Galerkin’s method. Secondly, the constrained electric field equation on the boundary in the frequency domain (FD-CEFEB) is also derived, which can solve the electric field intensity on the boundary and the dielectric interface with high accuracy. Thirdly, the calculation procedures of the FD-FEM with FD-CEFEB are introduced in detail. Finally, a numerical example of the press-packed insulated gate bipolar transistor under the working condition of the repetitive turn-on and turn-off is given.

The FD-CEFEB improves numerical accuracy of electric field intensity on the boundary and interfacial charge density, which can be achieved by modifying the existing FD-FEMs’ code in appropriate steps. Moreover, the proposed FD-FEM and the FD-CEFEB will only increase calculation costs by a little compared with the traditional FD-FEMs.

The FD-CEFEB can directly solve the electric field intensity on the boundary and the dielectric interface with high accuracy. This paper provides a new FD-FEM for the transient electric field in the non-sinusoidal steady state with high accuracy, which is suitable for combined insulation structure with a long time constant.

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.

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.

This paper aims to provide SIF calculation method for engineering application.

In this paper, the stress intensity factors (SIFs) calculation method is applied to the anisotropic Ni-based single crystal film cooling holes (FCHs) structure.

Based on contour integral, the anisotropic SIFs analysis finite element method (FEM) in Ni-based single crystal is proposed. The applicability and mesh independence of the method is assessed by comparing the calculated SIFs using mode of plate with an edge crack. Anisotropic SIFs can be calculated with excellent accuracy using the finite element contour integral approach. Then, the effect of crystal orientation and FCHs interference on the anisotropic SIFs is clarified. The SIFs of FCH edge crack in the [011] orientated Ni-based single crystal increases faster than the other two orientations. And the SIF of horizontal interference FCHs edge crack is also larger than that of the inclined interference one.

The SIFs of the FCH edge crack in the turbine air-cooled blade are innovatively computed using the sub-model method. Both the Mode I and II SIFs of FCHs edge crack in blade increase with crack growing.