Research Article Implementation of Structural Assemblies by Simultaneous Projection and Density-Based Topology Optimization

This article proposed a methodology that combines two well-known projections and modiﬁed density-based optimized techniques in one formulation methodology. This methodology contains an eﬀective explicit geometric entity identiﬁed by shape variables that provide easy control in desired particular regions; implicit density-based topological optimization entities utilized topological variables that oﬀer critical design elsewhere. Our main attractive key point of this combined formulation approach is structural assemblies. Structure always manufactures in many patches and joins them by utilizing the structural assemblies such as welding and riveting. It is not easy to execute the structural patches at the required region without acknowledging their dimensions. This proposed approach demonstrates the competence to impose the restrictions related to shape and topological variables of interfaces among the speciﬁc patches. Numerous standard numerical examples make sure the validation of the introduced methodology. It remarked that the optimal design could minimize compliance and the minimum number of iterations through numerically performed, concerning computational cost minimized without any kind loss of accuracy of the ﬁnal structure.


Introduction
Topology optimization is well-known computational structural methodology which we used for spawning different fundamental structural design. is methodology has acquired significant evolution in various engineering fields and a wide range of industrial areas. During the previous three decades, the inaugural study of Bendsøe and Kikuchi [1] proposes various numerical techniques and mathematical expressions. Topology optimization is more valuable as parallel to the size and shape optimization. It enables optimized designs of complex geometry that we are incapable of solving and shape optimization. Topology optimization increases with every passing day in industrial engineering, automobile, and aerospace industries sides to achieving minimum weight and maximum stiffness purpose. Despite remarkable advances in computer accomplishments, large-scale topology optimization still faces difficulty due to massive computational costs. In topology optimization, many researchers have put their efforts into decreasing the computational expense regarding their point of view. Nguyen et al. [2] advised the multiresolution topology optimization (MOTP) approach to obtain the optimized design at a lower computational cost. Lieu et al. and Filipov et al. [3,4] mixed up this approach in the framework of isogeometric analysis and polygonal elements' techniques to obtaining optimized structure at a lower cost. Wu et al. [5] suggested a new method of exploring optimization and structural analysis for producing bone-like pours' designs. Long et al. [6] indicated that a well-organized and effective quadratic formulation depends on the standard-type isotropic material with penalization technique utilized to solve transient heat conduction problems. Generally, in the solid isotropic material with penalization (SIMP) technique, material properties of every element are described explicitly corresponding to relative density. e structured domain discretized into the number of FE; we used these elements to design variables for the optimization structure. Numerous methods are discussed in depth in topology optimization, including formulations and techniques. For example, the level set technique [7], the (SIMP) approach [8], structural and bidirectional structural optimization [9,10], and metaheuristic methods enabled topology optimization to obtain pleasing layouts. Zhang et al. [11] introduced a fair number of holes in the preferred design domain. ese holes improve the inhomogeneity of the structure domain which is helpful to obtaining the optimized microstructure results. Keshavarzzadeh et al. [12] suggested Gaussian process-based methods for solving the stress-based topology optimized. Amir and Shakour [13] described a density-based technique for optimizing the structure of prestressed concrete beams.
Wang et al. [14] suggested an approach for gaining maximum stiffness of continuum design within different moveable hole sizes utilizing the moving morphable bars. Long et al. [15] proposed a reciprocal SIMP technique for stress-constrained optimization continuum structure containing the harmonic load. Identifying the solution of novel systems shows that the significant performance of topology optimization is more efficient than traditional techniques. In the numerous manufacturing technologies, some geometrical characteristics of the optimized structure are suffering from few drawbacks, just like the deficiency of direct geometrical control (indeed, some geometric characteristics of the optimized structure can be controlled implicitly). In comparison, implicit level set and density-based constraints proposed appropriate design freedom. According to the literature review, it is difficult to obtain direct control throughout the resulting geometry by utilizing the simple LSM or density-based approach. Implicitly, topology optimization achieves by using different techniques utilizing existing ways. In the following form, the associated problem is summarized by using the implicit method. First, to provide the structural feature sizes, it is difficult to precisely control the structure although it is imperative for manufacturing viewpoint, under the framework of implied topology optimization [16][17][18][19][20][21][22][23][24].
is is due to the unavailability of explicit geometry information implied in model optimization and embedded techniques to gain full-length scale control. Furthermore, unique methods must be suggested for convenience and scale control [25,26]. Additionally, it is challenging to link computer-aided design (CAD) and optimization models for modeling systems under the implied system or framework. In these two settings, structural geometries are presented differently [27]. Second, three-dimensional problems in the suggested topology optimization approach are more significant due to design variables relatively. More design variables are achieved through higher resolutions.
ird, the optimization model and analysis model, in implicit approaches, are always coupled strongly. e problems that occur are checkerboard pattern and spurious local vibration. Sometimes, it can lead to simple numerical glitches that may stop the optimization algorithm from approaching purposeful explanations, notably the multiphysics setting, which is considered a topology optimization problem.
Recently, in the topology optimization process, an alternative so-called geometric projection class is dependent on the explicit distribution of parameters of the structure. Norato et al. [28] suggested a density-based geometric projection technique with differentiable mapping from discrete design variables. Meng and Keshtegar [29] proposed the adaptive conjugate single-loop approach to upgrade the efficiency and robustness of reliability-based design optimization RBDO utilizing the single-loop system based on the conjugate gradient. Besides, Zhang et al. [30] improved this technique for 3D structure execution. Furthermore, Watts and Tortorelli [31] developed a geometric projection approach to lattice material from inverse homogenization. White and Stowel suggested a practical approach to apply Fourier representation in combination with density-based topology optimization. Meng et al. [32] proposed an innovative hybrid RBTO approach to sort out multisource ambiguities, material variations, and assembling errors of geometric dimensions and loads. According to the manufacturing point of view, most of the above-discussed methodology indicates that these are an explicit geometric representation. So, in this way, we got easy access among CAD and optimization models for modeling system under the implied system or framework.
In this article, we combined the projection and density-based optimized formulation methodology. e following main key points are an effective explicit geometric entity identified by shape variables. ese variables provide easy control in desired selective regions. On the contrary, implicit density-based topological optimization entities used topological variables that play a critical role in other design parts.
ere are numerous benefits of combining projection, and density-based optimized formulation methodology can be summarized. First, coupling the shape and topological variables is helpful for prestressed concrete optimization (Amir and Shakour) [13].
is combined formulation approach's main attractive key point is structural assemblies, the structure produced in several patches. After that, it joined them by utilizing structural assemblies such as welding and riveting. It is challenging to execute the structural patches at the required region without acknowledging their dimension.
is proposed approach demonstrates the competence to impose the restrictions related to shape and topological variables of interfaces among the specific patches.
e order of this article is as follows. In Section 2, the primary setting of geometric and density based according to sigmoid function entities is explained. Furthermore, we briefly describe the combination methodology for specific control and belonging to interface components in Section 3. Section 4 indicates the sensitivity analysis, and Section 5 briefly discusses the capability of the suggested methodology concerning numerical aspects. Nevertheless, the results are stated in Section 6, respectively.

Topology Optimization Formulation.
To elaborate on the combination of explicit geometric projection and implicit density-based topology optimization methodology in a unique formulation, a practical, precise geometric entity is identified by shape variables to provide easy control in desired selective regions. On the other side, implicit densitybased topological optimization entities are utilized topological variables that offer critical design elsewhere. Here, geometric positions of nodal coordinates of line segments are utilized. e general mathematical formulation for structural topology optimization problem is [33] min ρ,x where f(ρ, x) is an objective function, x and ρ are significant design variable vectors concerning nodal coordinates and densities, where N ele and N node are the number of FE besides nodes classifying the projection profile, x lb denotes the lower bound of the geometric coordinate variable of the node, and x ub corresponds to its upper bound, respectively. K represents the stiffness matrix, although the symbols u and f are so-called displacement and applied loads, respectively.

Density-Based Projection Parameterization.
Obtain an improved depiction of the optimized topological structure by filtering the finite element mesh to achieve a better representation of the optimized topological structure by filtering the finite element mesh. Low-order element discretized the design ached through the distinguished checkerboard pathway. To prevent such problems, we assume density field variables such as ρ ⟶ ρ ⟶ ρ; the field of design variable ρ is regularized to get rid of these checkerboard pathway and mesh dependency through the density filter, which is denoted by ρ i [34]. e density filter ρ i components calculate average weight of the neighboring values as Here, the range of design variable ρ enclosed between one and zero helps differentiate between the solid and void spaces. On the other side, N e, i is the neighborhood set of cells lying in the specified filter domain of cell i. ω(Δx ij ) is a linear weighting function that is the distance between the centroids of the cells i and j. e mathematically linear weighting function is defined as where R represents the specified filter radius. Taking the derivative of ρ i w.r.t. the density variable field ρ j , we obtain as follows: Equation (2) can be rewritten in the matrix form as follows: where entries (H s ) and (H) are vectors that are [N ele × 1] and [N ele × N ele ]. ese vectors are described as In this paper, at this stage, to project a parameterized smooth projection function proposed in Guest et al. [18,[35][36][37], the filtered density regarding smooth Heaviside function that needs to have physical characterizations is mathematically calculated as In equation (7), the parameter ƞ values are lying between zero to one. By introducing the gamma (c) parameter, the physical density can be improved as Here, the behavior of equation (8) represents that the smooth Heaviside function is continuous. For ƞ � 0, the estimated assessment of the smooth Heaviside function ρ is related to a filtered variable that helps to certify the length scale on the solid region. In contrast, similarly on the opposite side, it helps to ensure the length scale on the solid region at ƞ � 1. Here, the concise and innovative smooth Heaviside function expression that rests on sigmoid function σ notwithstanding is traditionally illustrated as tanh function. is naive expression provides more straightforward and faster computing despite the defined smooth Heaviside function in the previous equation (8). is alternative smooth Heaviside sigmoid-based function expression is given as where a, b, C, and d are parameters. e partial derivative of the sigmoid-based smooth Heaviside projection regarding the design variable is shown respectively:

Journal of Mathematics
By utilizing the robust topology optimization technique (RTO) along eroded ρ e , intermediate ρ i , and dilated ρ d density vectors which explicitly depend on the represented filter density field ρ although ρ filter density also explicitly depends upon on design variable ρ, respectively, Lazarov et al. utilize the proposed dilated ρ d and eroded ρ e density vectors in the concerning problem for quantifying the worst case of volume and stiffness. As mentioned above, eroded ρ e and dilated ρ d density fields are written as follows: where η e � 0.6, η i � 0.5, and η d � 0.4 are threshold projection parameters of the eroded intermediate and dilated designs and c governs the sharpness of the projected smooth projection. e minimum value of the gamma parameter is 0.1.
In the optimization procedure, the intended dilated field vector ρ d is utilized for evaluating the volume constraint whereas and the eroded ρ e density field vector is utilized for the calculation of the compliance and structural stiffness, respectively. e static equilibrium equation is obtained as the displacement field. is field is calculated due to utilizing the linear equation system:

Material Interpolation Schemes.
Material interpolation approaches make sure the acknowledgment of intermediate material throughout the solution procedure. In topology optimization, usually material properties utilized for intermediate densities and material densities may obtain any suitable value between 1 and 0. e material penalization mathematically performs to carry out black and white design. According to Young's Modulus E of each element, the eroded density field vector is computed with the help of the modified solid isotropic material penalization (SIMP) interpolation technique [38]. e idea of the scheme is that where E min is a slight stiffness which is assigned to void spaces concerning the stiffness matrix not to be singular, and usually, PE is used as a solid isotropic material penalization factor which should be between 1 and 3, respectively. Here, taking the derivative of equation (13), In [39], the so-called RAMP (rational approximation of material properties) approach was recommended as an alternative interpolation approach. Let E(ρ e i ) denote the mathematical function used in the topology optimization is given as follows: where E min and PE is penalization factor which should be PE ≥ 0 and E max is an actual value of Young's Modulus for distribution material in topology optimization: We demonstrate the behaviors of the SIMP and RAMP in Figure 1, regarding to values of PE, respectively. For compliance minimization, we utilize a RAMP penalization instead of SIMP penalization. e output of equations (14) and (16) is obtained at given input zero and 1 regarding the design variable, where p > 1 and q > 0 are penalty parameters for SIMP and RAMP, respectively. Table 1 illustrates the computational results corresponding to the SIMP method, RAMP method, and Halpin Tsai, respectively.
It can be observed in Table 1 that the above amount of material is less at minimum computational period as compared to others.
ere are a number of conclusions according to above layouts which can be sketched as well as.
e SIMP method has lesser iteration numbers in comparison to the RAMP method.
e results/outputs obtained by the RAMP method would have a smaller volume as compared to the SIMP method if the objective function has a minimum value. e results/outputs acquired through the RAMP method are superior than the SIMP method if the objective function is minimum volume. Moreover, the results of the SIMP method would be more accurate than the RAMP method if the objective function is minimum compliance. e iterations obtained through the RAMP method are more efficient than the SIMP method.

Super-Gaussian Projection.
A super-Gaussian process is an illustrious machine learning tool for uncertainty quantification projects and is adopted in this article. It is assumed as piecewise linear profiles, and each linear profile contains four segments. Configuration of ith profile and every jth point which has the distance parameter are well-described; the super-Gaussian function is utilized for projection in topology optimization as follows: where μ Φ governs the curvature of the super-Gaussian projection. If two profiles are perpendicular or parallel concerning the desired position, each ith point encloses by two distances d i,1 and d i,2 corresponding to two segments which are belonging to a profile. e vicinity of the two profiles is correlated with the predicted variables Φ, whose values are more significant than one if the two profiles are similar to each other. is scenario is an undesired activity that we want to discourage. We want to have a predicted vector Φ at most equal to one for points close to two segments with two different profiles: for all i: e expression of equation (19) allows choosing merely minimum distance keen on according to the classification of maximum distance: for all i: We utilize the distance filter d 2 i to obtain an improved depiction of required layouts, concerning the suggested distance field d 2 i , is given as Equation (20) can be rewritten in the matrix form as follows: In this paper, at the final stage, to project a parameterized smooth projection, the super-Gaussian function is as follows:

Control over Projected Entities
is section expresses a critical, quite flexible topology optimization that can control many design features to mix projection and density-based design variables. It is assumed the volume constraint resting arranged the structure domain can be written as where v i is representing the volume of the ith finite element, whereas g * 0,dil is the existing solid volume fraction inside the dilated design.
Control of the volume of material inside a particular region is significant during the manufacturing development of structural assemblies. It requires minimizing the volume of material for welded at the different subcomponent interfaces of the structural assembly. Such desired control may be obtained by utilizing the local volume constraints; it is useful merely for those underlying structures in the proposed projection region. In particular, we express the constraints mathematically as where Φ i (x) is the projection function as a explicitly coordinate profile variable x for configuration concerning the ith element and g * 0,dil is the existing solid volume fraction inside the dilated design in topology optimization. e above-written equation represents a volume constraint introduced merely on a subdomain classified through elements for which Φ � 1. Now, we present the density filter ρ i components and calculate the average local material configuration of the neighboring N i values as Equation (27) can be rewritten in the matrix form as follows: e abovementioned equation N ele constraints can be rested in terms of an average percentage of solid material inside a neighboring of every element on maximum Von Mises stress as follows: e above maximum function is not appropriate for numerically optimized layout concerning first-order information because it cannot be differential. So, we adopt the Kreisselmeier-Steinhauser (KS) function of Von Mises or p norm function [12] as follows: where p is p norm which is not similar as mentioned PE penalization in the modified solid isotropic material with penalization approach. However, the projected quantities and optimization evolution suppose nearby discrete 0-1 values. After organizing these terms, we move next to the following constraints' computation: Here, to introduce another constraint to standardize a projection profile lying on a piecewise segmented profile, assume a limitation is resting on the segmented profile. Limit on a maximum slope is calculated as According to equation (30), the suggested constraint is mathematically calculated as is obtained constraint seems like linear constraint for the shape variable x i . is single constraint is helpful along with an optimization algorithm which is not adequate for the computation of linear control. erefore, adopt a p norm differentiable calculation as follows: is resultant constraint is more effective to standardize to the optimized segmented projection profile. e proposed projection variables Φ i allow to delegate distinct material characteristics to the elements rested inside the projection area. So, many can consign Young's modulus E Φ to the element inside the projection region with the help of the modified solid isotropic material with penalization interpolation or RAMP technique as follows: where proposing r E parameter specify distinct Young's modulus to solid elements within the selected projection region, whereas the Φ i projection variable must be equal to 1. e predefined expression of equation (35) allows the delegate solid elements to the inner side beside the outer side of the projection region through utilizing Young's modulus E Φ and E max . is scenario would be helpful during structural assembly computing of distinct welded interfaces. Amir and Lazarov [40] described a minimum length scale executed on black-and-white spaces is depicted along with every optimized layout by utilizing a circle whose diameter signifies the length scale. To control the dimensional irregularity of the length scale, use the density filter whose radius variation depends on the existing position of an element and the projection profile. e significance density filter representation can be as where this expression ω(Φx ij , Φ i ) defines the modified weight function: In equation (38), the filter radius is described as Here, r min and α are the significance filter radius and filter radius amplification feature, respectively.
Equation (36) can be rewritten in the matrix form as follows: Here, the filter mentioned above the modified weight function interchange to another Gaussian weight proposed by [36] e suggested Gaussian weight has the same layout, and the modified weight function resting n must be equivalent to 2, whereas the exp function and n both are positive, respectively. With due respect to this significant development, we assume the aforementioned (27) can be rewritten as follows: where H represents the Heaviside function which gives 1 at nonzero input values of x; on the other side, it provides only 0. is filter weight seems like the Gaussian weight function as the equation mentioned above. e weight ω G further converges to homogenous through rising the nonzero values of n and almost equivalent to 1. Follow the improvement of the filter radius r min resting upon the argument of exponent n such that r min � r min (1 +(1 − 2/n)), r min ϕ i , n � r min 1 + αϕ i , Obtain the required length-scale control resting upon rising the nonzero values of n due to r min . Pollini and Amir [33] introduced another constraint to govern the length scale along varying maximum feature size; the mathematical formulation of this constraint can be written down as Here, exponent n must be greater than to 2, for example {4, 6}.

Sensitivity Analysis
e objective function variables are implicitly stated over the equilibrium equations in optimization problems. Utilizing the adjoint method, sensitivity analysis assumed to formulate the gradient of objective and constraint functions which rested upon the topological variables and segmented projection profile x such as g 1 , g 2 , and g 4 although objective function f. Besides, the constraint g 3 merely rested upon the shape variable. e following gradient derivation of constraints g 1 , g 2 , g 3 , and g 4 is already described in [37]. Utilizing the adjoint sensitivity method, to formulate the gradient of the objective can be found as According to the equilibrium equation, if K(ρ e , x)u − f � 0, then f � f. After satisfying the equilibrium condition, take the derivative of f regarding the generic ς variable: Removing the implicit derivative (zu/zε) in the equation mentioned above. Calculate the vector λ to utilize the following equation: en, In our case, the objective function assumed is structural compliance. en, the problem behaves like a self-adjoint. Consequently, we put λ � − u and to identify the derivative (zf/zε), regarding ρ, utilizing the chain rule: Here, A is known as a matrix assembly operator, and the stiffness matrix also rests upon the shape coordinate of segmented profile and objective function f variable:

Numerical Discussion
Two-dimensional numerical examples demonstrate the proficiency and efficiency of the recommended modified mixed projection density-based optimization method. As assumed, moving limits in topology optimization for the geometric variable, which identifies the shape of the projection profile besides density variables, are equivalent to 0.005 and 0.2, respectively. In topology optimization, the MBB beam is a famous numerical assessment and point of the reference problem. Parameter c governs the sharpness of projected smooth projection of design variable ρ within the eroded ρ e , intermediate ρ i , and dilated ρ d outcomes. Furthermore, Young's modulus E max is equivalent to 1 for distribution material, E min is considered 1 × 10 − 6 a small stiffness is allotted to void spaces for the stiffness matrix not to be becoming singular, and the Poisson ratio is equivalent to v � 0.3, respectively. ere are many descriptive numerical cases to show the well-intentioned significance of the proposed modified method. Here, we present the initial structure domain of the MBB beam in Figure 2, whereas the participating material characteristic, the applied external loads, and geometry data, respectively. Following case (a), left side is clamped and the other opposite side force is applied at the bottom. In case (b), the left-hand edge is completely fastened, and a force f is vertically operated on the righthand of the middle side. Utilize the method of moving asymptotes as an optimizer to sort out the optimization problem [41]. Example 1. In the first example, the optimized cantilever structure problem under examination is presented in Figure 3. We endow with a moderately good evaluation of the suggested modified and merge projection density-based optimization algorithm to treat along with numerous constraints. For the finite element structure domain dimensions, nelx and nely are 3 × 1 where elements are placed horizontally and vertically, respectively. e primary purpose of the representation of black and blue dashed lines is to differentiate between the chosen outer boundaries and the actual structure of the design domain. PE � 3 is used as a solid isotropic material penalization constant factor and substituted for rE � 0.5 which is half of Young's modulus E max . Volume constraint is 40 percent, density filter is r, generally, r min � 10 is so-called density filter radius utilized to obtain minimized length scale, r max is so-called density filter radius used to obtain maximized length scale, in projection profile, r Φ � 10 is filter distance field, and g * 0 is the suitable volume fraction for g 0 constraint. We follow the focus obtained in an impressive result of the MBB beam. Moreover, in the projection profile, there are some 6 nodes that are fixed vertically at ycoordinates and moveable x-coordinates besides 5 segments.
After modifying this suggested approach, we get hold of optimized stable convergence as a reference result, respectively.

Example 2.
is example demonstrates the opportunity to optimize a structure with an interface between substructures also concurrently optimized. e existence of welding materials and various mechanical characteristics joining the different substructures are considered. e variable shape of the project file is used to represent the connection among the substructures. e finite element structure domain dimensions nelx and nely are 3 × 1. ere are two volume constraints to the total volume constraint g 0 and g 1 ; the constraint is a local volume that is operated at the projection region, respectively. To acquire supplementary efficient shapes' layout of optimized projection profile, we utilized g 3 as a third constraint. e projection profile organized with physical representation of a node in the following layout figure is a red circle, where the horizontally coordinate x belongs to a geometric variable. Geometric variables x i , where i � 1, . . . , 6, are firstly synchronized as 0.5, and their bounds which represent the figure are a magenta-colored dash line. Assume moving limits in topology optimization for the geometric variable, which identifies the shape of the projection profile besides density variables equivalent to 0.005 and 0.2.
It can be observed in Figure 4(a) that the above amount of material is less at minimum computational period as compared to the reference optimized result and at the exact value of volume fraction, respectively. We initially rested on the value μ Φ � 1.25 along increment of 1.4142 as well as in Figure 5; also, we take μ Φ � 1.50 along with the increment of 1.4142 at every sequel setup, which stops at a maximum value of 5. Here, μ Φ governs the curvature of the super-Gaussian utilized for the projection. Here, case (b) shows the final shape of the optimized projection area and case (c) is a graphical representation between the final compliance result and the total number of iterations along with obtained optimized history.
Example 3. Here, we represent 2D numerically example to make evident the proficiency and efficiency of the recommended algorithm. We optimize the MBB beam utilizing the modified solid isotropic material with penalization; for the finite element structure domain dimensions, nelx and nely are 300 × 100, where elements are placed horizontally and vertically, respectively. ere are two volume constraints to the total volume constraints g 0 and g 1 ; the constraint is a local volume that is operated at the projection region, respectively.  To acquire supplementary efficient shapes' layout of optimized projection profile, we utilized g 3 as a third constraint. e projection profile organized with physical representation of a node in the following layout figure is a red circle, where the horizontally coordinate x belongs to a geometric variable. Geometric variables x i , where i � 1, . . . , 6, are firstly synchronized as 0.5, and their bounds, which are represented in figure as a magenta-colored dash line. Parameters such as density filter radius is r min ; generally, r max is for maximum length scale control, r ɸ filter is utilized for transition corners in segment cut, and c Φ is helped to define the specified region, where g 1 constraint is functional. e threshold projection parameters of the eroded intermediate are dilated, where η e � 0.52 and η d � 0.48, and c governs the sharpness of the projected smooth projection. In the optimization procedure, the intended dilated field vector ρ d is utilized for evaluating the volume constraint, whereas the eroded ρ e density field vector utilizes the calculation of the compliance and structural stiffness. e above focused numerical topology optimized example unique layouts which are depicted in Figure 6 at minimum computational period. e fundamental goal is to find an optimized topology by simulation. is topology imposes volume limitation, which is revised locally in the projection area. is condition may interpret the requirement for minimizing the dimension of the structural elements at the interface of the various parts created distinctly and later bring together. With the help of the proposed technique, this impressive remarkable optimized layout utilizing the minimized amount of material in a subportion of the design domain is gained. is obtained optimal result is more flexible than the optimized reference result, although both have approximately the same volume fraction. Despite that, in this focus example, we utilized the total amount of material as smaller than the reference layout at the interface area.

Conclusion
is article proposed a novel methodology that is a combination of two well-known implicit projections and modified explicit density-based formulation techniques and computational procedure. e said procedure is used to drop down the dimensions of structural elements. is dimension reduction is at the interface of the various parts that are manufactured distinctly and later assembled. ese two representations coupled with explicit geometric entities that make availability of the specific control are provided, and the assigned domain and modified implicit entity make availability of structure freedom. However, geometric representations are linearly piecewise segments. ere are numerous benefits of the combination of projection, and density-based optimized formulation methodology can be summarized. First, due to coupling, the shape and topological variables are helpful for prestressed concrete optimization.
is existing work inspires by the optimal strategy of 2-dimensional structural assemblies. e applicability of the suggested technique illustrates numerous numerical example tests. Section 5 obtained layouts which magnificently indicate such selective limitations in the assigned domain and provide the complete freedom of modified density-based factors elsewhere.

Data Availability
Data and supplementary materials are available upon request and are not available in public.

Conflicts of Interest
e authors declare that they have no conflicts of interest to report with respect to the present study.