Abstract

This work introduces the Proficient Simulated Annealing Design Agent Model (PSADA), a cognitively inspired, agent-based model of engineering configuration design. PSADA models different proficiency agents using move selection heuristics and problem space search strategies, both of which are identified and extracted from prior human subject studies. The model is validated with two design problems. Agents are compared to human designers and show the accurate simulation of the behaviors of the different proficiency designers. These behavior differences lead to significantly different performance levels, matching the human performance levels with just one exception. These validated heterogeneous agents are placed into teams and confirmed previous findings that the most proficient member of a configuration design team has the largest impact (positive or negative) on team performance. The PSADA model is introduced as a scalable platform to further explore configuration design proficiency’s role in design team performance and organizational behavior.

1 Introduction

Teams are the cornerstone of engineering organizations; however, much is unknown about how best to construct them. Several frameworks have been used to describe the team composition and creation process [1]. It is not disputed that the composition of individual team member attributes plays an important role in a team’s performance. Still, there are inconsistent results on whether an attribute’s minimum level, maximum level, variance, or average level is the best predictor of team success [2]. The relative contribution framework is one such framework and focuses on how individual team members affect team outcomes relative to the other members of the team [3]. A previous human subject study demonstrated how the proficiency of individual designers contributes to overall team performance in a configuration design problem [4]. It was found that the most proficient member had the most significant effect on team performance metrics, positive or negative. In short, the maximum level of proficiency on the team is the best predictor of team success. This means that increasing the proficiency of the most proficient member of a team has a larger effect than increasing the proficiency of any other team member by the same amount. While this result is insightful, the limitations of a human subject study constrain the conclusions that can be drawn. Because of the high time cost of a human subject study, a computational model of proficiency in engineering configuration design would be beneficial and complement the findings from human subject studies. With a validated computational model, the challenges of limited data and isolation of desired attributes of team members can be mitigated.

A single computational agent represents each designer in a multi-agent model. Autonomous agents are modeled such that they perceive their environment, exhibit goal-directed behavior, and are capable of interacting with other agents [5]. Multi-agent models have been used to model various aspects of design, including design organizations [6,7], design negotiation [8], optimal team structure [9], navigating disrupted design objectives [10], social learning in teams [11], and multidisciplinary design [12]. Recently, heterogeneous agents have been utilized to evaluate how varying characteristics of an individual agent affect team outcomes in multi-agent models of engineering design teams [1315]. While many other multi-agent models have been used to solve various engineering problems, they are designed to either solve the problem or propose new design methods, not to model and study engineering design teams and their problem-solving strategies [1618].

The objective of this work is to model proficiency in engineering configuration design. Configuration design is a class of design tasks where the artifact being designed is constructed from a set of predefined components and obeys a set of constraints [19,20]. A validated model can draw conclusions about how different proficiency designers should be utilized to construct design teams more effectively. Proficiency is a designer’s ability to deal with a specific range of problem [21]. It is related to a person’s general aptitude and the complexity of the task but is most strongly correlated with the person’s experience with the task [22]. To model an attribute, there must be behavior differences for varying levels of that attribute. In the previous human subject study, high-proficiency designers moved quickly from design initialization to fine-tuning and optimization and spent less time building the topology of their design. Meanwhile, low-proficiency designers spent significantly more time developing a topology of their design before moving onto the optimization phase of the design process [4]. This finding fits the existing literature on expert versus novice designer behavior, which can translate to high and low proficiency. A well-established element of design expertise is a solution-focused approach. Experts move to a problem solution faster than less experienced designers and spend less time exploring the problem domain [2325]. The other difference in behavior for various proficiency designers is that different proficiency designers utilize distinctive move selection heuristics [4]. These heuristics were explored using extracted hidden Markov models (HMMs), which are shown to model the sequential design process heuristics effectively [26]. These two findings can be modeled to create cognitively inspired computational agents that mimic the behavior of different proficiency designers and consequently mimic the performance of different proficiency designers.

In summary, this paper aims to introduce and validate a cognitively inspired computational model of proficiency in engineering configuration design. This model can serve as a platform for future experimentation to better understand individual team member proficiency and its effect on team performance. The model implementations described in this work are available to the public for educational and academic research purposes.2

2 The Proficient Simulated Annealing Design Agent Model

This work introduces the Proficient Simulated Annealing Design Agents (PSADA) model. The inspiration for this model is the cognitively inspired simulated annealing teams (CISAT) model introduced by McComb et al. [27]. Like CISAT, each PSADA agent uses the simulated annealing stochastic optimization method to dictate their design process [28]. Simulated annealing describes significant aspects of human problem-solving behavior [29] and has been used to computationally model human behavior on mechanical design problems [27,29]. PSADA differs from CISAT by modeling agents heterogeneously, with varying proficiency levels. To accomplish this goal, PSADA models agents using characteristics that focus on individual designer decisions in the design process. These characteristics are modeled after the behavior of individual designers identified from a previous human subject study. Specifically, two characteristics of individual designers are explicitly modeled to mimic those behavior differences for various proficiency designers. They are:

  1. Design space search strategy—Describes a designer’s tendency to utilize breadth or depth search strategies based on proficiency [4,24].

  2. Move selection heuristics—Designers utilize past relevant experience to guide sequential design decisions [4,26,30].

Both characteristics and how they are modeled are discussed in detail in Secs. 2.1 and 2.2.

The different proficiency PSADA agents are validated against their equivalently proficient human designers on two engineering configuration design problems. After validating the different proficiency PSADA agents, they are placed into teams and interact in the same method as the human subject study. This team structure for configuration design problems mimics the team structure of sub-teams within a larger design team. For example, when a company designs a car, several sub-teams work on specific components, like the suspension, chassis, and powertrain. Each sub-team member can collaborate with their sub-team teammates but could realistically develop a solution independently.

Figure 1 summarizes the structure of the PSADA model for both an individual agent and a team of agents. In the PSADA model, each agent can work on the entire problem. Each agent iteratively works to create a design, shares their solution with teammates at probabilistic, self-determined moments, and probabilistically chooses whether to adopt their teammates’ designs when they are shared. This structure allows for individual agents to be modeled independently from their teammates; thus, the characteristics of individual agents and their effect on team performance can be analyzed. Agents complete all actions and decisions independently, and therefore, the model is easily scaled to include different numbers of agents on a team. Section 2.3 discusses the team-based characteristics in detail.

Fig. 1
(a) Logic flowchart for an individual PSADA agent and (b) logic flowchart for PSADA team framework
Fig. 1
(a) Logic flowchart for an individual PSADA agent and (b) logic flowchart for PSADA team framework
Close modal

2.1 Design Space Search Strategy.

Analogous to high-proficiency designers, expert designers often move more quickly to a solution path [24]. Analogous to low-proficiency designers, novice designers often use a “trial-and-error” approach [31]. It is posited that these differences in design space search strategy are due to experience and understanding of the entire design space. High-proficiency designers can quickly make conjectures about solution paths. In contrast, low-proficiency designers must take time to attempt solution paths and make decisions only after they gain feedback from their attempts [4]. This key insight is one way in which proficiency is modeled in PSADA.

Human decision-makers rarely search until the optimal solution is found and instead search through available solutions only until relevant metrics are reached [32]. Simulated annealing is an optimization technique that does not, however, guarantee that the optimal solution is found. After every action the agents take, the simulated annealing algorithm accepts or rejects the move based on the quality of the new design versus the quality of the old design and the temperature parameter prescribed to the simulated annealing algorithm. If the design is better than the previous iteration, it is accepted. If the design is worse, the algorithm probabilistically accepts the design based on Eq. (1):
(1)

For this equation, q is the quality of the previous iteration, q′ is the quality of the new iteration, T is the temperature parameter, and it assumed that q′ > q. The temperature parameter in simulated annealing decreases based on a cooling schedule throughout the search process. The initial temperature for the algorithm depends on the objective function used.

Within PSADA, high-proficiency designers are modeled with a more aggressive simulated annealing cooling schedule to enable these designers to move more quickly to a deterministic solution path, which mimics how high-proficiency designers move more quickly to a solution path. Middle- and low-proficiency designers are modeled with decreasingly aggressive cooling schedules, which allows for an increasing amount of trial-and-error searches early in the design session. Because the length of the design session is equal for all designers, a more aggressive cooling schedule allows for more iterations on a deterministic path. A geometric cooling schedule is used for each and is shown by Eq. (2),
(2)
where T is the temperature at the current iteration, T' is the temperature at the next iteration, and cp is the cooling constant for the agents’ proficiency. In PSADA, the cooling constants used are 0.9, 0.95, and 0.99 for the high-, middle-, and low-proficiency agents, respectively.

2.2 Move Selection Heuristics.

In any design process, designers use heuristics to guide them in solving the design challenge [33]. These heuristics are learned from past relevant experience with the task at hand. The move selection heuristics can be extracted using HMMs for a configuration design task [4,26]. An HMM represents a first-order Markov process with hidden states [34]. An HMM is represented by a transition matrix, T, and an emission matrix, E. After each step, the T matrix shows the probability of the next step being in each state. At each step, the E matrix shows the probability distribution for taking each action given the current state at that step. Experimental data are required to train an HMM for each level of proficiency. A PSADA agent will probabilistically select their action at each iteration, given their proficiency and current state. After each iteration, a PSADA agent will probabilistically transition to the next state, given their proficiency and current state.

A designer’s proficiency at a design task is most strongly related to their experience with it [22]. In PSADA, each level of proficiency requires a unique HMM trained on past human behavior for that task. Therefore, the task-specific element of proficiency is captured by these move selection heuristics in the PSADA model.

2.3 Team-Based Characteristics.

The PSADA model adopts team-focused characteristics from previous models of design teams that contribute to how agents solve problems as part of a team [27]. The modeling of each characteristic in PSADA is designed to focus on individual agent decision-making instead of team decision-making as in other models. Specifically, the communication in the PSADA model is designed so that agents independently determine when to communicate their solutions with their teammates, and their teammates assess other designs only when they are shared. The team-specific characteristics that are used in this model are:

  1. Organic interaction timing—Interaction occurs at irregular and self-determined intervals [35,36].

  2. Quality bias reduction—Individuals in a team develop multiple solution concepts to avoid premature convergence [37].

  3. Self-bias—Designers are biased toward their own solutions [38].

Team interaction in real engineering teams does not occur at regular intervals [35]. After each iteration, agents in the PSADA model probabilistically interact by deciding to share their current solution with their teammates. The probability of sharing is based on the agents’ relative solution quality. Agents are more likely to share their current iteration if it is better than their previous iteration. This models organic interaction timing. A step function is used for determining the probability of sharing. This step function is shown by Eq. (3) and is determined from the previous human subject study where the overall rate of moves per share is matched. The step function is
(3)

In Eq. (3), Pshare is the probability of sharing the new solution, q′ is the quality of the new solution, and q is the quality of the previous solution.

Individuals in a team develop multiple solution concepts to avoid premature convergence [37]. They tend to focus on the most promising alternatives but do not do so greedily [10]. When a design is shared within a team, the other agents on the team evaluate the design against their own and determine whether to adopt the shared design. Designers must make a binary choice. Quality bias reduction is modeled using logistic regression, which finds the probability of switching designs based on the quality of each alternative. Logistic regression is commonly used for discrete choice modeling. The form that is used for this study is shown by Eqs. (4) and (5), where the relative difference, δr, is found by
(4)
and the probability of switching is found by
(5)

In Eq. (4), qs is the objective function value for the shared solution and qc is the current objective function value for the agent. In Eq. (5), Pswitch is the probability of adopting the new design, and β0 and β1 are the coefficient estimates.

Finally, designers have been found to favor their own designs, even when it is not the best option available to them [38]. Because PSADA models interaction as a binary choice using logistic regression, the probability of switching designs can determine if designers in the original human subject study favored their own designs and, therefore, whether the PSADA model incorporates self-bias.

3 Design Tasks

Proficiency is inherently task-specific. The PSADA model is tested on two different configuration design tasks to ensure that the model is generalizable for modeling proficiency across configuration design problems. Table 1 summarizes the parameters that are used for each problem.

Table 1

PSADA model design parameters for the truss and UAV design implementations

ParameterTruss implementationUAV implementationNotes
Initial temperature (Ti)0.12000Dependent on objective function,
Empirically informed
Cooling constants (cp)High proficiency: 0.99
Med proficiency: 0.95
Low proficiency: 0.90
High proficiency: 0.99
Low proficiency: 0.90
Determines design space search strategy,
Empirically informed
Action selectionHidden Markov model for each proficiency levelHidden Markov model for each proficiency levelDetermines move selection heuristics,
Empirically derived
Sharing probabilitySee Eq. (3)N/A,
Team design not tested
Step function,
Empirically informed
Switching probabilitySee Eq. (5): β1 = 0.74
β2 = 0.42
N/A,
Team design not tested
Logistic regression,
Empirically derived
ParameterTruss implementationUAV implementationNotes
Initial temperature (Ti)0.12000Dependent on objective function,
Empirically informed
Cooling constants (cp)High proficiency: 0.99
Med proficiency: 0.95
Low proficiency: 0.90
High proficiency: 0.99
Low proficiency: 0.90
Determines design space search strategy,
Empirically informed
Action selectionHidden Markov model for each proficiency levelHidden Markov model for each proficiency levelDetermines move selection heuristics,
Empirically derived
Sharing probabilitySee Eq. (3)N/A,
Team design not tested
Step function,
Empirically informed
Switching probabilitySee Eq. (5): β1 = 0.74
β2 = 0.42
N/A,
Team design not tested
Logistic regression,
Empirically derived

3.1 Truss Design.

The insights that inspired the PSADA model are from a previous human subject study conducted by Brownell et al., in which human designers (n = 78) were tasked with designing various truss structures that met specific constraints [4]. Previous studies have shown that truss design is a representative configuration design task, and the results can translate to other configuration design problems [26,39].

To design a truss, designers utilize a sequential design process of building a truss structure to support loads given anchor nodes. Figure 2 shows the truss design problem used in this PSADA work. This problem represents a familiar yet non-trivial problem to mechanical engineers. Any two-dimensional truss problem can be used in the PSADA model.

Fig. 2
Truss design problem that is used. The three anchor nodes and two load-bearing nodes are immovable. The arrows represent a load in the direction the arrow is pointed and the length of the arrow represents the relative load size.
Fig. 2
Truss design problem that is used. The three anchor nodes and two load-bearing nodes are immovable. The arrows represent a load in the direction the arrow is pointed and the length of the arrow represents the relative load size.
Close modal
Every truss design problem begins with building the initial topology of the truss design by adding nodes and members until the truss design is stable and statically determinant so that it can be evaluated. In PSADA, a truss design is instantiated by adding nodes and members sequentially until a stable and statically determinant truss design is reached. Upon instantiation, agents can take one of seven actions at each step in the sequential design process. Table 2 shows the possible actions, along with a brief description of each. This set of actions mimics that of the human subject study. The designers in the human subject study were asked to “design a truss to meet a goal factor of safety (FOS) = 1.25, with as low of a mass (m) as possible.” An objective function is developed for use in PSADA to be consistent with the stated objective for the human subject study, as shown by Eq. (6). The objective function is
(6)
Table 2

Actions and action descriptions for truss design

Action nameAction description
1. Add node and attachAdds a node to the design space and attaches it to the three nearest nodes with members. This ensures that the truss design is stable and statically determinant to evaluate the FOS of the design.
2. Add memberAdds a member between the two nearest unconnected nodes in the design space.
3. Delete nodeProbabilistically selects and deletes a node from the design space and all members attached to it. The probability is proportional to the sum of the factors of safety of each member attached to the joint.
4. Delete memberProbabilistically selects and deletes a single member from the design space. The probability is proportional to the factor of safety of the member.
5. Resize all membersIf a majority of members have a factor of safety below the target factor of safety, all members are increased in size. Otherwise, all members are decreased in size.
6. Resize single memberProbabilistically selects a member with probability proportional to |FOS − FOStarget|. If at least one member is failing, only members that are failing can be selected. Once a member is selected, if FOS > FOStarget the size is decreased, otherwise it is increased.
7. Move nodeRandomly selects a node in the design space and moves it. After a node is selected a direction is randomly chosen and a distance is chosen from a set of five predefined distances.
Action nameAction description
1. Add node and attachAdds a node to the design space and attaches it to the three nearest nodes with members. This ensures that the truss design is stable and statically determinant to evaluate the FOS of the design.
2. Add memberAdds a member between the two nearest unconnected nodes in the design space.
3. Delete nodeProbabilistically selects and deletes a node from the design space and all members attached to it. The probability is proportional to the sum of the factors of safety of each member attached to the joint.
4. Delete memberProbabilistically selects and deletes a single member from the design space. The probability is proportional to the factor of safety of the member.
5. Resize all membersIf a majority of members have a factor of safety below the target factor of safety, all members are increased in size. Otherwise, all members are decreased in size.
6. Resize single memberProbabilistically selects a member with probability proportional to |FOS − FOStarget|. If at least one member is failing, only members that are failing can be selected. Once a member is selected, if FOS > FOStarget the size is decreased, otherwise it is increased.
7. Move nodeRandomly selects a node in the design space and moves it. After a node is selected a direction is randomly chosen and a distance is chosen from a set of five predefined distances.

In Eq. (6), q is the quality of the solution, mt is the target mass based on the specific truss problem, ma is the mass of the current truss, FOSt is the target factor of safety, which is always 1.25, and FOSa is the lowest factor of safety for any single truss member. This objective function is the truss’s strength-to-weight ratio (SWR).

3.2 Unmanned Aerial Vehicle Design.

HyForm is a platform that was developed to study hybrid teams. This platform includes an unmanned aerial vehicle (UAV) design interface, shown in Fig. 3. UAV design in the HyForm platform is a configuration design problem, where the UAVs are comprised of motor–rotor pairs, battery structures, and foils connected in a grid. Designers can construct a UAV and simulate its validity and performance [40]. Table 3 shows the possible actions, along with a brief description of their UAV PSADA implementation.

Fig. 3
HyForm UAV design platform with sample UAV
Fig. 3
HyForm UAV design platform with sample UAV
Close modal
Table 3

Actions and action descriptions for UAV design

Action nameAction description
1. Add motor–rotorAdd a motor–rotor to empty grid location. Place x- and y-axis symmetric motors to maintain balance.
2. Add foilAdd a foil to an empty grid location. If x = y, place 1 additional foil at -x, -y to maintain balance. If xy, place three additional x- and y-axis symmetric foils to maintain balance.
3. Add battery structureAdd a structure to an empty grid location. If x = y, place 1 additional structure at -x, -y to maintain balance. If xy, place three additional x- and y-axis symmetric structures to maintain balance.
4. Delete motor–rotorRandomly select and delete a motor-rotor and the x- and y-axis symmetric counterparts.
5. Delete foilRandomly select and delete a foil and the x- and y-axis symmetric counterparts.
6. Delete battery structureRandomly select and delete a structure and the x- and y-axis symmetric counterparts.
7. Resize motor–rotorRandomly select and resize a motor–rotor and the x- and y-axis symmetric counterparts.
8. Resize foilRandomly select and resize a foil and the x- and y-axis symmetric counterparts.
9. Resize battery structureRandomly select and resize a structure and the x- and y-axis symmetric counterparts.
10. Switch designSwitch to a previously saved UAV. After a better UAV is found, agents have a small probability to save their current design.
Action nameAction description
1. Add motor–rotorAdd a motor–rotor to empty grid location. Place x- and y-axis symmetric motors to maintain balance.
2. Add foilAdd a foil to an empty grid location. If x = y, place 1 additional foil at -x, -y to maintain balance. If xy, place three additional x- and y-axis symmetric foils to maintain balance.
3. Add battery structureAdd a structure to an empty grid location. If x = y, place 1 additional structure at -x, -y to maintain balance. If xy, place three additional x- and y-axis symmetric structures to maintain balance.
4. Delete motor–rotorRandomly select and delete a motor-rotor and the x- and y-axis symmetric counterparts.
5. Delete foilRandomly select and delete a foil and the x- and y-axis symmetric counterparts.
6. Delete battery structureRandomly select and delete a structure and the x- and y-axis symmetric counterparts.
7. Resize motor–rotorRandomly select and resize a motor–rotor and the x- and y-axis symmetric counterparts.
8. Resize foilRandomly select and resize a foil and the x- and y-axis symmetric counterparts.
9. Resize battery structureRandomly select and resize a structure and the x- and y-axis symmetric counterparts.
10. Switch designSwitch to a previously saved UAV. After a better UAV is found, agents have a small probability to save their current design.

An experiment by Zhang et al. published data for individual designers (n = 20) building drones to meet specific objectives [41]. These data on individual designer behavior and performance are used to extract move selection heuristics via HMMs for different proficiency drone designers to build PSADA model agents for UAV design. The first set of objectives for the human subject study was to “construct a fleet of UAVs that [are]: (1) Capable of flying a total distance of 20 miles and (2) capable of flying a total distance of 50 miles.” The objective function for the UAV design problem is, therefore, simply the distance that the UAV can travel.

Two significant groupings are found in this dataset of human designer performance. Therefore, no medium proficiency designers were categorized in this dataset. Human designers are categorized into high- and low-proficiency designers. HMMs are derived from the action sequences of these high- and low-proficiency designers.

4 Results

The results aim to validate the model of engineering proficiency proposed in Sec. 2. This is done by comparing results from the computational model with the previously conducted human subject studies. The only parameters that change across design tasks or for different levels of proficiency are those that are summarized in Table 1.

4.1 Individual Agent Validation.

Proficiency is inherently an individual characteristic. Therefore, the proficiency of individual agents must be validated before they are put into teams. The human subjects designed trusses and UAVs in the same scenarios as the agents, which enables a comparison of PSADA and human designers. As discussed in Secs. 3.1 and 3.2, the human data from the truss design task had three distinct proficiency levels, and the human data from the UAV design task had two distinct proficiency levels. Therefore, the PSADA model implementation has high-, medium-, and low-proficiency agents for the truss design task and high- and low-proficiency agents for the UAV design task. One hundred individual PSADA model agents of each proficiency level are instantiated, and their results are analyzed.

As a reminder, the PSADA model aims to model characteristics of designers with different proficiency levels to create agents that mimic human designer behavior and performance. In other words, the agents’ behavior should match that of human designers with the expectation that if the agents behave like various proficiency designers, they will perform like various proficiency designers.

Therefore, the behavior of PSADA model agents is compared to human designers on the two design problems first. Figures 4 and 5 show the action frequency of the different proficiency agents and the different proficiency designers in the truss and UAV design problems, respectively. The different proficiency human designers exhibit different action frequencies, and move selection heuristics that are modeled in PSADA are able to capture these differences for both the truss and UAV design problems. By visual inspection, the different proficiency agents select actions with similar probabilities to the corresponding proficiency human designers. In the truss design problem, the KL divergence for each proficiency pairing (high–high KL: 0.0128, medium–medium KL: 0.0180, low–low KL: 0.0157) confirms that the distributions are closer than each proficiency non-pair (high agent–medium human KL: 0.1543, medium agent–high human KL: 0.0625, low agent–medium human KL: 0.0633). For succinctness, the KL divergence for the furthest apart distributions is not included. The same is true for the UAV design problem (high–high KL: 0.0050, high agent–low human KL: 0.0739, low–low KL: 0.0073, low agent–high human KL: 0.0646).

Fig. 4
Move frequencies for human designers and PSADA model agents separated by proficiency for the truss design problem
Fig. 4
Move frequencies for human designers and PSADA model agents separated by proficiency for the truss design problem
Close modal
Fig. 5
Move frequencies for human designers and PSADA model agents separated by proficiency for the UAV design problem
Fig. 5
Move frequencies for human designers and PSADA model agents separated by proficiency for the UAV design problem
Close modal

The design space search strategy is the other behavior that is explicitly modeled in PSADA. In the original human subject study, it was discovered that low-proficiency designers spent longer exploring the design space before selecting a design topology and refining it. This was uncovered by analyzing the designs being created, including the number of topology elements in the designs. For truss design, the topology elements are the nodes and the members of the truss. Figure 6 replicates this analysis for the PSADA model agents and compares the results to the human designers. Because the design space search differs for different design problems, the equivalent analysis does not hold the same meaning for the UAV design problem. The topology of a UAV interacts differently for a UAV design problem. The number of topology elements for low- and high-proficiency designers were compared and 300 two-sample t-tests resulted in zero significant differences for this analysis. Thus, the comparison between human and agent designers was not conducted.

Fig. 6
The number of topological elements (truss nodes and members) in the trusses that are being designed by different proficiency (a) PSADA model agents and (b) human designers. For both humans and computational agents, the low-proficiency designers create more topologically complex trusses.
Fig. 6
The number of topological elements (truss nodes and members) in the trusses that are being designed by different proficiency (a) PSADA model agents and (b) human designers. For both humans and computational agents, the low-proficiency designers create more topologically complex trusses.
Close modal

Figure 6(a) shows that the truss designs that are created by the different proficiency agents replicate the results from the different proficiency human designers (Fig. 6(b)). In both cases, all designers add topology elements to their design until approximately action 20, which is the required number of actions required to create the most common truss topology for the design scenario that is used. After building the initial truss topology, higher proficiency human designers and PSADA model agents do not add as many topology elements during the remainder of the design session. Low-proficiency human and agent designers, however, continue to add topology elements to their truss designs.

The different proficiency PSADA model agents were compared to their respective proficiency level human designers using two-sample t-tests. Only the iterations after the 20th action were compared. In total 204 two-sample t-tests were conducted and thus, the Holm–Bonferroni correction is used to set the significance levels. Of the 204 tests, only 11 of them yielded significant results. The tests for the high proficiency designers in action numbers 21–27 (p < 0.00026) account for seven of the significant differences. After this initial difference, the high-proficiency human and agent designed trusses do not have a statistically different number of topology elements. The final four iterations (p < 0.00091) for the low-proficiency designers yielded statistically significant results. In general, these results demonstrate that the agents are not designing trusses that are significantly different from the human designers and that all proficiency level agents are exhibiting the appropriate search behavior.

Upon validating the behavior, the averaged performances of the different proficiency human designers versus the averaged performance of the different proficiency PSADA model agents are compared and shown in Figs. 7 and 8 for both design problems. For the truss problem in Fig. 7, all agents and human designers have a mean quality of their first stable and statically determinant truss slightly below 0.3. The final values of each level of proficiency are consistent between the two as well. Beyond the initial and final quality, the trends for all three proficiency levels are reproduced by the PSADA model, with some minor exceptions. For each iteration, a two-sample t-test is conducted between the mean truss quality for the high, medium, and low proficiency PSADA agents and human designer pairings. One hundred and fifty-four significance tests are conducted, and the Bonferroni correction is used to set the significance level. For the high-proficiency and medium-proficiency designer and agent pairings, the tests failed to find a significant difference between the means of the two populations at any iteration. For the low-proficiency designer and agent pairings, only iterations 9, 12, and 15 resulted in statistically significant differences.

Fig. 7
Truss design quality progression by iteration for (a) human subject designers and (b) PSADA model agents. Shaded areas in the plots represent the 95% confidence interval for the mean of each condition.
Fig. 7
Truss design quality progression by iteration for (a) human subject designers and (b) PSADA model agents. Shaded areas in the plots represent the 95% confidence interval for the mean of each condition.
Close modal
Fig. 8
UAV design quality progression by iteration for (a) human subject designers and (b) PSADA model agents. Shaded areas in the plots represent the 95% confidence interval for the mean of each condition.
Fig. 8
UAV design quality progression by iteration for (a) human subject designers and (b) PSADA model agents. Shaded areas in the plots represent the 95% confidence interval for the mean of each condition.
Close modal

For the UAV design problem in Fig. 8, all human designers and agents are given a starting basic UAV design, and therefore all had a minimum best UAV range of ∼10. Early in the design session, the average design quality is not significantly different between the high- and low-proficiency human designers or agents when using a two-sample t-test at a p < 0.01 significance level. A stricter significance level was used due to the high number of significance tests being conducted. Only after the 224th action for the human designers (p = 0.00180) and the 228th action for the agents (p = 0.00979) do the high-proficiency designers have significantly better designs on average than their lower-proficiency counterparts. The average design quality is never significantly different when comparing the low-proficiency human designers and agents. However, the difference between the design quality for high-proficiency designers and agents is significant late in the design session. The human designers achieved a significantly better UAV quality after the 237th iteration. Despite the difference between high-proficiency agents and human designers, there is still a significant difference between low- and high-proficiency agents.

Although many of the high-proficiency agents achieve high ranges for their UAVs, which is defined as a range > 50 (the high proficiency human designer range criteria), not all are able to achieve this. The UAV design problem has a complex design space with large areas of unstable UAVs with an objective function value of zero. Some agents get trapped and spend long periods in zero-quality design space areas, limiting their ability to produce high-quality designs. The model is built first to mimic the human designer’s behavior and then to determine whether the behavioral differences modeled between different proficiency designers can account for the differences in performance. Ultimately, the PSADA agents are not able to escape these zero-quality regions in the same way as the human designers for this single proficiency level.

Another limiting factor for the UAV design problem agents from mimicking the high-proficiency human designer performance is the limited number of high-proficiency human designers. While there are enough designers to optimize HMMs for high- and low-proficiency agents, the HMMs would improve with additional human subjects. Because the dataset is from published data instead of the author’s work, additional data are unavailable. Despite these issues, the PSADA model can still produce high-proficiency agents with significantly better-quality UAVs than their low-proficiency counterparts.

While Figs. 48 demonstrate the validity of the PSADA model in two different design scenarios, it is also important to demonstrate that the individual characteristics that are modeling proficiency—design space search strategy and move selection heuristics—contribute to these results. To test this, each characteristic is “turned off,” one at a time, so the modified PSADA model can run without the characteristic to show how the results change. The results of this experiment with the truss design problem are in Fig. 9.

Fig. 9
Design quality progression for modified PSADA agents with (a) design space search strategy characteristic “turned off” and (b) move selection heuristics characteristic “turned off.” Shaded areas in the plots represent the 95% confidence interval for the 95% confidence interval for the mean of each condition.
Fig. 9
Design quality progression for modified PSADA agents with (a) design space search strategy characteristic “turned off” and (b) move selection heuristics characteristic “turned off.” Shaded areas in the plots represent the 95% confidence interval for the 95% confidence interval for the mean of each condition.
Close modal

The design space search strategy characteristic determines the cooling constant, cp in Eq. (2), for each level of proficiency. To “turn off” this characteristic, each proficiency level is given the same cooling constant. In the test, each is given a cooling constant of 0.95. The rest of the PSADA model remains unchanged (i.e., the different proficiency agents still had unique HHMs) during this test, and 100 instances of the modified model are run with individual agents. Figure 9(a) shows the results of this test.

The move selection heuristics use HMMs to probabilistically choose an action at every iteration. Each level of proficiency has a different HMM derived from human data. To “turn off” move selection heuristics, each level of proficiency selects a random action at every iteration instead of choosing an action based on the HMMs. The rest of the model remains unchanged (i.e., each proficiency level had a unique cp) during this test, and again 100 instances of the modified model are run with individual agents. Figure 9(b) shows the results of this test.

As Fig. 9 shows that both characteristics significantly impact the design quality progression for all proficiency levels when compared to Fig. 7(b). When both characteristics are used, PSADA agents closely follow human designers’ behavior and performance trends. When one of the characteristics is not modeled, the performance differences are not fully captured. These results demonstrate that the characteristics modeled in PSADA capture the behavioral and performance differences between the different proficiency human designers in configuration design tasks.

4.2 The Relative Contribution of Different Proficiency Designers in Configuration Design Teams.

The previous human subject study found that the most proficient team member of a configuration design team had the largest effect, positive or negative, on team performance. To validate this result, the team design portion of the truss design human subject study is reproduced using PSADA model agents. In the study, designers completed an assessment to determine their proficiency in truss design and then were placed into teams of three to design a truss collectively. The analysis is done with a multiple regression model. A nearly identical model is used to replicate this study with PSADA agents, as shown in Eq. (7). To test the effect of individual team member proficiency on team performance, the model
(7)
is used, where S denotes the set of teams, m1,i is the teams’ lowest proficiency members’ proficiency score, m2,i is the teams’ middle proficiency members’ proficiency score, and m3,i is the teams’ highest proficiency members’ proficiency score. In general, mn,i will represent the team’s nth least proficient member’s proficiency score. The βs are the unknown coefficients and qi is the teams’ final truss quality. The proficiency score for each group of agents is the final mean quality when designing individually, which is shown in Fig. 7(b). This ensures that the relative quality of the different proficiency levels is accounted for in the regression model.

With three different proficiency agents modeled in PSADA, 10 different team configurations are possible for a team of three agents, one with three low proficiency agents, one with two low and one medium proficiency agents, one with two low and one high proficiency agents, etc. Each condition is run for 250 instances, and each team acts for 70 iterations. The results are shown in Fig. 10. The coefficient estimates for β1 (the effect size of the lowest proficiency team member), β2 (the effect size of the medium proficiency team member), and β3 (the effect size of the highest proficiency team member), as shown in Fig. 10(b) show that improving any member of the design team will result in a positive effect on team performance. This result fits the logic that better team members will lead to better team performance. What is not as predictable is that raising the proficiency of the most proficient member has the largest effect on team performance. Therefore, the computational study with PSADA agents confirms the results from the previous human subject study that improving the highest proficiency member of a configuration design team will have the largest positive effect on team performance. This also validates the team characteristics of the PSADA model. With Welch’s t-test, this result is statistically significant when comparing β1 (lowest proficiency team member) and β3 (highest proficiency team member) (p = .0108). However, the differences in effect size are not statistically significant when comparing β2 (medium proficiency team member) and β3 (highest proficiency team member) (p = .7793) or comparing β1 (lowest proficiency team member) and β2 (medium proficiency team member) (p = .0757). The general trend shows that the more proficient a team member is relative to their teammates, the larger their effect on team performance.

Fig. 10
(a) SWR Progression for each team condition throughout the design session. In the legend, low proficiency is represented by “1,” medium proficiency is represented by “2,” and high proficiency is represented by “3.” (b) Coefficient estimates for Eq. (5), with 95% confidence intervals.
Fig. 10
(a) SWR Progression for each team condition throughout the design session. In the legend, low proficiency is represented by “1,” medium proficiency is represented by “2,” and high proficiency is represented by “3.” (b) Coefficient estimates for Eq. (5), with 95% confidence intervals.
Close modal

5 Conclusions

A computational model of configuration design proficiency is introduced in this work. The Proficient Simulated Annealing Design Agents model is shown to accurately model the behavior of human designers in two different configuration design problems. Additionally, the differences in behavior lead to significantly different performance levels, mimicking human performance with just one exception. In a truss design problem, three proficiency levels are validated against human subjects from a previously conducted human subject experiment. In a UAV design problem, two proficiency levels are modeled and compared to human subjects’ behavior and performance. The behaviors of action selection and design space search are accurately modeled in both problems. The performance of all three proficiency levels in a truss design problem is accurately matched. The performance of the low-proficiency designers accurately matches the UAV design problem. The high-proficiency agents perform significantly better than the low-proficiency agents but do not mimic the performance of their human-designer counterparts. The reason for this difference in performance is discussed.

Upon validation, PSADA agents can be placed into teams to collaboratively design configurations with the same effect as their human counterparts. For the truss design, it is confirmed that the highest proficiency team member has the largest effect on team performance. This shows the PSADA model’s ability to computationally model certain engineering design teams. In this introduction to the PSADA model, individual proficiency in configuration design problems is studied. However, now that the model has been validated to mimic human designers of different proficiency levels and design teams, PSADA can be utilized in future experiments to test other characteristics, other engineering design tasks, other team structures, and engineering design organizations that distribute designers to multiple teams simultaneously.

PSADA allows for an in-depth analysis of designer proficiency’s role in configuration design teams. With this model, it can be confirmed that the previously discovered behavior differences between various proficiency designers account for a large portion of the variance between different proficiency human designers in engineering configuration design. The PSADA model does not currently model other characteristics that interact with proficiency. Personality, operational learning, and leadership could affect how low-proficiency designers learn during a design session to become higher-proficiency designers and how they use their knowledge and expertise to influence lower-proficiency teammates positively. Incorporating new cognitive and behavioral characteristics into the PSADA structure is possible and is a fruitful area for future work. Additionally, the PSADA model requires human data for the move selection heuristics. Without a previously conducted human subject study to train HHMs for different proficiency agents, applying the PSADA model to new design tasks would not be possible.

It is impractical to assume that every designer will be of the highest proficiency and that understanding how best to utilize people of different proficiencies can lead to better outcomes. This work shows that the PSADA model is a foundation for a deeper understanding of the role that individual team member proficiency plays in engineering design teams, yielding novel insights into this area.

Footnote

Acknowledgment

This material is based on work supported by the United States Air Force Office of Scientific Research (ASFOR) through Grant No. FA9550-18-1-0088. Any opinions, findings, and conclusions or recommendations expressed in this paper are those of the authors and do not necessarily reflect the views of the sponsors.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

References

1.
Mathieu
,
J. E.
,
Tannenbaum
,
S. I.
,
Donsbach
,
J. S.
, and
Alliger
,
G. M.
,
2014
, “
A Review and Integration of Team Composition Models: Moving Toward a Dynamic and Temporal Framework
,”
J. Manage.
,
40
(
1
), pp.
130
160
.
2.
Sonesh
,
S. C.
,
Lacerenza
,
C.
,
Marlow
,
S.
, and
Salas
,
E.
,
2018
, “What Makes an Expert Team? A Decade of Research,”
The Cambridge Handbook of Expertise and Expert Performance
,
K. A.
Ericsson
,
R. R.
Hoffman
,
A.
Kozbelt
, and
A. M.
Williams
, eds.,
Cambridge University Press
,
Cambridge, UK
, pp.
506
532
.
3.
Steiner
,
I. D.
,
1972
, “Group Performance of Unitary Tasks,”
Group Process and Productivity
,
Academic Press
,
New York
, pp.
14
39
.
4.
Brownell
,
E.
,
Cagan
,
J.
, and
Kotovsky
,
K.
,
2021
, “
Only as Strong as the Strongest Link: The Relative Contribution of Individual Team Member Proficiency in Configuration Design
,”
ASME J. Mech. Des.
,
143
(
8
), p.
081402
.
5.
Wooldridge
,
M.
, and
Jennings
,
N. R.
,
1995
, “
Intelligent Agents: Theory and Practice
,”
Knowl. Eng. Rev.
,
10
(
2
), pp.
115
152
.
6.
Campbell
,
M. I.
,
Cagan
,
J.
, and
Kotovsky
,
K.
,
1999
, “
A-Design: An Agent-Based Approach to Conceptual Design in a Dynamic Environment
,”
Res. Eng. Des.
,
11
(
3
), pp.
172
192
.
7.
Jin
,
Y.
, and
Levitt
,
R. E.
,
1996
, “
The Virtual Design Team: A Computational Model of Project Organizations
,”
Comput. Math. Organ. Theory
,
2
(
3
), pp.
171
195
.
8.
Jin
,
Y.
, and
Lu
,
S. C.-Y.
,
2004
, “
Agent Based Negotiation for Collaborative Design Decision Making
,”
CIRP Ann.
,
53
(
1
), pp.
121
124
.
9.
McComb
,
C.
,
Cagan
,
J.
, and
Kotovsky
,
K.
,
2017
, “
Optimizing Design Teams Based on Problem Properties: Computational Team Simulations and an Applied Empirical Test
,”
ASME J. Mech. Des.
,
139
(
4
), p.
041101
.
10.
McComb
,
C.
,
Cagan
,
J.
, and
Kotovsky
,
K.
,
2015
, “
Rolling with the Punches: An Examination of Team Performance in a Design Task Subject to Drastic Changes
,”
Des. Stud.
,
36
(
C
), pp.
99
121
.
11.
Singh
,
V.
,
Dong
,
A.
, and
Gero
,
J. S.
,
2013
, “
Social Learning in Design Teams: The Importance of Direct and Indirect Communications
,”
Artif. Intell. Eng. Des. Anal. Manuf.
,
27
(
2
), pp.
167
182
.
12.
Hulse
,
D.
,
Tumer
,
K.
,
Hoyle
,
C.
, and
Tumer
,
I.
,
2019
, “
Modeling Multidisciplinary Design with Multiagent Learning
,”
Artif. Intell. Eng. Des. Anal. Manuf.
,
33
(
1
), pp.
85
99
.
13.
Singh
,
H.
,
Cascini
,
G.
, and
McComb
,
C.
,
2021
, “
Comparing Design Outcomes Achieved by Teams of Expert and Novice Designers Through Agent-Based Simulation
,”
Proc. Des. Soc.
,
1
, pp.
661
670
.
14.
Singh
,
H.
,
Cascini
,
G.
,
Casakin
,
H.
, and
Singh
,
V.
,
2019
, “
A Computational Framework for Exploring the Socio-Cognitive Features of Teams and Their Influence on Design Outcomes
,”
Proc. Des. Soc. Int. Conf. Eng. Des.
,
1
(
1
), pp.
1
10
.
15.
Perisic
,
M. M.
,
Martinec
,
T.
,
Storga
,
M.
, and
Gero
,
J. S.
,
2019
, “
A Computational Study of the Effect of Experience on Problem/Solution Space Exploration in Teams
,”
Proc. Des. Soc. Int. Conf. Eng. Des.
,
1
(
1
), pp.
11
20
.
16.
Soria Zurita
,
N. F.
,
Colby
,
M. K.
,
Tumer
,
I. Y.
,
Hoyle
,
C.
, and
Tumer
,
K.
,
2018
, “
Design of Complex Engineered Systems Using Multi-Agent Coordination
,”
ASME J. Comput. Inf. Sci. Eng.
,
18
(
1
), p.
011003
.
17.
Manion
,
C.
,
Soria
,
N. F.
,
Tumer
,
K.
,
Hoyle
,
C.
, and
Tumer
,
I. Y.
,
2016
, “
Designing a Self-Replicating Robotic Manufacturing Factory
,”
American Society of Mechanical Engineering
.
18.
Dimeas
,
A. L.
, and
Hatziargyriou
,
N. D.
,
2005
, “
Operation of a Multiagent System for Microgrid Control
,”
IEEE Trans. Power Syst.
,
20
(
3
), pp.
1447
1455
.
19.
Mittal
,
S.
, and
Frayman
,
F.
,
1989
, “
Towards a Generic Model of Configuration Tasks
,”
IJCAI'89: Proceedings of the 11th international joint conference on Artificial intelligence
, Vol.
2
, pp.
1395
1401
.
20.
Wielinga
,
B.
, and
Schreiber
,
G.
,
1997
, “
Configuration-Design Problem Solving
,”
IEEE Expert
,
49
(
2
), pp.
49
56
.
21.
Bass
,
B. M.
,
1980
, “
Team Productivity and Individual Member Competence
,”
Small Group Res.
,
11
(
4
), pp.
431
504
.
22.
Lance
,
C. E.
,
Hedge
,
J. W.
, and
Alley
,
W. E.
,
1989
, “
Joint Relationships of Task Proficiency With Aptitude, Experience, and Task Difficulty: A Cross-Level, Interactional Study
,”
Hum. Perform.
,
2
(
4
), pp.
249
272
.
23.
Atman
,
C. J.
,
Chimka
,
J. R.
,
Bursic
,
K. M.
, and
Nachtmann
,
H. L.
,
1999
, “
A Comparison of Freshman and Senior Engineering Design Processes
,”
Des. Stud.
,
20
(
2
), pp.
131
152
.
24.
Cross
,
N.
,
2004
, “
Expertise in Design: An Overview
,”
Des. Stud.
,
25
(
5
), pp.
427
441
.
25.
Cross
,
N.
,
2018
, “Expertise in Professional Design,”
The Cambridge Handbook of Expertise and Expert Performance
,
K. A.
Ericsson
,
R. R.
Hoffman
,
A.
Kozbelt
,
and A. M.
Williams
, eds.,
Cambridge University Press
,
Cambridge, UK
, pp.
372
388
.
26.
McComb
,
C.
,
Cagan
,
J.
, and
Kotovsky
,
K.
,
2017
, “
Mining Process Heuristics From Designer Action Data via Hidden Markov Models
,”
ASME J. Mech. Des.
,
139
(
11
), p.
111412
.
27.
McComb
,
C.
,
Cagan
,
J.
, and
Kotovsky
,
K.
,
2015
, “
Lifting the Veil: Drawing Insights About Design Teams From a Cognitively-Inspired Computational Model
,”
Des. Stud.
,
40
, pp.
119
142
.
28.
Kirkpatrick
,
S.
,
Gelatt
,
C. D.
, and
Vecchi
,
M. P.
,
1983
, “
Optimization by Simulated Annealing
,”
Science
,
220
(
4598
), pp.
671
680
.
29.
Cagan
,
J.
, and
Kotovsky
,
K.
,
1997
, “
Simulated Annealing and the Generation of the Objective Function: A Model of Learning During Problem Solving
,”
Comput. Intell.
,
13
(
4
), pp.
534
581
.
30.
Puentes
,
L.
,
Cagan
,
J.
, and
McComb
,
C.
,
2021
, “
Data-Driven Heuristic Induction From Human Design Behavior
,”
ASME J. Comput. Inf. Sci. Eng.
,
21
(
2
), p.
024501
.
31.
Ahmed
,
S.
,
Wallace
,
K. M.
, and
Blessing
,
L. T.
,
2003
, “
Understanding the Differences Between How Novice and Experienced Designers Approach Design Tasks
,”
Res. Eng. Des.
,
14
(
1
), pp.
1
11
.
32.
Simon
,
H. A.
,
1956
, “
Rational Choice and the Structure of the Environment
,”
Psychol. Rev.
,
63
(
2
), pp.
129
138
.
33.
Yilmaz
,
S.
,
Seifert
,
C. M.
, and
Gonzalez
,
R.
,
2010
, “
Cognitive Heuristics in Design: Instructional Strategies to Increase Creativity in Idea Generation
,”
Artif. Intell. Eng. Des. Anal. Manuf.
,
24
(
3
), pp.
335
355
.
34.
Leonard
,
E. B.
,
Ted
,
P.
,
George
,
S.
, and
Norman
,
W.
,
1970
, “
A Maximization Technique Occurring in the Statistical Analysis of Probabilistic Functions of Markov Chains
,”
Ann. Math. Stat.
,
41
(
1
), pp.
164
171
.
35.
Stempfle
,
J.
, and
Badke-Schaub
,
P.
,
2002
, “
Thinking in Design Teams—An Analysis of Team Communication
,”
Des. Stud.
,
23
(
5
), pp.
473
496
.
36.
Eckert
,
C.
,
Maier
,
A.
, and
McMahon
,
C.
,
2005
, “Communication in Design,”
Design Process Improvement
,
J.
Clarkson
, and
C.
Eckert
, eds.,
Springer
,
London
, pp.
232
261
.
37.
Dieter
,
G.
, and
Schmidt
,
L.
,
2009
,
Engineering Design
, 4th ed.,
McGraw-Hill Education
,
New York
.
38.
Nikander
,
J. B.
,
Liikkanen
,
L. A.
, and
Laakso
,
M.
,
2014
, “
The Preference Effect in Design Concept Evaluation
,”
Des. Stud.
,
35
(
5
), pp.
473
499
.
39.
McComb
,
C.
,
Cagan
,
J.
, and
Kotovsky
,
K.
,
2016
, “
Linking Properties of Design Problems to Optimal Team Characteristics
,”
Proceedings of the ASME Design Engineering Technical Conference
,
Charlotte, NC
.
40.
Song
,
B.
,
Soria Zurita
,
N. F.
,
Zhang
,
G.
,
Stump
,
G.
,
Balon
,
C.
,
Miller
,
S. W.
,
Yukish
,
M.
,
Cagan
,
J.
, and
McComb
,
C.
,
2020
, “
Toward Hybrid Teams: A Platform to Understnad Human-Computer Collaboration During the Design of Complex Engineered Systems
,”
Proc. Des. Soc. Des. Conf.
,
1
, pp.
1551
1560
.
41.
Zhang
,
G.
,
Soria Zurita
,
N. F.
,
Stump
,
G.
,
Song
,
B.
,
Cagan
,
J.
, and
McComb
,
C.
,
2021
, “
Data on the Design and Operation of Drones by Both Individuals and Teams
,”
Data Brief
,
36
, p.
107008
.