Next Article in Journal
The Bearing Capacity of Compressed Corrosion-Damaged Reinforced Concrete Elements under Lateral Pulse Loading
Next Article in Special Issue
A Memetic Algorithm for the Solution of the Resource Leveling Problem
Previous Article in Journal
The Geo-Distribution and Spatial Characteristics of Tulou Dwellings in Chaozhou, Guangdong, China
Previous Article in Special Issue
BIM-Based Automated Code Compliance Checking System in Malaysian Fire Safety Regulations: A User-Friendly Approach
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Weight Optimization of Discrete Truss Structures Using Quantum-Based HS Algorithm

1
School of Industrial Design & Architectural Engineering, Korea University of Technology & Education, 1600 Chungjeol-ro, Byeongcheon-myeon, Cheonan 31253, Republic of Korea
2
School of Liberal Arts, Korea University of Technology & Education, 1600 Chungjeol-ro, Byeongcheon-myeon, Cheonan 31253, Republic of Korea
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Buildings 2023, 13(9), 2132; https://doi.org/10.3390/buildings13092132
Submission received: 31 July 2023 / Revised: 18 August 2023 / Accepted: 21 August 2023 / Published: 22 August 2023

Abstract

:
Recently, a new field that combines metaheuristic algorithms and quantum computing has been created and is being applied to optimization problems in various fields. However, the application of quantum computing-based metaheuristic algorithms to the optimization of structural engineering is insufficient. Therefore, in this paper, we tried to optimize the weight of the truss structure using the QbHS (quantum-based harmony search) algorithm, which combines quantum computing and conventional HS (harmony search) algorithms. First, the convergence performance according to the parameter change of the QbHS algorithm was compared. The parameters selected for the comparison of convergence performance are QHMS, QHMCR, QPAR, ϵ , and θ r . The selected parameters were compared using six benchmark functions, and the range for deriving the optimal convergence performance was found. In addition, weight optimization was performed by applying it to a truss structure with a discrete cross-sectional area. The QbHS algorithm derived a lower weight than the QEA (quantum-inspired evolutionary algorithm) and confirmed that the convergence performance was better. A new algorithm that combines quantum computing and metaheuristic algorithms is required for application to various engineering problems, and this effort is essential for the expansion of future algorithm development.

1. Introduction

Quantum computers are rapidly emerging as a next-generation future technology and as one of the key technologies that will lead the fourth industrial revolution. Classical computers use the bit, expressed as 0 or 1, as the minimum information processing unit for computation. Quantum computers, on the other hand, use the qubit, or |1>, as the minimum information processing unit for computation. Due to these characteristics, the operation processing speed increases exponentially, attracting the interest of many researchers [1,2].
The possibility of a computational system based on quantum mechanics was first proposed by Feynman in 1982, and Deutsch proved in the same year that data processing was possible by applying quantum states [3]. In 1994, Shor’s algorithm and quantum searching algorithm were developed by Shor and Grover, and the full-scale development of quantum computing began [4,5,6]. Quantum computing is being developed based on the expectation that rapid computational processing is possible when quantum computer hardware is developed, and research is being steadily conducted in areas such as optimization, sensing, computing, and security [7,8].
In particular, quantum computing in the field of optimization is being combined with metaheuristic algorithms, and a new optimization algorithm based on quantum computing has been proposed [9,10]. Metaheuristics algorithms are applied and used in various engineering fields, such as weight optimization of structures, damage identification, optimal sensor placement, seismic collapse probability, life cycle cost, and smart dampers [11,12,13,14,15]. The quantum computing-based metaheuristic algorithm was first proposed by Narayanan and Moore in 1996. Narayanan and Moore applied quantum computing to genetic algorithms and tried to solve the traveling salesperson problem [16]. In 2000, Han and Kim expressed qubit probabilities and proposed the GQA (genetic quantum algorithm), which expresses the overlap of qubit states. GQA is expressed as a binary string by the probability of qubit, and the qubit rotates using the lookup table. In addition, a new termination condition was proposed using the convergence probability of qubit, and the possibility of GQA was confirmed by applying it to the knapsack problem [17]. In 2002, Han and Kim proposed the QEA (quantum-inspired evolutionary algorithm), which incorporated the evolutionary algorithm using the expression method of qubits used in GQA [18]. In 2004, Sun et al. proposed QDPSO (quantum delta-potential-well-based particle swarm optimization) using quantum wave functions and confirmed that it has a convergence performance similar to the results of conventional PSO (particle swarm optimization) algorithms using the benchmark function [19]. Since then, quantum computing has been applied to engineering problems and numerical problems in combination with algorithms such as the CSA (Cuckoo Search Algorithm), FA (Firefly Algorithm), GSA (Gravitational Search Algorithm), and TLBO (Teaching–Learning-Based Optimization) [20,21,22,23].
As explained earlier, new fields began to be created in the 1990s by combining quantum computing with various metaheuristic algorithms. The conventional HS algorithm was first proposed by Geem et al. [24] and is used to optimize many engineering problems because it is easy to apply to optimization problems [25]. The conventional HS algorithms, like other metaheuristic algorithms, underwent early attempts to combine them with quantum computing. In 2005, Geem proposed a BHS (binary HS) algorithm that expressed HM (harmony memory) in decimals in conventional HS algorithms [26], and in 2011, Wang et al. proposed a hybrid BHS algorithm using an ant system [27]. In 2013, Layeb proposed the QIHS (quantum-inspired HS) algorithm by combining quantum computing and conventional HS algorithms [28], and in 2016, Alfailakawi et al. tried to express the quantum gate as a two-dimensional circuit [29]. However, these attempts have the disadvantage of not being applicable to real-time problems because only binary problems determined by 0 or 1 can be solved, such as switch problems or knapsack problems. To solve this problem, in 2023, Lee et al. proposed a QbHS (quantum-based HS) algorithm that performs operations using probabilistic representations and overlapping qubit states and applied it to the weight optimization of truss structures with continuous cross-sectional areas [30]. However, it is not easy to determine the optimal parameter because the convergence performance according to changes in various parameters used in the QbHS algorithm is not comparable.
Attempts have been made to find the minimum weight by applying the conventional HS algorithm to the truss structure. In 2004, Lee et al. performed the weight optimization of 10-bar, 17-bar, 18-bar, 22-bar, 25-bar, 72-bar, and 200-bar truss structures, as well as 120-bar dome structures [31]. Lee et al. used a continuous cross-sectional area for weight optimization and performed size optimization. As constraints, the allowable stress of the elements, the maximum displacement of the nodes, and the buckling stress of the elements were considered. In 2005, Lee et al. performed the weight optimization of 25-bar, 52-bar, 72-bar, and 47-bar truss structures and used discrete cross-sectional areas [32]. As constraints, the allowable stress of the elements, the maximum displacement of the nodes, and the buckling stress of the elements were considered. In 2010, Srikanth et al. performed the weight optimization of a 22-bar truss structure and used the allowable stress of the elements, the maximum displacement of the nodes, and the buckling stress of the elements as constraints [33]. In 2012, Degetekin performed the weight optimization of 10-bar, 25-bar, 72-bar, and 200-bar truss structures using the IHS (improved HS) algorithm and used allowable stress, maximum displacement of the node, and buckling stress as constraints [34]. Since then, weight optimization of truss structures has been steadily performed using the HS algorithm [35,36,37]. Natural frequencies are widely used as constraints to avoid the resonance of structures in the weight optimization research of truss structures. However, there are few cases in which natural frequencies are included as constraints in studies that solve the weight optimization problem of truss structures using HS algorithms.
In order to apply the QbHS algorithm to various engineering problems, it is necessary to define the parameters with the best convergence performance, and it is necessary to apply them to various structure engineering problems using a quantum computing-based metaheuristic algorithm. Therefore, in this paper, we compare the convergence performance according to the parameter changes of the QbHS algorithm and perform weight optimization of the truss structure with discrete cross-sectional areas containing the natural frequency as a constraint. Examples adopted for the weight optimization problem are 20-bar, 24-bar, and 72-bar truss structures, each of which has a discrete cross-sectional area. Section 2 describes the QbHS algorithm, and Section 3 compares the convergence performance according to changes in parameters used in the QbHS algorithm. Section 4 performs weight optimization using example truss structures, and Section 5 concludes this paper.

2. Quantum-Based HS Algorithm

The QbHS algorithm first proposed by Lee et al. has a similar computational structure to the conventional HS algorithm and is classified into a total of five steps [30]. Although the conventional HS algorithm uses decimals, the QbHS algorithm is calculated using a binary, represented by the measurement of qubits. To express qubits, bracket notation is used and can be expressed as Equation (1). Here, α and β mean the probability amplitudes of |0> and |1>. α and β must satisfy Equation (2), and the single qubit state can be represented as a vector matrix, as shown in Equation (3).
| ψ > = α | 0 > + β | 1 >
| α | 2 + | β | 2 = 1
q = α β
In Step 1, an optimization problem is defined and parameters used in the QbHS algorithm are initialized. The parameters used in the QbHS algorithm are divided into parameters used in contentional HS algorithms such as QHMS (quantum harmony memory size), QHMCR (quantum harmony memory considering rate), and QPAR (quantum pitch adjusting rate), and parameters are added by combining them with quantum computing. Parameters added by combination with quantum computing include the number of qubits, ϵ , θ r , the number of measurements, tolBW, BWQ, q b w m a x , and q b w m i n .
In Step 2, QHM (quantum harmony memory) is initialized, and QHM is configured as shown in Figure 1. Here, N refers to the dimensions of the problem. For example, assuming that there are three qubits, each design variable is expressed as the probability information of qubits. The design variables of the qubit can be expressed as a binary through measurement. Since the information in the qubit consists of probability information, each measurement may have a different value.
Since the QbHS algorithm uses qubits, the process exists only when the initial qubit state is determined. Lee et al. used the QbHS H G algorithm when the qubit had the same probability of 0 or 1 being selected of 50%, and the QbHS R V algorithm when the probability of 0 or 1 being selected was random [30]. If the QbHS R V algorithm is used, it is evaluated once again with the H ϵ gate. The H ϵ gate serves to prevent the qubit from fully converging to 0 or 1 in the local minima state, and the convergence of the qubit is prevented by the size of ϵ . Equations (4)–(6) are used for the H ϵ gate. Equations (4) and (6) return the convergence probability to ϵ if | α | 2 or | β | 2 of the qubit converges above ϵ .
[ α i β i ] T = [ ϵ 1 ϵ ] T
[ α i β i ] T = [ 1 ϵ ϵ ] T
[ α i β i ] T = [ α i β i ] T
In Step 3, pitch adjusting is performed in the conventional HS algorithm; this is the most important step that determines the convergence performance of the algorithm. The QbHS algorithm is also performed by pitch adjusting the probabilities of QHMCR and QPAR. Lee et al. proposed performing sound control using the basic qubit state and expressed it as Equation (7) [30]. Here, r is a random number between 0 and 1, and pitch adjusting is performed around the current probability information. Q b w is calculated by Equation (8).
α i t + 1 = α i t 2 + r × Q b w r < 0.5 α i t + 1 = α i t 2 r × Q b w e l s e
Q b w = 0.7 × 0.9 × q b w m a x × exp l o g q b w m i n q b w m a x 0.7 × t t m a x
In addition, the QbHS algorithm was proposed to change the number of qubits that perform pitch adjusting according to the number of generations. Within a certain number of generations, all qubits perform pitch adjusting, but when the probability mean of qubits exceeds t o l B W , qubits, in addition to BWQ probabilities, are adopted to perform pitch adjusting. These characteristics improve the exploitation performance toward the end of the generation. The qubit performs rotation using the current generation and uses a rotation gate. A rotation gate is defined in Equation (9), and θ is defined in Equation (10). Here, Δ θ is determined using a lookup table, and θ r is used as a variable. Table 1 shows the lookup table.
α i t + 1 β i t + 1 = cos ( θ ) sin ( θ ) sin ( θ ) cos ( θ ) α i t β i t
θ = Δ θ × s i g n ( α i β i )
In Step 4, as in the conventional HS algorithm, whether or not to update is determined by comparing the qubit information of the previous generation with the qubit information in which pitch adjusting has been performed in Step 3. Since the QbHS algorithm constructs QHM using the qubit containing probability information, it compares fitness through the measurement of qubits. After comparison, the probability information of the qubit derived from better fitness is updated to the QHM. Through this process, the qubit converges to 0 or 1, and information accumulates.
In Step 5, the algorithm is terminated by the termination condition, and the optimization result is shown. As in the conventional HS algorithm, the QbHS algorithm mainly uses termination conditions using the number of generations. However, the QbHS algorithm can use the new termination condition proposed by Han et al. that uses the convergence of qubits [38], and the new termination condition can be expressed as Equations (11) and (12). Unlike the conventional HS algorithms, the new termination condition uses the convergence of qubits, so the algorithm can be terminated faster, and new results can be derived through measurement.
C b ( q ) = 1 m i = 1 m | 1 2 | α i | 2 | or C b ( q ) = 1 m i = 1 m | 1 2 | β i | 2 |
C a v = 1 N j = 1 N C b ( q j ) > ( 1 2 ϵ ) γ

3. Characteristics of the QbHS Algorithm

Table 2 shows the benchmark function used to compare the convergence performance according to parameter changes and the number of qubits used in each benchmark function [39,40]. Here, Min is the minimum value of the function, and t m a x is the maximum number of generations. The parameters selected for the comparison of convergence performance with changes in values are QHMS, QHMCR, QPAR, ϵ , and θ r . In addition, among the methods of initializing QHM for convergence performance comparison, the QbHS R V algorithm, which constructs the initial qubit state as a random number, was used.

3.1. QHMS

HMS (harmony memory size), one of the parameters for constructing HM (harmony memory) in conventional HS algorithms, is one of the parameters sensitive to convergence performance. Therefore, QHMS, with the same role as the HM of conventional HS algorithms, was adopted in QbHS algorithms, and the effect of the QHMS size change on convergence performance was compared. Table 3 presents the parameters for the interpretation of QHMS changes. QHMS was interpreted by changing it to 1, 5, 10, 20, 40, 60, and 100, and other parameters had fixed values. The interpretation was repeated 50 times.
Figure 2 is the convergence graph of best fitness according to the change in QHMS, and the interpretation results are summarized in Table A1. The smaller the size of the QHMS, the closer it was to green, and the larger the size of the QHMS, the closer it was to red. In all six functions, the larger the size of QHMS, the closer the result was to Min, and on the contrary, the smaller the size of QHMS, the farther the result was from Min. In terms of the average value using BF (best fitness) and MF (mean fitness), as presented in Table A1, QHMS was the worst at 6.92 when it had a value of 1 and the best at 1.17 when it had a value of 100.
Therefore, it was confirmed that the larger the QHMS, the better the convergence performance of the QbHS algorithm. This characteristic is because the larger the QHMS, the larger the QHM, so the probability of selecting a better value increases. It is also a characteristic similar to that of the conventional HS algorithms.

3.2. QHMCR and QPAR

The parameters that respond most sensitively to convergence performance in the conventional HS algorithms are known as HMCR and PAR [41]. The convergence performance of the conventional HS algorithms is the best when HMCR has values of 0.7–0.95 and PAR has values of 0.1–0.5. Therefore, the convergence performance of QHMCR and QPAR, which play the same role as HMCR and PAR in the conventional HS algorithms, was compared. Table 4 is a parameter for the interpretation of QHMCR and QPAR changes. QHMCR and QPAR were changed to 0.1, 0.3, 0.5, 0.7, and 0.9, respectively, and other parameters were fixed. The interpretation was repeated 50 times.
Figure 3 is a 3D bar graph expressing the analysis results of each benchmark function, and the analysis results are summarized in Table A2. In Figure 3, the closer to Min, the more blue, and the further away from Min, the more red. First, checking the change in QHMCR, the closer the value of QHMCR is to 0.9, the more blue it becomes. This characteristic means that it has a value close to Min, especially when it has a range of 0.7 to 0.9, where it has the closest value to Min. Second, the closer QPAR is to 0.1, the closer it is to Min. However, when QHMCR has a value close to 0.1, the effect on the change in QPAR does not occur clearly. That is, changes in QHMCR more dominantly affect convergence performance than do changes in QPAR.
Therefore, the closer QHMCR is to 1, the closer QPAR is to 0, and the better the convergence performance. These characteristics of the QbHS algorithm are similar to the range of HMCR and PAR, which are commonly used in the conventional HS algorithm.

3.3. ϵ

ϵ is a parameter used for the H ϵ gate, which prevents the qubit from fully converging to 0 or 1. That is, the smaller the value of ϵ , the closer the convergence of the qubit to 0 or 1, and the larger the value, the less the qubit converges. The values of ϵ changed to 0.00, 0.005, 0.01, 0.015, 0.02, and 0.03, and other parameters were used, as presented in Table 5. The interpretation was repeated 50 times.
Figure 4 shows the best or mean fitness according to the change in ϵ , and the analysis results are summarized in Table A3. The gray circles in Figure 4 are the results of a single analysis, and there are 50 gray circles depending on the size of ϵ . The red line means the best fitness, and the blue line means the mean fitness among 50 analyses. Checking the red lines, f 01 and f 03 were closest to Min when ϵ was 0.000, and f 02 and f 04 were closest to Min when ϵ was 0.005. However, if ϵ is 0.000, it is difficult to proactively escape from many problems with local minima. Therefore, 0.000 was excluded from the best ϵ range. In the average ranking using BF and MF, as shown in Table A3, it can be seen that the convergence performance is the best when ϵ is 0.005–0.015.
Figure 5 is the qubit probability according to the number of generations. Regardless of the size of ϵ , as the number of generations progresses, the qubit converges to one value. However, as ϵ increases, it converges at a value far from 1.0. Therefore, it was confirmed that an appropriate ϵ exists to increase convergence performance, and in this paper, the convergence performance was the best when the range was 0.005–0.015.

3.4. θ r

θ r is a parameter used in the rotation gate, and the rotation angle of the qubit is determined by the size of θ r . It can be predicted that if θ r has a large value, the qubit will converge quickly, and if θ r has a small value, the qubit will converge slowly. θ r was changed to 0.00, 0.05, 0.01, 0.02, 0.04, 0.06, 0.1, and 0.2, and other parameters were used as well, as shown in Table 6. The interpretation was repeated 50 times.
Figure 6 presents the change in the best or mean fitness according to the change in θ r , and the analysis results are summarized in Table A4. Checking for the best fitness, f 01 , f 02 , f 03 , and f 04 were closest to Min when θ r was 0.060, and f 05 and f 06 were closest when θ r was 0.040. In addition, when θ r was 0.000, the convergence performance was the worst for all functions because there was no rotation angle of the qubit. Checking the average values using BF and MF in Table A4, it can be seen that the convergence performance is the best when θ r has a range of 0.040–0.100.
Figure 7 is the qubit probability according to the number of generations. All functions converge to one value as the number of generations progresses, except when θ r is 0.000. In particular, the larger θ r is, the faster it converges to one value. However, when θ r is 0.2, the angle of rotation is so large that it converges near 0.9 and not any further.
Therefore, it is confirmed that θ r of an appropriate size exists to improve the convergence performance of the QbHS algorithm. In this paper, the convergence performance was the best when θ r was in the range 0.040–0.100.

4. Truss Structure Example

The weight optimization of truss structures is aimed at the minimum weight of the problem structure and can be defined as Equation (13). Equation (14) is a constraint for performing weight optimization [42].
Minimize F ( x ) = ρ i = 1 n B i A i L i + j = 1 m b j p e n a l t y
Subject   to g k ( x ) 0 , k = 1 , 2 , 3 , 4 , 5 , 6 , 7 g 1 ( x ) = | B i σ i | σ i m a x 0 g 2 ( x ) = | δ j | δ j m a x 0 g 3 ( x ) = | B i σ i c o m p | σ i c r 0 , σ i c r = k i A i E i L i 2 g 4 ( x ) = f r r r m a x 0 g 5 ( x ) = A m i n A i A m a x g 6 = C h e c k v a l i d i t y o f s t r u c t u r e g 7 = C h e c k k i n e m a t i c s t a b i l i t y
The cross-section that each truss structure element may have is discrete. Table A5 is the size of the cross-sectional area that each member can have and may have a total of 64 cross-sectional areas [43]. A total of 7 constraints were used for weight optimization of the truss structure. g 1 , g 2 , g 3 , and g 4 utilize the maximum stress of the member, the maximum displacement of the node, the buckling stress, and the natural frequency of the structure through FEA (finite element analysis) of the structure. A penalty of 10 4 is given if the constraint is not satisfied. g 5 is the range of cross-sectional areas that the element can have. g 6 evaluates the validity of the structure. That is, it determines whether there is a node serving as a support and a node acting as a load. A penalty of 10 9 is given if the constraint is not satisfied. g 7 evaluates the kinetic stability of the structure. To evaluate kinetic stability, check the degree of freedom and stiffness matrix of the structure. If the degree of freedom of the structure is greater than 0, a penalty of 10 8 is given. In addition, if eig(K) is less than 10 5 , a penalty of 10 7 is given.
The QEA algorithm was used to compare the weight optimization results of the QbHS algorithm. Table 7 presents parameters used in the QbHS and QEA algorithms to perform weight optimization, which was repeatedly interpreted 100 times.

4.1. 20-Bar Truss Structure

Figure 8 is the initial shape of a 20-bar truss structure, consisting of 9 nodes and 20 elements. E and ρ are 200,000 MPa and 7860 kg/m 3 , respectively, and the loading conditions are classified into two cases. The first case is F 1 = 500 kN and F 2 = 0 kN, and the second case is F 1 = 0 kN and F 2 = 500 kN. The allowable stress of each element is 180 MPa, and the maximum displacement that can occur on the Y-axis of the 4-node is 10 mm. Finally, the first natural frequency should be more than 60 Hz, and the second natural frequency should be more than 100 Hz.
Figure 9 is a convergence result graph for a 20-bar truss structure. For all three sets of algorithm results, the best and mean weights converge to one value. The best weights were derived as 332.503 kg for the QbHS H G algorithm, 331.211 kg for the QbHS R V algorithm, and 344.095 kg for the QEA algorithm. The mean weights were derived as 485.021 kg for the QbHS H G algorithm, 422.130 kg for the QbHS R V algorithm, and 478.228 kg for the QEA algorithm. The qubit probability shows that all algorithms converge to values close to 1.
Table 8 presents the weight optimization results of the 20-bar truss structure. Three algorithms contain elements 1, 5, 8, 11, 13, 15, 18, and 20, and the QbHS R V algorithm adds element 4. Therefore, the QbHS H G and QbHS R V algorithms adopted a total of eight elements, and the QEA algorithm adopted a total of nine elements. Among the results of the three algorithms, the QbHS R V algorithm had the best convergence performance by deriving the best, mean, and standard deviation (Std) of 331.211 kg, 422.130 kg, and 61.786, respectively.

4.2. 24-Bar Truss Structure

Figure 10 is the initial shape of the 24-bar truss structure, consisting of 8 nodes and 24 elements. E and ρ are 200,000 MPa and 7860 kg/m 3 , respectively, and the loading conditions are classified into two cases. The first case is F 1 = 100 kN and F 2 = 0 kN, and the second case is F 1 = 0 kN and F 2 = 100 kN. The allowable stress of each element is 180 MPa, and the maximum displacement that can occur on the Y-axis of the 5, 6-node is 10 mm. Finally, the first natural frequency should be more than 30 Hz.
Figure 11 is a convergence result graph of a 24-bar truss structure. In all three algorithm results, the best and mean weights converge to one value. The best weights were 243.762 kg for the QbHS H G algorithm, 250.718 kg for the QbHS R V algorithm, and 264.944 kg for the QEA algorithm. The mean weights were derived as 364.978 kg for the QbHS H G algorithm, 342.582 kg for the QbHS R V algorithm, and 364.060 kg for the QEA algorithm. The qubit probability shows that all algorithms converge to values close to 1.
Table 9 presents the weight optimization results of the 24-bar truss structure. As a result of weight optimization, all three algorithms contain elements 7, 8, 9, 12, and 15, and the QbHS H G algorithm adds elements 17, 20, 21, and 22. The QbHS R V algorithm adds elements 14, 16, 22, and 23, and the QEA algorithm adds elements 14, 16, and 24. Thus, the QbHS H G algorithm adopted a total of 10 elements, the QbHS R V algorithm a total of 9 elements, and the QEA algorithms a total of 8 elements. The QbHS R V algorithm derived 243.762 kg from the best weight, which was the best convergence performance, but the results of the QbHS R V algorithm, 342.582 kg and 72.289 for the mean and std, were the best. In addition, the results of the three algorithms satisfied all constraints.

4.3. 72-Bar Truss Structure

Figure 12 is the initial shape of a 72-bar truss structure, consisting of 20 nodes and 72 elements. The 72 elements were grouped into a total of 16 ( G 1 G 16 ) and are shown in Table A5. E and ρ are 68,950 MPa and 2767.99 kg/m 3 , respectively. The loading conditions are also classified into two cases. The first case has a load of 22.25 kN applied in the X, Y, and −Z directions at node 1. In the second case, 22.25 kN is applied in the −Z direction of nodes 1, 2, 3, and 4. The allowable stress of each element is 172.375 MPa, and the maximum displacement that can occur on the X or Y axis of nodes 1, 2, 3, and 4 is 6.35 mm. Finally, the first natural frequency should be more than 4 Hz, and the third natural frequency should be more than 6 Hz.
Figure 13 is a convergence result graph of a 72-bar truss structure. The results of the three algorithms show that the best and mean weight converge to one value. The weight optimization of the 72-bar truss structure resulted in the best weights of 445.833 kg for the QbHS H G algorithm, 449.190 kg for the QbHS R V algorithm, and 446.018 kg for the QEA algorithm. Mean weights were derived as 484.945 kg for the QbHS H G algorithm, 498.136 kg for the QbHS R V algorithm, and 522.369 kg for the QEA algorithm. The qubit probability shows that all algorithms converge to values close to 1.
Table 10 presents the weight optimization results of the 72-bar truss structure. The results of all three algorithms include groups 1, 2, 5, 6, 9, 10, 13, and 14, and the QbHS H G algorithm adds groups 8 and 11. The QbHS R V algorithm adds groups 4 and 15, and the QEA algorithm adds groups 8 and 15. Therefore, all three algorithms adopted a total of 10 groups. The best, mean, and std of the QbHS R V algorithm were 445.833 kg, 484.945 kg, and 21.306, respectively, showing the best convergence performance. In addition, the results of the three algorithms satisfied all constraints.

5. Conclusions

In this paper, the convergence performance of the QbHS algorithm, which combines quantum computing and conventional HS algorithms, was compared according to parameter changes. In addition, the weight optimization of 20-bar, 24-bar, and 72-bar truss structures with discrete cross-sectional areas was performed.
  • First, the convergence performance according to the size change of each parameter was compared. The convergence performance of the QbHS algorithm was better because the QHM increased as the QHMS increased. The larger the value of QHMCR, and the smaller the value of QPAR, the better the convergence performance of the QbHS algorithm. This aspect is judged to be similar to the conventional HS algorithm. The convergence performance of the QbHS algorithm was the best when ϵ had a range of 0.005–0.015 and θ r had a range of 0.040–0.100.
  • The weight optimization of 20-bar, 24-bar, and 72-bar truss structures with discrete cross-sectional areas was performed using the QbHS algorithm. For the 20-bar truss structure, the QbHS R V algorithm derived it as 331.211 kg, and for the 24-bar truss structure, the QbHS H G algorithm derived it as 243.762 kg. The 72-bar truss structure was derived as 549.954 kg by the QbHS H G algorithm.
Therefore, the convergence performance according to the changes in the parameters of the QbHS algorithm was compared using the six benchmark functions, and a parameter that could derive the best convergence performance was proposed. In addition, by applying it to the weight optimization of truss structures with discrete cross-sectional areas, a lower weight was derived than the QE algorithm, confirming that the convergence performance was better.
Research that combines quantum computing with existing metaheuristic algorithms, such as the QbHS algorithm, is creating new fields. However, it is extremely rare to apply quantum computing-based metaheuristics algorithms to engineering problems. Thus, quantum computing-based metaheuristic algorithms require a considerable amount of effort to solve optimization problems in engineering fields such as structures, machinery, and mechatronics. In addition, the truss structure applied in this paper is a basic structure, but it needs to be applied to optimize various structures such as large trusses and dome structures or complex buildings. Finally, since quantum computing-based metaheuristic algorithms are still in an early stage, it is necessary to combine them with various metaheuristic algorithms and develop various quantum operators. The application of various engineering problems of quantum computing-based metaheuristic algorithms and the development of quantum operators are expected to expand the field of computer information.

Author Contributions

Conceptualization, S.S., S.L. and J.H.; methodology, S.S. and J.H.; programming, S.S. and D.L.; validation, D.L. and S.S.; formal analysis, D.L.; investigation, D.L.; data curation, S.S. and D.L.; writing original draft preparation, D.L. and S.L.; visualization, D.L.; supervision, S.S. and D.L.; project administration, S.L. and D.L.; funding acquisition, D.L. All authors have read and agreed to the published version of the manuscript.

Funding

This paper was supported by the Education and Research Promotion Program of KOREATECH in 2023. This research was also supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF), funded by the Ministry of Education (RS-2023-00244008). Finally, this research was also supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT (NRF-2019R1A2C2010693).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Benchmark Function Results with Parameter Size Changes

Appendix A contains the benchmark function results according to the parameter size change of the QbHS algorithm. The parameters selected for convergence performance comparison are QHMS, QHMCR, QPAR, ϵ , and θ r , and six benchmark functions were used. The BF (best fitness), MF (mean fitness), and Std (standard deviation) of the interpretation results are indicated, in addition to the average ranking of BF and MF and the total average ranking of BF and MF.
Table A1. Benchmark function results according to QHMS.
Table A1. Benchmark function results according to QHMS.
Funct.IndexQHMS
1510204060100
f 01 BF3.162  × 10 + 2 1.387  × 10 0 1.200  × 10 1 4.785  × 10 3 6.863  × 10 4 1.403  × 10 4 3.085  × 10 5
MF1.048  × 10 + 3 8.719  × 10 0 8.578  × 10 1 6.648  × 10 2 4.966  × 10 3 1.149  × 10 3 1.558  × 10 4
Std4.461  × 10 + 2 7.235  × 10 0 8.215  × 10 1 4.487  × 10 2 3.712  × 10 2 7.524  × 10 4 1.158  × 10 4
f 02 BF3.091  × 10 0 2.955  × 10 1 2.507  × 10 2 4.300  × 10 3 9.549  × 10 4 4.891  × 10 4 6.105  × 10 5
MF5.528  × 10 0 2.104  × 10 0 1.706  × 10 0 1.603  × 10 0 8.970  × 10 1 5.772  × 10 1 1.394  × 10 1
Std1.344  × 10 0 4.933  × 10 1 9.041  × 10 1 8.460  × 10 1 9.534  × 10 1 8.727  × 10 1 4.135  × 10 1
f 03 BF1.076  × 10 0 2.910  × 10 2 1.429  × 10 3 3.956  × 10 5 2.842  × 10 6 3.060  × 10 7 1.477  × 10 7
MF1.849  × 10 0 2.260  × 10 1 1.637  × 10 1 1.112  × 10 1 7.753  × 10 2 6.532  × 10 2 5.948  × 10 2
Std9.218  × 10 1 1.322  × 10 1 1.352  × 10 1 9.429  × 10 2 9.152  × 10 2 5.555  × 10 2 5.274  × 10 2
f 04 BF2.238  × 10 + 1 1.460  × 10 + 1 8.689  × 10 0 4.945  × 10 0 3.709  × 10 0 3.236  × 10 0 1.236  × 10 0
MF4.297  × 10 + 1 3.006  × 10 + 1 2.310  × 10 + 1 1.808  × 10 + 1 1.499  × 10 + 1 1.164  × 10 + 1 9.036  × 10 0
Std9.897  × 10 0 9.566  × 10 0 6.573  × 10 0 6.608  × 10 0 6.888  × 10 0 5.747  × 10 0 4.374  × 10 0
f 05 BF4.890  × 10 + 2 4.588  × 10 + 2 1.719  × 10 + 2 3.477  × 10 + 1 7.286  × 10 1 6.248  × 10 1 4.186  × 10 1
MF1.426  × 10 + 3 9.323  × 10 + 2 5.291  × 10 + 2 2.149  × 10 + 2 4.734  × 10 + 1 8.553  × 10 0 1.102  × 10 0
Std3.708  × 10 + 2 2.738  × 10 + 2 2.615  × 10 + 2 1.399  × 10 + 2 7.457  × 10 + 1 2.469  × 10 + 1 2.629  × 10 1
f 06 BF1.772  × 10 + 1 1.713  × 10 + 1 1.718  × 10 + 1 1.638  × 10 + 1 1.805  × 10 + 1 1.441  × 10 + 1 1.617  × 10 + 1
MF1.944  × 10 + 3 3.597  × 10 + 2 2.778  × 10 + 2 2.087  × 10 + 2 1.746  × 10 + 2 1.079  × 10 + 2 1.080  × 10 + 2
Std1.051  × 10 + 4 5.556  × 10 + 2 4.280  × 10 + 2 4.588  × 10 + 2 3.433  × 10 + 2 1.376  × 10 + 2 2.186  × 10 + 2
RankingBF6.835.675.003.833.671.831.17
MF7.006.005.004.003.001.831.17
AVG(BF, MF)6.925.835.003.923.331.831.17
Table A2. Benchmark function results according to QHMCR and QPAR.
Table A2. Benchmark function results according to QHMCR and QPAR.
Funct.QPARIndexQHMCR
0.10.30.50.70.9
f 01 0.1BF1.279  × 10 + 4 7.571  × 10 + 3 4.188  × 10 + 3 1.260  × 10 + 3 5.714  × 10 0
MF1.864  × 10 + 4 1.354  × 10 + 4 8.048  × 10 + 3 2.544  × 10 + 3 2.023  × 10 + 1
Std2.287  × 10 + 3 2.652  × 10 + 3 1.555  × 10 + 3 6.756  × 10 + 2 1.119  × 10 + 1
0.3BF1.245  × 10 + 4 8.951  × 10 + 3 3.540  × 10 + 3 1.561  × 10 + 3 1.495  × 10 + 1
MF1.853  × 10 + 4 1.376  × 10 + 4 8.812  × 10 + 3 2.765  × 10 + 3 4.104  × 10 + 1
Std2.044  × 10 + 3 1.959  × 10 + 3 1.616  × 10 + 3 6.516  × 10 + 2 1.846  × 10 + 1
0.5BF9.470  × 10 + 3 9.985  × 10 + 3 6.326  × 10 + 3 2.109  × 10 + 3 3.052  × 10 + 1
MF1.871  × 10 + 4 1.425  × 10 + 4 8.938  × 10 + 3 3.252  × 10 + 3 8.475  × 10 + 1
Std2.794  × 10 + 3 1.942  × 10 + 3 1.304  × 10 + 3 7.582  × 10 + 2 3.833  × 10 + 1
0.7BF1.295  × 10 + 4 9.541  × 10 + 3 3.935  × 10 + 3 1.832  × 10 + 3 8.140  × 10 + 1
MF1.905  × 10 + 4 1.423  × 10 + 4 8.885  × 10 + 3 3.340  × 10 + 3 1.614  × 10 + 2
Std2.719  × 10 + 3 1.990  × 10 + 3 1.583  × 10 + 3 8.188  × 10 + 2 4.966  × 10 + 1
0.9BF1.360  × 10 + 4 9.186  × 10 + 3 5.664  × 10 + 3 1.823  × 10 + 3 1.035  × 10 + 2
MF1.878  × 10 + 4 1.409  × 10 + 4 8.785  × 10 + 3 3.822  × 10 + 3 2.707  × 10 + 2
Std2.301  × 10 + 3 1.856  × 10 + 3 1.438  × 10 + 3 9.882  × 10 + 2 8.303  × 10 + 1
f 02 0.1BF1.739  × 10 + 1 1.628  × 10 + 1 1.464  × 10 + 1 9.621  × 10 0 1.030  × 10 0
MF1.875  × 10 + 1 1.782  × 10 + 1 1.604  × 10 + 1 1.163  × 10 + 1 2.542  × 10 0
Std3.967  × 10 + 1 4.695  × 10 1 6.112  × 10 1 8.795  × 10 1 4.854  × 10 1
0.3BF1.706  × 10 + 1 1.644  × 10 + 1 1.465  × 10 + 1 9.175  × 10 0 2.442  × 10 0
MF1.872  × 10 + 1 1.796  × 10 + 1 1.609  × 10 + 1 1.200  × 10 + 1 3.238  × 10 0
Std4.467  × 10 1 4.182  × 10 1 6.520  × 10 1 1.126  × 10 0 4.882  × 10 1
0.5BF1.802  × 10 + 1 1.477  × 10 + 1 1.448  × 10 + 1 1.107  × 10 + 1 2.950  × 10 0
MF1.874  × 10 + 1 1.781  × 10 + 1 1.611  × 10 + 1 1.245  × 10 + 1 3.908  × 10 0
Std3.539  × 10 1 6.198  × 10 1 7.660  × 10 1 8.004  × 10 1 5.072  × 10 1
0.7BF1.755  × 10 + 1 1.679  × 10 + 1 1.274  × 10 + 1 1.075  × 10 + 1 3.618  × 10 0
MF1.873  × 10 + 1 1.788  × 10 + 1 1.615  × 10 + 1 1.286  × 10 + 1 4.864  × 10 0
Std3.785  × 10 1 4.855  × 10 1 9.103  × 10 1 8.175  × 10 1 5.925  × 10 1
0.9BF1.760  × 10 + 1 1.698  × 10 + 1 1.510  × 10 + 1 1.035  × 10 + 1 4.313  × 10 0
MF1.867  × 10 + 1 1.794  × 10 + 1 1.650  × 10 + 1 1.314  × 10 + 1 5.716  × 10 0
Std3.507  × 10 1 4.265  × 10 1 6.957  × 10 1 9.032  × 10 1 7.792  × 10 1
f 03 0.1BF1.218  × 10 + 2 6.916  × 10 + 1 4.396  × 10 + 1 6.147  × 10 0 7.966  × 10 1
MF1.592  × 10 + 2 1.148  × 10 + 2 6.585  × 10 + 1 1.835  × 10 + 1 1.038  × 10 0
Std1.866  × 10 + 1 1.932  × 10 + 1 9.703  × 10 0 5.688  × 10 0 5.307  × 10 2
0.3BF6.040  × 10 + 1 7.250  × 10 + 1 3.162  × 10 + 1 9.727  × 10 0 1.021  × 10 0
MF1.586  × 10 + 2 1.168  × 10 + 2 6.771  × 10 + 1 1.983  × 10 + 1 1.082  × 10 0
Std2.277  × 10 + 1 1.526  × 10 + 1 1.238  × 10 + 1 5.482  × 10 0 3.705  × 10 2
0.5BF1.227  × 10 + 2 6.784  × 10 + 1 3.018  × 10 + 1 1.224  × 10 + 1 1.048  × 10 0
MF1.649  × 10 + 2 1.159  × 10 + 2 7.145  × 10 + 1 2.244  × 10 + 1 1.209  × 10 0
Std1.672  × 10 + 1 1.769  × 10 + 1 1.523  × 10 + 1 5.062  × 10 0 1.118  × 10 1
0.7BF1.129  × 10 + 2 6.054  × 10 + 1 4.580  × 10 + 1 1.551  × 10 + 1 1.138  × 10 0
MF1.601  × 10 + 2 1.146  × 10 + 2 7.000  × 10 + 1 2.532  × 10 + 1 1.499  × 10 0
Std1.996  × 10 + 1 1.869  × 10 + 1 1.237  × 10 + 1 5.446  × 10 0 2.556  × 10 1
0.9BF1.042  × 10 + 2 6.623  × 10 + 1 4.747  × 10 + 1 1.716  × 10 + 1 1.266  × 10 0
MF1.626  × 10 + 2 1.216  × 10 + 2 7.440  × 10 + 1 2.787  × 10 + 1 1.831  × 10 0
Std2.052  × 10 + 1 1.931  × 10 + 1 1.325  × 10 + 1 6.395  × 10 0 3.024  × 10 1
f 04 0.1BF1.416  × 10 + 2 1.098  × 10 + 2 7.243  × 10 + 1 3.336  × 10 + 1 6.330  × 10 3
MF1.708  × 10 + 2 1.421  × 10 + 2 9.851  × 10 + 1 4.619  × 10 + 1 9.242  × 10 1
Std1.129  × 10 + 1 9.687  × 10 0 9.558  × 10 0 7.032  × 10 0 9.045  × 10 1
0.3BF1.440  × 10 + 2 1.190  × 10 + 2 7.421  × 10 + 1 3.444  × 10 + 1 1.637  × 10 0
MF1.753  × 10 + 2 1.436  × 10 + 2 1.034  × 10 + 2 5.306  × 10 + 1 5.301  × 10 0
Std9.039  × 10 0 1.069  × 10 + 1 9.931  × 10 0 7.812  × 10 0 1.853  × 10 0
0.5BF1.510  × 10 + 2 1.257  × 10 + 2 8.469  × 10 + 1 4.214  × 10 + 1 8.086  × 10 0
MF1.731  × 10 + 2 1.460  × 10 + 2 1.051  × 10 + 2 5.907  × 10 + 1 1.176  × 10 + 1
Std9.748  × 10 0 9.287  × 10 0 9.593  × 10 0 8.551  × 10 0 2.495  × 10 0
0.7BF1.451  × 10 + 2 1.256  × 10 + 2 8.393  × 10 + 1 5.166  × 10 + 1 9.087  × 10 0
MF1.748  × 10 + 2 1.461  × 10 + 2 1.123  × 10 + 2 6.500  × 10 0 2.147  × 10 + 1
Std1.244  × 10 + 1 1.017  × 10 + 1 1.153  × 10 + 1 9.367  × 10 0 4.274  × 10 0
0.9BF1.432  × 10 + 2 1.157  × 10 + 2 8.962  × 10 + 1 5.383  × 10 + 1 1.825  × 10 + 1
MF1.769  × 10 + 2 1.491  × 10 + 2 1.131  × 10 + 2 7.220  × 10 + 1 3.194  × 10 + 1
Std1.037  × 10 + 1 1.025  × 10 + 1 1.109  × 10 + 1 9.904  × 10 0 5.732  × 10 0
f 05 0.1BF3.511  × 10 + 3 2.715  × 10 + 3 1.553  × 10 + 3 4.516  × 10 + 2 2.155  × 10 1
MF4.245  × 10 + 3 3.400  × 10 + 3 2.196  × 10 + 3 7.358  × 10 + 2 6.138  × 10 1
Std2.498  × 10 + 2 2.479  × 10 + 2 2.562  × 10 + 2 1.514  × 10 + 2 2.722  × 10 1
0.3BF3.402  × 10 + 3 2.852  × 10 + 3 1.445  × 10 + 3 3.188  × 10 + 2 3.710  × 10 1
MF4.280  × 10 + 3 3.439  × 10 + 3 2.285  × 10 + 3 8.625  × 10 + 2 9.603  × 10 1
Std2.432  × 10 + 2 2.064  × 10 + 2 2.861  × 10 + 2 1.973  × 10 + 2 2.974  × 10 1
0.5BF3.396  × 10 + 3 2.796  × 10 + 3 1.652  × 10 + 3 4.463  × 10 + 2 7.093  × 10 1
MF4.312  × 10 + 3 3.439  × 10 + 3 2.376  × 10 + 3 9.559  × 10 + 2 1.708  × 10 0
Std2.443  × 10 + 2 1.949  × 10 + 2 2.522  × 10 + 2 2.075  × 10 + 2 9.339  × 10 1
0.7BF3.094  × 10 + 3 2.923  × 10 + 3 1.855  × 10 + 3 7.097  × 10 + 2 2.129  × 10 0
MF4.239  × 10 + 3 3.516  × 10 + 3 2.442  × 10 + 3 1.096  × 10 + 3 1.369  × 10 + 1
Std2.796  × 10 + 2 2.477  × 10 + 2 2.275  × 10 + 2 1.950  × 10 + 2 1.450  × 10 + 1
0.9BF3.374  × 10 + 3 3.012  × 10 + 3 1.615  × 10 + 3 6.972  × 10 + 2 9.921  × 10 0
MF4.277  × 10 + 3 3.539  × 10 + 3 2.490  × 10 + 3 1.248  × 10 + 3 6.615  × 10 + 1
Std2.328  × 10 + 2 2.200  × 10 + 2 2.551  × 10 + 2 2.397  × 10 + 2 4.256  × 10 + 1
f 06 0.1BF6.279  × 10 + 6 3.299  × 10 + 6 5.584  × 10 + 5 4.194  × 10 + 3 7.815  × 10 0
MF1.556  × 10 + 7 6.417  × 10 + 6 1.232  × 10 + 6 1.160  × 10 + 4 7.341  × 10 + 1
Std3.617  × 10 + 6 1.836  × 10 + 6 3.954  × 10 + 5 6.609  × 10 + 3 4.125  × 10 + 1
0.3BF7.724  × 10 + 6 3.085  × 10 + 6 4.562  × 10 + 5 2.498  × 10 + 3 1.222  × 10 + 1
MF1.647  × 10 + 7 6.917  × 10 + 6 1.353  × 10 + 6 2.014  × 10 + 4 9.499  × 10 + 1
Std4.446  × 10 + 6 1.905  × 10 + 6 3.870  × 10 + 5 1.307  × 10 + 4 8.885  × 10 + 1
0.5BF6.871  × 10 + 6 3.031  × 10 + 6 6.777  × 10 + 5 8.852  × 10 + 3 1.792  × 10 + 1
MF1.573  × 10 + 7 7.109  × 10 + 6 1.479  × 10 + 6 2.865  × 10 + 4 1.455  × 10 + 2
Std4.255  × 10 + 6 1.911  × 10 + 6 4.602  × 10 + 5 1.583  × 10 + 4 1.259  × 10 + 2
0.7BF6.470  × 10 + 6 2.280  × 10 + 6 8.091  × 10 + 5 1.692  × 10 + 4 2.201  × 10 + 1
MF1.706  × 10 + 7 6.661  × 10 + 6 1.712  × 10 + 6 4.240  × 10 + 4 1.425  × 10 + 2
Std4.091  × 10 + 6 1.997  × 10 + 6 6.701  × 10 + 5 2.043  × 10 + 4 1.204  × 10 + 2
0.9BF7.309  × 10 + 6 1.616  × 10 + 6 7.881  × 10 + 5 1.734  × 10 + 4 3.822  × 10 + 1
MF1.508  × 10 + 7 7.260  × 10 + 6 1.856  × 10 + 6 6.825  × 10 + 4 2.557  × 10 + 2
Std3.890  × 10 + 6 2.255  × 10 + 6 6.431  × 10 + 5 3.201  × 10 + 4 2.838  × 10 + 2
Table A3. Benchmark function results according to ϵ .
Table A3. Benchmark function results according to ϵ .
Funct.Index ϵ
0.0000.0050.0100.0150.0200.030
f 01 BF7.217  × 10 2 8.412  × 10 2 1.303  × 10 1 3.919  × 10 1 1.058  × 10 0 1.546  × 10 + 1
MF1.536  × 10 + 1 5.910  × 10 1 8.094  × 10 1 1.843  × 10 0 7.131  × 10 0 4.015  × 10 + 1
Std8.826  × 10 + 1 4.960  × 10 1 6.828  × 10 1 1.230  × 10 0 4.676  × 10 0 2.086  × 10 + 1
f 02 BF6.849  × 10 2 1.067  × 10 2 1.778  × 10 2 8.569  × 10 2 2.483  × 10 1 1.519  × 10 0
MF2.260  × 10 0 2.034  × 10 0 1.636  × 10 0 1.606  × 10 0 1.837  × 10 0 2.422  × 10 0
Std4.630  × 10 1 5.163  × 10 1 9.638  × 10 1 8.638  × 10 1 7.380  × 10 1 3.304  × 10 1
f 03 BF5.667  × 10 5 3.471  × 10 4 1.427  × 10 3 3.185  × 10 2 1.675  × 10 1 1.001  × 10 0
MF2.037  × 10 1 1.398  × 10 1 1.408  × 10 1 1.760  × 10 1 5.378  × 10 1 1.051  × 10 0
Std1.530  × 10 1 1.202  × 10 1 1.214  × 10 1 1.150  × 10 1 1.671  × 10 1 1.959  × 10 2
f 04 BF1.017  × 10 + 1 6.538  × 10 0 8.177  × 10 0 1.316  × 10 + 1 8.749  × 10 0 1.613  × 10 + 1
MF3.558  × 10 + 1 2.621  × 10 + 1 2.274  × 10 + 1 2.235  × 10 + 1 2.176  × 10 + 1 2.835  × 10 + 1
Std1.001  × 10 + 1 7.451  × 10 0 6.573  × 10 0 7.118  × 10 0 7.489  × 10 0 6.294  × 10 0
f 05 BF4.090  × 10 + 2 2.375  × 10 + 2 6.962  × 10 + 1 3.549  × 10 + 1 1.379  × 10 + 2 1.128  × 10 + 2
MF1.056  × 10 + 3 6.778  × 10 + 2 5.269  × 10 + 2 3.984  × 10 + 2 4.456  × 10 + 2 4.938  × 10 + 2
Std3.048  × 10 + 2 2.705  × 10 + 2 2.587  × 10 + 2 1.790  × 10 + 2 1.697  × 10 + 2 2.098  × 10 + 2
f 06 BF1.569  × 10 + 1 1.678  × 10 + 1 1.487  × 10 + 1 1.717  × 10 + 1 1.809  × 10 + 1 1.786  × 10 + 1
MF3.305  × 10 + 2 1.301  × 10 + 2 2.149  × 10 + 2 2.403  × 10 + 2 2.450  × 10 + 2 1.560  × 10 + 2
Std5.148  × 10 + 2 2.575  × 10 + 2 3.705  × 10 + 2 4.587  × 10 + 2 4.609  × 10 + 2 2.434  × 10 + 2
RankingBF2.832.332.173.674.675.33
MF5.332.672.672.333.334.67
AVG(BF, MF)4.082.502.423.004.005.00
Table A4. Benchmark function results according to θ r .
Table A4. Benchmark function results according to θ r .
Funct.Index θ r
0.0000.0050.0100.0200.0400.0600.1000.200
f 01 BF1.167  × 10 + 4 6.822  × 10 + 2 1.341  × 10 + 2 3.490  × 10 0 8.767  × 10 1 1.613  × 10 1 4.959  × 10 1 7.588  × 10 + 1
MF2.188  × 10 + 4 1.912  × 10 + 3 3.654  × 10 + 2 3.751  × 10 + 1 2.909  × 10 0 1.088  × 10 0 2.196  × 10 0 2.977  × 10 + 2
Std2.723  × 10 + 3 6.883  × 10 + 2 1.985  × 10 + 2 2.118  × 10 + 1 1.545  × 10 0 9.910  × 10 1 1.631  × 10 0 1.196  × 10 + 2
f 02 BF1.738  × 10 + 1 4.392  × 10 0 1.744  × 10 0 4.369  × 10 1 7.336  × 10 2 3.799  × 10 2 8.735  × 10 2 3.850  × 10 0
MF1.890  × 10 + 1 6.274  × 10 0 2.994  × 10 0 2.079  × 10 0 1.883  × 10 0 1.939  × 10 0 1.309  × 10 0 5.037  × 10 0
Std3.998  × 10 1 1.069  × 10 0 4.198  × 10 1 5.096  × 10 1 7.852  × 10 1 6.986  × 10 1 8.848  × 10 1 6.141  × 10 1
f 03 BF1.144  × 10 + 2 1.387  × 10 0 8.973  × 10 1 1.770  × 10 1 3.009  × 10 3 1.972  × 10 3 9.346  × 10 2 1.238  × 10 0
MF1.738  × 10 + 2 2.242  × 10 0 1.083  × 10 0 4.677  × 10 1 1.530  × 10 1 1.338  × 10 1 3.247  × 10 1 2.196  × 10 0
Std2.184  × 10 + 1 4.968  × 10 1 5.778  × 10 2 1.670  × 10 1 1.127  × 10 1 1.122  × 10 1 1.638  × 10 1 5.735  × 10 1
f 04 BF1.536  × 10 + 2 1.249  × 10 + 1 1.146  × 10 + 1 1.218  × 10 + 1 1.167  × 10 + 1 1.068  × 10 + 1 1.319  × 10 + 1 2.416  × 10 + 1
MF1.860  × 10 + 2 2.738  × 10 + 1 2.569  × 10 + 1 2.239  × 10 + 1 2.419  × 10 + 1 2.279  × 10 + 1 2.324  × 10 + 1 4.346  × 10 + 1
Std1.141  × 10 + 1 7.019  × 10 0 8.376  × 10 0 6.382  × 10 0 6.914  × 10 0 6.321  × 10 0 6.296  × 10 0 8.413  × 10 0
f 05 BF3.581  × 10 + 3 1.999  × 10 + 2 1.303  × 10 + 2 1.378  × 10 + 2 1.035  × 10 + 2 1.577  × 10 + 2 1.536  × 10 + 2 3.001  × 10 + 2
MF4.287  × 10 + 3 5.863  × 10 + 2 5.391  × 10 + 2 5.645  × 10 + 2 4.983  × 10 + 2 5.890  × 10 + 2 4.393  × 10 + 2 7.602  × 10 + 2
Std2.986  × 10 + 2 2.028  × 10 + 2 1.849  × 10 + 2 2.168  × 10 + 2 2.387  × 10 + 2 2.041  × 10 + 2 1.775  × 10 + 2 2.223  × 10 + 2
f 06 BF6.081  × 10 + 6 1.714  × 10 + 1 1.729  × 10 + 1 1.767  × 10 + 1 1.584  × 10 + 1 1.809  × 10 + 1 1.716  × 10 + 1 4.437  × 10 + 1
MF1.844  × 10 + 7 2.872  × 10 + 2 2.439  × 10 + 2 2.139  × 10 + 2 2.914  × 10 + 2 2.171  × 10 + 2 2.286  × 10 + 2 4.657  × 10 + 2
Std5.343  × 10 + 6 5.463  × 10 + 2 4.703  × 10 + 2 3.005  × 10 + 2 5.546  × 10 + 2 4.317  × 10 + 2 4.018  × 10 + 2 6.146  × 10 + 2
RankingBF8.004.673.674.002.673.503.835.67
MF8.006.174.673.003.172.502.176.33
AVG(BF, MF)8.005.424.173.502.923.003.006.00

Appendix B. Discrete Area of Truss Structure Element

Each element of the truss structure may have a discrete cross-sectional area and is adopted as one of a total of 64 cross-sectional areas.
Table A5. Discrete area of truss structure.
Table A5. Discrete area of truss structure.
No.Area (cm 2 )Thickness (cm)No.Area (cm 2 )Thickness (cm)
10.71610.15103324.77410.8880
20.90970.17023424.96770.8915
31.26450.20063525.03220.8926
41.61290.22663626.96770.9265
51.98060.25113727.22580.9309
62.52260.28343828.96770.9602
72.85160.30133929.61280.9709
83.63230.34004030.96770.9928
93.88390.35164132.06451.0103
104.94190.39664233.03221.0254
115.06450.40154337.03221.0857
126.41290.45184446.58061.2177
136.45160.45324551.41931.2793
147.92260.50224651.41931.2793
158.16770.50994759.99991.3820
169.40000.54704869.99991.4927
1710.08390.56664974.19431.5368
1810.45160.57685087.09661.6650
1911.61290.60805189.67721.6895
2012.83870.63935291.61271.7077
2113.74190.66145399.99981.7841
2215.35480.699154103.22561.8127
2316.90320.733555109.03201.8630
2416.96770.734956121.29011.9649
2518.58060.769157128.38682.0216
2618.90320.775758141.93522.1255
2719.93540.796659147.74162.1686
2820.19350.801760158.06422.2431
2921.80640.833161170.96742.3328
3022.38710.844262180.64482.3979
3122.90320.853863193.54802.4821
3223.41930.863464216.12862.6229

References

  1. Morsch, O. Quantum Bits and Quantum Secrets: How Quantum Physics Is Revolutionizing Codes and Computers; John Wiley & Sons: Hoboken, NJ, USA, 2008. [Google Scholar]
  2. Jones, J.A. Quantum computing with NMR. arXiv 2010, arXiv:1011.1382. [Google Scholar]
  3. Berthiaume, A.; Feynman, R.P. Quantum computation. Complex. Theory Retrosp. II 1997, 2, 23. [Google Scholar]
  4. Shor, P.W. Algorithms for quantum computation: Discrete logarithms and factoring. In Proceedings of the 35th Annual Symposium on Foundations of Computer Science, Santa Fe, NM, USA, 20–22 November 1994; pp. 124–134. [Google Scholar]
  5. Monz, T.; Nigg, D.; Martinez, E.A.; Brandl, M.F.; Schindler, P.; Rines, R.; Wang, S.X.; Chuang, I.L.; Blatt, R. Realization of a scalable Shor algorithm. Science 2016, 351, 1068–1070. [Google Scholar] [CrossRef]
  6. Bernstein, D.J. Grover vs. mceliece. In Proceedings of the Post-Quantum Cryptography: 3rd International Workshop, PQCrypto, Darmstadt, Germany, 25–28 May 2010; pp. 73–80. [Google Scholar]
  7. National Academies of Sciences, Engineering, and Medicine and others. Quantum Computing: Progress and Prospects; The National Academies Press: Washington, DC, USA, 2019. [Google Scholar]
  8. Zhang, J.; Tian, Q.; Tang, C.; Wang, L.; Xu, J.; Fang, J. Study on Worldwide Development and Trends of Quantum Technologies Based on Patent Data. Int. J. Inf. Educ. Technol. 2020, 10, 239–244. [Google Scholar] [CrossRef]
  9. Ross, O.H.M. A review of quantum-inspired metaheuristics: Going from classical computers to real quantum computers. IEEE Access 2019, 8, 814–838. [Google Scholar] [CrossRef]
  10. Gharehchopogh, F.S. Quantum-inspired metaheuristic algorithms: Comprehensive survey and classification. Artif. Intell. Rev. 2023, 56, 5479–5543. [Google Scholar] [CrossRef]
  11. Xiao, F.; Hulsey, J.L.; Chen, G.S.; Xiang, Y. Optimal static strain sensor placement for truss bridges. Int. J. Distrib. Sens. Netw. 2017, 13, 1550147717707929. [Google Scholar] [CrossRef]
  12. Xiao, F.; Zhu, W.; Meng, X.; Chen, G.S. Parameter identification of structures with different connections using static responses. Appl. Sci. 2022, 12, 5896. [Google Scholar] [CrossRef]
  13. Xiao, F.; Sun, H.; Mao, Y.; Chen, G.S. Damage identification of large-scale space truss structures based on stiffness separation method. In Structures; Elsevier: Amsterdam, The Netherlands, 2023; Volume 53, pp. 109–118. [Google Scholar]
  14. Rayegani, A.; Nouri, G. Seismic collapse probability and life cycle cost assessment of isolated structures subjected to pounding with smart hybrid isolation system using a modified fuzzy based controller. In Structures; Elsevier: Amsterdam, The Netherlands, 2022; Volume 44, pp. 30–41. [Google Scholar]
  15. Rayegani, A.; Nouri, G. Application of smart dampers for prevention of seismic pounding in isolated structures subjected to near-fault earthquakes. J. Earthq. Eng. 2022, 26, 4069–4084. [Google Scholar] [CrossRef]
  16. Narayanan, A.; Moore, M. Quantum-inspired genetic algorithms. In Proceedings of the IEEE International Conference on Evolutionary Computation, Nagoya, Japan, 20–22 May 1996; pp. 61–66. [Google Scholar]
  17. Han, K.H.; Kim, J.H. Genetic quantum algorithm and its application to combinatorial optimization problem. In Proceedings of the Congress on Evolutionary Computation, CEC00 (Cat. No. 00TH8512), La Jolla, CA, USA, 16–19 July 2000; Volume 2, pp. 1354–1360. [Google Scholar]
  18. Han, K.H.; Kim, J.H. Quantum-inspired evolutionary algorithm for a class of combinatorial optimization. IEEE Trans. Evol. Comput. 2002, 6, 580–593. [Google Scholar] [CrossRef]
  19. Sun, J.; Feng, B.; Xu, W. Particle swarm optimization with particles having quantum behavior. In Proceedings of the Congress on Evolutionary Computation (IEEE Cat. No. 04TH8753), Portland, OR, USA, 19–23 June 2004; Volume 1, pp. 325–331. [Google Scholar]
  20. Boussalia, S.R.; Chaoui, A. Optimizing QoS-based web services composition by using quantum inspired cuckoo search algorithm. In Mobile Web Information Systems: Proceedings of the 11th International Conference, MobiWIS 2014, Barcelona, Spain, 27–29 August 2014; Springer: Berlin/Heidelberg, Germany, 2014; pp. 41–55. [Google Scholar]
  21. Shareef, H.; Ibrahim, A.; Salman, N.; Mohamed, A.; Ai, W.L. Power quality and reliability enhancement in distribution systems via optimum network reconfiguration by using quantum firefly algorithm. Int. J. Electr. Power Energy Syst. 2014, 58, 160–169. [Google Scholar] [CrossRef]
  22. Soleimanpour-Moghadam, M.; Nezamabadi-Pour, H.; Farsangi, M.M. A quantum inspired gravitational search algorithm for numerical function optimization. Inf. Sci. 2014, 267, 83–100. [Google Scholar] [CrossRef]
  23. Gao, H.; Zhang, X.; Du, Y.; Diao, M. Quantum-inspired teaching-learning-based optimization for linear array pattern synthesis. In Communications, Signal Processing, and Systems: Proceedings of the International Conference on Communications, Signal Processing, and Systems; Springer: Berlin/Heidelberg, Germany, 2019; pp. 2106–2115. [Google Scholar]
  24. Geem, Z.W.; Kim, J.H.; Loganathan, G.V. A new heuristic optimization algorithm: Harmony search. Simulation 2001, 76, 60–68. [Google Scholar] [CrossRef]
  25. Dubey, M.; Kumar, V.; Kaur, M.; Dao, T.P. A systematic review on harmony search algorithm: Theory, literature, and applications. Math. Probl. Eng. 2021, 2021, 5594267. [Google Scholar] [CrossRef]
  26. Geem, Z.W. Harmony search in water pump switching problem. In Advances in Natural Computation: Proceedings of the 1st International Conference, Changsha, China, 27–29 August 2005; Proceedings, Part III 1; Springer: Berlin/Heidelberg, Germany, 2005; pp. 751–760. [Google Scholar]
  27. Wang, L.; Zhou, P.; Fang, J.; Niu, Q. A hybrid binary harmony search algorithm inspired by ant system. In Proceedings of the IEEE 5th International Conference on Cybernetics and Intelligent Systems (CIS), Qingdao, China, 17–19 September 2011; pp. 153–158. [Google Scholar]
  28. Layeb, A. A hybrid quantum inspired harmony search algorithm for 0–1 optimization problems. J. Comput. Appl. Math. 2013, 253, 14–25. [Google Scholar] [CrossRef]
  29. Alfailakawi, M.G.; Ahmad, I.; Hamdan, S. Harmony-search algorithm for 2D nearest neighbor quantum circuits realization. Expert Syst. Appl. 2016, 61, 16–27. [Google Scholar] [CrossRef]
  30. Lee, D.; Shon, S.; Lee, S.; Ha, J. Size and Topology Optimization of Truss Structures Using Quantum-Based HS Algorithm. Buildings 2023, 13, 1436. [Google Scholar] [CrossRef]
  31. Lee, K.S.; Geem, Z.W. A new structural optimization method based on the harmony search algorithm. Comput. Struct. 2004, 82, 781–798. [Google Scholar] [CrossRef]
  32. Lee, K.S.; Geem, Z.W.; Lee, S.H.; Bae, K.W. The harmony search heuristic algorithm for discrete structural optimization. Eng. Optim. 2005, 37, 663–684. [Google Scholar] [CrossRef]
  33. Srikanth, D.; Barai, S. Structural optimization using harmony search algorithm. In Soft Computing in Industrial Applications: Algorithms, Integration, and Success Stories; Springer: Berlin/Heidelberg, Germany, 2010; pp. 61–69. [Google Scholar]
  34. Degertekin, S. Improved harmony search algorithms for sizing optimization of truss structures. Comput. Struct. 2012, 92, 229–241. [Google Scholar] [CrossRef]
  35. Kim, B.I.; Kown, J.H. Optimum Design of Truss on Sizing and Shape with Natural Frequency Constraints and Harmony Search Algorithm. J. Ocean Eng. Technol. 2013, 27, 36–42. [Google Scholar] [CrossRef]
  36. Cheng, M.Y.; Prayogo, D.; Wu, Y.W.; Lukito, M.M. A Hybrid Harmony Search algorithm for discrete sizing optimization of truss structure. Autom. Constr. 2016, 69, 21–33. [Google Scholar] [CrossRef]
  37. Talatahari, S.; Goodarzimehr, V.; Shojaee, S. Symbiotic organisms search and harmony search algorithms for discrete optimization of structures. Int. J. Optim. Civ. Eng. 2021, 11, 177–194. [Google Scholar]
  38. Han, K.H.; Kim, J.H. Quantum-inspired evolutionary algorithms with a new termination criterion, H/sub/spl epsi//gate, and two-phase scheme. IEEE Trans. Evol. Comput. 2004, 8, 156–169. [Google Scholar] [CrossRef]
  39. Campos, M.; Krohling, R.A.; Enriquez, I. Bare bones particle swarm optimization with scale matrix adaptation. IEEE Trans. Cybern. 2013, 44, 1567–1578. [Google Scholar] [CrossRef]
  40. Lee, D.; Kim, J.; Shon, S.; Lee, S. An Advanced Crow Search Algorithm for Solving Global Optimization Problem. Appl. Sci. 2023, 13, 6628. [Google Scholar] [CrossRef]
  41. Mahdavi, M.; Fesanghary, M.; Damangir, E. An improved harmony search algorithm for solving optimization problems. Appl. Math. Comput. 2007, 188, 1567–1579. [Google Scholar] [CrossRef]
  42. Savsani, V.J.; Tejani, G.G.; Patel, V.K. Truss topology optimization with static and dynamic constraints using modified subpopulation teaching–learning-based optimization. Eng. Optim. 2016, 48, 1990–2006. [Google Scholar] [CrossRef]
  43. Savsani, V.J.; Tejani, G.G.; Patel, V.K.; Savsani, P. Modified meta-heuristics using random mutation for truss topology optimization with static and dynamic constraints. J. Comput. Des. Eng. 2017, 4, 106–130. [Google Scholar] [CrossRef]
Figure 1. Concept of QHM.
Figure 1. Concept of QHM.
Buildings 13 02132 g001
Figure 2. Comparison of convergence performance according to changes in QHMS: (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Figure 2. Comparison of convergence performance according to changes in QHMS: (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Buildings 13 02132 g002aBuildings 13 02132 g002b
Figure 3. 3D bar graph according to changes in QHMCR and QPAR: (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Figure 3. 3D bar graph according to changes in QHMCR and QPAR: (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Buildings 13 02132 g003
Figure 4. Scatter plot according to changes in ϵ : (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Figure 4. Scatter plot according to changes in ϵ : (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Buildings 13 02132 g004aBuildings 13 02132 g004b
Figure 5. Probability according to changes in ϵ : (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Figure 5. Probability according to changes in ϵ : (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Buildings 13 02132 g005aBuildings 13 02132 g005b
Figure 6. Scatter plot according to changes in θ r : (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Figure 6. Scatter plot according to changes in θ r : (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Buildings 13 02132 g006aBuildings 13 02132 g006b
Figure 7. Probability according to changes in θ r : (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Figure 7. Probability according to changes in θ r : (a) f 01 ; (b) f 02 ; (c) f 03 ; (d) f 04 ; (e) f 05 ; (f) f 06 .
Buildings 13 02132 g007aBuildings 13 02132 g007b
Figure 8. Shape of a 20-bar truss structure.
Figure 8. Shape of a 20-bar truss structure.
Buildings 13 02132 g008
Figure 9. Convergence graph of 20-bar truss structures: (a) QbHS H G ; (b) QbHS R V ; (c) QE.
Figure 9. Convergence graph of 20-bar truss structures: (a) QbHS H G ; (b) QbHS R V ; (c) QE.
Buildings 13 02132 g009
Figure 10. Shape of the 24-bar truss structure.
Figure 10. Shape of the 24-bar truss structure.
Buildings 13 02132 g010
Figure 11. Convergence graph of the 24-bar truss structure: (a) QbHS H G ; (b) QbHS R V ; (c) QE.
Figure 11. Convergence graph of the 24-bar truss structure: (a) QbHS H G ; (b) QbHS R V ; (c) QE.
Buildings 13 02132 g011
Figure 12. Shape of the 72-bar truss structure.
Figure 12. Shape of the 72-bar truss structure.
Buildings 13 02132 g012
Figure 13. Convergence graph of the 72-bar truss structure: (a) QbHS H G ; (b) QbHS R V ; (c) QE.
Figure 13. Convergence graph of the 72-bar truss structure: (a) QbHS H G ; (b) QbHS R V ; (c) QE.
Buildings 13 02132 g013
Table 1. Lookup table for rotation gate.
Table 1. Lookup table for rotation gate.
x i b i f ( x ) < f ( b ) Δ θ s i g n ( α i β i )
α i β i > 0 α i β i < 0 α i = 0 β i = 0
00True00000
00False00000
01True θ P 1−10 ± 1
01False00000
10True θ N 11−1 ± 1 0
10False00000
11True00000
11False00000
1 θ N = θ P
Table 2. Benchmark function for comparison.
Table 2. Benchmark function for comparison.
FunctionNameBoundaryMin t max Qubit
f 01 Sphere function[−100 100]050018
f 02 Ackley’s function[−32 32]080018
f 03 Griewank’s function[−600 600]0100021
f 04 Rastrigin’s function[−5.12 5.12]0200017
f 05 Schwefel’s function 2.26[−500 500]0400022
f 06 Rosenbrock’s function[−30 30]0800018
Table 3. Parameters for QHMS analysis.
Table 3. Parameters for QHMS analysis.
dQHMSQHMCRQPAR ϵ θ r Mea. qbw
tolBW qbw max qbw min
201–1000.90.10.010.0610.951.00.01
Table 4. Parameters for QHMCR and QPAR analysis.
Table 4. Parameters for QHMCR and QPAR analysis.
dQHMSQHMCRQPAR ϵ θ r Mea. qbw
tolBW qbw max qbw min
20100.1–0.90.1–0.90.010.0680.951.00.01
Table 5. Parameters for ϵ analysis.
Table 5. Parameters for ϵ analysis.
dQHMSQHMCRQPAR ϵ θ r Mea. qbw
tolBW qbw max qbw min
20100.90.10.00–0.030.0610.951.00.01
Table 6. Parameters for θ r analysis.
Table 6. Parameters for θ r analysis.
dQHMSQHMCRQPAR ϵ θ r Mea. qbw
tolBW qbw max qbw min
20100.90.10.010.00–0.210.951.00.01
Table 7. Parameter for weight optimization of truss structures.
Table 7. Parameter for weight optimization of truss structures.
AlgorithmParameters
QbHSA Q H M S = 10, Q H M C R = 0.9, Q P A R = 0.1, q u b i t = 18, ϵ = 0.01, θ r = 0.06, M e a . = 2, t o l B W = 0.95, B W Q = 0.3, q b w m a x = 1.0, q b w m i n = 0.01
QEALocal group size = 10, Global migration period = 100, q u b i t = 18, ϵ = 0.01, θ r = 0.06, M e a . = 2
Table 8. Results of the 20-bar truss structure.
Table 8. Results of the 20-bar truss structure.
VariableQbHS HG QbHS RV QE
A 1 262320
A 2 ---
A 3 ---
A 4 --21
A 5 242532
A 6 ---
A 7 ---
A 8 243321
A 9 ---
A 10 ---
A 11 242522
A 12 ---
A 13 272732
A 14 ---
A 15 272725
A 16 ---
A 17 ---
A 18 322725
A 19 ---
A 20 322733
Best (kg)332.503331.211344.095
Mean (kg)485.021422.130478.228
Std92.24761.786103.993
σ m a x  (MPa)177.35177.35179.61
σ m a x c r  (MPa)339.35495.48468.39
δ 4 y m a x  (mm)9.4389.6849.829
f 1  (Hz)80.20279.22678.543
f 2  (Hz)100.004100.141134.305
Table 9. Results of the 24-bar truss structure.
Table 9. Results of the 24-bar truss structure.
VariableQbHS HG QbHS RV QE
A 1 ---
A 2 ---
A 3 ---
A 4 ---
A 5 ---
A 6 ---
A 7 253232
A 8 121212
A 9 12116
A 10 ---
A 11 ---
A 12 849
A 13 ---
A 14 -91
A 15 91212
A 16 -2728
A 17 17--
A 18 ---
A 19 ---
A 20 17--
A 21 17--
A 22 19-
A 23 -8-
A 24 17-8
Best (kg)243.762250.718264.944
Mean (kg)364.978342.582364.060
Std91.22372.28973.645
σ m a x  (MPa)155.94162.50175.43
σ m a x c r  (MPa)129.07111.44111.44
δ 5 y m a x  (mm)3.2671.6573.048
δ 6 y m a x  (mm)3.039.6849.829
f 1  (Hz)32.02430.54933.925
Table 10. Results of the 72-bar truss structure.
Table 10. Results of the 72-bar truss structure.
VariableQbHS HG QbHS RV QE
G 1 686
G 2 888
G 3 988
G 4 101010
G 5 988
G 6 889
G 7 ---
G 8 ---
G 9 999
G 10 888
G 11 ---
G 12 ---
G 13 999
G 14 888
G 15 ---
G 16 ---
Best (kg)549.954551.654551.729
Mean (kg)806.250816.971900.185
Std177.718238.728260.150
σ m a x  (MPa)86.5981.9186.78
σ m a x c r  (MPa)133.78133.78133.78
δ m a x  (mm)2.9422.9682.932
f 1  (Hz)4.0084.0134.008
f 3  (Hz)6.8836.8836.940
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lee, S.; Ha, J.; Shon, S.; Lee, D. Weight Optimization of Discrete Truss Structures Using Quantum-Based HS Algorithm. Buildings 2023, 13, 2132. https://doi.org/10.3390/buildings13092132

AMA Style

Lee S, Ha J, Shon S, Lee D. Weight Optimization of Discrete Truss Structures Using Quantum-Based HS Algorithm. Buildings. 2023; 13(9):2132. https://doi.org/10.3390/buildings13092132

Chicago/Turabian Style

Lee, Seungjae, Junhong Ha, Sudeok Shon, and Donwoo Lee. 2023. "Weight Optimization of Discrete Truss Structures Using Quantum-Based HS Algorithm" Buildings 13, no. 9: 2132. https://doi.org/10.3390/buildings13092132

APA Style

Lee, S., Ha, J., Shon, S., & Lee, D. (2023). Weight Optimization of Discrete Truss Structures Using Quantum-Based HS Algorithm. Buildings, 13(9), 2132. https://doi.org/10.3390/buildings13092132

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop