CN106407544B - A kind of method for building up of IC10 unidirectional solidification material Stiffness prediction models - Google Patents

A kind of method for building up of IC10 unidirectional solidification material Stiffness prediction models Download PDF

Info

Publication number
CN106407544B
CN106407544B CN201610814449.0A CN201610814449A CN106407544B CN 106407544 B CN106407544 B CN 106407544B CN 201610814449 A CN201610814449 A CN 201610814449A CN 106407544 B CN106407544 B CN 106407544B
Authority
CN
China
Prior art keywords
modulus
alloy
young
affected zone
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN201610814449.0A
Other languages
Chinese (zh)
Other versions
CN106407544A (en
Inventor
张宏建
肖健峰
温卫东
崔海涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201610814449.0A priority Critical patent/CN106407544B/en
Publication of CN106407544A publication Critical patent/CN106407544A/en
Application granted granted Critical
Publication of CN106407544B publication Critical patent/CN106407544B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention discloses the method for building up that a kind of IC10 unidirectional solidification material considers the rigidity model that crystal boundary and loading direction influence, comprising steps of 1) monotonic tension of the IC10 monocrystalline at room temperature along [001], [010] and [011] direction is tested;2) IC10 directional solidificating alloy is tested along the monotonic tension of different loading directions at room temperature;3) scanning electron microscope and transmission electron microscopy observation of IC10 directional solidificating alloy;4) based single crystal performance, establishes the Stiffness prediction models of IC10 directional solidificating alloy, and model parameter is calculated by software;5) after the completion of model foundation, the Young's modulus to IC10 directional solidificating alloy along different loading directions is verified.The present invention can accurately predict the Young's modulus that IC10 unidirectional solidification material is loaded along different directions, provide accurate elastic material parameter for the further intensity of material and fatigue study, be of great significance to the engineering design of material.

Description

A kind of method for building up of IC10 unidirectional solidification material Stiffness prediction models
Technical field
The present invention relates to aeronautical material fields, and in particular to turbine disk directionally solidified superalloy, DS superalloy material stiffness model Method for building up.
Background technique
Directionally solidified superalloy, DS superalloy is developed to meet the turbine entrance temperature temperature of aero-engine continuous improvement One of the novel high-temperature alloy come.So-called directional solidification, when being exactly that high temperature alloy melt solidifies in casting mold, by controlling crystal grain The direction of growth, generate the column crystal that is almost parallel to each other.The grain growth direction of directional solidificating alloy and the maximum of material are main Axis direction is parallel, and mechanical property is generally better than the polycrystalline material with general crystal boundary.Currently, directional solidificating alloy elastic constant Acquisition be to be replaced based on test or with the elastic constant approximation of monocrystalline, first method needs a large number of experiments, higher cost, Second method precision of prediction is not high, is in most cases unable to satisfy requirement of engineering precision.
Currently, the prediction simplest method of polycrystalline elastic constant is exactly Voigt the and Reuss method of average, Voigt is based on often answering Become it is assumed that giving the upper limit that polycrystal effective modulus really solves, Reuss is based on iso-stress hypothesis and gives the effective mould of polycrystal The lower limit really solved is measured, though this method predicts that the elastic constant effect of general polycrystal is good, since it does not consider that orientation is solidifying Gu the influence of alloy crystal boundary is less suitable for direct use in prediction directional solidificating alloy.There is researcher to be established with self-consistency theory The prediction technique of directional solidificating alloy elastic constant, but it realizes that process is complex, need to iterate, and engineering application has one It settles finally sex-limited.
Summary of the invention
Goal of the invention: being directed to the above-mentioned prior art, propose a kind of rigidity model method for building up of IC10 unidirectional solidification material, It being capable of Young's modulus under Accurate Prediction directional solidificating alloy difference loading direction.
Technical solution: a kind of rigidity model method for building up of IC10 unidirectional solidification material includes the following steps:
1) monotonic tension test, is carried out along [001], [010] and [011] direction to IC10 single crystal alloy, obtains it respectively Monocrystalline elastic material constant D in three directions11、D12、D44
2) monotonic tension test, is carried out along [001] and [010] direction to IC10 directional solidificating alloy, obtains it along [001] [010] elasticity modulus in direction determines crystal boundary affecting parameters f (T) in rigidity model in conjunction with matlab nonlinear fitting module And n;Wherein, f (T) indicates that crystal boundary is model relevant to a temperature system to the limited degree perpendicular to crystal boundary Direction distortion Number, n are the coefficient for indicating GB affected zone size;
3) it, is seen with the microstructure that scanning electron microscope and transmission electron microscope obtain IC10 directional solidificating alloy Mapping obtains the diameter D and grain boundary width d of each crystal grain in IC10 directional solidificating alloy;
4) elastic property, based on IC10 single crystal alloy considers the influence of loading direction and crystal boundary, and it is solidifying to establish IC10 orientation Gu the rigidity model of alloy, includes the following steps:
4-1), the rigidity model of intra-die is established:
In formula, E1It is the Young's modulus of single intra-die, D11、D12、D44It is normal that monocrystalline elastic material is obtained for step 1) Number, α1、β1、γ1It is coefficient of angularity relevant to loading direction, is defined as follows:
α1=-sin (ψ)
β1=-cos (ψ) sin (θ)
γ1=cos (ψ) cos (θ)
In formula, ψ and θ are Eulerian angles;
4-2), according to the Young's modulus of single intra-die and crystal boundary affecting parameters, the rigidity in GB affected zone is established Model:
In formula, E2It is the Young's modulus in GB affected zone;
4-3), it is assumed that the width d' of GB affected zone is
D'=nd
In formula, n is the coefficient for indicating GB affected zone size, and d is grain boundary width;
According to Voigt and Reuss averaging method, in conjunction with intra-die rigidity model and GB affected zone in rigidity mould Type obtains
In formula, EvoigtAnd EReussRespectively indicate the bound of directional solidificating alloy Young's modulus, fv1Indicate single crystal grain Volume fraction, fv2Indicate the volume fraction of GB affected zone;fv1And fv2It is the function of crystal grain diameter D and grain boundary width d:
fv2=1-fv1
The then Young's modulus E of directional solidificating alloy are as follows:
The utility model has the advantages that showing that crystal boundary has a certain impact to the performance tool of material by experimental study, in monocrystal material not Containing crystal boundary, and contain crystal boundary in unidirectional solidification material, therefore, in the Stiffness prediction models for establishing unidirectional solidification material with regard to necessary Consider the influence of crystal boundary.The performance of monocrystal material and unidirectional solidification material is anisotropy, and experimental study shows different loads The rigidity of monocrystalline is different under direction, and the rigidity of monocrystal material and unidirectional solidification material is also different, so to establish IC10 Rigidity model of the directional solidificating alloy along different loading directions, it is necessary to fully consider the influence of loading direction and crystal boundary to rigidity. The method of the present invention quantitatively considers crystal boundary in deformation process on the basis of original monocrystalline rigidity model and particle near its is become The restriction effect of shape, the actual response flexible deformation behavior of directional solidificating alloy, rigidity model prediction effect of the invention compared with It is good, it can accurately predict the Young's modulus that IC10 unidirectional solidification material is loaded along different directions, be the further intensity of material Accurate elastic material parameter is provided with fatigue study, is of great significance to the further engineering design of material.
Detailed description of the invention
Fig. 1 is the flow chart of the method for the present invention;
Fig. 2 is Young's modulus test result and prediction result figure of the IC10 directional solidificating alloy along different loading directions.
Specific embodiment
Further explanation is done to the present invention with reference to the accompanying drawing.
As shown in Figure 1, a kind of rigidity model method for building up of IC10 unidirectional solidification material, includes the following steps:
1) it, in order to obtain the elastic constant of IC10 single crystal alloy at room temperature, draws materials from IC10 single crystal alloy masterbatch, edge [001], the standard tensile specimen that φ 5mm is processed into [010] and [011] three direction carries out the test of static(al) monotonic tension, respectively It is obtained in the monocrystalline elastic material constant D in three directions11、D12、D44, experimental condition is shown in Table 1.
2) monotonic tension test, is carried out along [001] and [010] direction at room temperature to IC10 directional solidificating alloy, is obtained It determines that crystal boundary influences in rigidity model in conjunction with matlab nonlinear fitting module along the elasticity modulus in [001] and [010] direction Parameter f (T) and n;Wherein, f (T) indicates that crystal boundary is one related to temperature to the limited degree perpendicular to crystal boundary Direction distortion Model coefficient, n is the coefficient for indicating GB affected zone size, and n characterizes the size of GB affected zone.And edge [025], [011] and monotonic tension test is carried out in [052] direction, and the elasticity modulus that these directions obtain establishes model to the present invention Model verifying.Experimental condition is shown in Table 1.
Table 1
Experimental condition IC10 monocrystalline tension test IC10 directional solidificating alloy tension test
Specimen size φ5mm φ5mm
Strain rate 10-3/s 10-3/s
Loading direction [001]、[010]、[011] [001]、[025]、[011]、[052]、[010]
Temperature Room temperature Room temperature
Range of strain It stretches until sample fracture It stretches until sample fracture
Testing equipment SDS-50 electro-hydraulic servo static and dynamic test machine SDS-50 electro-hydraulic servo static and dynamic test machine
3) sheet metal of 1.5mm thickness, is cut out on the sample of not test (N.T.) IC10 directional solidificating alloy respectively, respectively according to sweeping The production specification for retouching Electronic Speculum and transmission electron microscope sample processes scanning electron microscope and transmission electron microscope sample for observation, with sweeping It retouches electron microscope and transmission electron microscope obtains the microstructure observation chart of IC10 directional solidificating alloy, obtain IC10 orientation The diameter D and grain boundary width d of each crystal grain in solidified superalloy.
4) elastic property, based on IC10 single crystal alloy considers the influence of loading direction and crystal boundary, and it is solidifying to establish IC10 orientation Gu the rigidity model of alloy, includes the following steps:
4-1), the rigidity model of intra-die is established:
In formula, E1It is the Young's modulus of single intra-die, D11、D12、D44It is normal that monocrystalline elastic material is obtained for step 1) Number, α1、β1、γ1It is coefficient of angularity relevant to loading direction, is defined as follows:
α1=-sin (ψ)
β1=-cos (ψ) sin (θ)
γ1=cos (ψ) cos (θ)
In formula, ψ and θ are Eulerian angles.
4-2), according to the Young's modulus of single intra-die and crystal boundary affecting parameters, the rigidity in GB affected zone is established Model:
In formula, E2It is the Young's modulus in GB affected zone;F (T) is obtained by step 2).
4-3), it is assumed that the width d' of GB affected zone is
D'=nd
In formula, n is the coefficient for indicating GB affected zone size, and d is grain boundary width;
According to Voigt and Reuss averaging method, in conjunction with intra-die rigidity model and GB affected zone in rigidity mould Type obtains
In formula, EvoigtAnd EReussRespectively indicate the bound of directional solidificating alloy Young's modulus, fv1Indicate single crystal grain Volume fraction, fv2Indicate the volume fraction of GB affected zone;fv1And fv2It is the function of crystal grain diameter D and grain boundary width d:
fv2=1-fv1
The then Young's modulus of directional solidificating alloyAre as follows:
In the present embodiment, each parameter is as shown in table 2 in model:
Table 2
Model parameter D11/GPa D12/GPa D44/GPa D/mm d/mm f n
Numerical value 290.92 188.76 133.62 400 4.4 0.35 7
The present invention quantitatively considers in deformation process crystal boundary for matter near its on the basis of original monocrystalline rigidity model The restriction effect of point deformation, the actual response flexible deformation behavior of directional solidificating alloy, so the new rigidity model of exploitation Prediction effect is preferable.After the completion of the rigidity model to IC10 directional solidificating alloy is established, loaded with IC10 along different directions The Young's model that monotonic tension test measures is verified, and the Young's modulus of prediction is compared with test result, sees Fig. 2, It was found that prediction result and test result are coincide well, the reliability of model is demonstrated.
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered It is considered as protection scope of the present invention.

Claims (1)

1. a kind of rigidity model method for building up of IC10 unidirectional solidification material, which comprises the steps of:
1) monotonic tension test, is carried out along [001], [010] and [011] direction to IC10 single crystal alloy, obtains it respectively at this The monocrystalline elastic material constant D in three directions11、D12、D44
2) monotonic tension test, is carried out along [001] and [010] direction to IC10 directional solidificating alloy, obtain its edge [001] and [010] elasticity modulus in direction, in conjunction with matlab nonlinear fitting module determine in rigidity model crystal boundary affecting parameters f (T) and n;Wherein, f (T) indicates that crystal boundary is a model coefficient relevant to temperature to the limited degree perpendicular to crystal boundary Direction distortion, N is the coefficient for indicating GB affected zone size;
3) it, is observed with the microstructure that scanning electron microscope and transmission electron microscope obtain IC10 directional solidificating alloy Figure, obtains the diameter D and grain boundary width d of each crystal grain in IC10 directional solidificating alloy;
4) elastic property, based on IC10 single crystal alloy considers the influence of loading direction and crystal boundary, establishes the conjunction of IC10 directional solidification The rigidity model of gold, includes the following steps:
4-1), the rigidity model of intra-die is established:
In formula, E1It is the Young's modulus of single intra-die, D11、D12、D44Monocrystalline elastic material constant, α are obtained for step 1)1、 β1、γ1It is coefficient of angularity relevant to loading direction, is defined as follows:
α1=-sin (ψ)
β1=-cos (ψ) sin (θ)
γ1=cos (ψ) cos (θ)
In formula, ψ and θ are Eulerian angles;
4-2), according to the Young's modulus of single intra-die and crystal boundary affecting parameters, the rigidity model in GB affected zone is established:
In formula, E2It is the Young's modulus in GB affected zone;
4-3), it is assumed that the width d' of GB affected zone is
D'=nd
In formula, n is the coefficient for indicating GB affected zone size, and d is grain boundary width;
According to Voigt and Reuss averaging method, in conjunction with intra-die rigidity model and GB affected zone in rigidity model obtain It arrives
In formula, EvoigtAnd EReussRespectively indicate the bound of directional solidificating alloy Young's modulus, fv1Indicate the volume of single crystal grain Score, fv2Indicate the volume fraction of GB affected zone;fv1And fv2It is the function of crystal grain diameter D and grain boundary width d:
fv2=1-fv1
The then Young's modulus of directional solidificating alloyAre as follows:
CN201610814449.0A 2016-09-12 2016-09-12 A kind of method for building up of IC10 unidirectional solidification material Stiffness prediction models Expired - Fee Related CN106407544B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610814449.0A CN106407544B (en) 2016-09-12 2016-09-12 A kind of method for building up of IC10 unidirectional solidification material Stiffness prediction models

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610814449.0A CN106407544B (en) 2016-09-12 2016-09-12 A kind of method for building up of IC10 unidirectional solidification material Stiffness prediction models

Publications (2)

Publication Number Publication Date
CN106407544A CN106407544A (en) 2017-02-15
CN106407544B true CN106407544B (en) 2019-05-14

Family

ID=57999279

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610814449.0A Expired - Fee Related CN106407544B (en) 2016-09-12 2016-09-12 A kind of method for building up of IC10 unidirectional solidification material Stiffness prediction models

Country Status (1)

Country Link
CN (1) CN106407544B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111103316B (en) * 2018-10-25 2021-05-25 中国科学院上海硅酸盐研究所 Calculation method for non-conductive ceramic material non-charge balance voltage
CN114807635B (en) * 2022-05-06 2024-03-05 云南锡业集团(控股)有限责任公司研发中心 Method for preparing high-purity indium by multi-channel array directional solidification based on machine learning

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103406666A (en) * 2013-06-16 2013-11-27 北京工业大学 IC10 alloy connecting and repairing method for controlling dendrite growth direction through laser

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103406666A (en) * 2013-06-16 2013-11-27 北京工业大学 IC10 alloy connecting and repairing method for controlling dendrite growth direction through laser

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
IC10合金的力学性能试验及本构模型研究;张宏建;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20110615(第6期);第C028-1页
Recrystallization behaviors of alloy IC10 at elevated temperature:experiments and modeling;Hongjian Zhang 等;《Journal of Materials Science》;20110909;第46卷(第4期);第1076-1082页
Temperature dependence of faceting in Σ5(310)[001] grain boundary of SrTiO3;Sung Bo Lee 等;《Acta Materialia》;20030225;第51卷(第4期);第975-981页
双向应力状态下IC10高温合金的屈服行为研究;陈雷 等;《航空学报》;20120125;第33卷(第1期);第77-84页

Also Published As

Publication number Publication date
CN106407544A (en) 2017-02-15

Similar Documents

Publication Publication Date Title
JP6828476B2 (en) Edge fracture prediction method, program and recording medium
CN106407544B (en) A kind of method for building up of IC10 unidirectional solidification material Stiffness prediction models
CN112036057B (en) Hole position control method for hole making of aircraft structural part
JP5403256B2 (en) Apparatus and method for evaluating fatigue strength of polycrystalline metal materials
Wan et al. Improved inverse filter for the correction of distorted measured cutting forces
CN106096163A (en) A kind of method for building up of carbon fibre composite cutting Model
CN105067406A (en) Argon ion polishing/scanning electron microscope analyzing and sampling method for shale samples
Litak et al. Cutting force response in milling of Inconel: Analysis by wavelet and Hilbert-Huang transforms
CN103049585B (en) Automobile Chassis Parts and vehicle body be load extracting method and device under limiting condition
CN109635365A (en) A kind of process of control casing cutting parts deformation
CN110096776B (en) Parametric modeling and orientation optimization method for special-shaped air film holes on single-crystal turbine blade
CN107326314B (en) A method of δ phase resolving in dynamic state volume fraction in prediction nickel-base alloy containing niobium
CN106777773B (en) Diamond circular saw blade design and sawing process parameter making method
CN105510523A (en) Failure analysis method for triode transistor
CN106442407B (en) Reliability verification method for additive manufacturing product
CN116205112A (en) Grain boundary stress concentration factor calculation method, apparatus, medium, and program
CN110465584A (en) A kind of progressive molding limit test method of mesh plate
CN105510550B (en) A method of it is promoted and increases material titanium alloy material fatigue behaviour
CN113182565B (en) Weak-rigidity molded surface regional milling method
CN115544687A (en) Method for predicting cutting performance of precise complex cutter
CN111209657B (en) Solid deformation interface calculation method considering liquid surface tension
CN109117547B (en) Method for optimizing dynamic characteristics of large part structure of machine tool
Peng et al. Analysis of micro/mesoscale sheet forming process by strain gradient plasticity and its characterization of tool feature size effects
Schäfer et al. Stage-I fatigue crack studies in order to validate the dislocation-free zone model of fracture for bulk materials
Sun et al. Study on the effect of diamond wire saw-cutting process parameters on brittle cracks in monocrystalline silicon

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190514