CN101614272B - Reshaping method of automobile torque converter blade lattice system - Google Patents

Abstract

A reshaping method of an automobile torque converter blade lattice system comprises the following steps: 1) acquiring geometric information of the existing torque converter blade lattice system; 2) calculating target inlet and outlet angles of blades and the inner and outer ring intersecting lines of a circulating circle; 3) entering a genetic algorithm program to start iterative computation; 4) randomly modifying each point coordinate to generate an initial group; 5) calculating the inlet and outlet angles of the blades and the inner and outer ring intersecting lines of the circulating circle after being modified; 6) judging the matching degree of the inlet and outlet angles of the blades and the inner and outer ring intersecting lines of the circulating circle after being modified with the target inlet and outlet angles obtained by calculation in step 2; and 7) if the matching degree is high, inheriting the individual to the next generation for further intersection and mutation to obtain point coordinate combination with higher matching degree; and if the matching degree is low, eliminating the individual. The invention has the advantages of modification of the existing torque converter CAD model, thus realizing the optimization of the torque converter.

Description

The automobile torque converter blade system change the shape method

Technical field

The present invention relates to a kind of design method of automobile torque converter, particularly a kind of method that existing fluid torque converter blade system is changed shape.

Background technique

Fluid torque converter is one of critical component of automatic transmission system, and its performance plays material impact to vehicle performance.Because the complexity of fluid torque converter flow field, one dimension beam model and posterior infromation that its design is based on simplification are substantially manufactured experimently test repeatedly.The shortcoming of this method is that error is big, and the design cycle is long, and possibly can't reach optimal performance.

Summary of the invention

The invention provides and a kind ofly utilize existing torque converter CAD model that the fluid torque converter blade system is optimized, under the prerequisite that keeps original flow path features, is guide parameters with blade into and out of bicker, realization is to the modification of blade shape, thus reach the automobile torque converter blade system that improves the fluid torque converter performance change the shape method.

In order to solve above technical problem, what the invention provides a kind of automobile torque converter blade system changes the shape method, and this changes the shape method and comprises the steps:

1. obtain existing fluid torque converter blade system geological information;

2. calculate target on blade and the circulate circle inner and outer rings intersection into and out of bicker;

3. enter genetic algorithm program and begin iterative computation;

4. each point coordinates of random modification generates initial population;

5. calculate to revise on rear blade and the circulate circle inner and outer rings intersection into and out of bicker;

6. judge the matching degree of revising on rear blade and the circulate circle inner and outer rings intersection into and out of bicker and the 2. middle calculating gained of step target turnover bicker;

7. the matching degree height carries out the genetic operator operation, if reach maximum evolutionary generation, then changes next step; Matching degree is low, and then this individuality is eliminated;

8. obtain two curves being asked, make up new cad model.

Described obtain existing fluid torque converter blade system geological information be: according to existing torque converter CAD model, obtain the coordinate information of some well-distributed pointses on each impeller blade inner and outer rings streamline respectively.

Target on described calculating blade and the circulate circle inner and outer rings intersection into and out of the method for bicker is:

Blade center line of flow path upper blade is enterprising into and out of bicker and inner and outer rings, satisfy the described relation of anti-potential flow theories between the exit angle, that is:

$\left\{\begin{array}{c}{\cot \mathrm{\β}}_i=\frac{r_i\×\cot \mathrm{\β}}{r}\\ \cot {\mathrm{\β}}_o=\frac{r_o\×\cot \mathrm{\β}}{r}\end{array}\right.---\left(I\right)$

Wherein, β represents that the blade center line of flow path is into and out of bicker;

β _iEnterprising, the exit angle of ring in the expression blade;

β _oEnterprising, the exit angle of expression blade outer shroud;

R represents blade center line of flow path inlet/outlet place radius;

r _iEnterprising, the outlet port radius of ring in the expression blade;

r _oEnterprising, the outlet port radius of expression blade outer shroud.

The method of described each point coordinates of modification is:

Still being positioned on the circulate circle inner and outer rings curved surface for guarantee revising the back each point, realizing modification on the plane of vertical axis each point coordinates by each point is rotated to an angle around central shaft; Revising each point coordinates of back is:

$\left\{\begin{array}{c}{x}_i^{\′}=x_i\\ y_i^{\′}=\sqrt{y_i^2+z_i^2}\×\cos \left(\right|\arctan (z_i/y_i)+{\mathrm{\θ}}_i\left|\right)\\ z_i^{\′}=-\sqrt{y_i^2+z_i^2}\×\sin \left(\right|\arctan (z_i/y_i)+{\mathrm{\θ}}_i\left|\right)\end{array}\right.---\left(\mathrm{II}\right)$

Wherein, (x _i, y _i, z _i) for revising the coordinate of preceding each point; θ _iBe the each point angle of swing.

In the present invention, according to the blade inlet/outlet curved surface scope that changes, fixed intermediate portion is divided some point coordinates, and the spoon of blade of institute's fixed point will not change.

The method into and out of bicker that described calculating is revised on rear blade and the circulate circle inner and outer rings intersection is:

At first use method of least squares will revise rear blade and become 3 times space curve with inner and outer rings intersection each point coordinate fitting:

$\left\{\begin{array}{c}{x}^{\′}{=a}_0+a_{1}i+a_2{i}^2+a_3{i}^3\\ y^{\′}=b_0+b_{1}i+b_2{i}^2+b_3{i}^3\\ z^{\′}=c_0+c_{1}i+c_2{i}^2+c_3{i}^3\end{array}\right.---\left(\mathrm{III}\right)$

Wherein, i is the each point sequence number;

Intersection as blade and inner and outer rings is described by N point respectively, and definition inlet's point sequence number is 1, and some sequence number in outlet port is N.

Wherein the calculating formula of every coefficient is:

$\underset{k=0}{\overset{n}{\mathrm{\Σ}}}\left(\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{i}^{j+k}\right)a_k=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{i}^j{x}_j^{\′}$

$\underset{k=0}{\overset{n}{\mathrm{\Σ}}}\left(\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{i}^{j+k}\right)b_k=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{i}^j{y}_j^{\′},j=\mathrm{0,1},\·\·\·N$

$\underset{k=0}{\overset{n}{\mathrm{\Σ}}}\left(\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{i}^{j+k}\right)c_k=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{i}^j{z}_j^{\′}$

Inner and outer rings upper blade turnover bicker is:

$\mathrm{\θ}=\arccos \left(\frac{T\·N}{\left|T\right|\left|N\right|}\right)---\left(\mathrm{IV}\right)$

T is in blade and the circulate circle, the tangent vector of outer shroud intersection inlet, and the method for asking according to the space curve tangential equation has following result:

If the parametric equation of certain space curve is:

$\left\{\begin{array}{c}x=f_x\left(t\right)\\ y=f_y\left(t\right)\\ z=f_z\left(t\right)\end{array}\right.$

Then it is at (x ₀, y ₀, z ₀) point tangential equation be:

$\frac{x-f_x\left(t_0\right)}{f_x^{\′}\left(t_0\right)}=\frac{y-f_y\left(t_0\right)}{f_y^{\′}\left(t_0\right)}=\frac{z-f_z\left(t_0\right)}{f_z^{\′}\left(t_0\right)}$

The curvilinear equation that obtains according to match:

$\left\{\begin{array}{c}x=a_0+a_{1}t+a_2{t}^2+a_3{t}^3\\ y=b_0+b_{1}t+b_2{t}^2+b_3{t}^3\\ z=c_0+c_{1}t+c_2{t}^2+c_3{t}^3\end{array}\right.$

Can get:

Inlet, T=(a ₁+ 2a ₂+ 3a ₃, b ₁+ 2b ₂+ 3b ₃, c ₁+ 2c ₂+ 3c ₃);

The outlet port, T=(a ₁+ 2a ₂N+3a ₃N ², b ₁+ 2b ₂N+3b ₃N ², c ₁+ 2c ₂N+3c ₃N ²);

N is blade and circulate circle inside and outside ring intersection inlet peripheral velocity tangent vector, and the method for asking according to the space curve tangential equation has following result:

The parametric equation of circle is:

$\left\{\begin{array}{c}x=x_o\\ y=R\sin \mathrm{\θ}\\ z=R\cos \mathrm{\θ}\end{array}\right.$

Wherein, o≤θ≤2 π.

It is at (x ₀, y ₀, z ₀) point, promptly curve is at θ=θ ₀The tangential equation at place is:

$\left\{\begin{array}{c}x=x_0\\ \frac{y-R{\sin \mathrm{\θ}}_0}{{R\cos \mathrm{\θ}}_0}=\frac{z-R\cos {\mathrm{\θ}}_0}{-R\sin {\mathrm{\θ}}_0}\end{array}\right.$

Tangent vector is:

N＝(0，Rcosθ ₀，-Rsinθ ₀)。

Wherein, ${R\cos \mathrm{\θ}}_0=\sqrt{y_0^2+z_0^2}{\cos \mathrm{\θ}}_0;$ ${-R\sin \mathrm{\θ}}_0=-\sqrt{y_0^2+z_0^2}{\sin \mathrm{\θ}}_0,$ (x when calculating inlet angle ₀, y ₀, z ₀) be the inlet point coordinate, θ ₀Be this angle; (x when calculating bicker ₀, y ₀, z ₀) be the exit point coordinate, θ ₀Be this angle.

Described genetic algorithm comprises following step:

1) sets up the objective function collection of optimization problem;

2) coding and initialization population;

3) calculate each individual fitness;

4) judge whether to satisfy condition?

Matching degree is low, then eliminates; Matching degree is high, then enters next step;

5) genetic operator operation reaches maximum evolutionary generation;

6) finish.

Wherein: 1, optimization aim is:

$\left\{\begin{array}{c}\underset{x\&Element;\mathrm{\&Omega;}}{\max F\left(\mathrm{\&theta;}\right)}=1/|{\mathrm{\&beta;}}_{\mathrm{cal}}-{\mathrm{\&beta;}}_{t\arg \mathrm{et}}|\\ \mathrm{st}.t_{j\min }<t_j\left(x\right)<t_{j\max },j=\mathrm{1,2}.......n\end{array}\right.$

Wherein, Ω represents solution space;

θ represents the optimization variable vector;

β _CalExpression is revised rear blade inner and outer rings design path into and out of bicker;

β _TargetExpression blade inner and outer rings streamline target is into and out of bicker;

t _j(θ) vector of constraint is optimized in expression;

t _JminVector is optimized in expression needs satisfied minimum value;

t _JmaxVector is optimized in expression needs satisfied maximum value.

2, optimization variable is:

θ＝[θ ₁，θ ₂，…，θ _i] ^T

Wherein, θ _iCurve is enterprising in order to be asked, outlet section is modified the each point rotation angle value.

Optimize constraint:

${\mathrm{\&theta;}}_i=\frac{\mathrm{temp}}{\mathrm{chrom}\_\mathrm{value}\&times;({\mathrm{param}}_b-{\mathrm{param}}_s)+{\mathrm{param}}_s}---\left(V\right)$

Wherein, param _bThe maximum value of expression optimization range;

Param _sThe minimum value of expression optimization range;

Chrom_value is determined by chromosome quantity, gets 63.

3, coding and colony's initialization:

Adopt binary coding, that is:

θ _i＝binary _i

Wherein, binary is 6 bits;

Adopt consistent random fashion to choose initial population, the chromosome number is 18, and colony's quantity is 1000.

Fitness function is:

$f\left({\mathrm{\&theta;}}_i\right)=\frac{1}{|{\mathrm{\&beta;}}_{\mathrm{cal}}-{\mathrm{\&beta;}}_{t\arg \mathrm{et}}|}---\left(\mathrm{VI}\right)$

When not satisfying constraint conditio, promptly for inlet angle θ ₁＞θ ₂＞...＞θ _i＞0 or θ ₁＜θ ₂＜...＜θ _i＜0, for exit angle θ _N＞θ _N-1＞...＞θ _N-i＞0 or θ _N＜θ _N-1＜...＜θ _N-iFitness function was in＜0 o'clock:

$f\left({\mathrm{\&theta;}}_i\right)=\frac{1}{|{\mathrm{\&beta;}}_{\mathrm{cal}}-{\mathrm{\&beta;}}_{t\arg \mathrm{et}}|+10}---\left(\mathrm{VII}\right)$

4, genetic operator:

Select operator, adopt the Monte Carlo to select.Each individual selecteed probability and its fitness value are proportional, the individuality that the fitness function value is high, and its covering of the fan ratio shared on wheel disc is then high, and selected probability is also big.Mathematical expression is:

$P_{\mathrm{ri}}=f\left(X_i\right)/\underset{j=1}{\overset{\mathrm{PopSize}}{\mathrm{\&Sigma;}}}f\left(X_j\right)$

Wherein, i=1,2 ... M.

Crossover operator adopts non-consistent arithmetic crossover operator.Mathematical expression is:

binary ₁′＝λ _ibinary ₁+(1-λ _i)binary ₂

binary ₂′＝(1-λ _i)binary ₁+λ _ibinary ₂i∈{1，2，…10}

Wherein: binary ₁' be the situation before the gene conversion for the first time;

Binary ₁Be the situation after the gene conversion for the first time;

Binary ₂' be the situation after the gene conversion for the second time;

Binary ₂Situation after the conversion for the second time of branch gene;

λ _iFor obeying equally distributed random numbers.

Mutation operator, even mutation operation is adopted in the variation computing.Select to produce the individual θ of variation during variation _i, be binary corresponding to the binary form in the colony _i, the position that produces mutant gene then at random changes 1 according to the binary system characteristics into 0, perhaps changes 0 into 1, produces the individuality after the new variation.

5, the Operational Limits of genetic operator:

Code length is 18;

Group size is 1000;

Crossover probability is 0.6;

The variation probability is 0.02;

Stopping algebraically is 200.

Superior effect of the present invention is: the present invention helps reducing the actual loading test number of times in conjunction with the CFD simulation calculation, shortens the fluid torque converter development and Design cycle, reduces development cost, improves fluid torque converter enterprise product seriation degree.

Description of drawings

Fig. 1 is the shape method flow diagram that changes of the present invention;

Fig. 2 is enterprising for blade inside and outside ring of the present invention, the definition figure of exit angle;

Fig. 3 is the flow chart of genetic algorithm of the present invention;

Fig. 4 is a turbine new model of the present invention;

Fig. 5 is a runner model of the present invention.

Embodiment

See also shown in the accompanying drawing, the invention will be further described.

As shown in Figure 1, what the invention provides a kind of automobile torque converter blade system changes the shape method, and this changes the shape method and comprises the steps:

1. obtain existing fluid torque converter blade system geological information;

2. calculate target on blade and the circulate circle inner and outer rings intersection into and out of bicker;

3. enter genetic algorithm program and begin iterative computation;

4. each point coordinates of random modification generates initial population;

5. calculate to revise on rear blade and the circulate circle inner and outer rings intersection into and out of bicker;

6. judge the matching degree of revising on rear blade and the circulate circle inner and outer rings intersection into and out of bicker and the 2. middle calculating gained of step target turnover bicker;

7. matching degree height carries out the genetic operator operation, and this individual inheritance to of future generation, further intersected and make a variation, in the hope of obtaining the higher point coordinates combination of matching degree, as if reaching maximum evolutionary generation, then changes next step;

Matching degree is low, and then this individuality is eliminated;

8. obtain two curves being asked, make up new cad model.

Described obtain existing fluid torque converter blade system geological information be: according to existing torque converter CAD model, obtain the coordinate information of some well-distributed pointses on each impeller blade inner and outer rings streamline respectively.

Target on described calculating blade and the circulate circle inner and outer rings intersection into and out of the method for bicker is:

Blade center line of flow path upper blade is enterprising into and out of bicker and inner and outer rings, satisfy the described relation of anti-potential flow theories between the exit angle, that is:

$\left\{\begin{array}{c}{\cot \mathrm{\&beta;}}_i=\frac{r_i\&times;\cot \mathrm{\&beta;}}{r}\\ \cot {\mathrm{\&beta;}}_o=\frac{r_o\&times;\cot \mathrm{\&beta;}}{r}\end{array}\right.---\left(I\right)$

Wherein, β represents that the blade center line of flow path is into and out of bicker;

β _iEnterprising, the exit angle of ring in the expression blade;

β _oEnterprising, the exit angle of expression blade outer shroud;

R represents blade center line of flow path inlet/outlet place radius;

r _iEnterprising, the outlet port radius of ring in the expression blade;

r _oEnterprising, the outlet port radius of expression blade outer shroud.

Described blade inside and outside ring is enterprising, exit angle is defined as relative velocity and satellite speed between angle, can represent to form by two plane angle beta_xy and beta_yz, as shown in Figure 2.

The method of described each point coordinates of modification is:

Still being positioned on the circulate circle inner and outer rings curved surface for guarantee revising the back each point, realizing modification on the plane of vertical axis each point coordinates by each point is rotated to an angle around central shaft; Revising each point coordinates of back is:

$\left\{\begin{array}{c}{x}_i^{\&prime;}=x_i\\ y_i^{\&prime;}=\sqrt{y_i^2+z_i^2}\&times;\cos \left(\right|\arctan (z_i/y_i)+{\mathrm{\&theta;}}_i\left|\right)\\ z_i^{\&prime;}=-\sqrt{y_i^2+z_i^2}\&times;\sin \left(\right|\arctan (z_i/y_i)+{\mathrm{\&theta;}}_i\left|\right)\end{array}\right.---\left(\mathrm{II}\right)$

Wherein, (x _i, y _i, z _i) for revising the coordinate of preceding each point; θ _iBe the each point angle of swing.

In the present invention, according to the blade inlet/outlet curved surface scope that changes, fixed intermediate portion is divided some point coordinates, and the spoon of blade of institute's fixed point will not change.

The method into and out of bicker that described calculating is revised on rear blade and the circulate circle inner and outer rings intersection is:

At first use method of least squares will revise rear blade and become 3 times space curve with inner and outer rings intersection each point coordinate fitting:

$\left\{\begin{array}{c}{x}^{\&prime;}{=a}_0+a_{1}i+a_2{i}^2+a_3{i}^3\\ y^{\&prime;}=b_0+b_{1}i+b_2{i}^2+b_3{i}^3\\ z^{\&prime;}=c_0+c_{1}i+c_2{i}^2+c_3{i}^3\end{array}\right.---\left(\mathrm{III}\right)$

Wherein, i is the each point sequence number;

Intersection as blade and inner and outer rings is described by N point respectively, and definition inlet's point sequence number is 1, and some sequence number in outlet port is N.

Wherein the calculating formula of every coefficient is:

$\underset{k=0}{\overset{n}{\mathrm{\&Sigma;}}}\left(\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}^{j+k}\right)a_k=\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}^j{x}_j^{\&prime;}$

$\underset{k=0}{\overset{n}{\mathrm{\&Sigma;}}}\left(\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}^{j+k}\right)b_k=\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}^j{y}_j^{\&prime;},j=\mathrm{0,1},\&CenterDot;\&CenterDot;\&CenterDot;N$

$\underset{k=0}{\overset{n}{\mathrm{\&Sigma;}}}\left(\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}^{j+k}\right)c_k=\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{i}^j{z}_j^{\&prime;}$

Inner and outer rings upper blade turnover bicker is:

$\mathrm{\&theta;}=\arccos \left(\frac{T\&CenterDot;N}{\left|T\right|\left|N\right|}\right)---\left(\mathrm{IV}\right)$

T is in blade and the circulate circle, the tangent vector of outer shroud intersection inlet, and the method for asking according to the space curve tangential equation has following result:

If the parametric equation of certain space curve is:

$\left\{\begin{array}{c}x=f_x\left(t\right)\\ y=f_y\left(t\right)\\ z=f_z\left(t\right)\end{array}\right.$

Then it is at (x ₀, y ₀, z ₀) point tangential equation be:

$\frac{x-f_x\left(t_0\right)}{f_x^{\&prime;}\left(t_0\right)}=\frac{y-f_y\left(t_0\right)}{f_y^{\&prime;}\left(t_0\right)}=\frac{z-f_z\left(t_0\right)}{f_z^{\&prime;}\left(t_0\right)}$

The curvilinear equation that obtains according to match:

$\left\{\begin{array}{c}x=a_0+a_{1}t+a_2{t}^2+a_3{t}^3\\ y=b_0+b_{1}t+b_2{t}^2+b_3{t}^3\\ z=c_0+c_{1}t+c_2{t}^2+c_3{t}^3\end{array}\right.$

Can get:

Inlet, T=(a ₁+ 2a ₂+ 3a ₃, b ₁+ 2b ₂+ 3b ₃, c ₁+ 2c ₂+ 3c ₃);

The outlet port, T=(a ₁+ 2a ₂N+3a ₃N ², b ₁+ 2b ₂N+3b ₃N ², c ₁+ 2c ₂N+3c ₃N ²);

N is blade and circulate circle inside and outside ring intersection inlet peripheral velocity tangent vector, and the method for asking according to the space curve tangential equation has following result:

The parametric equation of circle is:

$\left\{\begin{array}{c}x=x_o\\ y=R\sin \mathrm{\&theta;}\\ z=R\cos \mathrm{\&theta;}\end{array}\right.$

Wherein, o≤θ≤2 π.

It is at (x ₀, y ₀, z ₀) point, promptly curve is at θ=θ ₀The tangential equation at place is:

$\left\{\begin{array}{c}x=x_0\\ \frac{y-R{\sin \mathrm{\&theta;}}_0}{{R\cos \mathrm{\&theta;}}_0}=\frac{z-R\cos {\mathrm{\&theta;}}_0}{-R\sin {\mathrm{\&theta;}}_0}\end{array}\right.$

Tangent vector is:

N＝(0，Rcosθ ₀，-Rsinθ ₀)。

Wherein, ${R\cos \mathrm{\&theta;}}_0=\sqrt{y_0^2+z_0^2}{\cos \mathrm{\&theta;}}_0;$ ${-R\sin \mathrm{\&theta;}}_0=-\sqrt{y_0^2+z_0^2}{\sin \mathrm{\&theta;}}_0,$ (x when calculating inlet angle ₀, y ₀, z ₀) be the inlet point coordinate, θ ₀Be this angle; (x when calculating bicker ₀, y ₀, z ₀) be the exit point coordinate, θ ₀Be this angle.

As shown in Figure 3, described genetic algorithm comprises following step:

1, sets up the objective function collection of optimization problem;

2, coding and initialization population;

3, calculate each individual fitness;

4, judge whether to satisfy condition?

Matching degree is low, then eliminates; Matching degree is high, then enters next step;

5, genetic operator operation reaches maximum evolutionary generation;

6, finish.

By above step, change the shape design and optimization in conjunction with what the analysis result of computation fluid dynamics (CFD) just can be finished fluid torque converter, optimizing the cad model that obtains can generate prototype fast by auxiliary make (CAM) of modern computer, if verify that by experiment satisfying performance requirement just can carry out volume production.Therefore the present invention has higher utility.

Wherein: 1, optimization aim is:

$\left\{\begin{array}{c}\underset{x\&Element;\mathrm{\&Omega;}}{\max F\left(\mathrm{\&theta;}\right)}=1/|{\mathrm{\&beta;}}_{\mathrm{cal}}-{\mathrm{\&beta;}}_{t\arg \mathrm{et}}|\\ \mathrm{st}.t_{j\min }<t_j\left(x\right)<t_{j\max },j=\mathrm{1,2}.......n\end{array}\right.$

Wherein, Ω represents solution space;

θ represents the optimization variable vector;

β _CalExpression is revised rear blade inner and outer rings design path into and out of bicker;

β _TargetExpression blade inner and outer rings streamline target is into and out of bicker;

t _j(θ) vector of constraint is optimized in expression;

t _{J min}Vector is optimized in expression needs satisfied minimum value;

t _{J max}Vector is optimized in expression needs satisfied maximum value.

2, optimization variable is:

θ＝[θ ₁，θ ₂，…，θ _i] ^T

Wherein, θ _iCurve is enterprising in order to be asked, outlet section is modified the each point rotation angle value.

Optimize constraint:

${\mathrm{\&theta;}}_i=\frac{\mathrm{temp}}{\mathrm{chrom}\_\mathrm{value}\&times;({\mathrm{param}}_b-{\mathrm{param}}_s)+{\mathrm{param}}_s}---\left(V\right)$

Wherein, param _bThe maximum value of expression optimization range;

Param _sThe minimum value of expression optimization range;

Chrom_value is determined by chromosome quantity, gets 63.

3, coding and colony's initialization:

Adopt binary coding, that is:

θ _i＝binary _i

Wherein, binary is 6 bits;

Adopt consistent random fashion to choose initial population, the chromosome number is 18, and colony's quantity is 1000.

Fitness function is:

$f\left({\mathrm{\&theta;}}_i\right)=\frac{1}{|{\mathrm{\&beta;}}_{\mathrm{cal}}-{\mathrm{\&beta;}}_{t\arg \mathrm{et}}|}---\left(\mathrm{VI}\right)$

When not satisfying constraint conditio, promptly for inlet angle θ ₁＞θ ₂＞...＞θ _i＞0 or θ ₁＜θ ₂＜...＜θ _i＜0, for exit angle θ _N＞θ _N-1＞...＞θ _N-i＞0 or θ _N＜θ _N-1＜...＜θ _N-iFitness function was in＜0 o'clock:

$f\left({\mathrm{\&theta;}}_i\right)=\frac{1}{|{\mathrm{\&beta;}}_{\mathrm{cal}}-{\mathrm{\&beta;}}_{t\arg \mathrm{et}}|+10}---\left(\mathrm{VII}\right)$

4, genetic operator:

Select operator, adopt the Monte Carlo to select.Each individual selecteed probability and its fitness value are proportional, the individuality that the fitness function value is high, and its covering of the fan ratio shared on wheel disc is then high, and selected probability is also big.Mathematical expression is:

$P_{\mathrm{ri}}=f\left(X_i\right)/\underset{j=1}{\overset{\mathrm{PopSize}}{\mathrm{\&Sigma;}}}f\left(X_j\right)$

Wherein, i=1,2 ... M.

Crossover operator adopts non-consistent arithmetic crossover operator.Mathematical expression is:

binary ₁′＝λ _ibinary ₁+(1-λ _i)binary ₂

binary ₂′＝(1-λ _i)binary ₁+λ _ibinary ₂i∈{1，2，…10}

Wherein: binary ₁' be the situation before the gene conversion for the first time;

Binary ₁Be the situation after the gene conversion for the first time;

Binary ₂' be the situation after the gene conversion for the second time;

Binary ₂Situation after the conversion for the second time of branch gene;

λ _iFor obeying equally distributed random numbers.

Mutation operator, even mutation operation is adopted in the variation computing.Select to produce the individual θ of variation during variation _i, be binary corresponding to the binary form in the colony _i, the position that produces mutant gene then at random changes 1 according to the binary system characteristics into 0, perhaps changes 0 into 1, produces the individuality after the new variation.

5, the Operational Limits of genetic operator:

Code length is 18;

Group size is 1000;

Crossover probability is 0.6;

The variation probability is 0.02;

Stopping algebraically is 200.

Existing a hydraulic pitch device cad model based on Pro/E software is decided to be turbine blade and changes the shape object.Former turbine blade center line of flow path is respectively 35 ° and 147 ° into and out of bicker, changes the center line of flow path inlet angle into 30 ° in the present embodiment, and the center line of flow path exit angle changes 130 ° into.

1. measure point coordinates

Utilize Pro/E to measure the instrument of point coordinates, get each point coordinates.

With the turbine is example, and optional blade on the curve that its pressure side and circulate circle inside and outside ring intersect, roughly measures 9 points (is example with 9 points) equably, writes down its three-coordinate value, attention, and these 9 points must comprise inlet summit and summit, outlet port.

2. caller calculates the blade and the circulate circle inside and outside ring intersection that change after the shape, promptly is respectively the curvilinear equation of blade and inside and outside ring intersection.

3. the rear blade inner and outer rings of will retrofiting streamline imports among the Pro/E, finishes new torque-converters cad model and makes up, and finally obtains the turbine new model, as shown in Figure 4.

4. set up the runner model

Set up new fluid torque converter runner model,, carry out the CFD simulation calculation, provide the simulation calculation foundation for next round changes shape as Fig. 5.

Claims (2)

1. An automobile torque converter blade system change the shape method, it is characterized in that: this changes the shape method and comprises the steps:

   1. obtain existing fluid torque converter blade system geological information;

   2. calculate target on blade and the circulate circle inner and outer rings intersection into and out of bicker;

   3. enter genetic algorithm program and begin iterative computation;

   4. each point coordinates of random modification generates initial population;

   5. calculate to revise on rear blade and the circulate circle inner and outer rings intersection into and out of bicker;

   6. judge the matching degree of revising on rear blade and the circulate circle inner and outer rings intersection into and out of bicker and the 2. middle calculating gained of step target turnover bicker;

   7. the matching degree height carries out the genetic operator operation, if reach maximum evolutionary generation, then changes next step; Matching degree is low, and then this individuality is eliminated;

   8. obtain two curves being asked, make up new cad model;

   Target on described calculating blade and the circulate circle inner and outer rings intersection into and out of the method for bicker is:

   Blade center line of flow path upper blade is enterprising into and out of bicker and inner and outer rings, satisfy the described relation of anti-potential flow theories between the exit angle, that is:

   

   Wherein, β represents that the blade center line of flow path is into and out of bicker;

   β _iEnterprising, the exit angle of ring in the expression blade;

   β _oEnterprising, the exit angle of expression blade outer shroud;

   R represents blade center line of flow path inlet/outlet place radius;

   r _iEnterprising, the outlet port radius of ring in the expression blade;

   r _oEnterprising, the outlet port radius of expression blade outer shroud;

   The method of described each point coordinates of modification is:

   Still being positioned on the circulate circle inner and outer rings curved surface for guarantee revising the back each point, realizing modification on the plane of vertical axis each point coordinates by each point is rotated to an angle around central shaft; Revising each point coordinates of back is:

   Wherein, (x _i, y _i, z _i) for revising the coordinate of preceding each point; θ _iBe the each point angle of swing;

   The method into and out of bicker that described calculating is revised on rear blade and the circulate circle inner and outer rings intersection is:

   At first use method of least squares will revise rear blade and become 3 times space curve with inner and outer rings intersection each point coordinate fitting:

   

   Wherein, i is the each point sequence number;

   Intersection as blade and inner and outer rings is described by N point respectively, and definition inlet's point sequence number is 1, and some sequence number in outlet port is N;

   Wherein the calculating formula of every coefficient is:

   

   

   

   Inner and outer rings upper blade turnover bicker is:

   

   T is in blade and the circulate circle, the tangent vector of outer shroud intersection inlet, and the method for asking according to the space curve tangential equation has following result:

   If the parametric equation of certain space curve is:

   

   Then it is at (x ₀, y ₀, z ₀) point tangential equation be:

   The curvilinear equation that obtains according to match:

   

   Can get:

   Inlet, T=(a ₁+ 2a ₂+ 3a ₃, b ₁+ 2b ₂+ 3b ₃, c ₁+ 2c ₂+ 3c ₃);

   The outlet port, T=(a ₁+ 2a ₂N+3a ₃N ², b ₁+ 2b ₂N+3b ₃N ², c ₁+ 2c ₂N+3c ₃N ²);

   N is blade and circulate circle inside and outside ring intersection inlet peripheral velocity tangent vector, and the method for asking according to the space curve tangential equation has following result:

   The parametric equation of circle is:

   

   Wherein, o≤θ≤2 π;

   It is at (x ₀, y ₀, z ₀) point, promptly curve is at θ=θ ₀The tangential equation at place is:

   

   Tangent vector is:

   N＝(0，Rcosθ ₀，-Rsinθ ₀)?;

   Wherein,

   

   

   (x when calculating inlet angle ₀, y ₀, z ₀) be the inlet point coordinate, θ ₀Be this angle; (x when calculating bicker ₀, y ₀, z ₀) be the exit point coordinate, θ ₀Be this angle.
2. By the described automobile torque converter blade system of claim 1 change the shape method, it is characterized in that:

   Described genetic algorithm comprises following step:

   1) sets up the objective function collection of optimization problem;

   2) coding and initialization population;

   3) calculate each individual fitness;

   4) judge whether to satisfy condition?

   The individuality that matching degree is low can be eliminated in the next round loop iteration;

   Matching degree is high, then enters next step;

   5) genetic operator operation reaches maximum evolutionary generation;

   6) finish.

Images