Abstract: The mechanical behavior of porous media is governed by the interaction between its solid skeleton and the fluid existing inside its pores. The interaction occurs through the interface of gains and fluid. The traditional analysis methods of porous media, based on the effective stress and Darcy's law, are unable to account for these interactions. For an accurate analysis, the porous media is represented in a fluid-filled porous solid on the basis of the Biot theory of wave propagation in poroelastic media. In Biot formulation, the equations of motion of the soil mixture are coupled with the global mass balance equations to describe the realistic behavior of porous media. Because of irregular geometry, the domain is generally treated as an assemblage of fmite elements. In this investigation, the numerical formulation for the field equations governing the dynamic response of fluid-saturated porous media is analyzed and employed for the study of transient wave motion. A finite element model is developed and implemented into a computer code called DYNAPM for dynamic analysis of porous media. The weighted residual method with 8-node elements is used for developing of a finite element model and the analysis is carried out in the time domain considering the dynamic excitation and gravity loading. Newmark time integration scheme is developed to solve the time-discretized equations which are an unconditionally stable implicit method Finally, some numerical examples are presented to show the accuracy and capability of developed model for a wide variety of behaviors of porous media.
Abstract: In this study the elastic-plastic stress distribution in
weld-bonded joint, fabricated from austenitic stainless steel (AISI
304) sheet of 1.00 mm thickness and Epoxy adhesive Araldite 2011,
subjected to axial loading is investigated. This is needed to improve
design procedures and welding codes, and saving efforts in the
cumbersome experiments and analysis. Therefore, a complete 3-D
finite element modelling and analysis of spot welded, bonded and
weld-bonded joints under axial loading conditions is carried out. A
comprehensive systematic experimental program is conducted to
determine many properties and quantities, of the base metals and the
adhesive, needed for FE modelling, such like the elastic – plastic
properties, modulus of elasticity, fracture limit, the nugget and heat
affected zones (HAZ) properties, etc. Consequently, the finite
element models developed, for each case, are used to evaluate
stresses distributions across the entire joint, in both the elastic and
plastic regions. The stress distribution curves are obtained,
particularly in the elastic regions and found to be consistent and in
excellent agreement with the published data. Furthermore, the
stresses distributions are obtained in the weld-bonded joint and
display the best results with almost uniform smooth distribution
compared to spot and bonded cases. The stress concentration peaks at
the edges of the weld-bonded region, are almost eliminated resulting
in achieving the strongest joint of all processes.
Abstract: Stress analysis of functionally graded composite plates
composed of ceramic, functionally graded material and metal layers is
investigated using 3-D finite element method. In FGM layer, material
properties are assumed to be varied continuously in the thickness
direction according to a simple power law distribution in terms of the
volume fraction of a ceramic and metal. The 3-D finite element model
is adopted by using an 18-node solid element to analyze more
accurately the variation of material properties in the thickness
direction. Numerical results are compared for three types of materials.
In the analysis, the tensile and the compressive stresses are
summarized for various FGM thickness ratios, volume fraction
distributions, geometric parameters and mechanical loads.
Abstract: In the present work, behavior of inoxydable steel as
reinforcement bar in composite concrete is being investigated. The
bar-concrete adherence in reinforced concrete (RC) beam is studied
and focus is made on the tension stiffening parameter. This study
highlighted an approach to observe this interaction behavior in
bending test instead of direct tension as per reported in many
references. The approach resembles actual loading condition of the
structural RC beam. The tension stiffening properties are then
applied to numerical finite element analysis (FEA) to verify their
correlation with laboratory results. Comparison with laboratory
shows a good correlation between the two. The experimental settings
is able to determine tension stiffening parameters in RC beam and
the modeling strategies made in ABAQUS can closely represent the
actual condition. Tension stiffening model used can represent the
interaction properties between inoxydable steel and concrete.
Abstract: Human middle-ear is the key component of the
auditory system. Its function is to transfer the sound waves through
the ear canal to provide sufficient stimulus to the fluids of the inner
ear. Degradation of the ossicles that transmit these sound waves from
the eardrum to the inner ear leads to hearing loss. This problem can
be overcome by replacing one or more of these ossicles by middleear
prosthesis. Designing such prosthesis requires a comprehensive
knowledge of the biomechanics of the middle-ear. There are many
finite element modeling approaches developed to understand the
biomechanics of the middle ear. The available models in the
literature, involve high computation time. In this paper, we propose a
simplified model which provides a reasonably accurate result with
much less computational time. Simulation results indicate a
maximum sound pressure gain of 10 dB at 5500 Hz.
Abstract: To date, theoretical studies concerning the Carbon
Fiber Reinforced Polymer (CFRP) strengthening of RC beams with
openings have been rather limited. In addition, various numerical
analyses presented so far have effectively simulated the behaviour of
solid beam strengthened by FRP material. In this paper, a two
dimensional nonlinear finite element analysis is presented to validate
against the laboratory test results of six RC beams. All beams had the
same rectangular cross-section geometry and were loaded under four
point bending. The crack pattern results of the finite element model
show good agreement with the crack pattern of the experimental
beams. The load midspan deflection curves of the finite element
models exhibited a stiffer result compared to the experimental beams.
The possible reason may be due to the perfect bond assumption used
between the concrete and steel reinforcement.
Abstract: The present paper presents a finite element model and
analysis for the interaction between a piezoresistive tactile sensor and
biological tissues. The tactile sensor is proposed for use in minimally
invasive surgery to deliver tactile information of biological tissues to
surgeons. The proposed sensor measures the relative hardness of soft
contact objects as well as the contact force. Silicone rubbers were
used as the phantom of biological tissues. Finite element analysis of
the silicone rubbers and the mechanical structure of the sensor were
performed using COMSOL Multiphysics (v3.4) environment. The
simulation results verify the capability of the sensor to be used to
differentiate between different kinds of silicone rubber materials.
Abstract: A transient finite element model has been developed
to study the heat transfer and fluid flow during spot Gas Tungsten
Arc Welding (GTAW) on stainless steel. Temperature field, fluid
velocity and electromagnetic fields are computed inside the cathode,
arc-plasma and anode using a unified MHD formulation. The
developed model is then used to study the influence of different
helium-argon gas mixtures on both the energy transferred to the
workpiece and the time evolution of the weld pool dimensions. It is
found that the addition of helium to argon increases the heat flux
density on the weld axis by a factor that can reach 6.5. This induces
an increase in the weld pool depth by a factor of 3. It is also found
that the addition of only 10% of argon to helium decreases
considerably the weld pool depth, which is due to the electrical
conductivity of the mixture that increases significantly when argon is
added to helium.
Abstract: In this study the behavior of interlaminar fracture of
carbon-epoxy thermoplastic laminated composite is investigated
numerically and experimentally. Tests are performed with Arcan
specimens. Testing with Arcan specimen gives the opportunity of
utilizing just one kind of specimen for extracting fracture properties
for mode I, mode II and different mixed mode ratios of materials with
exerting load via different loading angles. Variation of loading angles
in range of 0-90° made possible to achieve different mixed mode
ratios. Correction factors for various conditions are obtained from
ABAQUS 2D finite element models which demonstrate the finite
shape of Arcan specimens used in this study. Finally, applying the
correction factors to critical loads obtained experimentally, critical
interlaminar fracture toughness of this type of carbon- epoxy
composite has been attained.
Abstract: In the present research, a finite element model is
presented to study the geometrical and material nonlinear behavior of
reinforced concrete plane frames considering soil-structure
interaction. The nonlinear behaviors of concrete and reinforcing steel
are considered both in compression and tension up to failure. The
model takes account also for the number, diameter, and distribution
of rebar along every cross section. Soil behavior is taken into
consideration using four different models; namely: linear-, nonlinear
Winkler's model, and linear-, nonlinear continuum model. A
computer program (NARC) is specially developed in order to
perform the analysis. The results achieved by the present model show
good agreement with both theoretical and experimental published
literature. The nonlinear behavior of a rectangular frame resting on
soft soil up to failure using the proposed model is introduced for
demonstration.
Abstract: Development, calibration and validation of a threedimensional
model of the Legform impactor for pedestrian crash with
bumper are presented. Lower limb injury is becoming an increasingly
important concern in vehicle safety for both occupants and
pedestrians. In order to prevent lower extremity injuries to a
pedestrian when struck by a car, it is important to elucidate the
loadings from car front structures on the lower extremities and the
injury mechanism caused by these loadings. An impact test
procedure with a legform addressing lower limb injuries in car
pedestrian accidents has been proposed by EEVC/WG17. In this
study a modified legform impactor is introduced and validated
against EEVC/WG17 criteria. The finite element model of this
legform is developed using LS-DYNA software. Total mass of
legform impactor is 13.4 kg.Technical specifications including the
mass and location of the center of gravity and moment of inertia
about a horizontal axis through the respective centre of gravity in
femur and tibia are determined. The obtained results of legform
impactor static and dynamic tests are as specified in the
EEVC/WG17.
Abstract: In the analysis of structures, the nonlinear effects due to large displacement, large rotation and materially-nonlinear are very important and must be considered for the reliable analysis. The non-linear fmite element analysis has potential as usable and reliable means for analyzing of civil structures with the availability of computer technology. In this research the large displacements and materially nonlinear behavior of shear wall is presented with developing of fmite element code using the standard Galerkin weighted residual formulation. Two-dimensional plane stress model was carried out to present the shear wall response. Total Lagangian formulation, which is computationally more effective, is used in the formulation of stiffness matrices and the Newton-Raphson method is applied for the solution of nonlinear transient equations. The details of the program formulation are highlighted and the results of the analyses are presented, along with a comparison of the response of the structure with Ansys software results. The presented model in this paper can be developed for nonlinear analysis of civil engineering structures with different material behavior and complicated geometry.
Abstract: Resistance spot welding process comprises of electric,
thermal and mechanical phenomenon, which makes this process
complex and highly non-linear and thus, it becomes difficult to model
it. In order to obtain good weld nugget during spot welding, hit and
trial welds are usually done which is very costly. Therefore the
numerical simulation research has been conducted to understand the
whole process. In this paper three different cases were analyzed by
varying the tip contact area and it was observed that, with the
variation of tip contact area the nugget formation at the faying
surface is affected. The tip contact area of the welding electrode
becomes large with long welding cycles. Therefore in order to
maintain consistency of nugget formation during the welding process,
the current compensation in control feedback is required. If the
contact area of the welding electrode tip is reduced, a large amount of
current flows through the faying surface, as a result of which
sputtering occurs.
Abstract: Reliable information about tool temperature
distribution is of central importance in metal cutting. In this study,
tool-chip interface temperature was determined in cutting of ST37
steel workpiece by applying HSS as the cutting tool in dry turning.
Two different approaches were implemented for temperature
measuring: an embedded thermocouple (RTD) in to the cutting tool
and infrared (IR) camera. Comparisons are made between
experimental data and results of MSC.SuperForm and FLUENT
software.
An investigation of heat generation in cutting tool was performed
by varying cutting parameters at the stable cutting tool geometry and
results were saved in a computer; then the diagrams of tool
temperature vs. various cutting parameters were obtained. The
experimental results reveal that the main factors of the increasing
cutting temperature are cutting speed (V ), feed rate ( S ) and depth
of cut ( h ), respectively. It was also determined that simultaneously
change in cutting speed and feed rate has the maximum effect on
increasing cutting temperature.
Abstract: The integral form of equations of motion of composite
beams subjected to varying time loads are discretized using a
developed finite element model. The model consists of a straight five
node twenty-two degrees of freedom beam element. The stability
analysis of the beams is studied by solving the matrix form
characteristic equations of the system. The principle of virtual work
and the first order shear deformation theory are employed to analyze
the beams with large deformation and small strains. The regions of
dynamic instability of the beam are determined by solving the
obtained Mathieu form of differential equations. The effects of nonconservative
loads, shear stiffness, and damping parameters on
stability and response of the beams are examined. Several numerical
calculations are presented to compare the results with data reported
by other researchers.
Abstract: The purpose of this paper is to perform a multidisciplinary design and analysis (MDA) of honeycomb panels used in the satellites structural design. All the analysis is based on clamped-free boundary conditions. In the present work, detailed finite element models for honeycomb panels are developed and analysed. Experimental tests were carried out on a honeycomb specimen of which the goal is to compare the previous modal analysis made by the finite element method as well as the existing equivalent approaches. The obtained results show a good agreement between the finite element analysis, equivalent and tests results; the difference in the first two frequencies is less than 4% and less than 10% for the third frequency. The results of the equivalent model presented in this analysis are obtained with a good accuracy. Moreover, investigations carried out in this research relate to the honeycomb plate modal analysis under several aspects including the structural geometrical variation by studying the various influences of the dimension parameters on the modal frequency, the variation of core and skin material of the honeycomb. The various results obtained in this paper are promising and show that the geometry parameters and the type of material have an effect on the value of the honeycomb plate modal frequency.
Abstract: Solidification cracking and hydrogen cracking are some defects generated in the fusion welding of ultrahigh carbon steels. However, friction stir welding (FSW) of such steels, being a solid-state technique, has been demonstrated to alleviate such problems encountered in traditional welding. FSW include different process parameters that must be carefully defined prior processing. These parameters included but not restricted to: tool feed, tool RPM, tool geometry, tool tilt angle. These parameters form a key factor behind avoiding warm holes and voids behind the tool and in achieving a defect-free weld. More importantly, these parameters directly affect the microstructure of the weld and hence the final mechanical properties of weld. For that, 3D finite element (FE) thermo-mechanical model was developed using DEFORM 3D to simulate FSW of carbon steel. At points of interest in the joint, tracking is done for history of critical state variables such as temperature, stresses, and strain rates. Typical results found include the ability to simulate different weld zones. Simulations predictions were successfully compared to experimental FSW tests. It is believed that such a numerical model can be used to optimize FSW processing parameters to favor desirable defect free weld with better mechanical properties.
Abstract: Stick models are widely used in studying the
behaviour of straight as well as skew bridges and viaducts subjected
to earthquakes while carrying out preliminary studies. The
application of such models to highly curved bridges continues to
pose challenging problems. A viaduct proposed in the foothills of the
Himalayas in Northern India is chosen for the study. It is having 8
simply supported spans @ 30 m c/c. It is doubly curved in horizontal
plane with 20 m radius. It is inclined in vertical plane as well. The
superstructure consists of a box section. Three models have been
used: a conventional stick model, an improved stick model and a 3D
finite element model. The improved stick model is employed by
making use of body constraints in order to study its capabilities. The
first 8 frequencies are about 9.71% away in the latter two models.
Later the difference increases to 80% in 50th mode. The viaduct was
subjected to all three components of the El Centro earthquake of May
1940. The numerical integration was carried out using the Hilber-
Hughes-Taylor method as implemented in SAP2000. Axial forces
and moments in the bridge piers as well as lateral displacements at
the bearing levels are compared for the three models. The maximum
difference in the axial forces and bending moments and
displacements vary by 25% between the improved and finite element
model. Whereas, the maximum difference in the axial forces,
moments, and displacements in various sections vary by 35%
between the improved stick model and equivalent straight stick
model. The difference for torsional moment was as high as 75%. It is
concluded that the stick model with body constraints to model the
bearings and expansion joints is not desirable in very sharp S curved
viaducts even for preliminary analysis. This model can be used only
to determine first 10 frequency and mode shapes but not for member
forces. A 3D finite element analysis must be carried out for
meaningful results.
Abstract: Finite element method was applied to model damage
development in the femoral neck during a sideways fall. The femoral
failure was simulated using the maximum principal strain criterion.
The evolution of damage was consistent with previous studies. It was
initiated by compressive failure at the junction of the superior aspect
of the femoral neck and the greater trochanter. It was followed by
tensile failure that occurred at the inferior aspect of the femoral neck
before a complete transcervical fracture was observed. The estimated
failure line was less than 50° from the horizontal plane (Pauwels type
II).
Abstract: In this paper, a new formulation for acoustics coupled with linear elasticity is presented. The primary objective of the work is to develop a three dimensional hp adaptive finite element method code destinated for modeling of acoustics of human head. The code will have numerous applications e.g. in designing hearing protection devices for individuals working in high noise environments. The presented work is in the preliminary stage. The variational formulation has been implemented and tested on a sequence of meshes with concentric multi-layer spheres, with material data representing the tissue (the brain), skull and the air. Thus, an efficient solver for coupled elasticity/acoustics problems has been developed, and tested on high contrast material data representing the human head.