Subdomain hybrid cellular automata method for solving car body thickness optimization

A subdomain hybrid cellular automata method for solving car body thickness optimization, comprising the following steps:

• S1. building an initially designed crash finite element model for thickness optimization of a car body structure;
• S2. building a cellular automata model for the thickness optimization of the car body structure in subdomains, wherein a cell internal energy density distribution for car body framework parts approaches a step target internal energy density function, the car body framework parts are selected from the group consisting of A-pillar, B-pillar, sill, roof-rail, front door, rear door, rear side member, seat crossbeam, front side member rear section, seat rear crossbeam, rear floor front crossbeam, roof front crossbeam, roof middle crossbeam and roof rear crossbeam, and defining thickness variables and field variables;
• S3. executing an outer loop:
• obtaining a cell internal energy density and a constraint function value at a current design point for each of the car body framework parts through simulation analysis, and updating a target mass for each of the car body framework parts using a penalty function method according to an extent, wherein the current design point violates a constraint boundary to the extent;
• S4. executing an inner loop:
• S4.1. constructing a step target internal energy density function, and updating a target internal energy density for each of the car body framework parts;
• wherein a process of constructing the step target internal energy density function is:
• S4.1.1. according to subscripts i and j of a cell Ωi,j, defining a sequence number for the cell using id (i, j)={circumflex over (N)}Ωi−1*(i−1)+j, j∈[1, {circumflex over (N)}Ωi], {circumflex over (N)}Ω0, =0, wherein {circumflex over (N)}Ωi−1is a number of cells in an (i−1)th subdomain;
• S4.1.2. traversing all the cells, and calculating a difference between internal energy densities Sid(k) of all the cells and an average S(k) of the internal energy densities in a kth outer loop:
• ΔSid(k)=Sid(k)−S(k), wherein
• S
• _
• (
• k
• )
• =
• 1
• ∑
• i
• =
• 1
• l
• ⁢
• ⁢
• N
• ^
• Ω
• i
• ⁢
• ∑
• i
• =
• 1
• l
• ⁢
• ⁢
• ∑
• j
• =
• 1
• N
• ^
• Ω
• i
• ⁢
• ⁢
• S
• Ω
• i
• ,
• j
• (
• k
• )
• is the average of the internal energy densities of all the cells in the kth outer loop;
• S4.1.3. determining “step points” and “step intervals”:
• traversing all the cells, and when ΔSid(k)*ΔSid+1(k)<0 is established, defining a subscript id of ΔSid(k) as a “step point,” wherein m “step points” form m+1 “step intervals”;
• S4.1.4. updating the “step points” and the “step intervals”:
• if idi+1−idi+1<Hthreshod is established, when i=1, deleting a “step point” id1, and a “step interval” is updated from [id0,id1] to [id0,id2];
• when i>1, deleting a “step point” idi−1, and the “step interval” is updated from [idi−1,idi] to [idi−2,idi];
• if idi+1−idi+1<Hthreshold is not established, retaining the “step points” and “step intervals” of S4.1.3;
• S4.1.5 constructing the step target internal energy density function:
• S
• *
• (
• h
• ,
• k
• )
• =
• {
• S
• 1
• *
• (
• h
• ,
• k
• )
• ,
• 1
• ≤
• i
• ⁢
• ⁢
• d
• ≤
• i
• ⁢
• ⁢
• d
• 1
• S
• 2
• *
• (
• h
• ,
• k
• )
• ,
• i
• ⁢
• ⁢
• d
• 1
• <
• i
• ⁢
• ⁢
• d
• ≤
• i
• ⁢
• ⁢
• d
• 2
• …
• S
• i
• *
• (
• h
• ,
• k
• )
• ,
• i
• ⁢
• ⁢
• d
• i
• -
• 1
• <
• i
• ⁢
• ⁢
• d
• ≤
• i
• ⁢
• ⁢
• d
• i
• …
• S
• m
• ′
• *
• (
• h
• ,
• k
• )
• ,
• id
• m
• ′
• -
• 1
• <
• i
• ⁢
• ⁢
• d
• ≤
• i
• ⁢
• ⁢
• d
• m
• ′
• S
• m
• ′
• +
• 1
• *
• (
• h
• ,
• k
• )
• ,
• id
• m
• ′
• <
• i
• ⁢
• ⁢
• d
• ≤
• ∑
• i
• =
• 1
• l
• ⁢
• ⁢
• N
• ^
• Ω
• i
• wherein Si*(h,k) is a target internal energy density in the “step interval” [idi−1,idi] in the kth outer loop and an hth inner loop,
• S4.2. updating a cell thickness using cell thickness update rules based on a PID control strategy;
• S4.3. judging whether the inner loop converges:
• if a convergence condition is not satisfied, returning to S4.1;
• if the convergence condition is satisfied, exiting the inner loop and performing S5;
• S5. judging whether the outer loop converges:
• if global convergence conditions are not satisfied, returning to S3 and updating the cell thickness using the inner loop;
• if the global convergence conditions are satisfied, outputting a crashworthiness optimization design of a car body based upon internal energy densities of the car body framework parts.