• Journal of the Korea Furniture Society
• /
• v.21 no.1
• /
• pp.26-39
• /
• 2010

This is a leading study to replace the structural analysis methodology on the specific traditional joint by a numerical analysis. Tests were carried out to test the compressive methodologies with the numerical results. The Japanese larch was used as a sample. The Orthotropic property of wood was specifically considered for the finite element numerical analysis. Linear numerical analysis and non-linear numerical analysis for the BEAM element and the two SOLID elements of ANSYS were used to analyze the compressive performance. In addition, more finely divided elements were used to raise the accuracy of the numerical result. Finally, the statistically significant differences were tested between that of the analytical and numerical results. It could be concluded that the SOLID 64 element shows the most optimum result when the non-linear analysis with the more finely divided element was used. However, finely dividing of the element is a considerable time consuming process, and it is quite difficult to raise the accuracy of the non-linear numerical analysis. Therefore, if considering the vertical displacement to be of the only interest, the BEAM element is more efficient than the SOLID element because the BEAM element is reflected as a simple line, which is less time consuming and difficult in dividing the elements. But, the BEAM element cannot accurately model the knot as a strength defect factor which is an important property in the orthotropic property of wood. Therefore, the SOLID element should be used to model the strength defect factor, knot, as it can be efficiently applied on the structural size flexure member which could be more strongly effected by the knot. In addition, it is useful at times when the failure types of members are to be more closely investigated, as the SOLID element is able to examine the local stress distribution of the member. The conclusion drawn by this study is of the good concordance between analytical results and numerical results of compressive wood members, but how orthotropic properties should only be considered. The numerical analysis on the specific Korean traditional joints will be based on the current study results.

• Structural Engineering and Mechanics
• /
• v.72 no.5
• /
• pp.617-629
• /
• 2019

Based on the finite element method of traditional straight Euler-Bernoulli beams and the coupled relations between linear displacement and angular displacement of a pre-twisted Euler-Bernoulli beam, the shape functions and stiffness matrix are deduced. Firstly, the stiffness of pre-twisted Euler-Bernoulli beam is developed based on the traditional straight Euler-Bernoulli beam. Then, a new finite element model is proposed based on the displacement general solution of a pre-twisted Euler-Bernoulli beam. Finally, comparison analyses are made among the proposed Euler-Bernoulli model, the new numerical model based on displacement general solution and the ANSYS solution by Beam188 element based on infinite approach. The results show that developed numerical models are available for the pre-twisted Euler-Bernoulli beam, and which provide more accurate finite element model for the numerical analysis. The effects of pre-twisted angle and flexural stiffness ratio on the mechanical property are investigated.

• Transactions of the Korean Society of Automotive Engineers
• /
• v.14 no.2
• /
• pp.174-183
• /
• 2006

The resistance spot-welded region in most current finite element crash models is characterized as a rigid beam at the location of the welded spot. The region is modeled to fail with a failure criterion which is a function of the axial and shear load at the rigid beam. The role of this rigid beam is simply to transfer the load across the welded components. The calculation of the load acting on the rigid beam is important to evaluate the failure of the spot-weld. In this paper, numerical simulation is carried out to evaluate the calculation of the load at the rigid beam. The load calculated from the precise finite element model of the spot-welded region considering the residual stress due to the thermal history during the spot welding procedure is regarded as the reference value and the value of the load is compared with the one obtained from the spot-welded model using the rigid beam with respect to the element size, the element shape and the number of imposed constraints. Analysis results demonstrate that the load acting on the spot-welded element is correctly calculated by the change of the element shape around the welded region and the location of welded constrains. The results provide a guideline for an accurate finite element modeling of the spot-welded region in the crash analysis of vehicles.

• International Journal of Automotive Technology
• /
• v.7 no.3
• /
• pp.329-336
• /
• 2006

The resistance spot-welded region in most current finite element crash models is characterized as a rigid beam at the location of the welded spot. The region is modeled to fail with a failure criterion which is a function of the axial and shear load at the rigid beam. The calculation of the load acting on the rigid beam is important to evaluate the failure of the spot-weld. In this paper, numerical simulation is carried out to evaluate the calculation of the load at the rigid beam. At first, the load on the spot-welded region is calculated with the precise finite element model considering the residual stress due to the thermal history during the spot welding procedure. And then, the load is compared with the one obtained from the model used in the crash analysis with respect to the element size, the element shape and the number of imposed constraints. Analysis results demonstrate that the load acting on the spot-welded element is correctly calculated by the change of the element shape around the welded region and the location of welded constrains. The results provide a guideline for an accurate finite element modeling of the spot-welded region in the crash analysis of vehicles.

• Journal of Korean Association for Spatial Structures
• /
• v.3 no.2 s.8
• /
• pp.85-90
• /
• 2003

If we use a third order approximation for the displacement function of beam element in finite element methods, finite element solutions of beams yield nodal displacement values matching to beam theory results to have no connection with the number increasing of elements of beams. It is assumed that, as the member displacement value at beam nodes are correct, the calculation procedure of beam element stiffness matrix have no numerical errors. A the member forces are calculated by the equations of $\frac{-M}{EI}=\frac{{d^2}{\omega}}{dx^2}\;and\;\frac{dM}{dx}=V$, the member forces at nodes of beams have errors in a moment and a shear magnitudes in the case of smaller number of element. The nodal displacement value of plate subject to the lateral load converge to the exact values according to the increase of the number of the element. So it is assumed that the procedures of plate element stiffness matrix calculations has a error in the fundamental assumptions. The beam methods for the high accuracy ratio solution Is also applied to the plate analysis. The method of reducing a error ratio of member forces and element stiffness matrix in the finite element methods is studied. Results of study were as follows. 1. The matrixes of EI[B] and [K] in the equations of M(x)=EI[B]{q} and M(x) = [K]{q}+{Q} of beams are same. 2. The equations of $\frac{-M}{EI}=\frac{{d^2}{\omega}}{dx^2}\;and\;\frac{dM}{dx}=V$ for the member forces have a error ratio in a finite element method of uniformly loaded structures, so equilibrium node loads {Q} must be substituted in the equation of member forces as the numerical examples of this paper revealed.

• Journal of Advanced Navigation Technology
• /
• v.16 no.2
• /
• pp.367-374
• /
• 2012

A dynamic simulation and test of frame with thin walled closed section beams considering warping conditions have been performed. When a beam is subjected under torsional moment, the cross section will deform an warping as well as twist. For some thin-walled sections warping will be large, and accompanying warping restraint will induce axial and shear stresses and reduce the twist of beam which stiffens the beam in torsion. This paper presents that an warping restraint factor in finite element model effects the behavior of beam deformation and dynamic mode shape. The computer modelling of frame is discussed in linear beam element model and linear thin shell element model, also presents a correlation between computer predicted and actual experimental results for static deflection, natural frequencies and mode shapes of frame.

• Proceedings of the Korea Committee for Ocean Resources and Engineering Conference
• /
• 2004.11a
• /
• pp.271-276
• /
• 2004

In this paper. dynamic behavior of the cracked beam under a moving mass is presented using the finite element method (FEM). Model accuracy is improved with the following consideration: (1) FE model with Timoshenko beam element (2) Additional flexibility matrix due to crack presence (3) Interaction forces between the moving mass and supported beam. The Timoshenko bean model with a two-node finite element is constructed based on Guyan condensation that leads to the results of classical formulations. but in a simple and systematic manner. The cracked section is represented by local flexibility matrix connecting two unchanged beam segments and the crack as modeled a massless rotational spring. The inertia force due to the moving mass is also involved with gravity force equivalent to a moving load. The numerical tests for various mass levels. crack sizes. locations and boundary conditions were performed.

• Structural Engineering and Mechanics
• /
• v.39 no.1
• /
• pp.21-46
• /
• 2011

The natural frequencies of continuous systems depend on the governing partial differential equation and can be numerically estimated using the finite element method. The accuracy and convergence of the finite element method depends on the choice of basis functions. A basis function will generally perform better if it is closely linked to the problem physics. The stiffness matrix is the same for either static or dynamic loading, hence the basis function can be chosen such that it satisfies the static part of the governing differential equation. However, in the case of a rotating beam, an exact closed form solution for the static part of the governing differential equation is not known. In this paper, we try to find an approximate solution for the static part of the governing differential equation for an uniform rotating beam. The error resulting from the approximation is minimized to generate relations between the constants assumed in the solution. This new function is used as a basis function which gives rise to shape functions which depend on position of the element in the beam, material, geometric properties and rotational speed of the beam. The results of finite element analysis with the new basis functions are verified with published literature for uniform and tapered rotating beams under different boundary conditions. Numerical results clearly show the advantage of the current approach at high rotation speeds with a reduction of 10 to 33% in the degrees of freedom required for convergence of the first five modes to four decimal places for an uniform rotating cantilever beam.

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

This paper describes the time rate of change of dynamic response of a cantilever beam inserted with a damping element, such as bonding, which is excited under a general force at various locations. A sensitivity analysis was performed in a finite element model to show that two types of second-order algebraic governing equations were used to predict the rate of change of dynamic displacement: one is related to the modal coordinate linked to a physical coordinate, and the other to the design parameter of the time rate of change of displacement. The sensitivity differential equation formulation includes more complicated terms compared with that of the undamped cantilever beam. The sensitivities of the dynamic response were observed by changing the location of the excitation force, displacement extraction, and cross-sectional area of the beam. The analytical results obtained by this suggested theory showed a relatively good agreement when compared with those obtained using the commercial finite element program. The suggested analysis procedure enables the prediction of the response sensitivity for any finite element model of the dynamic system.

• Proceedings of the Korea Concrete Institute Conference
• /
• 2004.05a
• /
• pp.510-513
• /
• 2004

The objectives of this research are to evaluate the effect of the compressive strength of concrete, reinforcing bar size, spacing of column transverse bars related to the concrete confinement effects on anchorage bond strength and bond behavior of beam-column joints subjected to cyclic loading and to predict the bond behavior of beam-column joints according to the variables by Finite Element Analysis appling the interface element between concrete and reinforced bar surface in a three-dimensional configuration. This paper shows that to verify the results by three-dimensional nonlinear finite element analysis appling a interface element, the test results that were already conducted are compared with analytic results. The behavior of bond and anchorage of beam bar is expressed by a local bond stress-slip relationship and the failure mode of bond is predicted by principal stress contour.