CN103455674A - Geometrical reconstruction method of topological optimization result - Google Patents
Geometrical reconstruction method of topological optimization result Download PDFInfo
- Publication number
- CN103455674A CN103455674A CN2013103910235A CN201310391023A CN103455674A CN 103455674 A CN103455674 A CN 103455674A CN 2013103910235 A CN2013103910235 A CN 2013103910235A CN 201310391023 A CN201310391023 A CN 201310391023A CN 103455674 A CN103455674 A CN 103455674A
- Authority
- CN
- China
- Prior art keywords
- broken line
- contour
- contour curve
- point
- line
- 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
Images
Landscapes
- Processing Or Creating Images (AREA)
Abstract
The invention provides a geometrical reconstruction method of a topological optimization result. The method comprises the following steps of: firstly, acquiring the surface geometrical structure IGES (Initial Graphics Exchange Specification) file and data format of the topological optimization result; secondly, ordering contour lines to obtain a parameterized image fitted by a straight line and a circular arc; lastly, outputting an enclosed profile line fitted by a long broken line and a circular arc to CAD (Computer-Aided Design) software to obtain a corresponding modified geometrical model. According to the method, a complex profile image is fitted by the simple segmented straight line and the circular arc, and fitting of different complexity degrees is realized by controlling a plurality of adjustable parameters. The method has the advantages of high interactivity, simple algorithm, effect joint of CAD and CAE (Computer Aided Engineering) software and rapid realization of reconstruction and ordering of a CAD model.
Description
Technical field
The present invention relates to a kind of geometry reconstruction method of topological optimization result.
Background technology
Structural Topology Optimization can, in starting stage of Structural Design for the deviser provides an alternate conceptual design, make structure adopt optimal case on layout.But commercial CAD and CAE software interface are also not too perfect at present, because topological optimization is used finite element model, optimum results can not be directly changed into the model of Parametric CAD, need the designer again to create a Parametric CAD model that approaches the topological optimization result, repeating to import CAE software verifies and further optimizes, simultaneously, the profile as a result of topological optimization is usually also more complicated, is not easy to processing and manufacturing.Be very loaded down with trivial details and consuming time to topological optimization reconstruct as a result and finishing, increased to a certain extent design cycle and design cost.
Domestic to how realizing automatic reconfiguration, research is not deeply.In some Primary Study, proposed several method for transformation at home and abroad, wherein a kind of is that the figure of result after topological optimization is created to the B-spline surface cross section by identification figure gray scale, forms the cad model of 3D by a plurality of B battens cross section; Also having a kind of is that the complexity point cloud file of optimum results is simply put to the cloud file by one of the way formation that reduces grid node, retain its elementary contour information, utilize some cloud file to form cad model, above two kinds of methods have obtained success to a certain extent, but a lot of simple graphs do not need to go to describe with complicated B batten, and shortage logic function operation when the refined model structure, these two kinds of methods of while lack the interactivity with the designer, can not guarantee that result is exactly that the deviser is needed; The Shane Larsen and C.Greg Jensen of Brigham Young university has proposed the shape template method, and some cross section profile geometric coordinate of structure is changed into to polar form, and contrasts with the polar coordinates of template, approaches structure outline.The method allows the complexity of deviser's Definition Model matching, freer and more meet deviser's demand, but weak point is that the template complexity is various, and travelling speed is slow.
Summary of the invention
The technical matters solved
In order to overcome the deficiencies in the prior art, the invention provides a kind of geometry reconstruction method of topological optimization result.
Technical scheme
Technical scheme of the present invention is:
The geometry reconstruction method of described a kind of topological optimization result is characterized in that: adopt following steps:
Step 1: surface geometry IGES file and the data of obtaining the topological optimization result transform:
Step 1.1: by CAE software, generate topological optimization IGES file as a result;
Step 1.2: how much of the elements of IGES file including and positional information are converted into to the discernible data matrix of MATLAB, thereby obtain the outline line geometric parameter;
Step 2: outline line is carried out regular, obtain the parametrization image by straight line and circular fitting:
Step 2.1: contour segmentation: to forming any contour curve AB of closed contour, adopt following methods to carry out segmentation: to connect contour curve and put A point and B point the whole story, the A point and the B point coordinate that according to step 1, obtain, obtain straight line AB equation, and the upper each point of contour curve AB is to the distance of straight line AB, wherein most distant point is C, and maximum distance is CH; If judgement CH/AB>α sets up, mean that contour curve AB can divide, the waypoint using a C as curve A B, contour curve AB is divided into two sections of contour curve AC and contour curve CB; Again to contour curve AC and contour curve CB segmentation, until contour curve is inseparable, inseparable contour curve is put to straight line the whole story and be connected to form broken line; Wherein α is fragmentation threshold; Obtain the new outline line that some sections broken line endings are connected after the closed contour segmentation;
Step 2.2: profile parameters matching: according to the cross-sectional data obtained in step 1, the overall length that obtains closed contour is L
0; Broken line is distinguished: wherein n section broken line length is L
nif, L
n/ L
0<β, n section broken line is short broken line, otherwise is long broken line, β is correction threshold;
For any short broken line PQ, matching obtains the circular arc PMQ of too short broken line two-end-point P point and Q point and upper 1 M of contour curve PQ, with circular arc PMQ, replaces short broken line PQ; If the distance of center circle of the matching circular arc of adjacent many short broken lines is from being less than error threshold, and the difference of the radius of the matching circular arc of adjacent many short broken lines is less than error threshold, matching must be the circular arc of three end points in these adjacent many short broken lines, with this circular arc, replaces this adjacent many short broken lines; Whole closed contour is formed by long broken line and circular fitting;
Step 3: the closed contour that will be formed by long broken line and circular fitting exports in CAD software and obtains corresponding amended geometric model.
Beneficial effect
The contour images of simple segmentation straight line and circular fitting complexity for the present invention, realize the matching of differing complexity by controlling a plurality of adjustable parameters.Interactivity is good, algorithm is succinct, effectively is connected CAD and CAE software, realizes fast reconstruct and the finishing of cad model.
The accompanying drawing explanation
Fig. 1. the technical scheme process flow diagram that the present invention proposes;
Fig. 2. original contour line curve;
Fig. 3. the curve after segmentation;
Fig. 4. the three points arc matching;
Fig. 5. the example ELEMENT MESH GRAPH;
Fig. 6. the planimetric map of deriving by the finite element file;
Fig. 7. the reconstruction result figure of different threshold values.
Embodiment
Below in conjunction with drawings and Examples, the present invention is further described.
In the present embodiment, the finite element grid file of the topological optimization result of a certain planar part, as Fig. 5, has a plurality of irregular holes on plane, and appearance profile is more complicated also.The planar graph of by topological optimization software, deriving the IGES form as shown in Figure 6, has eight closed curves, seven inner outlines wherein, and an outer contour, every outline line is joined and is formed by multistage B-spline curves first place.
The step of the topological optimization result being carried out to the geometry reconstruction method is:
Step 1: surface geometry IGES file and the data of obtaining the topological optimization result transform:
Step 1.1: by CAE software, generate topological optimization IGES file as a result;
Step 1.2: how much of the elements of IGES file including and positional information are converted into to the discernible data matrix of MATLAB, thereby obtain the outline line geometric parameter; Store with a cell array data mode in MATLAB as the B-spline curves in nurbs surface, comprise the geometric parameter information such as its curve type and each point coordinate.
Step 2: outline line is carried out regular, obtain the parametrization image by straight line and circular fitting:
Step 2.1: contour segmentation:
Closed contour is joined end to end and is formed by many curves, to forming any contour curve AB of closed contour, as shown in Figure 2, adopt following methods to carry out segmentation: to connect contour curve and put A point and B point the whole story, the A point and the B point coordinate that according to step 1, obtain, obtain straight line AB equation, and the upper each point of contour curve AB is to the distance of straight line AB, wherein most distant point is C, and maximum distance is CH; If judgement CH/AB>α sets up, mean that contour curve AB can divide, the waypoint using a C as curve A B, contour curve AB is divided into two sections of contour curve AC and contour curve CB; To contour curve AC and contour curve CB segmentation, until contour curve is inseparable, as shown in Figure 3, former contour curve AB is divided into curve A C, CD, DE, EF, FB again.Travel through all contour curves, complete the segmentation of whole profile.Inseparable contour curve is put to straight line the whole story and be connected to form broken line; Obtain the new outline line that some sections broken line endings are connected after the closed contour segmentation.Wherein α is fragmentation threshold, and this value is less, and broken line is more, to portraying also more accurately of profile.
Step 2.2: profile parameters matching:
If only larger with the segmented fitting curve error, introduced thus the circular fitting of outline line, when guaranteeing pattern precision, avoid figure too complicated.
At first the outline line of segmentation screened, according to the cross-sectional data obtained in step 1, the overall length that obtains closed contour is L
0; Broken line is distinguished: wherein n section broken line length is L
nif, L
n/ L
0<β, n section broken line is short broken line, otherwise, be long broken line, β is correction threshold, is the value that can revise according to actual geometric model, and this value is larger, and the zone that is identified as short broken line is more.
For any short broken line PQ, as Fig. 4, matching obtains the circular arc PMQ of too short broken line two-end-point P point and Q point and upper 1 M of contour curve PQ, with circular arc PMQ, replaces short broken line PQ; If the distance of center circle of the matching circular arc of adjacent many short broken lines is from being less than error threshold, and the difference of the radius of the matching circular arc of adjacent many short broken lines is less than error threshold, matching must be the circular arc of three end points in these adjacent many short broken lines, with this circular arc, replaces this adjacent many short broken lines; Whole closed contour is formed by long broken line and circular fitting;
For long broken line, assert that the contour curve that this broken line is corresponding can be replaced by straight-line segment, now this section contour curve only needs the geological information of record start point and terminating point, generates straight-line equation.
Through above step process, whole geometric profile is formed by straight-line segment and circular fitting, original profile B-spline curves surface information is converted into straight line and circular parameters equation information, greatly reduces the complexity of data, for the automatic reconfiguration of cad model provides Information base.
Step 3: but by the data format file of parameter read-in CAD software identifying processing: the closed contour that will be formed by long broken line and circular fitting exports in CAD software and obtains corresponding amended geometric model.
Commercial CAD software generally has its corresponding identified documentation form, during parameterized model output, and be by corresponding form by parameter read-in.Can write statement by the file of MATLAB according to corresponding language format, straight line and circular parameters are write to text with order, generate revisable geometric model by corresponding input method in commercial CAD software.
Business software UG take in the present embodiment as example, write statement according to the Grip language format by MATLAB, straight line after matching and circular parameters are write to text with order, generate the grs file, after this document is compiled by the NX Open Grip in the UG instrument, can directly import in UG and generate revisable cad model, as shown in Figure 7.In the process of contour segmentation and matching, α, β are adjustable values, α is less, broken line is more, β is larger, and the zone that is identified as short broken line is more, also more accurate to the matching of profile, but processability may be poorer, the deviser can regulate parameter to obtain expected result according to actual conditions.Fig. 7 (a) and (b), (c) are respectively different α, geometry reconstruction model corresponding to β value.After geometry reconstruction, the deviser can also continue the element in cad model is adjusted on demand easily.
Claims (1)
1. the geometry reconstruction method of a topological optimization result is characterized in that: adopt following steps:
Step 1: surface geometry IGES file and the data of obtaining the topological optimization result transform:
Step 1.1: by CAE software, generate topological optimization IGES file as a result;
Step 1.2: how much of the elements of IGES file including and positional information are converted into to the discernible data matrix of MATLAB, thereby obtain the outline line geometric parameter;
Step 2: outline line is carried out regular, obtain the parametrization image by straight line and circular fitting:
Step 2.1: contour segmentation: to forming any contour curve AB of closed contour, adopt following methods to carry out segmentation: to connect contour curve and put A point and B point the whole story, the A point and the B point coordinate that according to step 1, obtain, obtain straight line AB equation, and the upper each point of contour curve AB is to the distance of straight line AB, wherein most distant point is C, and maximum distance is CH; If judgement CH/AB>α sets up, mean that contour curve AB can divide, the waypoint using a C as curve A B, contour curve AB is divided into two sections of contour curve AC and contour curve CB; Again to contour curve AC and contour curve CB segmentation, until contour curve is inseparable, inseparable contour curve is put to straight line the whole story and be connected to form broken line; Wherein α is fragmentation threshold; Obtain the new outline line that some sections broken line endings are connected after the closed contour segmentation;
Step 2.2: profile parameters matching: according to the cross-sectional data obtained in step 1, the overall length that obtains closed contour is L
0; Broken line is distinguished: wherein n section broken line length is L
nif, L
n/ L
0<β, n section broken line is short broken line, otherwise is long broken line, β is correction threshold;
For any short broken line PQ, matching obtains the circular arc PMQ of too short broken line two-end-point P point and Q point and upper 1 M of contour curve PQ, with circular arc PMQ, replaces short broken line PQ; If the distance of center circle of the matching circular arc of adjacent many short broken lines is from being less than error threshold, and the difference of the radius of the matching circular arc of adjacent many short broken lines is less than error threshold, matching must be the circular arc of three end points in these adjacent many short broken lines, with this circular arc, replaces this adjacent many short broken lines; Whole closed contour is formed by long broken line and circular fitting;
Step 3: the closed contour that will be formed by long broken line and circular fitting exports in CAD software and obtains corresponding amended geometric model.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310391023.5A CN103455674B (en) | 2013-08-31 | A kind of geometry reconstruction method of topological optimization result |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310391023.5A CN103455674B (en) | 2013-08-31 | A kind of geometry reconstruction method of topological optimization result |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103455674A true CN103455674A (en) | 2013-12-18 |
| CN103455674B CN103455674B (en) | 2016-11-30 |
Family
ID=
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104492081A (en) * | 2014-12-31 | 2015-04-08 | 苏州乐米信息科技有限公司 | Method for realizing curve track in mobile games and curve track equation setting device |
| CN106023283A (en) * | 2016-05-25 | 2016-10-12 | 山东省科学院海洋仪器仪表研究所 | Drawing method of smooth coast line |
| CN107194885A (en) * | 2017-05-02 | 2017-09-22 | 拓卡奔马机电科技有限公司 | A kind of CAD diagram Intelligent Portable Correcting System and method |
| CN109029421A (en) * | 2018-06-22 | 2018-12-18 | 南京擎天科技有限公司 | A kind of regional location coordinate acquisition method based on LBS technology |
| CN110276095A (en) * | 2019-05-05 | 2019-09-24 | 广东省汉兴科技有限公司 | A kind of part model contour line form adaptive decomposition method and device |
| CN110688704A (en) * | 2019-10-17 | 2020-01-14 | 广东三维家信息科技有限公司 | Home decoration design method and system and electronic equipment |
| CN112104515A (en) * | 2020-11-19 | 2020-12-18 | 中国人民解放军国防科技大学 | Network reconstruction method and device, computer equipment and storage medium |
| CN115004188A (en) * | 2019-08-28 | 2022-09-02 | 西门子工业软件有限公司 | System and method for processing topology optimized geometries |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102308298A (en) * | 2008-08-18 | 2012-01-04 | 空客运营有限公司 | Computer-assisted method for the advanced design of folded pieces of composite material |
| CN102955874A (en) * | 2011-08-11 | 2013-03-06 | 利弗莫尔软件技术公司 | Methods and systems for designing addendum section of die in sheet metal forming |
| JP2013114464A (en) * | 2011-11-29 | 2013-06-10 | Asahi Kasei Chemicals Corp | Rib design method, rib design device and rib design program |
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102308298A (en) * | 2008-08-18 | 2012-01-04 | 空客运营有限公司 | Computer-assisted method for the advanced design of folded pieces of composite material |
| CN102955874A (en) * | 2011-08-11 | 2013-03-06 | 利弗莫尔软件技术公司 | Methods and systems for designing addendum section of die in sheet metal forming |
| JP2013114464A (en) * | 2011-11-29 | 2013-06-10 | Asahi Kasei Chemicals Corp | Rib design method, rib design device and rib design program |
Non-Patent Citations (2)
| Title |
|---|
| 葛文杰，黄杰，杨方: "拓扑优化技术及其在汽车设计中的应用", 《机床与液压》 * |
| 黄杰，葛文杰，杨方: "实现机翼前缘形状连续变化柔性机构的拓扑优化", 《航空学报》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104492081A (en) * | 2014-12-31 | 2015-04-08 | 苏州乐米信息科技有限公司 | Method for realizing curve track in mobile games and curve track equation setting device |
| CN106023283A (en) * | 2016-05-25 | 2016-10-12 | 山东省科学院海洋仪器仪表研究所 | Drawing method of smooth coast line |
| CN106023283B (en) * | 2016-05-25 | 2018-12-25 | 山东省科学院海洋仪器仪表研究所 | A kind of method for drafting in smooth coastline |
| CN107194885A (en) * | 2017-05-02 | 2017-09-22 | 拓卡奔马机电科技有限公司 | A kind of CAD diagram Intelligent Portable Correcting System and method |
| CN107194885B (en) * | 2017-05-02 | 2020-07-14 | 拓卡奔马机电科技有限公司 | Intelligent CAD (computer-aided design) drawing correction system and method |
| CN109029421A (en) * | 2018-06-22 | 2018-12-18 | 南京擎天科技有限公司 | A kind of regional location coordinate acquisition method based on LBS technology |
| CN110276095A (en) * | 2019-05-05 | 2019-09-24 | 广东省汉兴科技有限公司 | A kind of part model contour line form adaptive decomposition method and device |
| CN110276095B (en) * | 2019-05-05 | 2023-03-24 | 广东省汉兴科技有限公司 | Workpiece model contour line shape self-adaptive decomposition method and device |
| CN115004188A (en) * | 2019-08-28 | 2022-09-02 | 西门子工业软件有限公司 | System and method for processing topology optimized geometries |
| CN110688704A (en) * | 2019-10-17 | 2020-01-14 | 广东三维家信息科技有限公司 | Home decoration design method and system and electronic equipment |
| CN112104515A (en) * | 2020-11-19 | 2020-12-18 | 中国人民解放军国防科技大学 | Network reconstruction method and device, computer equipment and storage medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105787226A (en) | Reconstruction of parameterization model of quadrilateral finite element grid model | |
| CN106373184B (en) | A kind of 3 D-printing model puts required amount of support Method of fast estimating | |
| CN105426628B (en) | A kind of complex product isomery Model Fusion method | |
| JP2005038219A (en) | Method for generating volume data from boundary expression data and program therefor | |
| CN104156546B (en) | The shape face redesign method of the car panel die based on T battens | |
| CN110688792B (en) | Design method of flat-bottom symmetrical Y-shaped steel bifurcated pipe based on CATIA and ANSYS | |
| CN102682172B (en) | Numerous-parameter optimization design method based on parameter classification for supercritical aerofoil | |
| CN111709092A (en) | Radial tire tread pattern finite element automatic modeling method based on finite element unit nodes | |
| CN105118091A (en) | Method and system for constructing multi-precision non-uniform geological grid curved surface model | |
| CN112926207B (en) | Method and system for automatically constructing editable model of equal-geometry topological optimization result | |
| CN106874610A (en) | A kind of sloping core earth and rockfill dam modeling method based on VB and CATIA | |
| CN103544343A (en) | Serialized and parameterized modeling method of power chuck | |
| CN102201017A (en) | Optimal design method for automotive body | |
| CN103177167A (en) | Cupboard design method | |
| CN104361185A (en) | Wire distribution space automatic generating method for cable virtual design | |
| CN105354881A (en) | Mesh distortion optimization algorithm based on discrete attribute data | |
| CN103310050A (en) | Method for defining curve surface features by user based on existing models | |
| CN110795835A (en) | Three-dimensional process model reverse generation method based on automatic synchronous modeling | |
| CN114241087B (en) | Building plane layout generation method based on bubble diagrams | |
| CN101546351B (en) | Geometric parameterization modeling method for optimizing variable-complexity shape | |
| CN111539064B (en) | Neutral line and thickness extraction method and finite element grid automatic generation method | |
| CN103412985A (en) | Parametrization design method for trailing edge wedge slit of gas cooling blade | |
| CN105931297A (en) | Data processing method applied to three-dimensional geological surface model | |
| CN108804791B (en) | Aircraft parameterization method suitable for embedded air inlet layout | |
| CN107644139B (en) | Attribute mapping method from CAD model to CAE model |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C41 | Transfer of patent application or patent right or utility model | ||
| CB03 | Change of inventor or designer information |
Inventor after: Li Zhanwu Inventor after: Xu An Inventor after: Ou Jianjun Inventor after: Lv Yue Inventor before: Fang Zongde Inventor before: Bian Xiang Inventor before: Qin Kun Inventor before: Lian Lifei Inventor before: Gao Yang Inventor before: Hou Xiangying Inventor before: Zhang Baoyu |
|
| COR | Change of bibliographic data | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20160601 Address after: 710051 Baqiao Shaanxi District, Changle East Road, a word No. 1 Applicant after: Air Engineering Univ., PLA Address before: 710072 Xi'an friendship West Road, Shaanxi, No. 127 Applicant before: Northwestern Polytechnical University |
|
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20161130 Termination date: 20180831 |