• Structural Engineering and Mechanics
• /
• v.48 no.1
• /
• pp.41-56
• /
• 2013

Initiation and growth of delamination is a great concern of designers of composite structures. Interface elements with de-cohesive constitutive law in the content of continuum damage mechanics can be used to predict initiation and growth of delamination in single and mixed mode conditions. In this paper, an interface element based on the cohesive zone method has been developed to simulate delaminatoin growth of post-buckled laminate under fatigue loading. The model was programmed as the user element and user material by the "User Programmable Features" in ANSYS finite element software. The interface element is a three-dimensional 20 node brick with small thickness. Because of mixed-mode condition of stress field at the delamination-front of post-buckled laminates, a mixed-mode bilinear constitutive law has been used as user material in this model. The constitutive law of interface element has been verified by modelling of a single element. A composite laminate with initial delamination under quasi-static compressive Loading available from literature has been remodeled with the present approach. Moreover, it will be shown that, the closer the delamination to the free surface of laminate, the slower the delamination growth under compressive fatigue loading. The effects of laminate configuration on delamination growth are also investigated.

• Composites Research
• /
• v.12 no.6
• /
• pp.65-73
• /
• 1999

Crippling stress and failure behavior of Z-section graphite/epoxy composite laminated stringers are investigated by the nonlinear finite element method. Stringers are idealized using 9-node laminated shell element. The complete unloading model is introduced into the finite element method for the progressive failure analysis. A modified Riks method is used to trace the post-failure equilibrium path after local buckling. Finite element results are validated with previous experimental results. The results show that the most important parameter affecting the crippling stress of Z-section stringers is the flange width. In terms of stacking sequence. the highest cripping stress is found at the stringer with $[{\pm}45/0/90]s$ lamination.

• Advances in Computational Design
• /
• v.5 no.3
• /
• pp.305-321
• /
• 2020

In this paper is presented the solution method for three-dimensional problem of transversely isotropic body's elastoplastic deformation by the finite element method (FEM). The process of problem solution consists of: determining the effective parameters of a transversely isotropic medium; construction of the finite element mesh of the body configuration, including the determination of the local minimum value of the tape width of non-zero coefficients of equation systems by using of front method; constructing of the stiffness matrix coefficients and load vector node components of the equation for an individual finite element's state according to the theory of small elastoplastic deformations for a transversely isotropic medium; the formation of a resolving symmetric-tape system of equations by summing of all state equations coefficients summing of all finite elements; solution of the system of symmetric-tape equations systems by means of the square root method; calculation of the body's elastoplastic stress-strain state by performing the iterative process of the initial stress method. For each problem solution stage, effective computational algorithms have been developed that reduce computational operations number by modifying existing solution methods and taking into account the matrix coefficients structure. As an example it is given, the problem solution of fibrous composite straining in the form of a rectangle with a system of circular holes.

• Ocean Systems Engineering
• /
• v.8 no.1
• /
• pp.57-77
• /
• 2018

The present paper deals with the finite element analysis of water tanks with rigid baffle. Fluid is discretized by two dimensional eight-node isoparametric elements and the governing equation is simulated by pressure based formulation to reduce the degrees of freedom in the domain. Both free vibration and force vibration analysis are carried out for different sizes and positions of block at tank bottom. The fundamental frequency depends on block height and it reduces with the increase of block height. The variation of hydrodynamic pressure on tank walls not only depends of the exciting frequency but also on the size and position of rigid block at tank bottom. The hydrodynamic pressure has higher value when the exciting frequency is equal and lower than the fundamental frequency of the water in the tank. Similarly, the hydrodynamic pressure increases with the increase of width of the block for all exciting frequencies when the block is at the centre of tank. The left and right walls of tank have experienced different hydrodynamic pressure when the block is placed at off-centre. However, the increase in hydrodynamic pressure on nearest tank wall becomes insignificant after a certain value of the distance between the wall and the rigid block.

• Journal of Korean Society of Industrial and Systems Engineering
• /
• v.26 no.3
• /
• pp.58-66
• /
• 2003

The cut tree approach of Montreuil and Ratliff [16] and eigenvector approach [10] are used to automatically draw a feasible facility layout with aisle structure. The department arrangement can minimize an aisle distance criterion considering door locations and dimension constraints. The aisle distance is measured by the door to door distance between departments. An eigenvector and cut tree approaches [1] are implemented based on the branch and bound technique in Kim et al. [2] in order to obtain feasible layouts. Then, the algorithm to fix the door location of each department is developed. After the door locations are determined, the factory layout is evaluated in terms of aisle distance. The aisle structure is obtained by expanding the original layout. The solution is kept until we will find better factory layout. The proposed approach based on the branch and bound technique, in theory, will provide the optimal solution. If the runs are time and/or node limited, the proposed method is a strong heuristic. The technique is made further practical by the fact that the solution is constrained such that the rectangular shape dimensions length(l) and width(w) are fixed and a perfect fit is generated if a fit is possible.

• Steel and Composite Structures
• /
• v.29 no.6
• /
• pp.749-758
• /
• 2018

This article derived a hybrid coupling technique using the higher-order displacement polynomial and three soft computing techniques (teaching learning-based optimization, particle swarm optimization, and artificial bee colony) to predict the optimal stacking sequence of the layered structure and the corresponding frequency values. The higher-order displacement kinematics is adopted for the mathematical model derivation considering the necessary stress and stain continuity and the elimination of shear correction factor. A nine noded isoparametric Lagrangian element (eighty-one degrees of freedom at each node) is engaged for the discretisation and the desired model equation derived via the classical Hamilton's principle. Subsequently, three soft computing techniques are employed to predict the maximum natural frequency values corresponding to their optimum layer sequences via a suitable home-made computer code. The finite element convergence rate including the optimal solution stability is established through the iterative solutions. Further, the predicted optimal stacking sequence including the accuracy of the frequency values are verified with adequate comparison studies. Lastly, the derived hybrid models are explored further to by solving different numerical examples for the combined structural parameters (length to width ratio, length to thickness ratio and orthotropicity on frequency and layer-sequence) and the implicit behavior discuss in details.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.20 no.5
• /
• pp.8-16
• /
• 2021

This study develops a shaft analysis model for the optimization of the power transmission system. The finite element method was used for the shaft analysis model. The shaft and gear were assumed Timoshenko beams. Strength was evaluated according to DIN 743, and gear misalignment was calculated through ISO 6336 and the coordinate system rotation. The analysis software for a power transmission system was developed using Visual Studio 2019. The analysis results of the developed program were compared with those of commercial software (MASTA, KISSsoft, and Romax). We confirmed that the force, deformation, and safety factors at each node were the same as those of the commercial software. The absolute value of the gear misalignment of the developed program and commercial software was different. However, the gear misalignment tended to increase with increasing the displacement in the tooth width direction.

• Restorative Dentistry and Endodontics
• /
• v.19 no.2
• /
• pp.345-371
• /
• 1994

Restorative procedures can lead to weakening tooth due to reduction and alteraton of tooth structure. It is essential to prevent fractures to conserve tooth. Among the several parameters in cavity designs, cavity isthmus and depth are very important. In this study, MO amalgam cavity was prepared on maxillary first premolar. Three dimensional. finite element models were made by serial photographic method and cavity depth(1.7mm, 2.4mm) and isthmus (11 4, 1/3, 1/2 of intercuspal distance) were varied. linear, eight and six-nodal, isoparametric brick elements were used for the three dimensional finite element model. The periodontal ligament and alveolar bone surrounding the tooth were excluded in these models. Three types model(B, G and R model) were developed. B model was assumed perfect bonding between the restoration and cavity wall. Both compressive and tensile forces were distributed directly to the adjacent regions. G model(Gap Distance: 0.000001mm) was assumed the possibility of play at the interface simulated the lack of real bonding between the amalgam and cavity wall (enamel and dentin). When compression occurred along the interface, the forces were transferred to the adjacent regions. However, tensile forces perpendicular to the interface were excluded. R model was assumed non-connection between the restoration and cavity wall. No force was transferred to the adjacent regions. A load of 500N was applied vertically at the first node from the lingual slope of the buccal cusp tip. This study analysed the displacement, von Mises stress, 1 and 2 direction normal stress and strain with FEM software ABAQUS Version 5.2 and hardware IRIS 4D/310 VGX Work-station. The results were as follows: 1. G model showed stress and strain patterns between Band R model. 2. B model and G model showed the bending phenomenon in the displacement. 3. R model showed the greatest amount of the displacement of the buccal cusp followed by G and B model in descending order. G model showed the greatest amount of the displacement of the lingual cusp followed by B and R model in descending order. 4. B model showed no change of the displacement as increasing depth and width of the cavity. G and R model showed greater displacement of the buccal cusp as increasing depth and width of the cavity, but no change in the displacement of the lingual cusp. 5. As increasing of the width of the cavity, stress and strain were not changed in B model. Stress and strain were increased on the distal marginal ridge and buccopulpal line angle in G and R model. The possibility of the tooth fracture was increased. 6. As increasing of the depth of the cavity, stress and strain were not changed in B and G model. Stress and strain were increased on the distal marginal ridge and buccopulpal line angle in R model. The possibility of the tooth fracture was increased.

• Journal of the Korea Institute of Information and Communication Engineering
• /
• v.20 no.2
• /
• pp.306-316
• /
• 2016

In this paper, an NMOS-diode eFuse OTP (One-Time Programmable) memory cell is proposed using a parasitic junction diode formed between a PW (P-Well), a body of an isolated NMOS (N-channel MOSFET) transistor with the small channel width, and an n+ diffusion, a source node, in a DNW (Deep N-Well) instead of an NMOS transistor with the big channel width as a program select device. Blowing of the proposed cell is done through the parasitic junction formed in the NMOS transistor in the program mode. Sensing failures of '0' data are removed because of removed contact voltage drop of a diode since a NMOS transistor is used instead of the junction diode in the read mode. In addition, a problem of being blown for a non-blown eFuse from a read current through the corresponding eFuse OTP cell is solved by limiting the read current to less than $100{\mu}A$ since a voltage is transferred to BL by using an NMOS transistor with the small channel width in the read mode.

• Journal of Radiation Industry
• /
• v.6 no.1
• /
• pp.31-40
• /
• 2012

There are two methods to fabricate the readout electronic to a large-area CMOS image sensor (LACIS). One is to design and manufacture the sensor part and signal processing electronics in a single chip and the other is to integrate both parts with bump bonding or wire bonding after manufacturing both parts separately. The latter method has an advantage of the high yield because the optimized and specialized fabrication process can be chosen in designing and manufacturing each part. In this paper, LACIS chip, that is optimized design for the latter method of fabrication, is presented. The LACIS chip consists of a 3-TR pixel photodiode array, row driver (or called as a gate driver) circuit, and bonding pads to the external readout ICs. Among 4 types of the photodiode structure available in a standard CMOS process, $N_{photo}/P_{epi}$ type photodiode showed the highest quantum efficiency in the simulation study, though it requires one additional mask to control the doping concentration of $N_{photo}$ layer. The optimized channel widths and lengths of 3 pixel transistors are also determined by simulation. The select transistor is not significantly affected by channel length and width. But source follower transistor is strongly influenced by length and width. In row driver, to reduce signal time delay by high capacitance at output node, three stage inverter drivers are used. And channel width of the inverter driver increases gradually in each step. The sensor has very long metal wire that is about 170 mm. The repeater consisted of inverters is applied proper amount of pixel rows. It can help to reduce the long metal-line delay.