Engineering Management and Modular Design: A Path to Robust Manufacturing Processes

1. Introduction

To manufacture a product, effective and efficient utilization of technical personnel, physical resources, and information is crucial (Chang, 2005; Sharon et al., 2011). The responsibility for this crucial task is generally assigned to engineering management (Savage et al., 2019). Engineering management is a diverse field that combines the principles of engineering, business administration, and management to plan, design, organize, supervise, and lead technical projects, teams, and organizations. It applies engineering and management principles to ensure that engineering projects are completed on time, within budget, and to the required quality standards. The responsibilities of an engineering manager include coordinating and collaborating with different stakeholders involved in a project, such as engineers, clients, vendors, and regulatory bodies, while managing people, resources, and risk. Engineering management encompasses a broad range of topics, including technology management, quality control, project management, cost analysis, design, procurement, construction, commissioning, project control, quality assurance/quality control, project service, and more. Ultimately, engineering management aims to effectively and efficiently utilize technical personnel, physical resources, and information to produce products or services that meet an organization's goals and objectives (Lannes, 2001; Farr et al., 2001; Daughton, 2017; Savage et al., 2019).

The resource allocation, decision-making processes, and the overall organizational strategy are handled by manufacturing management to ensure effective operations (Womack and Jones, 2003; Chang, 2005; Jacobs et al., 2011). In this regard, manufacturing management can significantly impact engineering management (Ahmad et al., 2017).

Based on the provided definition, manufacturing management refers to the planning, coordination, and control of the production process to ensure that goods are produced efficiently, on time, and to the required quality standards. It involves managing the resources and activities involved in the manufacturing process, including people, materials, equipment, and information. (Stevenson, 2019).

Inefficient manufacturing management can lead to suboptimal use of resources, delays in production schedules, and reduced productivity, negatively impacting engineering management. On the other hand, by implementing effective manufacturing management systems, organizations can optimize resource allocation and streamline decision-making processes. This can result in improved production processes and increased efficiency in utilizing technical personnel and physical resources. (Chase et al., 1998; Jacobs et al., 2004; Slack et al., 2010; Egilmez et al., 2018).

Often, manufacturing environments may experience various rapid and sudden changes, a lot of uncertainty, unforeseen risks, unpredictable events, and/or uncontrollable fluctuations (Mollajan et al., 2023). Considering that manufacturing management assumes the role of overseeing resource allocation, decision-making processes, and the overarching organizational strategy, it is crucial that the manufacturing management is capable enough of identifying, analyzing, and mitigating potential risks that may arise during the production process. This includes risks related to equipment failure, supply chain disruptions, safety hazards, quality issues, and other factors that may impact the production process or product quality (Pallagst et al., 2013; Selviaridis et al., 2014; Egilmez et al., 2018).

The poor design of manufacturing processes can significantly contribute to the occurrence and propagation of operational risks (Rahmani and Thomson, 2011; Mollajan et al., 2021; Mollajan et al., 2022). A poorly designed manufacturing process architecture lacks alignment with good design principles, leading to complexity, low modularity, inefficient resource allocation, and limited flexibility. This deficiency increases the risk of errors, defects, and safety hazards, while hindering adaptation to evolving requirements and unforeseen challenges. Inefficient processes drive up production costs, waste materials, and introduce more errors. They also impede problem identification and resolution, resulting in delays and complications. Poor process design also limits visibility and hampers process improvement over time (Rahmani and Thomson, 2011). On the other hand, A well-designed manufacturing process architecture offers benefits to manufacturing and engineering management. It facilitates bottleneck detection and improvement identification, supported by comprehensive documentation and efficient communication channels. Real-time monitoring enables prompt issue detection and resolution, enhancing overall efficiency (Cheng et al., 2012). Effective design is crucial in minimizing the likelihood of operational risks and their potential impact on manufacturing processes (Mollajan et al., 2023).

The design of manufacturing processes has a cascading effect on engineering management in addition to manufacturing management. In fact, design and management of manufacturing processes have a significant impact on engineering management. Engineering teams bear the responsibility of designing, implementing, and improving these processes, which necessitates specific skills and expertise aligned with the complexity and level of automation employed. Furthermore, the effectiveness of engineering management systems, including quality control and continuous improvement programs, relies on the underlying architecture of manufacturing. Therefore, the design of manufacturing processes architecture falls within the scope of engineering management responsibilities, offering an opportunity for them to make a substantial impact.

Achieving a well-designed architecture for manufacturing processes is crucial, and manufacturing management plays a pivotal role in contributing to engineering management. One of the most effective and critical ways in which manufacturing management can make this contribution is by actively participating in the design phase of a product. By involving manufacturing experts from the outset, potential production and assembly challenges can be identified and addressed proactively, resulting in an optimized manufacturing process. This collaborative approach ensures the development of an efficient and streamlined process, ultimately leading to the creation of high-quality products.

2. Design of Robust Architecture for Manufacturing Processes

When it comes to manufacturing processes, particularly in the face of environmental disturbances, the significance of robustness cannot be overstated. A robust design for manufacturing processes is key to ensuring that a manufacturing process can withstand, adapt to, and recover from any unforeseen disruptions that may arise.

Al-Refaie (2011) contributed to the understanding and modeling of robustness, highlighting the need for comprehensive methodologies that consider the various factors influencing performance. Koc et al. (2011) conducted a study on process parameter variation in worm hydroforming, aiming to optimize quality, robustness, and capability. They employed response surface analysis with a Box-Behnken design to optimize parameters and evaluate process capability using univariate indices such as Cp and Cpk. However, the study did not consider the interrelationships among response variables and noise variables, indicating the need to account for these factors in robustness assessments. The challenges associated with noise variables in robustness studies were highlighted by Mevik et al. (2001) and Robinson, Borror, and Myers (2004). Their work emphasized the limitations of certain designs, such as BBD, in effectively accounting for noise variables, which are crucial for studying robustness.

Scholz-Reiter et al. (2011) suggested improving robustness by enhancing production capacities through resource utilization optimization. They explore how optimizing resource utilization can contribute to robustness in production systems. Karimi and Djokoto (2012) also discussed the topic of robust stability analysis, further emphasizing the importance of accounting for uncertainties and disturbances in manufacturing processes. Kusumoto et al. (2012), they also discuss robust management in manufacturing systems, focusing on product and process robustness. They emphasize the importance of state-space analysis and structured analysis tools in achieving robustness. Alem and Morabito (2012) implement a robust optimization scheduling approach in furniture manufacturing. Their study surpasses deterministic worst-case methods in most situations, demonstrating the effectiveness of robust optimization techniques in improving robustness in manufacturing processes. Meyer et al. (2013) conducted research to analyze the influence of capacity adjustments on the performance robustness of different operational key figures in manufacturing systems. Their study investigated the trade-off between the robustness of operational performance and cost-efficiency, emphasizing the importance of systems being able to withstand fluctuations and disturbances to maintain high performance and profitability. Mondal et al. (2013) focused on the measurement and modeling of robustness in serial multi-stage manufacturing processes. They employed sensitivity analysis and sensitivity matrices to evaluate and measure the robustness of these processes. The authors highlighted the complexities and interactions present in multi-stage processes, emphasizing the need for a comprehensive understanding of robustness measurement. However, they acknowledged that their methods were more suitable for large batch production scenarios and may not be effective in situations such as one-kilogram production. Additionally, Mondal (2013) introduced a modeling approach that incorporated input-output and in-process parameters to assess and model process robustness in manufacturing systems. They demonstrated the application of their approach through a numerical example involving a two-stage worm gear manufacturing process. Malmström et al. (2013) explored the relationship between resource constraints and robustness in smaller companies. Their findings contribute to understanding the challenges faced by these companies in achieving and maintaining high levels of robustness. Montgomery (2013) emphasized the importance of interaction effects and noise variables in achieving robustness in design and manufacturing. He highlighted the need to understand their impact and efficiently incorporate them in experiments. Montgomery also discussed the relationship between process capability and robustness, emphasizing the need for further modeling and understanding. These insights contribute to the understanding of robust design and process robustness. Xiong et al. (2013) conducted research on evaluating the robustness of production systems. They introduced indicators related to material variances, free capacity, and the overlap between slack periods and machine failure. Their studies explored the system's ability to handle disruptions and maintain smooth operations. Additionally, they focused on evaluating robustness in a flexible job shop problem using the concept of float time. Their research provided valuable insights into the impact of machine and tool breakdowns on the robustness of production systems. Becker et al. (2013) defined MS robustness as the ratio of feasible operations with disruptions of a specific machine to the total feasible operations without disruptions. The specific details or context of their research are not mentioned in the text. Sharda and Banerjee (2013) presented a mathematical optimization model for designing a robust manufacturing system (MS). Their model incorporated factors such as processing times, equipment failure and repairs, and product demand to improve performance in the face of uncertainties. The model aimed to minimize makespan, mean work-in-progress (WIP), and the number of machines, considering these uncertainties during the design process. Mondal et al. (2014) critically reviewed robust approaches for measuring and modeling process robustness in product and manufacturing processes. This research contributed to understanding the effectiveness of existing methodologies and techniques in assessing and modeling process robustness. Taguchi's concept of process robustness emphasized designing processes less sensitive to variations and disruptions for consistent performance. Stricker and Lanza (2014) conducted research that supports the importance of robustness in production systems. They utilized the Value-at-Risk (VaR) measure as an indicator of robustness and explored its implications in uncertain environments. Their work contributes to the understanding of robustness and its role in the success of production systems. Bebera and Becker (2014) focused on exploring the relationship between the topological characteristics and dynamic behavior in manufacturing networks. They highlighted the importance of the network's static structure, represented by its topology, and the material flow within the network in influencing factors such as performance and robustness. Varas et al. (2014) contributed to the field of robust optimization by demonstrating its benefits in reducing variability and uncertainty in production schedules. Their work specifically focused on a sawmill scenario, highlighting the effectiveness of robust optimization strategies. Benderbal et al. (2015) emphasized the significance of technical disruptions in production systems. They analyzed the impact of these disruptions on system robustness and focus on optimizing schedules and designs. Tian et al. (2015) introduced sensitivity analysis to analyze the effects of machining parameters in the electrical discharge machining process. Their research focused on understanding the sensitivity of these parameters and their impact on the overall process. Putnik et al. (2015) studied the impact of system size and network structure on the robustness of production networks under changing demand. They defined robustness as the ratio of production rate with perturbations to production rate without perturbations, exploring the relationship between system characteristics and robustness. Their research assessed robustness in manufacturing systems, without specifying the approach or measure used. Benderbal et al. (2015) explored robustness in a reconfigurable production system using an average-based measure that quantified potential delays from machine failure. Their study highlighted various approaches to measuring robustness, including time span, float time, and average-based measures, providing insights into the system's ability to handle variations and disruptions. Stricker et al. (2015), their work reinforces the notion that high robustness is associated with costs. They provide insights into the trade-off between robustness and cost considerations in manufacturing systems. Giannetti and Ransing (2016) improved output quality in the steel foundry industry using tolerance synthesis to assess process robustness. Their robustness indicator, a likelihood ratio, compared the likelihood of achieving the expected output using modified input parameter regions. Validation of their methodology demonstrated its potential in improving output quality. Zhao et al. (2016) investigated the performance evaluation of the axial ventilator in the turbofan engine lubrication system. Their study aimed to assess the sensitivity of the ventilator's performance to different factors. Meyer (2016) examined the role of redundancy in achieving robustness in manufacturing system performance. The study derived specific definitions of robustness and redundancy for the manufacturing context and explored indicators such as nestedness and elementary flux modes. Analysis of different manufacturing system configurations revealed a significant correlation between redundancy and robustness. These findings emphasized the potential of incorporating robustness in the design of manufacturing systems. Boorla and Howard (2016) contribute to the discussion of robust management in manufacturing systems. They highlight the significance of robustness measures and analysis in ensuring the resilience of products and manufacturing processes. Sakhaii et al. (2016) further contribute to the topic of technical disruptions in production systems. Their research explores the optimization of system characteristics and the impact of these disruptions. Moslemi et al. (2017) focus on modeling robustness in multistage manufacturing systems, introducing process capability as a surrogate measure of robustness. They utilize response surfaces and a robust coefficient estimation method to analyze quality characteristics. The results demonstrate the superiority of the robust multistage process and highlight the relationship between robustness and capability in multistage manufacturing. This research contributes to understanding robustness in multistage processes and the use of robustness indices as surrogates. Zeng and Yen (2017), their work focuses on the impact of global value chains on production systems. They highlight how optimizing costs in these chains can increase vulnerability to factors like geo-political disputes and transportation issues, further emphasizing the need for robustness.

Zhang et al. (2017), their research primarily focuses on ensuring robustness in manufacturing systems through technology control. They address the importance of technological advancements and control mechanisms in achieving robustness. Himmiche et al. (2018) conducted a study on the impact of disturbances and uncertainties on meeting deadlines in a production system. They evaluated various machine capacity plans and assessed the likelihood of meeting deadlines under these conditions. Additionally, their research focused on addressing technical disruptions, such as machine and tool breakdowns, and optimizing system characteristics to enhance robustness. Efthymiou et al. (2018) explore robustness in the aerospace domain, applying their metric in a supplier flow workshop. Their approach quantifies capability and enhances controllability across scenarios using likelihood-based indicators. The case study showcases the practical application of their approach, providing valuable insights into production system robustness in the aerospace industry. Pagone et al. (2019) focused on the association between high robustness and costs. They contribute to the understanding of the economic implications of achieving and maintaining robustness in manufacturing operations.

Frederico et al. (2020) addressed the challenges posed by digitization in production systems. They highlight the need for rapid adaptation to new technologies, business models, and competitors, underscoring the importance of robustness in these dynamic environments. Martín et al. (2020) illustrated the benefits of a robust optimization approach in achieving high robustness at lower costs in an automotive supplier context. Their study compared the performance of the robust optimization approach with heuristic planning approaches for production schedules, emphasizing the cost-effectiveness of the robust optimization strategy.

Hyder et al. (2021) tackled uncertainty and sustainability in cloud manufacturing, exploring tools and proposing a Collaborative Cloud Service Platform (CCSP) to enhance supplier-customer-manufacturer relations while considering robustness and sustainability in service composition. Stockmann et al. (2021) emphasized a holistic approach to assess robustness in production systems, considering tradeoffs between performance, resource requirements, and output. Their methodology evaluated robustness from an input-throughput-output perspective, connecting it with adjacent capabilities for a comprehensive understanding of production system robustness. Liang and Li (2022) explored cooperation uncertainty's impact on robustness in manufacturing service systems (MSS). They proposed a methodology treating MSS as a network of service subsystems to identify critical paths and standardize cooperation mechanisms for enhanced robustness.

Exhibit 1 presents a comprehensive overview of the significant aspects of the present study, in comparison to the relevant studies found in the literature. The table showcases the diverse range of studies in the literature, which primarily concentrate on various aspects such as developing mathematical/statistical methods, presenting informative process capability indices, incorporating robust design principles, introducing useful robustness measures, and considering the complexity of interrelationships among manufacturing processes to tackle the issue of robustness in manufacturing process structure. Additionally, a select few of these studies also approach the problem from systems engineering perspective, specifically addressing the concept of system architecture. To achieve such robustness, a system engineering approach is adopted in the present study. This approach recognizes that all activities involved in forming the manufacturing process should be viewed as an integrated system, working together harmoniously to enhance its robustness. It is worth noting that within this category of studies, only a limited number of works propose modularization as an effective approach for achieving robustness in manufacturing processes. Exhibit 1 provides an overview of studies that utilize the Axiomatic Design (AD) theory to modularize the system architecture and contribute to the design of robust manufacturing processes. However, the literature review reveals that the focus has mainly been on the allocated architecture of the system, neglecting the significance of the process architecture in system robustness. According to the AD theory (Suh 2001), the physical components of the system, which are part of the physical architecture, are influenced by a set of process variables (such as men, machines, methods, and materials) belonging to the process architecture. Therefore, this study emphasizes the importance of considering modularity in the process architecture, as it indirectly supports modularity in the allocated architecture. In other words, examining the modularity of the process domain architecture is crucial for the effective allocation of components within the physical system. The literature currently lacks a clear understanding of how modularization of manufacturing processes architecture can effectively enhance the robustness of the system architecture while minimizing system complexity. This gap emphasizes the need for a mathematical model that elucidates the relationship between manufacturing processes architecture robustness and system architecture modularization. Such a model would provide valuable support to systems engineers in designing robust manufacturing processes architecture.

The purpose of this study is to address this significant gap in the literature and develop a mathematical model that precisely defines the proposed approach. By clarifying the relationship between manufacturing processes architecture and system architecture through modularization, the model aims to facilitate the creation of a robust characteristic in the manufacturing processes architecture while ensuring controlled system complexity.

To modularize the system architecture and consider the identified modules more effectively, this study employs the DSM methodology.

The paper is organized as follows: The introduction provides an overview of the study. Section 2 discusses the robust architecture for manufacturing processes, including a literature review and identification of research gaps. Section 3 presents background materials related to the topic. In Section 4, the research method and mathematical model development are described, specifically focusing on the relationship between "modularity" and the property of "robustness of manufacturing processes." Section 5 offers a case study that applies the robustness analysis method to manufacturing processes in Barez Industrial Group, an Iranian company. Finally, Section 6 concludes the paper by discussing the conclusions drawn from the robustness analysis method and providing practical recommendations.

Exhibit 1. Features of the present research comparing with the relevant researches.

Features
	Development of a Statistical/Mathematical Optimization Model	Process Capability Indices Consideration	Utilization of Robust Design Principle	Utilization of Response Surface Methodology	Development of Measurement of Robustness	Complexity of Interrelationship in Manufacturing Processes	Consideration of Manufacturing Processes Architecture	Modularization of Manufacturing Processes Architecture
This Research	✓	✓	✓	✓	✓	✓	✓	✓
Al-Refaie (2011)	✓		✓	✓	✓	✓
Koc et al. (2011)	✓		✓	✓	✓	✓
Mevik et al. (2001)	✓		✓	✓	✓	✓
Robinson, Borror, and Myers (2004)	✓		✓	✓	✓	✓
Nourelfath (2011)	✓		✓	✓	✓	✓
Scholz-Reiter et al. (2011)	✓		✓	✓	✓	✓
Karimi and Djokoto (2012)	✓		✓	✓	✓	✓
Kusumoto et al. (2012)	✓		✓	✓	✓	✓
Alem and Morabito (2012)	✓		✓	✓	✓	✓
Meyer et al. (2013)	✓		✓	✓	✓	✓
Mondal et al. (2013)	✓	✓	✓	✓	✓	✓
Malmström et al. (2013)	✓		✓	✓	✓	✓
Montgomery (2013)	✓	✓	✓	✓	✓	✓
Mondal et al. (2013)	✓		✓	✓	✓	✓
Becker et al. (2013)	✓		✓	✓	✓	✓
Sharda and Banerjee (2013)	✓		✓	✓	✓	✓
Xiong et al. (2013)	✓	✓	✓	✓	✓	✓
Mondal et al. (2014)	✓		✓	✓	✓	✓
Stricker and Lanza (2014)	✓		✓	✓	✓	✓
Bebera and Becker (2014)	✓		✓	✓	✓	✓
Varas et al. (2014)	✓	✓	✓	✓	✓	✓
Benderbal et al. (2015)	✓	✓	✓	✓	✓	✓
Tian et al. (2015)	✓		✓	✓	✓	✓
Putnik et al. (2015)	✓		✓	✓	✓	✓
Benderbal et al. (2015)	✓		✓	✓	✓	✓
Stricker et al. (2015)	✓		✓	✓	✓	✓
Giannetti and Ransing (2016)	✓		✓	✓	✓	✓
Zhao et al. (2016)	✓		✓	✓	✓	✓
Meyer (2016)	✓		✓	✓	✓	✓
Boorla and Howard (2016)	✓		✓	✓	✓	✓
Sakhaii et al. (2016)	✓		✓	✓	✓	✓
Moslemi et al. (2017)	✓		✓	✓	✓	✓
Zeng and Yen (2017)	✓		✓	✓	✓	✓
Zhang et al. (2017)	✓		✓	✓	✓	✓
Himmiche et al. (2018)	✓		✓	✓	✓	✓
Efthymiou et al. (2018)	✓		✓	✓	✓	✓
Pagone et al. (2019)	✓		✓	✓	✓	✓
Frederico et al. (2020)	✓		✓	✓	✓	✓
Martín et al. (2020)	✓		✓	✓	✓	✓
Hyder et al. (2021)	✓		✓	✓	✓	✓
Stockmann et al. (2021)	✓		✓	✓	✓	✓
Liang and Li (2022)	✓		✓	✓	✓	✓

Exhibit 1. Features of the present research comparing with the relevant researches.

Features
	Development of a Statistical/Mathematical Optimization Model	Process Capability Indices Consideration	Utilization of Robust Design Principle	Utilization of Response Surface Methodology	Development of Measurement of Robustness	Complexity of Interrelationship in Manufacturing Processes	Consideration of Manufacturing Processes Architecture	Modularization of Manufacturing Processes Architecture
This Research	✓	✓	✓	✓	✓	✓	✓	✓
Al-Refaie (2011)	✓		✓	✓	✓	✓
Koc et al. (2011)	✓		✓	✓	✓	✓
Mevik et al. (2001)	✓		✓	✓	✓	✓
Robinson, Borror, and Myers (2004)	✓		✓	✓	✓	✓
Nourelfath (2011)	✓		✓	✓	✓	✓
Scholz-Reiter et al. (2011)	✓		✓	✓	✓	✓
Karimi and Djokoto (2012)	✓		✓	✓	✓	✓
Kusumoto et al. (2012)	✓		✓	✓	✓	✓
Alem and Morabito (2012)	✓		✓	✓	✓	✓
Meyer et al. (2013)	✓		✓	✓	✓	✓
Mondal et al. (2013)	✓	✓	✓	✓	✓	✓
Malmström et al. (2013)	✓		✓	✓	✓	✓
Montgomery (2013)	✓	✓	✓	✓	✓	✓
Mondal et al. (2013)	✓		✓	✓	✓	✓
Becker et al. (2013)	✓		✓	✓	✓	✓
Sharda and Banerjee (2013)	✓		✓	✓	✓	✓
Xiong et al. (2013)	✓	✓	✓	✓	✓	✓
Mondal et al. (2014)	✓		✓	✓	✓	✓
Stricker and Lanza (2014)	✓		✓	✓	✓	✓
Bebera and Becker (2014)	✓		✓	✓	✓	✓
Varas et al. (2014)	✓	✓	✓	✓	✓	✓
Benderbal et al. (2015)	✓	✓	✓	✓	✓	✓
Tian et al. (2015)	✓		✓	✓	✓	✓
Putnik et al. (2015)	✓		✓	✓	✓	✓
Benderbal et al. (2015)	✓		✓	✓	✓	✓
Stricker et al. (2015)	✓		✓	✓	✓	✓
Giannetti and Ransing (2016)	✓		✓	✓	✓	✓
Zhao et al. (2016)	✓		✓	✓	✓	✓
Meyer (2016)	✓		✓	✓	✓	✓
Boorla and Howard (2016)	✓		✓	✓	✓	✓
Sakhaii et al. (2016)	✓		✓	✓	✓	✓
Moslemi et al. (2017)	✓		✓	✓	✓	✓
Zeng and Yen (2017)	✓		✓	✓	✓	✓
Zhang et al. (2017)	✓		✓	✓	✓	✓
Himmiche et al. (2018)	✓		✓	✓	✓	✓
Efthymiou et al. (2018)	✓		✓	✓	✓	✓
Pagone et al. (2019)	✓		✓	✓	✓	✓
Frederico et al. (2020)	✓		✓	✓	✓	✓
Martín et al. (2020)	✓		✓	✓	✓	✓
Hyder et al. (2021)	✓		✓	✓	✓	✓
Stockmann et al. (2021)	✓		✓	✓	✓	✓
Liang and Li (2022)	✓		✓	✓	✓	✓

3. Background Material

3.1. Robustness

There is a wide variety of definitions for robustness. According to Carlson and Doyle (2002), robustness is “the maintenance of some desired system characteristic despite fluctuations in the behavior of its component parts or its environment”. Fricke and Schulz (2005) define robustness as “a system’s ability to be insensitive towards changing environments. Robust systems deliver their intended functionality under varying operating conditions without being changed”. We use the definition by Barber and Salido (2011): robustness is the ability of a system to withstand stress, pressure, perturbations, unpredictable changes or variability in its operating environment without loss of functionality.

3.2. Axiomatic Design (AD) Theory

The Axiomatic Design (AD) Theory is a coherent framework rooted in techniques and logical thought processes, aimed at conducting design activities in a systematic and scientific manner (Suh, 1990). The primary elements employed within AD encompass 'Domains', 'Mapping and Zigzagging', the 'Independence (First) Axiom', and the 'Information (Second) Axiom'.

The Axiomatic Design (AD) method is a well-structured approach established to guide design activities across four distinct design domains: the 'Customer Domain', the 'Functional Requirements (FRs) Domain', the 'Physical or Design Parameters (DPs) Domain', and the 'Process Variables (PVs) Domain' (as depicted in Exhibit 2). The customer domain encapsulates customer attributes (CAs). These attributes' elements need to be transformed into functional requirements (FRs) and constraints (Cs) within the functional domain. To fulfil these FRs, design parameters (DPs) are established within the physical domain. Ultimately, to bring the system into production, process variables (PVs) are defined within the process domain. The design process operates as a dynamic interplay between two adjacent domains, cascading in a top-to-bottom manner, commencing at the system level and progressively delving into sublevels of finer detail until a comprehensive representation of the system of interest emerges. The linkage between two contiguous domains is aptly named the 'zigzagging' process, achieved by methodically breaking down the system into smaller design components (as illustrated in Exhibit 2).

Exhibit 2. Design domains in Axiomatic Design (AD) theory

Exhibit 2. Design domains in Axiomatic Design (AD) theory

The AD theory offers two design axioms for gauging the quality of solutions and provides guidance on how to assess and choose from various alternatives (Suh, 2001).

• Independence Axiom: Preserve the autonomy of the FRs (DPs). This axiom emphasizes that each FR (DP) should maintain its autonomy. Therefore, system designers face a pivotal challenge in disassembling intricate systems into subsystems characterized by independent FRs (DPs).
• Information Axiom: Minimize the informational complexity of a design. This axiom asserts that among designs adhering to the Independence Axiom, the one with the lowest informational complexity is regarded as the optimal design choice. The informational complexity associated with a particular FR (DP) quantifies the likelihood of achieving a given FR (DP) successfully (Suh, 2001).

In allocation between functional and physical domains, the informational complexity of a design with a single FR and DP can be expressed as the logarithm of the reciprocal of the probability of system success in fulfilling the FR (P, Equation 1);

$$I={\mathit{log}}_2\frac{1}{P}$$

(1)

Where ‘I’ denotes the information content, measured in bits (logarithm to base 2), and signifies the extent of effort needed to accomplish the given FR.

With the assumption of a uniform probability density function (pdf), the information content ‘I’ can be formulated as depicted in Equation (2);

$$I={\mathit{log}}_X\frac{(AreaofthesystemRange)}{(Areaofthecommonrange)}$$

(2)

The system range's extent is represented as SR, calculated from the probability density function of the FR. In essence, the SR of a particular FR signifies the performance range, PR, linked with the FR in question. The shared range's size, CR, constitutes the segment of the earlier mentioned area within the design range, DR, as illustrated in Exhibit 3. The DR defines the admissible scope linked to the relevant FR.

In the context of an uncoupled design involving n FRs (DPs), the aggregate information is formulated according to Equation (3);

$$I=\sum _{i=1}^n-{\mathit{log}}_2{P}_i=\sum _{i=1}^n{I}_i$$

(3)

where ‘Pi’ represents the probability of FRi (DPi) being fulfilled by DPi (PVi).

Exhibit 3. Probability of success in fulfilling a functional requirement (FR).

Exhibit 3. Probability of success in fulfilling a functional requirement (FR).

3.3. Design Structure Matrix (DSM) Methodology

A Design Structure Matrix (DSM) stands as a potent method for illustrating intricate interdependencies among components within a system. Serving as a versatile and straightforward modelling tool, the DSM aids in the design, depiction, and management of systems. The DSM encapsulates a system's architecture by denoting interactions and dependencies among its elements. Forming a binary square matrix with m columns and rows, along with n non-zero elements, where m symbolizes the nodes and n signifies the directed connections linking them, the DSM captures the intrinsic structure. Within a set of m elements within a system, the corresponding DSM takes shape as an m×m matrix, wherein the system elements occupy row and column headings, following the same order. Precedence relationships manifest in the matrix's off-diagonal entries. In a binary DSM, the value of entry ij (column i, row j) is one (or marked by "X" or "●") when element j depends on element i (i.e., i feeds j), while otherwise, it remains zero (or left vacant). Sequentially aligning the elements as they appear in the DSM, marks below the diagonal denote forward information from element i to j, and those above the diagonal embody feedback information from element j to i.

Central to DSM analysis is the aspiration to minimize feedback loops, achieved by orchestrating the sequencing of processes/elements, thereby facilitating a lower-triangular DSM formation. The marks above the diagonal signify feedback couplings, contrasting with the sub-diagonal marks denoting feed-forward couplings. The above-diagonal markings connote iterative interactions, which should ideally be reduced as much as feasible. However, when feedback marks cannot be eliminated, the implicated elements can be categorized into iterative sub-cycles. In Exhibit 4, for example, items (1, 2, 3) and items (6, 7, 8, 9, 10) fall into two such iterative sub-cycles.

Typically, the realization of this objective necessitates a two-step approach: partitioning and tearing (Eppinger and Browning, 2012; Huang et al., 2015), holds promise for achieving optimal results.

Exhibit 4. Design Structure Matrix (DSM) illustrating the interconnections among elements within the system.

Exhibit 4. Design Structure Matrix (DSM) illustrating the interconnections among elements within the system.

3.4. Modularity in System Architecture

Definition

The literature review underscores the wide array of definitions associated with the term 'Modularity' (Gershenson et al., 2003; Lameche et al., 2017). Generally, the definition of modularity emanates from the characterization of the constituent modules within a system. In essence, these definitions can primarily be categorized into two groups: 'Structural' definitions (These pertain to the relationships and interactions among the system components (Browning, 2001)). 'Architectural' definitions (These address the connections between system components and the functions established at a specific abstraction level within the functional domain (Stone et al., 2000), or through the relationships between models of physical and/or logical representations of systems (Faulconbridge and Ryan, 2014; Kruchten, 1995)). Regarding the notion of modularization within each of these classes, a comprehensive elucidation has been provided by Efatmaneshnik et al. (2018). In general, the process of modularization aims to maximize the dependency within modules while minimizing relations or interactions between modules (Baldwin and Clark, 2003; Gershenson et al., 2003).

3.5. Measures of Modularity for an Individual Architecture (Functional/Physical / Process Architecture)

Various modularity indices assess system architecture. Jung and Simpson (2017), Bonvoisin et al. (2016) reviewed many. Tools like modularity matrix (Dahmus et al., 2001), heuristics (Stone et al., 2000), DSM (Steward, 1981), axiomatic design (Suh, 2001), etc., identify modules. This study chooses DSM for module identification, considering its benefits.

To gauge modularity using the Design Structure Matrix (DSM) methodology, the literature offers a range of metrics (Yu, Yassine, and Goldberg, 2007; Hölttä-Otto and de Weck, 2007; Hölttä-Otto et al., 2012). The research emphasizes studying important interactions between element pairs. Among various modularity measures, focus is on those developed from interactions among the elements. Notably, five key metrics stand out: 'Grouping Capability Index (GCI)', 'Grouping Efficacy', 'Grouping Efficiency', 'Utilization Rate (U)', and 'Modularization Function' (Mikkola, 2006; Mikkola and Gassmann, 2003), which have been frequently utilized (Gupta, 2019; Al-Bashir et al., 2016).

Before proceeding to the measures, establishing notation is beneficial. Let N1 be the count of 1s in diagonal blocks, N2 be the count of 0s in diagonal blocks, N3 be the count of 0s in off-diagonal blocks, and N4 be the count of 1s in off-diagonal blocks. Moreover, N1+N2 indicates the total number of elements (0s and 1s) in diagonal blocks, and N3+N4 represents the total number of elements (0s and 1s) in off-diagonal blocks. The following section provides a detailed outline of the measures.

• Grouping capability index

The Grouping Capability Index (GCI) signifies the capacity of a collection of interdependent entities to establish a module (Equation (4)).

$$GCI=1-\frac{N_4}{N_1+N_4}$$

(4)

2.

Grouping efficacy

Grouping efficacy, denoted as τ, quantifies the ability of a module or group to perform a task. It is contingent on the successful clustering of elements within a group. (Equation (5)).

$$\tau =\frac{1-\frac{N_4}{{N_1+N}_4}}{1-\frac{N_2}{{N_1+N}_4}}$$

(5)

3.

Grouping efficiency

Grouping efficiency, represented as η, serves as a measure of the fit quality for a modular solution within a module. A higher η value indicates a more suitable modular arrangement. η (Equation (6)) symbolizes the relative capacity of components to form a module or group, as well as the compactness of a module.

$$\eta =q.{\eta }_1+(1-q){\eta }_2$$

(6)

${\eta }_1$ and ${\eta }_2$ are expressed by Equations (7) and (8), respectively.

$${\eta }_1=\frac{N_1}{N_1+N_2}$$

(7)

$${\eta }_2=\frac{N_3}{N_3+N_4}$$

(8)

Where, q is a weighting factor (0<q<1);

4.

Utilization rate

Utilization rate, denoted as U, indicates the extent to which components within a given system are dedicated to forming modules and how effectively they contribute to building modular systems (Equation (9)).

$$U=\frac{N_1}{{N_1+N}_2}$$

(9)

5.

Modularization Function

The Modularization Function, denoted as M(u), was initially developed by Mikkola and Gassmann (2003), and subsequently refined by Mikkola (2006).

$$M\left(u\right)=e^{-\frac{u^2}{2Ns\delta }}$$

(10)

Where ‘u’ signifies the count of distinctive elements/components encompassed within an architecture; "N" denotes the total number of standard and unique components, calculated as N = nSTD + u; 𝛿 represents the average intensity of coupling among elements/components and modules; "s" indicates the substitutability factor of unique elements/components within other systems. Additionally, there exists a relationship between 𝛿 (coupling intensity) and "k" (number of interfaces in a module), expressed as 𝛿∼ (𝑘/𝑛), where "n" represents the count of components in a module.

3.6. Measure of Modularity between Two Adjacent Domains

It is possible to determine a quantitative measure of independence for a given design matrix. Suh (2001) proposes two metrics of independence: reangularity, R, and semangularity, S. Reangularity evaluates the orthogonality between DPs and serves as a measure of interdependence.

$$R=sin(\theta )={(1-{Cos}^2(\theta ))}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}$$

(11)

The angle between the DP axes is denoted by θ (Equation 11). To provide additional clarity on the concept of reangularity, let us examine a two-dimensional case.

$$\left\{\begin{array}{c}{FR}_1\\ {FR}_2\end{array}\right\}=\left[\begin{array}{cc}{A}_{11}& A_{12}\\ A_{21}& A_{22}\end{array}\right]·\left\{\begin{array}{c}{DP}_1\\ {DP}_2\end{array}\right\}=\left\{\begin{array}{c}{A}_{11}\\ A_{12}\end{array}\right\}{·DP}_1+\left\{\begin{array}{c}{A}_{21}\\ A_{22}\end{array}\right\}{·DP}_2$$

(12)

where ${\mathrm{A}}_{\mathrm{i}\mathrm{j}}$ represents the entry (i, j) of the DM and relates DP_j to FR_i (i, j = 1, 2). Equation (12) can be rewritten as

$$\left\{\begin{array}{c}{FR}_1\\ {FR}_2\end{array}\right\}=\left[\begin{array}{cc}{A}_{11}& A_{12}\\ A_{21}& A_{22}\end{array}\right]·\left\{\begin{array}{c}{DP}_1\\ {DP}_2\end{array}\right\}=C_1·{(DP}_1)+C_2·{(DP}_2)$$

(13)

where

$C_1=\left\{\begin{array}{c}{A}_{11}\\ A_{12}\end{array}\right\}‚C_2=\left\{\begin{array}{c}{A}_{21}\\ A_{22}\end{array}\right\}$

From linear algebra,

$$\begin{gathered} cos(\theta ) \\ =\frac{C_1{·C}_2}{\left|C_1\right|·\left|C_2\right|} \end{gathered}$$

(14)

For the nth dimensional case, reangularity is

$$R=\prod _{\begin{array}{c}i=1,n-1\\ j=i+1,n\end{array}}{(1-\frac{{(\sum _{k=1}^n{A}_{ki}·A_{kj})}^2}{\left(\sum _{k=1}^n{A}_{ki}^2\right)·(\sum _{k=1}^n{A}_{kj}^2)})}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}$$

(15)

The metric S quantifies the angular relationship between the axes of DPs and FRs, evaluating the correlation between a given FR and any pair of DPs. In terms of coupling, S represents the magnitude of the diagonal elements in a normalized DM. In the case of an nth-dimensional scenario,

$$S=\prod _{j=1}^n(\frac{\left|A_{jj}\right|}{{\left(\sum _{k=1}^n{A}_{kj}^2\right)}^{\raisebox{1ex}{$1$}\!\left/ \!\raisebox{-1ex}{$2$}\right.}})$$

(16)

Both R and S have a maximum value of 1, corresponding to an uncoupled (ideal) design. As the degree of coupling increases, both reangularity and semangularity decrease. In particular, when R or S is equal to 0, the design is considered fully coupled.

3.7. Metrics for System Robustness

System robustness reflects the degree of system coupling. According to the mathematical equations for R and S, it is clear that, as both R and S belong to the interval [0,1], the larger R and S are, the greater R*S is. Concerning the information content of a system, I, the smaller the information content is, the smaller is the variability of the system where (R*S)/I becomes larger and can be employed to quantify the robustness of a system (Gao and He, 2015).

3.8. Robustness Analysis Method

The robustness analysis method looks at the relationship among contract clauses and evaluates alternative sequences. Managers can then choose the best organization of clauses to improve modularization and overall contract robustness. The steps of the method are described below and presented schematically in Exhibit 5:

• Steps 1 Consider manufacturing processes as an integrated system
• Step 2, 3 Develop functional and physical architecture of the system
• Step 4 Design manufacturing processes architecture for robustness
• Steps 5, 6, 7 use DSM to modularize the physical architecture of the system
• Step 8 Map the physical to the functional architecture of the system
• Step 9, 10 Calculate system architecture's information content and robustness
• Steps 11, 12, 13 use AD theory to modularize the mapping between the new physical architecture and the functional architecture
• Step 14, 15 Recompute system architecture's information content and robustness
• Step 16 compute the overall robustness of the system.

Exhibit 5. Methodology for the Robust Analysis Method.

Exhibit 5. Methodology for the Robust Analysis Method.

3.9. Contribution of the Study: Modularity and Robustness

Once robustness is analysed, mapping between a contract’s functional and physical architectures needs to be done, and determination of how well the physical elements satisfy the functional elements needs to be evaluated. A literature review of design robustness has indicated that both physical component redundancy and limiting the degradation of physical components have been mentioned by several authors (Jackson and Ferris, 2013; Jackson et al., 2015). In the present study, the focus is on avoiding propagation of failure/disruption/ disturbance, i.e., suppression of disturbances and redundancy in some of physical component in the physical architecture. To avoid propagation failures in a contract, the functional and physical components should not be coupled and the physical components should be independent. Robustness is achieved by having a significant level of modularity that avoids propagation of noise from an affected/impaired module or clause to other modules. Following is a mathematical investigation that demonstrates that modularity suppresses disturbances in the physical architecture as shown in Exhibit 6, where better modularity leads to better contact clauses, and ultimately, greater project success. The mathematical analysis underpins the use of AD Theory and DSM in the robustness analysis method.

Exhibit 6. Modularity and Robustness.

Exhibit 6. Modularity and Robustness.

3.10. Modularity and System Variance

System robustness can be improved by minimizing variance. In a multi-FR design, the prerequisite for minimizing variance is the satisfaction of the independence axiom. To mathematically consider the relationship for coupling, i.e., the absence of modularity and system variance, consider an n dimensional case

$$\left\{\begin{array}{c}{M}_1^{FR}\\ M_2^{FR}\\ ⋮\\ M_n^{FR}\end{array}\right\}=\left[\begin{array}{cccc}{A}_{11}& A_{12}& \dots & A_{1n}\\ A_{21}& A_{22}& \dots & A_{2n}\\ ⋮& ⋮& ⋮& ⋮\\ A_{n1}& A_{n2}& \dots & A_{nn}\end{array}\right]·\left\{\begin{array}{c}{M}_1^{DP}\\ M_2^{DP}\\ ⋮\\ M_n^{DP}\end{array}\right\}$$

(17)

where ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ (i=1, 2, …, n) represents the i^th module in FRs domain,

${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ (i=1, 2, …, n) represents the i^th module in DPs domain,

$\left[\mathrm{A}\right]$ represents a DM which associates DP-based modules with FR-based modules, and

${\mathrm{A}}_{\mathrm{i}\mathrm{j}}$ (i=1, 2, .., n and j=1, 2, .., n) are DM submatrices that have one or more dimensions where ${\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ (j=1, 2, .., n) is allocated to the ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ (i=1, 2, …, n).

In addition, $[{\mathrm{A}]}^{\mathrm{T}}$ represents the transposition of the matrix.

$${\left[A\right]}^T=\left[\begin{array}{cccc}{A}_{11}& A_{21}& \dots & A_{n1}\\ A_{21}& A_{22}& \dots & A_{2n}\\ ⋮& ⋮& ⋮& ⋮\\ A_{1n}& A_{2n}& \dots & A_{nn}\end{array}\right]$$

(18)

${\mathbf{M}}^{\mathrm{F}\mathrm{R}}$ is expressed as

$${\mathit{M}}^{FR}=\left[A\right]·{\mathit{M}}^{DP}$$

(19)

where ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$represents the domain of FR-based modules (the functional domain), and ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$represents the domain of DP-based modules (the physical domain).

The random variability of ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$, i.e., $\mathsf{\delta }{\mathrm{M}}^{\mathrm{F}\mathrm{R}}$ can be expressed as

$$\mathit{\delta }{\mathit{M}}^{FR}=\left[A\right]·\mathit{\delta }{\mathit{M}}^{DP}$$

(20)

where $\mathsf{\delta }{\mathbf{M}}^{\mathrm{D}\mathrm{P}}$represents random variability of ${\mathbf{M}}^{\mathrm{D}\mathrm{P}}$. It is clear that $\mathsf{\delta }{\mathbf{M}}^{\mathrm{F}\mathrm{R}}$ and $\mathsf{\delta }{\mathbf{M}}^{\mathrm{D}\mathrm{P}}$ are

$$\mathit{\delta }{\mathit{M}}^{FR}=\left[\begin{array}{c}{\delta M}_1^{FR}\\ {\delta M}_2^{FR}\\ ⋮\\ {\delta M}_m^{FR}\end{array}\right]$$

(21)

$$\mathit{\delta }{\mathit{M}}^{DP}=\left[\begin{array}{c}{\delta M}_1^{DP}\\ {\delta M}_2^{DP}\\ ⋮\\ {\delta M}_m^{DP}\end{array}\right]$$

(22)

where${\mathsf{\delta }\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ represents the random variability in ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ (i=1, 2, …, m), and${\mathsf{\delta }\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ represents the random variability in ${\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ (j=1, 2, …, m).

The transposition of equation (20) is

$${(\mathit{\delta }{\mathit{M}}^{FR})}^{\mathit{T}}={(\mathit{\delta }{\mathit{M}}^{DP})}^T·{\left[A\right]}^T$$

(23)

Therefore, the variance of the ${\mathbf{M}}^{\mathrm{F}\mathrm{R}}$ is

$$E\left\{\left(\mathit{\delta }{\mathit{M}}^{FR}\right)·{(\mathit{\delta }{\mathit{M}}^{FR})}^{\mathit{T}}\right\}=E\left\{\left(\left[A\right]·\mathit{\delta }{\mathit{M}}^{DP}\right)·({(\mathit{\delta }{\mathit{M}}^{DP})}^T·{\left[A\right]}^T)\right\}$$

(24)

$$=\left[A\right]·E\left\{\left(\delta {\mathit{M}}^{DP}\right)·(\delta {{(\mathit{M}}^{DP})}^{\mathit{T}})\right\}·{\left[A\right]}^T$$

(25)

where E expresses the mathematical expectation. Hence, the variance of ${\mathbf{M}}^{\mathrm{F}\mathrm{R}}$, $\mathrm{V}\mathrm{a}\mathrm{r}\left({\mathbf{M}}^{\mathrm{F}\mathrm{R}}\right)$, is

$$Var\left({\mathit{M}}^{FR}\right)=\left[A\right]·{\sum }_{{\mathit{M}}^{DP}}·{\left[A\right]}^T$$

(26)

where${\sum }_{{\mathbf{M}}^{\mathrm{D}\mathrm{P}}}$ represents the variance-covariance matrix for ${\mathbf{M}}^{\mathrm{D}\mathrm{P}}$.

Assume that the DP-based modules, ${\mathbf{M}}^{\mathrm{D}\mathrm{P}}$, are statistically independent of each other.

Hence,

$${\Sigma }_{\mathit{D}\mathit{P}}=\left[\begin{array}{cccc}{\Sigma }_{11}& 0& \dots & 0\\ 0& {\Sigma }_{22}& \dots & 0\\ ⋮& ⋮& ⋮& ⋮\\ 0& 0& \dots & {\Sigma }_{mm}\end{array}\right]$$

(27)

where${\mathsf{\Sigma }}_{\mathrm{i}\mathrm{i}}$ represents the variance or variance-covariance matrix of module ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$and (i=1, 2, …, m) when the module is one dimensional or multi-dimensional. Then,

$$\begin{gathered} Var({\mathbf{M}}^{\mathrm{F}\mathrm{R}})=\left[\begin{array}{cccc}{A}_{11}& A_{12}& \dots & A_{1n}\\ A_{21}& A_{22}& \dots & A_{2n}\\ ⋮& ⋮& ⋮& ⋮\\ A_{n1}& A_{n2}& \dots & A_{nn}\end{array}\right].\left[\begin{array}{cccc}{\Sigma }_{11}& 0& \dots & 0\\ 0& {\Sigma }_{22}& \dots & 0\\ ⋮& ⋮& ⋮& ⋮\\ 0& 0& \dots & {\Sigma }_{mm}\end{array}\right].\left[\begin{array}{cccc}{A}_{11}& A_{21}& \dots & A_{n1}\\ A_{21}& A_{22}& \dots & A_{2n}\\ ⋮& ⋮& ⋮& ⋮\\ A_{1n}& A_{2n}& \dots & A_{nn}\end{array}\right] \\ =\left[\begin{array}{cccc}{A}_{11}^2.{\Sigma }_{11}+\dots +A_{1n}^2.{\Sigma }_{mm}& A_{11}.A_{21}.{\Sigma }_{11}+\dots +A_{1n}.A_{2n}.{\Sigma }_{mm}& \dots & A_{11}.A_{n1}.{\Sigma }_{11}+\dots +A_{1n}.A_{nn}.{\Sigma }_{mm}\\ A_{21}.A_{11}.{\Sigma }_{11}+\dots +A_{2n}.A_{1n}.{\Sigma }_{mm}& A_{21}^2.{\Sigma }_{11}+\dots +A_{2n}^2.{\Sigma }_{mm}& \dots & A_{21}.A_{n1}.{\Sigma }_{11}+\dots +A_{2n}.A_{nn}.{\Sigma }_{mm}\\ ⋮& ⋮& ⋮& ⋮\\ A_{n1}.A_{11}.{\Sigma }_{11}+\dots +A_{nn}.A_{1n}.{\Sigma }_{mm}& A_{n1}.A_{21}.{\Sigma }_{11}+\dots +A_{nn}.A_{2n}.{\Sigma }_{mm}& \dots & A_{n1}^2.{\Sigma }_{11}+\dots +A_{nn}^2.{\Sigma }_{mm}\end{array}\right] \end{gathered}$$

(28)

Hence,

$$Var\left(M_1^{FR}\right)=A_{11}^2·{\Sigma }_{11}+A_{12}^2·{\Sigma }_{22}+\dots +A_{1n}^2·{\Sigma }_{mm}$$

(29)

$$Var\left(M_2^{FR}\right)=A_{21}^2·{\Sigma }_{11}+A_{22}^2·{\Sigma }_{22}+\dots +A_{2n}^2·{\Sigma }_{mm}$$

(30)

$$⋮$$

$$Var\left(M_n^{FR}\right)=A_{n1}^2·{\Sigma }_{11}+A_{n2}^2·{\Sigma }_{22}+\dots +A_{nn}^2·{\Sigma }_{mm}$$

(31)

Moreover, the system variance,$\mathrm{V}\mathrm{a}\mathrm{r}\left(\mathrm{S}\mathrm{y}\mathrm{s}\right)$

$$\begin{gathered} =A_{11}^2·{\Sigma }_{11}+A_{12}^2·{\Sigma }_{22}+\dots +A_{1n}^2·{\Sigma }_{mm}+A_{21}^2·{\Sigma }_{11}+A_{22}^2·{\Sigma }_{22}+\dots +A_{2n}^2 \\ ·{\Sigma }_{mm}+A_{n1}^2·{\Sigma }_{11}+A_{n2}^2·{\Sigma }_{22}+\dots +A_{nn}^2·{\Sigma }_{mm} \end{gathered}$$

(32)

$$\begin{gathered} {=(A}_{11}^2+A_{21}^2+\dots +A_{n1}^2){·\Sigma }_{11}+\left(A_{12}^2+A_{22}^2+\dots +A_{n2}^2\right){·\Sigma }_{22}+\dots +(A_{1n}^2 \\ +A_{2n}^2+\dots +A_{nn}^2){·\Sigma }_{mm} \end{gathered}$$

(33)

The larger the off-diagonal entries become, the higher is the system variance, i.e., a higher system variance indicates a more coupled system design. The development of a simple, well-defined equation to clarify the relationship between reangularity/semangularity and the lack of modularity may not be easy. However, the behavior of semangularity and reangularity with respect to the magnitude of off-diagonal entries of the DM are relatively similar. To demonstrate the degree of coupling, i.e., the lack of system modularity, for DMs with higher dimensions, computer simulation was used for four square DMs with sizes of 2×2, 3×3, 5×5, and 10×10 that were randomly generated 1000 times. The magnitude of the off-diagonal elements was gradually increased as the number of runs increased, which showed the effect of coupling on semangularity/reangularity and system variance. Exhibits 7a-d show reangularity versus lack of modularity where reangularity becomes smaller as the size of off-diagonal elements increases. Exhibits 7e-h show system variance versus lack of modularity where system variance becomes larger as the size of off-diagonal elements increases. This underlines the need to reduce the coupling among contract clauses as done in the robustness analysis method.

3.11. Modularization of the Allocated Architecture (Mapping between Two Adjacent Domains)

The allocated architecture of a contract refers to the organization of physical elements that yields a design solution which satisfies the functional requirements of a contract. Therefore, modularity of the allocated architecture is essential to achieving a robust contract. Sometimes, due to possible disruptions in project conditions, redundant clauses are used to increase a contract’s reliability to have a better chance of project success. When creating modules, then, it is important to determine whether there is any clause or module that is redundant since the use of redundant elements may result in a loss of contract modularity.

Concerning the improvement of contract reliability, the present study adopts the notion of redundancy from AD theory where in redundant designs there are more DPs than FRs to satisfy the FRs (Suh, 2001; Park, 2007; Heo and Lee, 2007; Gonçalves-Coelho et al., 2012). For instance, a contract objective may be described by a few clauses in such a way that if the most appropriate clause is not able to satisfy the objective, other associated clauses, as redundant elements, can be used to meet the objective. Nevertheless, contract modularity must be ensured. The following two subsections provide more detail about the modularization of the allocated architecture in the absence and presence of redundant physical modules.

3.12. Modularization of the Allocated Architecture in the Absence of Redundant Physical Modules

To map between DP-based and FR-based domains, consider a multi-FR design

$$M_i^{FR}=f(M_1^{DP}‚M_2^{DP}‚\dots ‚M_n^{DP})$$

(34)

where${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ represents the i^th module in the functional domain (i=1, 2, .., n), and${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ represents the i^th module in the physical domain (i=1, 2, .., n).

Note that ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ or ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$can be a module which can encompass only one element or several fully/partially coupled elements. Using equation (34), the derivative of ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$, i.e., $\mathrm{d}({\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}})$

$$=\frac{\partial f(M_1^{DP}‚M_2^{DP}‚\dots ‚M_n^{DP})}{\partial M_1^{DP}}·d(M_1^{DP})+\dots +\frac{\partial f(M_1^{DP}‚M_2^{DP}‚\dots ‚M_n^{DP})}{\partial M_n^{DP}}·d(M_n^{DP})$$

(35)

Or

$$d{(M}_i^{FR})=\sum _{j=1}^n\frac{\partial f\left(M_1^{DP}‚M_2^{DP}‚\dots ‚M_n^{DP}\right)}{\partial M_j^{DP}}·d(M_j^{DP})=\sum _{j=1}^n\frac{\partial M_i^{FR}}{\partial M_j^{DP}}·d(M_j^{DP})$$

(36)

where, $\mathrm{d}({\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}})$ represents the derivation of ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$. If ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ or ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ is a module of only one element, equation (36) can be utilized straightforwardly. However, when${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ or ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ is a vector which includes several fully/partially coupled elements, equation (36) needs to be expressed as the following. Let ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ and ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$(i=1, 2, …, n) be a vector that includes a specific subset of FRs and DPs, respectively. Then,

$$M_i^{FR}=\left[\begin{array}{c}{FR}_1^i\\ {FR}_2^i\\ ⋮\\ {FR}_{l_1}^i\end{array}\right]$$

(37)

And

$$M_j^{DP}=\left[\begin{array}{c}{DP}_1^j\\ {DP}_2^j\\ ⋮\\ {DP}_{k_1}^j\end{array}\right]$$

(38)

It is assumed that FRs are independent of each other whereas the associated DPs are partially/fully coupled to each other. Therefore, ${\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ can be expressed as

$$M_j^{DP}=\left[\begin{array}{c}{DP}_1^j\\ {DP}_2^j\\ ⋮\\ {DP}_{k_1}^j\end{array}\right]=\left[\begin{array}{cccc}0& x_{1‚2}& \dots & x_{1‚k_1}\\ x_{2‚1}& 0& \dots & x_{2‚k_1}\\ ⋮& ⋮& \dots & ⋮\\ x_{k_{1‚}1}& x_{k_1‚2}& \dots & 0\end{array}\right]·\left[\begin{array}{c}{DP}_1^{M_j}\\ {DP}_2^{M_j}\\ ⋮\\ {DP}_{k_1}^{M_j}\end{array}\right]$$

(39)

where ${\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ represents the j^th module in the physical domain with k₁ DP elements,

${\mathrm{D}\mathrm{P}}_{\mathrm{p}}^{\mathrm{j}}$ represents the p^th member of j^th module (p=1, 2, …, k₁),

${\mathrm{D}\mathrm{P}}_{\mathrm{q}}^{{\mathrm{M}}_{\mathrm{j}}}$ represents the q^th member of the j^th module (q=1, 2, …, k₁), and

${\mathrm{x}}_{\mathrm{p}‚\mathrm{q}}$ represents the entry (p,q) in matrix X and signifies the degree of dependency between${\mathrm{D}\mathrm{P}}_{\mathrm{p}}^{\mathrm{j}}$ and ${\mathrm{D}\mathrm{P}}_{\mathrm{q}}^{{\mathrm{M}}_{\mathrm{j}}}$ (p=1, 2, …, k₁ ; q=1, 2, …, k₁).

Matrix X in which all diagonal elements are zero is employed to express that every element of ${\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$is coupled with the remaining elements. Then,

$$M_j^{DP}=\mathit{X}·{\mathit{D}\mathit{P}}^{M_j};j=1,2,\dots ,n$$

(40)

In addition, assume that A is the design matrix (${\mathrm{l}}_1\times {\mathrm{k}}_1)$ which allocates DPs into FRs that belong to the module ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$. Therefore,

$$M_i^{FR}=\mathit{A}·M_j^{DP};i=1,2,\dots ,n\mathit{and}i=1,2,\dots ,n$$

(41)

Incorporating equation (40) into equation (41)

$$M_i^{FR}=\mathit{A}·\mathit{X}·{\mathit{D}\mathit{P}}^{M_j};i=1,2,\dots ,n\mathit{and}i=1,2,\dots ,n$$

(42)

Let

$${\mathit{B}}_{l_1\times k_1}={\mathit{A}}_{l_1\times k_1}\times {\mathit{X}}_{k_1\times k_1}$$

(43)

Then, equation (43) is

$$M_i^{FR}=\mathit{B}·{\mathit{D}\mathit{P}}^{M_j};i=1,2,\dots ,n\mathit{and}i=1,2,\dots ,n$$

(44)

Using equation (44), the derivative of ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$, i.e., ${\mathrm{d}(\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}})$ (i=1, 2, …, n), is

$${d(M}_i^{FR})=d(\mathit{B}·{\mathit{D}\mathit{P}}^{M_j});i=1,2,\dots ,n\mathit{and}i=1,2,\dots ,n$$

(45)

Where

$$d\left(\mathit{B}·{\mathit{D}\mathit{P}}^{M_j}\right)=\frac{\partial \left(\mathit{B}·{\mathit{D}\mathit{P}}^{M_j}\right)}{\partial ({\mathit{D}\mathit{P}}^{M_j})}·d({\mathit{D}\mathit{P}}^{M_j});i=1,2,\dots ,n\mathit{and}i=1,2,\dots ,n$$

(46)

In order to compute $\frac{\partial \left(\mathbf{B}·{\mathbf{D}\mathbf{P}}^{{\mathrm{M}}_{\mathrm{j}}}\right)}{\partial ({\mathbf{D}\mathbf{P}}^{{\mathrm{M}}_{\mathrm{j}}})}$ , consider $\frac{\partial {\left(\mathbf{B}·{\mathbf{D}\mathbf{P}}^{{\mathrm{M}}_{\mathrm{j}}}\right)}_{\mathbf{s}}}{\partial {({\mathbf{D}\mathbf{P}}^{{\mathrm{M}}_{\mathrm{j}}})}_{\mathrm{t}}}$. Now, ${\left(\mathit{B}·{\mathit{D}\mathit{P}}^{M_j}\right)}_{\mathrm{s}}$is

$${\left(\mathit{B}·{\mathit{D}\mathit{P}}^{M_j}\right)}_s={\sum }_{z=1}^{k_1}{\mathit{B}}_{s‚z}·{{\mathit{D}\mathit{P}}^{M_j}}_z;i=1,2,\dots ,n\mathit{and}i=1,2,\dots ,n$$

(47)

Therefore,

$$\frac{\partial {\left(\mathit{B}·{\mathit{D}\mathit{P}}^{M_j}\right)}_{\mathit{s}}}{\partial {({\mathit{D}\mathit{P}}^{M_j})}_t}=\frac{\partial }{\partial {({\mathit{D}\mathit{P}}^{M_j})}_t}{\sum }_{z=1}^{k_1}{\mathit{B}}_{s‚z}·{{\mathit{D}\mathit{P}}^{M_j}}_z{\mathit{B}}_{s‚t};i=1,2,\dots ,n\mathit{and}i=1,2,\dots ,n$$

(48)

And,

$${(\frac{\partial \mathit{B}·{\mathit{D}\mathit{P}}^{M_j}}{\partial {({\mathit{D}\mathit{P}}^{M_j})}_t})}_s=\frac{\partial {\left(\mathit{B}·{\mathit{D}\mathit{P}}^{M_j}\right)}_{\mathit{s}}}{\partial {({\mathit{D}\mathit{P}}^{M_j})}_t}={\left(\frac{\partial \mathit{B}·{\mathit{D}\mathit{P}}^{M_j}}{\partial {\mathit{D}\mathit{P}}^{M_j}}\right)}_{s‚t}={\mathit{B}}_{s‚t};i=1,2,\dots ,n\mathit{and}i=1,2,\dots ,n$$

(49)

Hence,

$$\frac{\partial (\mathit{B}·{\mathit{D}\mathit{P}}^{M_j})}{\partial ({\mathit{D}\mathit{P}}^{M_j})}=\mathit{B};i=1,2,\dots ,n\mathit{and}i=1,2,\dots ,n$$

(50)

And

$$d\left(\mathit{B}·{\mathit{D}\mathit{P}}^{M_j}\right)=\mathit{B}·d({\mathit{D}\mathit{P}}^{M_j};i=1,2,\dots ,n\mathit{and}i=1,2,\dots ,n$$

(51)

Or

$${d(M}_i^{FR})=\mathit{B}·d({\mathit{D}\mathit{P}}^{M_j});i=1,2,\dots ,n\mathit{and}i=1,2,\dots ,n$$

(52)

where (${\mathit{B})}_{l_1\times k_1}=({\mathit{A})}_{l_1\times k_1}\times {(\mathit{X})}_{k_1\times k_1}$.

$\mathrm{d}{(\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}})$must be within the required tolerance; thus,

$$d{(M}_i^{FR})={(d{(M}_i^{FR}))}_0+\delta {(M}_i^{FR})$$

(53)

where ${\mathrm{d}({(\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}))}_0$ is the desired change in ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ due to a change in ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$, and $\mathsf{\delta }({(\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}))$ is the required tolerance.

In an uncoupled design, $\mathrm{d}{(\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}})$ must be equal to ${\mathrm{d}{((\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}))}_0$. The deviation from this ideal condition, $\mathrm{d}{(\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}})-{\mathrm{d}{((\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}))}_0$, is given by ${\sum }_{\begin{array}{c}j\ne i\\ j=1\end{array}}^{\mathrm{n}}\frac{\partial {\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}}{\partial {\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}}·\mathrm{d}{\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ or $\mathrm{B}.\mathrm{d}({\mathrm{D}\mathrm{P}}^{{\mathrm{M}}_{\mathrm{j}}})$ where ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ or ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ is a module of only one element (case I) or a vector (case II) that includes several fully/partially coupled elements. For an uncoupled system, the required tolerance, $\mathsf{\delta }\left({\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}\right)$ must be ≥${\sum }_{\begin{array}{c}j\ne i\\ j=1\end{array}}^{\mathrm{n}}\frac{\partial {\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}}{\partial {\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}}·\mathrm{d}{\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ or ≥$\mathrm{B}.\mathrm{d}\left({\mathrm{D}\mathrm{P}}^{{\mathrm{M}}_{\mathrm{j}}}\right)$ in case I or II.

For case I,

$$\delta ({(M}_i^{FR}))>\sum _{\begin{array}{c}j=1\\ j\ne i\end{array}}^n\frac{\partial M_i^{FR}}{\partial M_j^{DP}}·d{M}_j^{DP}$$

(54)

and for case II,

$$\delta ({(M}_i^{FR}))>\mathit{B}·d({\mathit{D}\mathit{P}}^{M_j})$$

(55)

Specifically, for case I, consider ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}=\sum _{\mathrm{j}=1}^{\mathrm{n}}{\mathrm{A}}_{\mathrm{i}\mathrm{j}}·{\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$, then,

$$d{(M}_i^{FR})=\sum _{j=1}^n(\frac{\partial }{\partial M_j^{DP}}\left(\sum _{k=1}^n{A}_{ik}·M_k^{DP}\right))·d{M}_j^{DP}$$

(56)

Or,

$$d{(M}_i^{FR})=\sum _{j=1}^n(\sum _{k=1}^n\frac{\partial }{\partial M_j^{DP}}(A_{ik}·M_k^{DP}))·d{M}_j^{DP}$$

(57)

Using equation (54)

$$d{(M}_i^{FR})=\sum _{j=1}^n(\sum _{k=1}^n(A_{ik}·\frac{\partial M_k^{DP}}{M_j^{DP}}+\frac{\partial A_{ik}}{M_j^{DP}}·M_k^{DP}))d{M}_j^{DP}$$

(58)

Specifically, for case II, consider ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}=\sum _{\mathrm{j}=1}^{{\mathrm{l}}_1}{\mathrm{B}}_{\mathrm{i}\mathrm{j}}·{\mathrm{D}\mathrm{P}}_{\mathrm{j}}^{{\mathrm{M}}_{\mathrm{j}}}$, then,

$${d(M}_i^{FR})=\sum _{j=1}^{l_1}{B}_{ij}·{d(DP}_j^{M_j})$$

(59)

Or,

$$d{(M}_i^{FR})=\sum _{j=1}^n(\sum _{k=1}^n\frac{\partial }{\partial M_j^{DP}}(A_{ik}·M_k^{DP}))·d{M}_j^{DP}$$

(60)

Therefore,

$$d{(M}_i^{FR})=\sum _{j=1}^n(\sum _{k=1}^n(A_{ik}·\frac{\partial M_k^{DP}}{M_j^{DP}}+\frac{\partial A_{ik}}{M_j^{DP}}·M_k^{DP}))d{M}_j^{DP}$$

(61)

Equation (61) provides a criterion for determining when the off-diagonal elements of a basic matrix can be disregarded.

3.13. Independence and Tolerance Theorem

An independence and tolerance theorem can be formulated as follows: In case I, a design is considered uncoupled when the specified tolerance, $\mathsf{\delta }{(\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}})$, exceeds...

$$d{(M}_i^{FR})=\sum _{j=1}^n(\sum _{k=1}^n(A_{ik}·\frac{\partial M_k^{DP}}{M_j^{DP}}+\frac{\partial A_{ik}}{M_j^{DP}}·M_k^{DP}))d{M}_j^{DP}$$

(62)

Therefore, the non-diagonal elements of the design matrix can be disregarded. If the $A_{ik}$ (k= 1, 2, ..., n), which represents the stiffnesses relating changes of $M_i^{FR}$ to the changes of $M_k^{DP}$, are constant, then...

$$d{(M}_i^{FR})=\sum _{j=1}^n(\sum _{k=1}^n(A_{ik}·\frac{\partial {(M}_k^{DP}}{{(M}_j^{DP}}))d{(M}_j^{DP})$$

(63)

If $\sum _{\mathrm{j}=1}^{\mathrm{n}}(\sum _{\mathrm{k}=1}^{\mathrm{n}}({\mathrm{A}}_{\mathrm{i}\mathrm{k}}·\frac{\partial {(\mathrm{M}}_{\mathrm{k}}^{\mathrm{D}\mathrm{P}}}{{(\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}}))\mathrm{d}{(\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}})<\mathsf{\delta }({(\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}})$, then, the design is uncoupled.

3.14. B. Modularization of the Allocated Architecture in Presence of Redundant Physical Modules

$$M^{FR}=f\left(M_1^{DP}‚M_2^{DP}‚\dots ‚M_n^{DP}\right)$$

(64)

In this case, there are more ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ than ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}\mathrm{s}$, resulting in random variability (or noise) introduced by the excess ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ that can affect the variability of ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$. It's crucial for the variability to remain within the design range to ensure proper functionality. Therefore,

$$\mathsf{\delta }{\mathrm{M}}^{\mathrm{F}\mathrm{R}}=\frac{\partial f(M_1^{DP}‚M_2^{DP}‚\dots ‚M_n^{DP})}{\partial M_i^{DP}}∆M_i^{DP}+\sum _{\begin{array}{c}j=1\\ i\ne j\end{array}}^n\frac{\partial f(M_1^{DP}‚M_2^{DP}‚\dots ‚M_n^{DP})}{\partial M_i^{DP}}·\delta M_j^{DP}$$

(65)

The random variability in ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$ is denoted as $\mathsf{\delta }{\mathrm{M}}^{\mathrm{F}\mathrm{R}}$, and the variability of ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ is represented by $∆{\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$In an ideal one-${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$ design, only a single ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$ is necessary, as any additional ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ introduce the possibility of random variability. The first term on the right-hand side (RHS) of equation (65) accounts for the compensation effect of the primary ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$, which is selected to vary ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$, while the second term represents the cumulative random variability or noise introduced by any extra ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$. To ensure proper functionality, it's crucial to adjust the error term $∆{\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ in two possible ways to achieve $∆{\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$.

• Setting the coefficient ($\frac{\partial \mathrm{f}}{\partial {\mathrm{M}}^{\mathrm{D}\mathrm{P}}}$) that corresponds to any additional ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ to zero involves eliminating all random variability from every ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$, which is not practical given the variability that arises during manufacturing and assembly.
• To address the issue of random variability in ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$, it's necessary to set all values of ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ except for the chosen ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ to be fixed during operations, allowing that specific ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$ to vary ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$. This step is equivalent to setting the second term on the RHS of equation (65) as a constant. Once this adjustment is made, the next step involves adjusting $∆{\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ to compensate for any accumulated errors by setting the first and second terms on the RHS equal to one another. Then,

  $$\frac{\partial f(M_1^{DP}‚M_2^{DP}‚\dots ‚M_n^{DP})}{\partial M_i^{DP}}∆M_i^{DP}=-\sum _{\begin{array}{c}j=1\\ i\ne j\end{array}}^n\frac{\partial f(M_1^{DP}‚M_2^{DP}‚\dots ‚M_n^{DP})}{\partial M_i^{DP}}·\delta M_j^{DP}$$

  (66)

This approach eliminates the random fluctuations that may arise from additional ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$ units used during operations, while ensuring that the tolerance range specified in the design is met for ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$. The primary ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ unit must be capable of compensating for any random variations resulting from the use of extra ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$ units.

3.15. B. 2. Multi-${\mathit{M}}^{\mathit{F}\mathit{R}}$ Design

In a multi-$M^{FR}$ design where the number of $M^{DP}$ units exceeds that of $M^{FR}$ units, the design becomes redundant and each $M^{FR}$ unit may be affected by multiple $M^{DP}$ units. Consequently,

$$\{M^{FR}\}=[\mathit{Square}\mathit{DM}\left]\right\{M^{DP}\}+[\mathit{Extra}\mathit{Matrix}\left]\right\{M^{DP}{\}}^{\mathit{Extra}}$$

(67)

Let's consider the following scenario: {${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$} represents the ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ selected to fulfill the vector {${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$}, with the number of ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ in {${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$} being equal to the number of ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}\mathrm{s}$ in {${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$}. The vector {${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$}^Extra, which appears as the second term on the right-hand side of equation (67), contains redundant ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ that remain after {${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$} satisfies {${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$}, but these extra ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ may not be explicitly known. The square matrix [Square DM] represents the relationship between {${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$} and {${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$}, while [Extra Matrix] describes how the extra ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ affect the ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}\mathrm{s}$. To ensure compliance with the independence axiom, [Square DM] should be either diagonal or triangular. On the other hand, [Extra Matrix] can take any form, including a full matrix, and still satisfy the independence axiom as long as [Square DM] is diagonal or triangular. Specifically, let's look at a special case of equation (67) where we are dealing with a multi-${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$ design consisting of k ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}\mathrm{s}$ and numerous ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$.

$$\left\{\begin{array}{c}{M}_1^{FR}\\ M_2^{FR}\\ ⋮\\ M_k^{FR}\end{array}\right\}=\left[\begin{array}{cccc}{A}_{11}& 0& 0& 0\\ 0& A_{22}& 0& 0\\ ⋮& ⋮& ⋮& ⋮\\ 0& 0& 0& A_{kk}\end{array}\right]·\left\{\begin{array}{c}{M}_1^{DP}\\ M_2^{DP}\\ ⋮\\ M_k^{DP}\end{array}\right\}+\left[\begin{array}{cccccc}X& X& X& X& \dots & X\\ X& X& X& X& \dots & X\\ ⋮& ⋮& ⋮& ⋮& \dots & ⋮\\ X& X& X& X& \dots & X\end{array}\right]·\left\{\begin{array}{c}{M}_{k+1}^{DP}\\ M_{k+2}^{DP}\\ ⋮\\ M_n^{DP}\end{array}\right\}$$

(68)

where ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ represents the i^th (i=1, 2, …, k) module in the functional domain,

${\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ represents the j^th (j=1, 2, …, n) module in the physical domain,

${\mathrm{D}\mathrm{P}}_{\mathrm{p}}^{\mathrm{j}}$ represents the p^th member of j^th module in the physical domain for (p=1, 2, …, k₁), and${\mathrm{A}}_{\mathrm{i}\mathrm{i}}$ represents the sub-design matrix within DM and relates DPs of ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ to the associated FRs of ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ to satisfy them.

The first term on the right-hand side (RHS) of Equation (68) represents an uncoupled design, while the second term depicts a fully coupled design. However, even with the presence of the second term, this design can still be treated as an uncoupled design if we establish a fixed value for ${\mathrm{M}}_{\mathrm{k}+1}^{\mathrm{D}\mathrm{P}}$ through ${\mathrm{M}}_{\mathrm{n}}^{\mathrm{D}\mathrm{P}}$. We can then use ${\mathrm{M}}_1^{\mathrm{D}\mathrm{P}}$, ${\mathrm{M}}_2^{\mathrm{D}\mathrm{P}}$,… and ${\mathrm{M}}_{\mathrm{n}}^{\mathrm{D}\mathrm{P}}$ to ensure the independence of ${\mathrm{M}}_1^{\mathrm{F}\mathrm{R}}$, ${\mathrm{M}}_2^{\mathrm{F}\mathrm{R}}$, … and ${\mathrm{M}}_{\mathrm{k}}^{\mathrm{F}\mathrm{R}}$. In such cases, the design range of ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$ would be:

$$∆M_i^{FR}=\frac{\partial M_i^{FR}}{\partial M_i^{DP}}·∆M_i^{DP}+{\sum }_{\begin{array}{c}j=k+1\end{array}}^{j=n}\frac{\partial M_i^{FR}}{\partial M_j^{DP}}·\delta M_j^{DP}‚$$

(69)

and

$$∆M_i^{FR}=\frac{\partial M_i^{FR}}{\partial M_i^{DP}}·∆M_i^{DP}+\sum _{\begin{array}{c}j=k+1\end{array}}^{j=n}\left[ExtraTerms\right]$$

(70)

When the diagonal design matrix of the first term on the right-hand side (RHS) of Equation (68) is substituted with a triangular matrix in Equation (71), the resultant design can be treated as a decoupled design

$$\left\{\begin{array}{c}{M}_1^{FR}\\ M_2^{FR}\\ ⋮\\ M_k^{FR}\end{array}\right\}=\left[\begin{array}{cccc}{A}_{11}& 0& \dots & 0\\ A_{21}& A_{22}& \dots & 0\\ ⋮& ⋮& ⋮& ⋮\\ A_{k1}& A_{k2}& A_{k3}& A_{kk}\end{array}\right]·\left\{\begin{array}{c}{M}_1^{DP}\\ M_2^{DP}\\ ⋮\\ M_k^{DP}\end{array}\right\}+\left[\begin{array}{cccccc}X& X& X& X& \dots & X\\ X& X& X& X& \dots & X\\ ⋮& ⋮& ⋮& ⋮& ⋮& ⋮\\ X& X& X& X& X& X\end{array}\right]·\left\{\begin{array}{c}{M}_{k+1}^{DP}\\ M_{k+2}^{DP}\\ ⋮\\ M_n^{DP}\end{array}\right\}$$

(71)

Random variability in ${\mathrm{M}}_{\mathrm{k}+1}^{\mathrm{D}\mathrm{P}}$ through ${\mathrm{M}}_{\mathrm{n}}^{\mathrm{D}\mathrm{P}}$ can be counteracted by using fixed values for ${\mathrm{M}}_1^{\mathrm{D}\mathrm{P}}$, ${\mathrm{M}}_2^{\mathrm{D}\mathrm{P}}$, …, and ${\mathrm{M}}_{\mathrm{k}}^{\mathrm{D}\mathrm{P}}$, provided that any extra terms are fixed first and remain the same across all locations where they are used. In this manner, the design range $∆{\mathrm{M}}^{\mathrm{F}\mathrm{R}}$ may be achieved using ${\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ (j=1, 2, …, k), which is:

$$∆M_i^{FR}=\frac{\partial M_i^{FR}}{\partial M_i^{DP}}·∆M_i^{DP}+\sum _{\begin{array}{c}j=1\\ j\ne i\end{array}}^k\frac{\partial M_i^{FR}}{\partial M_j^{DP}}·∆M_j^{DP}+\sum _{\begin{array}{c}j=k+1\end{array}}^n\frac{\partial M_i^{FR}}{\partial M_k^{DP}}·\delta M_k^{DP}$$

(72)

$$=\frac{\partial M_i^{FR}}{\partial M_i^{DP}}·∆M_i^{DP}+\sum _{\begin{array}{c}j=1\\ j\ne i\end{array}}^k\frac{\partial M_i^{FR}}{\partial M_j^{DP}}·∆M_j^{DP}+(ExtraTerms)$$

(73)

${\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$ is the primary controller for ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$. ${\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ are secondary controllers chosen to satisfy ${\mathrm{M}}_{\mathrm{j}}^{\mathrm{F}\mathrm{R}}$, with a secondary effect on ${\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$ due to off-diagonal elements. ${\mathrm{M}}_{\mathrm{k}}^{\mathrm{D}\mathrm{P}}$ is not a primary controller and is the source of random variability. $∆{\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ (j=1, 2, …, k) should be set first according to the sequence defined by the triangular matrix to control the corresponding ${\mathrm{M}}_{\mathrm{j}}^{\mathrm{F}\mathrm{R}}$ (j=1, 2, …, k). The variability caused by the extra terms can be fixed by eliminating $\mathsf{\delta }{\mathrm{M}}_{\mathrm{j}}^{\mathrm{D}\mathrm{P}}$ (j=k, k+1, …, n) through fixing ${\mathrm{M}}_{\mathrm{k}}^{\mathrm{D}\mathrm{P}}$. Finally, $∆{\mathrm{M}}_{\mathrm{i}}^{\mathrm{F}\mathrm{R}}$can be controlled using $∆{\mathrm{M}}_{\mathrm{i}}^{\mathrm{D}\mathrm{P}}$. If Square DM, the first term on the right-hand side of equation (73), is a full matrix, the design is coupled, which violates the independence axiom. In this case, it is not possible to satisfy the independence of each ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$ when unknown ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ introduce random variability.

The previous discussion holds significant implications for designing a robust system that satisfies multiple ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}\mathrm{s}$ simultaneously. In a multi-${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$ design, the selection of ${\mathrm{M}}_1^{\mathrm{D}\mathrm{P}}$, ${\mathrm{M}}_2^{\mathrm{D}\mathrm{P}}$, and ${\mathrm{M}}_3^{\mathrm{D}\mathrm{P}}$ must meet the same set of conditions. However, in a decoupled multi-${\mathrm{M}}^{\mathrm{F}\mathrm{R}}$ design, we must also satisfy an additional condition for the elements of a triangular design matrix.

$${(\frac{\partial f}{\partial M^{DP}})}_{diagonal}\gg {(\frac{\partial f}{\partial M^{DP}})}_{off-diagonal}$$

(74)

If equation (74) is not satisfied, it indicates an unacceptable design resulting from choosing the wrong ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$. For instance, consider a scenario where diagonal elements have a smaller magnitude than off-diagonal elements. This implies that ${\mathrm{M}}_2^{\mathrm{D}\mathrm{P}}$ and ${\mathrm{M}}_3^{\mathrm{D}\mathrm{P}}$, selected to control ${\mathrm{M}}_2^{\mathrm{F}\mathrm{R}}$ and ${\mathrm{M}}_3^{\mathrm{F}\mathrm{R}}$, respectively, have a greater impact on ${\mathrm{M}}_1^{\mathrm{F}\mathrm{R}}$ than ${\mathrm{M}}_1^{\mathrm{D}\mathrm{P}}$ has on ${\mathrm{M}}_1^{\mathrm{F}\mathrm{R}}$. In extreme cases, when one of the diagonal elements is significantly smaller than the off-diagonal elements, the design becomes coupled because one ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}$ can influence more than two ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}\mathrm{s}$. Additionally, if the number of ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$ is less than the number of ${\mathrm{M}}^{\mathrm{F}\mathrm{R}}\mathrm{s}$ due to one or more diagonal elements being zero, the design becomes coupled due to an insufficient number of ${\mathrm{M}}^{\mathrm{D}\mathrm{P}}\mathrm{s}$

4. Application

This study focuses on verifying the proposed methodology for analysing the robustness of manufacturing processes in the Iranian tire industry, specifically within 'Barez Industrial Group,' an enterprise based in Tehran, Iran, specializing in Tire Manufacturing since its establishment in 1984.

The manufacturing processes under investigation are the 'system of interest' (SOI), and their architecture is analysed in this study. Considering the dynamic nature of manufacturing environments, the study emphasizes the system's ability to cope with rapid changes. The existing architecture is examined from the perspective of robustness, aiming to identify any shortcomings or limitations. If significant weaknesses are detected, the study proposes a new architecture design based on the engineering design methodology presented in this research, ensuring the system's capability for robustness. By addressing the critical point of verifying the proposed methodology and emphasizing the importance of robustness in manufacturing processes, this study contributes valuable insights to the field of tire manufacturing and engineering design practices within 'Barez Industrial Group.'

4.1. Functional Architecture of the SOI

In step 2 of the research method, the functional architecture of the system of interest (SOI) is addressed and its constituent elements are studied. Applying the Axiomatic Design (AD) theory, the elements of the SOI within the functional domain of the system are identified, considering the elements as 'Functional Requirements' and tagging them with FRs. Exhibit 8 presents and describes these FRs of the SOI.

Exhibit 8. Description of the leaf functional requirements (FRs) of the SOI.

Code	Description	Code	Description
FR1111	Organize meetings for the research committee.	FR1139	Request the required tests and follow up on the test.
FR1112	Make a decision about the annual budget.	FR113,10	Prepare and archive APQP file documents.
FR1113	Define new projects for the offices.	FR113,11	Prepare documents and archive SQA documents.
FR1114	Decide on the duration of projects.	FR113,12	Check and prepare a product requirements chart, and turn it into a technical specification.
FR1115	Complete the project charter.	FR113,13	Perform the process of improvement and change.
FR1121	Conduct phenomenological studies of the project.	FR1141	Check test requests.
FR1122	Select the initial formulation.	FR1142	Coordinate to receive tires from the production line.
FR1123	Produce the experimental mixture.	FR1143	Plan to take the test.
FR1124	Evaluate the results.	FR1144	Run the test.
FR1125	Optimize the formulation.	FR1145	Review and analyze results.
FR1126	Produce a prototype.	FR1146	Generate test reports.
FR1127	Acknowledge the research objectives.	FR1151	Manage knowledge of the combined design process and the final product in all design aspects.
FR1128	Design the production process for in-line mixing.	FR1152	Support the knowledge management process to fix product line defects.
FR1129	Prepare mixed pre-production samples.	FR1153	Manage previous knowledge within the organization.
FR112,10	Prepare pre-production stage tires.	FR1154	Coordinate the holding of internal seminars.
FR112,11	Validate to achieve project goals.	FR1155	Coordinate and disseminate new sciences in the field of tire production.
FR112,12	Perform daily line production.	FR1156	Receive and archive lessons learned related to research and development projects.
FR112,13	Prepare the mixture reference recipe and deliver it to the Research and Development Documentation Center.	FR1157	Publish lessons learned at the organization level.
FR112,14	Prepare project documentation and store a copy in the archive of the mixture design office and deliver another copy to the documentation center.	FR1221	Select the size and brand of comparable tires.
FR1131	Design the template profile and pattern.	FR1222	Request to buy.
FR1132	Simulate and certify the design.	FR1223	Determine the test and its standard.
FR1133	Prepare template drawings.	FR1231	Issue processes based on design department ORDERS.
FR1134	Prepare the Common Technical Document (CTD).	FR1232	Check the production notifications sent by the process department to decide on issuing or not issuing the process.
FR1135	Build the ORDER design.	FR1233	Prepare PFMEA for all radial tire production processes.
FR1136	Issue the manufacturing process.	FR1234	Investigate and eliminate the reasons for the non-functioning of the radial tire in the test center tests.
FR1137	Issue the baking process.	FR1235	Carry out projects defined in meetings.
FR1138	Experimentally construct and cook, then examine defects and uniformity.	FR1236	Prepare and update work standards for different parts of radial riding.

Exhibit 8. Description of the leaf functional requirements (FRs) of the SOI.

Code	Description	Code	Description
FR1111	Organize meetings for the research committee.	FR1139	Request the required tests and follow up on the test.
FR1112	Make a decision about the annual budget.	FR113,10	Prepare and archive APQP file documents.
FR1113	Define new projects for the offices.	FR113,11	Prepare documents and archive SQA documents.
FR1114	Decide on the duration of projects.	FR113,12	Check and prepare a product requirements chart, and turn it into a technical specification.
FR1115	Complete the project charter.	FR113,13	Perform the process of improvement and change.
FR1121	Conduct phenomenological studies of the project.	FR1141	Check test requests.
FR1122	Select the initial formulation.	FR1142	Coordinate to receive tires from the production line.
FR1123	Produce the experimental mixture.	FR1143	Plan to take the test.
FR1124	Evaluate the results.	FR1144	Run the test.
FR1125	Optimize the formulation.	FR1145	Review and analyze results.
FR1126	Produce a prototype.	FR1146	Generate test reports.
FR1127	Acknowledge the research objectives.	FR1151	Manage knowledge of the combined design process and the final product in all design aspects.
FR1128	Design the production process for in-line mixing.	FR1152	Support the knowledge management process to fix product line defects.
FR1129	Prepare mixed pre-production samples.	FR1153	Manage previous knowledge within the organization.
FR112,10	Prepare pre-production stage tires.	FR1154	Coordinate the holding of internal seminars.
FR112,11	Validate to achieve project goals.	FR1155	Coordinate and disseminate new sciences in the field of tire production.
FR112,12	Perform daily line production.	FR1156	Receive and archive lessons learned related to research and development projects.
FR112,13	Prepare the mixture reference recipe and deliver it to the Research and Development Documentation Center.	FR1157	Publish lessons learned at the organization level.
FR112,14	Prepare project documentation and store a copy in the archive of the mixture design office and deliver another copy to the documentation center.	FR1221	Select the size and brand of comparable tires.
FR1131	Design the template profile and pattern.	FR1222	Request to buy.
FR1132	Simulate and certify the design.	FR1223	Determine the test and its standard.
FR1133	Prepare template drawings.	FR1231	Issue processes based on design department ORDERS.
FR1134	Prepare the Common Technical Document (CTD).	FR1232	Check the production notifications sent by the process department to decide on issuing or not issuing the process.
FR1135	Build the ORDER design.	FR1233	Prepare PFMEA for all radial tire production processes.
FR1136	Issue the manufacturing process.	FR1234	Investigate and eliminate the reasons for the non-functioning of the radial tire in the test center tests.
FR1137	Issue the baking process.	FR1235	Carry out projects defined in meetings.
FR1138	Experimentally construct and cook, then examine defects and uniformity.	FR1236	Prepare and update work standards for different parts of radial riding.

4.2. Originate the Functional Architecture of the SOI in DSM Representation.

Following steps 3-6 of the method, the DSM tool portrays the functional architecture of the system of interest (SOI). Exhibit 9 shows the original architecture design. Positive values above the diagonal in the original DSM signify substantial coupling, suggesting changes in one element affect others. To curtail change propagation, a beneficial solution called 'modularization' is suggested.

Exhibit 9. Displaying the initial SOI design using a Numerical Design Structure Matrix

Exhibit 9. Displaying the initial SOI design using a Numerical Design Structure Matrix

4.3. Significance of Each System Element within the Functional Architecture of the SOI

4.3.1. Initial Significance of Each System Element

Each SOI element's significance is assessed by 'C₁' (system architecture dependency) and 'C₂' (element's info dependency on architecture). 'C₁' gauges system's reliance on an element, while 'C₂' measures element's info contribution and input from the system.

For insights into element significance, percentage average dependencies are calculated from the original Numerical Design Structure Matrix (NDSM) of the architecture. Exhibits 10 and 11 display this, providing a comprehensive view of individual element importance and role in the SOI;

Exhibit 10. Percentage of information provided by each system element for the entire system (NDSM rows).

Code	Percentage (%)	Code	Percentage (%)	Code	Percentage (%)
FR1121	0.57	FR1231	1.36	FR1157	1.87
FR1115	0.6	FR1138	1.38	FR1152	1.9
FR1129	0.69	FR1235	1.38	FR112,11	1.96
FR112, 10	0.73	FR1113	1.44	FR1223	2.04
FR1124	0.83	FR1112	1.52	FR1136	2.25
FR1154	0.87	FR1142	1.55	FR1137	2.26
FR1123	0.88	FR1131	1.55	FR1146	2.44
FR1236	0.96	FR112,13	1.55	FR1134	2.49
FR1233	0.99	FR1232	1.56	FR113,13	2.64
FR1111	1.04	FR1221	1.56	FR1156	2.75
FR1234	1.11	FR112,12	1.59	FR1145	2.87
FR1155	1.12	FR1222	1.6	FR1144	3.11
FR1125	1.19	FR1133	1.64	FR1141	3.17
FR1139	1.25	FR1128	1.68	FR113,11	3.72
FR1132	1.27	FR1153	1.74	FR113,12	3.92
FR1122	1.29	FR1143	1.77	FR112,14	3.94
FR1127	1.31	FR1114	1.79	FR1151	4.34
FR1126	1.32	FR1135	1.82	FR113,10	5.85

Exhibit 10. Percentage of information provided by each system element for the entire system (NDSM rows).

Code	Percentage (%)	Code	Percentage (%)	Code	Percentage (%)
FR1121	0.57	FR1231	1.36	FR1157	1.87
FR1115	0.6	FR1138	1.38	FR1152	1.9
FR1129	0.69	FR1235	1.38	FR112,11	1.96
FR112, 10	0.73	FR1113	1.44	FR1223	2.04
FR1124	0.83	FR1112	1.52	FR1136	2.25
FR1154	0.87	FR1142	1.55	FR1137	2.26
FR1123	0.88	FR1131	1.55	FR1146	2.44
FR1236	0.96	FR112,13	1.55	FR1134	2.49
FR1233	0.99	FR1232	1.56	FR113,13	2.64
FR1111	1.04	FR1221	1.56	FR1156	2.75
FR1234	1.11	FR112,12	1.59	FR1145	2.87
FR1155	1.12	FR1222	1.6	FR1144	3.11
FR1125	1.19	FR1133	1.64	FR1141	3.17
FR1139	1.25	FR1128	1.68	FR113,11	3.72
FR1132	1.27	FR1153	1.74	FR113,12	3.92
FR1122	1.29	FR1143	1.77	FR112,14	3.94
FR1127	1.31	FR1114	1.79	FR1151	4.34
FR1126	1.32	FR1135	1.82	FR113,10	5.85

Exhibit 11. Percentage of information input received by each system element from the entire system (NDSM columns).

Code	Percentage (%)	Code	Percentage (%)	Code	Percentage (%)
FR1111	0.85	FR112,14	1.56	FR1145	2.13
FR1112	0.94	FR1131	1.77	FR1146	1.76
FR1113	1.13	FR1132	1.35	FR1151	1.82
FR1114	1.3	FR1133	1.88	FR1152	1.84
FR1115	0.63	FR1134	2.83	FR1153	1.72
FR1121	1.62	FR1135	2.18	FR1154	0.99
FR1122	2.08	FR1136	2.36	FR1155	1.42
FR1123	1.92	FR1137	2.22	FR1156	1.98
FR1124	1.94	FR1138	1.44	FR1157	1.97
FR1125	2.01	FR1139	1.28	FR1221	1.57
FR1126	1.79	FR113,10	3.53	FR1222	2.07
FR1127	1.64	FR113,11	3.83	FR1223	1.74
FR1128	1.9	FR113,12	4.13	FR1231	1.33
FR1129	1.68	FR113,13	2.93	FR1232	1.97
FR112,10	1.59	FR1141	2.41	FR1233	1.43
FR112,11	2.38	FR1142	1.04	FR1234	1.51
FR112,12	2.01	FR1143	1.8	FR1235	1.4
FR112,13	1.46	FR1144	2.51	FR1236	1.43

Exhibit 11. Percentage of information input received by each system element from the entire system (NDSM columns).

Code	Percentage (%)	Code	Percentage (%)	Code	Percentage (%)
FR1111	0.85	FR112,14	1.56	FR1145	2.13
FR1112	0.94	FR1131	1.77	FR1146	1.76
FR1113	1.13	FR1132	1.35	FR1151	1.82
FR1114	1.3	FR1133	1.88	FR1152	1.84
FR1115	0.63	FR1134	2.83	FR1153	1.72
FR1121	1.62	FR1135	2.18	FR1154	0.99
FR1122	2.08	FR1136	2.36	FR1155	1.42
FR1123	1.92	FR1137	2.22	FR1156	1.98
FR1124	1.94	FR1138	1.44	FR1157	1.97
FR1125	2.01	FR1139	1.28	FR1221	1.57
FR1126	1.79	FR113,10	3.53	FR1222	2.07
FR1127	1.64	FR113,11	3.83	FR1223	1.74
FR1128	1.9	FR113,12	4.13	FR1231	1.33
FR1129	1.68	FR113,13	2.93	FR1232	1.97
FR112,10	1.59	FR1141	2.41	FR1233	1.43
FR112,11	2.38	FR1142	1.04	FR1234	1.51
FR112,12	2.01	FR1143	1.8	FR1235	1.4
FR112,13	1.46	FR1144	2.51	FR1236	1.43

(NDSM)Exhibits 12 and 13 are included to provide a clearer visual representation of the information presented in Tables 10 and 11, respectively;

Exhibit 12. Prioritization of the proportion of information contribution from each system element to the entire SOI (NDSM rows).

Exhibit 12. Prioritization of the proportion of information contribution from each system element to the entire SOI (NDSM rows).

Exhibit 13. Prioritization of the percentage of information input received by each system element from the entire SOI (NDSM columns).

Exhibit 13. Prioritization of the percentage of information input received by each system element from the entire SOI (NDSM columns).

Furthermore, Exhibit 14 offers a comparison between the 'percentage of information feed that each system element can provide' and the 'percentage of information inputs that each system element receives' for/from the entire SOI;

Exhibit 14. Comparing the information contribution and input percentages for each system element within the entire system.

Exhibit 14. Comparing the information contribution and input percentages for each system element within the entire system.

As observed in the aforementioned figures, it is evident that the criteria C₁ and C₂ for each system element are not identical. For instance, while elements 'FR_113,10', 'FR_113,11', and 'FR_113,12' demonstrate maximal percentage of information feed they can provide for the entire system (i.e., C₁), they also receive the highest percentage of information inputs from the entire system (i.e., C₂). However, it is essential to avoid potential misinterpretation, as the significance of each process (i.e., each node of the network representing the system architecture of the SOI) should be evaluated based on both criteria concurrently.

To achieve this, the present study utilizes a weighting method known as the 'compound measure,' proposed within this research. The subsequent section will delve into further details about this approach.

4.3.2. Compound Weights of each system element

To establish a singular weight for each element, two criteria, namely 'C₁' and 'C₂' are introduced with equal values (i.e., WC₁ = 0.5, WC₂ = 0.5). Subsequently, two Multi-Criteria Decision Making (MCDM) weighting methods, 'SAW' and 'TOPSIS,' are utilized. The resulting weights obtained through these methods are collectively referred to as 'Compound Weights' in this study.

Exhibit 15 displays the final compound weights of each element within the SOI, based on the outcomes derived from the SAW and TOPSIS methods.

Exhibit 15. Compound weight of each element of the SOI.

	SAW	TOPSIS	FRs	SAW	TOPSIS
FR1111	0.009	0.005	FR1139	0.013	0.010
FR1112	0.012	0.010	FR113,10	0.045	0.059
FR1113	0.013	0.010	FR113,11	0.038	0.044
FR1114	0.015	0.014	FR113,12	0.040	0.047
FR1115	0.006	0.000	FR113,13	0.028	0.031
FR1121	0.012	0.010	FR1141	0.027	0.032
FR1122	0.018	0.016	FR1142	0.013	0.011
FR1123	0.015	0.013	FR1143	0.018	0.017
FR1124	0.015	0.013	FR1144	0.028	0.032
FR1125	0.017	0.015	FR1145	0.024	0.028
FR1126	0.016	0.014	FR1146	0.020	0.022
FR1127	0.015	0.013	FR1151	0.029	0.037
FR1128	0.018	0.017	FR1152	0.019	0.019
FR1129	0.013	0.011	FR1153	0.017	0.017
FR112,10	0.012	0.010	FR1154	0.009	0.005
FR112,11	0.022	0.023	FR1155	0.013	0.010
FR112,12	0.018	0.018	FR1156	0.023	0.026
FR112,13	0.015	0.013	FR1157	0.019	0.019
FR112,14	0.025	0.033	FR1221	0.016	0.014
FR1131	0.017	0.015	FR1222	0.019	0.018
FR1132	0.013	0.010	FR1223	0.019	0.019
FR1133	0.018	0.017	FR1231	0.013	0.011
FR1134	0.027	0.030	FR1232	0.018	0.017
FR1135	0.020	0.020	FR1233	0.013	0.009
FR1136	0.023	0.025	FR1234	0.013	0.011
FR1137	0.022	0.024	FR1235	0.014	0.012
FR1138	0.014	0.012	FR1236	0.012	0.009

Exhibit 15. Compound weight of each element of the SOI.

	SAW	TOPSIS	FRs	SAW	TOPSIS
FR1111	0.009	0.005	FR1139	0.013	0.010
FR1112	0.012	0.010	FR113,10	0.045	0.059
FR1113	0.013	0.010	FR113,11	0.038	0.044
FR1114	0.015	0.014	FR113,12	0.040	0.047
FR1115	0.006	0.000	FR113,13	0.028	0.031
FR1121	0.012	0.010	FR1141	0.027	0.032
FR1122	0.018	0.016	FR1142	0.013	0.011
FR1123	0.015	0.013	FR1143	0.018	0.017
FR1124	0.015	0.013	FR1144	0.028	0.032
FR1125	0.017	0.015	FR1145	0.024	0.028
FR1126	0.016	0.014	FR1146	0.020	0.022
FR1127	0.015	0.013	FR1151	0.029	0.037
FR1128	0.018	0.017	FR1152	0.019	0.019
FR1129	0.013	0.011	FR1153	0.017	0.017
FR112,10	0.012	0.010	FR1154	0.009	0.005
FR112,11	0.022	0.023	FR1155	0.013	0.010
FR112,12	0.018	0.018	FR1156	0.023	0.026
FR112,13	0.015	0.013	FR1157	0.019	0.019
FR112,14	0.025	0.033	FR1221	0.016	0.014
FR1131	0.017	0.015	FR1222	0.019	0.018
FR1132	0.013	0.010	FR1223	0.019	0.019
FR1133	0.018	0.017	FR1231	0.013	0.011
FR1134	0.027	0.030	FR1232	0.018	0.017
FR1135	0.020	0.020	FR1233	0.013	0.009
FR1136	0.023	0.025	FR1234	0.013	0.011
FR1137	0.022	0.024	FR1235	0.014	0.012
FR1138	0.014	0.012	FR1236	0.012	0.009

Exhibit 16 presents a comparison of the compound weight of each element obtained from both the TOPSIS method and the SAW method, alongside the weights originally computed (i.e., weights obtained from C₁ and C₂ only).

Exhibit 16. Comparing the compound weights of each element with their original individual weights.

Exhibit 16. Comparing the compound weights of each element with their original individual weights.

As the figure illustrates, there are considerable differences between the original weights and the weights obtained from applying the MCDM methods. However, the latter weights (i.e., weights given by the SAW and TOPSIS methods) are nearly similar, while also being more reasonable and acceptable.

Given the considerable similarity between the results obtained from the SAW and TOPSIS methods, for the sake of simplicity and without losing generality, the weights given by the TOPSIS method are used to continue the analyses. Exhibit 17 presents the updated prioritization of the elements based on the compound weights derived from the TOPSIS method;

Exhibit 17. Prioritization of elements according to their compound weights determined by the TOPSIS method.

System Element	Priority	Compound Weight	System Element	Priority	Compound Weight
FR113,10	29	0.0449	FR1131	20	0.0168
FR113,12	31	0.0404	FR1125	10	0.0167
FR113,11	30	0.0378	FR1126	11	0.0160
FR1151	39	0.0286	FR1221	46	0.0157
FR113,13	32	0.0281	FR1127	12	0.0151
FR1144	36	0.0276	FR1114	4	0.0151
FR1141	33	0.0272	FR112,13	18	0.0150
FR1134	23	0.0269	FR1123	8	0.0149
FR112,14	19	0.0255	FR1124	9	0.0148
FR1145	37	0.0243	FR1138	27	0.0141
FR1136	25	0.0231	FR1235	53	0.0139
FR1156	44	0.0230	FR1234	52	0.0135
FR1137	26	0.0224	FR1231	49	0.0134
FR112,11	16	0.0220	FR1132	21	0.0132
FR1146	38	0.0204	FR1155	43	0.0130
FR1135	24	0.0203	FR1129	14	0.0127
FR1157	45	0.0193	FR1139	28	0.0127
FR1222	47	0.0187	FR1113	3	0.0126
FR1152	40	0.0187	FR1233	51	0.0125
FR1223	48	0.0187	FR1142	34	0.0125
FR112,12	17	0.0184	FR1236	54	0.0124
FR1128	13	0.0181	FR112,10	15	0.0123
FR1232	50	0.0180	FR1121	6	0.0118
FR1143	35	0.0179	FR1112	2	0.0118
FR1133	22	0.0178	FR1154	42	0.0094
FR1122	7	0.0175	FR1111	1	0.0093
FR1153	41	0.0173	FR1115	5	0.0062

Exhibit 17. Prioritization of elements according to their compound weights determined by the TOPSIS method.

System Element	Priority	Compound Weight	System Element	Priority	Compound Weight
FR113,10	29	0.0449	FR1131	20	0.0168
FR113,12	31	0.0404	FR1125	10	0.0167
FR113,11	30	0.0378	FR1126	11	0.0160
FR1151	39	0.0286	FR1221	46	0.0157
FR113,13	32	0.0281	FR1127	12	0.0151
FR1144	36	0.0276	FR1114	4	0.0151
FR1141	33	0.0272	FR112,13	18	0.0150
FR1134	23	0.0269	FR1123	8	0.0149
FR112,14	19	0.0255	FR1124	9	0.0148
FR1145	37	0.0243	FR1138	27	0.0141
FR1136	25	0.0231	FR1235	53	0.0139
FR1156	44	0.0230	FR1234	52	0.0135
FR1137	26	0.0224	FR1231	49	0.0134
FR112,11	16	0.0220	FR1132	21	0.0132
FR1146	38	0.0204	FR1155	43	0.0130
FR1135	24	0.0203	FR1129	14	0.0127
FR1157	45	0.0193	FR1139	28	0.0127
FR1222	47	0.0187	FR1113	3	0.0126
FR1152	40	0.0187	FR1233	51	0.0125
FR1223	48	0.0187	FR1142	34	0.0125
FR112,12	17	0.0184	FR1236	54	0.0124
FR1128	13	0.0181	FR112,10	15	0.0123
FR1232	50	0.0180	FR1121	6	0.0118
FR1143	35	0.0179	FR1112	2	0.0118
FR1133	22	0.0178	FR1154	42	0.0094
FR1122	7	0.0175	FR1111	1	0.0093
FR1153	41	0.0173	FR1115	5	0.0062

Consequently, the updated prioritization of the elements can be of significant assistance to managers, enabling them to identify the most and least crucial elements with greater reliability. This, in turn, offers valuable practical insights for making decisions that prioritize robustness.

Exhibit 18. Displaying the initial SOI design in the NDSM with updated weights.

Exhibit 18. Displaying the initial SOI design in the NDSM with updated weights.

4.4. Modularization of the Functional Architecture of the SOI

As highlighted in the preceding section, reducing the coupling identified within the functional architecture of the SOI is crucial. According to the principles of DSM methodology, the existence of coupling is inferred from positive values of upper-diagonal entries in the original DSM. However, it is important to note that not all positive values necessarily indicate a significant strength between the involved elements. In certain cases, some interactions between two elements may be negligible, and as a result, the associated values for those entries can be replaced by zero instead. This approach allows us to focus on the more substantial interactions while mitigating the impact of insignificant ones.

Based on the aforementioned logic, elite members of the SOI were interviewed to determine a set of candidate cut-off points for ignoring the strength of interaction between a given pair of elements. The cut-off point refers to a specific value below which the interactions can be disregarded. As a result, the values '0.10', '0.11', '0.11', '0.12', '0.13', '0.14', and '0.15'—each reflecting the strength of interaction between system elements—are introduced as acceptable candidates for identifying an appropriate cut-off point.

For each candidate, the associated NDSM representation of the functional architecture of interest is presented in Appendix A (Exhibits A₁-A₆). To diminish the coupling, the partitioning technique is applied for each of the resulted architecture. Exhibits 19–24 display the outcomes.

Exhibit 19. Developing a modular design for the SOI architecture with a cutting-point set at 0.10 (Design I).

Exhibit 19. Developing a modular design for the SOI architecture with a cutting-point set at 0.10 (Design I).

Exhibit 20. Developing a modular design for the SOI architecture with a cutting-point set at 0.11 (Design II).

Exhibit 20. Developing a modular design for the SOI architecture with a cutting-point set at 0.11 (Design II).

Exhibit 21. Developing a modular design for the SOI architecture with a cutting-point set at 0.12 (Design III)

Exhibit 21. Developing a modular design for the SOI architecture with a cutting-point set at 0.12 (Design III)

Exhibit 22. Developing a modular design for the SOI architecture with a cutting-point set at 0.13 (Design IV).

Exhibit 22. Developing a modular design for the SOI architecture with a cutting-point set at 0.13 (Design IV).

Exhibit 23. Developing a modular design for the SOI architecture with a cutting-point set at 0.14 (Design V).

Exhibit 23. Developing a modular design for the SOI architecture with a cutting-point set at 0.14 (Design V).

Exhibit 24. Developing a modular design for the SOI architecture with a cutting-point set at 0.15 (Design VI).

Exhibit 24. Developing a modular design for the SOI architecture with a cutting-point set at 0.15 (Design VI).

Exhibit 25 presents all the fundamental information (i.e., N₁-N₄) necessary for computing the criteria, introduced in section 2.3.2, used to measure the architecture modularity;

Exhibit 25. Fundamental details regarding the introduced measures of modularity.

Ns Designs	N1	N2	N3	N4
Design I	131	499	2226	60
Design II	101	295	2469	51
Design III	86	144	2645	41
Design IV	80	146	2656	34
Design V	54	86	2742	34
Design VI	30	44	2810	32

Exhibit 25. Fundamental details regarding the introduced measures of modularity.

Ns Designs	N1	N2	N3	N4
Design I	131	499	2226	60
Design II	101	295	2469	51
Design III	86	144	2645	41
Design IV	80	146	2656	34
Design V	54	86	2742	34
Design VI	30	44	2810	32

Using the information provided in Exhibit 25, Exhibit 26 compares the modularity of all six presented designs (i.e., Designs I-VI). Based on the applied criteria, it is concluded that 'Design VI' for the functional architecture of the SOI exhibits the highest modularity. Therefore, in terms of modularity, this specific design is considered the most optimal for the architecture. Notably, 'Design VI' introduces a new sequence to enhance the architecture while effectively removing significant coupling;

Exhibit 26. Assessing the degree of modularity for Designs I through VI.

Presented Designs	Modularity Measures
		Grouping Capability Index (GCI)	Grouping Efficacy	Grouping Efficiency	Utilization Rate (U)
	Design I	0.68	-0.42	0.60	0.21
	Design II	0.66	-0.71	0.62	0.25
	Design III	0.68	-5.05	0.68	0.37
	Design IV	0.70	-2.50	0.67	0.35
	Design V	0.61	27.00	0.69	0.38
	Design VI	0.48	26.67	0.70	0.40

Exhibit 26. Assessing the degree of modularity for Designs I through VI.

Presented Designs	Modularity Measures
		Grouping Capability Index (GCI)	Grouping Efficacy	Grouping Efficiency	Utilization Rate (U)
	Design I	0.68	-0.42	0.60	0.21
	Design II	0.66	-0.71	0.62	0.25
	Design III	0.68	-5.05	0.68	0.37
	Design IV	0.70	-2.50	0.67	0.35
	Design V	0.61	27.00	0.69	0.38
	Design VI	0.48	26.67	0.70	0.40

Exhibits 27a-27d display the scores of each design concerning the modularity measures utilized in Exhibit 26.

Using the modularity measure of the 'modularization function,' the thirty-eight hierarchical sub-systems are compared based on the strength of coupling within the involved modules (Exhibit 28). The modules included in this analysis are referred to as the unique components of the modularization function.;

Based on Exhibit 28, the SOI in the best design for the physical architecture is optimally organized into thirty-eight hierarchical subsystems. The achieved modularity for the functional architecture of the SOI is considered acceptable.

Exhibit 28. Reorganizing the functional architecture into 38 hierarchical sub-systems, inspired by 'Design VI'

	Modules	Components	Number of Interfaces within Module	$\mathbf{Coupling}\mathbf{within}\mathbf{Module}({\mathsf{\delta }}_{{\mathbf{M}}_{\mathbf{i}}}$)	Number of Interfaces within Sub-System	$\mathbf{Coupling}\mathbf{within}\mathbf{Sub}-\mathbf{System}({\mathsf{\delta }}_{{\mathbf{M}\mathbf{S}\mathbf{S}}_1}$
Sub-System1	M1	FR1111	0	0	0	0
Sub-System2	M2	FR1112	13	13	0	0
		FR1122			0	0
		FR1121
		FR112,12
		FR1126
		FR1123
		FR1124
		FR1125
		FR113,10
		FR1128
		FR1127
		FR113,12
		FR113,11
Sub- System3	M3	FR1113	0	0	0	0
Sub- System4	M4	FR1114	0	0	0	0
Sub-System5	M5	FR1115	0	0	0	0
Sub- System6	M6	FR1129	0	0	0	0
Sub-System7	M7	FR112,10	0	0	0	0
Sub- System8	M8	FR112,11	0	0	0	0
Sub-System9	M9	FR112,13	0	0	0	0
Sub- System10	M10	FR112,14	0	0	0	0
Sub-System11	M11	FR1131	0	0	0	0
Sub- System12	M12	FR1132	0	0	0	0
Sub-System13	M13	FR1133	0	0	0	0
Sub- System14	M14	FR1134	1	1	0	0
		FR1136
Sub- System15	M15	FR1135	0	0	0	0
Sub-System16	M16	FR1137	0	0	0	0
Sub- System17	M17	FR1138	0	0	0	0
Sub-System18	M18	FR1139	0	0	0	0
Sub- System19	M19	FR113,13	0	0	0	0
Sub-System20	M20	FR1144	0	0	0	0
Sub- System21	M21	FR1145	0	0	0	0
Sub- System22	M22	FR1146	0	0	0	0
Sub-System23	M23	FR1151	0	0	0	0
Sub- System24	M24	FR1152	0	0	0	0
Sub-System25	M25	FR1153	0	0	0	0
Sub- System26	M26	FR1154	0	0	0	0
Sub-System27	M27	FR1155	0	0	0	0
Sub- System28	M28	FR1156	0	0	0	0
Sub-System29	M29	FR1157	0	0	0	0
Sub- System30	M30	FR1221	0	0	0	0
Sub-System31	M31	FR1222	0	0	0	0
Sub-System32	M32	FR1223	3	3	0	0
		FR1141
		FR1143
		FR1142
Sub-System33	M33	FR1231	0	0	0	0
Sub-System34	M34	FR1232	0	0	0	0
Sub-System35	M35	FR1233	0	0	0	0
Sub-System36	M36	FR1234	0	0	0	0
Sub-System37	M37	FR1235	0	0	0	0
Sub-System38	M38	FR1236	0	0	0	0

Exhibit 28. Reorganizing the functional architecture into 38 hierarchical sub-systems, inspired by 'Design VI'

	Modules	Components	Number of Interfaces within Module	$\mathbf{Coupling}\mathbf{within}\mathbf{Module}({\mathsf{\delta }}_{{\mathbf{M}}_{\mathbf{i}}}$)	Number of Interfaces within Sub-System	$\mathbf{Coupling}\mathbf{within}\mathbf{Sub}-\mathbf{System}({\mathsf{\delta }}_{{\mathbf{M}\mathbf{S}\mathbf{S}}_1}$
Sub-System1	M1	FR1111	0	0	0	0
Sub-System2	M2	FR1112	13	13	0	0
		FR1122			0	0
		FR1121
		FR112,12
		FR1126
		FR1123
		FR1124
		FR1125
		FR113,10
		FR1128
		FR1127
		FR113,12
		FR113,11
Sub- System3	M3	FR1113	0	0	0	0
Sub- System4	M4	FR1114	0	0	0	0
Sub-System5	M5	FR1115	0	0	0	0
Sub- System6	M6	FR1129	0	0	0	0
Sub-System7	M7	FR112,10	0	0	0	0
Sub- System8	M8	FR112,11	0	0	0	0
Sub-System9	M9	FR112,13	0	0	0	0
Sub- System10	M10	FR112,14	0	0	0	0
Sub-System11	M11	FR1131	0	0	0	0
Sub- System12	M12	FR1132	0	0	0	0
Sub-System13	M13	FR1133	0	0	0	0
Sub- System14	M14	FR1134	1	1	0	0
		FR1136
Sub- System15	M15	FR1135	0	0	0	0
Sub-System16	M16	FR1137	0	0	0	0
Sub- System17	M17	FR1138	0	0	0	0
Sub-System18	M18	FR1139	0	0	0	0
Sub- System19	M19	FR113,13	0	0	0	0
Sub-System20	M20	FR1144	0	0	0	0
Sub- System21	M21	FR1145	0	0	0	0
Sub- System22	M22	FR1146	0	0	0	0
Sub-System23	M23	FR1151	0	0	0	0
Sub- System24	M24	FR1152	0	0	0	0
Sub-System25	M25	FR1153	0	0	0	0
Sub- System26	M26	FR1154	0	0	0	0
Sub-System27	M27	FR1155	0	0	0	0
Sub- System28	M28	FR1156	0	0	0	0
Sub-System29	M29	FR1157	0	0	0	0
Sub- System30	M30	FR1221	0	0	0	0
Sub-System31	M31	FR1222	0	0	0	0
Sub-System32	M32	FR1223	3	3	0	0
		FR1141
		FR1143
		FR1142
Sub-System33	M33	FR1231	0	0	0	0
Sub-System34	M34	FR1232	0	0	0	0
Sub-System35	M35	FR1233	0	0	0	0
Sub-System36	M36	FR1234	0	0	0	0
Sub-System37	M37	FR1235	0	0	0	0
Sub-System38	M38	FR1236	0	0	0	0

Regarding the optimal design VI, it is crucial to mention that among all developed subsystems (i.e., subsystems 1-38), subsystem₂, subsystem₃₂, and subsystem₁₄ are particularly susceptible to any adjustments due to the presence of coupled elements. For these subsystems, additional techniques like the 'tearing technique' may prove more helpful if all the necessary information for making required assumptions is available.

4.5. Optimal Modularized Design Properties for SOI Functional Architecture and Its Robust Contribution

Through the development of the optimal modular design for the functional architecture of the SOI, the following properties have been achieved:

• The information flow in the optimal design for the functional architecture follows a top-down approach, where elements placed at a higher priority require little to no information from other elements, while elements at a lower priority may require more information.
• The information flow in the optimal design is optimally organized into thirty-eight hierarchical subsystems, ensuring that all necessary information to fulfill the architecture and adapt to new situations is readily available.
• The information flow in the optimal design is now clearer, simpler, more agile, and more traceable.
• The achieved modularity in the architecture provides numerous possibilities for robustness. Among the presented subsystems (i.e., subsystems 1-38), many places can be easily identified and adjusted to enhance the architecture's robustness in response to new situations.
• The optimal design highlights places in the architecture where modifications are challenging due to the presence of coupled elements. These places are susceptible to any change or modification, and as a result, they remain unchanged in all possible architectures achieved through adaptation operations. These places are referred to as common platforms.

4.6. Investigating Functional Architecture Robustness and Optimal Design Selection

According to steps 7-10 of the proposed methodology, a simulation study aimed at assessing the robustness of the modularized architecture is conducted, encompassing Designs I to VI. This analysis entails evaluating the capacity of each design to withstand varying levels of complexity, which is indicative of the information content that can be accommodated by each architectural element. The outcomes of this evaluation are outlined in Exhibit 29 and Exhibit 30.

Exhibit 29 provides an in-depth examination of the functional integration within the System of Interest (SOI), yielding insights into the cohesive functioning of its constituent components. This analysis facilitates the comprehension of the restorative or adaptive effort necessitated in response to damage or modifications. This observation holds true for both the inherent design and the novel modular configurations embraced by the SOI. Of particular interest is the comparative evaluation illuminated by Exhibit 30. This comparison underscores a discernible dissimilarity between the exigencies of the original design and the alternative modular constructs. A quantitative assessment, quantified in temporal and resource-intensive terms (Bits), reveals a pronounced differential in the demands imposed by the original versus modular design paradigms. This divergence is further elucidated in Exhibit 30, a visual representation that vividly accentuates these variances. Notably, Design VI emerges as a salient exemplar within this analysis, distinctly showcasing its heightened robustness relative to the original design. Exhibit 30 supplements this discernment by delineating the correlation between system resilience and complexity. The graph emphasizes a positive link between simple design and effective problem-solving and robustness. In short, avoiding excessive complexity boosts the system's ability to adeptly handle challenges.

Exhibit 29. Complexity (work volume) of each SOI element under perturbation (Designs I-VI).

Elements of System Architecture	Design I	Design II	Design III	Design IV	Design V	Design VI
FR1111	6.49	4.57	2.52	2.43	1.92	1.19
FR1112	5.18	4.74	2.99	2.01	1.98	1.19
FR1113	5.33	3.38	3.21	2.36	1.93	1.24
FR1114	6.38	4.87	3.11	2.41	1.73	1.19
FR1115	6.38	4.44	3.02	2.31	1.83	1.18
FR1121	9.79	4.54	2.67	2.18	1.89	1.24
FR1122	9.77	3.45	2.83	2.17	1.87	1.24
FR1123	5.54	3.8	2.77	2.01	1.77	1.24
FR1124	7.78	3.42	2.57	2.03	1.68	1.2
FR1125	7.74	3.81	2.85	2.33	1.78	1.2
FR1126	5.25	4.42	3.04	2.09	1.7	1.22
FR1127	8.56	3.53	2.73	2.03	1.83	1.23
FR1128	6.65	4.19	3.1	2.09	1.9	1.18
FR1129	5.3	4.58	3.25	2.01	1.98	1.18
FR112,10	5.59	3.96	2.52	2.33	1.84	1.2
FR112,11	5.99	3.6	3.15	2.37	1.72	1.21
FR112,12	8.94	4.54	2.77	2.24	1.74	1.19
FR112,13	6.48	3.68	2.8	2.36	1.81	1.24
FR112,14	5.26	4.84	2.55	2.21	1.71	1.22
FR1131	5.36	4.06	2.95	2.3	1.87	1.19
FR1132	5.41	4.47	3.13	2.11	1.7	1.18
FR1133	5.31	3.47	2.9	2.13	1.71	1.2
FR1134	9.04	3.66	2.71	2	1.81	1.19
FR1135	6.38	4.29	2.88	2.11	1.93	1.22
FR1136	6.35	3.44	2.97	2.04	1.83	1.21
FR1137	6.61	3.41	2.75	2.19	1.68	1.19
FR1138	6.99	3.53	2.77	2.17	1.75	1.2
FR1139	5.25	4.09	2.73	2.2	1.75	1.24
FR113,10	7.47	3.54	3.32	2.41	1.91	1.23
FR113,11	6.17	3.64	2.98	2.09	1.94	1.2
FR113,12	5.93	4.9	2.92	2.15	1.95	1.19
FR113,13	7.98	3.97	2.52	2.03	1.88	1.22
FR1141	8.94	4.26	3.3	2.46	1.92	1.24
FR1142	5.85	3.51	3.03	2.27	1.75	1.18
FR1143	6.84	4.2	2.67	2.28	1.95	1.19
FR1144	7.13	3.44	2.51	2.06	1.95	1.24
FR1145	5.73	3.34	2.62	2.31	1.81	1.23
FR1146	9.2	3.77	2.83	2.06	1.77	1.23
FR1151	8.47	4.28	2.75	2.1	1.89	1.18
FR1152	6.91	3.94	2.72	2.17	1.73	1.23
FR1153	8.97	3.72	2.75	2.03	1.71	1.23
FR1154	9.1	3.82	3.08	2.26	1.81	1.22
FR1155	6.38	4.36	3.06	2.31	1.89	1.21
FR1156	7.33	4.16	3.07	2.12	1.73	1.22
FR1157	7.75	4.8	2.69	2.39	1.89	1.23
FR1221	6.94	4.15	3.16	2.4	1.88	1.18
FR1222	7.89	4.16	2.77	2.34	1.73	1.2
FR1223	9.29	3.58	2.97	2.26	1.96	1.18
FR1231	6.94	4.66	2.51	2.29	1.85	1.2
FR1232	7.02	4.97	3.06	2.48	1.97	1.21
FR1233	9.98	4.2	3.02	2.36	1.8	1.25
FR1234	5.53	3.74	2.6	2.14	1.91	1.23
FR1235	8.84	4.69	3.09	2.17	1.97	1.22
FR1236	6	3.7	3.32	2.24	1.77	1.2

Exhibit 29. Complexity (work volume) of each SOI element under perturbation (Designs I-VI).

Elements of System Architecture	Design I	Design II	Design III	Design IV	Design V	Design VI
FR1111	6.49	4.57	2.52	2.43	1.92	1.19
FR1112	5.18	4.74	2.99	2.01	1.98	1.19
FR1113	5.33	3.38	3.21	2.36	1.93	1.24
FR1114	6.38	4.87	3.11	2.41	1.73	1.19
FR1115	6.38	4.44	3.02	2.31	1.83	1.18
FR1121	9.79	4.54	2.67	2.18	1.89	1.24
FR1122	9.77	3.45	2.83	2.17	1.87	1.24
FR1123	5.54	3.8	2.77	2.01	1.77	1.24
FR1124	7.78	3.42	2.57	2.03	1.68	1.2
FR1125	7.74	3.81	2.85	2.33	1.78	1.2
FR1126	5.25	4.42	3.04	2.09	1.7	1.22
FR1127	8.56	3.53	2.73	2.03	1.83	1.23
FR1128	6.65	4.19	3.1	2.09	1.9	1.18
FR1129	5.3	4.58	3.25	2.01	1.98	1.18
FR112,10	5.59	3.96	2.52	2.33	1.84	1.2
FR112,11	5.99	3.6	3.15	2.37	1.72	1.21
FR112,12	8.94	4.54	2.77	2.24	1.74	1.19
FR112,13	6.48	3.68	2.8	2.36	1.81	1.24
FR112,14	5.26	4.84	2.55	2.21	1.71	1.22
FR1131	5.36	4.06	2.95	2.3	1.87	1.19
FR1132	5.41	4.47	3.13	2.11	1.7	1.18
FR1133	5.31	3.47	2.9	2.13	1.71	1.2
FR1134	9.04	3.66	2.71	2	1.81	1.19
FR1135	6.38	4.29	2.88	2.11	1.93	1.22
FR1136	6.35	3.44	2.97	2.04	1.83	1.21
FR1137	6.61	3.41	2.75	2.19	1.68	1.19
FR1138	6.99	3.53	2.77	2.17	1.75	1.2
FR1139	5.25	4.09	2.73	2.2	1.75	1.24
FR113,10	7.47	3.54	3.32	2.41	1.91	1.23
FR113,11	6.17	3.64	2.98	2.09	1.94	1.2
FR113,12	5.93	4.9	2.92	2.15	1.95	1.19
FR113,13	7.98	3.97	2.52	2.03	1.88	1.22
FR1141	8.94	4.26	3.3	2.46	1.92	1.24
FR1142	5.85	3.51	3.03	2.27	1.75	1.18
FR1143	6.84	4.2	2.67	2.28	1.95	1.19
FR1144	7.13	3.44	2.51	2.06	1.95	1.24
FR1145	5.73	3.34	2.62	2.31	1.81	1.23
FR1146	9.2	3.77	2.83	2.06	1.77	1.23
FR1151	8.47	4.28	2.75	2.1	1.89	1.18
FR1152	6.91	3.94	2.72	2.17	1.73	1.23
FR1153	8.97	3.72	2.75	2.03	1.71	1.23
FR1154	9.1	3.82	3.08	2.26	1.81	1.22
FR1155	6.38	4.36	3.06	2.31	1.89	1.21
FR1156	7.33	4.16	3.07	2.12	1.73	1.22
FR1157	7.75	4.8	2.69	2.39	1.89	1.23
FR1221	6.94	4.15	3.16	2.4	1.88	1.18
FR1222	7.89	4.16	2.77	2.34	1.73	1.2
FR1223	9.29	3.58	2.97	2.26	1.96	1.18
FR1231	6.94	4.66	2.51	2.29	1.85	1.2
FR1232	7.02	4.97	3.06	2.48	1.97	1.21
FR1233	9.98	4.2	3.02	2.36	1.8	1.25
FR1234	5.53	3.74	2.6	2.14	1.91	1.23
FR1235	8.84	4.69	3.09	2.17	1.97	1.22
FR1236	6	3.7	3.32	2.24	1.77	1.2

Upon close examination of Exhibit 30, Design VI emerges as the optimal candidate among the array of six distinct modular designs. Notably, Design VI exhibits a pronounced proclivity for effectively navigating challenging scenarios. This proclivity underscores its pre-eminence within the spectrum of designs under scrutiny, particularly in comparison to the original design.

Exhibit 30. Comparison of robustness factor of each system element among the Designs I–VI.

Exhibit 30. Comparison of robustness factor of each system element among the Designs I–VI.

Analogous to a formidable sentinel within the realm of design configurations, Design VI assumes the role of an architectural vanguard. Its singular aptitude lies in the augmentation of the System of Interest's (SOI) structural robustness and stability. Evidently, this revelation signifies Design VI's unique potential to substantially elevate the quality of the SOI's system architecture, thereby endowing it with a heightened capacity to navigate exigent circumstances.

Exhibit 30 provides empirical support for this observation, wherein the comparative assessment of robustness operations across all six modular designs showcases Design VI as being intrinsically associated with the most modest complexities. This distinct trait bestows upon modularized Design VI a superior capability in effectuating robustness operations within the broader system architecture of the SOI. This assertion holds particularly true when juxtaposed against the original SOI design, further validating Design VI's exceptional competence in orchestrating structural fortitude amidst deviations from anticipated circumstances.

5. Discussion and Conclusion

In this study, we addressed the challenges faced by manufacturing processes in dynamic environments characterized by rapid changes, uncertainties, risks, and uncontrollable fluctuations. Our objective was to enhance the overall robustness of manufacturing processes through an innovative approach that integrates engineering management principles and modular design. By adopting a systems engineering perspective, we treated all manufacturing process activities as interrelated elements of an integrated system, aiming to prevent the spread of risks and disruptions among various processes.

Our research revealed a profound relationship between manufacturing management and engineering management. Manufacturing management plays a pivotal role in overseeing resource allocation, decision-making processes, and the overall organizational strategy for effective operations. Inefficient manufacturing management can lead to suboptimal use of resources, delays in production schedules, and reduced productivity, negatively impacting engineering management. On the other hand, effective manufacturing management systems can optimize resource allocation, streamline decision-making processes, and result in improved production processes and increased efficiency in utilizing technical personnel and physical resources.

Engineering management, in turn, is responsible for overseeing the design, development, and improvement of manufacturing processes. The effectiveness of engineering management systems, including quality control and continuous improvement programs, relies on the underlying architecture of manufacturing. A well-designed manufacturing process architecture is crucial in minimizing the likelihood of operational risks and their potential impact on manufacturing processes.

Our findings showcased a significant improvement in modularization, leading to enhanced overall robustness in manufacturing processes. This highlights the importance of engineering management's proactive involvement in the design phase of manufacturing processes. By engaging manufacturing experts from the outset, potential production and assembly challenges can be addressed proactively, resulting in an optimized and robust manufacturing process.

Furthermore, our research emphasizes the importance of modularity within the manufacturing process architecture. Considering modularity in the process domain architecture indirectly supports modularity in the allocated architecture, fostering a comprehensive approach that optimizes component allocation within the physical system.

Additionally, we highlighted the significance of redundancy in certain system elements, as it enhances system reliability while preserving modularity. Engineering management's consideration of redundancy offers another avenue to enhance the overall robustness of manufacturing processes.

The successful integration of engineering management principles with modular design strategies introduces new avenues for enhancing robustness and efficiency in manufacturing processes. By proactively engaging engineering management in the design phase and applying our modularization approach, manufacturing organizations can enhance their overall system robustness and improve their capacity to respond to changing circumstances effectively.

Looking ahead, our research calls for further investigations into the interplay between manufacturing management and engineering management. Understanding the dynamics of this relationship will allow organizations to optimize resource allocation and streamline decision-making processes while fostering a culture of collaboration and knowledge exchange between engineering and manufacturing teams.

Moreover, our study opens avenues for future research in the development of comprehensive strategies that address other essential dimensions of manufacturing processes. Studies focused on process legality, reliability, adaptability, and recoverability will contribute to a more comprehensive understanding of resilient manufacturing processes.

Develop a resilience framework for manufacturing processes that encompasses not only robustness but also other dimensions such as flexibility, adaptability, and recoverability. This framework would provide a holistic view of manufacturing process resilience and guide decision-makers in making informed choices to enhance overall system resilience.