• Geomechanics and Engineering
• /
• v.8 no.4
• /
• pp.559-594
• /
• 2015

A nonlinear creep function incorporated into the elastic visco-plastic model may describe the long-term soil deformation more accurately. However, by applying the conventional procedure, there are challenges to determine the model parameters due to limitation of suitable data points. This paper presents a numerical solution to obtain several parameters simultaneously for a nonlinear elastic visco-plastic (EVP) model using the available consolidation data. The finite difference scheme using the Crank-Nicolson procedure is applied to solve a set of coupled partial differential equations of the time dependent strain and pore water pressure dissipation. The model parameters are determined by applying the algorithm of trust-region reflective optimisation in conjunction with the finite difference solution. The proposed method utilises all available consolidation data during dissipation of the excess pore water pressure to determine the required model parameters. Moreover, the reference time in the elastic visco-plastic model can readily be adopted as a unit of time; denoting creep is included in the numerical predictions explicitly from the very first time steps. In this paper, the settlement predictions of thick soft clay layers are presented and discussed to evaluate and compare the accuracy and reliability of the proposed method against the graphical procedure to obtain the model parameters. In addition, comparison of the available experimental results to the numerical predictions confirms the accuracy of the numerical procedure.

• Smart Structures and Systems
• /
• v.16 no.4
• /
• pp.701-723
• /
• 2015

Friction isolators are one of the most important types of bearings used to mitigate damages of earthquakes. The adaptive behavior of these isolators allows them to achieve multiple levels of performances and predictable seismic behavior during different earthquake hazard levels. There are three main types of friction isolators. The first generation with one sliding surface is known as Friction Pendulum System (FPS) isolators. The double concave friction pendulum (DCFP) with two sliding surfaces is an advanced form of FPS, and the third one, with fully adaptive behavior, is named as triple concave friction pendulum (TCFP). The current study has been conducted to investigate and compare seismic responses of these three types of isolators. The structure is idealized as a two-dimensional single degree of freedom (SDOF) resting on isolators. The coupled differential equations of motion are derived and solved using state space formulation. Seismic responses of isolated structures using each one of these isolators are investigated under seven near fault earthquake motions. The peak values of bearing displacement and base shear are studied employing the variation of essential parameters such as superstructure period, effective isolation period and effective damping of isolator. The results demonstrate a more efficient seismic behavior of TCFP isolator comparing to the other types of isolators. This efficiency depends on the selected effective isolation period as well as effective isolation damping. The investigation shows that increasing the effective isolation period or decreasing the effective isolation damping improves the seismic behavior of TCFP compared to the other isolators. The maximum difference in seismic responses, the base shear and the bearing displacement, for the TCFP isolator are calculated 26.8 and 13.4 percent less than the DCFP and FPS in effective isolation damping equal to10%, respectively.

• Transactions of the Korean Society of Mechanical Engineers B
• /
• v.29 no.9 s.240
• /
• pp.1022-1031
• /
• 2005

Temporal behavior of the laser induced incandescence (LII) signal is often used for soot particle sizing, which is possible because the cooling behavior of a laser heated particle is dependent on the particle size. In present study, LII signals of soot particles are modeled using two non-linear coupled differential equations deduced from the energy- and mass-balance of the process. The objective of this study is to obtain an appropriate calibration curve for determining primary particle size by comparing the gated signal ratio and double-exponential curve fitting methods. Not only the effects of laser fluence and gas temperature on the cooling behavior but also heat transfer mechanisms of heated soot particle have been investigated. The second-order exponential curve fitting showed better agreements with the LII signals than the gated signal ratio method which was based on the lust-order exponential curve fit. And the temporal decay rate of the LII signal and primary particle size showed nearly linear relationship, which was little dependent on the laser fluence. And it also could be reconfirmed that vaporization was dominant process of heat loss during first loons after laser pulse, then heat conduction played most important role while thermal radiation had little influence all the time.

• Journal of Environmental Science International
• /
• v.5 no.5
• /
• pp.593-603
• /
• 1996

A dual-porosity filmed agglomerate model for the porous cathode of the molten carbonate fuel has been investigated to predict the cell performance. A phenomenological treatment of molecular, kinetic and electrode parameters has been given. The major physical and chemical phenomena being modeled include mass transfer, ohmic losses and reaction kinetics at the electrode- electrolyte interface. The model predicts steady-state cell performance, given the above conditions that characterize the state of the electrode. Quasi-linearization and finite difference techniques are used to solve the coupled nonlinear differential equations. Also, the effective surface area of electrode pore was obtained by mercury porosimeter. The results of the investigation are presented in the form of plots of overpotential vs. current density with varied the electrode material, gas composition and mechanism. The predicted polarization curves are compared with the empirical data from 1 c$m^2$ cell. A fair correspondence is observed.

• Structural Engineering and Mechanics
• /
• v.51 no.2
• /
• pp.277-297
• /
• 2014

Since the isolation bearings undergo large displacements in base-isolated structures, impact with adjacent structures is inevitable. Therefore, in this investigation, the effect of impact on seismic response of isolated structures mounted on double concave friction pendulum (DCFP) bearings subjected to near field ground motions is considered. A non-linear viscoelastic model of collision is used to simulate structural pounding more accurately. 2-, 4- and 8-story base-isolated buildings adjacent to fixed-base structures are modeled and the coupled differential equations of motion related to these isolated systems are solved in the MATLAB environment using the SIMULINK toolbox. The variation of seismic responses such as base shear, displacement in the isolation system and superstructure (top floor) is computed to study the impact condition. Also, the effects of variation of system parameters: isolation period, superstructure period, size of seismic gap between two structures, radius of curvature of the sliding surface and friction coefficient of isolator are contemplated in this study. It is concluded that the normalized base shear, bearing and top floor displacement increase due to impact with adjacent structure. When the distance between two structures decreases, the base shear and displacement increase comparing to no impact condition. Besides, the increase in friction coefficient difference also causes the normalized base shear and displacement in isolation system and superstructure increase in comparison with bi-linear hysteretic behavior of base isolation system. Totally, the comparison of results indicates that the changes in values of friction coefficient have more significant effects on 2-story building than 4- and 8-story buildings.

• Journal of Biomedical Engineering Research
• /
• v.16 no.4
• /
• pp.415-424
• /
• 1995

Combined models, specified by two or more modeling formalisms, can represent a wide variety of complex systems. This paper describes a methodology for the development of combined models in two model types of discrete event and continuous process. The methodology is based on transformation of continuous state space into discrete one to homomorphically represent dynamics of continuous processes in discrete events. This paper proposes a formal structure which can combine model of the DES and the CS within a framework. The structure employs the DEVS formalism for the DES models and differential or polynomial equations for the CS models. To employ the proposed structure to specify a DEVS/CS combined model, a modeler needs to take the following steps. First, a modeler should identify events in the CS and transform the states of the CS into the DES. Second, a modular employs the formalism to specify the system as the DES. Finally, a moduler developes sub-models for the CS and continguos states of the DES and establishs one-to-one correspondence between the sub-models and such states. The proposed formal structre has been applied to develop a DEVS/CS combined model for the human cardiovascular system. For this, the cardiac cycle is partitioned into a set of phases based on events identified through observation. For each phase, a CS model has been developed and associated with the phase. To validate the DEVS/CS combined model developed, then simulate the model in the DEVSIM + + environment, which is a model simulation results with the results obtained from the CS model simulation using SPICE. The comparison shows that the DEVS/CS combined model adequately represents dynamics of the human heart system at each phase of cardiac cycle.

• Ocean Systems Engineering
• /
• v.7 no.4
• /
• pp.435-460
• /
• 2017

The main objective of this study is to investigate the turning and zig-zag maneuvering performance of the well-known naval surface combatant DTMB (David Taylor Model Basin) 5415 hull with URANS (Unsteady Reynolds-averaged Navier-Stokes) method. Numerical simulations of static drift tests have been performed by a commercial RANS solver based on a finite volume method (FVM) in an unsteady manner. The fluid flow is considered as 3-D, incompressible and fully turbulent. Hydrodynamic analyses have been carried out for a fixed Froude number 0.28. During the analyses, the free surface effects have been taken into account using VOF (Volume of Fluid) method and the hull is considered as fixed. First, the code has been validated with the available experimental data in literature. After validation, static drift, static rudder and drift and rudder tests have been simulated. The forces and moments acting on the hull have been computed with URANS approach. Numerical results have been applied to determine the hydrodynamic maneuvering coefficients, such as, velocity terms and rudder terms. The acceleration, angular velocity and cross-coupled terms have been taken from the available experimental data. A computer program has been developed to apply a fast maneuvering simulation technique. Abkowitz's non-linear mathematical model has been used to calculate the forces and moment acting on the hull during the maneuvering motion. Euler method on the other hand has been applied to solve the simultaneous differential equations. Turning and zig-zag maneuvering simulations have been carried out and the maneuvering characteristics have been determined and the numerical simulation results have been compared with the available data in literature. In addition, viscous effects have been investigated using Eulerian approach for several static drift cases.

• Coupled systems mechanics
• /
• v.9 no.3
• /
• pp.237-264
• /
• 2020

This article investigates the static behaviour of functionally graded (FG) plates sometimes declared as advanced composite plates by using a simple and accurate quasi-3D and 2D hyperbolic higher-order shear deformation theories. The properties of functionally graded materials (FGMs) are assumed to vary continuously through the thickness direction according to exponential law distribution (E-FGM). The kinematics of the present theories is modeled with an undetermined integral component and satisfies the free transverse shear stress conditions on the top and bottom surfaces of the plate; therefore, it does not require the shear correction factor. The fundamental governing differential equations and boundary conditions of exponentially graded plates are derived by employing the static version of principle of virtual work. Analytical solutions for bending of EG plates subjected to sinusoidal distributed load are obtained for simply supported boundary conditions using Navier'is solution procedure developed in the double Fourier trigonometric series. The results for the displacements and stresses of geometrically different EG plates are presented and compared with 3D exact solution and with other quasi-3D and 2D higher-order shear deformation theories to verify the accuracy of the present theory.

• Structural Engineering and Mechanics
• /
• v.90 no.4
• /
• pp.325-343
• /
• 2024

Triple-tower double-cable suspension bridges have increased confinement stiffness imposed by the main cable on the middle tower, which has bright application prospects. However, vertical bending and torsional vibrations of the double-cable and the girder are coupled in such bridges due to the hangers. In particular, the bending vibration of the towers in the longitudinal direction and torsional vibrations about the vertical axis influence the vertical bending and torsional vibrations of the stiffening girders, respectively. The conventional analytical algorithm for assessing the dynamic features of the suspension bridge is not directly applicable to this type of bridge. This study attempts to mitigate this problem by introducing an analytical algorithm for solving the triple-tower double-cable suspension bridge's natural frequencies and mode shapes. D'Alembert's principle is employed to construct the differential equations of the vertical bending and torsional vibrations of the stiffening girder continuum in each span. Vibrations of stiffening girders in each span are interrelated via the vibrations of the main cables and the bridge towers. On this basis, the natural frequencies and mode shapes are derived by separating variables. The proposed algorithm is then applied to an engineering example. The natural frequencies and mode shapes of vertical bending and torsional vibrations derived by the analytical algorithm agreed well with calculations via the finite element method. The fundamental frequency of vertical bending and first- and second-order torsion frequencies of double-cable suspension bridges are much higher than those of single-cable suspension bridges. The analytical algorithm has high computational efficiency and calculation accuracy, which can provide a reference for selecting appropriate structural parameters to meet the requirements of dynamics during the preliminary design.

• Korean Chemical Engineering Research
• /
• v.53 no.3
• /
• pp.295-301
• /
• 2015

A theoretical analysis was carried out to examine the concentration behavior of methane from a biogas using a polysulfone membrane. After the governing equations were derived for the cocurrent flow mode in a membrane module, the coupled nonlinear differential equations were numerically solved with the Compaq Visual Fortran 6.6 software. At the typical operating condition of mole fraction of 0.7 in a feed stream, the mole fraction of methane in the retentate increased to 0.76 while the normalized retentate flow rate to the feed flow rate decreased from 1 to 0.79. When either the mole fraction of methane in a feed increased or the pressure of the feed stream increased, the methane mole fraction in the retentate increased. On the other hand, it was found that as either the membrane area decreased or the ratio of the permeate pressure to the feed pressure increased, the methane mole fraction in the retentate decreased. In case that the stage cut increased, the methane mole fraction in the retentate increased while the recovery of methane slightly decreased.