CN107491627A - Composite drilling process finite element Global sensitivity analysis method - Google Patents
Composite drilling process finite element Global sensitivity analysis method Download PDFInfo
- Publication number
- CN107491627A CN107491627A CN201710961870.9A CN201710961870A CN107491627A CN 107491627 A CN107491627 A CN 107491627A CN 201710961870 A CN201710961870 A CN 201710961870A CN 107491627 A CN107491627 A CN 107491627A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- composite
- sample
- finite element
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 98
- 239000002131 composite material Substances 0.000 title claims abstract description 55
- 238000005553 drilling Methods 0.000 title claims abstract description 46
- 238000010206 sensitivity analysis Methods 0.000 title claims abstract description 31
- 238000004458 analytical method Methods 0.000 claims abstract description 14
- 238000005520 cutting process Methods 0.000 claims abstract description 12
- 239000000463 material Substances 0.000 claims abstract description 11
- 230000035945 sensitivity Effects 0.000 claims abstract description 10
- 238000009826 distribution Methods 0.000 claims abstract description 5
- 230000001419 dependent effect Effects 0.000 claims abstract description 4
- 238000007476 Maximum Likelihood Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 claims description 2
- 239000004744 fabric Substances 0.000 claims 1
- 238000004364 calculation method Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 3
- 229910052799 carbon Inorganic materials 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 3
- 239000000835 fiber Substances 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000009833 condensation Methods 0.000 description 2
- 230000005494 condensation Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005304 joining Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
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)
- Complex Calculations (AREA)
Abstract
A kind of composite drilling process finite element Global sensitivity analysis method, it is characterized in that it comprises the following steps：First, multigroup random number is produced according to the distribution of the relevant parameter of composite to input as sample；Secondly, drilling process is emulated by finite element software and produces corresponding cutting force as sample output by the use of least squares method；3rd, according to " state-contingent parametric method "（States Dependent Parameter,SDP）Conditional expectation of the parameters for output is sought respectively；4th, according to Global sensitivity analysis method --- " sobol methods " analyzes level of sensitivity of the different parameters for multiple material drill thrust.The present invention gives composite drilling process finite element Global sensitivity analysis method, embodies composite each parameter in drilling and, for the influence degree of drill thrust, is advantageous to, to Simplified analysis process, save the time, improves efficiency.
Description
Technical field
The present invention relates to a kind of composite technology, especially a kind of composite is subject to technology, specifically a kind of
Composite drilling process finite element Global sensitivity analysis method, belong to the finite element Sensitivity Analysis Method of composite,
The particularly composite Global sensitivity analysis method of material parameter to drill thrust in drilling process.
Background technology
Composite laminated structures have been obtained extensively with high specific strength, specific stiffness and low quality in aerospace field
Application.The design of many empty day structures, using the special Buckling of composite and its stronger designability, with up to
To the purpose for improving stability and reduction architecture quality.At present, Joining Technology main in composite assembling is that machinery connects
Connect, this just needs the drilling processing that pre-manufactured hole is carried out to it.But because composite materials property has larger dispersiveness
And randomness, it is therefore necessary to consider these uncertain influences to composite suffered drill thrust in drilling, i.e. composite wood
Expect level of sensitivity of each parameter for drill thrust, and then Simplified analysis process, improve efficiency, shorten analysis time.
The sensitivity analysis of current research is only to export response directly to seek input variable local derviation, then tries to achieve and is inputting
Derivative value at mean variable value point.The sensitivity analysis only considers that the uncertain of input variable responds not to output from part
Deterministic influence, therefore referred to as Local sensitivity analysis, due to the dispersiveness and randomness of Compound Material Engineering characteristic, the method
In analysis with certain limitation.Although having had some Global sensitivity analysis methods at present, such as Saltellid and
The nonparametric technique of the propositions such as Hehon, the square independent analysis that Borgonovo is proposed, and Sobol propose based on variance
Analysis method, but due to being started late to the research of composite drilling, composite is ground in drilling process sensitivity analysis
Study carefully deficiency, Global sensitivity analysis method is when composite drilling is analyzed using seldom.
The content of the invention
The purpose of the present invention is to be directed to influence of the composite parameters to composite suffered drill thrust in drilling
The needs of analysis, there is dispersed and randomness answer to overcome existing Local sensitivity analysis method not to be suitable for mechanical characteristic
The deficiency of condensation material, existing Global sensitivity analysis method is used for reference, invent a kind of global spirit of composite drilling process finite element
Basis of sensitivity analysis method, it considers composite each parameter in drilling and, for the influence degree of drill thrust, given based on side
The composite drilling process finite element Global sensitivity analysis method of difference.
The technical scheme is that：
A kind of composite drilling process finite element Global sensitivity analysis method, its feature are as follows：
First, multigroup random number is produced according to the distribution of composite relevant parameter to input as sample；
Secondly, the three-dimensional Drilling model of composite is established using finite element software；
3rd, above-mentioned three-dimensional Drilling model is called, emulating drilling process by finite element software produces drill thrust-time
Figure, obtain corresponding cutting force by the use of least squares method and exported as sample；
4th, each ginseng is asked according to " state-contingent parametric method " (States Dependent Parameter, SDP) respectively
Conditional expectation of the number for output；
5th, according to the Global sensitivity analysis method based on variance --- " sobol methods " analyzes different parameters for multiple
The size of the influence degree of material drill thrust.
The relevant parameter of described composite is all to obey the stochastic variable being just distributed very much, and the coefficient of variation is 0.06,
And it is uncorrelated, the random sample X that sample size is N is produced according to the distributed constant of each variable using Matlab softwarest(x1,x2,
x3…xi) (t=1,2,3 ... N) inputted as sample.
Described finite element software is ABAQUS/Explicit.
Three-dimensional Drilling model n times are called, corresponding cutting force Y is produced by Abaqus software emulations drilling processt(t=
1,2,3 ... N) exported as sample.
Parameters X is asked according to SDP methods respectively1, X2, X3,…XiFor sample output Y conditional expectation E (Y | X1), E (Y |
X2) ... E (Y | Xi) step is as follows：
According to control theory, establishing state space SDP models is：
observation equation:Yt=Pi,t+ei,t；
SDP model:Pi,t=Pi,t-1+di,t-1；
di,t=di,t-1+ηi,t, ηi,t~N (0, σ2)；
In formula：Pi,tFor with state variable Xi,tThe SDP state relevant parameters of change, it is basic variable Xi,tFunction, et
(observation interference), ηi,t(system disturbance) be average be zero, variance σ2And ση 2White Gaussian noise input, its value can use
Maximum likelihood optimization is estimated to obtain.
Then recursive Kalman filtering (kalam Filtreing, KF) is used in circulation backward fit procedure and is accordingly passed
Smooth (Fixed Interval Smoothing, the FIS) rule of fixed interval is returned to estimate to obtain each state relevant parameter Pi,t。
Finally according to formula：
E(Yt|Xi,t)=Yt-Pi,t；
Calculate each parameter for output conditional expectation E (Y | Xi,t)；
According to Global sensitivity analysis method --- " sobol methods " calculates each parameter XiGlobal sensitivity to exporting Y
Process is as follows：
According to formula:
Calculate foregoing gained sample output YtVariance V (Y) and each parameter XiFor Y conditional expectation E (Y | Xi) variance
V(E(Y|Xi))；
The formula carried according to sobol methods：
Composite parameters are calculated in drilling process for the global sensitivity index S of cutting forcei。
Beneficial effects of the present invention：
The Global sensitivity analysis method of the present invention can meet composite parameters to composite in drilling when institute
By the demand of the impact analysis of drill thrust, existing Sensitivity Analysis Method is overcome not to be suitable for mechanical characteristic dispersiveness and random
The deficiency of composite, when have studied the uncertainty of different input variables (Compound Material Engineering parameter) to composite drilling
Export the influence size of response quautity (drill thrust).
The present invention is advantageous to, to Simplified analysis process, save the time, improves efficiency.
Brief description of the drawings：
Fig. 1 is the analysis calculation process of the present invention.
Fig. 2 is sample calculation analysis finite element simulation figure of the present invention.
Fig. 3 is the analysis of cases sensitivity results figure of the present invention.
Embodiment：
The present invention is further illustrated with reference to the accompanying drawings and detailed description.
As Figure 1-3.
A kind of composite drilling process finite element Global sensitivity analysis method, it comprises the following steps：
First, multigroup random number is produced according to the distribution of composite relevant parameter to input as sample；
Secondly, the three-dimensional Drilling model of composite is established using finite element software ABAQUS/Explicit；
3rd, the corresponding cutting force of drilling process generation is emulated by finite element software and is used as sample output；
4th, each ginseng is asked according to " state-contingent parametric method " (States Dependent Parameter, SDP) respectively
Conditional expectation of the number for output；
5th, according to the Global sensitivity analysis method based on variance --- " sobol methods " analyzes different parameters for multiple
The size of the influence degree of material drill thrust.
Comprise the following steps that：
First, it is assumed that each parameter of composite is all to obey equally distributed stochastic variable, and it is uncorrelated, utilize
Matlab softwares produce the random sample X that sample size is N according to the distributed constant of each variablet(x1,x2,x3... xi) (t=1,2,
3 ... N) inputted as sample；
Secondly, the three-dimensional Drilling model of composite is established using finite element analysis software ABAQUS/Explicit；
3rd, call above-mentioned composite material model (i.e. the three-dimensional Drilling model of composite) n times, fortran language controls
Abaqus software emulations drilling process processed obtains the curve Y of corresponding cutting force time changet- t schemes, using in least square
Collimation method obtains the numerical value Y represented by the center line of curve in collimation method or arithmetic meant(t=1,2,3 ... N) exported as sample；
4th, parameters X is asked according to SDP methods respectively1, X2, X3,…XiFor sample output Y conditional expectation E (Y |
X1), E (Y | X2) ... E (Y | Xi) step is as follows：
According to control theory related content, state space SDP models are：
observation equation:Yt=Pi,t+ei,t；
SDP model：Pi,t=Pi,t-1+di,t-1
di,t=di,t-1+ηi,t,ηi,t~N (0, σ2)；
In formula：Pi,tFor with state variable Xi,tThe SDP state relevant parameters of change, it is basic variable Xi,tFunction, et
(observation interference), ηi,t(system disturbance) be average be zero, variance σ2And ση 2White Gaussian noise input, its value can use
Maximum likelihood optimization is estimated to obtain.
Then recursive Kalman filtering (kalam Filtreing, KF) is used in circulation backward fit procedure and is accordingly passed
Smooth (Fixed Interval Smoothing, the FIS) rule of fixed interval is returned to estimate to obtain each state relevant parameter Pi,t。
Finally substitute into formula：
E(Yt|Xi,t)=Yt-Pi,t；
Try to achieve each parameter for output conditional expectation E (Y | Xi,t)；
5th：According to Global sensitivity analysis method --- " sobol methods " calculates each parameter XiTo the global sensitive of output Y
Degree, its process are as follows：
Substitute into formula:
The gained sample of above-mentioned steps 2 output Y is tried to achieve respectivelytVariance V (Y), and each parameter X in step 3iFor Y condition
It is expected E (Y | Xi) variance V (E (and Y | Xi))。
The formula carried according to sobol methods：
Composite parameters are calculated in drilling process for the global sensitivity index S of cutting forcei。
With reference to T700 carbon fibre reinforced composites, details are as follows：
A kind of composite drilling process finite element Global sensitivity analysis method, its flow chart is as shown in figure 1, specific bag
Include following steps：
First, multigroup random number is produced according to the distribution of the relevant parameter of composite to input as sample：With T700 carbon
Exemplified by unidirectional (0 °) plate of fibre reinforced composites individual layer, such as Fig. 2.Its longitudinal stretching ultimate strength Y is obtained by testing1=
1900MPa, longitudinal ultimate compression strength Y2=1000MPa, cross directional stretch limiting chi1=84MPa, transverse compression limiting chi2=
250MPa, longitudinal Young's modulus E2=8.2GPa, horizontal Young's modulus E1=112GPa, Poisson's ratio ν=0.3, G=4500MPa.
Using these parameters as input variable, remaining parameter keeps constant, it is assumed that and each parameter is all to obey equally distributed stochastic variable,
And it is uncorrelated, produce the random sample X that sample size is 100 using uniform designt(t=1,2,3 ... 100) it is used as sample
Input, partial data such as table 1；
Secondly, the three-dimensional Drilling model of composite is established using finite element analysis software ABAQUS/Explicit, is such as schemed
Shown in 2, the cutter is that the cutter is 6 millimeters of drill bits, 118 ° of drift angle, 30 ° of helical angle.Kinematic constraint is in X and Y-direction, feeding speed
Degree is 1.5mm/s, drilling speed 314rad/s.Carbon fibre composite is multidirectional laminate, is 12mm*12mm*5mm's by size
40 are formed according to 0 °/- 45 °/90 °/45 ° directions stackings, and the laminate border is fixed.
3rd, control Abaqus finite element analysis softwares to call above-mentioned composite material model (i.e. multiple using fortran language
The three-dimensional Drilling model of condensation material) 100 times, pass through the song that software emulation drilling process produces corresponding cutting force time change
Line Yt- t schemes, and the numerical value Y represented by the center line of curve is obtained using collimation method in least squares method or arithmetic meant(t=1,
2,3 ... 100) export as sample, partial results such as table 1；
The material parameter of table 1. exports with drill thrust
4th, according to SDP methods ask respectively parameters for sample export Y conditional expectation E (Y | Xt) step is as follows：
According to control theory related content, state space SDP models are：
observation equation:Yt=Pi,t+ei,t；
SDP model：Pi,t=Pi,t-1+di,t-1
di,t=di,t-1+ηi,t；ηi,t~N (0, σ2)；
In formula：Pi,tFor with state variable Xi,tThe SDP state relevant parameters of change, it is basic variable Xi,tFunction, et
(observation interference), ηi,t(system disturbance) be average be zero, variance σ2And ση 2White Gaussian noise input, its value can use
Maximum likelihood optimization is estimated to obtain.
Then recursive Kalman filtering (kalam Filtreing, KF) is used in circulation backward fit procedure and is accordingly passed
Smooth (Fixed Interval Smoothing, the FIS) rule of fixed interval is returned to estimate to obtain each state relevant parameter Pi,t。
Finally substitute into formula：
E(Yt|Xi,t)=Yt-Pi,t
Try to achieve each parameter for output conditional expectation E (Y | Xi,t)；
5th：According to Global sensitivity analysis method --- " sobol methods " calculates each parameter XiTo the global sensitive of output Y
Degree, its process are as follows：
Substitute into formula:
The gained sample of above-mentioned steps 2 output Y is tried to achieve respectivelytVariance V (Y), and each parameter X in step 3iFor Y condition
It is expected E (Y | Xi) variance V (E (and Y | Xi))。
The formula carried according to sobol methods：
Composite parameters are calculated in drilling process for the global sensitivity index S of cutting forcei.Result of calculation
Such as table 2.
The Calculation of Sensitivity result of table 2.
For ease of analyzing to obtain block diagram, such as Fig. 3.Can with it is concluded that：The Young's modulus of material influences the most on drill thrust
Significantly, Poisson's ratio has negatively influencing to drill thrust size, and the coefficient of rigidity on drill thrust almost without influence.
Part that the present invention does not relate to is same as the prior art or can be realized using prior art.
Claims (6)
1. a kind of composite drilling process finite element Global sensitivity analysis method, its feature are as follows：
First, multigroup random number is produced according to the distribution of composite relevant parameter to input as sample；
Secondly, the three-dimensional Drilling model of composite is established using finite element software；
3rd, above-mentioned three-dimensional Drilling model is called, emulating drilling process by finite element software produces drill thrust-time diagram, profit
Corresponding cutting force is obtained by the use of least squares method to export as sample；
4th, parameters pair are asked according to " state-contingent parametric method " (States Dependent Parameter, SDP) respectively
In the conditional expectation of output；
5th, according to the Global sensitivity analysis method based on variance --- " sobol methods " analysis different parameters bore for multiple material
Cut the size of the influence degree of power.
2. according to the method for claim 1, it is characterized in that the relevant parameter of described composite is all to obey just to divide very much
The stochastic variable of cloth, the coefficient of variation are 0.06, and uncorrelated, are produced using Matlab softwares according to the distributed constant of each variable
The random sample X that raw sample size is Nt(x1,x2,x3…xi) (t=1,2,3 ... N) inputted as sample.
3. according to the method for claim 1, it is characterized in that described finite element software is ABAQUS/Explicit.
4. according to the method for claim 1, it is characterized in that calling three-dimensional Drilling model n times, pass through Abaqus software emulations
Drilling process produces corresponding cutting force Yt(t=1,2,3 ... N) exported as sample.
5. according to the method for claim 1, it is characterized in that seeking parameters X respectively according to SDP methods1, X2, X3,…XiFor
Sample output Y conditional expectation E (Y | X1), E (Y | X2) ... E (Y | Xi) step is as follows：
According to control theory, establishing state space SDP models is：
observation equation:Yt=Pi,t+ei,t；
SDP model:Pi,t=Pi,t-1+di,t-1；
di,t=di,t-1+ηi,t, ηi,t~N (0, σ2)；
In formula：Pi,tFor with state variable Xi,tThe SDP state relevant parameters of change, it is basic variable Xi,tFunction, et(observation is dry
Disturb), ηi,t(system disturbance) be average be zero, variance σ2And ση 2White Gaussian noise input, its value can use maximum likelihood
Optimal estimating obtains.
Then consolidated in circulation backward fit procedure with recursive Kalman filtering (kalam Filtreing, KF) and corresponding recurrence
Determine smooth (Fixed Interval Smoothing, the FIS) rule in section to estimate to obtain each state relevant parameter Pi,t。
Finally according to formula：
E(Yt|Xi,t)=Yt-Pi,t；
Calculate each parameter for output conditional expectation E (Y | Xi,t)；
6. according to the method for claim 1, it is characterized in that according to Global sensitivity analysis method --- " sobol methods " calculates
Each parameter XiThe process of global sensitivity to exporting Y is as follows：
According to formula:
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>Y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
<mo>-</mo>
<mi>E</mi>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mi>t</mi>
</msub>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>;</mo>
</mrow>
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mo>(</mo>
<mrow>
<mi>Y</mi>
<mo>|</mo>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mi>N</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>t</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mo>(</mo>
<mi>Y</mi>
<mo>|</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>-</mo>
<mi>E</mi>
<mo>(</mo>
<mi>E</mi>
<mo>(</mo>
<mrow>
<mi>Y</mi>
<mo>|</mo>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
<mo>&rsqb;</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>;</mo>
</mrow>
Calculate foregoing gained sample output YtVariance V (Y) and each parameter XiFor Y conditional expectation E (Y | Xi) variance V (E
(Y|Xi))；
The formula carried according to sobol methods：
<mrow>
<msub>
<mi>S</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>E</mi>
<mo>(</mo>
<mrow>
<mi>Y</mi>
<mo>|</mo>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>Y</mi>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
Composite parameters are calculated in drilling process for the global sensitivity index S of cutting forcei。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710961870.9A CN107491627B (en) | 2017-10-16 | 2017-10-16 | Finite element global sensitivity analysis method in composite material drilling process |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710961870.9A CN107491627B (en) | 2017-10-16 | 2017-10-16 | Finite element global sensitivity analysis method in composite material drilling process |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107491627A true CN107491627A (en) | 2017-12-19 |
| CN107491627B CN107491627B (en) | 2020-10-09 |
Family
ID=60653716
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710961870.9A Active CN107491627B (en) | 2017-10-16 | 2017-10-16 | Finite element global sensitivity analysis method in composite material drilling process |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107491627B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110232255A (en) * | 2019-06-24 | 2019-09-13 | 南京航空航天大学 | A kind of method of carbon fibre reinforced composite drilling axial force prediction |
| CN113111539A (en) * | 2021-05-03 | 2021-07-13 | 西北工业大学 | Fuzzy variable-based global sensitivity analysis method for series thermoelectric generator |
| CN114741792A (en) * | 2022-04-20 | 2022-07-12 | 合肥工业大学 | Method for calculating material parameter sensitivity of solid rocket engine nozzle |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103455679A (en) * | 2013-09-06 | 2013-12-18 | 湘潭电机股份有限公司 | Finite element analysis method of composite material |
| CN106202734A (en) * | 2016-07-14 | 2016-12-07 | 中国人民解放军国防科学技术大学 | global sensitivity analysis method based on Gaussian radial basis function |
-
2017
- 2017-10-16 CN CN201710961870.9A patent/CN107491627B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103455679A (en) * | 2013-09-06 | 2013-12-18 | 湘潭电机股份有限公司 | Finite element analysis method of composite material |
| CN106202734A (en) * | 2016-07-14 | 2016-12-07 | 中国人民解放军国防科学技术大学 | global sensitivity analysis method based on Gaussian radial basis function |
Non-Patent Citations (2)
| Title |
|---|
| 张屹尚 等: "复合随机振动系统的动力可靠性矩独立重要性测度及态相关参数解法", 《振动与冲击》 * |
| 阮文斌 等: "不确定条件下复合材料结构的全局灵敏度分析", 《复合材料学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110232255A (en) * | 2019-06-24 | 2019-09-13 | 南京航空航天大学 | A kind of method of carbon fibre reinforced composite drilling axial force prediction |
| CN113111539A (en) * | 2021-05-03 | 2021-07-13 | 西北工业大学 | Fuzzy variable-based global sensitivity analysis method for series thermoelectric generator |
| CN114741792A (en) * | 2022-04-20 | 2022-07-12 | 合肥工业大学 | Method for calculating material parameter sensitivity of solid rocket engine nozzle |
| CN114741792B (en) * | 2022-04-20 | 2024-02-13 | 合肥工业大学 | Method for calculating sensitivity of nozzle material parameters of solid rocket engine |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107491627B (en) | 2020-10-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107491627A (en) | Composite drilling process finite element Global sensitivity analysis method | |
| Boscolo et al. | Dynamic stiffness elements and their applications for plates using first order shear deformation theory | |
| Liao et al. | Subspace identification for fractional order Hammerstein systems based on instrumental variables | |
| Ahmadianfar et al. | Extracting optimal policies of hydropower multi-reservoir systems utilizing enhanced differential evolution algorithm | |
| D’Amore et al. | HPC computation issues of the incremental 3D variational data assimilation scheme in OceanVar software | |
| CN107958126A (en) | A kind of computational methods of bi-material layers v-notch stress intensity factor | |
| CN103020709A (en) | Optimization calculation method based on artificial bee colony algorithm and quantum-behaved particle swarm optimization algorithm | |
| Yang et al. | A delayed diffusive predator–prey system with Michaelis–Menten type predator harvesting | |
| Lam et al. | Uncertainty quantification for eigensystem-realization-algorithm, a class of subspace system identification | |
| CN113705878B (en) | Method and device for determining water yield of horizontal well, computer equipment and storage medium | |
| Zuo et al. | Reanalysis method for second derivatives of static displacement | |
| Yuan et al. | Application of differential evolution cuckoo search algorithm in parameter optimization of VG equation | |
| CN105932675A (en) | Parallel coordination algorithm for power flow of power system | |
| Akremi et al. | L∞/H∞ Functional Interval Observers Design for Multivariable Systems | |
| Rabiei et al. | Fuzzy regression model with interval-valued fuzzy input-output data | |
| Wang et al. | A short note on the boundedness analysis and control of the spatial fractal set from a kind of chain coupling logistic type map | |
| He et al. | A fuzzy control system for fitness service based on genetic algorithm during COVID-19 pandemic | |
| Yang et al. | Solving truss topological optimization via Swarm Intelligence | |
| Hu et al. | Perfect Manipulation of Open Stochastic Quantum Systems with Imperfect Detection | |
| Falco et al. | Analysis and optimal design of plates and shells under dynamic loads–I: finite element and sensitivity analysis | |
| Wang | Uzawa method for fuzzy linear system | |
| BR N et al. | Simulation of the drilling process in GFRP composites using system dynamics and validation by ANN and RSM | |
| Liao et al. | Analysis of the Solution of Time-Varying Linear Equations Based on the Method of Neural Dynamics | |
| Yang et al. | Solving truss topological optimization with discrete design variables via swarm intelligence | |
| Kora et al. | Frequency-and dissipation-dependent entanglement advantage in spin-network quantum reservoir computing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20211227 Address after: No. 29, Qinhuai District, Qinhuai District, Nanjing, Jiangsu Patentee after: Nanjing University of Aeronautics and Astronautics Asset Management Co.,Ltd. Address before: No. 29, Qinhuai District, Qinhuai District, Nanjing, Jiangsu Patentee before: Nanjing University of Aeronautics and Astronautics |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220207 Address after: Room 703-29, block C, Kechuang headquarters building, No. 320, pubin Road, Nanjing, Jiangsu 211899 Patentee after: Jiangsu Hangpu Guochuang Composite Material Research Institute Co.,Ltd. Address before: No. 29, Qinhuai District, Qinhuai District, Nanjing, Jiangsu Patentee before: Nanjing University of Aeronautics and Astronautics Asset Management Co.,Ltd. |