Abstract: A formulation of postbuckling analysis of end supported rods under self-weight has been presented by the variational method. The variational formulation involving the strain energy due to bending and the potential energy of the self-weight, are expressed in terms of the intrinsic coordinates. The variational formulation is accomplished by introducing the Lagrange multiplier technique to impose the boundary conditions. The finite element method is used to derive a system of nonlinear equations resulting from the stationary of the total potential energy and then Newton-Raphson iterative procedure is applied to solve this system of equations. The numerical results demonstrate the postbluckled configurations of end supported rods under self-weight. This finite element method based on variational formulation expressed in term of intrinsic coordinate is highly recommended for postbuckling analysis of end-supported rods under self-weight.
Abstract: This paper focuses on a variational formulation of large amplitude free vibration behavior of a very sag marine cable. In the static equilibrium state, the marine cable has a very large sag configuration. In the motion state, the marine cable is assumed to vibrate in in-plane motion with large amplitude from the static equilibrium position. The total virtual work-energy of the marine cable at the dynamic state is formulated which involves the virtual strain energy due to axial deformation, the virtual work done by effective weight, and the inertia forces. The equations of motion for the large amplitude free vibration of marine cable are obtained by taking into account the difference between the Euler’s equation in the static state and the displaced state. Based on the Galerkin finite element procedure, the linear and nonlinear stiffness matrices, and mass matrices of the marine cable are obtained and the eigenvalue problem is solved. The natural frequency spectrum and the large amplitude free vibration behavior of marine cable are presented.
Abstract: Investigation of fracture of wood components can prevent from catastrophic failures. Created fracture process zone (FPZ) in crack tip vicinity has important effect on failure of cracked composite materials. In this paper, a failure criterion for fracture investigation of cracked wood specimens under mixed mode I/II loading is presented. This criterion is based on maximum strain energy release rate and material nonlinearity in the vicinity of crack tip due to presence of microcracks. Verification of results with available experimental data proves the coincidence of the proposed criterion with the nature of fracture of wood. To simplify the estimation of nonlinear properties of FPZ, a damage factor is also introduced for engineering and application purposes.
Abstract: This paper is concerned with the development of a
fully implicit and purely Eulerian fluid-structure interaction method
tailored for the modeling of the large deformations of elastic
membranes in a surrounding Newtonian fluid. We consider a
simplified model for the mechanical properties of the membrane, in
which the surface strain energy depends on the membrane stretching.
The fully Eulerian description is based on the advection of a modified
surface tension tensor, and the deformations of the membrane are
tracked using a level set strategy. The resulting nonlinear problem
is solved by a Newton-Raphson method, featuring a quadratic
convergence behavior. A monolithic solver is implemented, and we
report several numerical experiments aimed at model validation and
illustrating the accuracy of the presented method. We show that
stability is maintained for significantly larger time steps.
Abstract: The breakage of bonding between sand particles and their dislodgment from the borehole wall are among the main factors resulting in a borehole failure in poorly cemented granular formations. The grain debonding usually precedes the borehole failure and it can be considered as a sign that the onset of the borehole collapse is imminent. Detecting the bonding breakage point and introducing an appropriate failure criterion will play an important role in borehole stability analysis. To study the influence of different factors on the initiation of sand bonding breakage at the borehole wall, a series of laboratory tests was designed and conducted on poorly cemented sand samples. The total absorbed strain energy per volume of material up to the point of the observed particle debonding was computed. The results indicated that the particle bonding breakage point at the borehole wall was reached both before and after the peak strength of the thick-walled hollow cylinder specimens depending on the stress path and cement content. Three different cement contents and two borehole sizes were investigated to study the influence of the bonding strength and scale on the particle dislodgment. Test results showed that the stress path has a significant influence on the onset of the sand bonding breakage. It was shown that for various stress paths, there is a near linear relationship between the absorbed energy and the normal effective mean stress.
Abstract: By applying coating onto a structural component, the corrosion and/or wear resistance requirements of the surface can be fulfilled. Since the layer adhesion of the coating influences the mechanical integrity of the coat/substrate interface during the service time, it should be examined accurately. At the present work, the tensile bonding strength of the 316 stainless steel plasma sprayed coating on aluminum substrate was determined by using tensile adhesion test, TAT, specimen. The interfacial fracture toughness was specified using four-point bend specimen containing a saw notch and modified chevron-notched short-bar (SB) specimen. The coating microstructure and fractured specimen surface were examined by using scanning electron- and optical-microscopy. The investigation of coated surface after tensile adhesion test indicates that the failure mechanism is mostly cohesive and rarely adhesive type. The calculated value of critical strain energy release rate proposes relatively good interface status. It seems that four-point bending test offers a potentially more sensitive means for evaluation of mechanical integrity of coating/substrate interfaces than is possible with the tensile test. The fracture toughness value reported for the modified chevron-notched short-bar specimen testing cannot be taken as absolute value because its calculation is based on the minimum stress intensity coefficient value which has been suggested for the fracture toughness determination of homogeneous parts in the ASTM E1304-97 standard.
Abstract: In this study, a general approach to the reliability
based limit analysis of laterally loaded piles is presented. In
engineering practice the uncertainties play a very important role. The
aim of this study is to evaluate the lateral load capacity of free-head
and fixed-head long pile when plastic limit analysis is considered. In
addition to the plastic limit analysis to control the plastic behaviour
of the structure, uncertain bound on the complementary strain energy
of the residual forces is also applied. This bound has significant effect
for the load parameter. The solution to reliability-based problems is
obtained by a computer program which is governed by the reliability
index calculation.
Abstract: This paper deals with nonlinear vibration analysis
using finite element method for frame structures consisting of elastic
and viscoelastic damping layers supported by multiple nonlinear
concentrated springs with hysteresis damping. The frame is supported
by four nonlinear concentrated springs near the four corners. The
restoring forces of the springs have cubic non-linearity and linear
component of the nonlinear springs has complex quantity to represent
linear hysteresis damping. The damping layer of the frame structures
has complex modulus of elasticity. Further, the discretized equations in
physical coordinate are transformed into the nonlinear ordinary
coupled differential equations using normal coordinate corresponding
to linear natural modes. Comparing shares of strain energy of the
elastic frame, the damping layer and the springs, we evaluate the
influences of the damping couplings on the linear and nonlinear impact
responses. We also investigate influences of damping changed by
stiffness of the elastic frame on the nonlinear coupling in the damped
impact responses.
Abstract: In this paper, an analytical simplified method for
calculating elasto-plastic stresses strains of notched bodies subject to
non-proportional loading paths is discussed. The method was based
on the Neuber notch correction, which relates the incremental elastic
and elastic-plastic strain energy densities at the notch root and the
material constitutive relationship. The validity of the method was
presented by comparing computed results of the proposed model
against finite element numerical data of notched shaft. The
comparison showed that the model estimated notch-root elasto-plastic
stresses strains with good accuracy using linear-elastic stresses. The
prosed model provides more efficient and simple analysis method
preferable to expensive experimental component tests and more
complex and time consuming incremental non-linear FE analysis.
The model is particularly suitable to perform fatigue life and fatigue
damage estimates of notched components subjected to nonproportional
loading paths.
Abstract: Each year nearly nine hundred persons die in head injuries and over fifty thousand persons are severely injured due to non wearing of helmets. In motor cycle accidents, the human head is exposed to heavy impact loading against natural protection. In this work, an attempt has been made for analyzing the helmet with all the standard data. The simulation software ‘ANSYS’ is used to analyze the helmet with different conditions such as bottom fixed-load on top surface, bottom fixed -load on top line, side fixed –load on opposite surface, side fixed-load on opposite line and dynamic analysis. The maximum force of 19.5 kN is applied on the helmet to study the model in static and dynamic conditions. The simulation has been carried out for the static condition for the parameters like total deformation, strain energy, von-Mises stress for different cases. The dynamic analysis has been performed for the parameter like total deformation and equivalent elastic strain. The result shows that these values are concentrated in the retention portion of the helmet. These results have been compared with the standard experimental data proposed by the BIS and well within the acceptable limit.
Abstract: Shear testing is one of the most complex testing areas where available methods and specimen geometries are different from each other. Therefore, a modified shear test specimen (MSTS) combining the simple uniaxial test with a zone of interest (ZOI) is tested which gives almost the pure shear. In this study, material parameters of polypropylene (PP) and high density polyethylene (HDPE) are first measured by tensile tests with a dogbone shaped specimen. These parameters are then used as an input for the finite element analysis. Secondly, a specially designed specimen (MSTS) is used to perform the shear stress tests in a tensile testing machine to get the results in terms of forces and extension, crack initiation etc. Scanning Electron Microscopy (SEM) is also performed on the shear fracture surface to find material behavior. These experiments are then simulated by finite element method and compared with the experimental results in order to confirm the simulation model. Shear stress state is inspected to find the usability of the proposed shear specimen. Finally, a geometry correction factor can be established for these two materials in this specific loading and geometry with notch using Linear Elastic Fracture Mechanics (LEFM). By these results, strain energy of shear failure and stress intensity factor (SIF) of shear of these two polymers are discussed in the special application of the screw cap opening of the medical or food packages with a temper evidence safety solution.
Abstract: The stiffness of the workpiece is very important to reduce the errors in manufacturing process. The high stiffness of the workpiece can be achieved by optimal positioning of fixture elements in the fixture. The minimization of the sum of the nodal deflection normal to the surface is used as objective function in previous research. The deflection in other direction has been neglected. The 3-2-1 fixturing principle is not valid for metal sheets due to its flexible nature. We propose a new fixture layout optimization method N-3-2-1 for metal sheets that uses the strain energy of the finite elements. This method combines the genetic algorithm and finite element analysis. The objective function in this method is to minimize the sum of all the element strain energy. By using the concept of element strain energy, the deformations in all the directions have been considered. Strain energy and stiffness are inversely proportional to each other. So, lower the value of strain energy, higher will be the stiffness. Two different kinds of case studies are presented. The case studies are solved for both objective functions; element strain energy and nodal deflection. The result are compared to verify the propose method.
Abstract: The separation efficiency of a hydrocyclone has
extensively been considered on the rigid particle assumption. A
collection of experimental studies have demonstrated their
discrepancies from the modeling and simulation results. These
discrepancies caused by the actual particle elasticity have generally
led to a larger amount of energy consumption in the separation
process. In this paper, the influence of particle elasticity on the
separation efficiency of a hydrocyclone system was investigated
through the Finite Element (FE) simulations using crude oil droplets
as the elastic particles. A Reitema-s design hydrocyclone with a
diameter of 8 mm was employed to investigate the separation
mechanism of the crude oil droplets from water. The cut-size
diameter eter of the crude oil was 10 - Ðçm in order to fit with the
operating range of the adopted hydrocylone model. Typical
parameters influencing the performance of hydrocyclone were varied
with the feed pressure in the range of 0.3 - 0.6 MPa and feed
concentration between 0.05 – 0.1 w%. In the simulation, the Finite
Element scheme was applied to investigate the particle-flow
interaction occurred in the crude oil system during the process. The
interaction of a single oil droplet at the size of 10 - Ðçm to the flow
field was observed. The feed concentration fell in the dilute flow
regime so the particle-particle interaction was ignored in the study.
The results exhibited the higher power requirement for the separation
of the elastic particulate system when compared with the rigid
particulate system.
Abstract: When considering the development of constitutive
equations describing the behavior of materials under cyclic plastic
strains, different kinds of formulations can be adopted. The primary
intention of this study is to develop computer programming of
plasticity models to accurately predict the life of engineering
components. For this purpose, the energy or cyclic strain is computed
in multi-surface plasticity models in non-proportional loading and to
present their procedures and codes results.
Abstract: Displacement measurement was conducted on compact normal and shear specimens made of acrylic homogeneous material subjected to mixed-mode loading by digital image correlation. The intelligent hybrid method proposed by Nishioka et al. was applied to the stress-strain analysis near the crack tip. The accuracy of stress-intensity factor at the free surface was discussed from the viewpoint of both the experiment and 3-D finite element analysis. The surface images before and after deformation were taken by a CMOS camera, and we developed the system which enabled the real time stress analysis based on digital image correlation and inverse problem analysis. The great portion of processing time of this system was spent on displacement analysis. Then, we tried improvement in speed of this portion. In the case of cracked body, it is also possible to evaluate fracture mechanics parameters such as the J integral, the strain energy release rate, and the stress-intensity factor of mixed-mode. The 9-points elliptic paraboloid approximation could not analyze the displacement of submicron order with high accuracy. The analysis accuracy of displacement was improved considerably by introducing the Newton-Raphson method in consideration of deformation of a subset. The stress-intensity factor was evaluated with high accuracy of less than 1% of the error.
Abstract: This paper presents a novel method for prediction of
the mechanical behavior of proximal femur using the general
framework of the quantitative computed tomography (QCT)-based
finite element Analysis (FEA). A systematic imaging and modeling
procedure was developed for reliable correspondence between the
QCT-based FEA and the in-vitro mechanical testing. A speciallydesigned
holding frame was used to define and maintain a unique
geometrical reference system during the analysis and testing. The
QCT images were directly converted into voxel-based 3D finite
element models for linear and nonlinear analyses. The equivalent
plastic strain and the strain energy density measures were used to
identify the critical elements and predict the failure patterns. The
samples were destructively tested using a specially-designed gripping
fixture (with five degrees of freedom) mounted within a universal
mechanical testing machine. Very good agreements were found
between the experimental and the predicted failure patterns and the
associated load levels.
Abstract: An inflation–extension test with human vena cava
inferior was performed with the aim to fit a material model. The vein
was modeled as a thick–walled tube loaded by internal pressure and
axial force. The material was assumed to be an incompressible
hyperelastic fiber reinforced continuum. Fibers are supposed to be
arranged in two families of anti–symmetric helices. Considered
anisotropy corresponds to local orthotropy. Used strain energy
density function was based on a concept of limiting strain
extensibility. The pressurization was comprised by four pre–cycles
under physiological venous loading (0 – 4kPa) and four cycles under
nonphysiological loading (0 – 21kPa). Each overloading cycle was
performed with different value of axial weight. Overloading data
were used in regression analysis to fit material model. Considered
model did not fit experimental data so good. Especially predictions
of axial force failed. It was hypothesized that due to
nonphysiological values of loading pressure and different values of
axial weight the material was not preconditioned enough and some
damage occurred inside the wall. A limiting fiber extensibility
parameter Jm was assumed to be in relation to supposed damage.
Each of overloading cycles was fitted separately with different values
of Jm. Other parameters were held the same. This approach turned out
to be successful. Variable value of Jm can describe changes in the
axial force – axial stretch response and satisfy pressure – radius
dependence simultaneously.
Abstract: Rise/span ratio has been mentioned as one of the
reasons which contribute to the lower buckling load as compared to
the Classical theory buckling load but this ratio has not been quantified
in the equation. The purpose of this study was to determine a more
realistic buckling load by quantifying the effect of the rise/span ratio
because experiments have shown that the Classical theory
overestimates the load. The buckling load equation was derived based
on the theorem of work done and strain energy. Thereafter, finite
element modeling and simulation using ABAQUS was done to
determine the variables that determine the constant in the derived
equation. The rise/span was found to be the determining factor of the
constant in the buckling load equation. The derived buckling load
correlates closely to the load obtained from experiments.
Abstract: Dual phase steels (DPS)s have a microstructure
consisting of a hard second phase called Martensite in the soft Ferrite
matrix. In recent years, there has been interest in dual-phase steels,
because the application of these materials has made significant usage;
particularly in the automotive sector Composite microstructure of
(DPS)s exhibit interesting characteristic mechanical properties such
as continuous yielding, low yield stress to tensile strength
ratios(YS/UTS), and relatively high formability; which offer
advantages compared with conventional high strength low alloy
steels(HSLAS). The research dealt with the characterization of
damage in (DPS)s. In this study by review the mechanisms of failure
due to volume fraction of martensite second phase; a new method is
introduced to identifying the mechanisms of failure in the various
phases of these types of steels. In this method the acoustic emission
(AE) technique was used to detect damage progression. These failure
mechanisms consist of Ferrite-Martensite interface decohesion and/or
martensite phase fracture. For this aim, dual phase steels with
different volume fraction of martensite second phase has provided by
various heat treatment methods on a low carbon steel (0.1% C), and
then AE monitoring is used during tensile test of these DPSs. From
AE measurements and an energy ratio curve elaborated from the
value of AE energy (it was obtained as the ratio between the strain
energy to the acoustic energy), that allows detecting important
events, corresponding to the sudden drops. These AE signals events
associated with various failure mechanisms are classified for ferrite
and (DPS)s with various amount of Vm and different martensite
morphology. It is found that AE energy increase with increasing Vm.
This increasing of AE energy is because of more contribution of
martensite fracture in the failure of samples with higher Vm. Final
results show a good relationship between the AE signals and the
mechanisms of failure.
Abstract: In this paper delamination phenomenon in
Carbon-Epoxy laminated composite material is investigated
numerically. Arcan apparatus and specimen is modeled in ABAQUS
finite element software for different loading conditions and crack
geometries. The influence of variation of crack geometry on
interlaminar fracture stress intensity factor and energy release rate for
various mixed mode ratios and pure mode I and II was studied. Also,
correction factors for this specimen for different crack length ratios
were calculated. The finite element results indicate that for loading
angles close to pure mode-II loading, a high ratio of mode-II to
mode-I fracture is dominant and there is an opposite trend for loading
angles close to pure mode-I loading. It confirms that by varying the
loading angle of Arcan specimen pure mode-I, pure mode-II and a
wide range of mixed-mode loading conditions can be created and
tested. Also, numerical results confirm that the increase of the mode-
II loading contribution leads to an increase of fracture resistance in
the CF/PEI composite (i.e., a reduction in the total strain energy
release rate) and the increase of the crack length leads to a reduction
of interlaminar fracture resistance in the CF/PEI composite (i.e., an
increase in the total interlaminar strain energy release rate).