1. Optimization of elastic systems using absolute nodal coordinate finite element formulationBojan Vohar, Marko Kegl, Zoran Ren, 2006, short scientific article Abstract: An approach to a shape optimization of elastic dynamic multibody systems is presented. The proposed method combines an appropriate shape parameterization concept and recently introduced finite element type using absolute nodal coordinate formulation (ANCF). In ANCF, slopes and displacements are used as the nodal coordinates instead of infinitesimal or finite rotations. This way one avoids interpolation of rotational coordinates and problems with finite rotations. ANCF elements are able to describe nonlinear deformation accurately; therefore, this method is very useful for simulations of lightweight multibody structures, where large deformations have to be taken into account. The optimization problem is formulated as a nonlinear programming problem and a gradient-based optimization procedure is implemented. The introduced optimization design variables are related to the cross-sectional parameters of the element and to the shape of the whole structure. The shape parameterization is based on the design element techniqueand a rational B ezier body is used as a design element. A body-like design element makes possible to unify the shape optimization of both simple beams and beam-like (skeletal) structures. Keywords: mechanics, dynamics of material systems, multibody systems, elastic mechanical systems, manipulators, dynamically loaded beams, optimum shape design, absolute nodal coordinate formulation, design element technique, finite element method Published: 31.05.2012; Views: 1128; Downloads: 67 Link to full text |
2. Velocity vorticity-based large eddy simulation with the bounadr element methodJure Ravnik, Leopold Škerget, Matjaž Hriberšek, 2006, published scientific conference contribution (invited lecture) Abstract: A large eddy simulation using the velocity-vorticity formulation of the incompressible Navier-Stokes equations in combination with the turbulent heat transfer equation is proposed for the solution of the turbulent natural convection drive flow in a 1:4 enclosure. The system of equations is closed by an enthropy-based subgrid scale model.The Prandtl turbulent number is used to estimate turbulent diffusion in the heat transfer equation. The boundary element method is used to solve the kinematics equation and estimate the boundary vorticity values. The vorticity transport equation is solved by FEM. The numerical example studied in this paper is the onset of a turbulent flow regime occuring at high Rayleigh number values ▫$(Ra=10^7-10^10)$▫. The formation of vortices in the boundary layer is observed, along with buoyancy driven diffusive convective transport. Quantitative comparison with the laminar flow model and the worh of other authors is also presented in terms of Nusselt number value oscillations. Keywords: fluid mechanics, incompressible viscous fluid, turbulent flow, velocity vorticity formulation, finite element method, large eddy simulation Published: 31.05.2012; Views: 1218; Downloads: 17 Link to full text |
3. On critical buckling load estimation for slender transversely cracked beam-columns by the application of a simple computational modelMatjaž Skrinar, 2008, original scientific article Abstract: This paper brings new insights into the implementation of a simplified computational model in the prediction of buckling load Pcr for slender beam-type structures with a transverse crack. From among several approaches discussed, two of them produced applicable results exhibiting considerably good agreement with those values from more precise and complex computational models. In the first approach, the critical load value is obtained from numerical solutions of analytically expressed characteristic equations (obtained from governing differential equations). Although producing excellent results, this approach limits the application since an analytical solution of the governing differential equation can only be obtained for moderate structures. The second approach implements a new cracked beam-columnfinite element, derived at on the basis of a fairly accurate approximation of the governing differential equation's solution. It allows for flexible utilization and also yields the smallest compact computational model, thus exhibiting itself as very suitable for inverse identification problems. Numerical examples covering several structures with different boundary conditions are briefly presented in order to support the discussed approaches. The results obtained using the presented approaches are further compared with those values from either references or more complex models, thus clearly proving the quality of the presented compact FE model. Keywords: beam columns, transverse cracks, stability problems, buckling load, computational model, finite element method Published: 31.05.2012; Views: 1601; Downloads: 60 Link to full text |
4. Crack tip shielding or anti-shielding due to smooth and discontinuous material inhomogeneitiesN.K. Simha, Franz Dieter Fischer, Otmar Kolednik, Jožef Predan, G.X. Shan, 2005, original scientific article Abstract: This paper describes a theoretical model and related computational methods forexamining the influence of inhomogeneous material properties on the crack driving force in elastic and elastic-plastic materials. Following the configurational forces approach, the crack tip shielding or anti-shielding dueto smooth (e.g. graded layer) and discontinuous (e.g. bimaterial interface)distributions in material properties are derived. Computational post-processing methods are described to evaluate these inhomogeneity effects.The utility of the theoretical model and computational methods is demonstrated by examining a bimaterial interface perpendicular to a crack in elastic and elastic-plastic compact tension specimens. Keywords: fracture mechanics, fracture toughness, composite materials, layered material, inhomogeneity, cracks, finite element method, elastic bimaterials Published: 01.06.2012; Views: 1136; Downloads: 50 Link to full text |
5. A FEM method for eddy current analysis in laminated mediaOszkar Biró, Kurt Preis, Igor Tičar, 2005, original scientific article Abstract: A method is presented to compute three-dimensional eddy current distributions in laminated media by means of the finite element method. In a first step, the laminated medium is assumed to have an anisotropic conductivity with zero (or very low) value in the direction normal to the laminations. In a second step, the eddy currents within the laminates are computed by solving the quasistatic electromagnetic field individually in each sheet. In these essentially two-dimensional analyses, the boundary conditions are taken from the three dimensional field distribution determined in the first step. Comparisons with results obtained from a finite element model taking account of each laminate prove the validity of the method. Keywords: electromagnetic field, eddy currents, finite element method analysis, laminates Published: 01.06.2012; Views: 1153; Downloads: 37 Link to full text |
6. On the application of a simple computational model for slender transversely cracked beams in buckling problemsMatjaž Skrinar, 2007, original scientific article Abstract: This paper discusses the implementation of a simplified computational model that is widely used for the computation of transverse displacements in transversely cracked slender beams into the Euler's elastic flexural buckling theory. Two alternatives are studied instead of solving the corresponding differential equations to obtain exact analytical expressions for the buckling load ▫$P_{cr}$▫ due to the complexity of this approach. The first approach implements wisely selected polynomials to describe the behavior of the structure, which allows the derivation of approximate expressions for the critical buckling load. Although the relevance of the results strongly depends on the proper prime selection of the polynomial, it is shown that the later upgrading of the polynomials can lead to even more reliable results. The second approach is a purely numerical approach and presents the geometrical stiffness matrix for a beam finite element with a transverse crack. To support the discussed approaches, numerical examples covering several structures with different boundary conditions are briefly presented. The results obtained with the presented approaches are further compared with the values from enormous 2D finite elements models, where a detailed description of the crack was achieved with the discrete approach. It is evident that the drastic difference in the computational effort is not reflected in the significant differences in the results between the models. Keywords: columns, transverse cracks, stability problems, buckling load, computational model, polynomial solutions, finite element method, geometrical stiffness matrix Published: 01.06.2012; Views: 1315; Downloads: 56 Link to full text |
7. New finite element for transversely cracked slender beams subjected to transverse loadsMatjaž Skrinar, Tomaž Pliberšek, 2007, original scientific article Abstract: The paper covers the derivation of a new finite element for beams with transverse cracks. The derivation is based on a simplified computational model that has already proved itself suitable for the inverse identification of cracks. The model embodies all the necessary major information about the structure's response from the inverse identification point of view, where the presence and location and, if possible, the depth of the crack should be detected from the measured response, usually dynamic. In such situations the stress distributions obtained from 2D finite elements analysis are not as important as the computational model being capable of reliably describing the displacement of the structure. However, from numerical studies it also became evident that the relevance of the model decreases with element thickness. This indicated that shear forces should be included in the analysis process. Therefore, derivation of a new finite element with the inclusion of shear forces effect has been executed. The stiffness matrix for transversely cracked slender beams, as well as the derivation of interpolation functions is presented and all expressions are given in symbolic forms. The example shows that, although with significantly less computational effort than with 2D FE meshes, significant improvement in transverse displacements can be obtained with the presented beam finite element. Keywords: beams, transverse cracks, computational model, finite element method, transverse displacements Published: 01.06.2012; Views: 1188; Downloads: 63 Link to full text |
8. Solving exterior problems of wave propagation based on an iterative variation of local DtN operatorsMiroslav Premrov, Igor Špacapan, 2004, original scientific article Abstract: This paper discusses the problem of wave reflection from the fictitious boundary, with particular regard to the higher harmonic modes. This problem occurs when solving the wave equation in exterior domains using an asymptotic local low-order Dirichlet-to-Neumann (DtN) map for computational procedures applied to a finite domain. We demonstrate that the amplitudes of the reflected fictitious harmonics depend on the wave number, the location of the fictitious boundary, as well as on the local DtN operator used in the computations. Moreover, we show that a constant value of the asymptotic local low-order operator cannot sufficiently eliminate the amplitudes of all reflected waves, and that the results are poor especially for higher harmonics. We propose therefore an iterative method, which varies the tangential dependence of the local operator in each computational step. We only discuss some logical and interesting choices for the operators although this method permits several possibilities on how to vary the operator. The method is simple to apply and the presented examples demonstrate that the accuracy is considerably improved by iterations. Keywords: wave motion, infinite domains, fictitious boundary, radiation condition, finite element method, Dirichlet-to-Neumann map Published: 01.06.2012; Views: 960; Downloads: 46 Link to full text |
9. Symmetry preserving algorithm for large displacement frictionless contact by the pre-discretization penalty methodDušan Gabriel, Jiři Plešek, Miran Ulbin, 2004, original scientific article Abstract: A three-dimensional contact algorithm based on the pre-discretization penalty method is presented. Using the pre-discretization formulation gives rise to contact searching performed at the surface Gaussian integration points. It is shown that the proposed method is consistent with the continuum formulation ofthe problem and allows an easy incorporation of higher-order elements with midside nodes to the analysis. Moreover, a symmetric treatment of mutually contacting surfaces is preserved even under large displacement increments. Theproposed algorithm utilizes the BFGS method modified for constrained non-linear systems. The effectiveness of quadratic isoparametric elements in contact analysis is tested in terms of numerical examples verified by analytical solutions and experimental measurements. The symmetry of the algorithm is clearly manifested in the problem of impact of two elastic cylinders. Keywords: mechanics, numerical methods, contacting surfaces, contact problems, 3D contact algorithm, discretization, higher order elements, finite element method, Gauss point search, pre-discretization penalty method Published: 01.06.2012; Views: 1488; Downloads: 42 Link to full text |
10. Weakly coupled analysis of a blade in multiphase mixing vesselMatej Vesenjak, Zoran Ren, Matjaž Hriberšek, 2004, short scientific article Abstract: Two or more physical systems frequently interact with each other, where the independent solution of one system is impossible without a simultaneous solution of the others. An obvious coupled system is that of a dynamic fluid-structure interaction. š8đ In this paper a computational analysis of thefluid-structure interaction in a mixing vessel is presented. In mixing vessels the fluid can have a significant influence on the deformation of blades during mixing, depending on speed of mixing blades and fluid viscosity.For this purpose a computational weakly coupled analysis has been performed to determine the multiphase fluid influences on the mixing vessel structure. The multiphase fluid field in the mixing vessel was first analyzed with the computational fluid dynamics (CFD) code CFX. The results in the form of pressure were then applied to the blade model, which was the analysed with the structural code MSC.visualNastran forWindows, which is based on the finiteelement method (FEM). Keywords: fluid mechanics, solid mechanics, coupled problems, computational analysis, two-phase fluid, mixing blade, pressure distribution, finite volume method, finite element method Published: 01.06.2012; Views: 875; Downloads: 51 Link to full text |