• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.24 no.6 s.177
• /
• pp.1565-1574
• /
• 2000

Surface edge crack subjected to contact stresses is analysed. A punch with corner radii is considered to press the semi-infinite plane. Partial slip problem is solved when a shear force is applied to the punch. Dislocation density function method is used to solve the present mixed mode crack problem. The crack length of positive K1 is examined, which is affected by the ratio of the flat portion to the total width of the punch. Surface traction during one cycle of the shear force is evaluated to simulate the fretting condition. The compliance change of the contact surface is also investigated during the shear cycle. It is found that the crack grows during only a part of the cycle, which may be termed as effective period of crack growing. A design method for restraining the fretting failure is discussed, from which recommendable geometry of the punch is suggested.

• Steel and Composite Structures
• /
• v.26 no.4
• /
• pp.387-405
• /
• 2018

A novel higher shear deformation theory (HSDT) is proposed for the bending, buckling and free vibration investigations of isotropic and functionally graded (FG) sandwich plates. It contains only four variables, which is even less than the first shear deformation theory (FSDT) and the conventional HSDTs. The model accounts for a parabolic variation of transverse shear stress, respects the traction free boundary conditions and contrary to the conventional HSDTs, the present one presents a novel displacement field which incorporates undetermined integral terms. Equations of motion determined in this work are applied for three types of FG structures: FG plates, sandwich plates with FG core and sandwich plates with FG faces. Analytical solutions are given to predict the transverse displacements, stresses, critical buckling forces and natural frequencies of simply supported plates and a comparison study is carried out to demonstrate the accuracy of the proposed model.

• Computational Structural Engineering
• /
• v.9 no.1
• /
• pp.93-99
• /
• 1996

In the particular situations where the crack is terminated at an interface of two materials, the order of stress singularity depends on the elastic constants which specify the properties of two materials. A multidomain boundary element technique is used to solve a crack normal to bimaterial interface. A correct order of shape function is used for displacement by using the isoparametric elements by shifting adequately the side nodes adjacent to this crack tip. A shape function containing the same order of singularity as that in the interface crack is also used for the interpolation of traction. Numerical testing of a binaterial with a crack normal to the interface is carried out with three-node elements. The results obtained are compared with the previous solutions.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.3
• /
• pp.223-229
• /
• 2009

Most of structural analyses are concerned with the deformation and stress in a body subjected to external loads. In many fields, however, the interpretation of inverse problems is needed to determine surface tractions or internal stresses from measured displacements. In this study, the inverse processes by using the boundary element method are formulated for the evaluation of boundary tractions from displacements measured on a remote surface. Small errors in measured displacements often result in a substantial loss of accuracy of an inverse system. Numerical results show that the error in reconstructed tractions by using the inverse boundary element methods is sensitive to measurement location and noise.

• Journal of Mechanical Science and Technology
• /
• v.14 no.2
• /
• pp.197-206
• /
• 2000

A boundary integral equation method in the shape design sensitivity analysis is developed for the elasticity problems with axisymmetric non-homogeneous bodies. Functionals involving displacements and tractions at the zonal interface are considered. Sensitivity formula in terms of the interface shape variation is then derived by taking derivative of the boundary integral identity. Adjoint problem is defined such that displacement and traction discontinuity is imposed at the interface. Analytic example for a compound cylinder is taken to show the validity of the derived sensitivity formula. In the numerical implementation, solutions at the interface for the primal and adjoint system are used for the sensitivity. While the BEM is a natural tool for the solution, more generalization should be made since it should handle the jump conditions at the interface. Accuracy of the sensitivity is evaluated numerically by the same compound cylinder problem. The endosseous implant-bone interface problem is considered next as a practical application, in which the stress value is of great importance for successful osseointegration at the interface. As a preliminary step, a simple model with tapered cylinder is considered in this paper. Numerical accuracy is shown to be excellent which promises that the method can be used as an efficient and reliable tool in the optimization procedure for the implant design. Though only the axisymmetric problem is considered here, the method can be applied to general elasticity problems having interface.

• Journal of the Korean Society for Precision Engineering
• /
• v.33 no.1
• /
• pp.69-75
• /
• 2016

Polydimethylsiloxane (PDMS) is used as a scaffold for cell culture. Because both the stress and strain acting on the substrate and the hemodynamic environment are important for studying mechano-transduction of cellular function, the traction force of the surface of a substrate has been measured using fluorescence images of particle distribution. In this study, deformation of the cross-sectional plane of a PDMS block was measured by correlating particle image distributions to validate the particle image strain measurement technique. Deformation was induced by a cone indentor and a shearing parallel plate. Measured deformations from particle image distributions were in agreement with the results of a computational structure analysis using the finite-element method. This study demonstrates that the particle image correlation method facilitates measurement of deformation of a polymer scaffold in the cross-sectional plane.

• Computers and Concrete
• /
• v.14 no.1
• /
• pp.91-107
• /
• 2014

Using the principle of damage mechanics, zero-thickness pore pressure cohesive elements (PPCE) are used to simulate the casing-cement interface (CCI) and cement-rock interface (CRI). The traction-separation law describes the emergence and propagation of the PPCE. Mohr-coulomb criteria determines the elastic and plastic condition of cement sheath and rock. The finite element model (FEM) of delamination fractures emergence and propagation along the casing-cement-rock (CCR) interfaces during hydraulic fracturing is established, and the emergence and propagation of fractures along the wellbore axial and circumferential direction are simulated. Regadless of the perforation angle (the angle between the perforation and the max. horizontal principle stress), mirco-annulus will be produced alonge the wellbore circumferential direction when the cementation strength of the CCI and the CRI is less than the rock tensile strength; the delamination fractures are hard to propagate along the horizontal wellbore axial direction; emergence and propagation of delamination fractures are most likely produced on the shallow formation when the in-situ stresses are lower; the failure mode of cement sheath in the deep well is mainly interfaces seperation and body damange caused by cement expansion and contraction, or pressure testing and well shut-in operations.

• Structural Engineering and Mechanics
• /
• v.64 no.5
• /
• pp.591-610
• /
• 2017

This paper is concerned with the static analysis of variable thickness of two directional functionally graded porous materials (FGPM) circular plate resting on a gradient hybrid foundation (Horvath-Colasanti type) with friction force and subjected to compound mechanical loads (e.g., transverse, in-plane shear traction and concentrated force at the center of the plate).The governing state equations are derived in terms of displacements based on the 3D theory of elasticity, assuming the elastic coefficients of the plate material except the Poisson's ratio varying continuously throughout the thickness and radial directions according to an exponential function. These equations are solved semi-analytically by employing the state space method (SSM) and one-dimensional differential quadrature (DQ) rule to obtain the displacements and stress components of the FGPM plate. The effect of concentrated force at the center of the plate is approximated with the shear force, uniformly distributed over the inner boundary of a FGPM annular plate. In addition to verification study and convergence analysis, numerical results are displayed to show the effect of material heterogeneity indices, foundation stiffness coefficients, foundation gradient indices, loads ratio, thickness to radius ratio, compressibility, porosity and friction coefficient of the foundation on the static behavior of the plate. Finally, the responses of FG and FG porous material circular plates to compound mechanical loads are compared.

• Korea-Australia Rheology Journal
• /
• v.13 no.1
• /
• pp.37-45
• /
• 2001

Numerical simulation of coextrusion process of viscoelastic fluids within a die has been carried out. In the coextrusion process velocity profile at the outflow boundary is not known a priori, which makes it difficult to impose the proper boundary condition at the outflow boundary. This difficulty has been avoided by using the open boundary condition (OBC) method. In this study, elastic viscous stress splitting (EVSS) formulation with streamline upwind (SU) method has been used in the finite element method. In order to test the validity of the OBC method, comparison between the results of fully developed condition at the outlet and those of OBC has been made for a Newtonian fluid. In the case of upper convected Maxwell (UCM) fluid, the effect of outflow boundary condition on the interface position has been investigated by using two meshes having different downstream lengths. In both cases, the results with the OBC method showed reasonable interface shape. In particular, for the UCM fluid the interface shape calculated with OBC was independent of the downstream length, while the results with the zero traction condition showed oscillation of interface position close to the outlet. Viscosity difference was found to be more important than elasticity difference in determining the final interface position. However, the overshoot of interface position near the con-fluent point increased with elasticity.

• Computers and Concrete
• /
• v.1 no.1
• /
• pp.47-60
• /
• 2004

The influence of the fracture process zone (FPZ) on the fracture properties is one of the hottest topics in the field of fracture mechanics for cementitious materials. Within the FPZ in front of a traction free crack, cohesive forces are distributed in accordance with the softening stress-separation constitutive relation of the material. Therefore, further crack propagation necessitates energy dissipation, which is the work done by the cohesive forces. In this paper $g_f$, the local fracture energy characterizing the energy consumption due to the cohesive forces, is discussed. The computational expression of $g_f$ in the FPZ can be obtained for any stage during the material fracture process regarding the variation of FPZ, whether in terms of its length or width. $G_{fa}$, the average energy consumption along the crack extension region, has also been computed and discussed in this paper. The experimental results obtained from the wedge splitting tests on specimens with different initial notch ratios are employed to investigate the property of the local fracture energy $g_f$ and the average value $G_{fa}$ over the crack extension length. These results can be used to indicate the influence of the FPZ. Additionally, changes in the length of the FPZ during the fracture process are also studied.