E. J. Sellountos¹, D. Polyzos²

CMES-Computer Modeling in Engineering & Sciences, Vol.4, No.6, pp. 619-636, 2003, DOI:10.3970/cmes.2003.004.619

Abstract A new meshless local Petrov-Galerkin (MLPG) method for solving two dimensional frequency domain elastodynamic problems is proposed. Since the method utilizes, in its weak formulation, either the elastostatic or the frequency domain elastodynamic fundamental solution as test function, it is equivalent to the local boundary integral equation (LBIE) method. Nodal points spread over the analyzed domain are considered and the moving least squares (MLS) interpolation scheme for the approximation of the interior and boundary variables is employed. Two integral equations suitable for the integral representation of the displacement fields in the local sub- domains are used. The first utilizes the… More >

Kazunori Shinohara^{1, *}

CMES-Computer Modeling in Engineering & Sciences, Vol.118, No.3, pp. 599-647, 2019, DOI:10.31614/cmes.2019.04472

Abstract According to the wave power rule, the second derivative of a function x(t) with respect to the variable t is equal to negative n times the function x(t) raised to the power of 2n-1. Solving the ordinary differential equations numerically results in waves appearing in the figures. The ordinary differential equation is very simple; however, waves, including the regular amplitude and period, are drawn in the figure. In this study, the function for obtaining the wave is called the leaf function. Based on the leaf function, the exact solutions for the undamped and unforced Duffing equations are presented. In the… More >

Yanxin Wang^{1, *}, Li Zhu¹, Zhi Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.118, No.2, pp. 339-350, 2019, DOI:10.31614/cmes.2019.04575

Abstract An Euler wavelets method is proposed to solve a class of nonlinear variable order fractional differential equations in this paper. The properties of Euler wavelets and their operational matrix together with a family of piecewise functions are first presented. Then they are utilized to reduce the problem to the solution of a nonlinear system of algebraic equations. And the convergence of the Euler wavelets basis is given. The method is computationally attractive and some numerical examples are provided to illustrate its high accuracy. More >

Lei Peng^1,2, Zhiqiang Feng^1,2,*, Pierre Joli², Christine Renaud²

CMES-Computer Modeling in Engineering & Sciences, Vol.118, No.1, pp. 111-137, 2019, DOI:10.31614/cmes.2018.04556

Abstract In order to investigate the mechanical behavior of systems with complex architecture and a large number of contacting bodies, a finite element software, named LiToTac, has been developed by using the object-oriented programming technique. This software, with an interactive graphical user interface, is able to handle highly non-linear problems including multiple contacts and large deformation. More importantly, the contact detection based on a hybrid three-stages methodology can be performed automatically, which is more efficient than the common strategies of pre-defining contact zones in commercial FEM software like ANSYS, ABAQUS, etc. In addition, the contact solver in LiToTac is portable between… More >

Yingjun Wang^1,*, Zhenpei Wang^2,*, Zhaohui Xia³, Leong Hien Poh²

CMES-Computer Modeling in Engineering & Sciences, Vol.117, No.3, pp. 455-507, 2018, DOI:10.31614/cmes.2018.04603

Abstract Isogeometric analysis (IGA), an approach that integrates CAE into conventional CAD design tools, has been used in structural optimization for 10 years, with plenty of excellent research results. This paper provides a comprehensive review on isogeometric shape and topology optimization, with a brief coverage of size optimization. For isogeometric shape optimization, attention is focused on the parametrization methods, mesh updating schemes and shape sensitivity analyses. Some interesting observations, e.g. the popularity of using direct (differential) method for shape sensitivity analysis and the possibility of developing a large scale, seamlessly integrated analysis-design platform, are discussed in the framework of isogeometric shape… More >

Adrien Scheuer^{1, 3, *}, Amine Ammar², Emmanuelle Abisset-Chavanne³, Elias Cueto⁴, Francisco Chinesta⁵, Roland Keunings¹, Suresh G. Advani⁶

CMES-Computer Modeling in Engineering & Sciences, Vol.117, No.3, pp. 367-386, 2018, DOI:10.31614/cmes.2018.04278

Abstract Describing the orientation state of the particles is often critical in fibre suspen-sion applications. Macroscopic descriptors, the so-called second-order orientation tensor (or moment) leading the way, are often preferred due to their low computational cost. Clo-sure problems however arise when evolution equations for the moments are derived from the orientation distribution functions and the impact of the chosen closure is often unpre-dictable. In this work, our aim is to provide macroscopic simulations of orientation that are cheap, accurate and closure-free. To this end, we propose an innovative data-based approach to the upscaling of orientation kinematics in the context of fibre… More >

Yang Yang^1,*, Jinbao Duan¹, Haitao Yu¹, Zhipeng Gao¹, Xuesong Qiu¹

CMES-Computer Modeling in Engineering & Sciences, Vol.117, No.3, pp. 555-575, 2018, DOI:10.31614/cmes.2018.04249

Abstract The development of deep learning has revolutionized image recognition technology. How to design faster and more accurate image classification algorithms has become our research interests. In this paper, we propose a new algorithm called stochastic depth networks with deep energy model (SADIE), and the model improves stochastic depth neural network with deep energy model to provide attributes of images and analysis their characteristics. First, the Bernoulli distribution probability is used to select the current layer of the neural network to prevent gradient dispersion during training. Then in the backpropagation process, the energy function is designed to optimize the target loss… More >

Qi Zhang¹, Bin Gao¹, Yue Shi¹, Chang Yu^1,*

CMES-Computer Modeling in Engineering & Sciences, Vol.116, No.2, pp. 247-262, 2018, DOI: 10.31614/cmes.2018.04082

Abstract Although pulsatile ECMO, as novel kinds of ECMO, has been attracted more and more attention, the differences of the hemodynamic effects of the pulsatile ECMO on the aorta, the cerebral perfusion, and left ventricular work were still under-investigated. The aim of this study was to clarify the hemodynamic differences of the cardiovascular system between the pulsatile and non-pulsatile VA ECMO. In this study, three ECMO support modes, named as “constant flow mode”, “co-pulse mode” and “counter pulse mode”, were designed. The computational fluid dynamics (CFD) study was carried out. The distribution of the oxygenated blood, the blood velocity vector, the… More >

Wei Wei¹, Qingsheng Yang^2,*, Xia Liu², Xiaoqiao He^3,*, Kim-Meow Liew³

CMES-Computer Modeling in Engineering & Sciences, Vol.115, No.3, pp. 387-409, 2018, DOI: 10.3970/cmes.2018.01856

Abstract In this paper, by introducing a chemical field, the J-integral formulation is presented for the chemo-mechanical coupled medium based on the laws of thermodynamics. A finite element implementation of the J-integral was performed to study the mode I chemo-mechanical coupled fracture problem. For derivation of the coupled J-integral, the equivalent domain integral (EDI) method was applied to obtain the mode I J-integral, with expression of the area integrals based on constitutive relationships of a linear elastic small deformation for chemo-mechanical coupling, instead of the finite deformation problem. A finite element procedure is developed to compute the mode I J-integral, and… More >

A. M. E. Ramírez-Mendoza¹

CMES-Computer Modeling in Engineering & Sciences, Vol.115, No.3, pp. 295-311, 2018, DOI: 10.3970/cmes.2018.00239

Abstract A spike response model (SRM) based on the spikes generator circuit (SGC) of adaptive fuzzy spiking neurons (AFSNs) is developed. The SRM is simulated in MatlabTM environment. The proposed model is applied to a configuration of a fuzzy exclusive or (fuzzy XOR) operator, as an illustrative example. A description of the comparison of AFSNs with other similar methods is given. The novel method of the AFSNs is used to determine the value of the weights or parameters of the fuzzy XOR, first with dynamic weights or self-tuning parameters that adapt continuously, then with fixed weights obtained after training, finally with… More >