# Numerical Investigation and Mold Optimization of the Automobile Coat Rack Compression Molding

Authors: Youmin Wang; Xiangli Li; He Sui
Journal: Advances in Materials Science and Engineering
 (2021)

---

## Abstract

In order to have more accurate control over the compression molding of automobile coat rack, improve the quality of molding products, and achieve the goal of lightweight design, a novel mechanical model for the main two-layer composite structure of the coat rack is proposed. In this regard, the main factors affecting the mechanical properties of the composite structure are obtained. The hot air convection is selected for the sheet preheating. During the experiment, the hot air temperature, preheating time, molding pressure, and pressing holding time are set to 250°C, 110 s, 13 MPa, and 80 s, respectively. Moreover, the error compensation method is applied to compensate for the shrinkage of the product during solidification and cooling. The LS-DYNA finite element software is used to simulate the molding process of the main body of the coat rack, and the node force information with large deformation is obtained accordingly. The load mapping is used as the boundary condition of mold topology optimization, and the compression molding of the main body of the coat rack is optimized. A lightweight design process and method for the compression molding of automotive interior parts and a mathematical model for the optimization of the solid isotropic material penalty (SIMP) (power law) material interpolation of the concave and convex molds are established. Based on the variable density method, OptiStruct is used for the lightweight design of the convex and concave molds of the main body of the coat rack, which reduces the mold weight by 15.6% and meets the requirements of production quality and lightweight.


---

## Body

## 1. Introduction

Car coat rack is defined as a vehicle accessory to provide beautiful coverage and protection for the covered area or system. In addition, it is an important acoustic component for blocking noise from the luggage compartment. Accordingly, the car coat rack has been widely used in almost all vehicles. With the increasing integration of automobile functions, the coat rack has become a carrier for many other parts. Recently, improving the molding quality of automotive interior parts has focused on controlling the mechanical properties of products, the optimal combination of process parameters, and mold optimization. More specifically, Zheng et al. [1] performed CAE analysis and physical tests on coat racks to find control parameters that affect the product strength and the load capacity. Moreover, they considered the requirements of product appearance and obtained optimal control parameter values to achieve the design goals, including low-cost and high-performance products. However, the mechanical properties of the material and the mechanical model of the coat rack were not considered. Guo et al. [2] proposed a method for the recycling felt waste of PP glass fiber composite and studied the impact of the size of crushed particles and the proportion of crushed material on the formability and mechanical properties of the coat rack. However, the thickness impact of the glass fiber sheet on the mechanical properties of the material was not considered so that the formulation of process parameters was not reliable. Miao et al. [3] discussed the effect of the multicomponent fiber composition process on the performance of the car coat rack. Through investigating the ratio, pressure, temperature, and speed of multicomponent fibers, the orthogonal test was carried out to find the best component. Studies showed that the combination scheme is an effective way to produce car coat racks and investigate all aspects of their performance. However, the experimental method was not used to calculate the material shrinkage so that the solidification shrinkage and the product cooling were not compensated. Yu [4] simulated, analyzed, and performed experiments on the nonwoven fabric hot-molding mold and the molding process and showed that the molding temperature is the most influential factor in the hot-molding process. Moreover, the temperature distribution was studied, which solved the problem of repeatedly adjusting the molding temperature based on experience, but the number of research objects and simulation experiments was small. Shao [5] optimized the process parameters of glass mat reinforced thermoplastics (GMT) composite material compression molding and obtained the optimal combination from the theoretical point of view. In this optimization, no analysis software was applied to optimize the mold. Demirci and Yildiz [6] studied the effect of conventional steel, new generation DP-TRIP steels, AA7108-AA7003 aluminum alloys, AM60-AZ31 magnesium alloys, and crash-box cross-sections on crash performance of thin-walled energy absorbers which are investigated numerically for the lightweight design of vehicle structures. However, it does not involve numerical research on nonmetallic materials such as car coat racks. It is worth noting that none of the foregoing studies has simulated the compression molding, while the molding simulation can accurately evaluate the correlation between level factors and the mold quality. In order to resolve this shortcoming, it was intended to investigate the mechanical properties of thermoplastic composite materials and design the car coat structure. In this regard, the LS-DYNA finite element software was applied to simulate the molding process, and the load mapping was used as the boundary condition for mold topology optimization. Finally, mold structure optimization was established and the mold structure was optimized to improve the molding quality of the coat rack.


## 2. Investigating the Mechanical Properties and Structural Design of Thermoplastic Composites for Automobile Coat Racks

In this section, it is intended to establish a mechanical model for the thermoplastic composite of the car coat rack (the mechanical model is solved), analyze some test cases, determine the thickness of the main body, complete the formulation of the main body molding process parameters, determine the shrinkage compensation value, and calculate the heat transfer in the cooling system. These items are discussed in detail in the following.


### 2.1. Establishment of Structural Mechanical Model for GMT/PET Composite

In the present paper, the composite two-layer structure is taken as the research object. The upper layer is a decorative layer, the main specimen is made of needle-punched PET fabric, and the lower layer is a GMT sheet. It should be indicated that the mechanical properties of GMT sheet depend on the spatial distribution of the fiber and PP and the combination between them [7]. In the composite structure, the materials, thickness, and other physical parameters of the upper and lower layers are different. Accordingly, the composite plate is regarded as a thin plate, and then its mechanical model is established by considering variations of different physical parameters along the Z-axis. Figure 1 shows the mechanical model of the composite structure.Figure 1
Mechanical model of the composite structure.In order to simplify the mechanical model and the corresponding solving method, it is assumed that the glass fiber and PP materials are evenly dispersed into the space occupied by the sheet, and the overall material performance is linearly combined by the volume fraction of each component. Under a specific glass fiber mass fraction, the GMT sheet has a specific material constant, which can be simplified to an isotropic material in the analysis.The volume fraction of PP and glass fiber in the GMT sheet can be expressed as follows:(1)Vp=ρ1⋅Mpρp,Vb=ρ1⋅Mbρb,where ρ1, ρp, and ρb are the total density of GMT sheets, the bulk density of polypropylene, and glass fiber at room temperature, respectively. Moreover, Mp and Mb denote the mass ratio of polypropylene and glass fiber, respectively. Then the material constant of the GMT sheet can be calculated. It is worth noting that the sum of quality scores satisfies the following expression:(2)Mb+Mp=1.Since the GMT sheet may contain air and cavities during the molding process, the sum of the practical volume fraction of glass fiber and PP is(3)Vb+Vp≤1.Moreover, the void volume rate is defined as(4)Vc=1−Vp+Vb=1−ρ1⋅Mbρp−ρ1⋅Mbρb.Let the physical property parameter of the GMT sheet beP, where P can represent the elastic modulus of the material E, Poisson’s ratio μ, linear expansion coefficient α, and so on. Then this parameter in the composite is defined as follows [8]:(5)P=Pp⋅Vp+Pb⋅Vb,where Pp and Pb are the physical properties of polypropylene and glass fiber, respectively. Similarly, the elastic modulus and Poisson’s ratio of the GMT sheet can be calculated from equations (6) and (7), respectively.(6)E1=Ep⋅Vp+Eb⋅Vb,(7)μ1=μp⋅Vp+μb⋅Vb,where Ep and Eb are the elastic modulus of polypropylene and glass fiber, respectively. Moreover, μp and μb are Poisson’s ratio of polypropylene and glass fiber, respectively. The upper layer is a single material of needle-punched PET fabric, so the elastic modulus and Poisson’s ratio are E2 and μ2. Then the elastic modulus of the two-layer composite structure is divided into two layers:(8)Ez=E1,−h1,h2,E2,h2,h3.The same operation can be conducted for Poisson’s ratio in the following form:(9)μz=μ1,−h1,h2,μ2,h2,h3.The differential equation of the elastic surface of the thin plate can be expressed in the following form:(10)D∇4w=q,where D is the bending stiffness of the sheet, D=Eh3/121−μ2, and h is the total thickness of the sheet. The left side of equation (10) D∇4w denotes the elastic resistance per unit area, and w is the deflection. Moreover, the left side q is the lateral load per unit area. The stress differential equation can be simplified as follows:(11)∂σZσZ=Ez2⋅1−μz2⋅h24−z2⋅∇4w.Integrating both sides with respect toz from −h1 to h3 yields the following expressions:(12)∫−h1h3∂σZσZdz=σZ|h3−h1=0−q=−q,∫−h1h3Ez21−μz2h24−z2∇4wdz=∫−h1h3Ez21−μz2h24−z2dz∇4w.In order to simplify expressions, parameterD is defined as follows:(13)D=−∫−h1h3Ez21−μz2h24−z2dz.Then the integral can be decomposed into the sum of integrals of each layer:(14)−D=∫−h1h2Ez21−μz2h24−z2dz+∫h2h3Ez21−μz2h24−z2dz=−D1−D2.


### 2.2. The Mechanical Model of the Composite Structure Is Solved

When solving the thin plate deflection-bending problem by the Navier method, the deflectionw should be initially obtained from the elastic surface differential equation under the boundary conditions of the plate edge [9].During the unsupported subsidence, the length and width for a rectangular thin plate, which is simply supported on four sides, area and b, respectively. In this case, the boundary conditions are(15)wx=0=0,∂2w∂x2x=0=0wx=a=0,∂2w∂x2x=a=0wy=0=0,∂2w∂x2y=0=0wy=b=0,∂2w∂x2y=b=0.The deflectionw can be calculated from the following trigonometric series:(16)w=∑m=1∞∑n=1∞wm=∑m=1∞∑n=1∞Amnsinmπxasinnπyb,where m and n are positive integers. Substitute this expression into equation (15) which indicates that all boundary conditions are satisfied. In order to find the coefficient Amn, equation (16) is introduced into equation (10):(17)π4D∑m=1∞∑n=1∞m2a2+n2b2Amnsinmπxasinnπyb=q,and expanding the load q on the right side of equation (17) into a multiple triangle series results in the following expression:(18)q=∑m=1∞∑n=1∞Amnsinmπxasinnπyb,where A can be solved through the method of determining the triangle series, and the solution is(19)Amn=4ab∫0a∫0bqx,ysinnπybsinmπxbdxdy.Then equation (18) can be rewritten in the following form:(20)q=4ab∑m=1∞∑m=1∞4ab∫0a∫0bqx,ysinnπybsinmπxbdxdy⋅sinnπybsinmπxb,(21)Amn=4∫0a∫0bqx,ysinnπy/bsinmπx/bdxdyπ4abDm2/a2+n2/b2.When the thin plate is subjected to a uniform load,q becomes a constant q0. Under this circumstance, equation (21) can be rewritten in the following form:(22)Amn=16q0π6mnDm2/a2+n2/b22,m=1,3,5,…;n=1,3,5,….Then equation (16) can be applied to calculate the deflection:(23)w=16q0π2D∑m=1,3,5,…∞∑m=1,3,5,…∞sinmπx/asinnπy/bmnm2/a2+n2/b22.Internal forces can be expressed in the following form:(24)Mx=−D∂2w∂x2+μ∂2w∂y2,My=−D∂2w∂y2+μ∂2w∂x2,Mxy=Myx=−D1−μ∂2w∂x∂y,FS,x=−D1−μ∂2∂x∇2w,FS,y=−D1−μ∂2∂y∇2w.Based on the performed calculations and analyses, the two-layer composite structure is simply supported on four sides, where the lengtha, width b, and GMT sheet thickness are 0.1 m, 0.1 m, and 3.5 mm, respectively. Moreover, the composite density is 1.3 g/cm3 and the void ratio is assumed to be 0. The mass fraction and density of the glass fiber are 30% and 2.6 g/cm3, respectively. In the present study, alkali-free fibers (E) with an elastic modulus of 72 GPa and Poisson’s ratio of 0.2 are applied. Furthermore, the density, elastic modulus, and Poisson’s ratio of PP are 0.93 g/cm3, 1.4 GPa, and 0.4, respectively. The lateral uniform load is 10 MPa, the thickness of the transition layer is 1.5 mm, and the specimen is made of needle-punched PET with an elastic modulus of 4000 MPa.Based on the established mathematical model, the thickness impact on the mechanical properties of the thin plate can be investigated by changing the thickness of the GMT sheet under the condition that the glass fiber content is maintained 30%, while the applied load and constraints remain constant. The iterative calculations are performed in the MATLAB environment. The maximum deflection and stress value of the two-layer structure are shown in Table1. The stress value can be obtained from the internal force value through the generalized Hook definition of the thin plate.Table 1
Maximum deflection and stress values for different thicknesses.
GMT plate thickness (mm)Deflection (m)σx, σy (Pa)τxy, τyx (Pa)3.30.002289934.91424.5983.40.001938030.42821.4783.50.001654726.35718.7593.60.001424223.45816.1023.70.001198419.85713.478Table1 indicates that, for a constant load and constraint, as the thickness of the GMT sheet increases, the deflection of the two-layer composite structure gradually decreases, thereby reducing the normal stress and the shear force. It is inferred that, as the sheet thickness increases, flexural, tensile, and shear resistance of the whole plate improve.When the thickness is fixed and different glass fiber ratios are used, the deflection and internal force of the plate are shown in Table2.Table 2
Deflection and internal force values for different glass fiber contents.
GMT plate glass fiber content (%)Deflection (mm)Bending moment (N⋅m)Torque (N⋅m)250.001673225.43429.562300.001654725.29329.646350.001623825.17729.739400.001596625.05429.822450.001562424.96429.018Table2 indicates that, under the condition of unchanged load and restraint, as the content of glass fiber increases, the deflection value and bending moment value decrease, while the corresponding torque increases. Therefore, the optimal content of GMT sheet glass fiber is selected according to the required mechanical performance index.


### 2.3. Thickness of the Car Coat Rack

Based on the foregoing discussions, it is concluded that the rigidity design requirements can be achieved by increasing the thickness of the GMT sheet without changing the main structure of the car coat rack. However, the material consumption for too thick products cannot be justified. More specifically, the production cost increases, the product becomes prone to stress concentration, and the molding quality reduces.The stiffness test principle is as follows: as a support frame, the sample deflects continuously at the midpoint of the span until the fracture or deformation reaches a certain predefined value. In this process, the applied stress to the sample is measured. Figure2 illustrates the schematic configuration of composite layers, indicating that the composite consists of a PET fabric and the GMT sheet.Figure 2
Schematic configurations of the composite layers.Since the fabric and the substrate are subjected to stress in parallel, the mixed law of elastic modulus of the parallel composite structure can be used [9]:(25)E0=E1V1+E2V2,where E0, E1, and E2 denote the elastic modulus of the composite material, PP glass fiberboard, and PET fabric, respectively. Substituting equations (6) into (25) gives(26)E0=EpVp+EbVbV1+E2V2.Then the predicted value of the elastic modulus in the composite shown in Figure2 can be expressed in the following form:(27)E0=EpVp+EbVbd1d1+d2+E2d2d1+d2,where d1 and d2 denote the thickness of PP glass fiber and PET fabric, respectively. In order to meet the rigidity design requirements, the following conditions should be satisfied:(28)E0I≥σ,I=b⋅d1+d2312,where σ is the required stiffness; I is the moment of inertia; and b is the width of the coat rack. In order to solve the equations easily, a thick needle-punched PET fabric with a thickness of 1.5 mm, elastic modulus of 4000 MPa is considered in the calculations, and GMT sheet porosity is assumed to be 0. The glass fiber is made of alkali-free fiber (E) with an elastic modulus of 72 GPa, while the elastic modulus of PP is 1.4 GPa. Moreover, the coat rack width b is 600 mm and the required bending stiffness is σ ≥ 60 N/mm2.The combination of glass fiber content and thickness of the GMT sheet that satisfy the stiffness requirements can be calculated from equations (27) and (28). The proposed combination is shown in Table 3.Table 3
Proposed GMT sheet glass fiber ratio and thickness combination.
GMT sheet glass fiber content (%)GMT sheet thickness (mm)Round thickness (mm)204.8865254.6364.8304.3864.5354.1364.3403.8864Based on the obtained results, a glass fiber content of 30% and a thickness of 4.5 mm is proposed as the selected composition. To maximize the safety factor, the thickness of the GMT sheet is rounded to 4.5 mm. Since the thicknesses of the GMT sheet and PET fabric are 4.5 mm and 1.5 mm, the overall thickness of the composite is 6 mm.


### 2.4. Formulation of the Process Parameters to Mold the Main Body of the Automobile Coat Rack

By analyzing the common quality defects and causes of the main body of the coat rack, the main processing parameters that affect the molding process of the car coat rack are determined. These parameters are the molding pressure, holding time, and the molding temperature.


#### 2.4.1. Sheet Warm-Up Time

In order to heat the GMT sheet, hot air convection preheating is used. In this case, the sheet can be expanded to 2-3 times the thickness of the cold sheet. It should be indicated that the sheet temperature should be less than 230°C; otherwise, the PP polypropylene will degrade, thereby affecting the surface of the molded product and reducing the heat resistance and strength of the product [10, 11]. The preheating time of the GMT sheet can be calculated through the following expression [12] (29):(29)t=70.9521+0.2318d−0.2810T,where t is the appropriate preheating time, min; d is the sheet thickness, mm; and T is the hot air temperature, °C.During the experiment, hot air flow from an oven is used to preheat the sheet. The thickness of the GMT sheet is 4.5 mm and the preheating temperature of the hot air is 250°C. Based on equation (29), the preheating time is adjusted to 99 s. In the experiment, the sheet is taken out of the oven and transferred to the mold. Since the temperatures of the material and the oven reduce during the material transfer, the preheating time of the sheet is set to 110 s.


#### 2.4.2. Molding Pressure and Mold Clamping Speed

Appropriate molding pressure not only improved the product quality but also made up for the shortage if other industrial parameters were not adequate [13].(1)
The method of maximum gauge pressure:(30)ftable=f1f2ST,whereftable, f1, and f2 denote the gauge pressure, unit pressure required by the product, and the rated gauge pressure of the press, respectively. Meanwhile, S and T denote the projected area of the product in the compression direction and the pressing force, respectively.(31)x=f1ST,wherex is the ratio of the force acting on the pressure direction of the product to the pressing force. When x > 1, the pressing force is too small to meet the production requirements. On the other hand, when x < 1, the pressing force is sufficient to meet the production demand.(2)
The method of Piston area method:(32)ftable=f1SF,whereF is the cross-sectional area of the piston of the press.In the present study, for a four-post-servo-hydraulic machine (model; company) withF = 0.5 m2 and S = 0.84 m2, the molding pressure is 13 MPa and ftable= 22 MPa. For these parameters, equation (31) indicates that x > 1, so that the selected hydraulic pressure meets the molding requirements.In the molding process to produce the main body of the automobile coat rack, the mold clamping speed is not a fixed value. It should be indicated that the molding process can be divided into 5 main stages:①
When the sheet is placed in the mold, the concave mold decreases at a speed of 60 mm/s.②
The interval between the convex and concave mold gradually decreases. Consequently, when the guide column enters the guide sleeve, the speed becomes 8 mm/s. Moreover, when the convex and concave molds are closed, the sheet is molded under pressure.③
Cooling and solidification.④
After holding the pressure for 80 s, the concave mold is split at a speed of 8 mm/s.⑥
After the split mold is completed, the concave mold is 60 mm/s. Then the speed rises and when the die is reset, the product is removed.These five stages are also shown in Figure3.Figure 3
Speed variation graph during the mold clamping.


#### 2.4.3. Holding Time

The cooling time calculation results are utilized as a reference to select the holding time. For large products such as car coat racks, the length and width are far beyond their thickness, and they have two-dimensional characteristics. However, the heat is only transferred in the one-dimensional direction, that is, the thickness direction of the part. Therefore, the cooling of such products can be conducted as follows.Heat transfer can be calculated from the following expression:(33)∂T∂θ′=k∂2T∂2x2,where T and k denote the average temperature when the sheet is demolded and the thermal diffusion coefficient of the sheet, respectively.In order to separate the variables, the following equation is used:(34)θ′=−2π2x12klnπT−TwT0−Tw,k=λρ⋅cp,where x1 is half the thickness of the part, m; Tw is the mold cavity temperature, °C; T0 is the sheet molding temperature, °C; ρ is the density of the sheet, kg/m3; λ is the thermal conductivity of the sheet, kW/m⋅°C; and cp is the specific heat capacity of the sheet, kJ/kg⋅°C.The process parameter values are shown in Table4.Table 4
The value of process parameters.
Parameter nameValueThe thickness of the main body of the car coat rack (mm)4.5The average temperature of the mold cavity wallTw (°C)40The molding temperatureT0 (°C)220The average temperature when the sheet is demoldedT (°C)80The density of the sheetρ (kg/m3)1130The specific heat capacity of the sheetcp (kJ/kg·°C)1.55The thermal conductivity of the sheetλ (kW/m·°C)0.33 × 10−3Substituting the abovementioned parameters into equation (34), the cooling time of the molding of the main body of the automobile coat rack is mathematically expressed as follows:(35)k=0.33×10−31130×1550×10−3=1.88×10−7,t=−2π20.002252klnπ80−404220−40≈74s.Combining with the previous pressure holding time, the pressure holding time is determined to be 80 s.


### 2.5. Shrinkage Compensation Value

The molding shrinkage rate of the product is measured experimentally. The experimental process is as follows: the main model of the coat rack of the reference model is selected. Moreover, 20 pieces of GMT plates with a thickness of 4.5 mm glass fiber content of 30% and a needle-punched PET fabric of 1.5 mm thickness are selected. It should be indicated that, in the experiment, the sheet preheating time, the pressure holding time, molding pressure, and mold closing speed utilize the abovementioned parameters. When the mold is opened, the distance between the selected positioning points is measured, and each piece is averaged at one measurement in the length and width directions. After the product is left at room temperature for 24 hours, the selected positioning point is measured again, using equation (29). Finally, the average shrinkage rate of the composite two-layer structure of GMT sheet + needle-punched PET fabric is 0.354%.


### 2.6. Heat Transfer Calculation for the Cooling System of the Automobile Coat Rack Compression Molding Mold

The cooling system adopts the water cooling method, and a cooling water channel is laid on the mold surface to pass the cooling water. In this section, the design and calculation of the cooling system of the forming mold of the main body of the car coat rack are established.


#### 2.6.1. The Total Heat of the Sheet into the Cavity

In a molding process, the total heat transferred from the sheet to the mold is mathematically expressed as follows:(36)Q1=Cp⋅T0−T⋅V⋅ρ,where T0 and T denote the sheet molding temperature and the average temperature when the sheet is demolded, respectively. Moreover, ρ, cp, and V are the density of the sheet, the specific heat capacity of the sheet, and the volume of the sheet, respectively. It is worth noting that the volume of the sheet is set to 1.4 m × 0.6 m × 0.0045 m.


#### 2.6.2. Heat Taken away by the Cooling Water

The heat taken away by the mold cooling water in one press molding is expressed as follows:(37)Q2=Q1−Qc−QR−QL,where Qc, QR, and QL denote the convection heat transfer from the mold to the air, the radiation heat transfer from the mold to the air, and the heat transfer from the mold to the hydraulic machine through the upper and lower mold bases, respectively.Under normal circumstances, 90%–95% of the heat brought into the sheet can be taken away by the cooling water through the mold cooling channel. Therefore, during the design process,Q2 = Q1.Since the main body of the car coat rack is a large thin-walled product, it can be designed according to 50% of the heat taken by the concave mold and the convex mold; then,(38)Q2a=Q2t=0.5Q2.


#### 2.6.3. Cooling Water Hole Diameter

The volume flow of the cooling water can be calculated according to equation (39) during the molding process of the die and the die of the die in one molding process:(39)V2a=Q2aρ1c1θ1−θ2⋅t,V2t=Q2tρ1c1θ1−θ2⋅t,where θ1 and θ2 are the inlet and outlet temperatures of the mold cooling water, respectively. It should be indicated that the temperature difference should not be too large and its maximum value should not exceed 5°C. Moreover, V (m3⋅min−1) and ρ1 (kg/m3) denote the volume flow of cooling water and the density of water, respectively. Furthermore, c1 (kJ/kg⋅°C) and t (min) are the specific heat capacity of water and the cooling time of one molding, respectively.The density of water is 1000 kg/m3, and the specific heat capacity of water is 4.2kJ/kg⋅°C, which is calculated according to equations (36) and (38):(40)V2a=V2t=0.5Q2ρ1c1θ1−θ2⋅t=0.5⋅Cp⋅T0−T⋅V⋅ρρ1c1θ1−θ2⋅t=16.275m3/min.According to the correlation between the cooling water pipe volume flow and the pipe diameter, the cooling water pipe diameter is determined asD = 27 mm.The correlation between the distanceP of the cooling water pipe, the distance H from the profile, and the diameter D of the channel is P = (3∼5) D, H = (1∼2) D; from D = 27 mm, take the pipe distance P = 80 mm, H = 30 mm.


## 2.1. Establishment of Structural Mechanical Model for GMT/PET Composite

In the present paper, the composite two-layer structure is taken as the research object. The upper layer is a decorative layer, the main specimen is made of needle-punched PET fabric, and the lower layer is a GMT sheet. It should be indicated that the mechanical properties of GMT sheet depend on the spatial distribution of the fiber and PP and the combination between them [7]. In the composite structure, the materials, thickness, and other physical parameters of the upper and lower layers are different. Accordingly, the composite plate is regarded as a thin plate, and then its mechanical model is established by considering variations of different physical parameters along the Z-axis. Figure 1 shows the mechanical model of the composite structure.Figure 1
Mechanical model of the composite structure.In order to simplify the mechanical model and the corresponding solving method, it is assumed that the glass fiber and PP materials are evenly dispersed into the space occupied by the sheet, and the overall material performance is linearly combined by the volume fraction of each component. Under a specific glass fiber mass fraction, the GMT sheet has a specific material constant, which can be simplified to an isotropic material in the analysis.The volume fraction of PP and glass fiber in the GMT sheet can be expressed as follows:(1)Vp=ρ1⋅Mpρp,Vb=ρ1⋅Mbρb,where ρ1, ρp, and ρb are the total density of GMT sheets, the bulk density of polypropylene, and glass fiber at room temperature, respectively. Moreover, Mp and Mb denote the mass ratio of polypropylene and glass fiber, respectively. Then the material constant of the GMT sheet can be calculated. It is worth noting that the sum of quality scores satisfies the following expression:(2)Mb+Mp=1.Since the GMT sheet may contain air and cavities during the molding process, the sum of the practical volume fraction of glass fiber and PP is(3)Vb+Vp≤1.Moreover, the void volume rate is defined as(4)Vc=1−Vp+Vb=1−ρ1⋅Mbρp−ρ1⋅Mbρb.Let the physical property parameter of the GMT sheet beP, where P can represent the elastic modulus of the material E, Poisson’s ratio μ, linear expansion coefficient α, and so on. Then this parameter in the composite is defined as follows [8]:(5)P=Pp⋅Vp+Pb⋅Vb,where Pp and Pb are the physical properties of polypropylene and glass fiber, respectively. Similarly, the elastic modulus and Poisson’s ratio of the GMT sheet can be calculated from equations (6) and (7), respectively.(6)E1=Ep⋅Vp+Eb⋅Vb,(7)μ1=μp⋅Vp+μb⋅Vb,where Ep and Eb are the elastic modulus of polypropylene and glass fiber, respectively. Moreover, μp and μb are Poisson’s ratio of polypropylene and glass fiber, respectively. The upper layer is a single material of needle-punched PET fabric, so the elastic modulus and Poisson’s ratio are E2 and μ2. Then the elastic modulus of the two-layer composite structure is divided into two layers:(8)Ez=E1,−h1,h2,E2,h2,h3.The same operation can be conducted for Poisson’s ratio in the following form:(9)μz=μ1,−h1,h2,μ2,h2,h3.The differential equation of the elastic surface of the thin plate can be expressed in the following form:(10)D∇4w=q,where D is the bending stiffness of the sheet, D=Eh3/121−μ2, and h is the total thickness of the sheet. The left side of equation (10) D∇4w denotes the elastic resistance per unit area, and w is the deflection. Moreover, the left side q is the lateral load per unit area. The stress differential equation can be simplified as follows:(11)∂σZσZ=Ez2⋅1−μz2⋅h24−z2⋅∇4w.Integrating both sides with respect toz from −h1 to h3 yields the following expressions:(12)∫−h1h3∂σZσZdz=σZ|h3−h1=0−q=−q,∫−h1h3Ez21−μz2h24−z2∇4wdz=∫−h1h3Ez21−μz2h24−z2dz∇4w.In order to simplify expressions, parameterD is defined as follows:(13)D=−∫−h1h3Ez21−μz2h24−z2dz.Then the integral can be decomposed into the sum of integrals of each layer:(14)−D=∫−h1h2Ez21−μz2h24−z2dz+∫h2h3Ez21−μz2h24−z2dz=−D1−D2.


## 2.2. The Mechanical Model of the Composite Structure Is Solved

When solving the thin plate deflection-bending problem by the Navier method, the deflectionw should be initially obtained from the elastic surface differential equation under the boundary conditions of the plate edge [9].During the unsupported subsidence, the length and width for a rectangular thin plate, which is simply supported on four sides, area and b, respectively. In this case, the boundary conditions are(15)wx=0=0,∂2w∂x2x=0=0wx=a=0,∂2w∂x2x=a=0wy=0=0,∂2w∂x2y=0=0wy=b=0,∂2w∂x2y=b=0.The deflectionw can be calculated from the following trigonometric series:(16)w=∑m=1∞∑n=1∞wm=∑m=1∞∑n=1∞Amnsinmπxasinnπyb,where m and n are positive integers. Substitute this expression into equation (15) which indicates that all boundary conditions are satisfied. In order to find the coefficient Amn, equation (16) is introduced into equation (10):(17)π4D∑m=1∞∑n=1∞m2a2+n2b2Amnsinmπxasinnπyb=q,and expanding the load q on the right side of equation (17) into a multiple triangle series results in the following expression:(18)q=∑m=1∞∑n=1∞Amnsinmπxasinnπyb,where A can be solved through the method of determining the triangle series, and the solution is(19)Amn=4ab∫0a∫0bqx,ysinnπybsinmπxbdxdy.Then equation (18) can be rewritten in the following form:(20)q=4ab∑m=1∞∑m=1∞4ab∫0a∫0bqx,ysinnπybsinmπxbdxdy⋅sinnπybsinmπxb,(21)Amn=4∫0a∫0bqx,ysinnπy/bsinmπx/bdxdyπ4abDm2/a2+n2/b2.When the thin plate is subjected to a uniform load,q becomes a constant q0. Under this circumstance, equation (21) can be rewritten in the following form:(22)Amn=16q0π6mnDm2/a2+n2/b22,m=1,3,5,…;n=1,3,5,….Then equation (16) can be applied to calculate the deflection:(23)w=16q0π2D∑m=1,3,5,…∞∑m=1,3,5,…∞sinmπx/asinnπy/bmnm2/a2+n2/b22.Internal forces can be expressed in the following form:(24)Mx=−D∂2w∂x2+μ∂2w∂y2,My=−D∂2w∂y2+μ∂2w∂x2,Mxy=Myx=−D1−μ∂2w∂x∂y,FS,x=−D1−μ∂2∂x∇2w,FS,y=−D1−μ∂2∂y∇2w.Based on the performed calculations and analyses, the two-layer composite structure is simply supported on four sides, where the lengtha, width b, and GMT sheet thickness are 0.1 m, 0.1 m, and 3.5 mm, respectively. Moreover, the composite density is 1.3 g/cm3 and the void ratio is assumed to be 0. The mass fraction and density of the glass fiber are 30% and 2.6 g/cm3, respectively. In the present study, alkali-free fibers (E) with an elastic modulus of 72 GPa and Poisson’s ratio of 0.2 are applied. Furthermore, the density, elastic modulus, and Poisson’s ratio of PP are 0.93 g/cm3, 1.4 GPa, and 0.4, respectively. The lateral uniform load is 10 MPa, the thickness of the transition layer is 1.5 mm, and the specimen is made of needle-punched PET with an elastic modulus of 4000 MPa.Based on the established mathematical model, the thickness impact on the mechanical properties of the thin plate can be investigated by changing the thickness of the GMT sheet under the condition that the glass fiber content is maintained 30%, while the applied load and constraints remain constant. The iterative calculations are performed in the MATLAB environment. The maximum deflection and stress value of the two-layer structure are shown in Table1. The stress value can be obtained from the internal force value through the generalized Hook definition of the thin plate.Table 1
Maximum deflection and stress values for different thicknesses.
GMT plate thickness (mm)Deflection (m)σx, σy (Pa)τxy, τyx (Pa)3.30.002289934.91424.5983.40.001938030.42821.4783.50.001654726.35718.7593.60.001424223.45816.1023.70.001198419.85713.478Table1 indicates that, for a constant load and constraint, as the thickness of the GMT sheet increases, the deflection of the two-layer composite structure gradually decreases, thereby reducing the normal stress and the shear force. It is inferred that, as the sheet thickness increases, flexural, tensile, and shear resistance of the whole plate improve.When the thickness is fixed and different glass fiber ratios are used, the deflection and internal force of the plate are shown in Table2.Table 2
Deflection and internal force values for different glass fiber contents.
GMT plate glass fiber content (%)Deflection (mm)Bending moment (N⋅m)Torque (N⋅m)250.001673225.43429.562300.001654725.29329.646350.001623825.17729.739400.001596625.05429.822450.001562424.96429.018Table2 indicates that, under the condition of unchanged load and restraint, as the content of glass fiber increases, the deflection value and bending moment value decrease, while the corresponding torque increases. Therefore, the optimal content of GMT sheet glass fiber is selected according to the required mechanical performance index.


## 2.3. Thickness of the Car Coat Rack

Based on the foregoing discussions, it is concluded that the rigidity design requirements can be achieved by increasing the thickness of the GMT sheet without changing the main structure of the car coat rack. However, the material consumption for too thick products cannot be justified. More specifically, the production cost increases, the product becomes prone to stress concentration, and the molding quality reduces.The stiffness test principle is as follows: as a support frame, the sample deflects continuously at the midpoint of the span until the fracture or deformation reaches a certain predefined value. In this process, the applied stress to the sample is measured. Figure2 illustrates the schematic configuration of composite layers, indicating that the composite consists of a PET fabric and the GMT sheet.Figure 2
Schematic configurations of the composite layers.Since the fabric and the substrate are subjected to stress in parallel, the mixed law of elastic modulus of the parallel composite structure can be used [9]:(25)E0=E1V1+E2V2,where E0, E1, and E2 denote the elastic modulus of the composite material, PP glass fiberboard, and PET fabric, respectively. Substituting equations (6) into (25) gives(26)E0=EpVp+EbVbV1+E2V2.Then the predicted value of the elastic modulus in the composite shown in Figure2 can be expressed in the following form:(27)E0=EpVp+EbVbd1d1+d2+E2d2d1+d2,where d1 and d2 denote the thickness of PP glass fiber and PET fabric, respectively. In order to meet the rigidity design requirements, the following conditions should be satisfied:(28)E0I≥σ,I=b⋅d1+d2312,where σ is the required stiffness; I is the moment of inertia; and b is the width of the coat rack. In order to solve the equations easily, a thick needle-punched PET fabric with a thickness of 1.5 mm, elastic modulus of 4000 MPa is considered in the calculations, and GMT sheet porosity is assumed to be 0. The glass fiber is made of alkali-free fiber (E) with an elastic modulus of 72 GPa, while the elastic modulus of PP is 1.4 GPa. Moreover, the coat rack width b is 600 mm and the required bending stiffness is σ ≥ 60 N/mm2.The combination of glass fiber content and thickness of the GMT sheet that satisfy the stiffness requirements can be calculated from equations (27) and (28). The proposed combination is shown in Table 3.Table 3
Proposed GMT sheet glass fiber ratio and thickness combination.
GMT sheet glass fiber content (%)GMT sheet thickness (mm)Round thickness (mm)204.8865254.6364.8304.3864.5354.1364.3403.8864Based on the obtained results, a glass fiber content of 30% and a thickness of 4.5 mm is proposed as the selected composition. To maximize the safety factor, the thickness of the GMT sheet is rounded to 4.5 mm. Since the thicknesses of the GMT sheet and PET fabric are 4.5 mm and 1.5 mm, the overall thickness of the composite is 6 mm.


## 2.4. Formulation of the Process Parameters to Mold the Main Body of the Automobile Coat Rack

By analyzing the common quality defects and causes of the main body of the coat rack, the main processing parameters that affect the molding process of the car coat rack are determined. These parameters are the molding pressure, holding time, and the molding temperature.


### 2.4.1. Sheet Warm-Up Time

In order to heat the GMT sheet, hot air convection preheating is used. In this case, the sheet can be expanded to 2-3 times the thickness of the cold sheet. It should be indicated that the sheet temperature should be less than 230°C; otherwise, the PP polypropylene will degrade, thereby affecting the surface of the molded product and reducing the heat resistance and strength of the product [10, 11]. The preheating time of the GMT sheet can be calculated through the following expression [12] (29):(29)t=70.9521+0.2318d−0.2810T,where t is the appropriate preheating time, min; d is the sheet thickness, mm; and T is the hot air temperature, °C.During the experiment, hot air flow from an oven is used to preheat the sheet. The thickness of the GMT sheet is 4.5 mm and the preheating temperature of the hot air is 250°C. Based on equation (29), the preheating time is adjusted to 99 s. In the experiment, the sheet is taken out of the oven and transferred to the mold. Since the temperatures of the material and the oven reduce during the material transfer, the preheating time of the sheet is set to 110 s.


### 2.4.2. Molding Pressure and Mold Clamping Speed

Appropriate molding pressure not only improved the product quality but also made up for the shortage if other industrial parameters were not adequate [13].(1)
The method of maximum gauge pressure:(30)ftable=f1f2ST,whereftable, f1, and f2 denote the gauge pressure, unit pressure required by the product, and the rated gauge pressure of the press, respectively. Meanwhile, S and T denote the projected area of the product in the compression direction and the pressing force, respectively.(31)x=f1ST,wherex is the ratio of the force acting on the pressure direction of the product to the pressing force. When x > 1, the pressing force is too small to meet the production requirements. On the other hand, when x < 1, the pressing force is sufficient to meet the production demand.(2)
The method of Piston area method:(32)ftable=f1SF,whereF is the cross-sectional area of the piston of the press.In the present study, for a four-post-servo-hydraulic machine (model; company) withF = 0.5 m2 and S = 0.84 m2, the molding pressure is 13 MPa and ftable= 22 MPa. For these parameters, equation (31) indicates that x > 1, so that the selected hydraulic pressure meets the molding requirements.In the molding process to produce the main body of the automobile coat rack, the mold clamping speed is not a fixed value. It should be indicated that the molding process can be divided into 5 main stages:①
When the sheet is placed in the mold, the concave mold decreases at a speed of 60 mm/s.②
The interval between the convex and concave mold gradually decreases. Consequently, when the guide column enters the guide sleeve, the speed becomes 8 mm/s. Moreover, when the convex and concave molds are closed, the sheet is molded under pressure.③
Cooling and solidification.④
After holding the pressure for 80 s, the concave mold is split at a speed of 8 mm/s.⑥
After the split mold is completed, the concave mold is 60 mm/s. Then the speed rises and when the die is reset, the product is removed.These five stages are also shown in Figure3.Figure 3
Speed variation graph during the mold clamping.


### 2.4.3. Holding Time

The cooling time calculation results are utilized as a reference to select the holding time. For large products such as car coat racks, the length and width are far beyond their thickness, and they have two-dimensional characteristics. However, the heat is only transferred in the one-dimensional direction, that is, the thickness direction of the part. Therefore, the cooling of such products can be conducted as follows.Heat transfer can be calculated from the following expression:(33)∂T∂θ′=k∂2T∂2x2,where T and k denote the average temperature when the sheet is demolded and the thermal diffusion coefficient of the sheet, respectively.In order to separate the variables, the following equation is used:(34)θ′=−2π2x12klnπT−TwT0−Tw,k=λρ⋅cp,where x1 is half the thickness of the part, m; Tw is the mold cavity temperature, °C; T0 is the sheet molding temperature, °C; ρ is the density of the sheet, kg/m3; λ is the thermal conductivity of the sheet, kW/m⋅°C; and cp is the specific heat capacity of the sheet, kJ/kg⋅°C.The process parameter values are shown in Table4.Table 4
The value of process parameters.
Parameter nameValueThe thickness of the main body of the car coat rack (mm)4.5The average temperature of the mold cavity wallTw (°C)40The molding temperatureT0 (°C)220The average temperature when the sheet is demoldedT (°C)80The density of the sheetρ (kg/m3)1130The specific heat capacity of the sheetcp (kJ/kg·°C)1.55The thermal conductivity of the sheetλ (kW/m·°C)0.33 × 10−3Substituting the abovementioned parameters into equation (34), the cooling time of the molding of the main body of the automobile coat rack is mathematically expressed as follows:(35)k=0.33×10−31130×1550×10−3=1.88×10−7,t=−2π20.002252klnπ80−404220−40≈74s.Combining with the previous pressure holding time, the pressure holding time is determined to be 80 s.


## 2.4.1. Sheet Warm-Up Time

In order to heat the GMT sheet, hot air convection preheating is used. In this case, the sheet can be expanded to 2-3 times the thickness of the cold sheet. It should be indicated that the sheet temperature should be less than 230°C; otherwise, the PP polypropylene will degrade, thereby affecting the surface of the molded product and reducing the heat resistance and strength of the product [10, 11]. The preheating time of the GMT sheet can be calculated through the following expression [12] (29):(29)t=70.9521+0.2318d−0.2810T,where t is the appropriate preheating time, min; d is the sheet thickness, mm; and T is the hot air temperature, °C.During the experiment, hot air flow from an oven is used to preheat the sheet. The thickness of the GMT sheet is 4.5 mm and the preheating temperature of the hot air is 250°C. Based on equation (29), the preheating time is adjusted to 99 s. In the experiment, the sheet is taken out of the oven and transferred to the mold. Since the temperatures of the material and the oven reduce during the material transfer, the preheating time of the sheet is set to 110 s.


## 2.4.2. Molding Pressure and Mold Clamping Speed

Appropriate molding pressure not only improved the product quality but also made up for the shortage if other industrial parameters were not adequate [13].(1)
The method of maximum gauge pressure:(30)ftable=f1f2ST,whereftable, f1, and f2 denote the gauge pressure, unit pressure required by the product, and the rated gauge pressure of the press, respectively. Meanwhile, S and T denote the projected area of the product in the compression direction and the pressing force, respectively.(31)x=f1ST,wherex is the ratio of the force acting on the pressure direction of the product to the pressing force. When x > 1, the pressing force is too small to meet the production requirements. On the other hand, when x < 1, the pressing force is sufficient to meet the production demand.(2)
The method of Piston area method:(32)ftable=f1SF,whereF is the cross-sectional area of the piston of the press.In the present study, for a four-post-servo-hydraulic machine (model; company) withF = 0.5 m2 and S = 0.84 m2, the molding pressure is 13 MPa and ftable= 22 MPa. For these parameters, equation (31) indicates that x > 1, so that the selected hydraulic pressure meets the molding requirements.In the molding process to produce the main body of the automobile coat rack, the mold clamping speed is not a fixed value. It should be indicated that the molding process can be divided into 5 main stages:①
When the sheet is placed in the mold, the concave mold decreases at a speed of 60 mm/s.②
The interval between the convex and concave mold gradually decreases. Consequently, when the guide column enters the guide sleeve, the speed becomes 8 mm/s. Moreover, when the convex and concave molds are closed, the sheet is molded under pressure.③
Cooling and solidification.④
After holding the pressure for 80 s, the concave mold is split at a speed of 8 mm/s.⑥
After the split mold is completed, the concave mold is 60 mm/s. Then the speed rises and when the die is reset, the product is removed.These five stages are also shown in Figure3.Figure 3
Speed variation graph during the mold clamping.


## 2.4.3. Holding Time

The cooling time calculation results are utilized as a reference to select the holding time. For large products such as car coat racks, the length and width are far beyond their thickness, and they have two-dimensional characteristics. However, the heat is only transferred in the one-dimensional direction, that is, the thickness direction of the part. Therefore, the cooling of such products can be conducted as follows.Heat transfer can be calculated from the following expression:(33)∂T∂θ′=k∂2T∂2x2,where T and k denote the average temperature when the sheet is demolded and the thermal diffusion coefficient of the sheet, respectively.In order to separate the variables, the following equation is used:(34)θ′=−2π2x12klnπT−TwT0−Tw,k=λρ⋅cp,where x1 is half the thickness of the part, m; Tw is the mold cavity temperature, °C; T0 is the sheet molding temperature, °C; ρ is the density of the sheet, kg/m3; λ is the thermal conductivity of the sheet, kW/m⋅°C; and cp is the specific heat capacity of the sheet, kJ/kg⋅°C.The process parameter values are shown in Table4.Table 4
The value of process parameters.
Parameter nameValueThe thickness of the main body of the car coat rack (mm)4.5The average temperature of the mold cavity wallTw (°C)40The molding temperatureT0 (°C)220The average temperature when the sheet is demoldedT (°C)80The density of the sheetρ (kg/m3)1130The specific heat capacity of the sheetcp (kJ/kg·°C)1.55The thermal conductivity of the sheetλ (kW/m·°C)0.33 × 10−3Substituting the abovementioned parameters into equation (34), the cooling time of the molding of the main body of the automobile coat rack is mathematically expressed as follows:(35)k=0.33×10−31130×1550×10−3=1.88×10−7,t=−2π20.002252klnπ80−404220−40≈74s.Combining with the previous pressure holding time, the pressure holding time is determined to be 80 s.


## 2.5. Shrinkage Compensation Value

The molding shrinkage rate of the product is measured experimentally. The experimental process is as follows: the main model of the coat rack of the reference model is selected. Moreover, 20 pieces of GMT plates with a thickness of 4.5 mm glass fiber content of 30% and a needle-punched PET fabric of 1.5 mm thickness are selected. It should be indicated that, in the experiment, the sheet preheating time, the pressure holding time, molding pressure, and mold closing speed utilize the abovementioned parameters. When the mold is opened, the distance between the selected positioning points is measured, and each piece is averaged at one measurement in the length and width directions. After the product is left at room temperature for 24 hours, the selected positioning point is measured again, using equation (29). Finally, the average shrinkage rate of the composite two-layer structure of GMT sheet + needle-punched PET fabric is 0.354%.


## 2.6. Heat Transfer Calculation for the Cooling System of the Automobile Coat Rack Compression Molding Mold

The cooling system adopts the water cooling method, and a cooling water channel is laid on the mold surface to pass the cooling water. In this section, the design and calculation of the cooling system of the forming mold of the main body of the car coat rack are established.


### 2.6.1. The Total Heat of the Sheet into the Cavity

In a molding process, the total heat transferred from the sheet to the mold is mathematically expressed as follows:(36)Q1=Cp⋅T0−T⋅V⋅ρ,where T0 and T denote the sheet molding temperature and the average temperature when the sheet is demolded, respectively. Moreover, ρ, cp, and V are the density of the sheet, the specific heat capacity of the sheet, and the volume of the sheet, respectively. It is worth noting that the volume of the sheet is set to 1.4 m × 0.6 m × 0.0045 m.


### 2.6.2. Heat Taken away by the Cooling Water

The heat taken away by the mold cooling water in one press molding is expressed as follows:(37)Q2=Q1−Qc−QR−QL,where Qc, QR, and QL denote the convection heat transfer from the mold to the air, the radiation heat transfer from the mold to the air, and the heat transfer from the mold to the hydraulic machine through the upper and lower mold bases, respectively.Under normal circumstances, 90%–95% of the heat brought into the sheet can be taken away by the cooling water through the mold cooling channel. Therefore, during the design process,Q2 = Q1.Since the main body of the car coat rack is a large thin-walled product, it can be designed according to 50% of the heat taken by the concave mold and the convex mold; then,(38)Q2a=Q2t=0.5Q2.


### 2.6.3. Cooling Water Hole Diameter

The volume flow of the cooling water can be calculated according to equation (39) during the molding process of the die and the die of the die in one molding process:(39)V2a=Q2aρ1c1θ1−θ2⋅t,V2t=Q2tρ1c1θ1−θ2⋅t,where θ1 and θ2 are the inlet and outlet temperatures of the mold cooling water, respectively. It should be indicated that the temperature difference should not be too large and its maximum value should not exceed 5°C. Moreover, V (m3⋅min−1) and ρ1 (kg/m3) denote the volume flow of cooling water and the density of water, respectively. Furthermore, c1 (kJ/kg⋅°C) and t (min) are the specific heat capacity of water and the cooling time of one molding, respectively.The density of water is 1000 kg/m3, and the specific heat capacity of water is 4.2kJ/kg⋅°C, which is calculated according to equations (36) and (38):(40)V2a=V2t=0.5Q2ρ1c1θ1−θ2⋅t=0.5⋅Cp⋅T0−T⋅V⋅ρρ1c1θ1−θ2⋅t=16.275m3/min.According to the correlation between the cooling water pipe volume flow and the pipe diameter, the cooling water pipe diameter is determined asD = 27 mm.The correlation between the distanceP of the cooling water pipe, the distance H from the profile, and the diameter D of the channel is P = (3∼5) D, H = (1∼2) D; from D = 27 mm, take the pipe distance P = 80 mm, H = 30 mm.


## 2.6.1. The Total Heat of the Sheet into the Cavity

In a molding process, the total heat transferred from the sheet to the mold is mathematically expressed as follows:(36)Q1=Cp⋅T0−T⋅V⋅ρ,where T0 and T denote the sheet molding temperature and the average temperature when the sheet is demolded, respectively. Moreover, ρ, cp, and V are the density of the sheet, the specific heat capacity of the sheet, and the volume of the sheet, respectively. It is worth noting that the volume of the sheet is set to 1.4 m × 0.6 m × 0.0045 m.


## 2.6.2. Heat Taken away by the Cooling Water

The heat taken away by the mold cooling water in one press molding is expressed as follows:(37)Q2=Q1−Qc−QR−QL,where Qc, QR, and QL denote the convection heat transfer from the mold to the air, the radiation heat transfer from the mold to the air, and the heat transfer from the mold to the hydraulic machine through the upper and lower mold bases, respectively.Under normal circumstances, 90%–95% of the heat brought into the sheet can be taken away by the cooling water through the mold cooling channel. Therefore, during the design process,Q2 = Q1.Since the main body of the car coat rack is a large thin-walled product, it can be designed according to 50% of the heat taken by the concave mold and the convex mold; then,(38)Q2a=Q2t=0.5Q2.


## 2.6.3. Cooling Water Hole Diameter

The volume flow of the cooling water can be calculated according to equation (39) during the molding process of the die and the die of the die in one molding process:(39)V2a=Q2aρ1c1θ1−θ2⋅t,V2t=Q2tρ1c1θ1−θ2⋅t,where θ1 and θ2 are the inlet and outlet temperatures of the mold cooling water, respectively. It should be indicated that the temperature difference should not be too large and its maximum value should not exceed 5°C. Moreover, V (m3⋅min−1) and ρ1 (kg/m3) denote the volume flow of cooling water and the density of water, respectively. Furthermore, c1 (kJ/kg⋅°C) and t (min) are the specific heat capacity of water and the cooling time of one molding, respectively.The density of water is 1000 kg/m3, and the specific heat capacity of water is 4.2kJ/kg⋅°C, which is calculated according to equations (36) and (38):(40)V2a=V2t=0.5Q2ρ1c1θ1−θ2⋅t=0.5⋅Cp⋅T0−T⋅V⋅ρρ1c1θ1−θ2⋅t=16.275m3/min.According to the correlation between the cooling water pipe volume flow and the pipe diameter, the cooling water pipe diameter is determined asD = 27 mm.The correlation between the distanceP of the cooling water pipe, the distance H from the profile, and the diameter D of the channel is P = (3∼5) D, H = (1∼2) D; from D = 27 mm, take the pipe distance P = 80 mm, H = 30 mm.


## 3. Simulation of the Molding Process

In the present study, the Ls-Dyna finite element program is used to simulate the molding process of the main body of the car coat rack to obtain the nodal force information of the stress concentration area of the mold during the molding process [14]. Therefore, the boundary conditions for the topological optimization of the convex and concave mold are obtained. Moreover, the design of the mold is optimized.


### 3.1. The Main Influencing Factors of the Modeling Process

The molding process is a complicated nonlinear challenging process. The influencing factors in this process include the compression molding of the sheet, the elastic deformation of the mold under high temperature and high pressure, and the frictional contact between the mold and the sheet during the mold clamping process; it is a complex finite element problem of thermal and mechanical coupling [15].


#### 3.1.1. Grid Division

The density distribution of the grid is one of the important factors that affect the calculation accuracy [16]. In order to simplify the meshing and improve the quality and solution rate of meshing, the CATIA software is utilized to simplify the model in the early stage. The convex and concave molds use sweeping meshing, SOLID168 elements are used to mesh, and the length direction is divided into 60. It should be indicated that the width direction is divided into 3 equal parts with a total of 942,000 units. The sheet adopts shell163 thin shell unit and the sheet grid is a regular quadrilateral unit with 15000 units. Moreover, the unit has 11 different algorithms. In the present study, the Belytschko−Tsay algorithm is used because its calculation speed is high. Figure 4 illustrates the grid division diagram.Figure 4
Grid division diagram.


#### 3.1.2. Heat Transfer Model Settings

According to the characteristics of the molding process, the following two heat transfer methods are described:(1)
Radiation and air-cooling heat exchangeOwing to the small proportion of the air-cooled heat transfer, air-cooling and radiation can be written as a unified heat transfer equation [17]:(41)q=h1+h2T−T0,whereq is the heat flux, density, J/m2⋅s. Moreover, h1 and h2denote the convective heat transfer coefficient, W/m2⋅K, and the radiant heat transfer coefficient, W/m2⋅K, respectively.(42)h2=εδT+T0T2−T02,whereε and δ are the blackness coefficient and the Boltzmann constant, J/K, respectively. Moreover, T and T0 are the surface temperature of the sheet, °C, and the ambient temperature, °C, respectively.Table5 presents the specific values.(2)
Contact heat transferA contact heat transfer exists between the sheet and the mold, and calculation of the contact heat transfer is challenging. Empirical equation [18] (43) is used in this regard:(43)q=hcondT−T1.Table 5
Heat transfer constant.
h1 + h2εδTT0200 W/(m2·K)0.65.672×10−8J/K220°C20°CIn the abovementioned equation,hcond=K/lgap, where K = 10 and lgap is the contact gap between the workpiece and the mold, which is automatically calculated in the LS-DYNA program, and its unit is mm. Moreover, T1 denotes the mold cavity temperature, °C.


#### 3.1.3. Contact and Friction Settings

The keyword∗CONTACT_FORMING_ONE_WAY_SURFACE_TO_SURFACE_THERMAL defines the contact method of the compression molding, and FORMING_ONE_WAY is the type of contact specifically used for molding. It should be indicated that the penalty function algorithm is used to calculate the contact interface force. Since LS-DYNA is originally used for the collision analysis, if it is used for the forming analysis, the default parameter settings should be modified [19]. In the present study, the penalty function stiffness factor SLSFAC is modified to 0.1.In the compression molding process, there is friction between the mold and the sheet. Therefore, the friction coefficient should be determined. In LS-DYNA, the friction force is calculated by the Coulomb model. In this study, the average friction coefficient is 0.1.


#### 3.1.4. Hourglass Control

In order to verify the validity of the calculation result, the ratio of the hourglass energy to total energy should be less than 10%. On the other hand, the hourglass control coefficient cannot be too large. When the hourglass control coefficient is greater than 0.15, the calculation will be unstable. Therefore, the hourglass control sets IHQ = 4, and the hourglass energy coefficient is set to QH = 0.1.


#### 3.1.5. Energy Dissipation Control

In the numerical simulation of the compression molding, the main reference index to determine the accuracy of the simulation result is the change in energy. The control of energy is set by the keyword∗CONTROL_ENERGY. In this study, set HGEN = 2, RWEN = 2, SLNTEN = 2, and RYLEN = 2. The purpose of the abovementioned setting is to consider various energy dissipations in the energy balance.


#### 3.1.6. Control of the Time Step

In the LS-DYNA software, the convergence of the time integration algorithm is conditional, and it is stable only when the time step Δt is less than the critical time step, which is mathematically expressed as(44)Δt<Δlc+x˙,c=E1−v2⋅ρ,where Δl, c, and x denote the characteristic length of the element, m, the sound velocity of the material, m/s, and the velocity of the node, m/s, respectively.In the compression molding simulation, as time passes, the deformation increases. Then, Δt quickly decreases, indicating that the integration algorithm divergence interrupts the simulation. In order to avoid this phenomenon, the mass scaling technique is used to set a duration value ∆tsz through the keyword ∗CONTROL_TIMESTEP. In the simulation, if ∆t < ∆tsz, the program automatically increases the density of the structural material and the time step size. In the present study, Δtsz is set as 0.0000025, time step size scaling factor TSSFAC = 0.9, and DTINIT = 0. Moreover, the initial step size is set by the program.


### 3.2. Setting of Materials and Boundary Conditions

(1)
Material: the sheet material is selected from the constitutive equation of the thermoplastic material. The thickness of the sheet is defined by the actual thickness of the part. The material of the convex and concave mold is 45#. The sheet is made of PP-based GMT composite material with a glass fiber content of 30%.(2)
Initial boundary conditions: the initial process parameters of the automobile coat rack compression molding are set as follows. Moreover, the punch is set to fix the lower table of the hydraulic press, and the concave mold should be closed under the pressure of the press. The molding pressure is set to 13 MPa, the mold closing speed after the upper mold contacts the sheet is 8 mm/s, and the initial temperature of the sheet is 220°C. The mold temperature is 40°C, the cooling channel water temperature is 20°C, and the pressure holding time is 80 s.


### 3.3. Finite Element Numerical Calculation Results

Figures5–7 show the numerical simulation results in the molding process, which are obtained after the software calculation.Figure 5
The equivalent stress diagram of the forming process.
(a)(b)(c)(d)(e)(f)Figure 6
Cloud diagram of stress changes during the packing process.
(a)(b)(c)(d)(e)(f)Figure 7
The main body shape of the car coat rack.
(a)(b)


#### 3.3.1. Equivalent Stress during Forming

The simulation results of the equivalent stress during the molding process show that, during the molding process of the main body of the automobile coat rack, the stress concentration area is located at the outer contour of the mold and the slope is steeper.


#### 3.3.2. Equivalent Stress during the Holding Pressure

The simulation results of the equivalent stress during the packing process show that the node displacement is the largest when the packing time is 11 s. As the packing time passes, the mold temperature decreases and the deformation gradually decreases too. After 64 s, it remains unchanged until the end of the packing process.


#### 3.3.3. Appearance of the Molding Process

The abovementioned simulation results show that the nodal force information of the mold during the compression molding process is obtained. Therefore, the load mapping is used as the boundary condition for the next topology optimization.


## 3.1. The Main Influencing Factors of the Modeling Process

The molding process is a complicated nonlinear challenging process. The influencing factors in this process include the compression molding of the sheet, the elastic deformation of the mold under high temperature and high pressure, and the frictional contact between the mold and the sheet during the mold clamping process; it is a complex finite element problem of thermal and mechanical coupling [15].


### 3.1.1. Grid Division

The density distribution of the grid is one of the important factors that affect the calculation accuracy [16]. In order to simplify the meshing and improve the quality and solution rate of meshing, the CATIA software is utilized to simplify the model in the early stage. The convex and concave molds use sweeping meshing, SOLID168 elements are used to mesh, and the length direction is divided into 60. It should be indicated that the width direction is divided into 3 equal parts with a total of 942,000 units. The sheet adopts shell163 thin shell unit and the sheet grid is a regular quadrilateral unit with 15000 units. Moreover, the unit has 11 different algorithms. In the present study, the Belytschko−Tsay algorithm is used because its calculation speed is high. Figure 4 illustrates the grid division diagram.Figure 4
Grid division diagram.


### 3.1.2. Heat Transfer Model Settings

According to the characteristics of the molding process, the following two heat transfer methods are described:(1)
Radiation and air-cooling heat exchangeOwing to the small proportion of the air-cooled heat transfer, air-cooling and radiation can be written as a unified heat transfer equation [17]:(41)q=h1+h2T−T0,whereq is the heat flux, density, J/m2⋅s. Moreover, h1 and h2denote the convective heat transfer coefficient, W/m2⋅K, and the radiant heat transfer coefficient, W/m2⋅K, respectively.(42)h2=εδT+T0T2−T02,whereε and δ are the blackness coefficient and the Boltzmann constant, J/K, respectively. Moreover, T and T0 are the surface temperature of the sheet, °C, and the ambient temperature, °C, respectively.Table5 presents the specific values.(2)
Contact heat transferA contact heat transfer exists between the sheet and the mold, and calculation of the contact heat transfer is challenging. Empirical equation [18] (43) is used in this regard:(43)q=hcondT−T1.Table 5
Heat transfer constant.
h1 + h2εδTT0200 W/(m2·K)0.65.672×10−8J/K220°C20°CIn the abovementioned equation,hcond=K/lgap, where K = 10 and lgap is the contact gap between the workpiece and the mold, which is automatically calculated in the LS-DYNA program, and its unit is mm. Moreover, T1 denotes the mold cavity temperature, °C.


### 3.1.3. Contact and Friction Settings

The keyword∗CONTACT_FORMING_ONE_WAY_SURFACE_TO_SURFACE_THERMAL defines the contact method of the compression molding, and FORMING_ONE_WAY is the type of contact specifically used for molding. It should be indicated that the penalty function algorithm is used to calculate the contact interface force. Since LS-DYNA is originally used for the collision analysis, if it is used for the forming analysis, the default parameter settings should be modified [19]. In the present study, the penalty function stiffness factor SLSFAC is modified to 0.1.In the compression molding process, there is friction between the mold and the sheet. Therefore, the friction coefficient should be determined. In LS-DYNA, the friction force is calculated by the Coulomb model. In this study, the average friction coefficient is 0.1.


### 3.1.4. Hourglass Control

In order to verify the validity of the calculation result, the ratio of the hourglass energy to total energy should be less than 10%. On the other hand, the hourglass control coefficient cannot be too large. When the hourglass control coefficient is greater than 0.15, the calculation will be unstable. Therefore, the hourglass control sets IHQ = 4, and the hourglass energy coefficient is set to QH = 0.1.


### 3.1.5. Energy Dissipation Control

In the numerical simulation of the compression molding, the main reference index to determine the accuracy of the simulation result is the change in energy. The control of energy is set by the keyword∗CONTROL_ENERGY. In this study, set HGEN = 2, RWEN = 2, SLNTEN = 2, and RYLEN = 2. The purpose of the abovementioned setting is to consider various energy dissipations in the energy balance.


### 3.1.6. Control of the Time Step

In the LS-DYNA software, the convergence of the time integration algorithm is conditional, and it is stable only when the time step Δt is less than the critical time step, which is mathematically expressed as(44)Δt<Δlc+x˙,c=E1−v2⋅ρ,where Δl, c, and x denote the characteristic length of the element, m, the sound velocity of the material, m/s, and the velocity of the node, m/s, respectively.In the compression molding simulation, as time passes, the deformation increases. Then, Δt quickly decreases, indicating that the integration algorithm divergence interrupts the simulation. In order to avoid this phenomenon, the mass scaling technique is used to set a duration value ∆tsz through the keyword ∗CONTROL_TIMESTEP. In the simulation, if ∆t < ∆tsz, the program automatically increases the density of the structural material and the time step size. In the present study, Δtsz is set as 0.0000025, time step size scaling factor TSSFAC = 0.9, and DTINIT = 0. Moreover, the initial step size is set by the program.


## 3.1.1. Grid Division

The density distribution of the grid is one of the important factors that affect the calculation accuracy [16]. In order to simplify the meshing and improve the quality and solution rate of meshing, the CATIA software is utilized to simplify the model in the early stage. The convex and concave molds use sweeping meshing, SOLID168 elements are used to mesh, and the length direction is divided into 60. It should be indicated that the width direction is divided into 3 equal parts with a total of 942,000 units. The sheet adopts shell163 thin shell unit and the sheet grid is a regular quadrilateral unit with 15000 units. Moreover, the unit has 11 different algorithms. In the present study, the Belytschko−Tsay algorithm is used because its calculation speed is high. Figure 4 illustrates the grid division diagram.Figure 4
Grid division diagram.


## 3.1.2. Heat Transfer Model Settings

According to the characteristics of the molding process, the following two heat transfer methods are described:(1)
Radiation and air-cooling heat exchangeOwing to the small proportion of the air-cooled heat transfer, air-cooling and radiation can be written as a unified heat transfer equation [17]:(41)q=h1+h2T−T0,whereq is the heat flux, density, J/m2⋅s. Moreover, h1 and h2denote the convective heat transfer coefficient, W/m2⋅K, and the radiant heat transfer coefficient, W/m2⋅K, respectively.(42)h2=εδT+T0T2−T02,whereε and δ are the blackness coefficient and the Boltzmann constant, J/K, respectively. Moreover, T and T0 are the surface temperature of the sheet, °C, and the ambient temperature, °C, respectively.Table5 presents the specific values.(2)
Contact heat transferA contact heat transfer exists between the sheet and the mold, and calculation of the contact heat transfer is challenging. Empirical equation [18] (43) is used in this regard:(43)q=hcondT−T1.Table 5
Heat transfer constant.
h1 + h2εδTT0200 W/(m2·K)0.65.672×10−8J/K220°C20°CIn the abovementioned equation,hcond=K/lgap, where K = 10 and lgap is the contact gap between the workpiece and the mold, which is automatically calculated in the LS-DYNA program, and its unit is mm. Moreover, T1 denotes the mold cavity temperature, °C.


## 3.1.3. Contact and Friction Settings

The keyword∗CONTACT_FORMING_ONE_WAY_SURFACE_TO_SURFACE_THERMAL defines the contact method of the compression molding, and FORMING_ONE_WAY is the type of contact specifically used for molding. It should be indicated that the penalty function algorithm is used to calculate the contact interface force. Since LS-DYNA is originally used for the collision analysis, if it is used for the forming analysis, the default parameter settings should be modified [19]. In the present study, the penalty function stiffness factor SLSFAC is modified to 0.1.In the compression molding process, there is friction between the mold and the sheet. Therefore, the friction coefficient should be determined. In LS-DYNA, the friction force is calculated by the Coulomb model. In this study, the average friction coefficient is 0.1.


## 3.1.4. Hourglass Control

In order to verify the validity of the calculation result, the ratio of the hourglass energy to total energy should be less than 10%. On the other hand, the hourglass control coefficient cannot be too large. When the hourglass control coefficient is greater than 0.15, the calculation will be unstable. Therefore, the hourglass control sets IHQ = 4, and the hourglass energy coefficient is set to QH = 0.1.


## 3.1.5. Energy Dissipation Control

In the numerical simulation of the compression molding, the main reference index to determine the accuracy of the simulation result is the change in energy. The control of energy is set by the keyword∗CONTROL_ENERGY. In this study, set HGEN = 2, RWEN = 2, SLNTEN = 2, and RYLEN = 2. The purpose of the abovementioned setting is to consider various energy dissipations in the energy balance.


## 3.1.6. Control of the Time Step

In the LS-DYNA software, the convergence of the time integration algorithm is conditional, and it is stable only when the time step Δt is less than the critical time step, which is mathematically expressed as(44)Δt<Δlc+x˙,c=E1−v2⋅ρ,where Δl, c, and x denote the characteristic length of the element, m, the sound velocity of the material, m/s, and the velocity of the node, m/s, respectively.In the compression molding simulation, as time passes, the deformation increases. Then, Δt quickly decreases, indicating that the integration algorithm divergence interrupts the simulation. In order to avoid this phenomenon, the mass scaling technique is used to set a duration value ∆tsz through the keyword ∗CONTROL_TIMESTEP. In the simulation, if ∆t < ∆tsz, the program automatically increases the density of the structural material and the time step size. In the present study, Δtsz is set as 0.0000025, time step size scaling factor TSSFAC = 0.9, and DTINIT = 0. Moreover, the initial step size is set by the program.


## 3.2. Setting of Materials and Boundary Conditions

(1)
Material: the sheet material is selected from the constitutive equation of the thermoplastic material. The thickness of the sheet is defined by the actual thickness of the part. The material of the convex and concave mold is 45#. The sheet is made of PP-based GMT composite material with a glass fiber content of 30%.(2)
Initial boundary conditions: the initial process parameters of the automobile coat rack compression molding are set as follows. Moreover, the punch is set to fix the lower table of the hydraulic press, and the concave mold should be closed under the pressure of the press. The molding pressure is set to 13 MPa, the mold closing speed after the upper mold contacts the sheet is 8 mm/s, and the initial temperature of the sheet is 220°C. The mold temperature is 40°C, the cooling channel water temperature is 20°C, and the pressure holding time is 80 s.


## 3.3. Finite Element Numerical Calculation Results

Figures5–7 show the numerical simulation results in the molding process, which are obtained after the software calculation.Figure 5
The equivalent stress diagram of the forming process.
(a)(b)(c)(d)(e)(f)Figure 6
Cloud diagram of stress changes during the packing process.
(a)(b)(c)(d)(e)(f)Figure 7
The main body shape of the car coat rack.
(a)(b)


### 3.3.1. Equivalent Stress during Forming

The simulation results of the equivalent stress during the molding process show that, during the molding process of the main body of the automobile coat rack, the stress concentration area is located at the outer contour of the mold and the slope is steeper.


### 3.3.2. Equivalent Stress during the Holding Pressure

The simulation results of the equivalent stress during the packing process show that the node displacement is the largest when the packing time is 11 s. As the packing time passes, the mold temperature decreases and the deformation gradually decreases too. After 64 s, it remains unchanged until the end of the packing process.


### 3.3.3. Appearance of the Molding Process

The abovementioned simulation results show that the nodal force information of the mold during the compression molding process is obtained. Therefore, the load mapping is used as the boundary condition for the next topology optimization.


## 3.3.1. Equivalent Stress during Forming

The simulation results of the equivalent stress during the molding process show that, during the molding process of the main body of the automobile coat rack, the stress concentration area is located at the outer contour of the mold and the slope is steeper.


## 3.3.2. Equivalent Stress during the Holding Pressure

The simulation results of the equivalent stress during the packing process show that the node displacement is the largest when the packing time is 11 s. As the packing time passes, the mold temperature decreases and the deformation gradually decreases too. After 64 s, it remains unchanged until the end of the packing process.


## 3.3.3. Appearance of the Molding Process

The abovementioned simulation results show that the nodal force information of the mold during the compression molding process is obtained. Therefore, the load mapping is used as the boundary condition for the next topology optimization.


## 4. Optimized Mold Design

Currently, the lightweight design of molds is the development direction of mold optimization [20]. Based on the empirical design, CAE tools are used to perform mold structure analysis and optimization design before processing. Aiming at the problem of excessive quality of molds and due to the current empirical design of compression molding dies for automotive interior parts, in the present study, a lightweight design process and method is proposed for compression molding dies.


### 4.1. Lightweight Design Process of the Mold Structure

According to the numerical simulation method for the molding process of automotive interior parts and the topology optimization theory of the variable density method [21], a lightweight design method for the mold structure of interior parts molding is proposed.(1)
According to the mold design requirements, the structure of the compression molding convex and concave mold is designed in accordance with the traditional mold design guidelines.(2)
With the help of the LS-DYNA finite element program, the compression molding process of the main body of the coat rack is numerically simulated to obtain the force information of the convex and concave mold nodes in the actual working conditions.(3)
Load mapping is used to construct convex and concave boundary conditions of die structure optimization.(4)
The improvement space of the mold structure design is clarified, and the design area and nondesign area are determined.(5)
Based on the OptiStruct, the topological optimization of the convex-concave structure is iteratively solved.(6)
The final shape is optimized according to the topology of the convex and concave mold structure, and the three-dimensional software is utilized to redesign the mold structure.(7)
According to the redesigned mold structure for mold processing and manufacturing.


### 4.2. OptiStruct-Based Optimization Design of the Mold Topology

OptiStruct uses LS-DYNA for preprocessing and definition of structural optimization problems. After completing the finite element modeling in LS-DYNA, the optimization definition panel is used to define optimization variables, constraints and goals, and optimization parameters. Then, OptiStruct is submitted for structural analysis and optimization. Figure8 shows the OptiStruct structure optimization design flowchart [22].Figure 8
Topology optimization flowchart.


#### 4.2.1. The Establishment of the Topological Optimization Model of the Mold Convex and Concave Mold Structure

In the optimization process of the convex and concave dies, the problem of the mass minimization can be equated to the problem of the volume minimization. The real working conditions of the convex and concave molds are analyzed. Moreover, based on the solid isotropic material penalty (SIMP) method and according to the numerical simulation calculation results of the compression molding, the maximum displacement of the nodes in the stress concentration area is used as the constraint condition. Furthermore, the material element density is used as the design variable. The minimum volume of the convex-concave structure is the SIMP optimization model [23], and the density penalty factor P is set to 3. The model is formulas (45) and (46),(45)Ee=ρep⋅E0,0<ρmin≤ρe≤1,(46)min:Vsubjectto:KU=F,0<ρmin≤ρe≤1,Uj≤Uj∗,where ρi are V denote the design variable for the topology optimization with the variable density method and the overall structure volume, respectively. Moreover, U, K, and F are the overall displacement column vector, the overall structure stiffness matrix, and the external load vector, respectively.


#### 4.2.2. Topological Optimization of Convex and Concave Dies Structure

According to the analysis of the real working conditions of the convex and concave molds, based on the SIMP material interpolation optimization model of the convex and concave molds, OptiStruct is used to solve. OptiStruct describes the three elements of the optimal design through different types of information cards.(1)
Set optimization variablesIn the topology optimization of convex and concave molds, the tighter the optimized cell density is, the more important the material is. Moreover, this should be retained when redesigning the convex and concave mold structure [24]. As it tends to 0, it can be determined according to the analysis of the convex and concave structure. Therefore, the design area and nondesign area should be selected when designing optimization variables.(2)
Define optimization goals and constraintsThrough the numerical simulation of the molding process of the main body of the car coat rack, the node number of the stress concentration and deformation of the mold during the molding process are obtained. Moreover, the number and displacement values with time change are obtained by applying LS-Prepositi postprocessing. Export as a .txt file, and define the corresponding node displacement load in Hyperworks according to these results, in order to optimize the mold design. Figure9 shows the nodal displacement results derived from the concave and convex molds.Figures10 and 11 show that boundary conditions are applied and the optimization area is determined. The analysis of the mold structure shows that, in the optimized design, the concave and convex model surfaces and the location of the cooling water channel will remain unchanged. Moreover, the design area and the nondesign area are divided according to this, and the pink area of the concave mold is determined as the nondesign area. Furthermore, the blue part is the design area, the light color of the punch is the nondesign area, and the dark color is the design area. The boundary conditions are as described above.(3)
Trial calculationThe preliminary trial calculation is passed, as shown in Figure12, through the topology optimization algorithm and through the iterative calculation of this process to find the best path for load transfer, to obtain the optimal solution for topology optimization, and thus to obtain the best distribution of convex and concave mold materials.(4)
Topology optimization resultsAfter the lamination calculation, the topological results of the die and punch of the die are obtained. Figure13 illustrates a topographic cloud image with a density threshold of 0.2.Figure 9
Derived node displacement result graph. (a) Nodal displacement results derived from the die. (b) Nozzle displacement results derived from punch.
(a)(b)Figure 10
Boundary condition diagram of the upper die definition. (a) The upper die defines the boundary conditions. (b) The upper die defines the boundary conditions.
(a)(b)Figure 11
Boundary condition diagram of the lower die definition. (a) The lower die defines the boundary conditions. (b) The lower die defines the boundary conditions.
(a)(b)Figure 12
Preliminary calculation results. (a) Calculation result of upper die displacement. (b) Calculation result of upper mold equivalent stress. (c) Calculation result of lower die displacement. (d) Calculation result of equivalent stress of lower die.
(a)(b)(c)(d)Figure 13
Topology optimization morphology. (a) Top mold topography. (b) Topology of the lower die.
(a)(b)


#### 4.2.3. Redesign of the Convex and Concave Mold Structure

According to the topological optimization of the optimal path of the material distribution, the topological shape of the model can be obtained. Through the OptiStruc-OSSmooth program, the IGES general CAD file is exported and fed back to CATIA. The CATIA software is used to reverse the convex and concave mold structure to obtain the topology. Figure14 presents the optimized convex-concave structure.Figure 14
Optimal design of the convex and concave die. (a) Punch optimization design results. (b) Die optimization design results.
(a)(b)By optimizing the topology of the convex and concave mold structure of the automobile coat rack compression mold, the optimized design results of the convex and concave mold structures of the automobile coat rack mold are obtained. After the CATIA volume measurement, the optimized design of the convex and concave mold can meet the strength and on the premise of stiffness, weight loss was 15.6%.


## 4.1. Lightweight Design Process of the Mold Structure

According to the numerical simulation method for the molding process of automotive interior parts and the topology optimization theory of the variable density method [21], a lightweight design method for the mold structure of interior parts molding is proposed.(1)
According to the mold design requirements, the structure of the compression molding convex and concave mold is designed in accordance with the traditional mold design guidelines.(2)
With the help of the LS-DYNA finite element program, the compression molding process of the main body of the coat rack is numerically simulated to obtain the force information of the convex and concave mold nodes in the actual working conditions.(3)
Load mapping is used to construct convex and concave boundary conditions of die structure optimization.(4)
The improvement space of the mold structure design is clarified, and the design area and nondesign area are determined.(5)
Based on the OptiStruct, the topological optimization of the convex-concave structure is iteratively solved.(6)
The final shape is optimized according to the topology of the convex and concave mold structure, and the three-dimensional software is utilized to redesign the mold structure.(7)
According to the redesigned mold structure for mold processing and manufacturing.


## 4.2. OptiStruct-Based Optimization Design of the Mold Topology

OptiStruct uses LS-DYNA for preprocessing and definition of structural optimization problems. After completing the finite element modeling in LS-DYNA, the optimization definition panel is used to define optimization variables, constraints and goals, and optimization parameters. Then, OptiStruct is submitted for structural analysis and optimization. Figure8 shows the OptiStruct structure optimization design flowchart [22].Figure 8
Topology optimization flowchart.


### 4.2.1. The Establishment of the Topological Optimization Model of the Mold Convex and Concave Mold Structure

In the optimization process of the convex and concave dies, the problem of the mass minimization can be equated to the problem of the volume minimization. The real working conditions of the convex and concave molds are analyzed. Moreover, based on the solid isotropic material penalty (SIMP) method and according to the numerical simulation calculation results of the compression molding, the maximum displacement of the nodes in the stress concentration area is used as the constraint condition. Furthermore, the material element density is used as the design variable. The minimum volume of the convex-concave structure is the SIMP optimization model [23], and the density penalty factor P is set to 3. The model is formulas (45) and (46),(45)Ee=ρep⋅E0,0<ρmin≤ρe≤1,(46)min:Vsubjectto:KU=F,0<ρmin≤ρe≤1,Uj≤Uj∗,where ρi are V denote the design variable for the topology optimization with the variable density method and the overall structure volume, respectively. Moreover, U, K, and F are the overall displacement column vector, the overall structure stiffness matrix, and the external load vector, respectively.


### 4.2.2. Topological Optimization of Convex and Concave Dies Structure

According to the analysis of the real working conditions of the convex and concave molds, based on the SIMP material interpolation optimization model of the convex and concave molds, OptiStruct is used to solve. OptiStruct describes the three elements of the optimal design through different types of information cards.(1)
Set optimization variablesIn the topology optimization of convex and concave molds, the tighter the optimized cell density is, the more important the material is. Moreover, this should be retained when redesigning the convex and concave mold structure [24]. As it tends to 0, it can be determined according to the analysis of the convex and concave structure. Therefore, the design area and nondesign area should be selected when designing optimization variables.(2)
Define optimization goals and constraintsThrough the numerical simulation of the molding process of the main body of the car coat rack, the node number of the stress concentration and deformation of the mold during the molding process are obtained. Moreover, the number and displacement values with time change are obtained by applying LS-Prepositi postprocessing. Export as a .txt file, and define the corresponding node displacement load in Hyperworks according to these results, in order to optimize the mold design. Figure9 shows the nodal displacement results derived from the concave and convex molds.Figures10 and 11 show that boundary conditions are applied and the optimization area is determined. The analysis of the mold structure shows that, in the optimized design, the concave and convex model surfaces and the location of the cooling water channel will remain unchanged. Moreover, the design area and the nondesign area are divided according to this, and the pink area of the concave mold is determined as the nondesign area. Furthermore, the blue part is the design area, the light color of the punch is the nondesign area, and the dark color is the design area. The boundary conditions are as described above.(3)
Trial calculationThe preliminary trial calculation is passed, as shown in Figure12, through the topology optimization algorithm and through the iterative calculation of this process to find the best path for load transfer, to obtain the optimal solution for topology optimization, and thus to obtain the best distribution of convex and concave mold materials.(4)
Topology optimization resultsAfter the lamination calculation, the topological results of the die and punch of the die are obtained. Figure13 illustrates a topographic cloud image with a density threshold of 0.2.Figure 9
Derived node displacement result graph. (a) Nodal displacement results derived from the die. (b) Nozzle displacement results derived from punch.
(a)(b)Figure 10
Boundary condition diagram of the upper die definition. (a) The upper die defines the boundary conditions. (b) The upper die defines the boundary conditions.
(a)(b)Figure 11
Boundary condition diagram of the lower die definition. (a) The lower die defines the boundary conditions. (b) The lower die defines the boundary conditions.
(a)(b)Figure 12
Preliminary calculation results. (a) Calculation result of upper die displacement. (b) Calculation result of upper mold equivalent stress. (c) Calculation result of lower die displacement. (d) Calculation result of equivalent stress of lower die.
(a)(b)(c)(d)Figure 13
Topology optimization morphology. (a) Top mold topography. (b) Topology of the lower die.
(a)(b)


### 4.2.3. Redesign of the Convex and Concave Mold Structure

According to the topological optimization of the optimal path of the material distribution, the topological shape of the model can be obtained. Through the OptiStruc-OSSmooth program, the IGES general CAD file is exported and fed back to CATIA. The CATIA software is used to reverse the convex and concave mold structure to obtain the topology. Figure14 presents the optimized convex-concave structure.Figure 14
Optimal design of the convex and concave die. (a) Punch optimization design results. (b) Die optimization design results.
(a)(b)By optimizing the topology of the convex and concave mold structure of the automobile coat rack compression mold, the optimized design results of the convex and concave mold structures of the automobile coat rack mold are obtained. After the CATIA volume measurement, the optimized design of the convex and concave mold can meet the strength and on the premise of stiffness, weight loss was 15.6%.


## 4.2.1. The Establishment of the Topological Optimization Model of the Mold Convex and Concave Mold Structure

In the optimization process of the convex and concave dies, the problem of the mass minimization can be equated to the problem of the volume minimization. The real working conditions of the convex and concave molds are analyzed. Moreover, based on the solid isotropic material penalty (SIMP) method and according to the numerical simulation calculation results of the compression molding, the maximum displacement of the nodes in the stress concentration area is used as the constraint condition. Furthermore, the material element density is used as the design variable. The minimum volume of the convex-concave structure is the SIMP optimization model [23], and the density penalty factor P is set to 3. The model is formulas (45) and (46),(45)Ee=ρep⋅E0,0<ρmin≤ρe≤1,(46)min:Vsubjectto:KU=F,0<ρmin≤ρe≤1,Uj≤Uj∗,where ρi are V denote the design variable for the topology optimization with the variable density method and the overall structure volume, respectively. Moreover, U, K, and F are the overall displacement column vector, the overall structure stiffness matrix, and the external load vector, respectively.


## 4.2.2. Topological Optimization of Convex and Concave Dies Structure

According to the analysis of the real working conditions of the convex and concave molds, based on the SIMP material interpolation optimization model of the convex and concave molds, OptiStruct is used to solve. OptiStruct describes the three elements of the optimal design through different types of information cards.(1)
Set optimization variablesIn the topology optimization of convex and concave molds, the tighter the optimized cell density is, the more important the material is. Moreover, this should be retained when redesigning the convex and concave mold structure [24]. As it tends to 0, it can be determined according to the analysis of the convex and concave structure. Therefore, the design area and nondesign area should be selected when designing optimization variables.(2)
Define optimization goals and constraintsThrough the numerical simulation of the molding process of the main body of the car coat rack, the node number of the stress concentration and deformation of the mold during the molding process are obtained. Moreover, the number and displacement values with time change are obtained by applying LS-Prepositi postprocessing. Export as a .txt file, and define the corresponding node displacement load in Hyperworks according to these results, in order to optimize the mold design. Figure9 shows the nodal displacement results derived from the concave and convex molds.Figures10 and 11 show that boundary conditions are applied and the optimization area is determined. The analysis of the mold structure shows that, in the optimized design, the concave and convex model surfaces and the location of the cooling water channel will remain unchanged. Moreover, the design area and the nondesign area are divided according to this, and the pink area of the concave mold is determined as the nondesign area. Furthermore, the blue part is the design area, the light color of the punch is the nondesign area, and the dark color is the design area. The boundary conditions are as described above.(3)
Trial calculationThe preliminary trial calculation is passed, as shown in Figure12, through the topology optimization algorithm and through the iterative calculation of this process to find the best path for load transfer, to obtain the optimal solution for topology optimization, and thus to obtain the best distribution of convex and concave mold materials.(4)
Topology optimization resultsAfter the lamination calculation, the topological results of the die and punch of the die are obtained. Figure13 illustrates a topographic cloud image with a density threshold of 0.2.Figure 9
Derived node displacement result graph. (a) Nodal displacement results derived from the die. (b) Nozzle displacement results derived from punch.
(a)(b)Figure 10
Boundary condition diagram of the upper die definition. (a) The upper die defines the boundary conditions. (b) The upper die defines the boundary conditions.
(a)(b)Figure 11
Boundary condition diagram of the lower die definition. (a) The lower die defines the boundary conditions. (b) The lower die defines the boundary conditions.
(a)(b)Figure 12
Preliminary calculation results. (a) Calculation result of upper die displacement. (b) Calculation result of upper mold equivalent stress. (c) Calculation result of lower die displacement. (d) Calculation result of equivalent stress of lower die.
(a)(b)(c)(d)Figure 13
Topology optimization morphology. (a) Top mold topography. (b) Topology of the lower die.
(a)(b)


## 4.2.3. Redesign of the Convex and Concave Mold Structure

According to the topological optimization of the optimal path of the material distribution, the topological shape of the model can be obtained. Through the OptiStruc-OSSmooth program, the IGES general CAD file is exported and fed back to CATIA. The CATIA software is used to reverse the convex and concave mold structure to obtain the topology. Figure14 presents the optimized convex-concave structure.Figure 14
Optimal design of the convex and concave die. (a) Punch optimization design results. (b) Die optimization design results.
(a)(b)By optimizing the topology of the convex and concave mold structure of the automobile coat rack compression mold, the optimized design results of the convex and concave mold structures of the automobile coat rack mold are obtained. After the CATIA volume measurement, the optimized design of the convex and concave mold can meet the strength and on the premise of stiffness, weight loss was 15.6%.


## 5. Conclusion

The following conclusions are drawn from the present study:(1)
A mechanical model of the composite two-layer structure of the main body of the car coat rack is established. By solving the model, it is found that the main factors affecting the mechanical properties of the composite two-layer structure are the thickness of the PP glass fiber sheet and the content of glass fiber. The process parameters of the molding of the main body of the hat rack are prepared: the sheet preheating is performed by convection with hot air; the temperature of the hot air is 250°C; the preheating time is 110 s; the molding pressure is 13 MPa; and the holding time is 80 s. Using the experimental method, the average shrinkage rate of the material is calculated to be 0.354%, and the shrinkage rate of the solidification and cooling of the product is compensated by error compensation.(2)
Aiming at the current empirical design of mold designers, a lightweight design process for molding molds for interior parts is proposed. The LS-DYNA finite element program is used to simulate the molding process to obtain the actual working conditions. Under the stress concentration area of the mold, the nodal force information with large deformation is obtained. Therefore, the load mapping is used as the boundary condition of the next mold topology optimization.(3)
A mathematical model for SIMP material interpolation optimization of the male and female molds is established. Based on the variable density method, OptiStruct software is used to design the male and female molds of the main body of the automobile coat rack compression mold. The weight of the mold is reduced by 15.6%, while ensuring the structural rigidity and strength of the mold.



---

*Source: 6665753-2021-04-30.xml*