• Computers and Concrete
• /
• v.1 no.3
• /
• pp.227-248
• /
• 2004

A combined experimental-computational study of a double edge-notched stone specimen subjected to tensile loading is presented. In the experimental part, the load-deformation response and the displacement field around the crack tip are recorded. An Electronic Speckle Pattern Interferometer (ESPI) is used to obtain the local displacement field. The experimental results are used to validate a numerical model for the description of fracture using finite elements. The numerical model uses displacement discontinuities to model cracks. At the discontinuity, a plasticity-based cohesive zone model is applied for monotonic loading and a combined damage-plasticity cohesive zone model is used for cyclic loading. Both local and global results from the numerical simulations are compared with experimental data. It is shown that local measurements add important information for the validation of the numerical model. Consequently, the numerical models are enhanced in order to correctly capture the experimentally observed behaviour.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.24 no.6
• /
• pp.777-784
• /
• 2000

A Navier-Stokes code based on an unstructured finite volume method is used to simulate the MIT flapping foil experiment. A low Reynolds number ${\kappa}-{\varepsilon}$ turbulence model is used to close the Reynolds averaged Navier-Stokes equations. Computations are carried out for the whole experimental domain involving two flapping foils and a downstream hydrofoil. The computational domain is meshed with unstructured quadrilateral elements, partly structured. Numerical solutions show good agreement with experiment. The first harmonics of the velocity in the boundary layer shows local peak value inside the boundary layer and also local minimum near the edge of boundary layer. It is intensified as it develops along the blade surface. This is shown to be caused as the unsteadiness inside the boundary layer is being convected at a speed less than the free stream value. It is also shown that there is negligible mixing of the unsteadiness between the boundary layer and the free stream.

• Structural Engineering and Mechanics
• /
• v.3 no.1
• /
• pp.11-23
• /
• 1995

The examined model consists of two substructures linked by a right angle rigid line joint. One element is a wall loaded externally along its upper edge by an uneven vertical load. The other element, defined as a plate, is not loaded. Stresses and displacements in the vicinity of the joint are analysed, considering the lateral distribution which leads to three-dimensional effects. The proposed solution combines classical approach with numerical means, using appropriate stress distribution polynomial functions along the joint. Space structure constructions supply cases of interest.

• Proceedings of the KOSOMBE Conference
• /
• v.1995 no.11
• /
• pp.34-38
• /
• 1995

Using a 3-D finite element method (FEM), the biomechanical characteristics of a threaded truncated acetabular component and a porous coated hemispherical acetabular component were studied. The Von-Mises stress/strain patterns in the acetabulum reconstructed with these two different types of cementless acetabular cups were investigated. The geometry and dimensions of human hemi-pelvis used in the present shape modeling for finite element analysis were scanned with a 3-D laser scanner(TDS-9000, Cyberware, USA). The scanned data was numerically handled with a shape modelling software 'Pro-Engineer'. Using 19836, 16853 tetrahedral elements, respectively, the stress and displacement field of the acetabulum reconstructed with the two different types of the acetabular components were computed. While the hemi-sphere component was found to show a relatively similar stress/strain patterns to those in the normal hip, the results with the threaded cup showed a considerably different patterns from those in the normal condition. Several regions in cancellous bone near the threads and the edge of the truncated cup was found to be overstressed, especially in the superior-lateral part of the acetabulum. It was postulated that the excessive reaming-out of subchondral bone layer when the truncated cup was used can cause the presence of these overstressed regions of cancellous bone. This theoretical prediction for the implanted acetabulum appeared to consistent with the pathological observation of proximal/medial migration of the threaded truncated acetabular prostheses in the previous publications.

• Structural Engineering and Mechanics
• /
• v.28 no.1
• /
• pp.107-127
• /
• 2008

Wind turbine blades are increasing in magnitude without a proportional increase of stiffness for which reason geometrical and inertial nonlinearities become increasingly important. Often these effects are analysed using a nonlinear truncated expansion in undamped fixed base mode shapes of a blade, modelling geometrical and inertial nonlinear couplings in the fundamental flap and edge direction. The purpose of this article is to examine the applicability of such a reduced-degree-of-freedom model in predicting the nonlinear response and stability of a blade by comparison to a full model based on a nonlinear co-rotating FE formulation. By use of the reduced-degree-of-freedom model it is shown that under strong resonance excitation of the fundamental flap or edge modes, significant energy is transferred to higher modes due to parametric or nonlinear coupling terms, which influence the response and stability conditions. It is demonstrated that the response predicted by such models in some cases becomes instable or chaotic. However, as a consequence of the energy flow the stability is increased and the tendency of chaotic vibrations is reduced as the number of modes are increased. The FE model representing the case of infinitely many included modes, is shown to predict stable and ordered response for all considered parameters. Further, the analysis shows that the reduced-degree-of-freedom model of relatively low order overestimates the response near resonance peaks, which is a consequence of the small number of included modes. The qualitative erratic response and stability prediction of the reduced order models take place at frequencies slightly above normal operation. However, for normal operation of the wind turbine without resonance excitation 4 modes in the reduced-degree-of-freedom model perform acceptable.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.19 no.3 s.73
• /
• pp.259-269
• /
• 2006

In this paper, the stress intensity factor, energy release rate and crack opening displacement are computed using the finite element method for axisymmetric viscoelastic cylinders with the penny-shaped and circumferential cracks. The triangular elements with quarter point nodes are used to describe the stress singularity around the crack edge. The analytical solutions are also derived by using the elastic-viscoelastic correspondence principle and compared with the numerical results to show the validity and accuracy of the presented method. Viscoelastic materials are assumed to behave elastically in dilatation and like a three-parameter standard linear solid.

• Computers and Concrete
• /
• v.21 no.3
• /
• pp.321-333
• /
• 2018

Flat-slabs, being a significant structural component, not only reduce the dead load of the structure but also reduce the amount of concrete required for construction. Moreover the use of recycled aggregates lowers the impact of large scale construction to nearby ecosystems. Recycled aggregate based concrete being a quasi-brittle material shows enormous cracking during failure. Crack growth in flat-slabs is mostly in sliding mode (Mode II). Therefore sufficient sections need to be provided for resistance against such failure modes. The main objective of the paper is to numerically determine the ultimate load carrying capacity of two self-similar flat-slab specimens and validate the results experimentally for the natural aggregate as well as recycled aggregate based concrete. Punching shear experiments are carried out on circular flat-slab specimen on a rigid circular knife-edge support built out of both normal (NAC) and recycled aggregate concrete (RAC, with full replacement). Uniaxial compression and bending tests have been conducted on cubes, cylinders and prisms using both types of concrete (NAC and RAC) for its material characterization and use in the numerical scheme. The numerical simulations have been conducted in ABAQUS (a known finite element software package). Eight noded solid elements have been used to model the flat slab and material properties have been considered from experimental tests. The inbuilt Concrete Damaged Plasticity model of ABAQUS has been used to monitor crack propagation in the specimen during numerical simulations.

• Coupled systems mechanics
• /
• v.12 no.2
• /
• pp.183-197
• /
• 2023

Scarce research has been published on crack propagation fracture of flywheels manufactured with carbon fiber-reinforced polymers. The present work deals with a calculation method to determine the conditions for which a crack propagates in the axial direction of the flywheel. The assumptions are: flywheels made with just a single thick ply or ply clustering laminates, oriented following the hoop direction; a single crack is analyzed in the plane defined by the hoop and axial directions; the crack starts close to one of the free edges; its axial length is initially large enough so that its tip is far away from that free edge, and the crack expands the entire circumferential perimeter and keeps its concentric position. The developed method provides information for a good design of flywheels. It is concluded that a fracture-based crack propagation criterion generally occurs at a lower speed than a stress-based criterion. Also, that the evolution of failure with thickness using the fracture criterion is exponential, demonstrating that thin flywheels are relatively not sensitive to crack propagation, whereas thick ones are very prone.

• Proceedings of the Korean Society for Noise and Vibration Engineering Conference
• /
• 2005.11a
• /
• pp.788-791
• /
• 2005

Nodal displacements are referred to the Initial configuration in the total Lagrangian formulation and to the last converged configuration in the updated Lagrangian formulation. This research proposes a relative nodal displacement method to represent the position and orientation for a node in truss structures. Since the proposed method measures the relative nodal displacements relative to its adjacent nodal reference frame, they are still small for a truss structure undergoing large deformations for the small size elements. As a consequence, element formulations developed under the small deformation assumption are still valid fer structures undergoing large deformations, which significantly simplifies the equations of equilibrium. A structural system is represented by a graph to systematically develop the governing equations of equilibrium for general systems. A node and an element are represented by a node and an edge in graph representation, respectively. Closed loops are opened to form a spanning tree by cutting edges. Two computational sequences are defined in the graph representation. One is the forward path sequence that is used to recover the Cartesian nodal displacements from relative nodal displacements and traverses a graph from the base node towards the terminal nodes. The other is the backward path sequence that is used to recover the nodal forces in the relative coordinate system from the known nodal forces in the absolute coordinate system and traverses from the terminal nodes towards the base node. One closed loop structure undergoing large deformations is analyzed to demonstrate the efficiency and validity of the proposed method.

• Geophysics and Geophysical Exploration
• /
• v.7 no.2
• /
• pp.148-154
• /
• 2004

This article reviews the development of three-dimensional (3-D) magnetotelluric (MT) modeling. The 3-D modeling of electromagnetic fields is essential in understanding the physics of MT soundings, and in implementing an inversion method to reconstruct a 3-D resistivity image. Although various numerical schemes have been developed over the last two decades, practical methods have been quite limited. However, the recent rapid improvement in computer speed and memory, as well as the advance in iterative solution algorithms for a large system of equations, makes it possible to model the MT responses of complex 3-D structures, which have been very difficult to simulate before. The use of staggered grids in finite difference method has become popular, conserving a magnetic flux and an electric current and allowing for realistic discontinuous fields. The convergence of numerical solutions has been greatly accelerated by adopting Krylov subspace methods, proper preconditioning techniques, and static divergence corrections. The vector finite-element method using edge elements is also free from the discontinuity problem, and seems a natural choice for modeling complex structures including irregular topography because its flexibility allows one to capture full geometric complexity.