Search results for: finite element time stepping method

Authors: Clarinda Vitorino Nhangumbe, Fredericks Ebrahim, Betuel Canhanga

Abstract:

The price of an option weather derivatives can be approximated as a solution of the two-dimensional convection-diffusion dominant partial differential equation derived from the Ornstein-Uhlenbeck process, where one variable represents the weather dynamics and the other variable represent the underlying weather index. With appropriate financial boundary conditions, the solution of the pricing equation is approximated using a nonstandard finite difference method. It is shown that the proposed numerical scheme preserves positivity as well as stability and consistency. In order to illustrate the accuracy of the method, the numerical results are compared with other methods. The model is tested for real weather data.

Keywords: nonstandard finite differences, Ornstein-Uhlenbeck process, partial differential equations approach, weather derivatives

Procedia PDF Downloads 111

Authors: A. S. Ibrahim Mamane, S. Giljean, M.-J. Pac, G. L’Hostis

Abstract:

Fiber reinforced composite laminates are particularly subject to residual stresses due to their heterogeneity and the complex chemical, mechanical and thermal mechanisms that occur during their processing. Residual stresses are now well known to cause damage accumulation, shape instability, and behavior disturbance in composite parts. Many works exist in the literature on techniques for minimizing residual stresses in thermosetting and thermoplastic composites mainly. To study in-depth the influence of processing mechanisms on the formation of residual stresses and to minimize them by establishing a reliable correlation, it is essential to be able to measure very precisely the profile of residual stresses in the composite. Residual stresses are important data to consider when sizing composite parts and predicting their behavior. The incremental hole drilling is very effective in measuring the gradient of residual stresses in composite laminates. This method is semi-destructive and consists of drilling incrementally a hole through the thickness of the material and measuring relaxation strains around the hole for each increment using three strain gauges. These strains are then converted into residual stresses using a matrix of coefficients. These coefficients, called calibration coefficients, depending on the diameter of the hole and the dimensions of the gauges used. The reliability of the incremental hole drilling depends on the accuracy with which the calibration coefficients are determined. These coefficients are calculated using a finite element model. The samples’ features and the experimental conditions must be considered in the simulation. Any mismatch can lead to inadequate calibration coefficients, thus introducing errors on residual stresses. Several calibration coefficient correction methods exist for isotropic material, but there is a lack of information on this subject concerning composite laminates. In this work, a Python program was developed to automatically generate the adequate finite element model. This model allowed us to perform a parametric study to assess the influence of experimental errors on the calibration coefficients. The results highlighted the sensitivity of the calibration coefficients to the considered errors and gave an order of magnitude of the precisions required on the experimental device to have reliable measurements. On the basis of these results, improvements were proposed on the experimental device. Furthermore, a numerical method was proposed to correct the calibration coefficients for different types of materials, including thick composite parts for which the analytical approach is too complex. This method consists of taking into account the experimental errors in the simulation. Accurate measurement of the experimental errors (such as eccentricity of the hole, angular deviation of the gauges from their theoretical position, or errors on increment depth) is therefore necessary. The aim is to determine more precisely the residual stresses and to expand the validity domain of the incremental hole drilling technique.

Keywords: fiber reinforced composites, finite element simulation, incremental hole drilling method, numerical correction of the calibration coefficients, residual stresses

Procedia PDF Downloads 132

Authors: Talebi Aliasghar, Ebrahimpour Komeleh Hooman, Maghsoudi Ali Akbar

Abstract:

In reinforced concrete (RC) structures, beam-column connection region has a considerable effect on the behavior of structures. Using fiber reinforcement polymer (FRP) for the strengthening of connections in RC structures can be one of the solutions to retrofitting this zone which result in the enhanced behavior of structure. In this paper, these changes in behavior by using FRP for high strength concrete beam-column connection have been studied by finite element modeling. The concrete damage plasticity (CDP) model has been used to analyze the RC. The results illustrated a considerable development in load-bearing capacity but also a noticeable reduction in ductility. The study also assesses these qualities for several modes of strengthening and suggests the most effective mode of strengthening. Using FRP in flexural zone and FRP with 45-degree oriented fibers in shear zone of joint showed the most significant change in behavior.

Keywords: HSC, beam-column connection, Fiber Reinforcement Polymer, FRP, Finite Element Modeling, FEM

Procedia PDF Downloads 159

Authors: Naeim Akbari Shahkhosravi, Reza Kakavand, Helen M. S. Davies, Amin Komeili

Abstract:

Equine locomotion and performance are significantly affected by hoof health. One of the most critical diseases of the hoof is laminitis, which can lead to horse lameness in a severe condition. This disease exhibits the mechanical properties degradation of the laminar junction tissue within the hoof. Therefore, it is essential to investigate the biomechanics of the hoof, focusing specifically on excessive and cumulatively accumulated stresses within the laminar junction tissue. For this aim, the current study generated a novel equine hoof Finite Element (FE) model under dynamic physiological loading conditions and employing a hyperelastic material model. Associated tissues of the equine hoof were segmented from computed tomography scans of an equine forelimb, including the navicular bone, third phalanx, sole, frog, laminar junction, digital cushion, and medial- dorsal- lateral wall areas. The inner tissues were connected based on the hoof anatomy, and the hoof was under a dynamic loading over cyclic strides at the trot. The strain distribution on the hoof wall of the model was compared with the published in vivo strain measurements to validate the model. Then the validated model was used to study the development of laminitis. The ultimate stress tolerated by the laminar junction before rupture was considered as a stress threshold. The tissue damage was simulated through iterative reduction of the tissue’s mechanical properties in the presence of excessive maximum principal stresses. The findings of this investigation revealed how damage initiates from the medial and lateral sides of the tissue and propagates through the hoof dorsal area.

Keywords: horse hoof, laminitis, finite element model, continuous damage

Procedia PDF Downloads 182

Authors: Adriano Trono, Federico Pinto, Diego Turello, Marcelo A. Ceballos

Abstract:

Numerical modeling of dynamic soil-structure interaction problems requires an adequate representation of the unbounded characteristics of the ground, material non-linearity of soils, and geometrical non-linearities such as large displacements due to rocking of the structure. In order to account for these effects simultaneously, it is often required that the equations of motion are solved in the time domain. However, boundary conditions in conventional finite element codes generally present shortcomings in fully absorbing the energy of outgoing waves. In this sense, the Perfectly Matched Layers (PML) technique allows a satisfactory absorption of inclined body waves, as well as surface waves. However, the PML domain is inherently unstable, meaning that it its instability does not depend upon the discretization considered. One way to stabilize the PML domain is to use multiaxial stretching functions. This development is questionable because some Jacobian terms of the coordinate transformation are not accounted for. For this reason, the resulting absorbing layer element is often referred to as "uncorrected M-PML” in the literature. In this work, the strong formulation of the "corrected M-PML” absorbing layer is proposed using multiaxial stretching functions that incorporate all terms of the coordinate transformation. The results of the stable model are compared with reference solutions obtained from extended domain models.

Keywords: mixed finite elements, multiaxial stretching functions, perfectly matched layer, soil-structure interaction

Procedia PDF Downloads 70

Authors: Emre Ozaslan, Bulent Acar, Mehmet Ali Guler

Abstract:

Unidirectional composites are very popular structural materials used in aerospace, marine, energy and automotive industries thanks to their superior material properties. However, the mechanical behavior of composite materials is more complicated than isotropic materials because of their anisotropic nature. Also, a stress concentration availability on the structure, like a hole, makes the problem further complicated. Therefore, enormous number of tests require to understand the mechanical behavior and strength of composites which contain stress concentration. Accurate finite element analysis and analytical models enable to understand mechanical behavior and predict the strength of composites without enormous number of tests which cost serious time and money. In this study, unidirectional Carbon/Epoxy composite specimens with central circular hole were investigated in terms of stress concentration factor and strength prediction. The composite specimens which had different specimen wide (W) to hole diameter (D) ratio were tested to investigate the effect of hole size on the stress concentration and strength. Also, specimens which had same specimen wide to hole diameter ratio, but varied sizes were tested to investigate the size effect. Finite element analysis was performed to determine stress concentration factor for all specimen configurations. For quasi-isotropic laminate, it was found that the stress concentration factor increased approximately %15 with decreasing of W/D ratio from 6 to 3. Point stress criteria (PSC), inherent flaw method and progressive failure analysis were compared in terms of predicting the strength of specimens. All methods could predict the strength of specimens with maximum %8 error. PSC was better than other methods for high values of W/D ratio, however, inherent flaw method was successful for low values of W/D. Also, it is seen that increasing by 4 times of the W/D ratio rises the failure strength of composite specimen as %62.4. For constant W/D ratio specimens, all the strength prediction methods were more successful for smaller size specimens than larger ones. Increasing the specimen width and hole diameter together by 2 times reduces the specimen failure strength as %13.2.

Keywords: failure, strength, stress concentration, unidirectional composites

Procedia PDF Downloads 156

Authors: Hasan Husnu Korkmaz, Serra Zerrin Korkmaz

Abstract:

The reinforced concrete (RC) framed structures composed of beams, columns, shear walls and the slabs. The other members are assumed to be nonstructural. Especially the brick infill walls which are used to separate the rooms or spaces are just handled as dead loads. On the other hand, if these infills are constructed within the frame bays, they also have higher shear and compression capacities. It is a well-known fact that, brick infills increase the lateral rigidity of the structure and thought to be a reserve capacity in the design. But, brick infills can create unfavorable failure or damage modes in the earthquake action such as soft story or short columns. The increase in the lateral rigidity also causes an over estimation of natural period of the structure and the corresponding earthquake loads in the design are less than the actual ones. In order to obtain accurate and realistic design results, the infills must be modelled in the structural design and their capacities must be included. Unfortunately, in Turkish Earthquake Code, there is no design methodology for the engineers. In this paper, finite element modelling of infilled reinforced concrete structures are studied. The proposed or used method is compared with the experimental results of a previous study. The effect of infills on the structural response is expressed within the paper.

Keywords: seismic loading, brick infills, finite element analysis, reinforced concrete, earthquake code

Procedia PDF Downloads 314

Authors: Nisha Nandlal Sharma, Julaluk Carmai, Saiprasit Koetniyom, Bernd Markert

Abstract:

Whiplash Injuries are usually associated with car accidents. The sudden forward or backward jerk to head causes neck strain, which is the result of damage to the muscle or tendons. Neck pain and headaches are the two most common symptoms of whiplash. Symptoms of whiplash are commonly reported in studies but the Injury mechanism is poorly understood. Neck muscles are the most important factor to study the neck Injury. This study focuses on the development of finite element (FE) model of human neck muscle to study the whiplash injury mechanism and effect of active muscle contraction on occupant kinematics. A detailed study of Injury mechanism will promote development and evaluation of new safety systems in cars, hence reducing the occurrence of severe injuries to the occupant. In present study, an active human finite element (FE) model with 3D neck muscle model is developed. Neck muscle was modeled with a combination of solid tetrahedral elements and 1D beam elements. Muscle active properties were represented by beam elements whereas, passive properties by solid tetrahedral elements. To generate muscular force according to inputted activation levels, Hill-type muscle model was applied to beam elements. To simulate non-linear passive properties of muscle, solid elements were modeled with rubber/foam material model. Material properties were assigned from published experimental tests. Some important muscles were then inserted into THUMS (Total Human Model for Safety) 50th percentile male pedestrian model. To reduce the simulation time required, THUMS lower body parts were not included. Posterior to muscle insertion, THUMS was given a boundary conditions similar to experimental tests. The model was exposed to 4g and 7g rear impacts as these load impacts are close to low speed impacts causing whiplash. The effect of muscle activation level on occupant kinematics during whiplash was analyzed.

Keywords: finite element model, muscle activation, neck muscle, whiplash injury prevention

Procedia PDF Downloads 357

Authors: Swati Sharma, R. P. Sharma

Abstract:

We intend to study the nonlinear evolution of the parallel propagating finite frequency Alfvén wave (also called Dispersive Alfvén wave/Hall MHD wave) propagating in the solar wind regime of the solar region when a perpendicularly propagating magnetosonic wave is present in the background. The finite frequency Alfvén wave behaves differently from the usual non-dispersive behavior of the Alfvén wave. To study the nonlinear processes (such as filamentation) taking place in the solar regions such as solar wind, the dynamical equation of both the waves are derived. Numerical simulation involving finite difference method for the time domain and pseudo spectral method for the spatial domain is then performed to analyze the transient evolution of these waves. The power spectra of the Dispersive Alfvén wave is also investigated. The power spectra shows the distribution of the magnetic field intensity of the Dispersive Alfvén wave over different wave numbers. For DAW the spectra shows a steepening for scales larger than the proton inertial length. This means that the wave energy gets transferred to the solar wind particles as the wave reaches higher wave numbers. This steepening of the power spectra can be explained on account of the finite frequency of the Alfvén wave. The obtained results are consistent with the observations made by CLUSTER spacecraft.

Keywords: solar wind, turbulence, dispersive alfven wave

Procedia PDF Downloads 601

Authors: Subodh Nath Patel, Raja Raman, Mananshi Adhikary, Jitendra Mathur, Sandip Bhattacharyya

Abstract:

The composition adjustment by sealed argon bubbling–oxygen blowing (CAS-OB) process is a process designed for adjusting steel composition and temperature during secondary metallurgy. One of the equipment in the said process is a snorkel or bell, fixed to a movable bracket. Snorkel serves the purpose of feeding ferroalloys into the liquid metal simultaneously removing gases to the gas cleaning system through its port at its top. The bell-shaped snorkel consists of two parts. The upper part has an inside liner, and the lower part is lined on both side with high-alumina castable reinforced with 2% stainless steel needles. Both the parts are coupled with a flange bolt system. These flanges were found to get buckled during operation, and the gap was generating between them. This problem was chronic since its. It was expected to give a life of 80 heats, but it was failing within 45-50 heats. After every 25-30 heats, it had to be repaired by changing and/or tightening its nuts and bolts. Visual observation, microstructural analysis through optical microscopes and SEM, hardness measurement and thermal strain calculation were carried out to find out the root cause of this problem. The calculated thermal strain was compared with actual thermal strain; comparison of the two revealed that thermal strain was responsible for buckling. Finite element analysis (FEA) was carried out to reaffirm the effect temperature on the flanges. FEA was also used in the modification in the design of snorkel flange to accommodate thermal strain. Thermal insulation was also recommended which increased its life from 45 heats to 65 heats, impacting business process positively.

Keywords: CAS OB process, finite element analysis, snorkel, thermal strain

Procedia PDF Downloads 137

Authors: Riya Roy, Anjana Bhasi

Abstract:

Column-supported embankments provide a practical and efficient solution for construction on soft soil due to the low cost and short construction times. In the recent years, geosynthetic have been used in combination with column systems to support embankments. The load transfer mechanism in these systems is a combination of soil arching effect, which occurs between columns and membrane effect of the geosynthetic. This paper aims at the study of soil arching effect on columnar embankments using finite element software, ABAQUS. An axisymmetric finite element model is generated and using this model, parametric studies are carried out. Thus the effects of various factors such as height of embankment fill, elastic modulus of pile and tensile stiffness of geosynthetic, on soil arching have been studied. The development of negative skin friction along the pile-soil interface have also been studied and the results obtained from this study are compared with the current design methods.

Keywords: ABAQUS, geosynthetic, negative skin friction, soil arching

Procedia PDF Downloads 379

Authors: Tzu-Hsuan Lei, Shu-Chien Wu, Kuang-Che Lo, Shu-Chi Liu, Yu-Cheng Liu

Abstract:

The focus of this study was on the effects of air propagation associated with vented holes on acoustic resistance properties. A cylindrical hole with diameter and depth of 0.7 mm and 1.0 mm, respectively, was the research target. By constructing a finite element analytical model of its sound field properties, the acoustic-specific airflow resistance relationships were obtained for the differences in sound pressure and flow velocity at the two ends of this vented hole. In addition, the acoustic properties of this vented hole were included in the in-ear earphone simulation model to complete the sound pressure curve simulation analysis of the in-ear earphone system with a vented hole of corresponding size. Then, the simulation results were compared with actual measurements obtained from the standard system. Based on the results, when the in-ear earphone vented hole simulation model considered the simulated specific airflow resistance values of this cylindrical hole, the overall simulated sound pressure performance was highly consistent with that of measured values. The difference in the first peak values of sound pressure at mid-to-low frequencies was reduced from 5.64% when the simulation model did not consider the specific airflow resistance of the cylindrical hole to 1.18%, and the accuracy of the overall simulation was around 70%. This indicates the importance of the acoustic resistance properties of vented holes. Moreover, as specific airflow resistance values were able to be further quantified, the accuracy of the entire in-ear earphone simulation was ultimately and effectively elevated.

Keywords: specific airflow resistance, vented holes, in-ear earphone, finite element method

Procedia PDF Downloads 43

Authors: Thamer Kubat, Riadh Al Mahiadi, Ahmad Shayan

Abstract:

The alkali-aggregate reaction (AAR) in concrete has a negative influence on the mechanical properties and durability of concrete. Confinement by carbon fiber reinforced polymer (CFRP) is an effective method of treatment for some AAR-affected elements. Eighteen reinforced columns affected by different levels of expansion due to AAR were confined using CFRP to evaluate the effect of expansion level on confinement efficiency. Strength and strain capacities (axial and circumferential) were measured using photogrammetry under uniaxial compressive loading to evaluate the efficiency of CFRP wrapping for the rehabilitation of affected columns. In relation to uniaxial compression capacity, the results indicated that the confinement of AAR-affected columns by one layer of CFRP is sufficient to reach and exceed the load capacity of unaffected sound columns. Parallel to the experimental study, finite element (FE) modeling using ATENA software was employed to predict the behavior of CFRP-confined damaged concrete and determine the possibility of using the model in a parametric study by simulating the number of CFRP layers. A comparison of the experimental results with the results of the theoretical models showed that FE modeling could be used for the prediction of the behavior of confined AAR-damaged concrete.

Keywords: ATENA, carbon fiber reinforced polymer (CFRP), confinement efficiency, finite element (FE)

Procedia PDF Downloads 76

Authors: Chollawat Pookpienlert, Jintana Sanwong

Abstract:

An element a of a semigroup S is called left (right) regular if there exists x in S such that a=xa² (a=a²x) and said to be intra-regular if there exist u,v in such that a=ua²v. Let T(X) be the semigroup of all full transformations on a set X under the composition of maps. For a fixed nonempty subset Y of X, let Fix(X,Y)={α ™ T(X) : yα=y for all y ™ Y}, where yα is the image of y under α. Then Fix(X,Y) is a semigroup of full transformations on X which fix all elements in Y. Here, we characterize left regular, right regular and intra-regular elements of Fix(X,Y) which characterizations are shown as follows: For α ™ Fix(X,Y), (i) α is left regular if and only if Xα\Y = Xα²\Y, (ii) α is right regular if and only if πα = πα², (iii) α is intra-regular if and only if | Xα\Y | = | Xα²\Y | such that Xα = {xα : x ™ X} and πα = {xα⁻¹ : x ™ Xα} in which xα⁻¹ = {a ™ X : aα=x}. Moreover, those regularities are equivalent if Xα\Y is a finite set. In addition, we count the number of those elements of Fix(X,Y) when X is a finite set. Finally, we determine the maximal congruence ρ on Fix(X,Y) when X is finite and Y is a nonempty proper subset of X. If we let | X \Y | = n, then we obtain that ρ = (Fixn x Fixn) ∪ (H ε x H ε) where Fixn = {α ™ Fix(X,Y) : | Xα\Y | < n} and H ε is the group of units of Fix(X,Y). Furthermore, we show that the maximal congruence is unique.

Keywords: intra-regular, left regular, maximal congruence, right regular, transformation semigroup

Procedia PDF Downloads 229

Authors: Kent Salomonsson, Xuefang Zhao, Sara Kallin

Abstract:

Human soft tissue is loaded and deformed by any activity, an effect known as a stress-strain relationship, and is often described by a load and tissue elongation curve. Several advances have been made in the fields of biology and mechanics of soft human tissue. However, there is limited information available on in vivo tissue mechanical characteristics and behavior. Confident mechanical properties of human soft tissue cannot be extrapolated from e.g. animal testing. Thus, there is need for non invasive methods to analyze mechanical characteristics of soft human tissue. In the present study, the internal mechanical conditions of the lower limb, which is subject to an external load, is studied by use of the finite element method. A detailed finite element model of the lower limb is made possible by use of MRI scans. Skin, fat, bones, fascia and muscles are represented separately and the material properties for them are obtained from literature. Previous studies have been shown to address macroscopic deformation features, e.g. indentation depth, to a large extent. However, the detail in which the internal anatomical features have been modeled does not reveal the critical internal strains that may induce hypoxia and/or eventual tissue damage. The results of the present study reveals that lumped material models, i.e. averaging of the material properties for the different constituents, does not capture regions of critical strains in contrast to more detailed models.

Keywords: FEM, tissue, indentation, properties

Procedia PDF Downloads 358

Authors: Alexander P. Smirnov, Sergey A. Matveev, Dmitry A. Zheltkov, Eugene E. Tyrtyshnikov

Abstract:

We propose a new computational technique for multidimensional (multicomponent) Smoluchowski coagulation equation. Using low-rank approximations in Tensor Train format of both the solution and the coagulation kernel, we accelerate the classical finite-difference Runge-Kutta scheme keeping its level of accuracy. The complexity of the taken finite-difference scheme is reduced from O(N^2d) to O(d^2 N log N ), where N is the number of grid nodes and d is a dimensionality of the problem. The efficiency and the accuracy of the new method are demonstrated on concrete problem with known analytical solution.

Keywords: tensor train decomposition, multicomponent Smoluchowski equation, runge-kutta scheme, convolution

Procedia PDF Downloads 432

Authors: Win-Jin Chang, Haw-Long Lee, Yu-Ching Yang

Abstract:

Atomic finite element simulation is applied to investigate the vibration frequency of a single-layer gamma graphyne with an attached mass for the CCCC, SSSS, CFCF, SFSF boundary conditions using the commercial code ANSYS. The fundamental frequencies of the graphyne sheet are compared with the results of the previous study. The results of the comparison are very good in all considered cases. The attached mass causes a shift in the resonant frequency of the graphyne. The frequencies of the single-layer gamma graphyne with an attached mass for different boundary conditions are obtained, and the order based on the boundary condition is CCCC >SSSS > CFCF> SFSF. The highest frequency shift is obtained when the attached mass is located at the center of the graphyne sheet. This is useful for the design of a highly sensitive graphyne-based mass sensor.

Keywords: graphyne, finite element analysis, vibration analysis, frequency shift

Procedia PDF Downloads 212

Authors: Atousa Ataieyan, Salvador A. Gomez-Lopera, Gennaro Sepede

Abstract:

Today, due to the growing urban population and consequently, the increasing water demand in cities, the amount of contaminants entering the water resources is increasing. This can impose harmful effects on the quality of the downstream water. Therefore, predicting the concentration of discharged pollutants at different times and distances of the interested area is of high importance in order to carry out preventative and controlling measures, as well as to avoid consuming the contaminated water. In this paper, the concentration distribution of an injected conservative pollutant in a square reservoir containing four symmetric blocks and three sources using Finite Volume Method (FVM) is simulated. For this purpose, after estimating the flow velocity, classical Advection-Diffusion Equation (ADE) has been discretized over the studying domain by Backward Time- Backward Space (BTBS) scheme. Then, the discretized equations for each node have been derived according to the initial condition, boundary conditions and point contaminant sources. Finally, taking into account the appropriate time step and space step, a computational code was set up in MATLAB. Contaminant concentration was then obtained at different times and distances. Simulation results show how using BTBS differentiating scheme and FVM as a numerical method for solving the partial differential equation of transport is an appropriate approach in the case of two-dimensional contaminant transfer in an advective-diffusive flow.

Keywords: BTBS differentiating scheme, contaminant concentration, finite volume, mass transfer, water pollution

Procedia PDF Downloads 135

Authors: M. A. Askari, M. A. Nazari, P. Perrier, Y. Payan

Abstract:

Growth and remodeling of biological structures have gained lots of attention over the past decades. Determining the response of the living tissues to the mechanical loads is necessary for a wide range of developing fields such as, designing of prosthetics and optimized surgery operations. It is a well-known fact that biological structures are never stress-free, even when externally unloaded. The exact origin of these residual stresses is not clear, but theoretically growth and remodeling is one of the main sources. Extracting body organs from medical imaging, does not produce any information regarding the existing residual stresses in that organ. The simplest cause of such stresses is the gravity since an organ grows under its influence from its birth. Ignoring such residual stresses might cause erroneous results in numerical simulations. Accounting for residual stresses due to tissue growth can improve the accuracy of mechanical analysis results. In this paper, we have implemented a computational framework based on fixed-point iteration to determine the residual stresses due to growth. Using nonlinear continuum mechanics and the concept of fictitious configuration we find the unknown stress-free reference configuration which is necessary for mechanical analysis. To illustrate the method, we apply it to a finite element model of healthy human face whose geometry has been extracted from medical images. We have computed the distribution of residual stress in facial tissues, which can overcome the effect of gravity and cause that tissues remain firm. Tissue wrinkles caused by aging could be a consequence of decreasing residual stress and not counteracting the gravity. Considering these stresses has important application in maxillofacial surgery. It helps the surgeons to predict the changes after surgical operations and their consequences.

Keywords: growth, soft tissue, residual stress, finite element method

Procedia PDF Downloads 355

Authors: Chi Zhang, Jun Jiang

Abstract:

Constitutive equations are very important in finite element (FE) modeling, and the accuracy of the material constants in the equations have significant effects on the accuracy of the FE models. A wide range of constitutive equations are available; however, fitting the material constants in the constitutive equations could be complex and time-consuming due to the strong non-linearity and relationship between the constants. This work will focus on the development of a set of unified MATLAB programs for fitting the material constants in the constitutive equations efficiently. Users will only need to supply experimental data in the required format and run the program without modifying functions or precisely guessing the initial values, or finding the parameters in previous works and will be able to fit the material constants efficiently.

Keywords: constitutive equations, FE modelling, MATLAB program, non-linear curve fitting

Procedia PDF Downloads 99

Authors: Santiago Montoya-Ospina, Raúl E. Jiménez-Mejía, Rosa Elvira Correa Gutiérrez

Abstract:

An experimental methodology is proposed for determining hemi-anechoic characteristics of an engineering acoustic room built at the facilities of Universidad Nacional de Colombia to evaluate the free-field conditions inside the chamber. Experimental results were compared with theoretical ones in both, the source and the sound propagation inside the chamber. Acoustic source was modeled by using monopole radiation pattern from punctual sources and the image method was considered for dealing with the reflective plane of the room, that means, the floor without insulation. Finite-difference time-domain (FDTD) method was implemented to calculate the sound pressure value at every spatial point of the chamber. Comparison between theoretical and experimental data yields to minimum error, giving satisfactory results for the hemi-anechoic characterization of the chamber.

Keywords: acoustic impedance, finite-difference time-domain, hemi-anechoic characterization

Procedia PDF Downloads 163

Authors: V. Ložar, N. Hadžić, T. Opetuk, R. Keser

Abstract:

The development of digital twins strongly depends on efficient algorithms and their capability to mirror real-life processes. Nowadays, such efforts are required to establish factories of the future faced with new demands of custom-made production. The ship hull processes face these challenges too. Therefore, it is important to implement design and evaluation approaches based on production system engineering. In this study, the recently developed finite state method is employed to describe the stell hull process as a platform for the implementation of digital twinning technology. The application is justified by comparing the finite state method with the analytical approach. This method is employed to rebuild a model of a real shipyard ship hull process using a combination of serial and splitting lines. The key performance indicators such as the production rate, work in process, probability of starvation, and blockade are calculated and compared to the corresponding results obtained through a simulation approach using the software tool Enterprise dynamics. This study confirms that the finite state method is a suitable tool for digital twinning applications. The conclusion highlights the advantages and disadvantages of methods employed in this context.

Keywords: digital twin, finite state method, production system engineering, shipyard

Procedia PDF Downloads 99

Authors: Suman Paik, P. V. Durgaprasad, Bijan K. Dutta

Abstract:

A study combining experimental and numerical investigation on the deformation behaviour of single crystals of copper is presented in this paper. Cylindrical samples were cut in specific orientations from high purity copper single crystal and subjected to uniaxial compression loading at quasi-static strain rate. The stress-strain curves along two different crystallographic orientations were then extracted. In order to study and compare the deformation responses, a single crystal plasticity model incorporating non-Schmid effects was developed assuming cross-slip plays an important role in orientation of the material. By making use of crystal plasticity finite element method, the model was applied to investigate the orientation dependence of the stress-strain behaviour of two crystallographic orientations. Finally, details of slip activities of deformed crystals were investigated by linking the orientation of slip lines with the theoretical traces of possible crystallographic planes. The experimentally determined active slip modes were matched with those determined by simulations.

Keywords: crystal plasticity, modelling, non-Schmid effects, finite elements, finite strain

Procedia PDF Downloads 213

Authors: Bassim Shaheen Bachy, Jörg Franke

Abstract:

In this study, a three dimensional numerical heat transfer model has been used to simulate the laser structuring of polymer substrate material in the Three-Dimensional Molded Interconnect Device (3D MID) which is used in the advanced multi-functional applications. A finite element method (FEM) transient thermal analysis is performed using APDL (ANSYS Parametric Design Language) provided by ANSYS. In this model, the effect of surface heat source was modeled with Gaussian distribution, also the effect of the mixed boundary conditions which consist of convection and radiation heat transfers have been considered in this analysis. The model provides a full description of the temperature distribution, as well as calculates the depth and the width of the groove upon material removal at different set of laser parameters such as laser power and laser speed. This study also includes the experimental procedure to study the effect of laser parameters on the depth and width of the removal groove metal as verification to the modeled results. Good agreement between the experimental and the model results is achieved for a wide range of laser powers. It is found that the quality of the laser structure process is affected by the laser scan speed and laser power. For a high laser structured quality, it is suggested to use laser with high speed and moderate to high laser power.

Keywords: laser structuring, simulation, finite element analysis, thermal modeling

Procedia PDF Downloads 349

Authors: Malorzata Kmiotek, Anna Kucaba-Pietal, Robert Smusz

Abstract:

Due to the miniaturisation process, in many technical devices, microchannels are used in cooling systems. Because of the small size of microchannels, the flow inside is laminar, which caused a slow heat exchange. In order to intensify the heat exchange, the flow must be disturbed, for example, by introducing obstacles. We present results on the influence of a thin obstacle, placed on microchannel wall, on the fluid and heat flow in the aspect of their use by constructors of heat exchangers. The obstacle is called 'thin' when its geometrical parameter (o=w/h, w- width, h - height of the obstacle) satisfies inequality: o < 0.5. In this work, we report numerical results on heat and mass transfer in the microchannels of 400 micrometer height (H - height of the microchannel), where thin obstacles are immersed on the walls, to disturb the flow. The Reynolds number of the flow in microchannel varies between 20 and 200 and is typical for the flow in micro heat exchangers. The equations describing the fluid and heat flows in microchannels were solved numerically by using the finite element method with an application of CFD&FSI package of ADINA R&D, Inc. 9.4 solver. In the case of flows in the microchannels with sequences of thin rectangular obstacles placed on the bottom and the top wall of a microchannel, the influence of distances s (s is the distance between two thin obstacles) and heights of obstacles on the fluid and heat transfer was investigated. Thermal and flow conditions of the application area of microchannels in electronic cooling systems, i.e., wall temperature of 60 °C, the fluid temperature of 20°C were used to solve equations. Additionally, the distance s between the thin obstacles in microchannels as a multiple of the amount of the channel height was determined. Results show that placing thin obstacles on microchannel walls increase the length of recirculation zones of the flow and improves the heat transfer.

Keywords: Finite Element Method, heat transfer, mechanical engineering, microchannel

Procedia PDF Downloads 134

Authors: Sudip Sudhir Godbole

Abstract:

In transmission line surface potential gradient is a critical design parameter for planning overhead line, as it determines the level of corona loss (CL), radio interference (RI) and audible noise (AN).With increase of transmission line voltage level bulk power transfer is possible, using bundle conductor configuration used, it is more complex to find accurate surface stress in bundle configuration. The majority of existing models for surface gradient calculations are based on analytical methods which restrict their application in simulating complex surface geometry. This paper proposes a novel technique which utilizes both analytical and numerical procedure to predict the surface gradient. One of 400 kV transmission line configurations has been selected as an example to compare the results for different methods. The different strand shapes are a key variable in determining.

Keywords: surface gradient, Maxwell potential coefficient method, market and Mengele’s method, successive images method, charge simulation method, finite element method

Procedia PDF Downloads 538

Authors: Alicia Kamischke, Souhail Elhouar, Yasser Khodair

Abstract:

The American Institute of Steel Construction (AISC) Specification (360-10) provides equations for calculating the capacity of a W-shaped steel member to resist concentrated forces applied to its flange. In the case of flange local bending, the capacity equations were primarily formulated for an interior point along the member, which is defined to be at a distance larger than ten flange thicknesses away from the member’s end. When a concentrated load is applied within ten flange thicknesses from the member’s end, AISC requires a fifty percent reduction to be applied to the flange bending capacity. This reduction, however, is not supported by any research. In this study, finite element modeling is used to investigate the actual reduction in capacity near the end of such a steel member. The results indicate that the AISC equation for flange local bending is quite conservative for forces applied at less than ten flange thicknesses from the member’s end and a new equation is suggested for the evaluation of available flange local bending capacity within that distance.

Keywords: flange local bending, concentrated forces, column, flange capacity

Procedia PDF Downloads 689

Authors: Jeong-Min Lee, Hyunwoo Song, Yonseok Kim, Junghan Yun, Jungin Byun, Jae-Mean Koo, Chang-Sung Seok

Abstract:

The Thermal Barrier Coating (TBC) prevents heat directly transferring from the high-temperature flame to the substrate. Top coat and bond coat compose the TBC and top coat consists of a ceramic and bond coat increases adhesion between the top coat and the substrate. The TBC technology drops the substrate surface temperature by about 150~200°C. In addition, the TBC system has a cooling system to lower the blade temperature by the air flow inside the blade. Then, as a result, the thermal gradient occurs inside the blade by cooling. Also, the internal stress occurs due to the difference in thermal expansion. In this paper, the finite element analyses (FEA) were performed and stress changes were derived according to the thermal gradient of the TBC system. The stress was increased due to the cooling, but difference of the stress between the top coat and bond coat was decreased. So, delamination in the interface between top coat and bond coat.

Keywords: gas turbine blade, Thermal Barrier Coating (TBC), thermal gradient, Finite Element Analysis (FEA)

Procedia PDF Downloads 608

Authors: S. H. Bae, D. W. Cho, W. B. Jeong, J. R. Cho

Abstract:

Fatigue fracture of an aluminum welded structure is a phenomenon frequently occurring from pores in a weld. In order to grasp the state of the welded structure in operation in real time, the acceleration signal of the structure is measured. At this time, the vibration characteristic of the signal according to the fatigue load is an important parameter of the state diagnosis. This paper was an experimental study on the variation of vibration characteristics of welded beams with vibration aging (especially bending vibration). First simple beams were produced according to welding conditions. Each beam was vibrated and measured beam's PSD (power spectral density) according to the degree of aging. Also, modal testing was conducted to compare the transfer functions of welded beams. Testing result shows that the natural frequencies of the beam changed with the vibration aging due to the change of stiffness in welding part and its stiffness was estimated by the finite element method.

Keywords: modal testing, natural frequency, vibration aging, welded structure

Procedia PDF Downloads 483

Authors: Manabu Sasajima, Takao Yamaguchi, Yoshio Koike, Mitsuharu Watanabe

Abstract:

This paper discusses the propagation of sound waves in air, specifically in narrow rectangular pathways of an occluded-ear simulator for acoustic measurements. In narrow pathways, both the speed of sound and the phase of the sound waves are affected by the damping of the air viscosity. Herein, we propose a new finite-element method (FEM) that considers the effects of the air viscosity. The method was developed as an extension of existing FEMs for porous, sound-absorbing materials. The results of a numerical calculation for a three-dimensional ear-simulator model using the proposed FEM were validated by comparing with theoretical lumped-parameter modeling analysis and standard values.

Keywords: ear simulator, FEM, simulation, viscosity

Procedia PDF Downloads 444