CASTING SIMULATION METHOD
20180001380 · 2018-01-04
Assignee
Inventors
- Makoto YOSHIDA (Shinjuku-ku, Tokyo, JP)
- Yuichi MOTOYAMA (Tsukuba-shi, Ibaraki, JP)
- Toshimitsu OKANE (Tsukuba-shi, Ibaraki, JP)
- Yoya FUKUDA (Sunto-gun, Shizuoka, JP)
Cpc classification
B22D46/00
PERFORMING OPERATIONS; TRANSPORTING
C21D11/00
CHEMISTRY; METALLURGY
International classification
Abstract
Provided is a casting simulation method capable of expressing influence of different inelastic strains produced at different temperatures on strain hardenability at room temperature. The following amount of effective equivalent inelastic strain ε.sub.effective inelastic is substituted into a constitutive equation in which an amount of equivalent inelastic strain is used as a degree of work hardening:
an amount of effective equivalent inelastic strain ε.sub.effective inelastic=∫.sub.o.sup.t{h.sub.(T)/h.sub.(RT)}{(Δε.sub.inelastic/Δt)}dt
, where T denotes a temperature with inelastic strain, h.sub.(T) denotes an increment of yield strength at room temperature with respect to an amount of inelastic strain at the temperature with inelastic strain, h.sub.(RT) denotes an increment of yield strength at room temperature with respect to an amount of inelastic strain applied at room temperature, h.sub.(T)/h.sub.(RT) denotes an effective inelastic strain coefficient α(T), Δε.sub.inelastic/Δt denotes an equivalent inelastic strain rate, and t denotes a time from 0 second in analysis.
Claims
1. A casting simulation method capable of expressing influence of different inelastic strains produced at different temperatures on work hardening, namely on increase in yield stress, at room temperature, the influence varying with differences in recovery at the different temperatures, by introducing an amount of effective equivalent inelastic strain to an elasto-plastic constitutive equation and/or an elasto-viscoplastic constitutive equation in which an amount of equivalent inelastic strain is used as a degree of work hardening, namely an amount of increase in yield stress, such as an elasto-plastic constitutive equation in which a yield function is expressed as f=f(σ.sub.eff,ε.sub.eff.sup.p,T) or an elasto-viscoplastic constitutive equation in which a relation between equivalent stress, equivalent viscoplastic strain, viscoplastic strain rate, and temperature is expressed as σ.sub.eff.=F(ε.sub.eff..sup.vp,{dot over (ε)}.sub.eff..sup.vp,T), wherein an amount of effective equivalent inelastic strain ε.sub.effective inelastic obtained by Eq. (1) below is used:
the amount of effective equivalent inelastic strain ε.sub.effective inelastic=∫.sub.o.sup.t{h.sub.(T)/h.sub.(RT)}{(Δε.sub.inelastic/Δt)}dt (1) , where T denotes a temperature with inelastic strain, h.sub.(T) denotes an increment of yield strength at room temperature with respect to an amount of inelastic strain at the temperature with inelastic strain, h.sub.(RT) denotes an increment of yield strength at room temperature with respect to an amount of inelastic strain applied at room temperature, h.sub.(T)/h.sub.(RT) denotes an effective inelastic strain coefficient α(T), Δε.sub.inelastic/Δt denotes an equivalent inelastic strain rate, and t denotes a time from 0 second in analysis.
2. The casting simulation method according to claim 1, the effective inelastic strain coefficient α(T) is obtained by: applying different inelastic pre-strains to a test piece at different temperatures; cooling the test piece to room temperature; performing a tensile test or a compression test on the test piece at room temperature; and measuring influence of amounts of the inelastic pre-strains applied at the different temperatures on the increase in yield stress.
3. The casting simulation method according to claim 1, wherein a stress-equivalent inelastic strain curve is transformed into a stress-effective equivalent inelastic strain curve using α(T), and based on the stress-effective equivalent inelastic strain curve, a determination is made of a material constant in a constitutive equation to which the amount of effective equivalent inelastic strain ε.sub.effective inelastic is introduced.
4. The casting simulation method according to claim 1, wherein when α(T) is substituted into Eq. (1) in a temperature range in which α(T) is 0 or a negative value and a stress-equivalent inelastic strain curve at the temperature T indicates work hardening, α(T) is corrected from 0 or the negative value to a small positive value.
5. The casting simulation method according to claim 2, wherein a stress-equivalent inelastic strain curve is transformed into a stress-effective equivalent inelastic strain curve using α(T), and based on the stress-effective equivalent inelastic strain curve, a determination is made of a material constant in a constitutive equation to which the amount of effective equivalent inelastic strain ε.sub.effective inelastic is introduced.
6. The casting simulation method according to claim 2, wherein when α(T) is substituted into Eq. (1) in a temperature range in which α(T) is 0 or a negative value and a stress-equivalent inelastic strain curve at the temperature T indicates work hardening, α(T) is corrected from 0 or the negative value to a small positive value.
7. The casting simulation method according to claim 3, wherein when α(T) is substituted into Eq. (1) in a temperature range in which α(T) is 0 or a negative value and a stress-equivalent inelastic strain curve at the temperature T indicates work hardening, α(T) is corrected from 0 or the negative value to a small positive value.
8. The casting simulation method according to claim 5, wherein when α(T) is substituted into Eq. (1) in a temperature range in which α(T) is 0 or a negative value and a stress-equivalent inelastic strain curve at the temperature T indicates work hardening, α(T) is corrected from 0 or the negative value to a small positive value.
Description
BRIEF DESCRIPTION OF THE DRAWING
[0034] In the accompanying drawings:
[0035]
[0036]
[0037]
[0038]
[0039]
[0040]
[0041]
[0042]
[0043]
[0044]
[0045]
[0046]
[0047]
[0048]
[0049]
[0050]
DETAILED DESCRIPTION
[0051] The following describes the present disclosure in detail.
[0052] As used herein, an equivalent inelastic strain rate calculated by thermal stress analysis is multiplied by an “effective inelastic strain coefficient α(T) indicative of temperature dependency”, which is preferably experimentally determined, and the result is used as an effective inelastic strain rate. The effective inelastic strain rate is then integrated over a time from 0 second in analysis to determine the amount of effective equivalent inelastic strain, which in turn is, in place of the amount of equivalent inelastic strain, used as a measure of hardening in a constitutive equation.
[0053] The amount of effective equivalent inelastic strain is applicable to any constitutive equation, whether an elasto-plastic constitutive equation or an elasto-viscoplastic constitutive equation, as long as it is a constitutive equation using the amount of equivalent inelastic strain conventionally as a measure of work hardening as described in paragraph 0008.
[0054] Specific examples are an elasto-plastic constitutive equation and/or an elasto-viscoplastic constitutive equation in which the amount of equivalent inelastic strain is used as a degree of work hardening (namely an amount of increase in yield stress) such as an elasto-plastic constitutive equation in which a yield function is expressed as f=f(σ.sub.eff,ε.sub.eff.sup.p,T) or an elasto-viscoplastic constitutive equation in which a relation between equivalent stress, equivalent viscoplastic strain, viscoplastic strain rate, and temperature is expressed as σ.sub.eff.=F(ε.sub.eff..sup.vp,{dot over (ε)}.sub.eff..sup.vp,T).
[0055] The determination of the effective inelastic strain coefficient α(T) is carried out by applying inelastic pre-strains of different magnitudes at the corresponding temperatures at which α(T) is to be obtained, then cooling to room temperature, conducting a tensile test or a compression test, and measuring the increase in yield stress.
[0056] The following describes how to determine α(T) in detail.
[0057] A tensile test piece is heated in a way as presented in the temperature history in
[0058] Desirably, the same temperature history is set for all test conditions. The reason for this is to eliminate the influence of the temperature history of the test piece on the measured values.
[0059] In the context of the present disclosure, the above-described test is not limited to a particular test, and may be a tensile test or a compression test as long as it can provide a stress-strain curve and enables measurement of yield stress. For the tensile test in this disclosure, for example, a publicly-known and widely-used tensile test may be used, such as JIS Z 2241:2011. For the compression test, for example, a publicly-known and widely-used compression test may be used, such as JIS K 7181:2011.
[0060] The results obtained in the above test are conceptually illustrated in
[0061] In a temperature range in which an effective inelastic strain coefficient α(T) is experimentally determined to be 0 or a negative value, i.e., in which inelastic strain applied at the temperature T should not contribute to work hardening at room temperature, if the stress-equivalent inelastic strain curve at the temperature T indicates work hardening and if the effective inelastic strain coefficient α(T) is 0, then a constitutive equation to which the amount of effective equivalent inelastic strain is introduced involves no effective inelastic strain, and thus is not able to express work hardening in principle. In other words, the stress-strain curve displays elasto-perfectly plastic behavior or elasto-perfectly viscoplastic behavior. In that case, the reproducibility of the stress-strain curve deteriorates, resulting in lower accuracy of predictions on thermal stress and deformation.
[0062] An effective inelastic strain is produced as long as the effective inelastic strain coefficient α(T) is not 0, making it possible to express work hardening in the stress-inelastic strain curve at the temperature T. Accordingly, even in a temperature range with the inelastic strain coefficient α(T) being 0 or a negative value, if a positive small value, rather than 0 or a negative value, is corrected appropriately for use as an effective inelastic strain coefficient α(T), it is possible to express work hardening in a stress-equivalent inelastic strain curve with a slight ineffective inelastic strain. The value of α(T) at the time of correction is in a range of 0<α<0.5, and desirably 0<α<0.1, although it depends on the alloy type. As an example, α(T) is corrected by the following linear interpolation using an effective inelastic strain coefficient at maximum temperature α.sub.min in a temperature range in which α(T) is experimentally determined to be non-zero and a maximum temperature T.sub.max (which may alternatively be a solidus temperature) in a temperature range in which α is experimentally determined to be 0:
where
[0063] α.sub.extrapolation(T) denotes a value of a as corrected in a temperature range in which α is experimentally determined to be 0;
[0064] T.sub.max denotes the maximum temperature (which may alternatively be a solidus temperature) in a temperature range where a is experimentally determined to be 0,
[0065] T.sub.min denotes a maximum temperature in a temperature range in which α is experimentally determined to be non-zero;
[0066] α.sub.min denotes a value of α in T.sub.min; and
[0067] T denotes a temperature above T.sub.min
[0068] As an example, material constants in a constitutive equation to which the amount of effective inelastic strain is introduced are determined as explained below in the case of introducing the amount of effective inelastic strain to the constants (K(T), n(T), and m(T)) of the extended Ludwik's law.
[0069] The extended Ludwik's law is as follows. ε.sub.0 is a constant necessary for calculation and usually a small value of 1×10.sup.−6.
σ=K(T)(ε.sub.inelastic+ε.sub.0).sup.n(T)({dot over (ε)}.sub.inelastic).sup.m(T)
[0070] When introducing the amount of effective inelastic strain, it is expressed as:
σ=K(T)(ε.sub.effective inelastic+ε.sub.0).sup.n(T)({dot over (ε)}.sub.inelastic).sup.m(T)
[0071] The term with an index n.sub.(T) representing the degree of work hardening includes the amount of effective equivalent inelastic strain as a variable. In addition, the term with an index m.sub.(T) representing the strain rate dependence of the stress-strain curve includes an equivalent inelastic strain rate as a variable. Since the equation as a whole includes the amount of effective equivalent inelastic strain, for each temperature, by substituting the inelastic strain rate into the term with m.sub.(T), K.sub.(T), n.sub.(T), and m.sub.(T) are determined by numerical optimization to fit the stress-effective equivalent inelastic strain curve.
[0072] Following the above procedure, the amount of effective equivalent inelastic strain ε.sub.effective inelastic is determined. In the disclosure, the amount of effective equivalent inelastic strain is used to simulate the influence of inelastic strain applied at the temperature in question on work hardening at room temperature. Although details of the procedures for casting simulation will be described in the Examples, the points are summarized as follows.
[0073] In the simulation according to the disclosure, it is possible to adopt a constitutional expression that uses the amount of equivalent inelastic strain as a measure of work hardening conventionally used for casting simulation. An exemplary elasto-plastic constitutive equation is:
ε.sub.ij=ε.sup.e.sub.ij+ε.sup.p.sub.ij
σ.sub.ij=D.sub.ijkl(T)ε.sup.e.sub.kl
f=f(σ.sub.eff,ε.sub.eff.sup.p,T)
, where T is the temperature, σ.sub.ij is the stress, σ.sub.eff. is the equivalent stress, ε.sub.ij is the total strain, ε.sup.e.sub.ij is the elastic strain, ε.sup.p.sub.ij is the plastic strain, f is the yield function, and D.sub.ijkl is the fourth-order constitutive tensor.
[0074] Alternatively, an exemplary elasto-viscoplastic constitutive equation is:
ε.sub.ij=ε.sup.e.sub.ij+ε.sup.vp.sub.ij
σ.sub.ij=D.sub.ijkl(T)ε.sup.e.sub.kl
σ.sub.eff.=F(ε.sub.eff..sup.vp,{dot over (ε)}.sub.eff..sup.vp,T)
, where σ.sub.eff. is the equivalent stress, ε.sub.eff..sup.vp is the equivalent viscoplastic strain, {dot over (ε)}.sub.eff..sup.vp is the equivalent viscoplastic strain rate.
[0075] To these constitutive equations, the amount of effective equivalent inelastic strain ε.sub.effective inelastic defined by Eq. (1) according to the present disclosure, instead of the amount of equivalent inelastic strain conventionally used, may be introduced or substituted.
[0076] The well-known and widely-used procedures for casting simulation are:
[0077] (I) element creation step;
[0078] (II) element definition step;
[0079] (III) heat transfer analysis step;
[0080] (IV) thermal stress analysis step; and
[0081] (V) analysis result evaluation step.
[0082] In the present disclosure, material constants in a constitutive equation are determined in step (II) using an equivalent stress-effective equivalent inelastic strain curve, and are input to a constitutive equation to which the amount of effective inelastic strain is introduced. Then, in the thermal stress analysis step (IV), the amount of effective equivalent inelastic strain is calculated, and the result, instead of the amount of equivalent inelastic strain conventionally used, is used as a parameter representing the amount of work hardening to calculate thermal stress.
EXAMPLES
[0083] The following describes how the present disclosure enables prediction with high accuracy of the influence of amounts of inelastic pre-strains produced at different temperatures on the work hardening behavior at room temperature in a typical aluminum die-casting alloy, JIS ADC12, using an extended Ludwik equation, which is a typical elasto-viscoplastic constitutive equation.
[0084] A specific form of the equation before and after the introduction of the amount of effective inelastic strain is as presented in paragraph 0032.
[0085] Firstly, a typical aluminum die-casting alloy, JIS ADC12, is analyzed for an effective inelastic strain coefficient α(T) and material constants (K(T), n(T), and m(T)) at each temperature according to the procedures described in paragraphs 0026 to 0033.
[0086] Firstly, tensile tests were performed to obtain stress-equivalent inelastic strain curves required to determine material constants K(T), n(T), and m(T). Stress-inelastic strain curves were obtained under a set of conditions including: experimental strain rates of 10.sup.−3/s and 10.sup.−4/s and test temperatures of RT, 200° C., 250° C., 300° C., 350° C., 400° C., and 450° C. Each test pieces was obtained by casting JIS ADC12 in a copper mold and processing it into the shape of a tensile test piece.
[0087] In these tests, all the test pieces were heated from room temperature to 450° C., then subjected to heat treatment at 450° C. for 1 hour to cause precipitates to be re-dissolved, and cooled to the test temperature as soon as possible so that the mechanical properties at the time of cooling can be examined accurately. As soon as the test temperature was reached, the tensile test was carried out.
[0088] In addition, tests for determining an effective inelastic strain coefficient α(T) were carried out at RT, 200° C., 250° C., 300° C., 350° C., 400° C., and 450° C. After solution treatment at 450° C. for 1 hour, the temperature was lowered to a temperature at which the target inelastic pre-strain was to be applied to the test piece following the temperature history presented in
[0089] After application of pre-strains, each test piece was cooled to room temperature and quenched with dry ice to eliminate the effect of the increase in yield stress caused by natural aging. Then, the 0.2% offset yield stress was determined by conducting a tensile test on each test piece at room temperature. The results are presented in
[0090] By using the effective inelastic strain coefficient α.sub.(T) thus obtained, the stress-equivalent inelastic strain curve obtained in paragraph 0040 was transformed into a stress-effective equivalent inelastic strain curve, and the values of K.sub.(T), m.sub.(T), and n.sub.(T) were obtained as described in paragraphs 0032 and 0033. The values of K.sub.(T), m.sub.(T), and n.sub.(T) are presented in
[0091]
[0092] It is conceivable, however, that if α is corrected as described in paragraph 0031 with the temperature of α=0, T.sub.max in paragraph 0031, being temporarily set as a liquidus temperature, the reproducibility of the stress-inelastic strain curve also improves.
[0093] For extended Ludwik equations incorporating or not incorporating (in the case of a conventional example) the amount of effective equivalent inelastic strain according to the present disclosure, calculation was made to determine the effect of inelastic pre-strains applied at different temperatures on the yield stress at room temperature, and the calculation results were compared as presented in
[0094] In the conventional extended Ludwik equation not incorporating the amount of effective equivalent inelastic strain, in principle, inelastic strains produced in different temperature ranges are all considered as equivalent to one another and included as a measure of hardening. Accordingly, as is clear from
[0095] On the other hand, as can be seen from
[0096] To examine the influence of inelastic pre-strain at different temperatures on the yield stress at room temperature, a comparison was made between experimental values and calculated values according to the extended Ludwik equation to which the amount of effective equivalent inelastic strain is introduced, and the results are presented in
[0097] It can be seen from the figure that the analysis program incorporating the amount of effective equivalent inelastic strain could reproduce the behavior at 300° C. or higher at which an increase in yield stress is independent of the amount of inelastic pre-strain. In addition, this program could accurately reproduce the behavior even at 300° C. or lower at which an increase in yield stress depends on the amount of inelastic pre-strain.
[0098]
[0099] It can be seen from the figure that pre-strain applied at 700° C. does not contribute to work hardening at room temperature. In contrast, inelastic strain applied at 350° C. contributes to work hardening at room temperature and the amount of work hardening is proportional to the amount of inelastic strain applied. This behavior is identical to that observed in ADC12, and from this follows that the present disclosure is also applicable to FCD400.