COLLECTIONS

JOURNALS

Documents published in Scipedia

• Metal forming analysis by fourier series expansion and further uses of pseudo-concentrations

  H. Antúnez, S. Idelsohn, E. Dvorkin

  Computers and Structures (1992). Vol. 44 (1-2), pp. 435-451

  
  

  Abstract

  An analysis of metal forming processes with axial geometry is performed by an expansion in Fourier series on the circumferential direction. For that purpose, an incremental form is presented on the basis of the flow formulation. For every increment of the non-null imposed boundary conditions an algebraic system is formed in which the non-linearity due to the visco-plastic properties of the material is placed as part of the load term. In this way, the global stiffness matrix is block diagonal, the coupling due to material non-linearity being included on the right-hand side. Free surfaces are treated as an interface inside the elements, as an extension of the pseudo-concentration method.

  Abstract
  An analysis of metal forming processes with axial geometry is performed by an expansion in Fourier series on the circumferential direction. For that purpose, an incremental [...]

  • 7
  •  read

• An effective automatic incremental/iterative method for static nonlinear structural analysis

  J. Jeusette, G. Laschet, S. Idelsohn

  Computers and Structures (1989). Vol. 32 (1), pp. 125-135

  
  

  Abstract

  In the prebuckling range of the complete nonlinear response of a structure submitted to destabilizing loads, a linearized stability analysis is generally an interesting tool. Used in conjunction with an arc-length algorithm, incremental bifurcation analyses are able to give optimal definition for the increment predictor size. When the associated eigenvalue problem is solved by minimizing the finite element discretized form of the Rayleigh quotient via a preconditioned conjugate gradient technique, it can be very easily inserted in an automatic nonlinear equations solver. Such a procedure is presented in this paper after a review of existing arc-length schemes and step-size control methods. Moreover, for tracing the postbuckling range of the nonlinear response automatically, a new recurrence formula is proposed to control the step-size. The well behaviour of this algorithm is illustrated in two classical structural applications. 

  Abstract
  In the prebuckling range of the complete nonlinear response of a structure submitted to destabilizing loads, a linearized stability analysis is generally an interesting tool. [...]

  • 8
  •  read

• A load-dependent basis for reduced nonlinear structural dynamics

  S. Idelsohn, A. Cardona

  Computers and Structures (1985). Vol. 20 (1-3), pp. 203-210

  
  

  Abstract

  A computational algorithm for predicting the dynamical response of a nonlinear structure by means of a reduction scheme is described. In it, the nonlinear system of ordinary differential equations obtained from the finite element discretization is reduced, by employing a Rayleigh-Ritz technique. A new criterion for the computation of the basis vectors is proposed. It is based on a sequence of basis vectors proposed by Wilson et al. for linear dynamics problems, which was augmented by adding derivatives of these vectors with respect to the generalized displacement amplitudes. The later vectors allow the treatment of nonlinear problems.

  Abstract
  A computational algorithm for predicting the dynamical response of a nonlinear structure by means of a reduction scheme is described. In it, the nonlinear system of ordinary [...]

  • 19
  •  read

• Computational strategies for the solution of large nonlinear problems via quasi-newton methods

  M. Geradin, S. Idelsohn, M. Hogge

  Computers and Structures (1981). Vol. 13 (1-3), pp. 73-81

  
  

  Abstract

  The usefulness of quasi-Newton methods for the solution of nonlinear systems of equations is demonstrated. After a review of the Newton iterative method, several quasi-Newton updates are presented and tested. Special attention is devoted to the solution of large sparse systems of equations such as those issued from spatial discretization of continua by finite elements. The numerical examples presented comprise static and dynamic analyses of geometrical, material and combined nonlinear structural problems and a model fluid flow problem with different levels of nonlinearity. All the results are assorted with a complete discussion of the different methods used, of the convergence rates and of the associated computer costs. From the present studies, it can be concluded that computational costs for the solution of large nonlinear systems of equations can be reduced drastically by using convenient quasi-Newton updates or by adequate combined Newton/quasi-Newton strategies.

  Abstract
  The usefulness of quasi-Newton methods for the solution of nonlinear systems of equations is demonstrated. After a review of the Newton iterative method, several quasi-Newton [...]

  • 10
  •  read

• Failure internal pressure of spherical steel containments

  G. Sarmiento, S. Idelsohn, A. Cardona, V. Sonzogni

  Nuclear Engineering and Desing (1985). Vol. 90 (2), pp. 209-222

  
  

  Abstract

  An application of the British CEGB's R6 Failure Assessment Approach to the determination of failure internal pressure of nuclear power plant spherical steel containments is presented. The presence of hypothetical cracks both in the base metal and in the welding material of the containment, with geometrical idealizations according to the ASME Boiler and Pressure Vessel Code (Section XI), was taken into account in order to analyze the sensitivity of the failure assessment with the values of the material fracture properties. Calculations of the elastoplastic collapse load have been performed by means of the Finite Element System SAMCEF. The clean axisymmetric shell (neglecting the influence of nozzles and minor irregularities) and two major penetrations (personnel and emergency locks) have been taken separately into account. Large-strain elastoplastic behaviour of the material was considered in the Code, using lower bounds of true stress—true strain relations obtained by testing a collection of tensile specimens. Assuming the presence of cracks in non-perturbed regions, the reserve factor for test pressure and the failure internal pressure have been determined as a function of the flaw depth.

  Abstract
  An application of the British CEGB's R6 Failure Assessment Approach to the determination of failure internal pressure of nuclear power plant spherical steel containments [...]

  • 7
  •  read

• Nonlinear structural dynamics via Newton and quasi-Newton methods

  M. Geradin, S. Idelsohn, M. Hogge

  Nuclear Engineering and Design (1980). Vol. 58 (3), pp. 339-348;

  
  

  Abstract

  This paper is an attempt to compare Newton and quasi-Newton methods in nonlinear structural dynamics. After a review of the classical iterative methods, several quasi-Newton updates are presented and tested. Special attention is devoted to the solution of large sparse systems for which two original procedures are described: a substructure correction and a vectorial correction. The numerical examples presented include the dynamic analyses of geometrical, material and combined nonlinearities. All the results are assorted with a complete discussion of the different methods used, of the convergence rates and of the associated computer costs. From the present results, Newton's methods appear to exhibit the best convergence rates when an efficient computational strategy is adopted. Nevertheless computational costs for the solution of large systems can be reduced drastically by using convenient quasi-Newton updates.

  Abstract
  This paper is an attempt to compare Newton and quasi-Newton methods in nonlinear structural dynamics. After a review of the classical iterative methods, several quasi-Newton [...]

  • 21
  •  read

• OpenCL‐based implementation of an unstructured edge‐based finite element convection‐diffusion solver on graphics hardware

  F. Mossaiby, R. Rossi, P. Dadvand, S. Idelsohn

  Int. J. Numer. Meth. Engng. (2012). Vol. 89 (13), pp. 1635-1651

  
  

  Abstract

  The solution of problems in computational fluid dynamics (CFD) represents a classical field for the application of advanced numerical methods. Many different approaches were developed over the years to address CFD applications. Good examples are finite volumes, finite differences (FD), and finite elements (FE) but also newer approaches such as the lattice‐Boltzmann (LB), smooth particle hydrodynamics or the particle finite element method. FD and LB methods on regular grids are known to be superior in terms of raw computing speed, but using such regular discretization represents an important limitation in dealing with complex geometries. Here, we concentrate on unstructured approaches which are less common in the GPU world. We employ a nonstandard FE approach which leverages an optimized edge‐based data structure allowing a highly parallel implementation. Such technique is applied to the ‘convection‐diffusion’ problem, which is often considered as a first step towards CFD because of similarities to the nonconservative form of the Navier–Stokes equations. In this regard, an existing highly optimized parallel OpenMP solver is ported to graphics hardware based on the OpenCL platform. The optimizations performed are discussed in detail. A number of benchmarks prove that the GPU‐accelerated OpenCL code consistently outperforms the OpenMP version.

  Abstract
  The solution of problems in computational fluid dynamics (CFD) represents a classical field for the application of advanced numerical methods. Many different approaches were [...]

  • 9
  •  read

• All‐hexahedral element meshing: automatic elimination of self‐intersecting dual lines

  N. Calvo, S. Idelsohn

  Int. J. Numer. Meth. Engng. (2002). Vol. 55 (12), pp. 1439-1449

  
  

  Abstract

  There has been some degree of success in all‐hexahedral meshing. Standard methods start with the object geometry defined by means of an all‐quadrilateral mesh, followed by the use of the combinatorial dual to the mesh in order to define the internal connectivities among elements. For all of the known methods using the dual concept, it is necessary to first prevent or eliminate self‐intersecting (SI) dual lines of the given quadrilateral mesh. The relevant features of SI lines are studied, giving a method to remove them, which avoids deforming the original geometry. Some examples of resulting meshes are shown where the current meshing method has been successfully applied. 

  Abstract
  There has been some degree of success in all‐hexahedral meshing. Standard methods start with the object geometry defined by means of an all‐quadrilateral mesh, followed [...]

  • 3
  •  read

• All‐hexahedral mesh smoothing with a node‐based measure of quality

  N. Calvo, S. Idelsohn

  Int. J. Numer. Meth. Engng. (2001). Vol. 50 (8), pp. 1957-1967

  
  

  Abstract

  This research work deals with the analysis and test of a normalized‐Jacobian metric used as a measure of the quality of all‐hexahedral meshes. Instead of element qualities, a measure of node quality was chosen. The chosen metric is a bound for deviation from orthogonality of faces and dihedral angles. We outline the main steps and algorithms of a program that is successful in improving the quality of initially invalid meshes to acceptable levels. For node movements, the program relies on a combination of gradient‐driven and simulated annealing techniques. Some examples of the results and speed are also shown.

  Abstract
  This research work deals with the analysis and test of a normalized‐Jacobian metric used as a measure of the quality of all‐hexahedral meshes. Instead of element qualities, [...]

  • 7
  •  read

• Non‐reflective planar boundary condition based on Gauss filter

  R. Bonet Chaple, N. Nigro, M. Storti, S. Idelsohn

  Int. J. Numer. Meth. Engng. (2000). Vol. 47 (5), pp. 969-983

  
  

  Abstract

  A non‐reflecting boundary condition based on the Gauss filter is employed for the determination of scattered potential governed by the equation. A filtering layer is used for closing infinite domain calculations. An expression for the reflection coefficient is derived and an optimal filtering layer is designed. Numerical results validate the performance of this method for unbounded wave guide problems.

  Abstract
  A non‐reflecting boundary condition based on the Gauss filter is employed for the determination of scattered potential governed by the equation. A filtering layer is used [...]

  • 8
  •  read