• Steel and Composite Structures
• /
• v.12 no.1
• /
• pp.53-72
• /
• 2012

In this study, the physically-based failure models for matrix and fibers in compression and tension loading are introduced. For the 3D stress based fiber kinking model a modification is proposed for calculation of the fiber misalignment angle. All of these models are implemented into the finite element code by using the advantage of damage variable and the numerical results are discussed. To investigate the matrix failure model, purely in-plane transverse compression experiments are carried out on the specimens made by Glass/Epoxy to obtain the fracture surface angle and then a comparison is made with the calculated numerical results. Furthermore, shear failure of $({\pm}45)_s$ model is investigated and the obtained numerical results are discussed and compared with available experimental results. Some experiments are also carried out on the woven laminated composites to investigate the fracture pattern in the matrix failure mode and shown that the presented matrix failure model can be used for the woven composites. Finally, the obtained numerical results for stress based fiber kinking model and improved ones (strain based model) are discussed and compared with each other and with the available results. The results show that these models can predict the kink band angle approximately.

• Smart Structures and Systems
• /
• v.22 no.4
• /
• pp.459-468
• /
• 2018

In this paper, the effect of non-persistent joints was determined on the behavior of concrete specimens subjected to biaxial loading through numerical modeling using particle flow code in two dimensions (PFC2D). Firstly, a numerical model was calibrated by uniaxial, Brazilian and triaxial experimental results to ensure the conformity of the simulated numerical model's response. Secondly, sixteen rectangular models with dimension of 100 mm by 100 mm were developed. Each model contains two non-persistent joints with lengths of 40 mm and 20 mm, respectively. The angularity of the larger joint changes from $30^{\circ}$ to $90^{\circ}$. In each configuration, the small joint angularity changes from $0^{\circ}$ to $90^{\circ}$ in $30^{\circ}$ increments. All of the models were under confining stress of 1 MPa. By using of the biaxial test configuration, the failure process was visually observed. Discrete element simulations demonstrated that macro shear fractures in models are because of microscopic tensile breakage of a large number of bonded discs. The failure pattern in Rock Bridge is mostly affected by joint overlapping whereas the biaxial strength is closely related to the failure pattern.

• Journal of Korea Water Resources Association
• /
• v.34 no.2
• /
• pp.177-186
• /
• 2001

By using he total variation diminishing (TVD) condition, accurate and upwind based schemes are firstly introduced to develop numerical models free from nonphysical oscillations in the vicinity of large gradients. These models are then applied to both abruptly and smoothly varying initial conditions. By comparing computed predictions to analytical solutions, it is clearly shown that the first-order upwind scheme produces the numerical viscosity and the second-order Lax-Wendroff scheme produces the spurious oscillations. However, the TVD scheme gives the most reasonable results.

• Steel and Composite Structures
• /
• v.49 no.4
• /
• pp.431-440
• /
• 2023

The functionally graded (FG) porous plates are usually characterized by the non-symmetric elastic modulus distribution through the thickness so that the plate neutral surface does not coincide with the mid-surface. Nevertheless, the conventional analysis models were mostly based on the plate mid-surface so that the accuracy of resulting numerical results is questionable. In this context, this paper presents the neutral surface-based static and free vibration analysis of FG porous plates and investigates the differences between the mid- and neutral surface-based analysis models. The neutral surface-based numerical method is formulated using the (3,3,2) hierarchical model and approximated by the last introduced natural element method (NEM). The volume fractions of metal and ceramic are expressed by the power-law function and the cosine-type porosity distributions are considered. The proposed numerical method is demonstrated through the benchmark experiment, and the differences between two analysis models are parametrically investigated with respect to the thickness-wise material and porosity distributions. It is found from the numerical results that the difference cannot be negligible when the material and porosity distributions are remarkably biased in the thickness direction.

• The Journal of Engineering Geology
• /
• v.25 no.3
• /
• pp.369-376
• /
• 2015

Numerical models simulating groundwater flow are often used to estimate groundwater discharge into a tunnel. In designing numerical models, the grid size should be carefully considered to ensure that groundwater discharge is accurately predicted. However, several recent studies have employed various grid sizes without providing an adequate explanation for their choice. This paper suggests the optimal grid size based on a comparison of numerical models with analytical solutions. Discrepancies between numerical and analytical solutions result from the effect of model boundaries as well as the grid size. By nullifying boundary effects, the errors solely associated with the grid size could be analyzed. The optimal grid size yielding accurate numerical solutions was thus derived. The suggested relationship between tunnel radius and optimal grid size is analogous to the relationship between the equivalent well block radius and grid size.

• Computers and Concrete
• /
• v.7 no.3
• /
• pp.203-220
• /
• 2010

In the present study, exterior beam-column sub-assemblage from a regular reinforced concrete (RC) building has been considered. Two different types of beam-column sub-assemblages from existing RC building have been considered, i.e., gravity load designed ('GLD'), and seismically designed but without any ductile detailing ('NonDuctile'). Hence, both the cases represent the under-designed structure at different time frame span before the introduction of ductile detailing. For designing 'NonDuctile' structure, Eurocode and Indian Standard were considered. Non-linear finite element (FE) program has been employed for analysing the sub-assemblages under cyclic loading. FE models were developed using quadratic concrete brick elements with embedded truss elements to represent reinforcements. It has been found that the results obtained from the numerical analysis are well corroborated with that of experimental results. Using the validated numerical models, it was proposed to correlate the energy dissipation from numerical analysis to that from experimental analysis. Numerical models would be helpful in practice to evaluate the seismic performance of the critical sub-assemblages prior to design decisions. Further, using the numerical studies, performance of the sub-assemblages with variation of axial load ratios (ratio is defined by applied axial load divided by axial strength) has been studied since many researchers have brought out inconsistent observations on role of axial load in changing strength and energy dissipation under cyclic load.

• Coupled systems mechanics
• /
• v.7 no.1
• /
• pp.95-109
• /
• 2018

The relevance of turbulent mixing in estuarine numerical models for stratified two-layer shallow water flows is analysed in this paper. A one-dimensional numerical model was developed for this purpose by extending an immiscible two-layer model with an additional source term, which accounts for turbulent mixing effects, namely the entrainment of fluid from the lower to the upper layer. The entrainment rate is quantified by an empirical equation as a function of the bulk Richardson number. A finite volume method based on an approximated Roe solver was used to solve the governing coupled system of partial differential equations. A comparison of numerical results with and without entrainment is presented to illustrate the influence of entrainment on both the salt-water intrusion length and lower layer dynamics. Furthermore, one example is given to demonstrate how entrainment terms may help to stabilize the numerical scheme and prevent a possible loss of hyperbolicity. Finally, the model with entrainment is validated by comparing the numerical results to field measurements.

• Atmosphere
• /
• v.14 no.4
• /
• pp.19-23
• /
• 2004

Progresses in numerical prediction of weather and climate have been in parallel with those of computing resources, especially the development of supercomputers. Advanced techniques in numerical modeling, computational schemes, and data assimilation cloud not have been practically achieved without the aid of supercomputers. With such techniques and computing powers, the accuracy of numerical forecasts has been tremendously improved. Supercomputers are also indispensible in constructing and executing the synthetic Earth system models. In this study, a brief overview on numerical weather / climate prediction, Earth system modeling, and the values of supercomputing is provided.

• Journal of Industrial Technology
• /
• v.36
• /
• pp.23-31
• /
• 2016

This paper presents EMQ models in which some proportion of defective items are produced and some of them are converted to good items through rework process and items not converted are disposed. Numerical models are developed for three cases of disposal and optimal solution of each model is derived. In the first model, if a defective item is found during the production process, only re-workable items are stored and reworked after normal production is finished. Not re-workable items are disposed immediately during normal production. The second model deals with the case where all defective items are stored and items to be disposed are determined in rework process. In the third model, an additional inspection process exists before rework to determine rework or disposal. Numerical examples are presented to validate the proposed models.

• Korea-Australia Rheology Journal
• /
• v.14 no.4
• /
• pp.161-174
• /
• 2002

A new technique for numerical calculation of viscoelastic flow based on the combination of Neural Net-works (NN) and Brownian Dynamics simulation or Stochastic Simulation Technique (SST) is presented in this paper. This method uses a "universal approximator" based on neural network methodology in combination with the kinetic theory of polymeric liquid in which the stress is computed from the molecular configuration rather than from closed form constitutive equations. Thus the new method obviates not only the need for a rheological constitutive equation to describe the fluid (as in the original Calculation Of Non-Newtonian Flows: Finite Elements St Stochastic Simulation Techniques (CONNFFESSIT) idea) but also any kind of finite element-type discretisation of the domain and its boundary for numerical solution of the governing PDE's. As an illustration of the method, the time development of the planar Couette flow is studied for two molecular kinetic models with finite extensibility, namely the Finitely Extensible Nonlinear Elastic (FENE) and FENE-Peterlin (FENE-P) models.P) models.