ADINA Publications

Page 30

The Theory used in ADINA is richly documented in the following books by K.J. Bathe and co-authors

Finite Element Procedures

Finite Element Procedures in Engineering Analysis

Numerical Methods in Finite Element Analysis

The Mechanics of Solids and Structures — Hierarchical ...

The Finite Element Analysis of Shells — Fundamentals

Inelastic Analysis of Solids and Structures

To Enrich Life
(Sample pages here)


Following are more than 700 publications — that we know of — with reference to the use of ADINA. The pages give the Abstracts of some papers published since 1986 referring to ADINA. The most recent papers are listed first. All these papers may be searched using the box:

(latest on Page 43)
Page 43 Page 42 Page 41 Page 40 Page 39 Page 38 Page 37 Page 36 Page 35 Page 34
Page 33 Page 32 Page 31 Page 30 Page 29 Page 28 Page 27 Page 26 Page 25 Page 24
Page 23 Page 22 Page 21 Page 20 Page 19 Page 18 Page 17 Page 16 Page 15 Page 14
Page 13 Page 12 Page 11 Page 10 Page 9  Page 8  Page 7  Page 6  Page 5  Page 4 
Page 3  Page 2  Page 1 

Resonant pumping in a multilayer impedance pump

L. Loumes,1 I. Avrahami,2 and M. Gharib1

1Division of Engineering and Applied Sciences, California Institute of Technology, 1200 E, California Blvd., MC 205-45, Pasadena, California 91125, USA
2Medical Engineering, AFEKA ý Tel Aviv Academic College of Engineering, Israel

Physics Of Fluids 20, 023103, 2008

Abstract: This paper introduces the concept of multilayer impedance pump, a novel pumping mechanism inspired by the embryonic heart structure. The pump is a composite two-layer fluid-filled elastic tube featuring a thick gelatinous internal. Pumping is based on the impedance pumping mechanism. In an impedance pump, elastic waves are generated upon external periodic compressions of the elastic tube. These waves propagate along the tube’s walls, reflect at the tube’s extremities, and drive the flow in a preferential direction. The originality in the multilayer impedance pump design relies on the use of the thick internal gelatinous layer to amplify the elastic waves responsible for the pumping. As a consequence, only small excitations are needed to produce significant flow. This fully coupled fluid-structure interaction problem is solved for the flow and the structure using the finite element method over a relevant range of frequencies of excitation. Results show that the multilayer impedance pump is a complex system that exhibits a resonant response. Flow output and inner wall motion are maximal when the pump is actuated at the resonant frequency. The wave interaction mechanism present in an impedance pump is described here in details for the case of a multilayer impedance pump. Using energy balance for the passive portion of the elastic tube, we show that the elastic tube itself works as a pump and that at resonance maximum energy transmission between the elastic tube and the fluid occurs. Finally, the pump is especially suitable for many biomedical applications.


Numerical Simulation of Hydraulic Turbine Based on Fluid-Structure Coupling

Demin Liu1, Shuhong Liu1, Yulin Wu1 and Xiao-bing Liu2

1State Key Laboratory of Hydroscience & Engineering, Dept. of Thermal Engineering, Tsinghua University, China
2Schoolof Energy and Environment, Xihua University, Chengdu, China

Proc. The 4th International Symposium on Fluid Machinery and Fluid Engineering, 345-351, 2008

Abstract: By adopting the arbitrary Lagrange-Euler (ALE) method of software ADINA, fluid-structure coupling (FSC) calculation of a Francis turbine is conducted. The vibration frequency and mode-of the runner in the air and water are obtained. The calculation results show that runner frequency of the turbine is reduced in certain degree under the effect of water pressure and viscous force, and the mode is changed, too. By using dynamic fracture mechanics, the possible cracking damage of the runner is predicated.

Keywords: fluid-structure coupling - hydraulic turbine – vibration - numerical simulation


Three-Dimensional Finite Element Simulations of Compression Tests on Bimrock

M. Barbero, M. Bonini, M. Borri-Brunetto

Department of Structural and Geotechnical Engineering, Politecnico di Torino, Italy

Proc. 12th Int. Conf. IACMAG, 2008

Abstract: The paper deals with the mechanical behaviour of “bimrock” (block in matrix rock), characterized by a heterogeneous structure, constituted by rock blocks included in a small-grained, well-cemented material (matrix) with lower mechanical characteristics. The bimrock is a very peculiar material that may give rise to serious problems in engineering applications. To properly model this kind of problems it is necessary to suitably characterize the bimrock. In order to contribute to the understanding of the bimrock mechanical behaviour and to individuate strength and deformability laws suitable for the bimrock as an equivalent continuum, as the problem is clearly three-dimensional, 3D numerical analyses have been carried out by using the finite element method implemented in the Adina code. Uniaxial and triaxial compression tests on bimrock laboratory samples are simulated. Different volumetric proportions of the blocks are considered. The numerical results show the influence of the volume proportion of the blocks on the strength and deformability of bimrock, confirming the results obtained with previous 2D analyses. Furthermore, the results suggest that the interfaces between the blocks and the matrix play an important role, and may have an influence on the mechanical behaviour of bimrock.

Keywords: Bimrock — melange — heterogeneous materials — numerical modelling


Seismic Design Of Two New York Bridges

B. Yin

Hardesty & Hanover, LLP, New York, USA

Proc. The 14th World Conference on Earthquake Engineering, 2008

Abstract: The Third Avenue Bridge and Willis Avenue Bridge are two adjacent bridges each
carrying one directional traffic across Harlem River in New York City. The replacement of 100-year old Third Avenue Bridge is complete. The new 5-lane crossing consists of 17 approach spans and a 107-m long swing span for an overall bridge length of 457m. Including on-grade approaches, the overall project length is 1067m. Construction of the $118.8 million bridge replacement project began in 2001 and completed in 2005. The Willis Avenue Swing Bridge over the Harlem River is being replaced under a $612 million project which is massive in scope as it extends over a mile in length, passes over the Harlem River and an adjoining railyard and provides connections between two major highways as well as three major arterial streets. The new alignment not only dramatically improves the alignment from that of the 100 year old existing bridge but also facilitates maintaining 70,000 vehicles per day of roadway traffic as well as maintaining navigation on the river. The project centerpiece is a new four lane, 106 meter long swing span. This paper focuses on the seismic analysis and seismic design of the two bridges. Challenges faced by the designers of the new bridge, and solutions developed to meet these challenges will be introduced. Since similarities exist between these two bridges, more detailed description will be on one of them, i.e. Third Avenue Bridge, and the other one is only in brief.

Keywords: Seismic — Design — Swing Bridge — Harlem River — Pivot Pier — Pivot Bearing — Truss — New York — Non-linear


Modelling of Tracer Tests in a Shear Zone at the Grimsel Test Site

A. Pudewills

Institut für Nukleare Entsorgung

Forschungszentrum Karlsruhe in der Helmholtz-Gemeinschaft — Wissenschaftliche Berichte FZKA 7417, 2008

Abstract: In the Colloid Formation and Migration (CFM) Project at the Grimsel Test Site an experiment has been planned in order to study the generation of bentonite colloid and colloid facilitated radionuclides transport in a low gradient groundwater flow. During the first phase of this project a series of tracer tests have been performed to evaluate suitable flow fields for long-term colloid migration experiments within the fractured shear zone. The tracer tests with uranine have been accompanied by numerical modelling to support the in situ measurements. The overall objective of this modelling work is to understand the flow path and the relevant transport processes within the water-conducting shear zone, considering the low hydraulic gradients. The numerical analysis was performed with the ADINA-F finite element code. The mathematical model is based on the Darcy’s law for groundwater flow and the advection-dispersion equation for the uranine transport through a porous shear zone. For the numerical simulation a two dimensional model in the plane of the shear zone has been used. The computational domain considers the circular tunnel sealed with a surface packer, both injection wells and the extraction hole. Both field data and modelling results provide a fairly consistent picture of the flow and transport characteristics in the shear zone at the planned CFM location. Furthermore, the simulation results match the experimental breakthrough curves fairly well and give confidence in the numerical model used. The estimated hydraulic and transport parameters from calibration work for two dipole fields show somewhat differences caused by natural heterogeneities in the hydraulic conductivity of the shear zone.


Seismic performance of a long span arch bridge taking account of fluctuation of axial force

H. Chavez and J.J. Alvarez

Universidad Michoacana de San Nicolas de Hidalgo, Mexico

Proc. The 14th World Conference on Earthquake Engineering, 2008

Abstract: The seismic behavior of a long-span arch bridge is investigated based on a series of nonlinear dynamic analyses. The bridge was modeled with ADINA by means of elements with nonlinear behavior Since the axial force and bending moment interaction is significant in the arch rib under strong ground motions, this effect is included in the analyses by frame elements varying their moment-curvature relations. It was found from the analyses that large axial forces, as much as 1.5 times the dead axial force are developed in the arch rib during excitations, so moment-curvature relations change in the arch rib and hence should be considered in design.

Keywords: Arch bridges — nonlinear dynamic analysis — seismic performance


Experimental And Numerical Study On Dynamic Properties Of Friction-Resistant Dry-Masonry Structures

S. Dorjpalam1, H. Kawase2, K. Yamaguchi3, and S.O.Citak4

1 Structural Engineer
2 Disaster Prevention Research Institute, Kyoto University, Kyoto, Japan
3 Grad. School of Human-Environment Studies, Kyushu University, Fukuoka, Japan
4 Ohsaki Research Institute, Tokyo, Japan

Proc. The 14th World Conference on Earthquake Engineering, 2008

Abstract: Dynamic properties of friction-resistant dry-masonry SRB-DUP (Steel Reinforced Brick based on Distributed Unbond Prestress) structures and structural elements are studied through experimental tests and numerical simulation analyses. First, dynamic tests are conducted on masonry wall elements, and basic dynamic properties are found. Next, a numerical model is formulated using micro-modeling technique, where constituent materials (bricks and reinforcing steel elements) are modeled as separate entities. The model is verified by analyzing the experimental wall models. Then simulation models of various wall elements are built and analyzed. Based on the analyses results a macro-model of a real size structure is built and analyzed for strong ground motion records. We found that the SRB-DUP structure shows quite strong seismic resistance.

Keywords: Friction-resistant — Reinforced Masonry — umerical Modeling — Simulation — FEM Analysis — Environmentally Friendly Structural System


Old and new solutions for contact splices in columns

J. Lindner

TU Berlin, Department of Steel Structures, c/o Consulting office, Bismarckallee 4, 14193 Berlin, Germany

Journal of Constructional Steel Research, 64:833–844, 2008

Abstract: Compression forces in a splice can be transferred by contact. Accordingly additional imperfections can be noticed and have to be accounted for. If a contact splice is placed inside the length of a column these initial additional imperfections will have an influence on the buckling behaviour and thus the load carrying capacity of the column. Report on theoretical investigations on the member behaviour and on the design of a splice and its fastenings is given. Also full scale tests on different types of splices are referred to. The test results are compared to ultimate load calculations and respective design methods are proposed. Especially for contact splices in columns arranged directly one on top of the other the stability behaviour for a splice at midspan is examined and leads to a special buckling curve.

Keywords: Columns — Splices — Contact — Stability — Initial imperfections — Full scale tests — Ultimate load calculations — Steel structures


Stability analysis of gravity dams on sloping layered rock foundation against deep slide

L. Zhang, D. Wang, H. Zhang, and W. Wang

Hohai University, Nanjing, China, 210098

Proc. 11th ASCE Aerospace Division International Conference (Earth and Space 2008),

Abstract: In this study, the stability of gravity dam on a sloping layered rock foundation is analyzed, and the stability problem of gravity dam against deep slide is particularly focused. The commercial finite element software ADINA is used to develop the constitutive law of anisotropic layered materials, such as the sloping layered rock foundation. Based on a strength reduction method, the stability of gravity dams on sloping layered rock foundation is evaluated, and a criterion for evaluating the instability is proposed. As an application, the proposed methodology for stability analysis of gravity dams against deep slide is applied to a gravity dam in ArHai Hydraulic Power Station in Yunnan Province, China. The complexity of dam structures as well as the anisotropic nature of layered rock foundation is tackled with a finite element simulation, from which the deep sliding mechanism of gravity dam and its stability against slide are analyzed. By using the different vertical and horizontal parameters in the layered rock foundation model, the stability of gravity dam against deep slide is analyzed. From the thorough numerical parametric study, it indicates that the methodology is effective, feasible and generally applicable.The methodology and numerical analysis presented provide a better understanding of stability of gravity dam on layered rock foundation and can be used as a valuable reference for design analysis of gravity dam.


Modal Identification Study of Vincent Thomas Bridge Using Simulated Wind-Induced Ambient Vibration Data

X. He1, B. Moaveni1, J.P Conte2, A. Elgamal2, S.F. Masri2

1 Tufts University
2 University of California, San Diego

Postprints, UC San Diego, DOI: 10.1111/j.1467-8667.2008.00544.x, 2008

Abstract: In this paper, wind-induced vibration response of Vincent Thomas Bridge, a suspension bridge located in San Pedro near Los Angeles, California, is simulated using a detailed three-dimensional finite element model of the bridge and a state-of-the-art stochastic wind excitation model. Based on the simulated wind-induced vibration data, the modal parameters (natural frequencies,damping ratios, and mode shapes) of the bridge are identified using the data-driven stochastic subspace identification method. The identified modal parameters are verified by the computed eigenproperties of the bridge model. Finally, effects of measurement noise on the system identification results are studied by adding zero-mean Gaussian white noise processes to the simulated response data. Statistical properties of the identified modal parameters are investigated under increasing level of measurement noise. The framework presented in this paper will allow to investigate the effects of various realistic damage scenarios in long-span cable-supported (suspension and cable-stayed) bridges on changes in modal identification results. Such studies are required in order to develop robust and reliable vibration-based structural health monitoring methods for this type of bridges, which is a long-term research objective of the authors.

Keywords: Vincent Thomas Bridge — experimental modal analysis — ambient vibration —  system identification


Fluid-structure interaction for a jacket model structure with imposed displacement

A.C. Mendes, H.J.D. Correia

Laboratory of Fluid Mechanics, Universidade da Beira Interior, Covilhã, Portugal.

Proc. Eighteenth Int. Offshore and Polar Eng. Conf., 2008

Abstract: In this paper we describe a series of experimental tests performed with
a 1:45  scale jacket structure in a water basin. The model structure is assembled on a seismic table that performs a controlled motion, with amplitudes and frequencies that could correspond to actual earthquake conditions translated to the scale of the model. The structural response to the imposed displacement is analyzed in terms of the accelerations and reaction forces at the base of the structure, first in air and afterwards in water. The experimental measurements are equally correlated with the computer predictions of dynamic loading that have been obtained with a FEM commercial software - the ADINA.

Keywords: Offshore jacket structure — earthquakes — fluid-structure interaction —  displacement based design


Flexural Strengthening of RC Beams with Externally Bonded CFRP Systems: Test Results and 3D Nonlinear FE Analysis

Renata Kotynia1, Hussien Abdel Baky2, Kenneth W. Neale3, and Usama A. Ebead4

1 Assistant Professor, Dept. of Concrete Structures, Technical Univ. of Lodz, Al. Politechnik 6, 90-924 Lodz, Poland.
2 Ph.D. Candidate, Dept. of Civil Engineering, Univ. of Sherbrooke, Sherbrooke, Quebec, Canada J1K 2R1.
3 Canada Research Chair in Advanced Engineered Material Systems, Dept. of Civil Engineering, Univ. of Sherbrooke, Sherbrooke, Quebec, Canada J1K 2R1
4 Assistant Professor, Dept. of Civil and Environmental Engineering, College of Engineering, United Arab Emirates Univ., P.O. Box 17555, Al Ain, United Arab Emirates

Journal Of Composites For Construction, DOI 10.1061/(ASCE)1090-0268(2008)12:2(190),  2008

Abstract: This paper presents experimental results and a numerical analysis of the reinforced concrete (RC) beams strengthened in flexure with various externally bonded carbon fiber-reinforced polymer (CFRP) configurations. The aim of the experimental work was to investigate the parameters that may delay the intermediate crack debonding of the bottom CFRP laminate, and increase the load carrying capacity and CFRP strength utilization ratio. Ten rectangular RC specimens with a clear span of 4.2 m, categorized in two series, were tested to evaluate the effect of using the additional U-shaped CFRP systems on the intermediate crack debonding of the bottom laminate. Two different configurations of the additional systems were proposed, namely, continuous U-shaped wet layup sheets and spaced side-bonded CFRP L-shaped laminates. The fiber orientation effect of the side-bonded sheets was also investigated. A numerical analysis using an incremental nonlinear displacement-controlled 3D finite-element (FE) model was developed to investigate the flexural and CFRP/concrete interfacial responses of the tested beams. The finite-element model accounts for the orthotropic behavior of the CFRP laminates. An appropriate bond-slip model was adopted to characterize the behavior of the CFRP/concrete interface. Comparisons between the FE predictions and experimental results show very good agreement in terms of the load-deflection and load-strain relationships, ultimate capacities, and failure modes of the beams.

Keywords: Fiber reinforced polymers - Concrete beams – Concrete – reinforced - Flexural strength - Finite element method



Tensile radial stress in the spinal cord related to arachnoiditis or tethering: a numerical model

C.D. Bertram1, L.E. Bilston1 and M.A. Stoodley2

1 Biofluid Mechanics Laboratory, Faculty of Engineering, University of New South Wales, Sydney 2052, Australia
2 Prince of Wales Medical Research Institute, UNSW, Barker Street, Randwick 2031, New South Wales, Australia

Med Biol Eng Compu, 46:701–707, 2008

Abstract: Abstract Spinal arachnoiditis comprises fibrous scarring of the subarachnoid space, following spinal trauma or inflammation, and is often associated with syringomyelia. We hypothesised that cord-to-dura attachments could cause transient tensile cord radial stress, as pressure waves propagate. This was tested in a fluid-structure interaction model, simulating three types of cord tethering, with ‘arachnoiditis’ confined to a short mid-section of the cord. The annular system was excited abdominally with a short transient, and the resulting Young and Lamb waves and reflections were analysed. Radial mid-section tethering was less significant than axial tethering, which gave rise to tensile radial stress locally when the cord was not fixed cranially. Simulated as inextensible string connections to the dura, arachnoiditis caused both localised tensile radial stress and localised low pressure in the cord as the transient passed. The extent of these effects was sensitive to the relative stiffness of the dura and cord. Tensile radial stress may create a syrinx in previously normal cord tissue, and transiently lowered pressure may draw in interstitial fluid, causing the syrinx to enlarge if fluid exit is inhibited. The suggested mechanism could also explain the juxtaposition of syrinxes to regions of arachnoiditis.

Keywords: Fluid-structure interaction — Finite-element model — Cerebrospinal fluid —  Wave propagation — Syringomyelia


Modeling Of Electro Resistance Welding Tube Forming Process

Sebastián D’hers1, Pablo A. Bucello1, Martín N. Arregui1, Alberto P. Sibileau1, Miguel A. Cavaliere2 and Esteban Baralla3

1 Instituto Tecnológico de Buenos Aires. Av. E. Madero 399 C1106ACD, Buenos Aires, Argentina.
2 Tenaris Siderca, Dr. Simini 250, B2804MHA Campana, Buenos. Aires, Argentina.
3 Tenaris SIAT, Guatemala 3400, B1822AXZ Valentín Alsina, Buenos Aires, Argentina.

Mecánica Computacional, 27:627-638, 2008

Abstract: During the forming process in an electro resistance welding (ERW) forming line, a flat skelp evolves from a flat surface to a round shape, enforced by a series of roll stages. The formed round shape determines, together with welding parameters, the obtained weld and tube quality. In order to improve the process and achieve a better understanding of the skelp evolution towards the welding station, a finite element model has been developed. The process has been usually studied with analytical techniques and “trial and error” plant testing, while the use of finite element in modeling this type of process is quite recent. Process complexity imposes great computing efforts and heavy pre and post processing work. A finite element model has been proposed to predict the skelp behavior along the critical stages of the forming line. Large displacement, contact algorithms and J2 material models have been used in the problem description. Agreement with plant measurements, computational speed and robustness are the
main achieved goals.

Key Words: Metal forming - finite element method - roll forming - edge shape - electric resistance welding


Numerical Analysis Of Soil Response To Ice Scouring

E. Evgin and Z. Fu

Department of Civil Engineering, University of Ottawa, Ottawa, Ontario, Canada

Proceedings of the ASME 27th International Conference on Offshore Mechanics and Arctic Engineering, OMAE2008-57293, 2008

Abstract: Icebergs and ice ridges frequently scour the surface of seabed deposits. Ice scouring can be problematic for offshore pipelines and other seabed installations. In order to reduce the risk of failure, pipelines are buried in the seabed. However, a stationary or moving ice feature could cause a significant increase in stresses and deformation in the seabed soil deposits below the contact surface between the soil and the ice, and consequently, might result in structural failure of buried pipeline. Safe burial depth for pipelines has been the subject of both experimental and numerical studies. In this paper, two and three dimensional analyses are conducted using PLAXIS and ADINA. In these analyses, geometric and material nonlinearities are considered. In order to establish the validity of the finite element calculations, the experimental results reported in the literature and the numerical results obtained in the present study are compared. The emphasis is placed on determining the importance of (1) using interface elements between different materials such as soil and ice, soil and pipelines; (2) using the soil model correctly, and (3) using a three dimensional analysis rather than a two dimensional analysis. The changes taking placed in the deformation pattern and the stress states in the seabed soils resulting from ice scouring are determined. The effects of pipeline burial depth, the shape of the ice feature, and the characteristics of seabed soils on the stresses acting on the pipeline are evaluated.


Influence of different load models on gear crack path shapes and fatigue lives

S. Podrug1, D. Jelaska1 and S. Glodez2

1 University of Split, Faculty of Electrical Engineering, Mechanical Engineering and Naval Architecture, R. Boskovica bb, 21000 Split, Croatia
2 University of Maribor, Faculty of Mechanical Engineering, Smetanova 17, 2000 Maribor, Slovenia

Fatigue Fract Engng Mater Struct, 31:327–339, 2008

Abstract: A computational model for determination of the service life of gears with regard to bending fatigue at gear tooth root is presented. In conventional fatigue models of the gear tooth root, it is usual to approximate actual gear load with a pulsating force acting at the highest point of the single tooth contact. However, in actual gear operation, the magnitude as well as the position of the force changes as the gear rotates. A study to determine the effect of moving gear tooth load on the gear service life is performed. The fatigue process leading to tooth breakage is divided into crack-initiation and crack-propagation period. The critical plane damage model has been used to determine the number of stress cycles required for the fatigue crack initiation. The finite-element method and linear elastic fracture mechanics theories are then used for the further simulation of the fatigue crack growth.

Keywords: crack closure - crack initiation - crack propagation - gear tooth root - moving load


Numerical simulation of fluid-structure interaction for wind induced dynamic response of the 3rd Jinan Yellow River cable stayed bridge

Qi-Lin Zhang, Zhen-Hua Liu, Ying Zhou

Department of Building Engineering, Tongji University, Shanghai 200092, P.R. China

Proc. 6th Int. Conf. on Comp. of Shell and Spatial Structures, IASS-IACM 2008: “Spanning Nano to Mega”, Cornell University, Ithaca, NY, USA, 2008

Abstract: The bridge aerodynamic design has considerably progressed to a more and more analytical methodology by which the mechanical aspects concerning bridge stability are deeply understood and characterized. In this paper, a response analysis was carried out using ADINA to evaluate wind-induced dynamic response of the 3rd Jinan Yellow River cable-stayed steel box girder bridge. The wind induced dynamic coefficients and the wind pressure distribution coefficient of the steel box girder can be calculated.


Simulation of blood flow in a small-diameter vascular graft model with a swirl (spiral) flow guider

ZhiGuo Zhang1,2, YuBo Fan1, XiaoYan Deng1, GuiXue Wang3, He Zhang3 & Robert Guidoin4

1 Department of Bioengineering, Beihang University, Beijing 100083, China
2 School of Aerospace, Tsinghua University, Beijing 100084, China
3 College of Bioengineering, Chongqing University, Chongqing 400044, China
4 Department of Surgery, Laval University and Quebec Biomaterials Institute CHUQ, Quebec G1K7P4, Canada

Sci China Ser C-Life Sci, 51(10):913-921, 2008

Abstract: Small-diameter vascular grafts are in large demand for coronary and peripheral bypass procedures, but present products still fail in long-term clinical application. In the present communication, a new type of small-diameter graft with a swirl flow guider was proposed to improve graft patency rate. Flow pattern in the graft was simulated  numerically and compared with that in a conventional graft. The numerical results revealed that the swirl flow guider could indeed make the blood flow rotate in the new graft. The swirling flow distal to the flow guider significantly altered the flow pattern in the new graft and the velocity profiles were re-distributed. Due to the swirling flow, the blood velocity near the vessel wall and wall shear rate were greatly enhanced. We believe that the increased blood velocity near the wall and the wall shear rate can impede the occurrence of acute thrombus formation and intimal hyperplasia, hence can improve the graft patency rate for long-term clinical use.

Keywords: arterial prosthesis - acute thrombus formation - flow simulation - swirling (spiral) flow


Computational investigation of the tensile behaviour of the hard coated Ti-6Al-4V alloy

W. Ziaja

Department of Materials Science, Rzeszow University of Technology, ul. W. Pola 2, 35-959 Rzeszow, Poland

Journal of Achievements in Materials and Manufacturing Engineering, 26(2):175-178, 2008

Purpose: Modification of the surface layer of the titanium alloys is frequently applied in order to improve their tribological properties. Various surface engineering techniques can be used to produce hard coatings, e.g. composed of metallic carbides, nitrides or more recently DLC. The coating and substrate materials possess significantly different stiffness and strength properties. This can lead to premature failure of the usually elastic coating in case of plastic deformation of the substrate when the high stresses are encountered. Cracking of the hard coating leads to stress concentration and localized plastic deformation of the substrate that can modify macroscopic deformation behaviour of the system. In the paper the influence of coating and substrate properties on local plastic deformation of substrate material was numerically investigated.
Design/methodology/approach: Two dimensional finite element analysis of the process of tensile deformation of titanium alloy with hard elastic coating was carried out. Two cases were analyzed, i.e. with and without diffusion strengthened layer underlying the coating.
Findings: The influence of the difference in Young’s modulus between coating and substrate material, yield strength of substrate material, coating thickness and depth of the crack in the coating on local plastic deformation of substrate material was determined.
Research limitations/implications: Some extension of the numerical model should be pursued in order to take into account initiation of microcracks in surface layer of the coated material and process of coating
Practical implications: The results could be used in the element design process for selection of parameters of surface layer with complex structure for load bearing applications.
Originality/value: The mechanical behaviour of hard coated material was most frequently studied for indentation and friction conditions and much less investigations were carried out for coated systems under tension or compression.

Keywords: Computational material science - Titanium alloys - Surface layer


3D MR Flow Analysis in Realistic Rapid-Prototyping Model Systems of the Thoracic Aorta: Comparison with In Vivo Data and Computational Fluid Dynamics in Identical Vessel Geometries

C. Canstein,1 P. Cachot,2 A. Faust,2 A.F. Stalder,1 J. Bock,1 A. Frydrychowicz,1 J. Küffer,2 J. Hennig,1 and M. Markl1

1 Department of Diagnostic Radiology, Medical Physics, University Hospital, Freiburg, Germany.
2 Institute for Product and Production engineering, Department of Technics, University of Applied Sciences Northwestern Switzerland, Muttenz, Switzerland.

Magnetic Resonance in Medicine, 59:535–546, 2008

Abstract: The knowledge of local vascular anatomy and function in the human body is of high interest for the diagnosis and treatment of cardiovascular disease. A comprehensive analysis of the hemodynamics in the thoracic aorta is presented based on the integration of flow-sensitive 4D MRI with state-of-the-art rapid prototyping technology and computational fluid dynamics (CFD). Rapid prototyping was used to transform aortic geometries as measured by contrast-enhanced MR angiography into realistic vascular models with large anatomical coverage. Integration into a flow circuit with patient-specific pulsatile in-flow conditions and application of flow-sensitive 4D MRI permitted detailed analysis of local and global 3D flow dynamics in a realistic vascular geometry. Visualization of characteristic 3D flow patterns and quantitative comparisons of the in vitro experiments with in vivo data and CFD simulations in identical vascular geometries were performed to evaluate the accuracy of vascular model systems. The results indicate the potential of such patient-specific model systems for detailed experimental simulation of realistic vascular hemodynamics. Further studies are warranted to examine the influence of refined boundary conditions of the human circulatory system such as fluid-wall interaction and their effect on normal and pathological blood flow characteristics associated with vascular geometry.

Key words: vascular hemodynamics - phase contrast - aorta - blood flow - velocity mapping - model system - rapid prototyping - computational fluid dynamics


Evaluation of critical fracture energy parameter Gfr and assessment of its transferrability

S. Acharyya1, S. Dhar1, J. Chattopadhyay2

1 Department of Mechanical Engineering, Jadavpur University, Kolkata 700032, India
2 Reactor Safety Division, Bhabha Atomic Research Center, Mumbai 400085, India

Engineering Fracture Mechanics, 75:253–274, 2008

Abstract: Prediction of maximum load bearing capacity and crack growth for ductile materials using existing models like J-R curve approach has the problem of transferability and the use of micro-mechanical model (e.g. Gurson Tvergaard, and Needleman [Tvergaard V, Needleman A. Analysis of cup cone fracture in a round tensile bar. Acta Metall 1984;32:157–169]) are limited by the requirements of the huge computation time and large numbers of critical metallurgical parameters as input to analysis. Marie and Chapuliot [Marie S, Chapuliot S. Ductile tearing simulation based on local energy criterion. Fatigue Fract Engng Mater Struct 1998;21:215–227] of CEA, France, proposed a simple but convenient ductile crack growth model using critical fracture energy (Gfr) for crack growth and Ji for initiation, both of which are material parameters. They also proposed several schemes, namely, graphical and slope of modi.ed plastic J-integral vs crack growth, JM-pl _ Da methods for the evaluation of the value of Gfr from specimens as well as from components. In all these methods the role of non-crack displacement in the crack growth process was not considered. The necessary modifications due to non-crack displacement in the above methods to evaluate the values of Gfr was studied and published [Acharyya S, Dhar S, Chattopadhyay J. (2003). The e.ect of non-crack component on Critical fracture energy on ductile material. Int J Pressure Vessels Piping 2004;81:345–353] by the authors earlier. In this paper, the modi.ed methods and formulation have been applied to evaluate the values of Gfr from experimental and FE simulated results for compact tensile (CT), three point bend (TPB) specimens and also from components like pipes and elbows. Then statistical estimation is done from these Gfr values to assess whether Gfr can be accepted as constant value material parameter. Finally, the mean value of Gfr obtained from statistical computation is used as material constant along with crack initiation toughness parameter (Ji)SZW to consider crack growth for FE simulation of load vs load-line–displacement (LLD) and load vs crack growth curves for different specimens and components. Finite element simulated results are compared with the experimental results and good matching between the two for several components are found and maximum error in prediction of maximum load is found to be within 12%.

Keywords: Ductile fracture - Critical fracture energy (Gfr) - Ernest’s JM-pl - Graphical method - Non-crack compliance - Reverse simulation


Delamination behaviour of very high modulus carbon/epoxy marine composites

N. Barala, P. Daviesb, C. Baleyc, B. Bigourdanb

aTrimaran Groupama, 56100 Lorient, France
bMaterials and Structures Group, IFREMER Brest Centre, 29280 Plouzane, France
cUniversite de Bretagne Sud, L2PIC, BP 92116, 56321 Lorient Cedex, France

Composites Science and Technology 68 (2008) 995–1007

Abstract: This paper presents a study of the delamination resistance of unidirectional carbon .bre reinforced epoxy composites manufactured from prepreg for racing yacht mast applications. Fibre modulus has been varied from 290 to 640 GPa. Two manufacturing methods, oven cure under vacuum and autoclave, were used to produce samples. Results show that delamination resistance under mode I and mixed mode loading decreases as fibre modulus increases. The manufacturing route affects Gc values by modifying interlaminar layer geometry but does not change the trend of lower toughness with increasing fibre modulus. For the pitch composites crack propagation is observed within the fibres, in addition to more usual matrix and interface cracking. For the very high modulus fibre composites it is no longer possible to measure the GI/IIc fracture envelope with the MMB specimen, compression failure occurs before delamination. FE analysis and instrumented specimens have been used to define the limits to the current standard test method for these materials.

Keywords: Delamination - Carbon fibres - Polymers - Fracture toughness - Mixed mode


Influence of microcalcifications on vulnerable plaque mechanics using FSI modeling

Danny Bluesteina, Yared Alemua, Idit Avrahamib,e, Morteza Gharibb, Kris Dumonta, John J. Ricottac, Shmuel Einava,d

aDepartment of Biomedical Engineering, Stony Brook University, Stony Brook, NY 11794-8181, USA
bAeronautics and Bioengineering, California Institute of Technology, Pasadena, CA, USA
cDepartment of Surgery, Stony Brook University Hospital, Stony Brook University, Stony Brook, NY 11794-8181, USA
dDepartment of Bioengineering, Tel Aviv University, Tel Aviv, Israel
eAfeka College of Engineering, Tel Aviv, Israel

Journal of Biomechanics 41 (2008) 1111–1118

Abstract: Sudden heart attacks remain one of the primary causes of premature death in the developed world. Asymptomatic vulnerable plaques that rupture are believed to prompt such fatal heart attacks and strokes. The role of microcalcifications in the vulnerable plaque rupture mechanics is still debated. Recent studies suggest the microcalcifications increase the plaque vulnerability. In this manuscript we present a numerical study of the role of microcalcifications in plaque vulnerability in an eccentric stenosis model using a transient fluid–structure
interaction (FSI) analysis. Two cases are being compared (i) in the absence of a microcalcification (ii) with a microcalcification spot fully embedded in the fibrous cap. Critical plaque stress/strain conditions were affected considerably by the presence of a calcified spot, and were dependent on the timing (phase) during the flow cycle. The vulnerable plaque with the embedded calcification spot presented higher
wall stress concentration region in the fibrous cap a bit upstream to the calcified spot, with stress propagating to the deformable parts of the structure around the calcified spot. Following previous studies, this finding supports the hypothesis that microcalcifications increase the plaque vulnerability. Further studies in which the effect of additional microcalcifications and parametric studies of critical plaque cap
thickness based on plaque properties and thickness, will help to establish the mechanism by which microcalcifications weaken the plaque and may lead to its rupture.

Keywords: Vulnerable plaque - Atherosclerosis - Fibrous cap - Rupture - Microcalcification - Numerical modeling - FSI


Non-Newtonian Bile Flow in Elastic Cystic Duct: One- and Three-Dimensional Modeling

W.G. Li1, X.Y. Luo1,S.B. Chin2,N.A. Hill1, A.G. Johnson3, and N.C. Bird3

1Department of Mathematics, University of Glasgow, Glasgow G12 8QW, UK;
2Department of Mechanical Engineering, University of Sheffield, Sheffield S1 3JD, UK; and
3Academic Surgical Unit, Royal Hallamshire Hospital, Sheffield S10 2JF, UK

Annals of Biomedical Engineering, Vol. 36, No. 11, pp. 1893–1908, 2008

Abstract: Bile flow is thought to play an essential role in the pathophysiological genesis of cholelithiasis (gallstone formation) and in gallbladder pain. In this paper, we extend our previous study of the human biliary system (Li et al., 2007, J. Biomech. Eng., 129:164–173) to include two important factors: the non-Newtonian properties of bile, and elastic deformation of the cystic duct. A one-dimensional (1D) model is analyzed and compared with three-dimensional (3D) fluid–structure interaction simulations. It is found that non-Newtonian bile raises resistance to the flow of bile, which can be augmented significantly by the elastic deformation (collapse) of the cystic duct. We also show that the 1D model predicts the pressure drop of the cystic duct flow well for all cases considered (Newtonian or non-Newtonian flow, rigid or elastic ducts), when compared with the full 3D simulations.

Keywords: Gallbladder – Gallstone - Cystic duct - Biliary system - Non-Newtonian fluid - Pressure drop - Fluid–structure interaction - 1D modeling - Finite element methods


Cyclic Bending Contributes to High Stress in a Human Coronary Atherosclerotic Plaque and Rupture Risk: In Vitro Experimental Modeling and Ex Vivo MRI-Based Computational Modeling Approach

Chun Yang1,2, Dalin Tang1, Shunichi Kobayashi3, Jie Zheng4, Pamela K. Woodard4, Zhongzhao Teng1, Richard Bach5, and David N. Ku6

1Mathematical Sciences Department, Worcester Polytechnic Institute, Worcester, MA 01609
2School of Mathematical Sciences, Beijing Normal University, Beijing, China
3Division of Creative Engineering, Shinshu University, Ueda, Nagano, Japan
4Mallinkcrodt Inst. of Radiology, Washington University, St. Louis, MO 63110, USA
5Division of Cardiovascular Diseases, Washington University, St. Louis, MO 63110, USA
6School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, GA, 30332-0405, USA

Mol Cell Biomech. 2008; 5(4): 259–274.

Abstract: Many acute cardiovascular syndromes such as heart attack and stroke are caused by atherosclerotic plaque ruptures which often happen without warning. MRI-based models with fluid-structure interactions (FSI) have been introduced to perform flow and stress/strain analysis for atherosclerotic plaques and identify possible mechanical and morphological indices for accurate plaque vulnerability assessment. In this paper, cyclic bending was added to 3D FSI coronary plaque models for more accurate mechanical predictions. Curvature variation was prescribed using the data of a human left anterior descending (LAD) coronary artery. Five computational models were constructed based on ex vivo MRI human coronary plaque data to assess the effects of cyclic bending, pulsating pressure, plaque structure, and axial stretch on plaque stress/strain distributions. In vitro experiments using a hydrogel stenosis model with cyclical bending were performed to observe effect of cyclical bending on flow conditions. Our results indicate that cyclical bending may cause more than 100% or even up to more than 1000% increase in maximum principal stress values at locations where the plaque is bent most. Stress increase is higher when bending is coupled with axial stretch, non-smooth plaque structure, or resonant pressure conditions (zero phase angle shift). Effects of cyclic bending on flow behaviors are more modest (21.6% decrease in maximum velocity, 10.8% decrease in flow rate, maximum flow shear stress changes were < 5%). Computational FSI models including cyclic bending, plaque components and structure, axial stretch, accurate in vivo measurements of pressure, curvature, and material properties should lead to significant improvement on stress-based plaque mechanical analysis and more accurate coronary plaque vulnerability assessment.

Keywords: Coronary artery - cardiovascular - cyclic bending - fluid-structure interaction - blood flow - atherosclerotic plaque rupture


Computational studies of resonance wave pumping in compliant tube

Idit Avrahami1 And Morteza Gharib2

1Medical Engineering, AFEKA-Tel Aviv. Academic College of Engineering, Bney Efraim, Tel Aviv 69107, Israel
2Aeronautics and Bioengineering, California Institute of Technology, Pasadena, CA 91125, USA

J. Fluid Mech. (2008), vol. 608, pp. 139–160.

Abstract: The valveless impedance pump is a simple design that allows the producion or amplification of a flow without the requirement for valves or impellers. It is based on fluid-filled flexible tubing, connected to tubing of different impedances. Pumping is achieved by a periodic excitation at an off-centre position relative to the tube ends. This paper presents a comprehensive study of the fluid and structural dynamics in an impedance pump model using numerical simulations. An axisymmetric finite element model of both the fluid and solid domains is used with direct coupling at the interface. By examining a wide range of parameters, the pump’s resonance nature is described and the concept of resonance wave pumping is discussed. The main driving mechanism of the flow in the tube is the reflection of waves at the tube boundary and the wave dynamics in the passive tube. This concept is supported by three different analyses: (i) time-dependent pressure and flow wave dynamics along the tube, (ii) calculations of pressure–flow loop areas along the passive tube for a description of energy conversion, and (iii) an integral description of total work done by the pump on the fluid. It is shown that at some frequencies, the energy given to the system by the excitation is converted by the elastic tube to kinetic energy at the tube outlet, resulting in an efficient pumping mechanism and thus signi.cantly higher flow rate. It is also shown that pumping can be achieved with any impedance mismatch at one boundary and that the outlet con.guration does not necessarily need to be a tube.


Analysis of the heating process and development of a microstructure suitable for thixoforming of steel

K.P. Solek, A. Lukaszek-Solek
AGH University of Science and Technology – al. Mickiewicza 30, 30-059 Kraków, Poland

Int J Mater Form (2008) Suppl 1:1015–1018

ABSTRACT: The thixoforming processes of metal alloys are carried out in the semi-solid state. The thixotropy phenomenon is conditional upon the globular microstructure that occurs in shaped material. In order to design and optimise such processes, some procedures of the Gleeble® system could be applied. The main objective of this work is an analysis of the resistance heating of samples in different material tests which are possible on the Gleeble® 3800 simulator. An example of development of steel alloy microstructure during heating on the Gleeble® simulator will also be presented. M2 tool steel was submitted for analysis. In this case, obtaining the globular microstructure was carried out on the basis of the SIMA method (Strain Induced Melt Activated).

Key words: thixoforming - resistance heating - globular microstructure - numerical analysis


Identification of nonlinear dynamic coefficients in plain journal bearings

V. Meruane, R. Pascual

Mechanical Engineering Department, Universidad de Chile, Beauchef 850, Santiago, Chile

Tribology International 41 (2008) 743–754

Abstract: This work proposes a framework to the numerical identification of nonlinear fluid film bearing parameters from large journal orbital motion (20–60% of the bearing clearance). Nonlinear coefficients are defined by a third order Taylor expansion of bearing reaction forces and are evaluated through a least mean square in time domain technique. The journal response is obtained from a computational fluid dynamic (CFD) model of a plain journal bearing on high dynamic loading conditions. The model considers fluid–structure interaction between the fluid flow and the journal. The case in study considers a laboratory test rig. Results indicate that nonlinear coefficients have an important effect on stiffness and damping. It was found a change on nonlinear behavior occurred when the Oil Whirl phenomenon starts, which it is not seen in classical linear models.

Keywords: Journal bearing - Computational fluid dynamics - Transient - Identification - Large orbital motion - Nonlinear


Using traditional and innovative approaches in the seismic evaluation and strengthening of a historic unreinforced masonry synagogue

Terrence F. Paret, Sigmund A. Freeman, Gary R. Searer, Mahmoud Hachem, Una M. Gilmartin

Wiss, Janney, Elstner Associates, Inc., 2200 Powell Street, Suite 925, Emeryville, CA 94608, United States

Engineering Structures 30 (2008) 2114–2126

Abstract: A monumental 100-year-old multi-story brick masonry synagogue in San Francisco was threatened with closure due to non-compliance with an Unreinforced Masonry Building Ordinance. The Ordinance required demonstrating that the building would protect life safety in the event of a major earthquake or upgrading the building to do so, despite the historical record showing that the building survived the Great 1906 Earthquake with relatively little damage. Much of the minor damage that did occur is still visible today as the historic finishes have remained intact since that
earthquake. The structure consists primarily of thick, brick masonry perimeter walls, wood-framed diaphragms with diagonal and straight sheathing, and riveted structural steel trusses and columns that support the dome over the main sanctuary and a large portion of the floors and main roof. The structure was subjected to linear dynamic and non-linear static analyses to benchmark its behavior during the 1906 earthquake, to better understand its vulnerabilities in the event of potentially stronger shaking, and to develop seismic improvements that minimize disturbance of its ornately detailed stone and painted plaster finishes. The most appropriate structural solution consistent with preserving the historic fabric takes advantage of the dynamic separation between the modes predominated by in-plane and out-of-plane wall shaking. The out-of-plane wall modes have long periods and low modal mass participation. The in-plane wall modes have higher mass participation, but act rigidly with little amplification. The solution consisted of a combination of intervention techniques, each developed to minimize disturbance to the non-structural finishes and retain the original dynamic characteristics. These interventions include a system of tension ties in the attic that interconnect the four perimeter walls, yet circumvent the domed sanctuary; center-cored reinforcement of the masonry walls; and fiber-wrap of a few critical piers. The tension ties contain super-elastic nitinol wires and were designed to be lightweight, easy to install, and to restrain the walls from falling outward while maintaining the modal separation and inherent flexibility of the system that enabled the structure to survive the 1906 earthquake.

Keywords: Unreinforced masonry building - Seismic strengthening - Historic structure


Seismic vulnerability of historical masonry buildings: A case study in Ferrara

Vincenzo Mallardoa, Roberto Malvezzia, Enrico Milanib, Gabriele Milanib

aDepartment of Architecture, University of Ferrara, Via Quartieri 8, 44100 Ferrara, Italy
bDepartment of Engineering, University of Ferrara, Via Saragat 1, 44100 Ferrara, Italy

Engineering Structures 30 (2008) 2223–2241

Abstract: In this paper, the seismic behaviour of an important Renaissance Palace in Ferrara (Italy), Palazzo Renata di Francia, is investigated. In the first part, a full 3D nonlinear analysis is performed in order to obtain a first insight into its structural performance. Then, three different 2D nonlinear analyses of the main front are critically compared and discussed. The first model relies on an equivalent frame approach, the second is based on an elastic–plastic damaging Finite Element (F.E.) approach and the third stems from a kinematic collapse analysis by means of a homogenised limit surface. The facade exhibits an adequate in-plane strength, when analysed with 2D plane stress F.Es., whereas the whole building undergoes out-of-plane partial collapses even in the presence of very low intensity horizontal actions, when analysed with 3D F.Es.

Keywords: Masonry - Nonlinear - Seismic - Homogenisation


Influence of Meniscectomy and Meniscus Replacement on the Stress Distribution in Human Knee Joint

Ashkan Vaziri,1,2 Hamid Nayeb-Hashemi,2 Arvinder Singh,2 and Bashir A. Tafti3

1School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA;
2Department of Mechanical, Industrial and Manufacturing Engineering, Northeastern University, Boston, MA 02115, USA; and
3Department of Surgery and Regenerative Medicine, Division of Plastic Surgery, Stanford University, Stanford, CA 94305, USA

Annals of Biomedical Engineering, Vol. 36, No. 8, August 2008, pp. 1335–1344

Abstract: Studying the mechanics of the knee joint has direct implications in understanding the state of human health and disease and can aid in treatment of injuries. In this work, we developed an axisymmetric model of the human knee joint using finite element method, which consisted of separate parts representing tibia, meniscus and femoral, and tibial articular cartilages. The articular cartilages were modeled as three separate layers with different material characteristics: top superficial layer, middle layer, and calcified layer. The biphasic characteristic of both meniscus and cartilage layers were included in the computational model. The developed model was employed to investigate several aspects of mechanical response of the knee joint under external loading associated with the standing posture. Specifically, we studied the role of the material characteristic of the articular cartilage and meniscus on the distribution of the shear stresses in the healthy knee joint and the knee joint after meniscectomy. We further employed the proposed computational model to study the mechanics of the knee joint with an arti.cial meniscus. Our calculations suggested an optimal elastic modulus of about 110 MPa for the arti.cial meniscus which was modeled as a linear isotropic material. The suggested optimum stiffness of the artificial meniscus corresponds to the stiffness of the physiological meniscus in the circumferential direction.

Keywords: Knee mechanics - Meniscectomy - Artificial meniscus


Collapse arrestors for deepwater pipelines. Cross-over mechanisms

Rita G. Toscano, Luciano O. Mantovano, Pablo M. Amenta, Roberto F. Charreau, Daniel H. Johnson, Andrea P. Assanelli, Eduardo N. Dvorkin

Center for Industrial Research, TENARIS, Dr. Jorge A. Simini 250, Campana 2804, Argentina

Computers and Structures 86 (2008) 728–743

Abstract: Using finite element models it is possible to determine the cross-over external pressure of different pipeline arrestor designs. In this paper these finite element models are discussed and validated by comparing their results with experimental determinations. The flipping and flattening cross-over mechanisms, that were previously described in the literature, are considered in the experimental validation of the numerical models.

Keywords: Collapse - Pipelines - Deepwater - Finite elements


Using 3D fluid–structure interaction model to analyse the biomechanical properties of erythrocyte

C.Y. Cheea, H.P. Leea,b, C. Lua

aInstitute of High Performance Computing, 1 Science Park Road, Capricorn S117528, Singapore
bDepartment of Mechanical Engineering, National University of Singapore, Singapore 119260, Singapore

Physics Letters A 372 (2008) 1357–1362

Abstract: This Letter presents a newly developed three-dimensional fluid–structure interaction model of the red blood cell (RBC). The model consists of a deformable liquid capsule modelled as Newtonian fluid enclosed by a hyperelastic membrane with viscoelastic property. Numerical results show that viscosity in the cytoplasm affects the deformed shape of RBC under loading. This observation is contrary to the earlier belief that viscosity of the cytoplasm can be neglected. Numerical simulations carried out to investigate large deformation induced on the RBC model using direct tensile forces show significant improvement in terms of correlation with experimental results. The membrane shear modulus estimated from the model ranges between 3.7 to 9.0 µNm-1 compares well with results obtained from micropipette aspiration experiments.

Keywords:Cells biomechanics - Computational fluid–structure interaction


Numerical analyses of dynamic response of saturated porous seabed-pipeline interaction under seismic loading

Xiaoling Zhang, Maotian Luan

Institute of Geotechnical Engineering, School ofCivil and Hydraulic Engineering, Dalian University ofTechnology, Dalian JJ6023, China

Geotechnical Engineering for Disaster Mitigation and Rehabilitation. Liu, Deng and Chu (eds) 2008 Science Press Beijing and Springer-Verlag GmbH Berlin Heidelberg

Abstract: Based on the Biot's theory of consolidation, the model of the seabed-pipeline interaction is established using viscous-elastic artificial boundary in this paper. The distribution ofthe pore pressure in seabed soil and the dynamic response of submarine pipeline are studied under El Centro seismic wave by virtue of the general-purpose FEM analysis package ADINA. The effects of the elastic anisotropic soil and two-layered seabed behaviors on the seismic-induced pore pressure and internal stresses are discussed emphatically. It actually shows the validity and practicality ofdynamic response ofsubmarine pipeline under seismic loading.


Mechanical behavior of segment rebar of shield tunnel in construction stage

Jun-sheng Chen, Hai-hong Mo

Department of Civil Engineering, South China University of Technology, Guangzhou 510640, China
State Key Laboratory of Subtropical Building Science, South China University of Technology, Guangzhou 510640, China

Journal of Zhejiang University SCIENCE A, 2008 9(7):888-899

Abstract: In this paper, a 3D finite element (FE) program ADINA was applied to analyzing a tunnel with 9 segment rings. The loads acting on these segment rings included the squeezing action of tail brush of shield machine under attitude deflection, the jacking forces, the grouting pressure and the soil pressure. The analyses focused on the rebar stress in two statuses: (1) normal construction status without shield machine squeezing; (2) squeezing action induced by shield machine under attitude deflection. The analyses indicated that the rebar stress was evidently affected by the construction loads. In different construction status, the rebar stress ranges from -80 MPa to 50 MPa, and the rebar is in elastic status. Even some cracks appear on segments, the stress of segment rebar is still at a low level. It is helpful to incorporate a certain quantity of steel fiber to improve the anti-crack and shock resistance performance.

Key words: Shield tunnel - Construction - Segment - Rebar - Finite element (FE) analysis


Tensile radial stress in the spinal cord related to arachnoiditis or tethering: a numerical model

C. D. Bertram1, L. E. Bilston2 and M. A. Stoodley2

1Biofluid Mechanics Laboratory, Faculty of Engineering, University of New South Wales, Sydney 2052, Australia
2Prince of Wales Medical Research Institute, UNSW, Barker Street, Randwick 2031, New South Wales, Australia

Med Biol Eng Comput (2008) 46:701–707

Abstract: Spinal arachnoiditis comprises fibrous scarring of the subarachnoid space, following spinal trauma or inflammation, and is often associated with syringomyelia. We hypothesised that cord-to-dura attachments could cause transient tensile cord radial stress, as pressure waves propagate. This was tested in a fluid–structure interaction model, simulating three types of cord tethering, with ‘arachnoiditis’ confined to a short mid-section of the cord. The annular system was excited abdominally with a short transient, and the resulting Young and Lamb waves and reflections were analysed. Radial mid-section tethering was less significant than axial tethering, which gave rise to tensile radial stress locally when the cord was not fixed cranially. Simulated as inextensible string connections to the dura, arachnoiditis caused both localised tensile radial stress and localised low pressure in the cord as the transient passed. The extent of these effects was sensitive to the relative stiffness of the dura and cord. Tensile radial stress may create a syrinx in previously normal cord tissue, and transiently lowered pressure may draw in interstitial fluid, causing the syrinx to enlarge if fluid exit is inhibited. The suggested mechanism could also explain the juxtaposition of syrinxes to regions of arachnoiditis.

Keywords: Fluid–structure interaction - Finite-element model - Cerebrospinal fluid - Wave propagation - Syringomyelia


Fluid–solid interaction simulation of flow and stress pattern in thoracoabdominal aneurysms: A patient-specific study

A. Borghia, N.B. Wooda, R.H. Mohiaddinb, X.Y. Xua

aDepartment of Chemical Engineering, South Kensington Campus, Imperial College, London SW72AZ, UK
bRoyal Brompton and Harefield NHS Trust, Sydney Street, London, UK

Journal of Fluids and Structures 24 (2008) 270–280

Abstract: Thoracoabdominal aneurysm (TA) is a pathology that involves the enlargement of the aortic diameter in the inferior descending thoracic aorta and has risk factors including aortic dissection, aortitis or connective tissue disorders. Abnormal flow patterns and haemodynamic stress on the diseased aortic wall are thought to play an important role in the development of this pathology and the internal wall stress has proved to be more reliable as a predictor of rupture than the maximum diameter for abdominal aortic aneurysms; but this assumption has not been validated yet for aneurysms involving the thoracic aorta. In the present study, three patients with TAs of different maximum diameters were scanned using magnetic resonance imaging (MRI) techniques. Realistic models of the aneurysms were reconstructed from the in vivo MRI data acquired from the patients, and subject-specific flow conditions were applied as boundary conditions. The wall and thrombus were modelled as hyperelastic materials and their properties were derived from the literature. A normal descending aorta was also simulated to provide data for comparison. Fully coupled fluid–solid interaction (FSI) simulations as well as solid static simulations were performed using ADINA 8.2. The results show that the wall stress distribution and its magnitude are strongly dependent on the 3-D shape of the aneurysm and the distribution of thrombus. Maximum wall stresses in all TA models are higher than in the normal aorta, and values of maximum wall stress are not directly related to the maximum aneurysm diameter. Comparisons between the FSI and solid static simulation results showed no significant difference in maximum wall stress, supporting those previous studies which found that FSI simulations were not  necessary for wall stress prediction.


Frictional Heat-Induced Phase Transformation on Train Wheel Surface

Su Hang1, Pan Tao1, Li Lit1, Yang Cai-fu1, Cui Yin-hui2, Ji Huai-zhong2

1 Division of Structural Material Research, Central Iron and Steel Research Institute, Beijing 100081, China
2 Ma' anshan Iron and Steel Company Limited, Ma' anshan 243000, Anhui, China

Journal Of Iron And Steel Research, International, 15(5): 49-55, 2008.

Abstract: By combining thermomechanical coupling finite element analysis with the characteristics of phase transformation [continuous cooling transformation (CCT) curve], the thermal fatigue behavior of train wheel steel under high speed and heavy load conditions was analyzed. The influence of different materials on the formation of the phase transformation zone of the wheel tread was discussed. The result showed that the peak temperature of wheel/track friction zone could be higher than the austenitizing temperature for braking. The depth of the austenitized region could reach a point of 0.9 mm beneath the wheel tread surface. The supercooled austenite is transformed to a hard and brittle martensite layer during the following rapid cooling process, which may lead to cracking and then spalling on the wheel tread surface. The decrease in carbon contents of the train wheel steel helps inhibit the formation of martensite by increasing the austenitizing temperature of the train wheel steel. When the carbon contents decrease from 0.7% to 0.4 %, the Ac3 of the wheel steel is increased by 45°C, and the thickness of the martensite layer is decreased by 30%, which is helpful in reducing the thermal cycling fatigue of the train wheel tread such as spalling.

Keywords: train wheel steel - thermal cycling fatigue - friction - martensite transformation -  thermomechanical coupling


The Carrying Capacity of Conical Interference Fit Joints with Laser Reinforcement Zones

Lucjan Sniezek2, Jolanta Zimmerman1, Andrzej Zimmerman1

1 Warsaw University of Technology, Faculty of Production Engineering, Institute of Mechanics and Printing, 85 Narbutta St., 02-524 Warsaw, Poland
2 Military University of Technology, Faculty of Mechanical Engineering, Department of Machine Design, 2 Kaliskiego St., 00-908 Warsaw, Poland

Journal of Materials Processing Technology, doi:10.1016/j.jmatprotec.2010.02.004, 2008.

The subject of the paper is the comparative research on the carrying capacity of the conical interference fit joints and the interference fit joints with the laser reinforcement paths. The calculations of strain and stress distribution in the shaft neck and the bush were carried out for three stages: during the joining process, after the interference and after the application of torsional moment. The FEM was used. The dependences of the joint carrying capacity on the interference force were determined and experimental tests of the capacity of the conical interference fit joints were carried out. The numerical calculation showed the increase in the capacity of the conical interference fit joints with the laser reinforcement zones in comparison with the standard joints by 50% in the case of the same axial displacement. Greater capability (about 25%) of this joint was achieved from experimental tests in comparison to one without

Keywords: conical interference fit joints - stress distribution - torque capacity - FEM

Estimation of NonlinearMechanical Properties of Vascular Tissues via Elastography

Reza Karimi1,2, Ting Zhu1, Brett E. Bouma3, Mohammad R. Kaazempur Mofrad4

1 Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
2 Department of Bioengineering, University of California, 208A Stanley Hall #1762, Berkeley, CA 94720, USA
3 Wellman Center for Photomedicine, Massachusetts General Hospital, Harvard Medical School,
Boston, MA, USA 02114
4 Department of Bioengineering, University of California, 208A Stanley Hall #1762, Berkeley, CA 94720, USA, e-mail:

Cardiovasc Eng., 8(4):191–202, 2008.

Abstract: A new method is proposed for estimation of nonlinear elastic properties of soft tissues. The proposed approach involves a combination of nonlinear finite element methods with a genetic algorithm for estimating tissue stiffness profile. A multipoint scheme is introduced that satisfies the uniqueness condition, improves the estimation performance, and reduces the sensitivity to image noise. The utility of the proposed techniques is demonstrated using optical coherence tomography (OCT) images. The approach is, however, applicable to other imaging systems and modalities, as well, provided a reliable image registration scheme. The proposed algorithm is applied to realistic (2D) and idealized (3D) arterial plaque models, and proves promising for the estimation of intra-plaque distribution of nonlinear material properties.

Keywords: Soft tissue mechanics - Nonlinear material properties – Atherosclerosis - Elastography - Genetic algorithm - Optical coherence tomography

Nonspherical Half-Shells by Ultrasonic Cleavage of the Hollow Polyhedral Particles in Water

Shuhong Jiao1, Kai Jiang2, Yanhong Zhang3, Ming Xiao1, Lifen Xu1, and Dongsheng Xu1

1 Beijing National Laboratory for Molecular Science (BNLMS), State Key Laboratory for Structural Chemistry of Unstable and Stable Species, College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
2 College of Chemistry and Environmental Science, Henan Normal University, Xinxiang, Henan 453002, China
3 China Institute of Water Resources and Hydropower Research, Beijing 100044, China

J. Phys. Chem. C, 112:3358-3361, 2008.

Abstract: We present a facile route to preparing nonspherical half-shells of copper sulfide in which the octahedral mesocages of copper sulfide were grown on the substrate by a two-step electrochemical method first and then were cleaved into tetrahedral half-shells by ultrasonic vibration in water. The obtained half-shells are 1-3 µm in lengths, and the shell thickness is uniform in ˜100 nm. It is found that only centrosymmetric octahedral and six-pod shells can be equally split into two similar half-shells. Furthermore, a finite element analysis has been used to model the stress distribution in the octahedral CuS shell structure during ultrasonic vibration.


(latest on Page 43)
Page 43 Page 42 Page 41 Page 40 Page 39 Page 38 Page 37 Page 36 Page 35 Page 34
Page 33 Page 32 Page 31 Page 30 Page 29 Page 28 Page 27 Page 26 Page 25 Page 24
Page 23 Page 22 Page 21 Page 20 Page 19 Page 18 Page 17 Page 16 Page 15 Page 14
Page 13 Page 12 Page 11 Page 10 Page 9  Page 8  Page 7  Page 6  Page 5  Page 4 
Page 3  Page 2  Page 1 

Page 30

ADINA logo