CN102663163B - Topological optimization design method based on geometrical background meshes and with consideration of manufacturing constraints of drawing die - Google Patents
Topological optimization design method based on geometrical background meshes and with consideration of manufacturing constraints of drawing die Download PDFInfo
- Publication number
- CN102663163B CN102663163B CN201210071062.2A CN201210071062A CN102663163B CN 102663163 B CN102663163 B CN 102663163B CN 201210071062 A CN201210071062 A CN 201210071062A CN 102663163 B CN102663163 B CN 102663163B
- Authority
- CN
- China
- Prior art keywords
- finite element
- grid
- parting
- profile
- withdrawing pattern
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Complex Calculations (AREA)
- Moulds For Moulding Plastics Or The Like (AREA)
Abstract
The invention discloses a topological optimization design method based on geometrical background meshes and with consideration of manufacturing constraints of drawing die, which is used for solving technical problems of introducing drawing die constraints into topological optimization design of irregular mesh models. A technical scheme in the invention adopts a geometrical background mesh method, wherein mesh centroid of a finite element is projected to a geometrical background mesh along the direction of split dies or profiles of the drawing die, and pseudo-density of the finite element unit which constrains the centroid in a same geometrical background mesh decreases in order along the direction of split dies or profiles. The method disclosed by the invention effectively overcomes the problem that finite element models which is divided by free meshes can not apply drawing die constrains.
Description
Technical field
The present invention relates to a kind of method of topological optimization design, particularly relate to a kind of based on background geometry grid and consider the method for topological optimization design of withdrawing pattern manufacturing constraints.
Background technology
With reference to Fig. 1~3.In fields such as Aero-Space, automobile makings, there are a large amount of components of machine need to be with Mould Machining manufacture or sand casting, so need to consider withdrawing pattern problem.Mould withdrawing pattern is in order to guarantee mould product demoulding smoothly in the process of producing product; Sand casting withdrawing pattern is not destroy sand mold in order to take out wooden model from sand.The result of traditional Structural Topology Optimization Design, cannot create by the mode of Mould Machining or sand casting.Along with Topology Optimization Method is used widely in engineering, consider that the Topology Optimization Method of the various manufacturing constraints such as withdrawing pattern constraint meets the requirement of structural design more, so need to consider withdrawing pattern problem in topology optimization design.
Hinder several situations of withdrawing pattern: in structure, have cap holes 1, cap holes 1 hinders withdrawing pattern; The angle 2 of the parting/profile 4 of structure outward flange and withdrawing pattern is less than 90 degree, and angle 2 hinders withdrawing pattern.
With reference to Fig. 4.Document " Zhou M; Shyy YK; Thomas HL (2001) Topology optimization with manufacturing constraints.In:4th world congress of structural and multidisciplinary optimization, Dalian " discloses a kind of method of topological optimization design of considering withdrawing pattern manufacturing constraints.Document is introduced withdrawing pattern manufacturing constraints condition, and, in withdrawing pattern direction, the pseudo-density of topology optimization design variable meets constraint:
Wherein,
pseudo-density for the unit in same row 3 in withdrawing pattern direction.
Although the disclosed method of document is incorporated into withdrawing pattern manufacturing constraints in topology optimization design, the method require the finite element grid of design section must be in withdrawing pattern direction complete matching, finite element grid is necessary for the regular grid that size is identical.Shortcoming is: the method can only be applied to simple structure; For the structure of dividing with irregular finite element grid, because its finite element grid can not complete matching in withdrawing pattern direction, therefore can not apply the method.
Summary of the invention
In order to solve the technical matters of introducing withdrawing pattern manufacturing constraints in the topology optimization design of irregular grid model, the invention provides a kind of based on background geometry grid and consider the method for topological optimization design of withdrawing pattern manufacturing constraints.The method is by definition background geometry grid, the barycenter of finite element grid is mapped in background geometry grid, cell density by constraint barycenter in same background geometry grid diminishes successively along parting/profile, can realize and in the topology optimization design of irregular finite element grid model, introduce withdrawing pattern constraint.
The technical solution adopted for the present invention to solve the technical problems is: a kind of based on background geometry grid and consider the method for topological optimization design of withdrawing pattern manufacturing constraints, be characterized in comprising the following steps:
(a) by the cad model of structure, set up finite element model, definition load and boundary condition.
(b) parting/profile of withdrawing pattern 4 is divided into the background geometry grid 7 measure-alike with finite element grid 6.For any one the finite element grid i on design space, its barycenter is necessarily under the jurisdiction of some background geometry grid m along the projection of withdrawing pattern parting/profile 4 normal direction, m ∈ (1,2 ..., 27).
(c) set up Topological optimization model:
find X=(x
1,x
2,K,x
n)
min Φ(X)
s.t.KU=F
(0<x
p≤x
q...≤x
w≤1)
m,m=1,2,K,M
(d
p≥d
q...≥d
w)
m,m=1,2,K,M
In formula, X is the pseudo-intensity vector in the unit in design domain; N is design variable number; The objective function that Φ (X) is topological optimization; K is finite element model global stiffness matrix; F is node equivalent load vectors; U is node global displacement vector; G
j(X) be j constraint function;
it is the upper limit of j constraint function; J is the quantity of constraint; x
p, x
q... x
wfor barycenter drops on the pseudo-density of the finite element grid in background geometry grid m along the projection of withdrawing pattern parting/profile 4 normal direction; M is the total columns in withdrawing pattern direction.D
ifor unit x
ithe distance of centroid distance parting/profile 4, i.e. d
p, d
q, d
wbe respectively unit x
p, x
q, x
wdistance apart from parting/profile 4.For barycenter, along the projection of parting/profile 4 normal direction of withdrawing pattern, drop on the unit x of same background geometry grid
p, x
q... x
w, its barycenter is larger apart from the distance of parting/profile 4, and the pseudo-density of its unit is less.
(d) by finite element soft Ansys, model is carried out to a finite element analysis; By structure optimization platform Boss-Quattro, be optimized sensitivity analysis again, try to achieve the sensitivity of objective function and constraint condition, choose gradient optimal method GCMMA and be optimized design, result is optimized.
The dimension of described background geometry grid 7 is than the low one dimension of the dimension of finite element grid 6.
The invention has the beneficial effects as follows: owing to having adopted the method for background geometry grid in topological optimization, finite element grid barycenter is mapped on background geometry grid along withdrawing pattern parting/profile, the pseudo-density of the finite element unit of constraint barycenter in same background geometry grid reduces successively along parting/profile direction, and the finite element model that has effectively overcome free grid division can not be applied the problem of withdrawing pattern constraint.
Below in conjunction with drawings and Examples, the present invention is elaborated.
Accompanying drawing explanation
Fig. 1 is the structural representation that can manufacture by Mould Machining or sand casting in background technology.
Fig. 2 is one of the structure that cannot manufacture by Mould Machining or sand casting in background technology schematic diagram.
Fig. 3 is two schematic diagram that cannot pass through the structure of Mould Machining or sand casting manufacture in background technology.
Fig. 4 is the schematic diagram that the disclosed method of background technology document is introduced withdrawing pattern constraint.
Fig. 5 is the stress model of the inventive method.
Fig. 6 is finite element grid and the background grid schematic diagram of the inventive method.
Fig. 7 is the topology optimization design result that the inventive method is introduced withdrawing pattern constraint.
In figure, 1-hole, 2-angle, 3-is the unit in same row in withdrawing pattern direction, 4-parting/profile, 5-border, 6-finite element grid, 7-background geometry grid.
Embodiment
With reference to Fig. 5~7.Take two-dimentional semi-girder topology optimization design as example explanation the present invention.Two dimension semi-girder is of a size of long 500mm, high 100mm.Semi-girder left side is complete fixing; Base angle, right side is subject to concentrated force load F=100N vertically upward.Design cantilever beam structure, makes its rigidity maximum, and volume fraction is 60%.Withdrawing pattern direction is for vertically upward.Method concrete steps are as follows:
(a) finite element modeling.
Cad model by structure is set up finite element model, definition load and boundary condition.It is that base angle, right side is subject to concentrated force load F=100N vertically upward that model is subject to load; The boundary condition of model is that semi-girder left side is all fixing.
(b) divide background geometry grid.
Parting/the profile of withdrawing pattern 4 is divided into the background geometry grid 7 measure-alike with finite element grid 6.It can be both the division that meets certain rule that background geometry grid 7 is divided, and can be also to divide arbitrarily.Because the parting/profile 4 of a three-dimensional model is a two dimensional surface, so the background geometry grid 7 of three-dimensional finite element model is two dimension; The background geometry grid 7 of corresponding two-dimensional finite element model is one dimension.Be that the dimension of background geometry grid 7 is than the low one dimension of the dimension of model finite element grid 6.Background geometry grid 7 sizes are suitable with finite element unit size, suppose total M background geometry grid 7, and the present embodiment has been divided 27 background geometry grids 7 altogether.For any one the finite element grid i on design space, its barycenter is necessarily under the jurisdiction of some background geometry grid m along the projection of withdrawing pattern parting/profile 4 normal direction, (m ∈ (1,2 ..., 27)).
(c) setting up Topological optimization model is:
find X=(x
1,x
2,K,x
n)
minΦ(X)
s.t.KU=F
(0<x
p≤x
q...≤x
w≤1)
m,m=1,2,K,M
(d
p≥d
q...≥d
w)
m,m=1,2,K,M
Wherein, X is the pseudo-intensity vector in unit in design domain; N is design variable number; The objective function that Φ (X) is topological optimization; K is finite element model global stiffness matrix; F is node equivalent load vectors; U is node global displacement vector; G
j(X) be j constraint function;
it is the upper limit of j constraint function; J is the quantity of constraint; x
p, x
q... x
wfor barycenter drops on the pseudo-density of the finite element grid 6 (unit) in background geometry grid m along the projection of withdrawing pattern parting/profile 4 normal direction, wherein the border of model geometric background grid is border 5; M is the total columns in withdrawing pattern direction, the i.e. sum of background geometry grid 7.D
ifor unit x
ithe distance of centroid distance parting/profile 4, i.e. d
p, d
q, d
wbe respectively unit x
p, x
q, x
wdistance apart from parting/profile 4.For barycenter, along the projection of parting/profile 4 normal direction of withdrawing pattern, drop on the unit x of same background geometry grid 7
p, x
q... x
w, its barycenter is larger apart from the distance of parting/profile 4, and the pseudo-density of unit is less.
The Topological optimization model that the present embodiment is set up is:
find X=(x
1,x
2,K,x
n)
min C(X)
s.t. KU=F
V(X)-0.6≤0
(0<x
p≤x
q...≤x
w≤1)
m,m=1,2,K,27
(d
p≥d
q...≥d
w)
m,m=1,2,K,27
(d) finite element analysis and Optimization Solution.
Model is carried out to a finite element analysis; By optimizing sensitivity analysis, try to achieve the sensitivity of objective function and constraint condition, choose certain optimized algorithm and be optimized design, result is optimized.
By finite element soft Ansys, model is carried out to a finite element analysis; By structure optimization platform Boss-Quattro, be optimized sensitivity analysis again, try to achieve the sensitivity of objective function and constraint condition, choose gradient optimal method GCMMA (Globally Convergent Method of Moving Asymptotes) and be optimized design, result is optimized.
By optimum results, can be found out: by introducing background geometry gridding technique, can apply withdrawing pattern manufacturing constraints to the finite element model of any cell type, rather than only be confined to the finite element model of regular grid.
Claims (1)
1. based on background geometry grid and consider the method for topological optimization design of withdrawing pattern manufacturing constraints, it is characterized in that comprising the following steps:
(a) by the cad model of structure, set up finite element model, definition load and boundary condition;
(b) parting/profile of withdrawing pattern (4) is divided into the background geometry grid (7) measure-alike with finite element grid (6); For any one the finite element grid i on design space, its barycenter is necessarily under the jurisdiction of some background geometry grid m along the projection of parting/profile (4) normal direction of withdrawing pattern, m ∈ (1,2 ..., 27);
(c) set up Topological optimization model:
find X=(x
1,x
2,...,x
n)
min Φ(X)
s.t.KU=F
(0<x
p≤x
q...≤x
w≤1)
m,m=1,2,...,M
(d
p≥d
q...≥d
w)
m,m=1,2,...,M
In formula, X is the pseudo-intensity vector in the unit in design domain; N is design variable number; The objective function that Φ (X) is topological optimization; K is finite element model global stiffness matrix; F is node equivalent load vectors; U is node global displacement vector; G
j(X) be j constraint function;
it is the upper limit of j constraint function; J is the quantity of constraint; x
p, x
q... x
wfor barycenter drops on the pseudo-density of the finite element grid in background geometry grid m along the projection of parting/profile (4) normal direction of withdrawing pattern; M is the total columns in withdrawing pattern direction; d
ifor unit x
ithe distance of centroid distance parting/profile (4), i.e. d
p, d
q, d
wbe respectively unit x
p, x
q, x
wdistance apart from parting/profile (4); For barycenter, along the projection of parting/profile (4) normal direction of withdrawing pattern, drop on the unit x of same background geometry grid
p, x
q... x
w, its barycenter is larger apart from the distance of parting/profile (4), and the pseudo-density of its unit is less;
(d) by finite element soft Ansys, model is carried out to a finite element analysis; By structure optimization platform Boss-Quattro, be optimized sensitivity analysis again, try to achieve the sensitivity of objective function and constraint condition, choose gradient optimal method GCMMA and be optimized design, result is optimized.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210071062.2A CN102663163B (en) | 2012-03-17 | 2012-03-17 | Topological optimization design method based on geometrical background meshes and with consideration of manufacturing constraints of drawing die |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210071062.2A CN102663163B (en) | 2012-03-17 | 2012-03-17 | Topological optimization design method based on geometrical background meshes and with consideration of manufacturing constraints of drawing die |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102663163A CN102663163A (en) | 2012-09-12 |
CN102663163B true CN102663163B (en) | 2014-11-05 |
Family
ID=46772654
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210071062.2A Expired - Fee Related CN102663163B (en) | 2012-03-17 | 2012-03-17 | Topological optimization design method based on geometrical background meshes and with consideration of manufacturing constraints of drawing die |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102663163B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103207940B (en) * | 2013-04-23 | 2015-07-22 | 西北工业大学 | Topology optimization design method for cyclic symmetry cylindrical grid structure |
CN103559740B (en) * | 2013-08-26 | 2016-06-29 | 空气动力学国家重点实验室 | A kind of prism mess generation method realizing interactive operation |
CN109002598B (en) * | 2018-06-29 | 2020-09-18 | 华中科技大学 | Self-supporting microstructure topology optimization method considering overhanging angle and minimum size constraint |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339575A (en) * | 2008-08-07 | 2009-01-07 | 上海交通大学 | Three-dimensional visualized process design system and its design method |
CN101609479A (en) * | 2009-06-23 | 2009-12-23 | 北京理工大学 | A kind of optimal design method for trajectory robust |
-
2012
- 2012-03-17 CN CN201210071062.2A patent/CN102663163B/en not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339575A (en) * | 2008-08-07 | 2009-01-07 | 上海交通大学 | Three-dimensional visualized process design system and its design method |
CN101609479A (en) * | 2009-06-23 | 2009-12-23 | 北京理工大学 | A kind of optimal design method for trajectory robust |
Non-Patent Citations (2)
Title |
---|
Topology optimization with manufacturering constraints;zhou M et al.;《Proceeding of the Fourth World Congress of structural and multidisciplinary optimization》;20010630;参见摘要,正文第4节、图1 * |
zhou M et al..Topology optimization with manufacturering constraints.《Proceeding of the Fourth World Congress of structural and multidisciplinary optimization》.2001, * |
Also Published As
Publication number | Publication date |
---|---|
CN102663163A (en) | 2012-09-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102622486B (en) | Topological optimization design method taking drafting mold manufacturing constraint into consideration | |
Xia et al. | A level set based method for the optimization of cast part | |
Qian | Full analytical sensitivities in NURBS based isogeometric shape optimization | |
CN102842147B (en) | Three-dimensional integral wall panel expanding and modeling method | |
CN100557614C (en) | A kind of method that generates three-dimensional finite element mesh | |
CN110555263B (en) | Level set topology optimization method for curved shell structure optimization design | |
CN102663163B (en) | Topological optimization design method based on geometrical background meshes and with consideration of manufacturing constraints of drawing die | |
CN104156546B (en) | The shape face redesign method of the car panel die based on T battens | |
CN110929443B (en) | Two-dimensional flood simulation method based on high-precision terrain generalization | |
CN104441488A (en) | Predeformation method and system of unequal-wall-thickness plastic part mold parting surface | |
CN102637223B (en) | Topological optimization design method on basis of pseudo-density sequencing and with consideration of draft manufacturing constraints | |
CN113505408B (en) | Empty box retaining wall full-parametric three-dimensional model construction method based on feature model | |
CN102968518A (en) | Topology optimization design method for novel down-the-hole arc-shaped hydraulic steel gate | |
Tang et al. | Lattice-skin structures design with orientation optimization | |
Chakraborty et al. | Automatic determination of parting directions, parting lines and surfaces for two-piece permanent molds | |
CN103294861A (en) | Multi-arc curve shape optimization design method | |
CN103942365A (en) | Associated design method of aircraft assembly jig framework | |
CN113673186A (en) | Cartesian grid rapid generation method based on STL file | |
CN107515960B (en) | Feature modeling based topological optimization design method for circularly symmetric cylindrical support structure | |
JPH08315183A (en) | Method and system for automatic mesh generation | |
CN104573846B (en) | A kind of polymorphism job shop layout optimization method based on CA PSO hybrid optimization algorithms | |
CN101169645A (en) | Bending-torsional non-linear variable-diameter tube pressed mould manufacture method | |
Brito-Parada et al. | The advantages of using mesh adaptivity when modelling the drainage of liquid in froths | |
Wei et al. | Rapid development technique for drip irrigation emitters | |
CN107577894B (en) | Real-time binary encryption-decryption method for quadrilateral unit |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141105 Termination date: 20150317 |
|
EXPY | Termination of patent right or utility model |