Pub Date : 2018-01-01DOI: 10.2495/CMEM-V6-N6-1008-1018
A. Buchau, M. Jüttner
Nowadays, a variety of numerical methods and numerical formulations exits to solve complex or coupled field problems in three dimensions. Most of them are generally applicable to nearly arbitrary kind of field problems. On the other hand, some highly optimized methods are available, which are predestined for the solution of a specific kind of problem. Especially in the case of weakly coupled multiphysics problems, a mixture of several numerical methods is very advantages to benefit from different properties of numerical methods for diverse physical sub-problems. A very promising approach for a flexible coordination of the related solution process is the application of software agents. Then, the results of one sub-problem are converted into boundary values or volume source distributions for another sub-problem and software agents choose solution methods independently for each subproblem. Furthermore, two main aspects have to be considered in applications of numerical methods. First, the solution of a boundary value problem should be computed efficiently and second, the solution is evaluated for visualization and interpretation of obtained results. In practice, it is difficult to choose a single appropriate method, which is well suited both for the solution of a problem and its evaluation, since the demands differ in both cases. Here, a concept is presented to apply various numerical methods successfully to the solution and evaluation of complex field problems. Attention is mainly turned on the integration of boundary element methods into the concept of mixed numerical formulations.
{"title":"A concept of separated numerical formulations for the solution and evaluation of complex field problems","authors":"A. Buchau, M. Jüttner","doi":"10.2495/CMEM-V6-N6-1008-1018","DOIUrl":"https://doi.org/10.2495/CMEM-V6-N6-1008-1018","url":null,"abstract":"Nowadays, a variety of numerical methods and numerical formulations exits to solve complex or coupled field problems in three dimensions. Most of them are generally applicable to nearly arbitrary kind of field problems. On the other hand, some highly optimized methods are available, which are predestined for the solution of a specific kind of problem. Especially in the case of weakly coupled multiphysics problems, a mixture of several numerical methods is very advantages to benefit from different properties of numerical methods for diverse physical sub-problems. A very promising approach for a flexible coordination of the related solution process is the application of software agents. Then, the results of one sub-problem are converted into boundary values or volume source distributions for another sub-problem and software agents choose solution methods independently for each subproblem. Furthermore, two main aspects have to be considered in applications of numerical methods. First, the solution of a boundary value problem should be computed efficiently and second, the solution is evaluated for visualization and interpretation of obtained results. In practice, it is difficult to choose a single appropriate method, which is well suited both for the solution of a problem and its evaluation, since the demands differ in both cases. Here, a concept is presented to apply various numerical methods successfully to the solution and evaluation of complex field problems. Attention is mainly turned on the integration of boundary element methods into the concept of mixed numerical formulations.","PeriodicalId":22520,"journal":{"name":"THE INTERNATIONAL JOURNAL OF COMPUTATIONAL METHODS AND EXPERIMENTAL MEASUREMENTS","volume":"32 1","pages":"1008-1018"},"PeriodicalIF":0.0,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91540157","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 : 2018-01-01DOI: 10.2495/CMEM-V6-N6-1043-1056
R. Vodička, Filip Kšiňan
A model for numerical analysis of compound structures made of various materials is presented. The mathematical concept of solution is based on quasi-static evolution of debonding processes occurring along the interface. It is formulated in terms of energies considering the stored energy represented by the elastic energy of the structures and dissipation due to damage processes, plastic slip at the interface or friction. The numerical solution includes a semi-implicit time stepping procedure, relying on splitting of the whole problem at a current time step into two problems of variational nature solved recursively. The space discretisation includes Symmetric Galerkin Boundary Element Method used to obtain the stored energies, and, in combination with the variational character of the recursive problems, also to calculate its gradients to be utilized in non-linear programming algorithms for finding the timeevolving solution. Numerical results are demonstrated for a steel-concrete interface frequently met in civil engineering applications to assess the model applicability in engineering practice.
{"title":"A quasi-static interface damage model with frictional contact – applications to steel reinforced concrete structures","authors":"R. Vodička, Filip Kšiňan","doi":"10.2495/CMEM-V6-N6-1043-1056","DOIUrl":"https://doi.org/10.2495/CMEM-V6-N6-1043-1056","url":null,"abstract":"A model for numerical analysis of compound structures made of various materials is presented. The mathematical concept of solution is based on quasi-static evolution of debonding processes occurring along the interface. It is formulated in terms of energies considering the stored energy represented by the elastic energy of the structures and dissipation due to damage processes, plastic slip at the interface or friction. The numerical solution includes a semi-implicit time stepping procedure, relying on splitting of the whole problem at a current time step into two problems of variational nature solved recursively. The space discretisation includes Symmetric Galerkin Boundary Element Method used to obtain the stored energies, and, in combination with the variational character of the recursive problems, also to calculate its gradients to be utilized in non-linear programming algorithms for finding the timeevolving solution. Numerical results are demonstrated for a steel-concrete interface frequently met in civil engineering applications to assess the model applicability in engineering practice.","PeriodicalId":22520,"journal":{"name":"THE INTERNATIONAL JOURNAL OF COMPUTATIONAL METHODS AND EXPERIMENTAL MEASUREMENTS","volume":"63 1","pages":"1043-1056"},"PeriodicalIF":0.0,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"90708478","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 : 2018-01-01DOI: 10.2495/CMEM-V6-N6-1127-1137
K. Matsushima, H. Isakari, Toru Takahashi, Toshiro Matsumoto
In this study, we investigate the distribution of eigenfrequencies of boundary integral equations (BIEs) of two-dimensional elastodynamics. The corresponding eigenvalue problem is classified as a nonlinear eigenvalue problem. We confirm that the Burton-Miller formulation can properly avoid fictitious eigenfrequencies. The boundary element method (BEM) is expected as a powerful numerical tool for designing sophisticated devices related to elastic waves such as acoustic metamaterials. However, the BEM is known that it loses its accuracy for certain frequencies, called as fictitious eigenfrequencies, for problems defined in the infinite domain. Recent researches It has also been revealed that not only the real-valued eigenfrequencies but also the complex-valued ones may affect the accuracy of the BEM results. We examine the distribution of complex eigenvalues obtained by BIEs for time-harmonic elastodynamic problems with the help of the Sakurai-Sugiura method which is applicable to nonlinear eigenvalue problems. We also examine its relation to the accuracy of the BEM numerical results. We also discuss an appropriate choice of the coupling parameter from a viewpoint of the distribution of fictitious eigenfrequencies.
{"title":"An investigation of Eigenfrequencies of boundary integral equations and the Burton-Miller formulation in two-dimensional elastodynamics","authors":"K. Matsushima, H. Isakari, Toru Takahashi, Toshiro Matsumoto","doi":"10.2495/CMEM-V6-N6-1127-1137","DOIUrl":"https://doi.org/10.2495/CMEM-V6-N6-1127-1137","url":null,"abstract":"In this study, we investigate the distribution of eigenfrequencies of boundary integral equations (BIEs) of two-dimensional elastodynamics. The corresponding eigenvalue problem is classified as a nonlinear eigenvalue problem. We confirm that the Burton-Miller formulation can properly avoid fictitious eigenfrequencies. The boundary element method (BEM) is expected as a powerful numerical tool for designing sophisticated devices related to elastic waves such as acoustic metamaterials. However, the BEM is known that it loses its accuracy for certain frequencies, called as fictitious eigenfrequencies, for problems defined in the infinite domain. Recent researches It has also been revealed that not only the real-valued eigenfrequencies but also the complex-valued ones may affect the accuracy of the BEM results. We examine the distribution of complex eigenvalues obtained by BIEs for time-harmonic elastodynamic problems with the help of the Sakurai-Sugiura method which is applicable to nonlinear eigenvalue problems. We also examine its relation to the accuracy of the BEM numerical results. We also discuss an appropriate choice of the coupling parameter from a viewpoint of the distribution of fictitious eigenfrequencies.","PeriodicalId":22520,"journal":{"name":"THE INTERNATIONAL JOURNAL OF COMPUTATIONAL METHODS AND EXPERIMENTAL MEASUREMENTS","volume":"32 1","pages":"1127-1137"},"PeriodicalIF":0.0,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"78103187","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 : 2018-01-01DOI: 10.2495/CMEM-V6-N6-1182-1191
M. Cvetković, D. Poljak
ABSTRACT The paper revisits the use of a surface equivalence theorem in deriving the surface integral equation (SIE) based formulation for a homogeneous bio-electromagnetics problem. The vector analog of Green’s 2nd identity is used to obtain the expression for the electric field representing the mathematical foundation of the equivalence theorem. The particular emphasis is put on the treatment of boundary integral when the observation and source points, respectively, coincide. The boundary conditions at infinity are taken into account via the Sommerfeld radiation conditions. The derived coupled SIE set can be used in problems involving biological body exposed to electromagnetic field radiation.
{"title":"Surface equivalence principle and surface integral equation (SIE) revisited for bioelectromagnetics application","authors":"M. Cvetković, D. Poljak","doi":"10.2495/CMEM-V6-N6-1182-1191","DOIUrl":"https://doi.org/10.2495/CMEM-V6-N6-1182-1191","url":null,"abstract":"ABSTRACT The paper revisits the use of a surface equivalence theorem in deriving the surface integral equation (SIE) based formulation for a homogeneous bio-electromagnetics problem. The vector analog of Green’s 2nd identity is used to obtain the expression for the electric field representing the mathematical foundation of the equivalence theorem. The particular emphasis is put on the treatment of boundary integral when the observation and source points, respectively, coincide. The boundary conditions at infinity are taken into account via the Sommerfeld radiation conditions. The derived coupled SIE set can be used in problems involving biological body exposed to electromagnetic field radiation.","PeriodicalId":22520,"journal":{"name":"THE INTERNATIONAL JOURNAL OF COMPUTATIONAL METHODS AND EXPERIMENTAL MEASUREMENTS","volume":"32 1","pages":"1182-1191"},"PeriodicalIF":0.0,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"80299895","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 : 2018-01-01DOI: 10.2495/CMEM-V6-N6-1120-1126
L. Alexeyeva, B. Alipova
The dynamics of multi-connected thermoelastic semiplane with the non-stationary power source and thermal effects by using of a model of coupled thermoelasticity is investigated. Green’s tensor in the space of the Laplace transforms in time describes the displacements of medium under the effect of the impulse concentrated power and thermal sources. The generalized solution of the problem of the dynamics of thermoelastic semiplane with the free boundary under the effect of arbitrary mass forces and thermal sources in 2D-case is built.
{"title":"Boundary integral equations of dynamics problems for multi-connected thermoelastic semi-plane with a free boundary","authors":"L. Alexeyeva, B. Alipova","doi":"10.2495/CMEM-V6-N6-1120-1126","DOIUrl":"https://doi.org/10.2495/CMEM-V6-N6-1120-1126","url":null,"abstract":"The dynamics of multi-connected thermoelastic semiplane with the non-stationary power source and thermal effects by using of a model of coupled thermoelasticity is investigated. Green’s tensor in the space of the Laplace transforms in time describes the displacements of medium under the effect of the impulse concentrated power and thermal sources. The generalized solution of the problem of the dynamics of thermoelastic semiplane with the free boundary under the effect of arbitrary mass forces and thermal sources in 2D-case is built.","PeriodicalId":22520,"journal":{"name":"THE INTERNATIONAL JOURNAL OF COMPUTATIONAL METHODS AND EXPERIMENTAL MEASUREMENTS","volume":"48 19","pages":"1120-1126"},"PeriodicalIF":0.0,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"91455618","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 : 2018-01-01DOI: 10.2495/CMEM-V6-N6-1079-1086
Shin Yina, W. Pei, Zhao Qiang, F. Qijing
This work investigates the mass ejected from surface perturbations as the shockwave reaches the AL-vacuum interface, which originates from unstable Richtmyer–Meshkov (RMI) impulse phenomena. The main purpose is to explore the relationships between the shockwave impulse and the geometric properties of surface perturbations, and how those relationships drive the total ejected mass, directionality and velocity distribution. We discuss in detail different types of surface geometry (sinusoidal, square-wave, chevron and semicircle), as well as the wavelengths and amplitudes of surface perturbation. The time evolutions of micro-jet ejection are simulated using a hydrodynamic Lagrangian-Remapping Eulerian method. The calculated results show that primary jetting ejection can be formed from the different shapes, and with increasing wavelength, the ejection mass keeps an increase while the jet head-velocity decreases. However, not all periodic perturbations behave similarly, and masses ejected from irregular surface cannot be normalized to its cross-sectional areas. The square-wave surface may yield pronounced, velocity-enhanced secondary jetting, which is a result of collision of primary jets.
{"title":"Exploring Richtmyer–meshkov Instability Phenomena And The Links Between Surface Perturbations And Shocked-induced Mass Ejection","authors":"Shin Yina, W. Pei, Zhao Qiang, F. Qijing","doi":"10.2495/CMEM-V6-N6-1079-1086","DOIUrl":"https://doi.org/10.2495/CMEM-V6-N6-1079-1086","url":null,"abstract":"This work investigates the mass ejected from surface perturbations as the shockwave reaches the AL-vacuum interface, which originates from unstable Richtmyer–Meshkov (RMI) impulse phenomena. The main purpose is to explore the relationships between the shockwave impulse and the geometric properties of surface perturbations, and how those relationships drive the total ejected mass, directionality and velocity distribution. We discuss in detail different types of surface geometry (sinusoidal, square-wave, chevron and semicircle), as well as the wavelengths and amplitudes of surface perturbation. The time evolutions of micro-jet ejection are simulated using a hydrodynamic Lagrangian-Remapping Eulerian method. The calculated results show that primary jetting ejection can be formed from the different shapes, and with increasing wavelength, the ejection mass keeps an increase while the jet head-velocity decreases. However, not all periodic perturbations behave similarly, and masses ejected from irregular surface cannot be normalized to its cross-sectional areas. The square-wave surface may yield pronounced, velocity-enhanced secondary jetting, which is a result of collision of primary jets.","PeriodicalId":22520,"journal":{"name":"THE INTERNATIONAL JOURNAL OF COMPUTATIONAL METHODS AND EXPERIMENTAL MEASUREMENTS","volume":"104 1","pages":"1079-1086"},"PeriodicalIF":0.0,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"87661049","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 : 2018-01-01DOI: 10.2495/CMEM-V6-N6-989-999
Y. Ochiai
{"title":"Meshless large plastic deformation analysis considering with a friction coefficient by triple-reciprocity boundary element method","authors":"Y. Ochiai","doi":"10.2495/CMEM-V6-N6-989-999","DOIUrl":"https://doi.org/10.2495/CMEM-V6-N6-989-999","url":null,"abstract":"","PeriodicalId":22520,"journal":{"name":"THE INTERNATIONAL JOURNAL OF COMPUTATIONAL METHODS AND EXPERIMENTAL MEASUREMENTS","volume":"114 1","pages":"989-999"},"PeriodicalIF":0.0,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"77642662","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 : 2018-01-01DOI: 10.2495/CMEM-V6-N6-1067-1078
Z. Yan
{"title":"Simulation of sound structure interactions by the coupled FEM/BEM","authors":"Z. Yan","doi":"10.2495/CMEM-V6-N6-1067-1078","DOIUrl":"https://doi.org/10.2495/CMEM-V6-N6-1067-1078","url":null,"abstract":"","PeriodicalId":22520,"journal":{"name":"THE INTERNATIONAL JOURNAL OF COMPUTATIONAL METHODS AND EXPERIMENTAL MEASUREMENTS","volume":"38 1","pages":"1067-1078"},"PeriodicalIF":0.0,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"74045515","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 : 2018-01-01DOI: 10.2495/cmem-v6-n6-1173-1181
A. Galybin
This study investigates solvability of boundary value problems of plane elasticity formulated in terms of principal directions of the stress tensor and the orientations of the displacement vector. The analysis of solvability is performed by using the following approach. Firstly, boundary values of the complex potentials are represented by the Cauchy-type integrals with unknown density. Then a system of singular integral equations is obtained by satisfying particular boundary conditions. This system is further reduced to the system of the Riemann boundary value problems for the determination of sectionally holomorphic functions. Solvability of the Riemann problems is investigated by calculating their indexes. This allows one to determine the number of linearly independent solutions and hence the number of arbitrary parameters entering into the general solution. Two novel formulations have been investigated for the case of elastic half-planes. In both cases the initial system of equations has been reduced to the form that allow for successive solution of its equations.
{"title":"Boundary value problems for elastic half-planes posed in terms of stress and displacement orientations","authors":"A. Galybin","doi":"10.2495/cmem-v6-n6-1173-1181","DOIUrl":"https://doi.org/10.2495/cmem-v6-n6-1173-1181","url":null,"abstract":"This study investigates solvability of boundary value problems of plane elasticity formulated in terms of principal directions of the stress tensor and the orientations of the displacement vector. The analysis of solvability is performed by using the following approach. Firstly, boundary values of the complex potentials are represented by the Cauchy-type integrals with unknown density. Then a system of singular integral equations is obtained by satisfying particular boundary conditions. This system is further reduced to the system of the Riemann boundary value problems for the determination of sectionally holomorphic functions. Solvability of the Riemann problems is investigated by calculating their indexes. This allows one to determine the number of linearly independent solutions and hence the number of arbitrary parameters entering into the general solution. Two novel formulations have been investigated for the case of elastic half-planes. In both cases the initial system of equations has been reduced to the form that allow for successive solution of its equations.","PeriodicalId":22520,"journal":{"name":"THE INTERNATIONAL JOURNAL OF COMPUTATIONAL METHODS AND EXPERIMENTAL MEASUREMENTS","volume":"37 1","pages":"1173-1181"},"PeriodicalIF":0.0,"publicationDate":"2018-01-01","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"81287715","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}