1. Introduction
To improve the quality of machine products, applying topological optimization technology to actual product structure design has received much attention. Since the 1980s, in the industrial world such as aerospace and automobiles, the structural analysis technique using the finite element method and a mathematical programming method had been the norm [
1]. Research on structural optimization for obtaining the optimum structural composition has been actively conducted after satisfying restrictive conditions such as strength and rigidity [
2].
Since the beginning of the 2000s, researchers have improved the analysis technology using supercomputers and other hardware, collision analysis and noise analysis, along with the approximate analytical model, have introduced the experimental design method and the neural network, and the range of study of optimization and the scale of analysis model are further expanded.
Recently, with the progress of analytical technology, the need for structural topological optimization has been increasing more and more from the field of design development. Regarding conventional structural shape optimization, depending on the type of design variable, it can be classified into size optimization, phase optimization, and shape optimization.
During size optimization, the most general setting should consist of a determination for every point in space, whether there is material at that point or not, such Azegami, H. proposed as a solution to shape optimization problems of domain boundaries in which boundary value problems of partial differential equations are defined [
3]. Alternatively, for a FEM (Finite element
method) discretization, every element is a potential void or structural member. Since the design variables of dimension optimization are simple, it is currently the situation that there are relatively few cases that can be applied to complicated structural topological optimization problems.
During topological optimization, the element types and element densities of the analysis model are used as design variables, and the number of elements increases or decreases during the optimization analysis [
4]. He, D., et al., mentions only the critical substructure to which the internal force of the structure is transmitted remains, when the so-called optimum phase structure can be obtained. However, the shape of the optimum structure obtained by phase optimization is difficult to predict in advance, it is difficult to use at the stage of detailed design, and there are features to be applied to the initial design [
5,
6].
During topological optimization, by changing the boundary shape of the structure using a parameter of a combination of appropriate interpolation functions, and by directly changing the design variable via the design variable, the design variable is changed in the middle of the optimization analysis, therefore, it is possible to control the shape change [
7]. Lu, J., et al., find it is common to use the definition and numerical formulation of CAD (computer-aided design) shape data to parameterize the structural shape change [
8]. However, similar to dimension optimization, there is a remaking problem that appropriately adjusts the analysis mesh inside the structure under the shape change of the structure boundary.
Conversely, the researchers proposed a topological optimization method based on the reverse variation principle on the edge, since it is a method for finding the optimum boundary shape mathematically to realize the design target under the given load condition and constraint condition [
9,
10]. Although it is applied to structural design or independent part design, it is easy to see that it is difficult to apply to complicated topological optimization problems often found in the design site [
11].
Observations from earlier studies [
12] during the actual machine part topological optimization problem find, as the external boundary shape of the structure is changed, the internal shape of the hole or the uneven partial shape is always adjusted according to requirements such as the assembly relationship and the processing condition between the part and the part. It is common to have a certain partial shape dependency. Alternatively, it is required that a partial shape dependency relationship always exists between the boundary shape outside the structure and the internal partial shape in the course of shape optimization [
13,
14].
To solve this partial shape subordination problem, some researchers have applied the basis vector method based on the basic shape change pattern [
15,
16], like Tang, Y., and Mam, K., et al., who proposed a COBYLA optimization method, and Wang, C., et al., who proposed evaluation of manufacturing constraints from the early stages of the conceptual design to perform an optimization coherent with the manufacturing technology chosen. Actually, when optimizing the structure using the basis vector method, one first considers the structural processing condition of the design object and the assembling relation with other parts, from the designer [
17,
18], then some basic shape change patterns and their change range. Then, each new basic shape change pattern is multiplied by a weighting coefficient and added to form a new structure shape. Each weighting coefficient is used as a design variable as it is, thereby obtaining an optimal weighting factor [
19,
20]. Here, since each basic shape change pattern is expressed in the form of a node coordinate vector having a similar phase configuration, this method is called a basis vector method. To the authors’ knowledge, there are no studies of shape optimization considering the partial shape dependence problem. However, when the basic vector method is applied to complex structural, topological optimization problems, there is no theoretical research on the partial shape dependency problem, and the basic vector problems which depend on the relationship of multiple basic shape change patterns are not solved [
21]. This has a significant impact on the solution of complex structure topological optimization problems. While these simple approaches are adequate for many applications, it is nevertheless clear for practitioners that they can become insufficient to address the partial shape dependency optimization method [
22,
23].
Conversely, many researchers have been making efforts to add optimum design functions in general-purpose analysis codes such as ANSYS and HyperWorks but, due to the complexity and difficulty of partial shape subordination problems, many cases have stopped halfway.
Until now, the partial shape dependency problem in complicated structural shape optimization has not been solved very well, and there is no theoretical research on the partial shape dependency problem. Here, to solve the complicated structural shape optimization problems, some important issues are handled and applied in detail when applying the basis vector method.
First, considering the partial shape dependence problem and the basis vector method, the subordinate relationship of the partial shape inside the structure is divided into three basic cases with the change of the external shape of the structure. The correctness of the basic vector method is verified using mathematical theory, which provides a theoretical basis for the basic vector method used in complex shape optimization design. Additionally, the basic vector generation and mesh regeneration problems in the practical application of the basic vector method are discussed, and the corresponding algorithms are given. The basic vector method proposed in this paper is used to analyze the complex shape optimization design of a two-dimensional turbine blade and a three-dimensional engine base, and the optimization design results verify the local shape dependence problem.
2. Partial Shape Subordination Problem
Here, it will be described using a simple example shown in
Figure 1. The initial shape is shown in
Figure 1a, where we let the design variable
be the width of the rectangular plate, the diameter of the hole in the center is
, the distance from the hole to the left boundary is
, and the distance to the right boundary is
b0. Shown in
Figure 1, when changing the design variable
to
, the shape of the holes inside is considered to be four change patterns from Case1 to Case4.
Case1: The ratio of the lateral length of the hole to the side length is made constant, and the hole changes from a circle to an ellipse.
Case2: Positions of holes which are always fixed regardless of changes in the design variables; the holes are fixed in the center.
Case3: The distance from the hole to the left boundary is made constant, and the hole and the left boundary are fixed.
Case4: The distance from the hole and the hole to the right boundary is made constant, and the hole and the right boundary are fixed.
According to the change of shape design variables described above, there are many subordinate problems during the partial shapes change; this problem actually exists in the shape optimization problem of machine parts, which has a great influence on practical use.
3. Basis Vector Method Used for Shape Optimization
To deal with the partial shape dependency problem mentioned in the previous section, the basis vector method is applied. The basis vector method is a method in which a movement of a certain group of nodes is regarded as a pattern of one shape change, corresponding with one shape design variable for some shape change patterns and adds each weighting factor to obtain a new shape.
Here, the principle of the basis vector method will be explained using the example shown in
Figure 2.
Figure 2a shows the initial shape,
Figure 2b shows the basic shape change pattern 1 in which the entire shape falls forward.
Figure 2c shows the basic shape change pattern 2, in which the upper-end face of the shape is reduced by half in the front-rear direction.
Figure 2d shows the basic shape change pattern 3 in which the upper-end face of the shape is reduced by half in the lateral direction. However, the mesh structure (element connectivity) of the four analysis models shown in
Figure 2 are the same, and only the node coordinate values are different. Here, the node coordinate values of the four analysis models are arranged in the same order, represented by vectors
,
,
,
respectively,
is referred to as an original vector while
,
and
are called basis vectors.
Further, each coordinate vector
,
,
,
is multiplied by a weighting coefficient
x0,
,
,
respectively. Adding the results, the linear combination creates a new node coordinate vector as follows:
When applying Equation (5) to the topological optimization problem in practice, the original vector and the basis vectors , , do not always change while the reference data of shapes change at all times. By changing the design variable , , , , the corresponding node coordinates of the new shape can be calculated directly using Equation (5).
However, when giving the same value to the original vector
and the basis vectors
,
,
, this paper considers the geometric requirement that the new node coordinate vector obtained by Equation (5) becomes the same. The node coordinates of the new shape corresponding to that can be calculated directly using Equation (5):
Substituting Equation (6) into Equation (5), the coordinate calculation formula of the basis vector method used for shape optimization is obtained as follows:
Substituting
,
,
into Equation (7), the analysis model shown in
Figure 3 can be formed using the obtained coordinate vector. Comparing
Figure 1 and
Figure 3, the shape obtained by combining the base shape 1 with 0%, the basis shape 2 at 50% and the basis shape 3 at 50%, according to the value of the design variable, is shown in
Figure 3. Therefore, using the design variables
,
,
, it may be said to be the adoption rate of each basis shape.
Regarding this, since the design variable
,
,
is a continuous variable, it is conceivable that the number of shapes formed by calculation using Equation (7) is infinite. Furthermore, the coordinate calculation of the basis vector method expressed by Equation (7) can be generalized when there are several basis vectors.
However, it is an important research subject, whether it is possible to solve the partial shape subordination problem, as shown in
Figure 1 by using Equation (8).
4. Basis Vector Method and Partial Shape Dependency Relation
Analyzing the partial shape dependency relationship shown in
Figure 1, it can be summarized into the following three basic fundamental relationships:
- (1)
The corresponding relation of a line segment is composed of any two points.
- (2)
There is a proportional relationship of distance between three arbitrary points.
- (3)
The relationship of three arbitrary points is maintained.
Thus, if these three basic subordinations are set, the complex partial shape subordination formed by them also can be ensured.
4.1. Corresponding Relation of a Line Segment Composed of Any Two Points
Considering two arbitrary points a and b inside the structure, and taking direction vectors from point a to point b in the original shape, the basic shape and transformed new shape are as follows:
Considering if the shape transformation on the weighting factor is
ω0 and
ω1, the direction vector after transformation according to Equation (8) is represented by Equation (12),
Taking Equations (9) and (10) into consideration, the above equation becomes Equation (13)
Considering Equation (13), it is found that the direction vector consisting of two arbitrary points inside the structure has a linear combination relation with the same weighting coefficient as the shape transformation formula.
Furthermore, to investigate the change of the line segment
, Equation (14) is obtained by squaring both ends of Equation (13).
Here,
is the angle between the vector
and
. Using this formula, it can calculate the length of the line segment after shape transformation. As a special case, when the basis vector becomes parallel, the above equation becomes Equation (15).
To express the length of the line segment, Equation (16) is obtained.
Also, when
is substituted, Equation (18) is obtained.
Here, the length of the line segment whose shape has been converted indicates that the original line segment length is always maintained, regardless of the weighting factor.
When using the Equation (18) basis vector method for shape optimization, if the complex shape is decomposed into multiple line segments and holes, and convex-concave are set to the same shape, then in the process of shape optimization—even if the combination of any weight coefficients is converted—the shape will remain unchanged.
4.2. Proportional Relationship of Distance between Arbitrary Three Points
Consider the
x coordinates of three arbitrary nodes
a,
b,
c inside the structure as an example. Here, between the three nodes
a,
b,
c in the original shape and each basic shape, it is assumed that there is a proportional relationship concerning the same proportional coefficient, respectively.
These formulas are rewritten to obtain Equation (20),
Also, the x coordinate of the three nodes a, b, and c in the newly transformed shape are as follows:
Substituting Equation (20) into the third equation of Equation (21) and rearranging it, Equation (22) is obtained.
Substituting the first and the second equations of Equation (21) into Equation (22), Equation (23) is obtained.
This expression is consistent with Equation (20). That is, if the same proportional relationship is set between the original shape and the node coordinates in each basic shape, it shows that the same proportional relationship is always maintained regardless of the weighting coefficient between the nodes in the newly transformed shape.
During the actual shape optimization, if some shapes and their proportional relationship are set in the same way using Equation (23), then even if it is transformed into the combination of arbitrary weight coefficients in the process of shape optimization, the proportional relationship will always remain unchanged.
4.3. Holding Relationship of Arbitrary Three Points
During the case of shape optimization accompanying the rotation of the structure, to investigate the change in the angle of the structure, consider the angle formed by the three nodes
a, b, c inside the structure and make the angle between the original shape and the basic shape coincide, as shown in Equation (24).
Here, is the angle between the vector and . Similarly, is the angle between the vector and .
Since the aim is to investigate the influence on the angle due to structural rotation, and it is assumed that the side length constituting the angle does not change:
Also, the rotation angle
from the original shape to the basic shape is taken. The vector that constitutes the new angle
in the transformed shape is as shown in Equation (26).
Then, take an inner product on each side of Equation (26):
Considering Equation (25) and rearranging it, Equation (28) can be obtained:
Conversely, Equation (29) is obtained using the inner product calculation formula of the vector:
Making the first equation of Equation (26) square on both sides, Equation (30) can be obtained:
Considering Equation (29) and rearranging it, Equation (31) is obtained:
Likewise, Equation (32) is obtained from the second equation of Equation (26):
By substituting Equations (31) and (32) into Equation (29), Equation (33) is obtained:
By comparing Equations (28) and (33), Equation (34) is obtained:
The new angle in the shape transformed by Equation (34) always coincides with the angle in the original shape, irrespective of the weighting coefficient.
When actual shape optimization is performed, it is possible to control the change of the angular value composed of three arbitrary points of the structure using Equation (34).
When performing the actual shape optimization using Equation (23), by making the angle value consisting of any three points of the structure in the original shape and in the basis shape identical in a shape optimization way, it can be ensured that the angle value does not ever change.
6. Application Example
6.1. Shape Optimization of Steam Turbine Implantation Section
The movable wing of the steam turbine for thermal power generation equipment is composed of a blade and a rotor. To replace the blade easily, the root of the tree-type blade is directly inserted into the rotor for assembly, as shown in
Figure 5.
When the steam turbine operates, centrifugal force is created due to the rotation of the rotor action and the concentrated stress is generated in the mutually contacting parts of the tree-type implanted part, which connects the blades and the rotating rotor. When there is improper internal stress in parts, it will seriously reduce the fatigue strength of parts; on the contrary, it will increase the life or fatigue strength of parts.
Since most defects where cracks occur are seen partially, it is necessary to reduce the concentrated stress of the steam turbine as much as possible. How to obtain the optimum shape of the tree-type e implantation part is an important issue.
Here, to simplify the analysis problem, as shown in
Figure 5, one blade is cut out from a large number of blades installed along the circumferential direction of the rotating rotor. Further considering the periodic symmetry in the circumferential direction, as shown in the figure, half of the blade implanting portion is made into the FEM model, the boundary constraint condition is cyclically symmetrical in the circumferential direction, the centrifugal force is converted from the mass of the blade effective portion, and the rotational speed of the rotor is applied to the root portion of the blade.
The structural design of the implant portion, upon considering the assembled condition of the parts, etc., as the basic pattern of the shape changes, is shown in
Figure 6.
Changing the slope of the tooth tip
Length of implanted part
Width of implanted part
Width of implanted part
Change in the inclination of the tooth contact surface
The radius of curvature of the first hook
Change in the slope of the first hook
Considering the influence of fatigue stress, the constraint condition is the stress difference between the hooks being less than or equal to the allowable value of the stress in each hook, the stress difference between the hooks is added below the allowable value, and the number of constraint conditions is 27.
However, the seven shape variables shown in
Figure 6 cannot be changed freely, there is a partial shape dependency relationship among them, and some shapes should be the same as far as possible. To obtain simple shapes that are easy to process, the shape changes are shown in
Figure 6.
- (1)
Tooth tip points a, b, c are always on the same straight line
- (2)
The distance between and of a tooth tip point are always equal
- (3)
The curved contour line of each hook is always the same
An optimum shape will be determined on the premise that the requirements relating to these three shape changes are forcibly satisfied.
To realize this problem, following the shape indicated by the basic shape change pattern shown in
Figure 7 and using the smoothing method discussed in the previous section to generate a basis vector, shape optimization can be performed according to the algorithm shown in
Figure 4.
Aiming for the optimization constraint condition, the maximum mises stress value on the blade side and the rotor side in each hook is limited to the allowable amount or less. Furthermore, constraints are added to the stress difference between the hooks so stress is not concentrated on one hook, but dispersed over the whole. The number of constraint conditions is 27 in total.
Figure 7 shows the stress distribution of the optimum shape and initial shape obtained by optimization calculation. Considering
Figure 7, a relatively low distribution is seen in the upper part on the rotor side and the lower part on the blade side in the stress distribution of the initial shape. As a result of averaging the load sharing rate by optimization, the tooth tip inclination angle becomes larger with respect to the vertical plane, and the lateral width of the moving blade increases.
Since the shape and centrifugal force of the blade implant are limited, the average stress value is determined by the ratio of centrifugal force and cross-section, so it is difficult to reduce the comprehensive stress level only by adjusting the internal shape. Therefore, the purpose of the optimization problem is to balance the stress distribution in the implant.
Accompanying the increase of the width and the movement of the material, the stress on the blade side is lower. Since the stress on the rotor side tends to be higher, the angle of the contact surface becomes acute relative to the horizontal plane, and the angle of the tooth surface clamp becomes larger. Additionally, to reduce the stress concentration of the first hook, the stress balance should be maintained to increase the radius of curvature.
Table 1 shows the change in the basis vector coefficient corresponding to each shape change pattern. The change range of each coefficient is from −1 to 1, and the initial value is 0. Regarding the optimum shape, it coincides with the shape shown in
Figure 6 where the coefficient becomes 1, is adopted 100%, and the coefficient becomes minus changes in the direction opposite to the shape shown in
Figure 6.
Figure 8a is a distribution showing the Mises stress value for the initial shape and the optimum shape along with the surface layer of the blade. The vertical axis is the stress value, and the hook corresponding to the horizontal axis is shown in
Figure 3. Similarly,
Figure 8b shows the stress distribution along with the surface layer on the rotor side. Going from there, the concentrated stress of the first hook on the blade side and the fourth hook on the rotor side is reduced, in which stress is high in the initial shape. Additionally, it is found that the maximum stress at each hook is almost the same value, the stress difference of each hook is decreased, and the overall stress balance is improved.
Further, as shown in
Figure 7, it is understood that the assembling dimension of the blade implant portion (the distance between the hook and the hook) is not changed as compared with the original. After optimization, it is easy to confirm that the size of the position of the hook coincides with the original. That is, the position of the vertex of each hook becomes a straight line, and all three requirements for shape change are satisfied.
6.2. Stress Distribution Optimization of the Engine Mount Bracket
Figure 9 shows a model of an engine mount bracket for fixing a passenger car engine. The load conditions occur when applying sudden braking, and the uniform load in the oblique direction is applied along the lower half inner diameter of the upper hole. The value of the uniform load is 6000 N in the axial direction, 23,000 N in the lateral direction, and 45,000 N in the downward direction.
The dotted line in
Figure 10 represents the range in which the structural shape under design can be changed and, in the optimization calculation, these change ranges represent the basis vectors. Furthermore, by adjusting the ratio and combination of each basis vector, the optimum shape can be obtained.
The objective function of the optimization calculation minimizes the weight of the structure. Concerning the constraint condition, the maximum stress of each part is set to the allowable stress value of 140 MPa or less. Optimum solutions are obtained by optimization, and the optimal parameters corresponding to the shape change pattern represented by the dotted line in
Figure 10 are shown in
Table 2. Regarding the table, column 1 is the shape parameter, column 2 is the maximum value of the shape parameter, column 3 is the initial value of the shape parameter, column 4 is the optimum value of the shape parameter, and column 5 is the adoption rate of the shape change pattern.
Considering the stress distribution diagram shown in
Figure 11, in the optimization calculation, each of these change ranges is represented by a basis vector. Concerning the original shape, the maximum Mises stress is 98.4 MPa, and there is a significant difference from the allowable value, which shows that useless material is used more than necessary. Regarding the optimized shape, the maximum Mises stress is 136.8 MPa higher than that of the original shape. Still, the difference from the allowable value is small, the stress in each part is dispersed, and the performance of the material can be fully utilized.
Compared to the original in
Figure 11, concentrated stress is reduced and stress is dispersed. Weight decrease, as the objective function, is 49.9% lighter than the original shape. Additionally, it shows that the requirements for the relative relationship between the external boundary shape and the inner hole are all satisfied and the partial shape subordination problem is solved.
7. Discussion
Based on the results of the examination of the practical example in the previous section, it is shown that it is effective to solve the shape optimization problem having the partial shape subordinate relation using the method proposed in this report. Regarding the case applying to the problem, it also is necessary to examine several problems in detail.
7.1. Creation of Basis Vectors
Prior to shape optimization using the basis vector method, it is necessary to create a basic shape change pattern as a reference for the shape change and a basis vector represented by its node coordinates. Additionally, compared to the original shape, the basic shape is required to be the same as that of the original model, and only the coordinate values of the basis vector and the original vector having the same node number are different from each other.
Here, the paper proposed a basis vector creation method making full use of the finite element method analysis mesh smoothing technique. Specifically, as explained using the simple example shown in
Figure 12, a basis vector can be created according to the following procedure:
- (1)
According to the design scheme of the original shape, the analytical model used to generate the basic vector is divided into three parts: the fixed part (
a1,
a2,,,
a5), the specified part (
b1,
b2,,,
b5), and other shape-related changes as shown in
Figure 12a.
- (2)
A new node coordinate is given to the designated part (
b1,
b2,,,
b5) of the shape change, and the node position of the designated part is updated as shown in
Figure 12b.
- (3)
To enable the shape-related change part, excluding the fixed part (a1,a2,,,a5) whose shape is not changed and the designation part (b1,b2,,,b5) to forcibly change, a smoothing calculation of the analysis mesh is performed using Equation (35).
Considering Equation (35), xi, yi, zi are the new coordinate values of the nodes i of the shape-related change part, xi, yi, zi are the neighboring node i coordinate values adjacent to the node, and j is a surrounding node-set adjacent to nodes i.
- (4)
Smoothing calculation of the mesh, with respect to the repetitive shape-related change part, is performed using Expression (35) and, as a result, a smooth basis shape can be obtained. Here, the basis shape obtained by the number of repetitive smoothing calculations is also different. The basic shape can be finely adjusted according to the actual shape optimization request so the number of repetitive calculations is small in
Figure 12c, and the shape with many repetitive calculations is shown in
Figure 12d.
- (5)
Basis vectors can be created directly by sequentially arranging the node coordinate values of the basic shape obtained last.
When actually generating a basis vector there are two methods; one method is moving a shape designation portion directly, as shown in
Figure 13, and the other is a projection method to a CAD shape in which a shape designation portion is accurately created, as shown in
Figure 14.
Figure 13a shows an original shape in which a circular hole is drilled in the center of a thick square plate, and
Figure 13b is a basic shape change pattern that moves only the shape designation portion around the inner hole in the upward direction.
Figure 14a shows an original shape having a thick-walled cylindrical shape, the inner surface of which is a fixed portion, and the outer surface is a shape designating portion. Shown in
Figure 14b, a CAD shape to be the destination of the shape designation portion is prepared.
Figure 14c shows the outer surface of the cylinder is projected to the CAD shape, at the same time, the analysis mesh of the thick portion of the cylinder is subjected to smoothing processing and, finally, a basic shape change pattern is shown in
Figure 14d.
7.2. Mutual Independence Problem of Basis Vectors
When changing the structure shape using the shape transformation Equation (3), for the combination of arbitrary weighting coefficients , the shape vector obtained by calculation using Equation (7) is considered to be represented by a linear combination of the basis vectors. Here, the mutual relationship between each basis vector is considered.
Figure 15 shows a structure composed of two elements, and its nodes 1, 2, and 3 are defined as shape change nodes. The basis vectors corresponding to the original shape and the basic shape change pattern shown in
Figure 15a–d are expressed by Equations (36)–(39).
Further, the shape change vector represented by the difference between the basic shape change pattern and the basis vector of the original shape is expressed by Equation (37).
Here, to evaluate the correlation between the two shape change vectors
and
, it is necessary to evaluate using the MAC (Modal Assurance Criterion) value calculated using the
value calculation, Equation (38), used for evaluating the mutual relationship of the vibration modes.
Using Equation (38), it is understood that . When shape change vectors and have a linearly independent relationship, it indicates that there is no mutual influence between shape change patterns represented by and . When it indicates that there is a proportional relation between shape change vectors and , and eliminating one in and does not affect the shape optimization result.
Therefore, with respect to the shape change vector of Equation (37), the MAC value between each shape change vector is calculated using Equation (38), and the result becomes Equation (39).
Considering the case of shape optimization by the basis vector method, it is ideal if there is a mutual linear independent relationship between each shape change vector; however, as shown in
Figure 7, when a plurality of basic shape change patterns are set along the same direction with respect to one boundary shape, as shown in Equation (39), a certain degree of change—even though there is a linear relationship since it is possible to express shape change accurately by a combination of shape change vectors—is verified by the results of numerical analysis that the result of shape optimization is hardly affected.
However, if the linear combination between the shape change vectors is potent and the MAC value is too close to 1, then the design sensitivity components corresponding to the design variable may interact with each other. When the corresponding basic shape change modes are adjusted, for example, two basic shape change modes are changed to one. It can greatly improve the analytical efficiency and quality of shape optimization.
7.3. Remeshing Problem of Analysis Model
During the process of shape optimization, with the change of shape, the mesh of the analysis model will change greatly. The aspect ratio of some meshes, or the inner angle of the element edge, will be increased so the mesh quality will have significantly deteriorated, so it may not be able to continue the analysis in the process of shape optimization.
To solve this problem, when each basic shape change pattern is determined and the corresponding basis vector is created, the mesh quality corresponding to each basis vector can be checked using the following Equation (40):
Here and are the initial shape and the basis vector corresponding to the basic shape change pattern, respectively, is the weight coefficient corresponding to the th basic shape change pattern, and n is the number of basis vectors.
A concrete method of checking the analysis mesh will be described using the simple example shown in
Figure 16a–c which are structural shapes and analysis meshes obtained by setting weighting coefficients
to
, a maximum value (
), and a minimum value (
), respectively. Also, the weighting coefficient
may be considered as the adoption rate of the
ith basic shape change pattern, but
and
mean that the respective basic shape change pattern is adopted.
Here, only the analysis mesh of
Figure 16b,c with the largest shape change is checked with respect to the analysis mesh of the original shape and, if there is no problem in the quality of the analysis mesh of these two cases, all of the analysis meshes of the structural shapes formed in the optimization process linearly change between
Figure 16b,c, so there is no problem in the quality of the analysis mesh.
When there is a problem with the quality of the analytical grid in
Figure 16b,c, if the analysis fails, the maximum or minimum value of the weighting coefficient can be adjusted appropriately in the direction close to 0, or the basic shape change mode
i is modified, and the system vector can be recreated according to the situation.