An Anisotropic Hardening Model for Springback Prediction
|
|
- Kristopher Nicholson
- 5 years ago
- Views:
Transcription
1 An Anisotropic Harening Moel for Springback Preiction Danielle Zeng an Z. Ceric Xia Scientific Research Laboratories For Motor Company Dearborn, MI 48 Abstract. As more Avance High-Strength Steels (AHSS are heavily use for automotive boy structures an closures panels, accurate springback preiction for these components becomes more challenging because of their rapi harening characteristics an ability to sustain even higher stresses. In this paper, a moifie Mroz harening moel is propose to capture realistic Bauschinger effect at reverse loaing, such as when material passes through ie raii or rawbea uring sheet metal forming process. his moel accounts for material anisotropic yiel surface an nonlinear isotropic/kinematic harening behavior. Material tension/compression test ata are use to accurately represent Bauschinger effect. he effectiveness of the moel is emonstrate by comparison of numerical an experimental springback results for a DP6 straight U-channel test. INRODUCION Springback is an important issue in the sheet metal part an ie esign. In orer to obtain a final part shape matching its esign intent, elimination or correction of springback has to be mae uring esign stages. Otherwise it will prouce imensional eviations an result in assembly ifficulties. raitionally, ie esign engineers take great effort on trial-an-error to either reuce or compensate for springback. A number of ie re-cuts are often neee uring tryout in orer to obtain a imensionally accurate part. his approach is not only time consuming but also costly. Because of weight reuction efforts to meet fuel economy pressure an increase vehicle safety requirement, the use of Avance High Strength Steel (AHSS sheets such as ual-phase an RIP alloys becomes more wiesprea. However, their higher yiel strength an rapi work harening make springback more ifficult to anticipate in avance an harer to control in tryout. In recent years, researchers an stamping engineers start to make ie compensation by using computer simulations before cutting physical ies [] []. In such a process, springback analysis is carrie out to obtain irections an magnitues for ie compensation. Accoringly accurate springback preiction becomes a key to the success of any ie compensation whose algorithm relies on numerically preicte springback. Unfortunately accurate springback preiction remains a critical challenge for both FEA coe evelopers an en users. As we know, springback is ue to the release of resiual stresses accumulate uring the forming stage. here are lots of factors affecting the stress preiction accuracy in numerical simulations, such as time integration schemes (imicit vs. exicit, element formulations, contact algorithms, material moels etc. In this paper, the attention is focuse on the material moeling aspect with the aim to evelop a more realistic constitutive relationship which hopefully represents material eformation behavior more accurately. Isotropic harening is the simest an most popular harening moel use toay. It offers reasonable approximation for monotonic loaing cases an is easy to imement. However, when materials experience unloaing an reverse loaing, their yiel stresses in reverse loaing is usually lower than those in the case of monotonic loaing (as shown in Figure. his phenomenon is generally calle Bauchinger effect. A purely kinematic harening moel was first introuce by Prager an Ziegler where the yiel surface translates in the stress space as material yiels. However this sime moel ignores harening effects in other stress irections. In sheet metal forming process, materials experience very comicate eformation. Cyclic bening an unbening occur when the sheet passes through a rawbea or a ie raius. In orer to represent the cyclic behavior more realistically, harening rules combining both isotropic harening an kinematic harening are later evelope by various researchers [-4]. Noticeably among them is 4
2 the multi-yiel surface moel propose by Mroz [5-6] which successfully moels the non-linear harening behavior an the smooth transition from elastic to astic eformation. his moel introuces the concept of a fiel of work-harening mouli to moel the nonlinear harening behavior instea of a single moulus use in most of other kinematic moels. Chu [7] generalize Mroz s iscrete multie yiel surface concept into a continuous fiel of yiel surfaces, an later ang [8] an ang etc. [9] use the moel to analyze sheet metal formability an springback. MODIFIED MROZ MODEL he soli line OABCDEF in Figure shows the stress-strain relationship of a sime uniaxial loaing, unloaing an reverse loaing behavior accoring to the original Mroz moel. he moel was further evelope by Chu [7] where the iscrete multie yiel surfaces were generalize into a continuous fiel of yiel surfaces an is more suitable for FEA imementations. Intereste reaers shoul fin etails in [7]. Figure Illustration of Bauschinger effect he Morz's moel has several avantages over other more comex harening moels. It captures anisotropic harening behavior nicely uring reverse loaing an reuces to isotropic harening for monotonic loaing. One of its most appealing characteristics for sheet metal forming analysis an springback in particular is that the formulation oes not require any aitional experimental tests beyon stanar tensile curves. In fact there are no extra parameters neee to fit the moel. However, the moel also assumes that the material elastic unloaing an reverse loaing region is constant an twice as the initial yiel stress. While this might be true for some low yiel alloys, experimental evience suggests it is not the case in general, an Han etc. [] foun that the amount of Bauschinger effect actually epens on the harening magnitue before the loaing is reverse. Accoringly a moifie Mroz moel is propose in this paper to take into account this variable Bauschinger effect. he paper is organize as follows. First, the new moel incorporating more accurate Bauschinger effect is escribe. An then its constitutive integration algorithm is erive an imemente in commercial FEA coe. Last, the propose moel is apie to a sime test case for springback preiction. Experiments are conucte an results are compare with numerical preiction to emonstrate the new moel's apicability. Figure Uniaxial stress-strain curves in Mroz moel It can be seen from Figure that line CD represents the material elastic unloaing. Point D is the reverse yiel starting point which represents the compressive yiel strength when the material starts to unloa at point C. In Mroz moel, it assumes that the elastic region is constant an represente by the initial yiel stress. he elastic range uring reverse loaing equals to twice of the initial yiel strength. hus when the material starts to unloa at higher strength, lower compressive yiel strength is obtaine. As shown in Figure, the soli pink line is the reversal of compressive yiel strength curve at the corresponing reverse point accoring to Mroz moel for a typical DP6 material. It can be seen that as the tensile strength increases, the corresponing compression strength ecreases. For the material with higher work harening, it is conceivable that the reverse yiel point will go to positive region, i.e., the moel preicts that the material might go into reverse astic loaing when the material is release from uniaxial tension even there is no external forces apie. his is contraictory to what was observe in real material compression testing []. Figure 4 gives the compressive yiel strength at various reversal points obtaine from compression tests for ifferent graes of steels. he ot shows that when the loa reverses at higher strength point, the compressive yiel strength is higher or at least in the same level. Accoring to the testing ata, the 4
3 compressive yiel strength curve shoul follow the otte pink line instea of the soli line in Figure as preicte by Mroz moel. Clearly the Mroz moel as it is can not aequately represent material eformation behavior at reverse loaing. rue stress (Mpa.5% Reverse Yiel Stength (MPa tensile Mroz compression testing rue strain Figure Comparison of material stress-strain curves Uncertainty is ± stanar eviation Strength at Reversal Point (MPa LC HSLA DP AKDQ DDQ DQSK IF-Rephos HS44 DP5 DP6 DP8 Figure 4 Yiel strength (.5% offset on reversal as a function of steel strength before reversal (Courtesy of C. Van yne In this paper, the Mroz moel is moifie to aress its rawback as iscusse above while preserving its original formulation. In the moifie moel, material eformation follows the same rule as that of Mroz moel at monotonic loaing. However, when the unloaing is initiate, the elastic region is no long constant. he size of the elastic region is etermine through reverse compression test. hus, in Figure, instea of line CD, the elastic unloaing curve becomes CD', where D'is the compressive yiel strength at the reverse point C an can be etermine through stanar compression test. Base on this moification, the size of the elastic region enote by σ B is a function of equivalent astic strain ε an can be expresse as: σ B ( σ Y + σ c σ y = ( where σ Y is the flow stress an σ c is the compressive yiel stress at corresponing reversal point. σ B can be rewritten in the following generalize form: σ B = cσ + ( c σ Y ( where σ is the initial yiel stress of the material, an c is a material parameter reflecting the Bauschinger effect an is consiere to be a function of the effective astic strain. It is easy to note that isotropic harening is a special case of this moel when c (thus σ B = σ Y, an the moel is reuce to the original Mroz moel when c is taken to be (c, with σ B = σ. In its most general form, c is a function of ε, where it can be expresse as: σ y σ c when σ y σ c ( ε = σ σ ( y when σ y = σ FORMULAION AND CONSIUIVE INEGRAION OF HE MODIFIED MROZ MODEL he material constitutive relationship an its integration for numerical imementation will be erive in this section. Since most sheet metals exhibit ifferent yiel strengths along ifferent irections, Hill's anisotropic yiel criterion is use to account the material anisotropic behavior. For the case of combine isotropic-kinematic harening, the astic yiel criterion can be expresse as []: f ( ( P( Y ( ε (4 = σ, σ, σ, σ, σ, σ is the stress vector in most general D stress state, where { } { α, α, α, α, α α } =, is the yiel surface center, ( Y ε is the flow stress which is a function of equivalent astic strain ε, an P is a 6x6 anisotropic asticity matrix in general D case an can be expresse in terms of the anisotropic r values in three orientations as: 4
4 r r + + r r r ( r9 + r r + ( + ( + r9 r r9 r r + r r9 P = r + r9 ( r + r9 ( r + ( r45 + ( r + r 9 r9 ( r + (5 Accoring to Hooke's law, the stress vector can be obtaine through elastic strain as: el = Dt+ (6 where D is the material elastic moulus matrix. Accoring to incremental astic flow theory, all the variables at time step t+ can be calculate once the eformation history at last time step t is known. hus el the elastic strain t+ can be expresse as: el el = t + (7 Here el t is the elastic strain at time t, an are the total strain increment an astic strain increment vectors, respectively. he associate flow rule gives: f = t+ (8 where is the equivalent astic strain increment. Also, uring the astic loaing, ( f (9 where * is usually terme as stress preictor: * el = D ( t + ( he stress vector can be re-written as: = M (4 ( * t+ t+ where M = [ I + DP] Y ε (5 he position of active yiel surface center at time t+ can be etermine accoring to the rule of Mroz moel. As shown in Figure 5, F t is the active yiel surface with center O at time t. F I is the inactive yiel surface with center O I in memory an tangent to surface F t at point P. F t+ is the active yiel surface with center O'at time step t+. β is a unit vector representing the moving irection of the center of current yiel surface an can be expresse as: = (6 ( I I t t P ( I t hen the amount of active yiel surface center movement from time t to t+ can be expresse as: = Y (7 where Y is the increment of the active yiel surface raius from time t to t+ an is a function of. Substituting Equation (4 into (9, Equation (9 becomes hus, from Equations (4 an (9, we have f = P( Y t+ herefore t+ ( f * * ( [ ( ] [ ( ] ( ε = P M Y ε = M (8 t+ t+ is the only unknown in this non-linear equation an can be solve numerically. In this stuy, Newton- Raphson metho is use to solve the equation. ε = P( Y t+ ( Substituting Equations (7 an ( into (6, we get: * ε = DP( t+ Y t+ ( Figure 5 Movement of the active yiel surface 44
5 It is also necessary to obtain material stiffness matrix if imicit time integration FEA is aopte. It can be expresse as, after some algebra: = D + Y f f f ( + ' (9 Figure 7 Springback of DP6 U channel NUMERICAL EXAMPLE Base on the formulation erive above, a user subroutine was evelope to compute the elasticastic constitutive equation, an has been imemente in imicit commercial software ABAQUS/Stanar. he user subroutine can be use for D soli, ane stress shell, D ane stress an ane strain elements. A straight U channel raw test was conucte to test the ability of the propose moel in preicting springback. Figure 6 shows the geometry of tooling set up. he upper ie raius is mm an the lower ie raius is 6mm. he punch with is 5mm with a.5mm raius. he raius of the rawbea is 6mm. he material properties of DP6 are liste in able. he tensile curve, which is obtaine from uniaxial tensile test is shown in Figure 8. In the physical testing, the blank was wrappe by a teflon sheet to reuce the friction. In the numerical simulation,.8 is use as the friction coefficient. able DP6 material properties Material E ν R R 45 R 9 σ y DP6 GPa MPa rue Stress (Mpa tensile Mroz compression testing mo_mroz rue astic strain Figure 8 DP6 uniaxial tensile curve an compressive yiel strength on reversal as a function of true strain before reversal Figure 6 U channel test set up A variety of tests were conucte with variations of materials, blank holer forces, an with an without rawbeas. he case presente in this paper is a.5 mm DP6 blank with rawbea an blank holer force of KN. he blank size is 4mm x mm. he raw epth is 7mm. Figure 7 shows the picture of testing piece in this loaing case after springback. It can be seen that DP6 has very big springback both in the wall open an sie wall curl. he material compression properties use for simulation in ifferent moels are also shown in Figure 8. he soli pink an green curves are the compression yiel stress at corresponing reversal points use by Mroz moel an the moifie Mroz moel, respectively. In this exame, the compressive yiel stress use in moifie Mroz moel is obtaine by offsetting the testing reverse yiel curve (from Figure 4 own to the point that compression yiel stress is the same as the tensile yiel stress at zero strain reverse point. In this case, the c value in Equation ( is about.5. he testing curve is not irectly use in this exame since the compression yiel stress was obtaine by a.5% offset in the original testing ata an generally.% offset is use. 45
6 Figure 9 shows the comparison of springback among the testing result an the preictions by using ifferent material moels. he anisotropic yiel function was accounte in all three simulation material moels. In Mroz moel an moifie Mroz moel, the combine isotropic-kinematic harening rule was use an followe the reverse compression yiel curves shown in Figure 8. It is shown in Figure 9 that the preiction by Mroz moel gives largest eviation from testing ata while that by moifie Mroz moel goes between the isotropic harening moel an Mroz moel, an it is closest to the testing result. his inicates that the Bauschinger effect is actually not as big as assume by the original Mroz moel. testing isotropic mroz mo_mroz Figure 9 Comparison of springback results DISCUSSIONS A moifie Mroz moel is presente in this paper to preict the springback for Avance High Strength Steel. he moel combines the material isotropic harening an kinematic harening, an also consiere the anisotropic yiel criterion. he moeling of material reverse loaing behavior is incorporate to aequately represent actual Bauschinger effect uner cyclic loaing case. he springback preiction by the propose moel for the DP6 U channel correlates with testing ata very well. Work is uner way to apy the new moel to comex prouction parts where Bauschinger effect is expecte to have a larger influence on springback than the sime D exame presente in this paper. In aition, more reverse compression test ata are neee to valiate the propose moel. ACKNOWLEDGEMEN he authors woul like to thank Mr. Craig Miller for conucting the U channel raw testing an analyzing the testing ata, an Prof. C. Van yne of Colorao School of Mines for proviing compression testing ata in Figure 4. REFERENCES. L. Wu, 'ooling mesh generation technique for iterative FEM ie surface esign algorithm to compensate for springback in sheet metal stamping', Engineering Computations, Vol. 4, No.6, 997, Lumin Geng, homas Oetjens an Chung-Yeh Sa, 'Springback preiction with LS-DYNA an ie face compensation of Aluminum hoo inner', SAE J. L. Chaboche, 'ime-inepenent constitutive theories for cyclic asticity', International Journal of Plasticity, Vol., No., 986, K. Chung, M. Lee, D. Kim, C. Kim, M. L. Wenner, F. Barlat, 'Spring-back evaluation of automotive sheets base on isotropic-kinematic harening laws an nonquaratic anisotropic yiel functions, Part I: theory an formulation', International Journal an Plasticity, Vol., 5, Z. Mroz, 'On the escription of anisotropic work harening', J. Mech. Phys. Solis, Vol. 5, 967, Akhtar S. Khan an Sujian Huang, Continuum heory of Plasticity, Awiley-Interscience Publication, C.C. Chu, 'A three-imensional moel of anisotropic harening in metals an its apication to the analysis of sheet metal formability', J. Mech. Phys. Solis, Vol., No., 984, S.C. ang, 'An anisotropic harening rule for the analysis of sheet metal forming operations', Avance echnology of Plasticity, Vol., 99, S.C. ang, Z.C. Xia an F. Ren, 'Apication of the raial return metho to compute stress increments from Mroz's harening rule', J. Engr. Mat. ech., Vol.,, K. Han, C. J. Van yne an B.S. Levy, 'Bauschinger effect response of automotive sheet steel', SAE Z.C. Xia, 'A General Anisotropic Plasticity Moel for Sheet Metals', For echnical report 46
VUMAT for Fabric Reinforced Composites
VUMAT or Fabric Reinorce Composites. Introuction This ocument escribes a constitutive mo or abric reinorce composites that was introuce in Abaqus/Exicit 6.8. The mo has been imemente as a built-in VUMAT
More informationA simple model for the small-strain behaviour of soils
A simple moel for the small-strain behaviour of soils José Jorge Naer Department of Structural an Geotechnical ngineering, Polytechnic School, University of São Paulo 05508-900, São Paulo, Brazil, e-mail:
More informationStrength Analysis of CFRP Composite Material Considering Multiple Fracture Modes
5--XXXX Strength Analysis of CFRP Composite Material Consiering Multiple Fracture Moes Author, co-author (Do NOT enter this information. It will be pulle from participant tab in MyTechZone) Affiliation
More informationFinite element analysis of electromagnetic bulging of sheet metals
International Journal of Scientific & Engineering Research Volume 3, Issue 2, Febraury-212 1 Finite element analysis of electromagnetic bulging of sheet metals Ali M. Abelhafeez, M. M. Nemat-Alla, M. G.
More information'HVLJQ &RQVLGHUDWLRQ LQ 0DWHULDO 6HOHFWLRQ 'HVLJQ 6HQVLWLYLW\,1752'8&7,21
Large amping in a structural material may be either esirable or unesirable, epening on the engineering application at han. For example, amping is a esirable property to the esigner concerne with limiting
More information3-D FEM Modeling of fiber/matrix interface debonding in UD composites including surface effects
IOP Conference Series: Materials Science an Engineering 3-D FEM Moeling of fiber/matrix interface eboning in UD composites incluing surface effects To cite this article: A Pupurs an J Varna 2012 IOP Conf.
More informationOn Springback Prediction In Stamping Of AHSS BIW Components Utilizing Advanced Material Models
On Springback Prediction In Stamping Of AHSS BIW Components Utilizing Advanced Material Models Ming F. Shi and Alex A. Konieczny United States Steel Corporation Introduction Origin of Springback AHSS Springback
More informationOF CHS. associated. indicate. the need. Rio de Janeiro, Brazil. a) Footbridge Rio. d) Maria Lenk. CHS K joints
EUROSTEEL 2, August 3 September 2, 2, Buapest, Hungary A NUMERICAL EVALUATION OF CHS T JOINTS UNDER AXIAL LOADS Raphael S. a Silva a, Luciano R. O. e Lima b, Pero C. G. a S. Vellasco b, José G. S. a Silva
More informationSimulation of Angle Beam Ultrasonic Testing with a Personal Computer
Key Engineering Materials Online: 4-8-5 I: 66-9795, Vols. 7-73, pp 38-33 oi:.48/www.scientific.net/kem.7-73.38 4 rans ech ublications, witzerlan Citation & Copyright (to be inserte by the publisher imulation
More informationMULTISCALE FRICTION MODELING FOR SHEET METAL FORMING
MULTISCALE FRICTION MODELING FOR SHEET METAL FORMING Authors J. HOL 1, M.V. CID ALFARO 2, M.B. DE ROOIJ 3 AND T. MEINDERS 4 1 Materials innovation institute (M2i) 2 Corus Research Centre 3 University of
More informationEstimation of hardness by nanoindentation of rough surfaces
Journal of MATERIALS RESEARCH Welcome Comments Help Estimation of harness by nanoinentation of rough surfaces M. S. Bobji an S. K. Biswas Department of Mechanical Engineering, Inian Institute of Science,
More informationInternational Conference on Advances in Energy, Environment and Chemical Engineering (AEECE-2015)
International Conference on Avances in Energy, Environment an Chemical Engineering (AEECE-2015) Stuy on Damage Characteristic of Unergroun Cavern Blasting Excavation Base on Dynamic Damage Constitutive
More informationModelling dowel action of discrete reinforcing bars in cracked concrete structures
Title Moelling owel action of iscrete reinforcing bars in cracke concrete structures Author(s) Kwan, AKH; Ng, PL; Lam, JYK Citation The 2n International Symposium on Computational Mechanics an the 12th
More informationThermal conductivity of graded composites: Numerical simulations and an effective medium approximation
JOURNAL OF MATERIALS SCIENCE 34 (999)5497 5503 Thermal conuctivity of grae composites: Numerical simulations an an effective meium approximation P. M. HUI Department of Physics, The Chinese University
More informationinflow outflow Part I. Regular tasks for MAE598/494 Task 1
MAE 494/598, Fall 2016 Project #1 (Regular tasks = 20 points) Har copy of report is ue at the start of class on the ue ate. The rules on collaboration will be release separately. Please always follow the
More informationThe Phenomenon of Anomalous Rock Embrittlement
https://papers.acg.uwa.eu.au/p/574_29_tarasov/ B.G. Tarasov A.V. Dyskin School of Civil an Resource Engineering The University of Western Australia The paper analyses a phenomenon of rock behaviour - the
More informationChapter 2 Governing Equations
Chapter 2 Governing Equations In the present an the subsequent chapters, we shall, either irectly or inirectly, be concerne with the bounary-layer flow of an incompressible viscous flui without any involvement
More informationA PROCEDURE FOR DETERMINATION OF THE ALTERNAT MODEL PARAMETERS
1. Mohamme Y. FATTAH,. Omar. F. S. al DAMLUJI, 3. Yousif. J. al SHAKARCHI A PROCEDURE FOR DETERMINATION OF THE ALTERNAT MODEL PARAMETERS 1. BUILDING AND CONSTRUCTION DEPARTMENT, UNIVERSITY OF TECHNOLOGY,
More informationNonlinear Dielectric Response of Periodic Composite Materials
onlinear Dielectric Response of Perioic Composite aterials A.G. KOLPAKOV 3, Bl.95, 9 th ovember str., ovosibirsk, 639 Russia the corresponing author e-mail: agk@neic.nsk.su, algk@ngs.ru A. K.TAGATSEV Ceramics
More informationThe Principle of Least Action
Chapter 7. The Principle of Least Action 7.1 Force Methos vs. Energy Methos We have so far stuie two istinct ways of analyzing physics problems: force methos, basically consisting of the application of
More informationMath Notes on differentials, the Chain Rule, gradients, directional derivative, and normal vectors
Math 18.02 Notes on ifferentials, the Chain Rule, graients, irectional erivative, an normal vectors Tangent plane an linear approximation We efine the partial erivatives of f( xy, ) as follows: f f( x+
More informationPredictive Control of a Laboratory Time Delay Process Experiment
Print ISSN:3 6; Online ISSN: 367-5357 DOI:0478/itc-03-0005 Preictive Control of a aboratory ime Delay Process Experiment S Enev Key Wors: Moel preictive control; time elay process; experimental results
More informationVerification of cylindrical interference fits under impact loads with LS-Dyna
Verification of cylinrical interference fits uner impact loas with LS-Dyna Prof. Dr.-Ing. elmut Behler, Jan Göbel, M.Eng. ochschule Mannheim, Paul-Wittsack-Straße 10, D-68163 Mannheim, Germany Steffen
More informationResilient Modulus Prediction Model for Fine-Grained Soils in Ohio: Preliminary Study
Resilient Moulus Preiction Moel for Fine-Graine Soils in Ohio: Preliminary Stuy by Teruhisa Masaa: Associate Professor, Civil Engineering Department Ohio University, Athens, OH 4570 Tel: (740) 59-474 Fax:
More information3.2 Shot peening - modeling 3 PROCEEDINGS
3.2 Shot peening - moeling 3 PROCEEDINGS Computer assiste coverage simulation François-Xavier Abaie a, b a FROHN, Germany, fx.abaie@frohn.com. b PEENING ACCESSORIES, Switzerlan, info@peening.ch Keywors:
More informationSurvey Sampling. 1 Design-based Inference. Kosuke Imai Department of Politics, Princeton University. February 19, 2013
Survey Sampling Kosuke Imai Department of Politics, Princeton University February 19, 2013 Survey sampling is one of the most commonly use ata collection methos for social scientists. We begin by escribing
More informationAPPROXIMATE SOLUTION FOR TRANSIENT HEAT TRANSFER IN STATIC TURBULENT HE II. B. Baudouy. CEA/Saclay, DSM/DAPNIA/STCM Gif-sur-Yvette Cedex, France
APPROXIMAE SOLUION FOR RANSIEN HEA RANSFER IN SAIC URBULEN HE II B. Bauouy CEA/Saclay, DSM/DAPNIA/SCM 91191 Gif-sur-Yvette Ceex, France ABSRAC Analytical solution in one imension of the heat iffusion equation
More informationAssessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Subjected to In-Plane Shear
Assessment of the Buckling Behavior of Square Composite Plates with Circular Cutout Sujecte to In-Plane Shear Husam Al Qalan 1)*, Hasan Katkhua 1) an Hazim Dwairi 1) 1) Assistant Professor, Civil Engineering
More informationCrack onset assessment near the sharp material inclusion tip by means of modified maximum tangential stress criterion
Focuse on Mechanical Fatigue of Metals Crack onset assessment near the sharp material inclusion tip by means of moifie maximum tangential stress criterion Onřej Krepl, Jan Klusák CEITEC IPM, Institute
More informationOptimization of Geometries by Energy Minimization
Optimization of Geometries by Energy Minimization by Tracy P. Hamilton Department of Chemistry University of Alabama at Birmingham Birmingham, AL 3594-140 hamilton@uab.eu Copyright Tracy P. Hamilton, 1997.
More information6 General properties of an autonomous system of two first order ODE
6 General properties of an autonomous system of two first orer ODE Here we embark on stuying the autonomous system of two first orer ifferential equations of the form ẋ 1 = f 1 (, x 2 ), ẋ 2 = f 2 (, x
More informationThe effect of nonvertical shear on turbulence in a stably stratified medium
The effect of nonvertical shear on turbulence in a stably stratifie meium Frank G. Jacobitz an Sutanu Sarkar Citation: Physics of Fluis (1994-present) 10, 1158 (1998); oi: 10.1063/1.869640 View online:
More informationELECTRON DIFFRACTION
ELECTRON DIFFRACTION Electrons : wave or quanta? Measurement of wavelength an momentum of electrons. Introuction Electrons isplay both wave an particle properties. What is the relationship between the
More informationTOWARDS THERMOELASTICITY OF FRACTAL MEDIA
ownloae By: [University of Illinois] At: 21:04 17 August 2007 Journal of Thermal Stresses, 30: 889 896, 2007 Copyright Taylor & Francis Group, LLC ISSN: 0149-5739 print/1521-074x online OI: 10.1080/01495730701495618
More informationSituation awareness of power system based on static voltage security region
The 6th International Conference on Renewable Power Generation (RPG) 19 20 October 2017 Situation awareness of power system base on static voltage security region Fei Xiao, Zi-Qing Jiang, Qian Ai, Ran
More informationAdhesive Wear Theory of Micromechanical Surface Contact
International Journal Of Computational Engineering esearch ijceronline.com Vol. Issue. hesive Wear Theory of Micromechanical Surface Contact iswajit era Department of Mechanical Engineering National Institute
More informationMarine gears load capacity of involute parallel axis spur and helical gears
(1990) (Rev.1 1994/ Corr. 1996) (Rev. Oct 013) (Rev.3 Oct 015) Marine gears loa capacity of involute parallel axis spur an helical gears.1 Basic principles - introuction an general influence factors.1.1
More informationApplications of First Order Equations
Applications of First Orer Equations Viscous Friction Consier a small mass that has been roppe into a thin vertical tube of viscous flui lie oil. The mass falls, ue to the force of gravity, but falls more
More information6. Friction and viscosity in gasses
IR2 6. Friction an viscosity in gasses 6.1 Introuction Similar to fluis, also for laminar flowing gases Newtons s friction law hols true (see experiment IR1). Using Newton s law the viscosity of air uner
More informationAdvanced friction modeling for sheet metal forming
Avance friction moeling for sheet metal forming Authors J.Hol a, M.V. Ci Alfaro b, M.B. e Rooij c, T. Meiners a Materials innovation institute (M2i) b Tata Steel Research, Development & Technology c University
More informationExperiment 2, Physics 2BL
Experiment 2, Physics 2BL Deuction of Mass Distributions. Last Upate: 2009-05-03 Preparation Before this experiment, we recommen you review or familiarize yourself with the following: Chapters 4-6 in Taylor
More informationChapter 6: Energy-Momentum Tensors
49 Chapter 6: Energy-Momentum Tensors This chapter outlines the general theory of energy an momentum conservation in terms of energy-momentum tensors, then applies these ieas to the case of Bohm's moel.
More information2-7. Fitting a Model to Data I. A Model of Direct Variation. Lesson. Mental Math
Lesson 2-7 Fitting a Moel to Data I BIG IDEA If you etermine from a particular set of ata that y varies irectly or inversely as, you can graph the ata to see what relationship is reasonable. Using that
More informationImpact Experimental Analysis and Computer Simulation Yucheng Liu Department of Mechanical Engineering, University of Louisville
Impact Experimental Analysis an Computer Simulation Yucheng Liu Department of Mechanical Engineering, University of Louisville Abstract In this paper, an automotive bumper system (a bumper connecte to
More informationConservation laws a simple application to the telegraph equation
J Comput Electron 2008 7: 47 51 DOI 10.1007/s10825-008-0250-2 Conservation laws a simple application to the telegraph equation Uwe Norbrock Reinhol Kienzler Publishe online: 1 May 2008 Springer Scienceusiness
More informationThis section outlines the methodology used to calculate the wave load and wave wind load values.
COMPUTERS AND STRUCTURES, INC., JUNE 2014 AUTOMATIC WAVE LOADS TECHNICAL NOTE CALCULATION O WAVE LOAD VALUES This section outlines the methoology use to calculate the wave loa an wave win loa values. Overview
More information05 The Continuum Limit and the Wave Equation
Utah State University DigitalCommons@USU Founations of Wave Phenomena Physics, Department of 1-1-2004 05 The Continuum Limit an the Wave Equation Charles G. Torre Department of Physics, Utah State University,
More informationLATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION
The Annals of Statistics 1997, Vol. 25, No. 6, 2313 2327 LATTICE-BASED D-OPTIMUM DESIGN FOR FOURIER REGRESSION By Eva Riccomagno, 1 Rainer Schwabe 2 an Henry P. Wynn 1 University of Warwick, Technische
More informationV = Flow velocity, ft/sec
1 Drag Coefficient Preiction Chapter 1 The ieal force acting on a surface positione perpenicular to the airflow is equal to a ynamic pressure, enote by q, times the area of that surface. Dynamic pressure
More informationCONTROL CHARTS FOR VARIABLES
UNIT CONTOL CHATS FO VAIABLES Structure.1 Introuction Objectives. Control Chart Technique.3 Control Charts for Variables.4 Control Chart for Mean(-Chart).5 ange Chart (-Chart).6 Stanar Deviation Chart
More informationNonPAS: A Program for Nonlinear Analysis of Flexible Pavements
International Journal of Integrate Engineering, Vol. 7 No. (05) pp. -8 NonPAS: A Program for Nonlinear Analysis of Flexible Pavements Ali eza Ghanizaeh,*, Arash Ziaie, Department of Civil Engineering,
More informationA Simple Model for the Calculation of Plasma Impedance in Atmospheric Radio Frequency Discharges
Plasma Science an Technology, Vol.16, No.1, Oct. 214 A Simple Moel for the Calculation of Plasma Impeance in Atmospheric Raio Frequency Discharges GE Lei ( ) an ZHANG Yuantao ( ) Shanong Provincial Key
More informationCrack-tip stress evaluation of multi-scale Griffith crack subjected to
Crack-tip stress evaluation of multi-scale Griffith crack subjecte to tensile loaing by using periynamics Xiao-Wei Jiang, Hai Wang* School of Aeronautics an Astronautics, Shanghai Jiao Tong University,
More informationINDIAN REGISTER OF SHIPPING CLASSIFICATION NOTES
INDIAN REGISTER OF SHIPPING CLASSIFICATION NOTES Marine Gears Calculation of Loa Capacity of Involute Parallel Axis Spur an Helical Gears January 05 January 05 Page of 9 CLASSIFICATION NOTES Marine Gears
More informationCURRENT ELECTRICITY Q.1
CUENT EECTCTY Q. Define Electric current an its unit.. Electric Current t can be efine as the time rate of flow of charge in a conuctor is calle Electric Current. The amount of flow of charge Q per unit
More informationA Novel Decoupled Iterative Method for Deep-Submicron MOSFET RF Circuit Simulation
A Novel ecouple Iterative Metho for eep-submicron MOSFET RF Circuit Simulation CHUAN-SHENG WANG an YIMING LI epartment of Mathematics, National Tsing Hua University, National Nano evice Laboratories, an
More informationRecommendations: Part 7: Transient Creep for service and accident conditions
Materials an Structures/Matériaux et Constructions, Vol. 31, June 1998, pp 290-295 RILEM TECHNICAL COMMITTEES RILEM TC 129-MHT: TEST METHODS FOR MECHANICAL PROPERTIES OF CONCRETE AT HIGH TEMPERATURES Recommenations:
More informationTable of Common Derivatives By David Abraham
Prouct an Quotient Rules: Table of Common Derivatives By Davi Abraham [ f ( g( ] = [ f ( ] g( + f ( [ g( ] f ( = g( [ f ( ] g( g( f ( [ g( ] Trigonometric Functions: sin( = cos( cos( = sin( tan( = sec
More informationVibration Analysis of Railway Tracks Forced by Distributed Moving Loads
IJR International Journal of Railway Vol. 6, No. 4 / December 13, pp. 155-159 The Korean Society for Railway Vibration Analysis of Railway Tracks Force by Distribute Moving Loas Sinyeob Lee*, Dongkyu Kim*,
More informationOptimal LQR Control of Structures using Linear Modal Model
Optimal LQR Control of Structures using Linear Moal Moel I. Halperin,2, G. Agranovich an Y. Ribakov 2 Department of Electrical an Electronics Engineering 2 Department of Civil Engineering Faculty of Engineering,
More informationEssential Considerations for Buckling Analysis
Worlwie Aerospace Conference an Technology Showcase, Toulouse, France, Sept. 24-26, 2001 Essential Consierations for Buckling Analysis 2001-120 Sang H. Lee MSC.Software Corporation, 2 MacArthur Place,
More informationAngles-Only Orbit Determination Copyright 2006 Michel Santos Page 1
Angles-Only Orbit Determination Copyright 6 Michel Santos Page 1 Abstract This ocument presents a re-erivation of the Gauss an Laplace Angles-Only Methos for Initial Orbit Determination. It keeps close
More informationExperimental Robustness Study of a Second-Order Sliding Mode Controller
Experimental Robustness Stuy of a Secon-Orer Sliing Moe Controller Anré Blom, Bram e Jager Einhoven University of Technology Department of Mechanical Engineering P.O. Box 513, 5600 MB Einhoven, The Netherlans
More informationDissipative numerical methods for the Hunter-Saxton equation
Dissipative numerical methos for the Hunter-Saton equation Yan Xu an Chi-Wang Shu Abstract In this paper, we present further evelopment of the local iscontinuous Galerkin (LDG) metho esigne in [] an a
More informationSome New Thoughts on the Multipoint Method for Reactor Physics Applications. Sandra Dulla, Piero Ravetto, Paolo Saracco,
Jeju, Korea, April 16-20, 2017, on USB 2017 Some New Thoughts on the Multipoint Metho for Reactor Physics Applications Sanra Dulla, Piero Ravetto, Paolo Saracco, Politecnico i Torino, Dipartimento Energia,
More informationP. A. Martin b) Department of Mathematics, University of Manchester, Manchester M13 9PL, United Kingdom
Time-harmonic torsional waves in a composite cyliner with an imperfect interface J. R. Berger a) Division of Engineering, Colorao School of Mines, Golen, Colorao 80401 P. A. Martin b) Department of Mathematics,
More informationSummary. Introduction
Optimal survey esign for marine borehole seismics Darrell Coles*, Schlumberger-Doll Research, Yi Yang, WesternGeco, Hugues Djikpesse, Michael Prange, Schlumberger-Doll Research, an Konstantin Osypov, WesternGeco
More informationSmectic-C tilt under shear in smectic-a elastomers
Smectic-C tilt uner shear in smectic-a elastomers Olaf Stenull an T. C. Lubensky Department of Physics an Astronomy, University of Pennsylvania, Philaelphia, Pennsylvania 19104, USA J. M. Aams an Mark
More informationIntroduction to variational calculus: Lecture notes 1
October 10, 2006 Introuction to variational calculus: Lecture notes 1 Ewin Langmann Mathematical Physics, KTH Physics, AlbaNova, SE-106 91 Stockholm, Sween Abstract I give an informal summary of variational
More informationSYNCHRONOUS SEQUENTIAL CIRCUITS
CHAPTER SYNCHRONOUS SEUENTIAL CIRCUITS Registers an counters, two very common synchronous sequential circuits, are introuce in this chapter. Register is a igital circuit for storing information. Contents
More informationSparse Reconstruction of Systems of Ordinary Differential Equations
Sparse Reconstruction of Systems of Orinary Differential Equations Manuel Mai a, Mark D. Shattuck b,c, Corey S. O Hern c,a,,e, a Department of Physics, Yale University, New Haven, Connecticut 06520, USA
More informationHomework 7 Due 18 November at 6:00 pm
Homework 7 Due 18 November at 6:00 pm 1. Maxwell s Equations Quasi-statics o a An air core, N turn, cylinrical solenoi of length an raius a, carries a current I Io cos t. a. Using Ampere s Law, etermine
More informationd dx But have you ever seen a derivation of these results? We ll prove the first result below. cos h 1
Lecture 5 Some ifferentiation rules Trigonometric functions (Relevant section from Stewart, Seventh Eition: Section 3.3) You all know that sin = cos cos = sin. () But have you ever seen a erivation of
More informationModule 5 Couplings. Version 2 ME, IIT Kharagpur
Moule 5 Couplings Version ME, IIT Kharagpur Lesson Design proceures for rigi an flexible rubber-bushe couplings Version ME, IIT Kharagpur Instructional Objectives At the en of this lesson, the stuents
More informationEfficient Macro-Micro Scale Coupled Modeling of Batteries
A00 Journal of The Electrochemical Society, 15 10 A00-A008 005 0013-651/005/1510/A00/7/$7.00 The Electrochemical Society, Inc. Efficient Macro-Micro Scale Couple Moeling of Batteries Venkat. Subramanian,*,z
More informationStable and compact finite difference schemes
Center for Turbulence Research Annual Research Briefs 2006 2 Stable an compact finite ifference schemes By K. Mattsson, M. Svär AND M. Shoeybi. Motivation an objectives Compact secon erivatives have long
More informationTime-Optimal Motion Control of Piezoelectric Actuator: STM Application
Time-Optimal Motion Control of Piezoelectric Actuator: STM Application Yongai Xu, Peter H. Mecl Abstract This paper exaes the problem of time-optimal motion control in the context of Scanning Tunneling
More informationConnections Between Duality in Control Theory and
Connections Between Duality in Control heory an Convex Optimization V. Balakrishnan 1 an L. Vanenberghe 2 Abstract Several important problems in control theory can be reformulate as convex optimization
More informationOptimum design of tuned mass damper systems for seismic structures
Earthquake Resistant Engineering Structures VII 175 Optimum esign of tune mass amper systems for seismic structures I. Abulsalam, M. Al-Janabi & M. G. Al-Taweel Department of Civil Engineering, Faculty
More informationNon-Conservative Stability Analysis of Hauger Types of Columns with Different Boundary Conditions
Proceeings of the Worl Congress on Engineering 08 Vol II WCE 08, July 4-6, 08, onon, U.K. on-conservative Stability nalysis of Hauger Types of Columns with Different Bounary Conitions S.. Fazelzaeh,. Tashakorian,
More informationMoist Component Potential Vorticity
166 JOURNAL OF THE ATMOSPHERIC SCIENCES VOLUME 60 Moist Component Potential Vorticity R. MCTAGGART-COWAN, J.R.GYAKUM, AND M. K. YAU Department of Atmospheric an Oceanic Sciences, McGill University, Montreal,
More informationElectromagnet Gripping in Iron Foundry Automation Part II: Simulation
www.ijcsi.org 238 Electromagnet Gripping in Iron Founry Automation Part II: Simulation Rhythm-Suren Wahwa Department of Prouction an Quality Engineering, NTNU Tronheim, 7051, Norway Abstract This paper
More informationLeast-Squares Regression on Sparse Spaces
Least-Squares Regression on Sparse Spaces Yuri Grinberg, Mahi Milani Far, Joelle Pineau School of Computer Science McGill University Montreal, Canaa {ygrinb,mmilan1,jpineau}@cs.mcgill.ca 1 Introuction
More informationImage Denoising Using Spatial Adaptive Thresholding
International Journal of Engineering Technology, Management an Applie Sciences Image Denoising Using Spatial Aaptive Thresholing Raneesh Mishra M. Tech Stuent, Department of Electronics & Communication,
More informationEVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL
386 EVALUATION OF LIQUEFACTION RESISTANCE AND LIQUEFACTION INDUCED SETTLEMENT FOR RECLAIMED SOIL Lien-Kwei CHIEN 1, Yan-Nam OH 2 An Chih-Hsin CHANG 3 SUMMARY In this stuy, the fille material in Yun-Lin
More informationModeling the effects of polydispersity on the viscosity of noncolloidal hard sphere suspensions. Paul M. Mwasame, Norman J. Wagner, Antony N.
Submitte to the Journal of Rheology Moeling the effects of polyispersity on the viscosity of noncolloial har sphere suspensions Paul M. Mwasame, Norman J. Wagner, Antony N. Beris a) epartment of Chemical
More information19 Eigenvalues, Eigenvectors, Ordinary Differential Equations, and Control
19 Eigenvalues, Eigenvectors, Orinary Differential Equations, an Control This section introuces eigenvalues an eigenvectors of a matrix, an iscusses the role of the eigenvalues in etermining the behavior
More informationTHE ACCURATE ELEMENT METHOD: A NEW PARADIGM FOR NUMERICAL SOLUTION OF ORDINARY DIFFERENTIAL EQUATIONS
THE PUBISHING HOUSE PROCEEDINGS O THE ROMANIAN ACADEMY, Series A, O THE ROMANIAN ACADEMY Volume, Number /, pp. 6 THE ACCURATE EEMENT METHOD: A NEW PARADIGM OR NUMERICA SOUTION O ORDINARY DIERENTIA EQUATIONS
More informationHow the potentials in different gauges yield the same retarded electric and magnetic fields
How the potentials in ifferent gauges yiel the same retare electric an magnetic fiels José A. Heras a Departamento e Física, E. S. F. M., Instituto Politécnico Nacional, México D. F. México an Department
More informationX-ray Diffraction from Materials
X-ray Diffraction from Materials 8 Spring Semester Lecturer; Yang Mo Koo Monay an Wenesay 4:45~6: Plastically ll eforme metal -Macro strain: uniform elastic strain over relatively large istance iffraction
More informationCode_Aster. Detection of the singularities and computation of a card of size of elements
Titre : Détection es singularités et calcul une carte [...] Date : 0/0/0 Page : /6 Responsable : Josselin DLMAS Clé : R4.0.04 Révision : 9755 Detection of the singularities an computation of a car of size
More informationECE 422 Power System Operations & Planning 7 Transient Stability
ECE 4 Power System Operations & Planning 7 Transient Stability Spring 5 Instructor: Kai Sun References Saaat s Chapter.5 ~. EPRI Tutorial s Chapter 7 Kunur s Chapter 3 Transient Stability The ability of
More informationarxiv:physics/ v4 [physics.class-ph] 9 Jul 1999
AIAA-99-2144 PROPULSION THROUGH ELECTROMAGNETIC SELF-SUSTAINED ACCELERATION arxiv:physics/9906059v4 [physics.class-ph] 9 Jul 1999 Abstract As is known the repulsion of the volume elements of an uniformly
More informationCombined Isotropic-Kinematic Hardening Laws with Anisotropic Back-stress Evolution for Orthotropic Fiber-Reinforced Composites
Combined Isotropic-Kinematic Hardening Laws with Antropic Back-stress Evolution for Orthotropic Fiber- Reinforced Composites Combined Isotropic-Kinematic Hardening Laws with Antropic Back-stress Evolution
More informationarxiv:hep-th/ v1 3 Feb 1993
NBI-HE-9-89 PAR LPTHE 9-49 FTUAM 9-44 November 99 Matrix moel calculations beyon the spherical limit arxiv:hep-th/93004v 3 Feb 993 J. Ambjørn The Niels Bohr Institute Blegamsvej 7, DK-00 Copenhagen Ø,
More informationBoth the ASME B and the draft VDI/VDE 2617 have strengths and
Choosing Test Positions for Laser Tracker Evaluation an Future Stanars Development ala Muralikrishnan 1, Daniel Sawyer 1, Christopher lackburn 1, Steven Phillips 1, Craig Shakarji 1, E Morse 2, an Robert
More informationFinite element analysis of the dynamic behavior of concrete-filled double skin steel tubes (CFDST) under lateral impact with fixed-simply support
Structures Uner Shock an Impact XIII 383 Finite element analysis of the ynamic behavior of concrete-fille ouble skin steel tubes (CFDST) uner lateral impact with fixe-simply support M. Liu & R. Wang Architecture
More informationOn Using Unstable Electrohydraulic Valves for Control
Kailash Krishnaswamy Perry Y. Li Department of Mechanical Engineering, University of Minnesota, 111 Church St. SE, Minneapolis, MN 55455 e-mail: kk,pli @me.umn.eu On Using Unstable Electrohyraulic Valves
More informationSchöck Isokorb type KST
Schöck Isokorb type Schöck Isokorb type Contents Page Element arrangements/connection layouts 288-289 Views/Dimensions 290-293 Design an capacity table 294 Torsion spring strength/notes on calculations
More informationTHE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE
Journal of Soun an Vibration (1996) 191(3), 397 414 THE VAN KAMPEN EXPANSION FOR LINKED DUFFING LINEAR OSCILLATORS EXCITED BY COLORED NOISE E. M. WEINSTEIN Galaxy Scientific Corporation, 2500 English Creek
More information