CN109408893A - A kind of more thickness area dividing methods of inner plate of car door considering different operating conditions - Google Patents
A kind of more thickness area dividing methods of inner plate of car door considering different operating conditions Download PDFInfo
- Publication number
- CN109408893A CN109408893A CN201811122089.3A CN201811122089A CN109408893A CN 109408893 A CN109408893 A CN 109408893A CN 201811122089 A CN201811122089 A CN 201811122089A CN 109408893 A CN109408893 A CN 109408893A
- Authority
- CN
- China
- Prior art keywords
- thickness
- inner plate
- car door
- unit
- sensitivity
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Automobile Manufacture Line, Endless Track Vehicle, Trailer (AREA)
Abstract
The invention discloses a kind of more thickness area dividing methods of inner plate of car door for considering different operating conditions, comprising the following steps: by inner plate of car door it is uniformly discrete be multiple finite elements, and give corresponding working condition；It defines relevant parameter and all units of initial inner plate of car door is defined as maximum gauge unit；Finite element analysis is carried out to when front vehicle door inner plate, calculates the element sensitivity of all units；It is filtered and is averaged by element sensitivity of the mesh independent filter to all units；Calculate the target volume of following iteration step, and the threshold values of updating unit sensitivity；Reset the thickness design variable of each unit；Above step is repeated until reaching corresponding confined volume, and meets the condition of convergence: the more thickness of inner plate of car door obtained under different-stiffness operating condition are reconfigured and are matched.In the case where material thickness and given respective volume score, realizes that material distribution has more aggregation, quality can be greatly reduced under the premise of not weakening global stiffness.
Description
Technical field
The present invention relates to Mechanical Engineering Design technical fields, and in particular to how thick a kind of inner plate of car door for considering different operating conditions is
Spend region segmentation method.
Background technique
It is likely to be by not other than being formed by different materials grade welding for actual Tailor-welded Blank Structure
Same thickness component composition, such as tailor welded inner plate of car door.At this moment just it needs to be determined that different-thickness region line of demarcation, in other words
Say be exactly the more thickness areas of inner plate of car door reasonable layout again (thickness matching optimization) so that caused result is given
The gross mass that door device is reduced under load-up condition, obtains maximum global stiffness, and each interregional line of demarcation is naturally
It may be considered the position of weld seam, to instruct welding process (Fig. 2) before welding, so need to be to the more caliper zones of inner plate of car door
Regional partition is designed the research in method.Although having had the research much about Tailor-welded Blank Structure design aspect, phase
It can't be limited to a certain extent and promote and apply in the design for truly instructing vehicle body parts by closing achievement,
For example the location layout of weld line relies primarily on the experience of engineer, up to the present there are no a set of effective design sides
Method needs to select suitable design method to further increase the practicability of more thickness tailor welded Element Design results and reliable
Property.
It is exactly given in thickness for considering the more thickness area dividing methods of the inner plate of car door of different operating conditions intrinsically
In the case of, the region that every kind of thickness of redistribution is occupied, and then realize that inner plate of car door is structure lightened between performance design
Optimum balance, this is a kind of topological behavior in fact.Topological optimization technology is that have prospect in structure optimization technology, with novelty
Technology is to switch to structural topology problem mathematically to design optimal material distribution in certain space or region,
The basic ideas and purpose of design be exactly seek the rigidity of structure in design space optimal distribution form or structure it is best
Path of Force Transfer obtains the most light design of weight with certain performances of optimum results or under conditions of meeting various performances.So
And object applied by general Topology Optimization Method is the structure that there is dual extension-compression modulus to form, to different-thickness
The inner plate of car door structure that is composed simultaneously is not suitable for, so needing when requiring to be redistributed (piecemeal) to each thickness area
Upgrade primal algorithm and reconfigures the interpolating function for thickness and sensitivity analysis.Therefore the present invention selects two-way progressive open up
Optimization algorithm (Bi-directional Evolutionary Structural Optimization:BESO) is flutterred as with reference to blue
Algorithm is upgraded to carry out weight to various thickness the inner plate of car door structure being spliced by different-thickness by this
New layout.
Summary of the invention
The technical problem to be solved by the present invention is to, in view of the foregoing defects the prior art has, provide it is a kind of consideration not
With the more thickness area dividing methods of inner plate of car door of operating condition, in the case where material thickness and given respective volume score, realize
Material distribution has more aggregation, i.e. the line of demarcation in different-thickness region is more obvious, preferably solves poor thick welding type vehicle
The design problem of door inner plate structure weld seam layout, optimizes automobile door structure under the premise of not weakening global stiffness, can substantially drop
Low quality provides certain application foundation for further effective popularization of vehicle body tailor welded components.
Used technical solution is the present invention to solve above-mentioned technical problem:
A kind of more thickness area dividing methods of inner plate of car door considering different operating conditions, comprising the following steps:
1) by the design section of inner plate of car door it is uniformly discrete be multiple finite elements, and according to the rigidity property of inner plate of car door
Give corresponding working condition；
2) relevant parameter of Tailor-welded Blank Structure is defined, relevant parameter includes the plate for constituting the multi-thickness plate of Tailor-welded Blank Structure
Thick and respective thickness plate confined volume, and all units of initial inner plate of car door are defined as maximum gauge unit；
3) finite element analysis is carried out to when front vehicle door inner plate under corresponding rigidity operating condition, calculates all units and is based on thickness
The element sensitivity of attribute；
4) it is filtered and is averaged by element sensitivity of the mesh independent filter to all units；
5) target volume of the following iteration step of inner plate of car door is calculated, and according to new target volume updating unit sensitivity
Threshold values, following iteration step target volume are as follows:
In formula,For the current volume of jth kind thickness,For the confined volume of jth kind thickness, k is iteration step, and ER is
Progressive rate, progressive rate are that the quantity of increasing is deleted for control unit, obtain desired optimum structure and the Optimal Parameters that introduce；
6) according to the element sensitivity of all units and new element sensitivity threshold values comparison result, each unit is reset
Thickness design variable；
7) step 3) -6 is repeated) until thick plates various in car door inside reach corresponding confined volume, and under meeting
The condition of convergence in face:
In formula, N is positive integer, and η is convergence limit, Ck-i+1Tailor-welded Blank Structure is average soft when expression kth-i+1 iteration step
Degree, Ck-N-i+1Indicate that the average flexibility of Tailor-welded Blank Structure when kth-N-i+1 iteration step, N and the number of iterations do not have specific relationship,
The number of iterations is more than or equal to 2N, by 2N iteration until the average flexibility of structure reaches a stationary value, such as: N=5 is
Using the relative variation of nearest 10 iteration results as judgment criteria；Can also N=6, then be by nearest 12 iteration results
Relative variation is as judgment criteria；
8) it is distributed according to the finite elements of different-thickness, obtains the more thickness of the inner plate of car door group again under different-stiffness operating condition
It closes and matches, to obtain optimal different-thickness distributed areas.
According to above-mentioned technical proposal, in the step 1), corresponding work is given according to the rigidity property of inner plate of car door
Condition condition includes a variety of different working conditions, and specifically respectively sagging rigidity operating condition, upper torsion stiffness operating condition and lower torsion are rigid
Operating condition is spent, corresponding working condition is given and refers to any one given in a variety of different working conditions.
According to above-mentioned technical proposal, in the step 1), inner plate of car door is by n given thickness t1,t2,…,tnIt is single
Member composition, and t1>t2…>tn。
Further comprising the steps of in the step 2) according to above-mentioned technical proposal: relevant parameter further includes progressive rate
ER and thickness interpolation penalize index p；The confined volume of the various thick plates is the various thick plates in entire inner plate of car door
Respective volume score, maximum gauge unit refers to the thickness of most thick thick plates in relevant parameter that element thickness is definition.
According to above-mentioned technical proposal, optimal progressive rate ER=2%.
According to above-mentioned technical proposal, in the step 3), the specific calculating process of the element sensitivity includes
Following steps:
A) the average flexibility C of Tailor-welded Blank Structure is calculated:
In formula, K and u are the global stiffness matrix and motion vector of structure respectively, and f=Ku, f are load vectors；
B) by average flexibility C to design variable xijIt carries out differential and obtains element sensitivity:
In formula, design variable xijIt is engineer's variable of i-th of unit between+1 thickness of jth kind thickness and jth
Or density value,WithIt is by thickness T respectivelyjAnd Tj+1The element stiffness matrix for i-th of unit being calculated, p are thickness
Degree interpolation penalizes index, similar to the punishment index of material interpolation；When the value of p is smaller, occur in structural topology result figure
The thick middle angle value of major part, structure distribution are unintelligible；When p value is larger, the number of iterations increases, and is not easy to restrain, and p is most
Good value is 3, uiIt is the motion vector of i-th of unit.
According to above-mentioned technical proposal, in the step 4), element sensitivity is filtered and is put down by filter
Equal detailed process the following steps are included:
A) element sensitivity around each node is averaged, and then obtains node sensitivity:
In formula, M indicates the unit connecting with j-th of node sum, ωiIt is weighting or power of i-th of unit for node j
Weight coefficient, and Indicate j-th of node sensitivity；
B) it is found out by filter with the related each node of i-th of element sensitivity in the design section of structure, it will
Node sensitivity is mapped on unit, has element sensitivity of overall importance, filtered element sensitivity is calculated:
In formula, K indicates the sum of domain of influence interior nodes, ω (rij) it is weight coefficient, rijIndicate node from unit center
Distance；
C) element sensitivity before filtered element sensitivity and filtering is averaged:
In formula,For filtered element sensitivity,For the element sensitivity before filtering, k is iteration step.
According to above-mentioned technical proposal, in the step 5), once the obligatory point product value of each thick plates is determined, volume will
Constant is kept, then, the threshold values of updating unit sensitivity is determined by bisection method.
It is thick when the design section of inner plate of car door has 100 units in the step 5) according to above-mentioned technical proposal
Spend design variable xijCorresponding element sensitivity is ordered as α1,j>α2,j…>α100, j, as thickness tjVolume be in next step
When the y% of total volume, corresponding units sensitivity threshold value setting is αy,j, j is jth kind thickness, j=1,2 ..., n-1, design change
Measure xijIt is the engineer's variable or density value of i-th of unit between+1 thickness of jth kind thickness and jth.
According to above-mentioned technical proposal, in the step 5), the current volume of jth kind thicknessAre as follows:
In formula, NE indicates the unit sum in the design space region of inner plate of car door, ViFor the volume of i-th of unit, Vi *
For the confined volume of i-th of unit, the thickness design variable x of unitijIt is i-th between+1 thickness of jth kind thickness and jth
The engineer's variable or density value of unit, design variable xijIt is set to binary one or xmin, xminFor indicate close to
0 value.
According to above-mentioned technical proposal, in the step 6), the thickness design variable of unit is xij, it is thick in jth kind
The engineer's variable or density value of i-th of unit between degree and+1 thickness of jth, for xij=1 unit, if accordingly
Element sensitivity is smaller than threshold values, xijIt is turned into xmin；For xij=xminUnit, if element sensitivity is more than or equal to
When sensitivity threshold values, xijBecome 1；Meanwhile updating the new design variable of thickness profile.
The invention has the following advantages:
The method of the present invention is by the calculation formula based on thickness function and the sensitivity of thickness, in material thickness and corresponding body
In the case that fraction is given, realize that material distribution has more aggregation, i.e. the line of demarcation in different-thickness region is more obvious, compared with
Solves the design problem of poor thick welding type inner plate of car door construction joint layout well, practical value with higher and stronger
Engineering significance optimizes automobile door structure under the premise of not weakening global stiffness, and quality can be greatly reduced, and is vehicle body tailor welded zero
Further effective popularization of component provides certain application foundation.
Detailed description of the invention
Fig. 1 is the flow chart that the more thickness area dividing methods of the inner plate of car door of different operating conditions are considered in the embodiment of the present invention；
Fig. 2 is the schematic diagram of the continuous laser welding process of more thickness Tailor-welded Blank Structures in the embodiment of the present invention；
Fig. 3 is schematic diagram of the inner plate of car door under sagging rigidity operating condition in the embodiment of the present invention；
Fig. 4 is schematic diagram of the inner plate of car door under upper torsion stiffness operating condition in the embodiment of the present invention；
Fig. 5 is schematic diagram of the inner plate of car door under lower torsion stiffness operating condition in the embodiment of the present invention；
Fig. 6 is that the region of search formed in the embodiment of the present invention with the related each node of i-th of element sensitivity shows
It is intended to；
Fig. 7 is in the embodiment of the present invention in the progressive of the average flexibility of sagging rigidity operating condition exit door inner panel and volume fraction
History；
Fig. 8 is in the embodiment of the present invention in the average soft of upper torsion stiffness operating condition or lower torsion stiffness operating condition exit door inner panel
The progressive history of degree and volume fraction；
Fig. 9 is the inner plate of car door Optimum Design Results in the embodiment of the present invention under sagging rigidity operating condition；
Figure 10 is the inner plate of car door optimization in the embodiment of the present invention under upper torsion stiffness operating condition or lower torsion stiffness operating condition
Design result；
Specific embodiment
The present invention is described in detail with reference to the accompanying drawings and examples.
Shown in referring to Fig.1~Figure 10, the inner plate of car door of the different operating conditions of consideration in one embodiment provided by the invention is more
Thickness area dividing method, comprising the following steps:
1) by the design section of inner plate of car door it is uniformly discrete be multiple finite elements, and according to the rigidity property of inner plate of car door
Give different working conditions；
2) relevant parameter of Tailor-welded Blank Structure is defined, relevant parameter includes the plate for constituting the multi-thickness plate of Tailor-welded Blank Structure
Thick and respective thickness plate confined volume, and all units of initial inner plate of car door are defined as maximum gauge unit；
3) finite element analysis is carried out to when front vehicle door inner plate, calculates element sensitivity of all units based on thickness profile；
4) it is filtered and is averaged by element sensitivity of the mesh independent filter to all units；
5) target volume of the following iteration step of inner plate of car door is calculated, and according to new target volume updating unit sensitivity
Threshold values, following iteration step target volume are as follows:
In formula,For the current volume of jth kind thickness,For the confined volume of jth kind thickness, k is iteration step, and ER is
Progressive rate, progressive rate are that the quantity of increasing is deleted for control unit, obtain desired optimum structure and the Optimal Parameters that introduce；
6) according to the element sensitivity of all units and new element sensitivity threshold values comparison result, each unit is reset
Thickness design variable；
7) step 3) -6 is repeated) until thick plates each in Tailor-welded Blank Structure reach corresponding confined volumeAnd satisfaction
The following condition of convergence:
In formula, N is positive integer, and η is convergence limit, Ck-i+1Tailor-welded Blank Structure is average soft when expression kth-i+1 iteration step
Degree, Ck-N-i+1Indicate that the average flexibility of Tailor-welded Blank Structure when kth-N-i+1 iteration step, the number of iterations are more than or equal to 2N, pass through 2N
Secondary iteration is until the average flexibility of structure reaches a stationary value；
8) it is distributed according to the finite elements of different-thickness, obtains the more thickness of the inner plate of car door group again under different-stiffness operating condition
It closes and matches, to obtain optimal different-thickness distributed areas.
In one embodiment of the present of invention, selects the left front door device of body side wall structure as research object, consider
To the primary structure of car door, i.e. door inner panel plays main effect in keeping the static characteristics such as the rigidity of door device,
, using door inner panel as the research object of topological optimization, the purpose is to design to have minimum average B configuration flexibility under all load conditions for this
The door device of (or maximum rigidity).Under every kind of load condition, evaluated by applying displacement at load or amount of deflection
The rigidity of structure.
It is considered herein that two kinds of typical cases, it may be assumed that
(a) sagging rigidity: t1=1.3mm (Vf1=40%), t2=0.8mm (Vf2=30%) and t3=0.5mm (Vf3=
30%)；
(b) upper and lower torsion: t1=1.4mm (Vf1=50%), t2=1.0mm (Vf2=20%) and t3=0.6mm (Vf3=
30%).
The load that sagging rigidity is applied is concentrated mainly in plate, so it is considered plane stress problem；And upper torsion
Load with lower torsion stiffness is considered buckling problem perpendicular to plate.
Specific implementation step is as follows:
(1) design section of finite elements discrete uniform thickness inner plate of car door structure, and given boundary and load-up condition are used.
Numerical modeling is carried out by business finite element software ABAQUS, which defines shell unit, cell attribute, material properties, load
Lotus and boundary condition etc., it is considered herein that load working condition condition be sagging rigidity, it is upper torsion and lower torsion stiffness, finite element mould
Type is as shown in Fig. 3~Fig. 5.Load working condition corresponding to door device different-stiffness performance objective is different, such as under
It hangs down rigidity, (point P at hinge1And P2) six-freedom degree it is restrained, and (point P at door handle3) horizontal displacement, that is, direction x and y
On displacement it is restrained and the direction z apply F1The load of=800N；It is complete at hinge for upper torsion and lower torsion stiffness
Fixed, the freedom degree at door handle in addition to the rotation is restrained, two kinds of load (F2=F3=200N) it is applied along the direction y
In upper and lower position (the point P of door inner panel4And P5)。
(2) target volume and required parameter (volume constraint V are defined*With progressive rate ER=2%), will be used to determine next
The volume that every kind of thickness area is occupied in a iteration step.Inner plate of car door structure design under more rigidity conditions can regard in fact
For a multi-objective optimization design of power problem, i.e. minimum under different-stiffness is averaged flexibility.For this problem, weight can be tried out
Because of the influence degree of each objective function of sub-definite, that is, solve the problems, such as that this simple mode is exactly to pass through weighting flexibility to incite somebody to action
Multiple target is converted into the single-objective problem with given weight coefficient, in the more thickness layout optimization problems of inner plate of car door structure just
That corresponding target response is regarded as to the sum of corresponding flexibility weighting of every sub- rigidity operating condition, be consider more rigidity (or
Load) operating condition mentality of designing, the global response of actually one definition over the entire structure.Inner plate of car door structure is more
More thickness segmentation design problem under rigidity condition can give are as follows:
In formula, M is all rigidity situation sums；CkIt is the average flexibility under kth kind rigidity；wkIt is the power system of kth kind rigidity
Number, value depend on the relative importance of every kind of rigidity, andFor situation (b), in order to avoid optimizing inner plate of car door
There is biggish amount of deflection, i.e. rigidity is smaller, and upper torsion situation has more large deflection than lower torsion and has important engineering significance,
The weight factor of upper and lower torsion situation is simply taken as 0.67 and 0.33 respectively；K and u is global stiffness matrix and the position of structure respectively
The amount of shifting to；N indicates the unit sum in the region of design space, ViFor the volume of i-th of unit,For the volume of jth kind thickness,
Vi *For the confined volume of i-th of unit.
It is by unit " cell density " each in finite element model design space, i.e. x in topology optimization problemijAs
Design variable changes cell density during Optimized Iterative and couples it among stiffness matrix.It is assumed that in car door
Hardened structure is by n given thickness t1,t2,…,tn(t1>t2…>tn) composition.Binary design variable xijBe in jth kind thickness and
The engineer's variable or density value of i-th of unit between+1 thickness of jth, are generally set to binary one or 0, and
During Numerical implementation, number 0 uses xmin(such as 0.001) is unstable to avoid numerical value to replace.That is:
The principle of each units thick angle value selection of inner plate of car door structure is: if the unit of the object element after Optimization Solution
Density value is 1, then it represents that the thickness profile at the cell position is important；On the contrary, if cell density be or relatively close to
xmin, then illustrate the thickness at the cell position not and be especially important, need to assign lesser thickness value.The former unit needs
Retain (larger thickness value), the units thick angle value of the latter is smaller, can thus achieve the purpose that material efficiently utilizes, and then realize
The structure lightened optimum balance between performance design of inner plate of car door.
(3) finite element analysis and computing unit sensitivity are carried out to current design.When inner plate of car door structural material parameter ratio
When changing such as thickness profile, structural object response (such as averagely flexibility etc.) can also change therewith.Sensitivity is used
Analyze the relationship between expression structure Parameters variation and response variation, sensitivity coefficient can be expressed as inner plate of car door structure
Respond the gradient information to relevant parameter.In finite element analysis, the equilibrium equation of static structure can be provided:
Ku=f (3)
The single rigidity working condition exit door inner plate structure flexibility value that is averaged may be expressed as:
When inner plate of car door structure specified load, the sensitivity of average flexibility C can be by it to design variable xijIt carries out
Differential obtains, it may be assumed that
In order to determine the sensitivity of motion vector, the flexibility C that needs to be averaged, which passes through, introduces Lagrangian multiplier vector λ,
The flexibility C that is averaged in this way is corrected for:
The sensitivity of revised average flexibility are as follows:
Due to the presence of equilibrium equation, the Section 3 in above formula is just become zero.Assuming that the variation of a cell density is to applying
The load vectors added do not influence, soTherefore the sensitivity of average flexibility C can be corrected further are as follows:
By formula (5) it is found that because (f-Ku) is always zero, it is possible to it is any to select Lagrangian multiplier vector λ, be
Eliminate the unknown quantity in above formulaVector λ can choose are as follows:
Equilibrium equation is recombined, above formula can be expressed further are as follows:
Formula (10) are substituted into formula (8), the differential for the flexibility C that can must be averaged:
For problem designed by the present invention, that is, what is studied is that the more thickness area segmentation design methods of inner plate of car door are asked
Topic cannot provide unified interpolating function, but need because different functional relations are presented from thickness in element stiffness matrix
Carry out certain transformation.According to limited meta-knoeledge it is found that for any mechanics problem, element stiffness matrix can unify table
It reaches are as follows:
In formula:It is material matrix and strain displacement matrix respectively with B.It should be worth noting
, strain displacement matrix B is different plane stress problem and buckling problem.For the invention example, door device
Sagging rigidity be plane stress problem, it is upper torsion and lower torsion stiffness be buckling problem, therefore the function T (t) of thickness t respectively table
It reaches are as follows:
Based on the description above, the thickness of i-th of unit is designed variable xijControl.In order to which clear expression and progress are sensitive
Property analysis, interpolation directly is carried out to the function of interior thickness:
Element stiffness matrix can provide in this way:
In formula,WithIt is by thickness T respectivelyjAnd Tj+1The element stiffness matrix for i-th of unit being calculated.
If element sensitivity can be obtained by interpolating function formula (14) by carrying out optimal design to n kind thickness:
In this way, being directed to different mechanics problems, sensitivity can unify explicit expression are as follows:
Under the more rigidity operating conditions of convolution (1), sensitivity can be described as the formula with weight factor, it may be assumed that
It similarly, can be with explicit expression for the sensitivity of i-th of unit of jth kind thickness are as follows:
(4) element sensitivity is filtered and is averaged by mesh independent filter.By the list around each node
First sensitivity is averaged, and then obtains node sensitivity:
In formula, M indicates the unit connecting with j-th of node sum, ωiIt is weighting or power of i-th of unit for node j
Weight coefficient, and Indicate j-th of node sensitivity, wherein ωiDefinition such as formula (21):
In formula, rijFor the linear distance between j-th of node and i-th of unit center, physical significance be can be understood as
The center of unit illustrates that its influence to node sensitivity will be bigger closer to node.
Conventional topologies optimization method cannot be made up again after the rejecting of some unit, and two-way progressive Topology Optimization Method
Just avoid the problem, and new unit can be increased around structure thicker region, preferably solve such as mesh dependence,
Chessboard effect and being difficult to such as restrains at the problem in tradition optimization, and largely improves its computational efficiency.In order to make list
First sensitivity have it is of overall importance, node sensitivity must be mapped on unit, the given of filter is to find out design section
It is interior that the domain of influence center is exactly the center (Fig. 6) of i-th of unit with i-th of element sensitivity related each node, and i-th
A element sensitivity is given are as follows:
In formula, K indicates the sum of domain of influence interior nodes, ω (rij) it is weight coefficient, its calculation formula is:
ω(rij)=rmin-rij (23)
In formula, radius r is filteredmin=100mm；rijIndicate node with a distance from unit center.
According to above description it can be concluded that the element sensitivity in the entire domain of influence can be because of the definition of cell filter
And it obtains smooth.If there is the unit of more high sensitivity in the domain of influence, the very low unit of certain initial sensitivities is sensitive
Degree will obtain certain raising.The unit that is removed of early stage, by the filter they can by adding back again, thus
Achieve the purpose that add new unit around material " high-effect " region.In addition, during Optimized Iterative, due to averaging unit
The application (formula (24)) of sensitivity strategy, so preferably reducing the fluctuation between volume fraction and objective function, mentions significantly
The high convergent speed of optimization and robustness.
(5) target volume for defining following iteration step, for example, the current volume of jth kind thickness isHere k is iteration
Step, is decided to be in the volume of next step:
Once binding occurrence is determined, volume will keep constant.Then, the threshold values for updating design variable is carried out by bisection method
It determines.For example, having 100 units, sensitivity x in design sectioni1It is ordered as α1,1>α2,1…>α100,1.As thickness t1Volume exist
When being in next step the 20% of total volume, this threshold values can be set as α20,1。
(6) design variable of each unit is reset.For xij=1 unit, if remolding sensitivity threshold values is small, xijBecome
For xmin；For xij=xminUnit, if sensitivity be more than or equal to threshold values, xijBecome 1, updates thickness profile
New design just passes through finite element analysis next time and obtains.
(7) step (3)-(6) are repeated until reaching target volumeAnd meet the condition of convergence of formula (26).
In formula, N is positive integer, and η is convergence limit.
By above-mentioned solution procedure it can be concluded that the progressive history of the topology of the average flexibility of inner plate of car door and volume fraction
Such as Fig. 7 and Fig. 8.For situation (a) and (b), the average flexibility value of optimization design be respectively 3394.622Nmm and
300.126Nmm.The progressive history of volume fraction shows this paper design method from thickness t1Initial designs start；Next, by
Gradually from thickness t1It is converted to thickness t2, until thickness t2Reach the volume fraction that need to meet；Later, thickness t2Volume fraction protect
Hold constant, thickness t3Volume fraction be gradually increased until until meeting its volume constraint, objective function and each thickness area
Volume fraction when keeping stable convergence, entire Optimized Iterative process will stop.
The inner plate of car door shape of optimal topological optimization is as shown in Figure 9 and Figure 10, since inner plate of car door structure is more complicated
The optimal location of geometry, different-thickness region is more complicated than simple structure, boundary interface, that is, weld seam position of different-thickness
It sets indicated by the solid line., it is apparent that being all the largest in the thickness of body hinge bearing, i.e., regardless of situation
The position should have maximum stiffness characteristics.As a result, can be by laser welding technology according to the topological optimization, it will be different
Thickness and suitable hardened structure is welded into entire part, then tailor welded inner plate of car door structure is struck out, finally by resistance spot welding
Etc. programs i.e. constitute tailor welded door device.
In order to further evaluate optimum results, table 1 lists pair between initial designs and the numerical result of optimal design
Than.From maximum gauge t1Start to carry out initial designs, later by changing thickness from t1To t2T is arrived again3, it is total to gradually change reduction
Quality.(20% or more) is greatly reduced due to structure gross mass, obtained average flexibility only more by a small margin increase (less than
9%), especially situation a, only the 1% of loss structure rigidity, can be reduced close to 30% weight, it means that by mentioning
Method of topological optimization design out can be such that the gross mass of automobile door structure substantially drops under the premise of not sacrificing structure global stiffness
It is low, show that the design method is effective in door device light-weight design.
1 tailor welded inner plate of car door of table before optimization after average flexibility and quality versus
From the present embodiment it is found that design method according to the present invention can make the material point in inner plate of car door different-thickness region
Cloth has more aggregation, and chessboard effect and intermediate density region is not present.Different-thickness region is redistributed, and can protected
On the basis of holding rigidity property, reduce the gross mass of door device, and each interregional line of demarcation is more obvious, i.e. weld line
Position is high-visible, is capable of the layout problem of the thick tailor welded automobile door structure weld line of guide difference, to demonstrate
This method is in the feasibility for handling the different more thickness area segmentations of operating condition exit door inner panel.The present disclosure applies equally to other by not
The consitutional thickness matching optimization problem of stack pile component.There is preferable reference and reference price for Practical Project
Value.
In conclusion the invention proposes a kind of more thickness area dividing methods of inner plate of car door for considering different operating conditions.It should
Method, as design source, traditional interpolating function only for elasticity modulus is converted using two-way progressive topological optimization algorithm
For the interpolation model based on thickness function, it is deduced the Calculation of Sensitivity formula of thickness interpolating function, is had in design method
Certain novelty；The shortcomings that only " digging a hole " the present invention overcomes conventional topologies optimization method, according to actual engineering demand,
In the case where material thickness and given respective volume score, material distribution has more aggregation, i.e. point in different-thickness region
Boundary line is more obvious, preferably solves the design problem of poor thick welding type inner plate of car door construction joint layout, with higher
Practical value and stronger engineering significance.
Above is only presently preferred embodiments of the present invention, and of course, the scope of rights of the present invention cannot be limited by this,
Therefore according to equivalence changes made by scope of the present invention patent, still belong to protection scope of the present invention.
Claims (10)
1. a kind of more thickness area dividing methods of inner plate of car door for considering different operating conditions, which comprises the following steps:
1) by the design section of inner plate of car door it is uniformly discrete be multiple finite elements, and it is given according to the rigidity property of inner plate of car door
Corresponding working condition；
2) relevant parameter is defined, relevant parameter includes the final goal volume of various thick plates, and by the institute of initial inner plate of car door
There is unit to be defined as maximum gauge unit；
3) finite element analysis is carried out to when front vehicle door inner plate under corresponding rigidity operating condition, calculates all units and is based on thickness profile
Element sensitivity；
4) it is filtered and is averaged by element sensitivity of the mesh independent filter to all units；
5) target volume of the following iteration step of inner plate of car door is calculated, and according to the valve of new target volume updating unit sensitivity
Value, the target volume of following iteration step are as follows:
In formula,For the current volume of jth kind thickness,For the confined volume of jth kind thickness, k is iteration step, and ER is progressive
Rate；
6) according to the element sensitivity of all units and new element sensitivity threshold values comparison result, the thickness of each unit is reset
Design variable；
7) step 3) -6 is repeated) reach corresponding confined volume until thick plates various in inner plate of car door, and meet following
The condition of convergence:
In formula, N is positive integer, and η is convergence limit, Ck-i+1Indicate the average flexibility of Tailor-welded Blank Structure when kth-i+1 iteration step,
Ck-N-i+1Indicate that the average flexibility of Tailor-welded Blank Structure when kth-N-i+1 iteration step, the number of iterations are more than or equal to 2N, repeatedly by 2N times
In generation, is until the average flexibility of structure reaches a stationary value；
8) according to the finite elements of different-thickness be distributed, obtain different-stiffness operating condition under the more thickness of inner plate of car door reconfigure and
Matching, to obtain optimal different-thickness distributed areas.
2. the inner plate of car door more thickness area dividing methods according to claim 1 for considering different operating conditions, which is characterized in that
In the step 1), giving corresponding working condition according to the rigidity property of inner plate of car door includes a variety of different operating condition items
Part, specifically respectively sagging rigidity operating condition, upper torsion stiffness operating condition and lower torsion stiffness operating condition, give corresponding working condition and refer to
Give any one in a variety of different working conditions.
3. the inner plate of car door more thickness area dividing methods according to claim 1 for considering different operating conditions, which is characterized in that
Further comprising the steps of in the step 2): relevant parameter further includes that progressive rate ER and thickness interpolation penalize index p.
4. the inner plate of car door more thickness area dividing methods according to claim 1 for considering different operating conditions, which is characterized in that
In the step 3), the specific calculating process of the element sensitivity the following steps are included:
A) the average flexibility C of Tailor-welded Blank Structure is calculated:
In formula, K and u are the global stiffness matrix and motion vector of structure respectively, and f=Ku, f are load vectors；
B) by average flexibility C to design variable xijIt carries out differential and obtains element sensitivity:
In formula, design variable xijIt is engineer's variable or close of i-th of unit between+1 thickness of jth kind thickness and jth
Angle value,WithIt is by thickness T respectivelyjAnd Tj+1The element stiffness matrix for i-th of unit being calculated, p insert for thickness
Value penalizes index, and ui is the motion vector of i-th of unit.
5. the inner plate of car door more thickness area dividing methods according to claim 1 for considering different operating conditions, which is characterized in that
In the step 4), mesh independent filter refers to using i-th of unit center as dot, with rminFor the circle domain of radius,
Circle contains at least one unit in domain.
6. the inner plate of car door more thickness area dividing methods according to claim 1 for considering different operating conditions, which is characterized in that
In the step 4), element sensitivity is filtered by filter and average detailed process the following steps are included:
A) element sensitivity around each node is averaged, and then obtains node sensitivity:
In formula, M indicates the unit connecting with j-th of node sum, ωiIt is weighting or weight system of i-th of unit for node j
Number, and Indicate j-th of node sensitivity；
B) by filter find out in the design section of structure with the related each node of i-th of element sensitivity, by node
Sensitivity is mapped on unit, has element sensitivity of overall importance, filtered element sensitivity is calculated:
In formula, K indicates the sum of domain of influence interior nodes, ω (rij) it is weight coefficient, rijIndicate node with a distance from unit center；
C) element sensitivity before filtered element sensitivity and filtering is averaged:
In formula,For filtered element sensitivity,For the element sensitivity before filtering, k is iteration step.
7. the inner plate of car door more thickness area dividing methods according to claim 1 for considering different operating conditions, which is characterized in that
In the step 5), once the obligatory point product value of each thick plates is determined, volume will keep constant, then, updating unit spirit
The threshold values of sensitivity is determined by bisection method.
8. the inner plate of car door more thickness area dividing methods according to claim 7 for considering different operating conditions, which is characterized in that
In the step 5), when the design section of inner plate of car door has 100 units, thickness design variable xijCorresponding unit spirit
Sensitivity is ordered as α1,j>α2,j…>α100,j, as thickness tjVolume next step be total volume y% when, corresponding units sensitivity
Threshold value setting is αy,j, j is jth kind thickness, j=1,2 ..., n-1, design variable xijIt is in+1 thickness of jth kind thickness and jth
Between i-th of unit engineer's variable or density value.
9. the inner plate of car door more thickness area dividing methods according to claim 1 for considering different operating conditions, which is characterized in that
In the step 5), the current volume of jth kind thicknessAre as follows:
In formula, NE indicates the unit sum in the design space region of inner plate of car door, ViFor the volume of i-th of unit, Vi *It is i-th
The confined volume of a unit, the thickness design variable x of unitijIt is i-th of unit between+1 thickness of jth kind thickness and jth
Engineer's variable or density value, design variable xijIt is set to binary one or xmin, xminFor indicating close to 0
Value.
10. the inner plate of car door more thickness area dividing methods according to claim 9 for considering different operating conditions, feature exist
In in the step 6), for the thickness design variable x of unitij=1 unit, if corresponding element sensitivity ratio
Threshold values is small, xijIt is turned into xmin；For xij=xminUnit, if element sensitivity be more than or equal to sensitivity threshold values,
xijBecome 1；Meanwhile updating the new design variable of thickness profile.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811122089.3A CN109408893A (en) | 2018-09-26 | 2018-09-26 | A kind of more thickness area dividing methods of inner plate of car door considering different operating conditions |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811122089.3A CN109408893A (en) | 2018-09-26 | 2018-09-26 | A kind of more thickness area dividing methods of inner plate of car door considering different operating conditions |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109408893A true CN109408893A (en) | 2019-03-01 |
Family
ID=65465308
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811122089.3A Pending CN109408893A (en) | 2018-09-26 | 2018-09-26 | A kind of more thickness area dividing methods of inner plate of car door considering different operating conditions |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109408893A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102799704A (en) * | 2012-05-21 | 2012-11-28 | 湖南大学 | Optimal design method for spliced and welded plate automobile door |
| CN106372347A (en) * | 2016-09-08 | 2017-02-01 | 厦门大学嘉庚学院 | Dynamic response topological optimization method implemented by application of improved bi-directional evolutionary structural optimization (BESO) to equivalent static load method |
| KR101703450B1 (en) * | 2015-11-04 | 2017-02-09 | 한양대학교 산학협력단 | Structural topology optimization method using bigbang-bigcrunch algorithm |
-
2018
- 2018-09-26 CN CN201811122089.3A patent/CN109408893A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102799704A (en) * | 2012-05-21 | 2012-11-28 | 湖南大学 | Optimal design method for spliced and welded plate automobile door |
| KR101703450B1 (en) * | 2015-11-04 | 2017-02-09 | 한양대학교 산학협력단 | Structural topology optimization method using bigbang-bigcrunch algorithm |
| CN106372347A (en) * | 2016-09-08 | 2017-02-01 | 厦门大学嘉庚学院 | Dynamic response topological optimization method implemented by application of improved bi-directional evolutionary structural optimization (BESO) to equivalent static load method |
Non-Patent Citations (3)
| Title |
|---|
| X. HUANG ET AL.: "Convergent and mesh-independent solutions for the bi-directional evolutionary structural optimization method", 《FINITE ELEMENTS IN ANALYSIS AND DESIGN》 * |
| 徐峰祥等: "基于拓扑优化的拼焊板车门焊缝位置布局", 《同济大学学报(自然科学版)》 * |
| 陈刚强等: "拼焊板车门焊缝布置与尺寸优化设计方法研究", 《中国机械工程》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kohler et al. | PSO+: A new particle swarm optimization algorithm for constrained problems | |
| Deb | Nonlinear goal programming using multi-objective genetic algorithms | |
| CN105512385B (en) | The method that the more material topological optimizations of structure are realized using sequence power function interpolation method | |
| Van Zyl et al. | Operational optimization of water distribution systems using a hybrid genetic algorithm | |
| Cai et al. | A multi-objective chaotic particle swarm optimization for environmental/economic dispatch | |
| CN107844676A (en) | A kind of Structural Topology Optimization Design method based on more performance constraints | |
| CN106295880A (en) | A kind of method and system of power system multi-objective reactive optimization | |
| Kumar et al. | Multi-objective particle swarm optimization: an introduction | |
| CN109063357A (en) | The automobile component pad layout optimization method of adaptive multi-step variable domain based on topological optimization | |
| CN108829914B (en) | Structure and process integrated design method of FRP structural member | |
| CN113065268A (en) | Rapid multi-objective space reflector optimization design method based on proxy model | |
| CN110287512B (en) | Hybrid cellular automaton method for solving optimal thickness distribution of variable-thickness thin-wall structure under manufacturing constraint | |
| CN113159588A (en) | Intelligent logistics vehicle scheduling algorithm based on Internet of things technology | |
| CN107066663A (en) | A kind of truss structure Multidisciplinary systems Topology Optimization Method based on fully stress constraint criterion | |
| CN115907254A (en) | Vehicle route planning method based on evolutionary multi-objective reinforcement learning | |
| Lee et al. | Applying interactive genetic algorithms to disassembly sequence planning | |
| Ghatee et al. | Ranking function-based solutions of fully fuzzified minimal cost flow problem | |
| CN109408893A (en) | A kind of more thickness area dividing methods of inner plate of car door considering different operating conditions | |
| Reddy et al. | Capacitor placement using fuzzy and particle swarm optimization method for maximum annual savings | |
| Bui et al. | An introduction to multi-objective optimization | |
| CN109271714A (en) | A kind of thick welding type construction joint location determining method of difference | |
| Kogiso et al. | Reliability-based multiobjective optimization using the satisficing trade-off method | |
| Muhamediyeva et al. | Fuzzy model of transport demand | |
| Quagliarella et al. | GAs for aerodynamic shape design II: multiobjective optimization and multi-criteria design | |
| Lagos et al. | Combining Tabu search and genetic algorithms to solve the capacitated multicommodity network flow problem |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190301 |
|
| RJ01 | Rejection of invention patent application after publication |