Impact of finite element idealisation on the prediction of welded fuselage stiffened panel buckling

Lap joints are widely used in the manufacture of stiffened panels and influence local panel sub-component stability, defining buckling unit dimensions and boundary conditions. Using the Finite Element method it is possible to model joints in great detail and predict panel buckling behaviour with accuracy. However, when modelling large panel structures such detailed analysis becomes computationally expensive. Moreover, the impact of local behaviour on global panel performance may reduce as the scale of the modelled structure increases. Thus this study presents coupled computational and experimental analysis, aimed at developing relationships between modelling fidelity and the size of the modelled structure, when the global static load to cause initial buckling is the required analysis output. Small, medium and large specimens representing welded lap-joined fuselage panel structure are examined. Two element types, shell and solid-shell, are employed to model each specimen, highlighting the impact of idealisation on the prediction of welded stiffened panel initial skin buckling.

An adaptive finite element water column model using the Mellor-Yamada level 2.5 turbulence closure scheme

A one-dimensional water column model using the Mellor and Yamada level 2.5 parameterization of vertical turbulent fluxes is presented. The model equations are discretized with a mixed finite element scheme. Details of the finite element discrete equations are given and adaptive mesh refinement strategies are presented. The refinement criterion is an "a posteriori" error estimator based on stratification, shear and distance to surface. The model performances are assessed by studying the stress driven penetration of a turbulent layer into a stratified fluid. This example illustrates the ability of the presented model to follow some internal structures of the flow and paves the way for truly generalized vertical coordinates. (c) 2005 Elsevier Ltd. All rights reserved.

A finite element-spherical harmonics model for radiative transfer in inhomogeneous clouds. Part II. some applications

EVENT has been used to examine the effects of 3D cloud structure, distribution, and inhomogeneity on the scattering of visible solar radiation and the resulting 3D radiation field. Large eddy simulation and aircraft measurements are used to create realistic cloud fields which are continuous or broken with smooth or uneven tops. The values, patterns and variance in the resulting downwelling and upwelling radiation from incident visible solar radiation at different angles are then examined and compared to measurements. The results from EVENT confirm that 3D cloud structure is important in determining the visible radiation field, and that these results are strongly influenced by the solar zenith angle. The results match those from other models using visible solar radiation, and are supported by aircraft measurements of visible radiation, providing confidence in the new model.

Application of radial return mapping algorithm for finite strain elastoplasticity in a hybrid finite element and particle-in-cell model

The radial return mapping algorithm within the computational context of a hybrid Finite Element and Particle-In-Cell (FE/PIC) method is constructed to allow a fluid flow FE/PIC code to be applied solid mechanic problems with large displacements and large deformations. The FE/PIC method retains the robustness of an Eulerian mesh and enables tracking of material deformation by a set of Lagrangian particles or material points. In the FE/PIC approach the particle velocities are interpolated from nodal velocities and then the particle position is updated using a suitable integration scheme, such as the 4th order Runge-Kutta scheme[1]. The strain increments are obtained from gradients of the nodal velocities at the material point positions, which are then used to evaluate the stress increment and update history variables. To obtain the stress increment from the strain increment, the nonlinear constitutive equations are solved in an incremental iterative integration scheme based on a radial return mapping algorithm[2]. A plane stress extension of a rectangular shape J2 elastoplastic material with isotropic, kinematic and combined hardening is performed as an example and for validation of the enhanced FE/PIC method. It is shown that the method is suitable for analysis of problems in crystal plasticity and metal forming. The method is specifically suitable for simulation of neighbouring microstructural phases with different constitutive equations in a multiscale material modelling framework.

Finite element modeling and analysis of thrombosis in middle cerebral artery

Thrombotic stroke, which is caused by blood clot in the cerebral artery, is a major source of increased mortality and morbidity. Considering as efficient and fastest methods, mathematical approaches have gained significant importance for analyzing and understanding the biological events like thrombosis. This paper presents a computational model to analyze the effects of thrombosis using the theory of coupled fluid dynamics-structure interaction. The finite element method is used for the modeling of thrombosis (blood clot) of different stages in the middle cerebral artery with physiological compliance. The developed model is used to investigate the consequences that occur due to the various sizes of clots in the artery in the form of blood flow velocity, blood pressure, and artery wall stress. Such numerical assessment will facilitate better understanding of the biophysical process in case of thrombosis and thus would support medical practitioners to take faster curing steps.

A finite element free convection model for the side wall heated cavity

This paper presents a finite element numerical solution of free convection in a cavity with side walls maintained at constant but different temperatures. The predictions from the model and the method of solution were validated by comparison with the 'bench mark' solution and Vahl Davis' results and good agreement was found. The present model was used to obtain additional results over a wide range of Rayleigh number (10(3)-10(6)) and L/H ratios varying from 0.1 to 1.0. The predicted stream function patterns, temperature and velocity profiles as well as the mean Nusselt number were presented and discussed. (C) 2000 Elsevier B.V. Ltd. All rights reserved.

A coupling technique for non-matching finite element meshes

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)