An Innovative Scheme to Make an Initial Guess for Iterative Optimization Methods to Calibrate Material Parameters of Strain-Hardening Elastoplastic Models

Author(s):  
Manouchehr Sanei ◽  
Philippe R. B. Devloo ◽  
Tiago L. D. Forti ◽  
Omar Durán ◽  
Erick S. R. Santos
1981 ◽  
Vol 21 (05) ◽  
pp. 551-557 ◽  
Author(s):  
Ali H. Dogru ◽  
John H. Seinfeld

Abstract The efficiency of automatic history matchingalgorithms depends on two factors: the computationtime needed per iteration and the number of iterations needed for convergence. In most historymatching algorithms, the most time-consumingaspect is the calculation of the sensitivitycoefficientsthe derivatives of the reservoir variables(pressure and saturation) with respect to the reservoirproperties (permeabilities and porosity). This paper presents an analysis of two methodsthe direct andthe variationalfor calculating sensitivitycoefficients, with particular emphasis on thecomputational requirements of the methods.If the simulator consists of a set of N ordinary differential equations for the grid-block variables(e.g., pressures)and there are M parameters forwhich the sensitivity coefficients are desired, the ratioof the computational efforts of the direct to thevariational method is N(M + 1)R = .N(N + 1) + M Thus, for M less than N the direct method is moreeconomical, whereas as M increases, a point isreached at which the variational method is preferred. Introduction There has been considerable interest in thedevelopment of automatic history matching algorithms.Although automatic history matching can offer significant advantages over trial-and-errorapproaches, its adoption has been somewhatlower than might have been anticipated when thefirst significant papers on the subject appeared. Oneobvious reason for the persistence of thetrial-and-error approach is that it does not requireadditional code development beyond that already involvedin the basic simulator, whereas automatic routinesrequire the appendixing of an iterative optimization routine to the basic simulator. Nevertheless, theinvestment of additional time in code developmentfor the history matching algorithm may be returned many fold during the actual history matchingexercise. In spite of the inherent advantages ofautomatic history matching, however, the automatic adjustment of the number of reservoir parameterstypically unknown even in a moderately sizedsimulation can require excessive amounts ofcomputation time. Therefore, it is of utmost importancethat an automatic history matching algorithm be asefficient as possible. Setting aside for the moment the issue of code complexity, the efficiency of analgorithm depends on two factors, the computationtime needed per iteration and the number ofiterations needed for convergence (whereconvergence is usually defined in terms of reaching acertain level of incremental change in either theparameters themselves or the objective function). Formost iterative optimization methods, the speed ofconvergence increases with the complexity of thealgorithm. SPEJ P. 551^


2007 ◽  
Vol 561-565 ◽  
pp. 1869-1874
Author(s):  
Quan Lin Jin ◽  
Yan Shu Zhang

A hybrid global optimization method combining the Real-coded genetic algorithm and some classical local optimization methods is constructed and applied to develop a special program for parameter identification. Finally, the parameter identification for both 26Cr2Ni4MoV steel and AZ31D magnesium alloy is carried out by using the program. A comparison of deformation test and numerical simulation shows that the parameter identification and the obtained two sets of material parameters are all available.


2020 ◽  
Author(s):  
Ahmad Sedaghat ◽  
Amir Mosavi

AbstractThe SIR type models are built by a set of ordinary differential equations (ODE), which are strongly initial value dependant. To fit multiple biological data with SIR type equations requires fitting coefficients of these equations by an initial guess and applying optimization methods. These coefficients are also extremely initial value-dependent. In the vast publication of these types, we hardly see, among simple to highly complicated SIR type methods, that these methods presented more than a maximum of two biological data sets. We propose a novel method that integrates an analytical solution of the infectious population using Weibull distribution function into any SIR type models. The Weibull-SIRD method has easily fitted 4 set of COVID-19 biological data simultaneously. It is demonstrated that the Weibull-SIRD method predictions for susceptible, infected, recovered, and deceased populations from COVID-19 in Kuwait and UAE are superior compared with SIRD original ODE model. The proposed method here opens doors for new deeper studying of biological dynamic systems with realistic biological data trends than providing some complicated, cumbersome mathematical methods with little insight into biological data’s real physics.


Author(s):  
Yinhui Zhang ◽  
Jian Shuai

Abstract As the main transportation carrier of oil and gas, pipelines play a very important role in the petroleum industry. When the crack-containing pipelines subjected to external loads, the cracks may propagate gradually, and result in serious failure eventually. Therefore, accurately obtaining the fracture toughness is very essential for the safety assessment of the crack-containing pipelines. However, the fracture toughness is not a material intrinsic parameter, but heavily depends on the constraint. To obtain the accurate relationship between the constraint and the fracture toughness for different materials, it is necessary to determine the effects of different material parameters on the change characteristics of the constraint and the fracture toughness. In this work, the commonly used pipelines steels are selected as the research materials. The SENB specimens and the complete Gurson model are used to conduct the simulation with ABAQUS. The material parameters analyzed include strain hardening exponent, yield strength and initial void volume fraction. The results show that for the thinner specimen, the higher strain hardening capacity will result in lower constraint. The higher strain hardening capacity will result in higher constraint for the thicker specimen. For the thinner specimen, the constraint is approximately the same for the materials with different yield strength. The constraint will decrease with the increase of yield strength for the thicker specimen. In the middle range of the thickness of specimen, higher initial void volume fraction will result in higher constraint. For the thicker and thinner specimen, the effect of initial void volume is very weak. As the increase of strain hardening capacity and yield strength, the decreasing degree of the fracture toughness becomes higher in the increasing process of the constraint. A higher initial void volume will result in a lower decreasing degree of the fracture toughness. All of the results indicate that the strain hardening capacity is the main factor affecting the constraint and the fracture toughness. The initial void volume also has a significant effect on the fracture toughness. For the relationship between the constraint and the fracture toughness, the main affecting factor is the strain hardening capacity.


Author(s):  
Digendranath Swain ◽  
S Karthigai Selvan ◽  
Binu P Thomas ◽  
Ahmedul K Asraff ◽  
Jeby Philip

Ramberg-Osgood (R-O) type stress-strain models are commonly employed during elasto-plastic analysis of metals. Recently, 2-stage and 3-stage R-O variant models have been proposed to replicate stress-strain behavior under large plastic deformation. The complexity of these models increases with the addition of each stage. Moreover, these models have considered deformation till necking only. In this paper, a simplistic multi-stage constitutive model is proposed to capture the strain-hardening non-linearity shown by metals including its post necking behavior. The constitutive parameters of the proposed stress-strain model can be determined using only elastic modulus and yield strength. 3-D digital image correlation was used as an experimental tool for measuring full-field strains on the specimens, which were subsequently utilized to obtain the material parameters. Our constitutive model is demonstrated with an aerospace-grade stainless steel AISI 321 wherein deformation response averaged over the gauge length (GL) and at a local necking zone are compared. The resulting averaged and local material parameters obtained from the proposed model provide interesting insights into the pre and post necking deformation behavior. Our constitutive model would be useful for characterizing highly ductile metals which may or may not depict non-linear strain hardening behavior including their post necking deformations.


2012 ◽  
Vol 472-475 ◽  
pp. 1003-1008 ◽  
Author(s):  
Pei Pei Zhang ◽  
Mei Zhan ◽  
Tao Huang ◽  
He Yang

Spring-back is one of the key factors affecting the forming quality of the NC bending of high-strength TA18 tubes (TA18-HS tubes). Since material parameters have a direct influence on stress and strain fields during the bending and after unloading, the springback of TA18-HS tubes after NC bending depends on material properties to a great degree. In order to study the effect of material parameters, the sensitivity of material parameters on spring-back of TA18-HS tubes is analyzed in this study, using the numerical simulation and the multi-parameters sensitivity analysis method. The results show the following: (1) The springback angle has a positive correlation with the strength coefficient and initial yield stress, and has a negative correlation with the elastic modulus and strain hardening exponent. Besides, with the increase of elastic modulus, the fluctuation of springback goes gently; with the increase of the strength coefficient and initial yield stress, the fluctuation of springback goes abruptly; but with the variation of the strain hardening exponent, the springback fluctuates slightly; (2) The elastic modulus is the most sensitive material parameter on spring-back, the strength coefficient and initial yield stress rank the second and third, respectively, and the strain hardening exponent is the last. The achievement of the study is valuable to eliminate the non-sensitivity parameters, simplify the optimization project, and improve the spring-back prediction capability.


1985 ◽  
Vol 107 (3) ◽  
pp. 228-233 ◽  
Author(s):  
S. T. Clegg ◽  
R. B. Roemer

In cancer hyperthermia treatments, it is important to be able to predict complete tissue temperature fields from sampled temperatures taken at the limited number of locations allowed by clinical constraints. An initial attempt to do this automatically using unconstrained optimization techniques to minimize the differences between experimental temperatures and temperatures predicted from treatment simulations has been previously reported [1]. This paper reports on a comparative study which applies a range of different optimization techniques (relaxation, steepest descent, conjugate gradient, Gauss, Box-Kanemasu, and Modified Box-Kanemasu) to this problem. The results show that the Gauss method converges more rapidly than the others, and that it converges to the correct solution regardless of the initial guess for the unknown blood perfusion vector. A sensitivity study of the error space is also performed, and the relationships between the error space characteristics and the comparative speeds of the optimization techniques are discussed.


2020 ◽  
Author(s):  
Ahmad Sedaghat ◽  
Amir Mosavi

The SIR type models are built by a set of ordinary differential equations (ODE), which are strongly initial value dependant. To fit multiple biological data with SIR type equations requires fitting coefficients of these equations by an initial guess and applying optimization methods. These coefficients are also extremely initial value-dependent. In the vast publication of these types, we hardly see, among simple to highly complicated SIR type methods, that these methods presented more than a maximum of two biological data sets. We propose a novel method that integrates an analytical solution of the infectious population using Weibull distribution function into any SIR type models. The Weibull-SIRD method has easily fitted 4 set of COVID-19 biological data simultaneously. It is demonstrated that the Weibull-SIRD method predictions for susceptible, infected, recovered, and deceased populations from COVID-19 in Kuwait and UAE are superior compared with SIRD original ODE model. The proposed method here opens doors for new deeper studying of biological dynamic systems with realistic biological data trends than providing some complicated, cumbersome mathematical methods with little insight into biological data's real physics.


2008 ◽  
Vol 15 (3-4) ◽  
pp. 257-272 ◽  
Author(s):  
Felipe A.C. Viana ◽  
Valder Steffen Jr. ◽  
Marcelo A.X. Zanini ◽  
Sandro A. Magalhães ◽  
Luiz C.S. Góes

This work deals with the application of a nature-inspired optimization technique to solve an inverse problem represented by the identification of an aircraft landing gear model. The model is described in terms of the landing gear geometry, internal volumes and areas, shock absorber travel, tire type, and gas and oil characteristics of the shock absorber. The solution to this inverse problem can be obtained by using classical gradient-based optimization methods. However, this is a difficult task due to the existence of local minima in the design space and the requirement of an initial guess. These aspects have motivated the authors to explore a nature-inspired approach using a method known as LifeCycle Model. In the present formulation two nature-based methods, namely the Genetic Algorithms and the Particle Swarm Optimization were used. An optimization problem is formulated in which the objective function represents the difference between the measured characteristics of the system and its model counterpart. The polytropic coefficient of the gas and the damping parameter of the shock absorber are assumed as being unknown: they are considered as design variables. As an illustration, experimental drop test data, obtained under zero horizontal speed, were used in the non-linear landing gear model updating of a small aircraft.


Sign in / Sign up

Export Citation Format

Share Document