• Geomechanics and Engineering
• /
• v.29 no.1
• /
• pp.1-11
• /
• 2022

A semi-analytical solution to responses of overconsolidated (OC) unsaturated soils surrounding an expanding spherical cavity under constant suction condition is presented. To capture the elastoplastic hydro-mechanical property of OC unsaturated soils, the unified hardening (UH) model for OC unsaturated soil is adopted in corporation with a soil-water characteristic curve (SWCC) and two suction yield surfaces. Taking the specific volume, radial stress, tangential stress and degree of saturation as the four basic unknowns, the problem investigated is formulated by solving a set of first-order ordinary differential equations with the help of an auxiliary variable and an iterative algorithm. The present solution is validated by comparing with available solution based on the modified Cam Clay (MCC) model. Parametric studies reveal that the hydraulic and mechanical responses of spherical cavity expanding in unsaturated soils are not only coupled, but also affected by suction and overconsolidation ratio (OCR) significantly. More importantly, whether hydraulic yield will occur or not depends only on the initial relationship between suction yield stress and suction. The presented solution can be used for calibration of some insitu tests in OC unsaturated soil.

• Coupled systems mechanics
• /
• v.12 no.5
• /
• pp.419-428
• /
• 2023

We have developed a model for estimating the parameters of viscous materials from indirect tensile tests for asphalt. This is a simple Burger nonlinear rheological two-cell model or standard model. At the same time, we begin to develop a more versatile and complex multi-cell model. The simple model is validated using experimental load-displacement results from laboratory tests: The recorded displacements are used as input values and the measured force data are simulated with the model. The optimal model parameters are estimated using the Levenberg-Marquardt method and a very good agreement between the experimental results and the model calculations is shown. However, not all parts of the model are active in the loading phase of the experiment, so we extended the validation of the model to the simulation of the relaxation behaviour. In this stage, the other model parameters are activated and the simulation results are consistent with the literature. At this stage, we have estimated the parameters only for the two-cell uniaxial model, but further work will include results for the multi-cell model.

• Coupled systems mechanics
• /
• v.13 no.2
• /
• pp.171-186
• /
• 2024

In the present paper, thermal buckling characteristics of functionally graded rectangular plates made of porous material that are integrated with surface-bonded piezoelectric actuators subjected to the combined action of thermal load and constant applied actuator voltage are investigated by utilizing a Navier solution method. The uniform temperature rise loading is considered. Thermomechanical material properties of FGM plates are assumed to be temperature independent and supposed to vary through thickness direction of the constituents according to power-law distribution (P-FGM) which is modified to approximate the porous material properties with even and uneven distributions of porosities phases. The governing differential equations of stability for the piezoelectric FGM plate are derived based on higher order shear deformation plate theory. Influences of several important parameters on the critical thermal buckling temperature are investigated and discussed in detail.

• Structural Engineering and Mechanics
• /
• v.71 no.6
• /
• pp.627-641
• /
• 2019

The paper focuses on bending analysis of the functionally graded (FG) plates with arbitrary shapes and boundary conditions. The material property of FG plates is modelled by using the power law distribution. Based on the first order shear deformation plate theory (FSDT), the governing equations as well as boundary conditions are formulated and obtained by using the principle of virtual work. The coupled Boundary Element-Radial Basis Function (BE-RBF) method is established to solve the complex FG plates. The proposed methodology is developed by applying the concept of the analog equation method (AEM). According to the AEM, the original governing differential equations are replaced by three Poisson equations with fictitious sources under the same boundary conditions. Then, the fictitious sources are established by the application of a technique based on the boundary element method and approximated by using the radial basis functions. The solution of the actual problem is attained from the known integral representations of the potential problem. Therefore, the kernels of the boundary integral equations are conveniently evaluated and readily determined, so that the complex FG plates can be easily computed. The reliability of the proposed method is evaluated by comparing the present results with those from analytical solutions. The effects of the power index, the length to thickness ratio and the modulus ratio on the bending responses are investigated. Finally, many interesting features and results obtained from the analysis of the FG plates with arbitrary shapes and boundary conditions are demonstrated.

• Advances in nano research
• /
• v.13 no.3
• /
• pp.285-295
• /
• 2022

Creation of plastic deformation under seismic loads, is one of the most serious subjects in RC structures with steel bars which reduces the life threatening risks and increases dissipation of energy. Shape memory alloy (SMA) is one of the best choice for the relocating plastic hinges. In a challenge to study the seismic response of concrete moment resisting frame (MRF), this article investigates numerically a new type of concrete frames with nano fiber reinforced polymer (NFRP) and shape memory alloy (SMA) hinges, simultaneously. The NFRP layer is containing carbon nanofibers with agglomeration based on Mori-Tanaka model. The tangential shear deformation (TASDT) is applied for modelling of the structure and the continuity boundary conditions are used for coupling of the motion equations. In SMA connections between beam and columns, since there is phase transformation, hence, the motion equations of the structure are coupled with kinetic equations of phase transformation. The Hernandez-Lagoudas theory is applied for demonstrating of pseudoelastic characteristics of SMA. The corresponding motion equations are solved by differential cubature (DC) and Newmark methods in order to obtain the peak ground acceleration (PGA) and residual drift ratio for MRF-2%. The main impact of this paper is to present the influences of the volume percent and agglomeration of nanofibers, thickness and length of the concrete frame, SMA material and NFRP layer on the PGA and drift ratio. The numerical results revealed that the with increasing the volume percent of nanofibers, the PGA is enhanced and the residual drift ratio is reduced. It is also worth to mention that PGA of concrete frame with NFRP layer containing 2% nanofibers is approximately equal to the concrete frame with steel bars.

• Journal of the Korean Society of Propulsion Engineers
• /
• v.22 no.3
• /
• pp.46-52
• /
• 2018

A simple Crocco's $n-{\tau}$ time delay model and linear analysis of fluid flow coupled with acoustics are combined to investigate the high frequency combustion instability in the combustion chamber of LOX/hydrocarbon engines. The partial differential equation of the velocity potential is separated into ordinary differential equations, and eigenvalues that correspond to tangential resonance modes in the cylindrical chamber are determined. A general solution is obtained by solving the differential equation in the axial direction, and boundary conditions at the injector face and nozzle entrance are applied in order to calculate the chamber admittance. Frequency analysis of the transfer function is used to evaluate the stability of system. Stability margin is determined from the system gain and phase angle for the desired frequency range of 1T mode. The chamber model with variable baffle length and configurations are also considered in order to enhance the 1T mode stability of the combustion chamber.

• The Transactions of the Korea Information Processing Society
• /
• v.2 no.4
• /
• pp.535-542
• /
• 1995

In this paper we present two preconditioners for solving large sparse linear systems arising from elliptic partial differential equations on massively parallel machines, such as the CM-5. Most massively parallel machines do heavily rely on the message-passing for the interprocessor communications. but according to the current manufacturing standards the cost of communications is very high compared to that of floating point arithmetic computations. Due to this we need an algorithm which minimizes the amount of interprocessor communication on the massively parallel machines. We will show that Block SOR(Successive Over Relaxation) method coupled with the multi-coloring technique is one of such preconditioner on the massively parallel machines, by conducting experiments in the CM-5. Also, we implemented the ADI(Alternation Direction Implicit) method in the CM-5, which has been conventionally one of the most powerful parallel preconditioner. Our experiment shows that Block SOR method coupled with the multi-coloring technique could yield a speedup with 50% efficiency with the range of number of processors form 16 to 512 for a matrix with dimension 512x512. On the other hand, the ADI method shows a very poor performance.

• Smart Structures and Systems
• /
• v.18 no.2
• /
• pp.355-374
• /
• 2016

This paper presents and compares a one-dimensional (1D) bending theory for piezoelectric thin beam-type structures with resistive-inductive electrodes to ANSYS$^{(R)}$ three-dimensional (3D) finite element (FE) analysis. In particular, the lateral deflections and vibrations of slender piezoelectric beams are considered. The peculiarity of the piezoelectric beam model is the modeling of electrodes in such a manner that is does not fulfill the equipotential area condition. The case of ideal, perfectly conductive electrodes is a special case of our 1D model. Two-coupled partial differential equations are obtained for the lateral deflection and for the voltage distribution along the electrodes: the first one is an extended Bernoulli-Euler beam equation (second-order in time, forth order in space) and the second one the so-called Telegrapher's equation (second-order in time and space). Analytical results of our theory are validated by 3D electromechanically coupled FE simulations with ANSYS$^{(R)}$. A clamped-hinged beam is considered with various types of electrodes for the piezoelectric layers, which can be either resistive and/or inductive. A natural frequency analysis as well as quasi-static and dynamic simulations are performed. A good agreement between the extended beam theory and the FE results is found. Finally, the practical relevance of this type of electrodes is shown. It is found that the damping capability of properly tuned resistive or resistive-inductive electrodes exceeds the damping performance of beams, where the electrodes are simply linked to an optimized impedance.

• The Transactions of the Korea Information Processing Society
• /
• v.2 no.6
• /
• pp.960-968
• /
• 1995

In this paper we present a reordering scheme that could lead to efficient vectorization of the preconditioners for the large sparse linear systems arising from partial differential equations on the CRAY-2, This reordering scheme is a line version of the conventional red/black ordering. This reordering scheme, coupled with a variant of ILU(Incomplete LU) preconditioning, can overcome the poor rate of convergence of the conventional red/black reordering, if relatively large number of fill-ins were used. We substantiate our claim by conducting various experiments on the CRAY-2 machine. Also, the computation of the Frobenius norm of the error matrices agree with our claim.

• Smart Structures and Systems
• /
• v.16 no.6
• /
• pp.1023-1047
• /
• 2015

In this study, the optimal location of the MR fluid segments in a partially treated laminated composite sandwich plate has been identified to maximize the natural frequencies and the loss factors. The finite element formulation is used to derive the governing differential equations of motion for a partially treated laminated composite sandwich plate embedded with MR fluid and rubber material as the core layer and laminated composite plate as the face layers. An optimization problem is formulated and solved by combining finite element analysis (FEA) and genetic algorithm (GA) to obtain the optimal locations to yield maximum natural frequency and loss factor corresponding to first five modes of flexural vibration of the sandwich plate with various combinations of weighting factors under various boundary conditions. The proposed methodology is validated by comparing the natural frequencies evaluated at optimal locations of MR fluid pockets identified through GA coupled with FEA and the experimental measurements. The converged results suggest that the optimal location of MR fluid pockets is strongly influenced not only by the boundary conditions and modes of vibrations but also by the objectives of maximization of natural frequency and loss factors either individually or combined. The optimal layout could be useful to apply the MR fluid pockets at critical components of large structure to realize more efficient and compact vibration control mechanism with variable damping.