Pub Date : 1900-01-01DOI: 10.23967/wccm-apcom.2022.092
E. Tsuchida
. Computational materials design is an active area of research which aims at predicting phys-ical and chemical properties of various materials from first-principles electronic structure calculations. To keep the computational costs manageable, the Schr¨odinger equations are often approximated by Kohn-Sham equations within the framework of density-functional theory. These Kohn-Sham equations are solved numerically either by a basis set expansion or real-space discretization under given boundary conditions. In the case of a plane-wave basis set, it is common practice to apply periodic boundary conditions in all directions, while isolated boundary conditions are more common for the atomic basis set. However, there are many other options besides these standard boundary conditions. In this pre-sentation, we will explore several non-standard boundary conditions which exploit the characteristics of each system, such as surfaces, interfaces, and cyclic/helical structures, to minimize the computational costs of electronic structure calculations. Most of these boundary conditions are easily implemented by minor modifications of existing electronic structure codes. Numerical examples on a few model systems are also presented for the validation of these boundary conditions.
{"title":"Practical Boundary Conditions for Electronic Structure Calculations","authors":"E. Tsuchida","doi":"10.23967/wccm-apcom.2022.092","DOIUrl":"https://doi.org/10.23967/wccm-apcom.2022.092","url":null,"abstract":". Computational materials design is an active area of research which aims at predicting phys-ical and chemical properties of various materials from first-principles electronic structure calculations. To keep the computational costs manageable, the Schr¨odinger equations are often approximated by Kohn-Sham equations within the framework of density-functional theory. These Kohn-Sham equations are solved numerically either by a basis set expansion or real-space discretization under given boundary conditions. In the case of a plane-wave basis set, it is common practice to apply periodic boundary conditions in all directions, while isolated boundary conditions are more common for the atomic basis set. However, there are many other options besides these standard boundary conditions. In this pre-sentation, we will explore several non-standard boundary conditions which exploit the characteristics of each system, such as surfaces, interfaces, and cyclic/helical structures, to minimize the computational costs of electronic structure calculations. Most of these boundary conditions are easily implemented by minor modifications of existing electronic structure codes. Numerical examples on a few model systems are also presented for the validation of these boundary conditions.","PeriodicalId":429847,"journal":{"name":"15th World Congress on Computational Mechanics (WCCM-XV) and 8th Asian Pacific Congress on Computational Mechanics (APCOM-VIII)","volume":"34 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"1900-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"114774717","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 1900-01-01DOI: 10.23967/wccm-apcom.2022.120
G. Galik, V. Kutǐs, J. Paulech, V. Goga, M. Uličný
{"title":"Transient cooling of reactor vessel wall during LOCA","authors":"G. Galik, V. Kutǐs, J. Paulech, V. Goga, M. Uličný","doi":"10.23967/wccm-apcom.2022.120","DOIUrl":"https://doi.org/10.23967/wccm-apcom.2022.120","url":null,"abstract":"","PeriodicalId":429847,"journal":{"name":"15th World Congress on Computational Mechanics (WCCM-XV) and 8th Asian Pacific Congress on Computational Mechanics (APCOM-VIII)","volume":"90 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"1900-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"122036159","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 1900-01-01DOI: 10.23967/wccm-apcom.2022.045
M. Ichimiya
. Gas entrainment is one of the major defects caused in the casting filling process. Since the particle method is a Lagrangian method that does not use a lattice, it can easily analyze large deformations and boundary movements, so it has the potential to be applied to gas defect prediction as a methodology. However, in the two-phase flow simulation including gas entrainment, the analysis fails when the gas / liquid density ratio becomes smaller than about 1/10 in the conventional SPH. Therefore, two dimensional two-phase flow SPH methodology was developed. Then, authors extended the methodology to three dimensions that can be applied to gas entrainent during a die cast filling process.
{"title":"Development of two-phase flow simulation using SPH Method","authors":"M. Ichimiya","doi":"10.23967/wccm-apcom.2022.045","DOIUrl":"https://doi.org/10.23967/wccm-apcom.2022.045","url":null,"abstract":". Gas entrainment is one of the major defects caused in the casting filling process. Since the particle method is a Lagrangian method that does not use a lattice, it can easily analyze large deformations and boundary movements, so it has the potential to be applied to gas defect prediction as a methodology. However, in the two-phase flow simulation including gas entrainment, the analysis fails when the gas / liquid density ratio becomes smaller than about 1/10 in the conventional SPH. Therefore, two dimensional two-phase flow SPH methodology was developed. Then, authors extended the methodology to three dimensions that can be applied to gas entrainent during a die cast filling process.","PeriodicalId":429847,"journal":{"name":"15th World Congress on Computational Mechanics (WCCM-XV) and 8th Asian Pacific Congress on Computational Mechanics (APCOM-VIII)","volume":"30 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"1900-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"124680125","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 1900-01-01DOI: 10.23967/wccm-apcom.2022.021
M. Negi, M. Mukherjee
. While characterizing the mechanical behavior of granular assemblies through DEM simulations, various macro and micro level heterogeneities are often encountered. Such macro level heterogeneities may arise due to stress and void concentration near the wall boundaries; whereas, the micro level heterogeneities are attributed to consideration of only limited number of particles within the representative volume element (RVE). The present study assesses these macro and micro level heterogeneities in reference to the mechanical characterization of sand in DEM-based biaxial test simulation with both rigid and flexible lateral boundaries. In this regard, stresses and strains have been calculated using a wall-based global estimation and a representative area element (RAE)-based local estimation. It has been suggested that the RAE should occupy a maximum of 90% area of the specimen in order to avoid any macro level heterogeneity and can still be able to capture its overall mechanical behavior. For obtaining the spatial variation of field variables, RAE of smaller diameters are often employed. In such cases, depending on the average particle size of the granular assembly and the specimen dimensions, the diameter of the RAE should be selected ensuring that it is small enough to aptly capture the local variation of field variables and at the same time, large enough to avoid any micro level heterogeneity.
{"title":"Assessment of Macro and Micro level Heterogeneities for Characterizing Mechanical Behavior of Sand in Biaxial Test employing DEM","authors":"M. Negi, M. Mukherjee","doi":"10.23967/wccm-apcom.2022.021","DOIUrl":"https://doi.org/10.23967/wccm-apcom.2022.021","url":null,"abstract":". While characterizing the mechanical behavior of granular assemblies through DEM simulations, various macro and micro level heterogeneities are often encountered. Such macro level heterogeneities may arise due to stress and void concentration near the wall boundaries; whereas, the micro level heterogeneities are attributed to consideration of only limited number of particles within the representative volume element (RVE). The present study assesses these macro and micro level heterogeneities in reference to the mechanical characterization of sand in DEM-based biaxial test simulation with both rigid and flexible lateral boundaries. In this regard, stresses and strains have been calculated using a wall-based global estimation and a representative area element (RAE)-based local estimation. It has been suggested that the RAE should occupy a maximum of 90% area of the specimen in order to avoid any macro level heterogeneity and can still be able to capture its overall mechanical behavior. For obtaining the spatial variation of field variables, RAE of smaller diameters are often employed. In such cases, depending on the average particle size of the granular assembly and the specimen dimensions, the diameter of the RAE should be selected ensuring that it is small enough to aptly capture the local variation of field variables and at the same time, large enough to avoid any micro level heterogeneity.","PeriodicalId":429847,"journal":{"name":"15th World Congress on Computational Mechanics (WCCM-XV) and 8th Asian Pacific Congress on Computational Mechanics (APCOM-VIII)","volume":"91 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"1900-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"114974542","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 1900-01-01DOI: 10.23967/wccm-apcom.2022.025
S. Yamamoto, G. Shoji, M. Ohsumi
. This study aims to clarify the damage mechanism of a long-period bridge system — the Ohkirihata Bridge damaged in the 2016 Kumamoto earthquake — subjected to the combined effects of long-period pulsive ground motions and surface fault displacements. The target bridge’s site-specific waveforms at abutment A1 were estimated using the finite difference method. Linear dynamic analysis with a three-dimensional finite element model of the bridge structure-underground interconnected system was performed to examine the effects of long-period pulsive ground motions on the coupled responses of essential structural components: superstructure, rubber bearings, abutments, piers, foundations and underground.
{"title":"Clarification of the Damage Mechanism of the Long-Period Bridge System Damaged by the 2016 Kumamoto Earthquake","authors":"S. Yamamoto, G. Shoji, M. Ohsumi","doi":"10.23967/wccm-apcom.2022.025","DOIUrl":"https://doi.org/10.23967/wccm-apcom.2022.025","url":null,"abstract":". This study aims to clarify the damage mechanism of a long-period bridge system — the Ohkirihata Bridge damaged in the 2016 Kumamoto earthquake — subjected to the combined effects of long-period pulsive ground motions and surface fault displacements. The target bridge’s site-specific waveforms at abutment A1 were estimated using the finite difference method. Linear dynamic analysis with a three-dimensional finite element model of the bridge structure-underground interconnected system was performed to examine the effects of long-period pulsive ground motions on the coupled responses of essential structural components: superstructure, rubber bearings, abutments, piers, foundations and underground.","PeriodicalId":429847,"journal":{"name":"15th World Congress on Computational Mechanics (WCCM-XV) and 8th Asian Pacific Congress on Computational Mechanics (APCOM-VIII)","volume":"63 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"1900-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"130272291","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 1900-01-01DOI: 10.23967/wccm-apcom.2022.086
S. Wani, R. Samala, R. Kandasami, A. Chaudhuri
. Exploratory studies have been carried out to identify the potential natural gas hydrate reserves for commercially producing gas. While extracting the gas from the hydrate-bearing sediments using various dissociation techniques, there will be a significant loss of strength in these sediments. It is well known that the behavior of gas hydrate sediments is governed by Thermo Hydro Mechanical Chemical – THMC coupled process during the gas extraction. Thus, in this study, in order to understand the influence of depressurization at the well-bore and the permeability of the hydrate reservoir on the sediment deformation characteristics, a 2D (plane strain condition) hydrate reservoir is simulated (using a multiphase numerical schema). From the study, it is observed that the flow response, i.e., the rate of change of gas pressure near the well-bore, decreases with the increase in the duration of the extraction. The maximum settlement occurs for reservoirs having low well-bore pressure (higher amount of depressurization) and high intrinsic permeability. Additionally, these same reservoir conditions also lead to maximum cumulative gas production. Thus, the continuous gas
{"title":"Numerical study on the hydrate-rich sediment behaviour during depressurization","authors":"S. Wani, R. Samala, R. Kandasami, A. Chaudhuri","doi":"10.23967/wccm-apcom.2022.086","DOIUrl":"https://doi.org/10.23967/wccm-apcom.2022.086","url":null,"abstract":". Exploratory studies have been carried out to identify the potential natural gas hydrate reserves for commercially producing gas. While extracting the gas from the hydrate-bearing sediments using various dissociation techniques, there will be a significant loss of strength in these sediments. It is well known that the behavior of gas hydrate sediments is governed by Thermo Hydro Mechanical Chemical – THMC coupled process during the gas extraction. Thus, in this study, in order to understand the influence of depressurization at the well-bore and the permeability of the hydrate reservoir on the sediment deformation characteristics, a 2D (plane strain condition) hydrate reservoir is simulated (using a multiphase numerical schema). From the study, it is observed that the flow response, i.e., the rate of change of gas pressure near the well-bore, decreases with the increase in the duration of the extraction. The maximum settlement occurs for reservoirs having low well-bore pressure (higher amount of depressurization) and high intrinsic permeability. Additionally, these same reservoir conditions also lead to maximum cumulative gas production. Thus, the continuous gas","PeriodicalId":429847,"journal":{"name":"15th World Congress on Computational Mechanics (WCCM-XV) and 8th Asian Pacific Congress on Computational Mechanics (APCOM-VIII)","volume":"18 3","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"1900-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"114111636","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 1900-01-01DOI: 10.23967/wccm-apcom.2022.013
T. Asaka, K. Yoshimatsu, K. Schneider
. We develop a wavelet-based three-dimensional convolutional neural network (WCNN3d) for superresolution of coarse-grained data of homogeneous isotropic turbulence. The turbulent flow data are computed by high resolution direct numerical simulation (DNS), while the coarse-grained data are obtained by applying a Gaussian filter to the DNS data. The CNNs are trained with the DNS data and the coarse-grained data. We compare vorticity- and velocity-based approaches and assess the proposed WCNN3d method in terms of flow visualization, enstrophy spectra and probability density functions. We show that orthogonal wavelets enhance the efficiency of the learning of CNN. of isotropic turbulence in a periodic box and the coarse-grained data are obtained by the application of a Gaussian low-pass filter to the DNS data. We assessed the WCNN3d in terms of 3D visualization of vorticity, PDF of vorticity, and enstrophy spectra. We found that WCNN3d well reproduces vorticity statistics and the positions of the vortices from coarse-grained vorticity fields. For the vorticity-based approach, the use of wavelets enhances deep learning of turbulent flows considered here thanks to the sparsity of the wavelet representation which improves deep learning. For the velocity-based approach, we showed that weighting the wavelet coefficients of velocity, which yields velocity gradient information due to norm equivalence, improves the accuracy and yields results similar to the vorticity-based model. Furthermore, we assessed the divergence issue of the predicted fields and showed that its impact is negligible. We demonstrated the capability to predict a turbulent flow whose Reynolds number is higher than the flows used for the training.
{"title":"A wavelet-based three-dimensional Convolutional Neural Network for superresolution of turbulent vorticity","authors":"T. Asaka, K. Yoshimatsu, K. Schneider","doi":"10.23967/wccm-apcom.2022.013","DOIUrl":"https://doi.org/10.23967/wccm-apcom.2022.013","url":null,"abstract":". We develop a wavelet-based three-dimensional convolutional neural network (WCNN3d) for superresolution of coarse-grained data of homogeneous isotropic turbulence. The turbulent flow data are computed by high resolution direct numerical simulation (DNS), while the coarse-grained data are obtained by applying a Gaussian filter to the DNS data. The CNNs are trained with the DNS data and the coarse-grained data. We compare vorticity- and velocity-based approaches and assess the proposed WCNN3d method in terms of flow visualization, enstrophy spectra and probability density functions. We show that orthogonal wavelets enhance the efficiency of the learning of CNN. of isotropic turbulence in a periodic box and the coarse-grained data are obtained by the application of a Gaussian low-pass filter to the DNS data. We assessed the WCNN3d in terms of 3D visualization of vorticity, PDF of vorticity, and enstrophy spectra. We found that WCNN3d well reproduces vorticity statistics and the positions of the vortices from coarse-grained vorticity fields. For the vorticity-based approach, the use of wavelets enhances deep learning of turbulent flows considered here thanks to the sparsity of the wavelet representation which improves deep learning. For the velocity-based approach, we showed that weighting the wavelet coefficients of velocity, which yields velocity gradient information due to norm equivalence, improves the accuracy and yields results similar to the vorticity-based model. Furthermore, we assessed the divergence issue of the predicted fields and showed that its impact is negligible. We demonstrated the capability to predict a turbulent flow whose Reynolds number is higher than the flows used for the training.","PeriodicalId":429847,"journal":{"name":"15th World Congress on Computational Mechanics (WCCM-XV) and 8th Asian Pacific Congress on Computational Mechanics (APCOM-VIII)","volume":"19 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"1900-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"115168276","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 1900-01-01DOI: 10.23967/wccm-apcom.2022.123
P. Kumari, V. Kumar
: E-vehicles and light weight structural parts in automotive and aerospace industry has led to design and development of new structures where traditional materials aluminium/steel are joined with composite laminated materials. This approach has led the engineers/researchers to reduce weight and ultimately save fuel consumption and reduce carbon footprints. Moreover, prosthetic limbs are also designed to have varying material along the length for better suitability. The above problems cannot be analysed using functionally graded theories/concepts. In theory, material properties vary linearly, exponentially, or power-law-like along x-values but for the above cases, material property does not follow a particular variation. Further even, it is not always practical to produce or manufacture components having very smooth variations along the length. In an attempt is made to develop the 2D analytical solution for multi-segmented Al/ steel- composite panel under transverse loading. Extended Kantorovich method is used for developing governing equations. Continuity of displacement and stresses are satisfied at interface of each segment. Two segmented panels having aluminium/steel and Gr/Ep equal and unequal segment are considered. The deflection and stresses are compared with the finite element solution and found in good agreement.
{"title":"Two-dimensional analytical solution for multi-segmented Al/ steel- composite panel-An Aerospace Application","authors":"P. Kumari, V. Kumar","doi":"10.23967/wccm-apcom.2022.123","DOIUrl":"https://doi.org/10.23967/wccm-apcom.2022.123","url":null,"abstract":": E-vehicles and light weight structural parts in automotive and aerospace industry has led to design and development of new structures where traditional materials aluminium/steel are joined with composite laminated materials. This approach has led the engineers/researchers to reduce weight and ultimately save fuel consumption and reduce carbon footprints. Moreover, prosthetic limbs are also designed to have varying material along the length for better suitability. The above problems cannot be analysed using functionally graded theories/concepts. In theory, material properties vary linearly, exponentially, or power-law-like along x-values but for the above cases, material property does not follow a particular variation. Further even, it is not always practical to produce or manufacture components having very smooth variations along the length. In an attempt is made to develop the 2D analytical solution for multi-segmented Al/ steel- composite panel under transverse loading. Extended Kantorovich method is used for developing governing equations. Continuity of displacement and stresses are satisfied at interface of each segment. Two segmented panels having aluminium/steel and Gr/Ep equal and unequal segment are considered. The deflection and stresses are compared with the finite element solution and found in good agreement.","PeriodicalId":429847,"journal":{"name":"15th World Congress on Computational Mechanics (WCCM-XV) and 8th Asian Pacific Congress on Computational Mechanics (APCOM-VIII)","volume":"35 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"1900-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"132897085","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 1900-01-01DOI: 10.23967/wccm-apcom.2022.119
J. Wackers, G. Deng, C. Raymond, E. Guilmineau, A. Leroyer, P. Queutey, M. Visonneau
. Adaptive grid refinement is tested for routine, automated simulations of ship resistance in calm water. A simulation protocol for these computations is fine-tuned on one test case and then applied unchanged to three different cases. The solutions are numerically accurate and compare well with experiments. Effective numerical uncertainty estimation increases the trustworthiness of the solutions.
{"title":"Towards automated computation with uncertainty estimation for industrial simulation of ship flow","authors":"J. Wackers, G. Deng, C. Raymond, E. Guilmineau, A. Leroyer, P. Queutey, M. Visonneau","doi":"10.23967/wccm-apcom.2022.119","DOIUrl":"https://doi.org/10.23967/wccm-apcom.2022.119","url":null,"abstract":". Adaptive grid refinement is tested for routine, automated simulations of ship resistance in calm water. A simulation protocol for these computations is fine-tuned on one test case and then applied unchanged to three different cases. The solutions are numerically accurate and compare well with experiments. Effective numerical uncertainty estimation increases the trustworthiness of the solutions.","PeriodicalId":429847,"journal":{"name":"15th World Congress on Computational Mechanics (WCCM-XV) and 8th Asian Pacific Congress on Computational Mechanics (APCOM-VIII)","volume":"30 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"1900-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"125042827","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 1900-01-01DOI: 10.23967/wccm-apcom.2022.075
A. Manuello, J. Melchiorre, L. Sardone, G. Marano
. Over the past decades, different approaches, physical and geometrical, were implemented to identify the optimal shape, reducing the internal stresses, of grid shells and vaults. As far as their original organic shape is concerned, the design of grid shell structures inspired architects and structural engineers worldwide and in any time. The method, here presented, is developed and extended, from its original formulation, employing a self-made code based on the dynamic equilibrium, ensured by the d'Alembert principle, of masses interconnected by rope elements in the space-time domain. The equilibrium corresponding the optimized shape to be defined, is obtained through an iterative process in the falling masses connected by a net for the definition of the "catenary surface" coinciding with the best shape of the shell (form minimizing the bending moment) according to the conditions of zero velocities and accelerations of the nodes. The implementation of the method is realized in MATLAB and set up for Python in an interpreted high-level general-purpose programming language. By the use of this code as well as its object-oriented architecture the MRA Python code will be linked to the Grasshopper environment for the direct visualization of the shapes and their fast-parametrization phase.
{"title":"Multi-body Rope Approach for the Form-Finding of Shape Optimized Grid Shell Structures","authors":"A. Manuello, J. Melchiorre, L. Sardone, G. Marano","doi":"10.23967/wccm-apcom.2022.075","DOIUrl":"https://doi.org/10.23967/wccm-apcom.2022.075","url":null,"abstract":". Over the past decades, different approaches, physical and geometrical, were implemented to identify the optimal shape, reducing the internal stresses, of grid shells and vaults. As far as their original organic shape is concerned, the design of grid shell structures inspired architects and structural engineers worldwide and in any time. The method, here presented, is developed and extended, from its original formulation, employing a self-made code based on the dynamic equilibrium, ensured by the d'Alembert principle, of masses interconnected by rope elements in the space-time domain. The equilibrium corresponding the optimized shape to be defined, is obtained through an iterative process in the falling masses connected by a net for the definition of the \"catenary surface\" coinciding with the best shape of the shell (form minimizing the bending moment) according to the conditions of zero velocities and accelerations of the nodes. The implementation of the method is realized in MATLAB and set up for Python in an interpreted high-level general-purpose programming language. By the use of this code as well as its object-oriented architecture the MRA Python code will be linked to the Grasshopper environment for the direct visualization of the shapes and their fast-parametrization phase.","PeriodicalId":429847,"journal":{"name":"15th World Congress on Computational Mechanics (WCCM-XV) and 8th Asian Pacific Congress on Computational Mechanics (APCOM-VIII)","volume":"77 1","pages":"0"},"PeriodicalIF":0.0,"publicationDate":"1900-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"125098918","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":0,"RegionCategory":"","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}