Implementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis
|
|
- Buddy Daniels
- 6 years ago
- Views:
Transcription
1 CST0 191 October, 011, Krabi Imlementation and Validation of Finite Volume C++ Codes for Plane Stress Analysis Chakrit Suvanjumrat and Ekachai Chaichanasiri* Deartment of Mechanical Engineering, Faculty of Engineering, Mahidol University, Nakornathom, Thailand *Corresonding Author: Telehone ext. 6409, Fax ext. 649 Abstract The interaction between solid and fluid is an imortant and interesting subject for the resent. In general, to simulate the interaction between solid and fluid, finite element method (FEM) and finite volume method (FVM) are used for analysis of solid and fluid resectively. The transformation of results between FVM and FEM has a disadvantage since it takes time to transfer data between FEM and FVM. The individual using FVM for analysis of solid and fluid can avoid transferring data and it is more stable when simulating comlicated roblems. This research resented the FVM codes using C++ to erform a lane stress analysis. The code was imlemented using Oen Source Software that was OenFOAM. The written code was validated with lane stress roblems as the test cases, which comrised the tension of thin rectangular late with and without a circular hole at the center. The results of FVM were comared with the analytical solutions. From the results, the average difference between the FVM results and analytical solutions was less than 1.68%. These results suggested the otential of using FVM for stress analysis and it will be used for analyzing solidfluid interaction. Keywords: Finite volume method, Finite element method, C++, Stress analysis, OenFOAM 1. Introduction The interaction between solid and fluid is an interesting subject for the resent. For examle, a sloshing of liquid in vehicles [1,], an imact of the Tsunami wave on buildings [3]. In general, to simulate the interaction between solid and fluid, finite element method (FEM) and finite volume method (FVM) are used. The FEM is used for analysis of solid and the FVM is used for analysis of fluid. The transformation of results FEM has a disadvantage since it takes time to transfer data between FEM and FVM. The FVM is inherently good at treating comlicated, couled and nonlinear differential equations, widely resent in fluid flows. The individual using FVM for analysis of solid and fluid can avoid transferring data and it is more stable when simulating comlicated roblems. This research resented the FVM codes using C++ to erform a lane stress analysis. The code was imlemented using Oen Source Software that was OenFOAM.
2 . Materials and Methods.1 Governing Equation The Cauchy s first law of motion for solid body was written in the differential equation, ( D) t CST0 191 October, 011, Krabi (1) where D was the dislacement vector, was the stress tensor and was the body force, resectively. The stress tensor in an isotroic linear elasticity was written as tr ( ) I () 1 ( D D T (3) when ) where I was the identity tensor, and were the Lame s coefficients which could be written relating to the Young s modulus (E) and Poisson s ratio ( ) in lane stress, E (1 ) E ( 1)(1 ) (4) (5) Then Eq. (1) was arranged in the dislacement vector form that ( D) T { D ( D) t tr( D) I} (6) Equation (6) was comleted to use for the governing equation of the solid body deformation. The initial condition consist of the secified distribution of D and at t = 0. D t The boundary conditions could be of the fixed dislacement, fixed traction, fixed ressure and lane of symmetry. The solution of Eq. (1) could be transformed using Eqs. () and (3) to determine strain and stress, resectively.. Finite Volume Discretization The FVM consist of discretization of the time interval and sace. For the steadystate roblems, the time ste was effectively set to infinite. The sace discretization subdivides the satial domain into the of olyhedral control volume (CV). Every internal face was shared by two CVs. A tyical CV with a comutational oint in its centroid was shown in Fig. 1. A sharing face was also marked and the neighboring CV had oint in its cenroid [4]. Equation (6) was exressed in the other form that t ( D) T {( ) D} { ( D) tr( D) I ( ) D} (7) For this research, considering and as constants thus ( tr( D) I) ( D) (8) and ( D) T ( D) (9) Then Eq. (7) was rearranged that ( D) {( ) D} t (10) The discretization of FVM uses the integral form of Eq. (10) over the CV around oint with the volume V. Using the Gauss theorem the following equation was t V ( D) dv t Z Y X S f ( )( D) nds dv (11) V Fig. 1 Tyical control volume system S f f
3 CST0 191 October, 011, Krabi The volume integral terms were evaluated using the midoint rules [5]. The surface integral terms were slit into the sum of integrals over the cell faces and also evaluated using the midoint rules. The one equation assembled for each CV was a D a D S (1) when and where a a D ( ) S f, (13) d D n ( ) S f (14) n S V F (15) B d was the distance between the midoint and the midoint of neighbor CV (). Thus creating a system of algebric equations A D S (16) where [A] was the sarse matrix with coefficients a on the diagonal and a off the diagonal, [ D ] was the dislacement vector for all CVs..3 Validation of written C++ codes The discretization method described in section. had been imlemented in the Oen Source Software, OenFOAM, using C++ codes. The written code was validated with lane stress roblems as the test cases, which comrised the tension of a thin rectangular late with and without a circular hole at the center and the cross section of the late was uniform. Figure shows the geometry of the rectangular late. According to Timoshenko and Goodier [6], an analytical solution for infinite width late was justified for the finite width late. When the width of the late was not less than four times of the hole diameters, the error for max was less than 6 ercent. Based on Walter and Deborah [7], when the length of the late was longer than d n Fig. Rectangular late subjected to a uniform tension of magnitude T its width, the late could be considered as an infinite length and the effect from both end of the late was negligible. Therefore, the width was 4 m, the length was 8 m and the thickness was m. The hole situated at the center of the rectangular late and had a diameter of 1 m. All materials were assumed isotroic linear elastic. The Young s modulus was assumed Pa and Poisson s ratio was assumed 0.3. Three cases of a rectangular late subjected to a uniform tension of magnitude T in the x direction were analyzed. The first case was a rectangular late without a hole. The second case was a rectangular late with the hole. The third case was a quarter of a late with the hole. This case was modeled based on symmetry of the late. The magnitude of the tension used for all cases were Pa distributed uniformly along both sides of the rectangle as shown in Fig.. The late was discretized as shown in Fig. 3. The results of FVM were comared with the analytical solutions for all cases. 3. Results and Discussion 3.1 Rectangular late without the hole In this load case, the late was subjected to a uniform tension of Pa. Therefore, this was an analytical solution of for this case. The distribution of, and calculated by FVM was shown in Fig. 4.
4 Fig. 3 Discretized lates for FVM Fig. 4 The distribution of, CST0 191 October, 011, Krabi and for the first case calculated by FVM It could be seen that the stress distributions were slightly difference between the center and the edges of the late. The average values of, and from every center oint of cells were , and Pa, resectively. These were due to aroximation error, which was the nature of numerical methods. When comared with the analytical result, it was found that the average difference of from every oint between the FVM and analytical solution was only %. 3. Rectangular late with the hole The analytical solution for this case was based on a late with infinite width [6]. 4 T a T 3a 4a rr 1 1 cos 4 r r r 4 T a T 3a 1 1 cos 4 r r 4 T 3a a r 1 sin 4 r r (17) (18) (19), where a was a radius of a circle, r was a radius from center of the late to any oints on the late and was the angle between the x axis of the late to the radius (Fig. ). The rr, and r could be transformed to xy coordinates by the relation [6] rr cos sin r sin cos (0) sin cos sin cos (1) rr r rr sin cos r cos sin Therefore the analytical solution in terms of, and was [8] 4 a 3 3 a (3) () T 1 cos cos 4 cos 4 4 r r 4 a 1 3 a T cos cos 4 cos 4 4 r r 4 a 1 3 a T sin sin 4 sin 4 4 r r 1 where r x y and = tan y. x (4) (5) Figure 5 shows the distribution of calculated by FVM. As, and
5 exected, the local high stress was observed around the hole. The, and were comared with analytical solution using 16 oints as shown in Fig. 6. These oints were selected based on radial and angle osition around the hole at r = 0.5, 1, 1.5, m and = 0, 30, 60, 90. The order of oints was arranged by starting from r = 0.5 m to r = m and = 0 to 90. The first oint was the osition r = 0.5 m and = 0. CST0 191 October, 011, Krabi The comarison of the, and were shown in Tables. 1, and 3, resectively. It could be seen that the difference of stresses of oints at the boundary of the hole were higher than the other oints. This effect was received from the curvature of the cells at the boundary. Table. 1 Comarison of analytical solution for the second case FVM (Pa) Analytical (Pa) analytical (%) Table. Comarison of analytical solution for the second case Fig. 5 The distribution of, and for the second case calculated by FVM 16 4 Fig. 6 s for stress comarison for the second case 1 13 FVM (Pa) Analytical (Pa) analytical (%)
6 Table. 3 Comarison of analytical solution for the second case FVM (Pa) Analytical (Pa) CST0 191 October, 011, Krabi analytical (%) A quarter of the rectangular late with the hole For this case, the distribution of, and calculated by FVM was shown in Fig. 7. The localized high stress could be observed around the hole. The stress distribution attern for this case was similar to Fig. 5. Figure 8 shows the location of oints for stress comarison calculated by FVM and analytical solution. The analytical solution was calculated by the same equation as for the second case. From Tables. 4, 5 and 6, it could be seen that the difference of stresses were high at the boundary of the hole as same as the second case. 4. Conclusion In this research, the FVM that was develoed for comutational fluid dynamics was alied to analyze solid mechanics. The lane stress roblems were used as a test case. The code for the FVM was written using C++ language and imlemented in the Oen Source Software, OenFOAM. The obtained results were corresonded to the analytical solution. For the Fig. 7 The distribution of, and for the third case calculated by FVM Fig. 8 s for stress comarison for the third case Table. 4 Comarison of analytical solution for the third case 1 FVM (Pa) 13 Analytical (Pa) analytical (%)
7 Table. 5 Comarison of analytical solution for the third case FVM (Pa) Analytical (Pa) CST0 191 October, 011, Krabi analytical (%) Table. 6 Comarison of analytical solution for the third case FVM (Pa) Analytical (Pa) analytical (%) late without the hole, the results were slightly difference from the analytical solution. For the late with a hole, the FVM could cature the localized stress concentration. The quarter of a late with a hole had the results similar to the full late. Therefore, the symmetry model was justified. The overall difference of the stress calculated by FVM from analytical solution was less than 1.68 %. These results suggested the otential of using FVM for stress analysis and it will be used for analyzing solidfluid interaction. 5. References [1] Romero, J.A., Ramírez, O., Fortanell, J.M., Martinez, M. and Lozano, A. (006). Analysis of lateral sloshing forces within road containers with high fill levels, Proceedings of the Institution of Mechanical Engineers, Part D: Journal of Automobile Engineering, vol.0(3), [] Kim, H.S. and Lee, Y.S. (008). Otimization design technique for reduction of sloshing by evolutionary methods, Journal of Mechanical Science and Technology, vol.(1), [3] Thusyanthan, N.I. and Madabhushi, S.P.G. (008). Tsunami wave loading on coastal houses: A model aroach, Proceedings of the Institution of Civil Engineers: Civil Engineering, vol.161(), [4] OenFOAM. (009). User Guide Version 1.6, edition, OenCFD Limited, London. [5] Jasak, H. and Weller, H.G. (000). Alication of the finite volume method and unstructured meshes to linear elasticity, International Journal for Numerical Methods in Engineering, vol.48(), [6] Timoshenko, S.P. and Goodier, J.N. (198). Theory of Elasticity, 3 rd edition, McGrawHill, London. [7] Walter, D.P. and Deborah, F.P. (008). Peterson's Stress Concentration Factors, 3 rd edition, John Wiley & Sons, Inc., New Jersey. [8] Demirdzic, I. and Muzaferija, S. (1994). Finite volume method for stress analysis in comlex domains, International Journal for Numerical Methods in Engineering, vol.37(1),
Lower bound solutions for bearing capacity of jointed rock
Comuters and Geotechnics 31 (2004) 23 36 www.elsevier.com/locate/comgeo Lower bound solutions for bearing caacity of jointed rock D.J. Sutcliffe a, H.S. Yu b, *, S.W. Sloan c a Deartment of Civil, Surveying
More informationChurilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics
Churilova Maria Saint-Petersburg State Polytechnical University Deartment of Alied Mathematics Technology of EHIS (staming) alied to roduction of automotive arts The roblem described in this reort originated
More informationA Simple And Efficient FEM-Implementation Of The Modified Mohr-Coulomb Criterion Clausen, Johan Christian; Damkilde, Lars
Aalborg Universitet A Simle And Efficient FEM-Imlementation Of The Modified Mohr-Coulomb Criterion Clausen, Johan Christian; Damkilde, Lars Published in: Proceedings of the 9th Nordic Seminar on Comutational
More informationFE FORMULATIONS FOR PLASTICITY
G These slides are designed based on the book: Finite Elements in Plasticity Theory and Practice, D.R.J. Owen and E. Hinton, 1970, Pineridge Press Ltd., Swansea, UK. 1 Course Content: A INTRODUCTION AND
More informationA SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE
THE 19 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE K.W. Gan*, M.R. Wisnom, S.R. Hallett, G. Allegri Advanced Comosites
More informationNonlinear Static Analysis of Cable Net Structures by Using Newton-Raphson Method
Nonlinear Static Analysis of Cable Net Structures by Using Newton-Rahson Method Sayed Mahdi Hazheer Deartment of Civil Engineering University Selangor (UNISEL) Selangor, Malaysia hazheer.ma@gmail.com Abstract
More informationREFINED STRAIN ENERGY OF THE SHELL
REFINED STRAIN ENERGY OF THE SHELL Ryszard A. Walentyński Deartment of Building Structures Theory, Silesian University of Technology, Gliwice, PL44-11, Poland ABSTRACT The aer rovides information on evaluation
More informationZEEMAN EFFECT: p...(1). Eigenfunction for this Hamiltonian is specified by
ZEEMAN EFFECT: Zeeman Effect is a magneto-otical henomenon discovered by Zeeman in 1896. He observed that when an atom (light soce) is laced in an external magnetic field, the sectral lines it emits are
More informationCombining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)
Combining Logistic Regression with Kriging for Maing the Risk of Occurrence of Unexloded Ordnance (UXO) H. Saito (), P. Goovaerts (), S. A. McKenna (2) Environmental and Water Resources Engineering, Deartment
More information10 th Jubilee National Congress on Theoretical and Applied Mechanics, Varna September 2005
th Jubilee National Congress on Theoretical and Alied Mechanics, Varna 6 Setember 25 FEM APPLIED IN HYDRO-MECHANICAL COUPLED ANALYSIS OF A SLOPE A. Yanakieva, M. Datcheva, R. Iankov, F. Collin 2, A. Baltov
More informationImplementation of an Isotropic Elastic-Viscoplastic Model for Soft Soils using COMSOL Multiphysics
Imlementation of an Isotroic Elastic-Viscolastic Model for Soft Soils using COMSOL Multihysics M. Olsson 1,, T. Wood 1,, C. Alén 1 1 Division of GeoEngineering, Chalmers University of Technology, Gothenburg,
More informationSHAPE OPTOMIZATION OF H-BEAM FLANGE FOR MAXIMUM PLASTIC ENERGY DISSIPATION
The Fourth China-Jaan-Korea Joint Symosium on Otimization of Structural and Mechanical Systems Kunming, Nov. 6-9, 2006, China SHAPE OPTOMIZATION OF H-BEAM FLANGE FOR MAXIMUM PLASTIC ENERGY DISSIPATION
More informationDeformation Effect Simulation and Optimization for Double Front Axle Steering Mechanism
0 4th International Conference on Comuter Modeling and Simulation (ICCMS 0) IPCSIT vol. (0) (0) IACSIT Press, Singaore Deformation Effect Simulation and Otimization for Double Front Axle Steering Mechanism
More informationAnalysis of cold rolling a more accurate method
Analysis of cold rolling a more accurate method 1.1 Rolling of stri more accurate slab analysis The revious lecture considered an aroximate analysis of the stri rolling. However, the deformation zone in
More informationFlexible Pipes in Trenches with Stiff Clay Walls
Flexible Pies in Trenches with Stiff Clay Walls D. A. Cameron University of South Australia, South Australia, Australia J. P. Carter University of Sydney, New South Wales, Australia Keywords: flexible
More informationFEM simulation of a crack propagation in a round bar under combined tension and torsion fatigue loading
FEM simulation of a crack roagation in a round bar under combined tension and torsion fatigue loading R.Citarella, M.Leore Det. of Industrial Engineering University of Salerno - Fisciano (SA), Italy. rcitarella@unisa.it
More informationKeywords: pile, liquefaction, lateral spreading, analysis ABSTRACT
Key arameters in seudo-static analysis of iles in liquefying sand Misko Cubrinovski Deartment of Civil Engineering, University of Canterbury, Christchurch 814, New Zealand Keywords: ile, liquefaction,
More informationDomain Dynamics in a Ferroelastic Epilayer on a Paraelastic Substrate
Y. F. Gao Z. Suo Mechanical and Aerosace Engineering Deartment and Princeton Materials Institute, Princeton University, Princeton, NJ 08544 Domain Dynamics in a Ferroelastic Eilayer on a Paraelastic Substrate
More information8.7 Associated and Non-associated Flow Rules
8.7 Associated and Non-associated Flow Rules Recall the Levy-Mises flow rule, Eqn. 8.4., d ds (8.7.) The lastic multilier can be determined from the hardening rule. Given the hardening rule one can more
More informationA General Damage Initiation and Evolution Model (DIEM) in LS-DYNA
9th Euroean LS-YNA Conference 23 A General amage Initiation and Evolution Model (IEM) in LS-YNA Thomas Borrvall, Thomas Johansson and Mikael Schill, YNAmore Nordic AB Johan Jergéus, Volvo Car Cororation
More informationPrediction of the Excitation Force Based on the Dynamic Analysis for Flexible Model of a Powertrain
Prediction of the Excitation Force Based on the Dynamic Analysis for Flexible Model of a Powertrain Y.S. Kim, S.J. Kim, M.G. Song and S.K. ee Inha University, Mechanical Engineering, 53 Yonghyun Dong,
More informationFinite Element Analysis of V-Bending of Polypropylene Using Hydrostatic-Pressure-Dependent Plastic Constitutive Equation*
Materials Transactions, Vol. 48, No. 1 (7). 6 to 664 #7 The Jaan Society for Technology of Plasticity Finite Element Analysis of V-Bending of Polyroylene Using Hydrostatic-Pressure-Deendent Plastic onstitutive
More informationA PIEZOELECTRIC BERNOULLI-EULER BEAM THEORY CONSIDERING MODERATELY CONDUCTIVE AND INDUCTIVE ELECTRODES
Proceedings of the 6th International Conference on Mechanics and Materials in Design, Editors: J.F. Silva Gomes & S.A. Meguid, P.Delgada/Aores, 26-30 July 2015 PAPER REF: 5513 A PIEZOELECTRIC BERNOULLI-EULER
More informationCMSC 425: Lecture 4 Geometry and Geometric Programming
CMSC 425: Lecture 4 Geometry and Geometric Programming Geometry for Game Programming and Grahics: For the next few lectures, we will discuss some of the basic elements of geometry. There are many areas
More informationIsogeometric analysis based on scaled boundary finite element method
IOP Conference Series: Materials Science and Engineering Isogeometric analysis based on scaled boundary finite element method To cite this article: Y Zhang et al IOP Conf. Ser.: Mater. Sci. Eng. 37 View
More informationpp physics, RWTH, WS 2003/04, T.Hebbeker
1. PP TH 03/04 Accelerators and Detectors 1 hysics, RWTH, WS 2003/04, T.Hebbeker 2003-12-03 1. Accelerators and Detectors In the following, we concentrate on the three machines SPS, Tevatron and LHC with
More informationVIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES
Journal of Sound and Vibration (998) 22(5), 78 85 VIBRATION ANALYSIS OF BEAMS WITH MULTIPLE CONSTRAINED LAYER DAMPING PATCHES Acoustics and Dynamics Laboratory, Deartment of Mechanical Engineering, The
More informationUniversity of North Carolina-Charlotte Department of Electrical and Computer Engineering ECGR 4143/5195 Electrical Machinery Fall 2009
University of North Carolina-Charlotte Deartment of Electrical and Comuter Engineering ECG 4143/5195 Electrical Machinery Fall 9 Problem Set 5 Part Due: Friday October 3 Problem 3: Modeling the exerimental
More informationModelling a Partly Filled Road Tanker during an Emergency Braking
Proceedings of the World Congress on Engineering and Comuter Science 217 Vol II, October 25-27, 217, San Francisco, USA Modelling a Partly Filled Road Tanker during an Emergency Braking Frank Otremba,
More informationFINITE ELEMENT SIMULATION OF BACKWARD EXTRUSION DIES USING SHRINK FITTED RINGS
Vienna, AUSTRIA, 000, Setember 5-8 FINITE ELEMENT SIMULATION OF BACKWARD EXTRUSION DIES USING SHRINK FITTED RINGS M.E. Yurci 1, V. Taskin and N. Urkmez 1 Yildiz Technical University, Istanbul, Turkey Trakya
More informationMANDATORY APPENDIX 41 ELECTRIC IMMERSION HEATER ELEMENT SUPPORT PLATES
41-1 41-5 Page 1 of 5 No changes, age is included for reference. MANDATORY APPENDIX 41 ELECTRIC IMMERSION HEATER ELEMENT SUPPORT PLATES 41-1 SCOPE 41-3 41-1.1 The rules in this Mandatory Aendix cover the
More informationUnderstanding DPMFoam/MPPICFoam
Understanding DPMFoam/MPPICFoam Jeroen Hofman March 18, 2015 In this document I intend to clarify the flow solver and at a later stage, the article-fluid and article-article interaction forces as imlemented
More informationOptimization of Gear Design and Manufacture. Vilmos SIMON *
7 International Conference on Mechanical and Mechatronics Engineering (ICMME 7) ISBN: 978--6595-44- timization of Gear Design and Manufacture Vilmos SIMN * Budaest Universit of Technolog and Economics,
More informationPrinciples of Computed Tomography (CT)
Page 298 Princiles of Comuted Tomograhy (CT) The theoretical foundation of CT dates back to Johann Radon, a mathematician from Vienna who derived a method in 1907 for rojecting a 2-D object along arallel
More informationHydrogen-assisted stress corrosion cracking simulation using the stress-modified fracture strain model
Journal of Mechanical Science and Technology 26 (8) (2012) 2631~2638 www.sringerlink.com/content/1738-494x DOI 10.1007/s12206-012-0642-x Hydrogen-assisted stress corrosion cracking simulation using the
More informationPreliminary Uncertainty Estimation of the Pressure Distortion Coefficient of a Pressure. Balance by FEM Calculations
Preliminary Uncertainty Estimation of the Pressure Distortion Coefficient of a Pressure Balance by FEM Calculations G. Molinar*, M. Bergoglio*, G. Mosso*,G. Buonanno**, M. Dell Isola** * Istituto di Metrologia
More informationEstimating Laterally Loaded Pile Response
Estimating Laterally Loaded Pile Resonse C. Y. Lee, PhD Deartment of Civil Engineering College of Engineering Universiti Tenaga Nasional ABSTRACT A simlified and ractical aroach for estimating the laterally
More informationANALYTIC APPROXIMATE SOLUTIONS FOR FLUID-FLOW IN THE PRESENCE OF HEAT AND MASS TRANSFER
Kilicman, A., et al.: Analytic Aroximate Solutions for Fluid-Flow in the Presence... THERMAL SCIENCE: Year 8, Vol., Sul.,. S59-S64 S59 ANALYTIC APPROXIMATE SOLUTIONS FOR FLUID-FLOW IN THE PRESENCE OF HEAT
More informationAn-Najah National University Civil Engineering Departemnt. Fluid Mechanics. Chapter [2] Fluid Statics
An-Najah National University Civil Engineering Deartemnt Fluid Mechanics Chater [2] Fluid Statics 1 Fluid Statics Problems Fluid statics refers to the study of fluids at rest or moving in such a manner
More informationNumerical and experimental investigation on shot-peening induced deformation. Application to sheet metal forming.
Coyright JCPDS-International Centre for Diffraction Data 29 ISSN 197-2 511 Numerical and exerimental investigation on shot-eening induced deformation. Alication to sheet metal forming. Florent Cochennec
More informationA Special Case Solution to the Perspective 3-Point Problem William J. Wolfe California State University Channel Islands
A Secial Case Solution to the Persective -Point Problem William J. Wolfe California State University Channel Islands william.wolfe@csuci.edu Abstract In this aer we address a secial case of the ersective
More informationNumerical Methods for Particle Tracing in Vector Fields
On-Line Visualization Notes Numerical Methods for Particle Tracing in Vector Fields Kenneth I. Joy Visualization and Grahics Research Laboratory Deartment of Comuter Science University of California, Davis
More informationThe effect of dynamic bending moments on the ratchetting behavior of stainless steel pressurized piping elbows
International Journal of echanical Engineering and Alications 2014; 2(2): 31-37 Published online ay 30, 2014 (htt://www.scienceublishinggrou.com/j/ijmea) doi: 10.11648/j.ijmea.20140202.12 The effect of
More informationChapter 2 Introduction to the Stiffness (Displacement) Method. The Stiffness (Displacement) Method
CIVL 7/87 Chater - The Stiffness Method / Chater Introdction to the Stiffness (Dislacement) Method Learning Objectives To define the stiffness matrix To derive the stiffness matrix for a sring element
More informationOn split sample and randomized confidence intervals for binomial proportions
On slit samle and randomized confidence intervals for binomial roortions Måns Thulin Deartment of Mathematics, Usala University arxiv:1402.6536v1 [stat.me] 26 Feb 2014 Abstract Slit samle methods have
More informationMultivariable Generalized Predictive Scheme for Gas Turbine Control in Combined Cycle Power Plant
Multivariable Generalized Predictive Scheme for Gas urbine Control in Combined Cycle Power Plant L.X.Niu and X.J.Liu Deartment of Automation North China Electric Power University Beiing, China, 006 e-mail
More informationAutomatic Generation and Integration of Equations of Motion for Linked Mechanical Systems
Automatic Generation and Integration of Equations of Motion for Linked Mechanical Systems D. Todd Griffith a, John L. Junkins a, and James D. Turner b a Deartment of Aerosace Engineering, Texas A&M University,
More informationOn the elasticity of transverse isotropic soft tissues (L)
J_ID: JAS DOI: 10.1121/1.3559681 Date: 17-March-11 Stage: Page: 1 Total Pages: 5 ID: 3b2server Time: 12:24 I Path: //xinchnasjn/aip/3b2/jas#/vol00000/110099/appfile/ai-jas#110099 1 2 3 4 5 6 7 AQ18 9 10
More informationMODEL-BASED MULTIPLE FAULT DETECTION AND ISOLATION FOR NONLINEAR SYSTEMS
MODEL-BASED MULIPLE FAUL DEECION AND ISOLAION FOR NONLINEAR SYSEMS Ivan Castillo, and homas F. Edgar he University of exas at Austin Austin, X 78712 David Hill Chemstations Houston, X 77009 Abstract A
More informationNumerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model
Numerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model Q. Sun, Y.J. Liu, R.H, Jin School of Aeronautics, Northwestern Polytechnical University, Xi an 710072,
More informationFinite Element Solutions for Geotechnical Engineering
Release Notes Release Date: June, 2017 Product Ver.: GTSNX 2017(v1.1) Integrated Solver Otimized for the next generation 64-bit latform Finite Element Solutions for Geotechnical Engineering 1. Analysis
More informationINTRODUCING THE SHEAR-CAP MATERIAL CRITERION TO AN ICE RUBBLE LOAD MODEL
Symosium on Ice (26) INTRODUCING THE SHEAR-CAP MATERIAL CRITERION TO AN ICE RUBBLE LOAD MODEL Mohamed O. ElSeify and Thomas G. Brown University of Calgary, Calgary, Canada ABSTRACT Current ice rubble load
More informationSolved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.
Solved Problems Solved Problems P Solve the three simle classification roblems shown in Figure P by drawing a decision boundary Find weight and bias values that result in single-neuron ercetrons with the
More informationFinite Element Evaluation of J-integral in 3D for Nuclear Grade Graphite Using COMSOL-Multiphysics
Finite Element Evaluation of J-integral in 3D for Nuclear Grade Grahite Using COMSOL-Multihysics wani Bhushan *1, S.K. Panda 1, Sanjeev Kumar Singh 1, Debashis Khan 1 1 Indian Institute of Technology (BHU)
More informationSpin Diffusion and Relaxation in a Nonuniform Magnetic Field.
Sin Diffusion and Relaxation in a Nonuniform Magnetic Field. G.P. Berman, B. M. Chernobrod, V.N. Gorshkov, Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545 V.I. Tsifrinovich
More informationClosed-form solution for the size of plastic zone in an edge-cracked strip
University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Faculty Publications from the Deartment of Engineering Mechanics Mechanical & Materials Engineering, Deartment of May 2002
More informationMODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL
Technical Sciences and Alied Mathematics MODELING THE RELIABILITY OF CISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL Cezar VASILESCU Regional Deartment of Defense Resources Management
More informationAdam Paweł Zaborski. 8 Plasticity. reloading. 1. Bauschinger s effect. 2. unchanged yielding limit. 3. isotropic hardening
8 lasticity Introduction Definitions loading/unloading/reloading words commonly used in lasticity lastic strain a ermanent strain that doesn t vanish after unloading, lastically assive rocess a rocess
More informationCharacteristics of Beam-Based Flexure Modules
Shorya Awtar e-mail: shorya@mit.edu Alexander H. Slocum e-mail: slocum@mit.edu Precision Engineering Research Grou, Massachusetts Institute of Technology, Cambridge, MA 039 Edi Sevincer Omega Advanced
More informationarxiv: v1 [math-ph] 10 Sep 2010
Calculation of generalized ubbell rectangular source integral Jonathan Murley and Nasser Saad arxiv:9.7v [math-h] Se Deartment of Mathematics and Statistics, University of Prince Edward Island Charlottetown,
More informationAn Investigation on the Numerical Ill-conditioning of Hybrid State Estimators
An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators S. K. Mallik, Student Member, IEEE, S. Chakrabarti, Senior Member, IEEE, S. N. Singh, Senior Member, IEEE Deartment of Electrical
More informationMaximum Entropy and the Stress Distribution in Soft Disk Packings Above Jamming
Maximum Entroy and the Stress Distribution in Soft Disk Packings Above Jamming Yegang Wu and S. Teitel Deartment of Physics and Astronomy, University of ochester, ochester, New York 467, USA (Dated: August
More informationMeasurement of cyclone separator
Measurement of cyclone searator. Aim of the measurement Cyclones are widely used in industry (in food and chemical industry, in energy technology and in buildings) to remove dust and other articles from
More informationSession 5: Review of Classical Astrodynamics
Session 5: Review of Classical Astrodynamics In revious lectures we described in detail the rocess to find the otimal secific imulse for a articular situation. Among the mission requirements that serve
More informationUncorrelated Multilinear Principal Component Analysis for Unsupervised Multilinear Subspace Learning
TNN-2009-P-1186.R2 1 Uncorrelated Multilinear Princial Comonent Analysis for Unsuervised Multilinear Subsace Learning Haiing Lu, K. N. Plataniotis and A. N. Venetsanooulos The Edward S. Rogers Sr. Deartment
More informationCREEP AND SHRINKAGE EFFECT ON THE DYNAMIC ANALYSIS OF REINFORCED CONCRETE SLAB-AND-BEAM STRUCTURES
ECCM 99 Euroean Conference on Comutational Mechanics August 31 Setember 3 München, Germany CREEP AND SHRINKAGE EFFECT ON THE DYNAMIC ANALYSIS OF REINFORCED CONCRETE SLABANDBEAM STRUCTURES Evangelos J.
More informationHomogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers
Excert from the roceedings of the COMSOL Conference 1 aris Homogeneous and Inhomogeneous Model for Flow and Heat ransfer in orous Materials as High emerature Solar Air Receivers Olena Smirnova 1 *, homas
More information+++ Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic +++ Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic
Modeling of Structural-dynamic Systems by UML Statecharts in AnyLogic Daniel Leitner, Johannes Krof, Günther Zauner, TU Vienna, Austria, dleitner@osiris.tuwien.ac.at Yuri Karov, Yuri Senichenkov, Yuri
More informationMulti-Operation Multi-Machine Scheduling
Multi-Oeration Multi-Machine Scheduling Weizhen Mao he College of William and Mary, Williamsburg VA 3185, USA Abstract. In the multi-oeration scheduling that arises in industrial engineering, each job
More informationHigher order theory for analytic saddle point approximations to the Ρ Ρ and Ρ Ś reflected arrivals at a solid/solid interface
Higher order theory for analytic saddle oint aroximations to the Ρ Ρ and Ρ Ś reflected arrivals at a solid/solid interface P.F Daley ABSTACT The high frequency solution to the roblem of a Ρ P and Ρ S reflected
More informationWolfgang POESSNECKER and Ulrich GROSS*
Proceedings of the Asian Thermohysical Proerties onference -4 August, 007, Fukuoka, Jaan Paer No. 0 A QUASI-STEADY YLINDER METHOD FOR THE SIMULTANEOUS DETERMINATION OF HEAT APAITY, THERMAL ONDUTIVITY AND
More informationAn equivalent diameter in calculating borehole thermal resistance for spiral coil type GHE
The 01 World Congress on Advances in Civil, Environmental, and Materials esearch (ACEM 1) Seoul, Korea, August 6-30, 01 An equivalent diameter in calculating borehole thermal resistance for siral coil
More informationUniformly best wavenumber approximations by spatial central difference operators: An initial investigation
Uniformly best wavenumber aroximations by satial central difference oerators: An initial investigation Vitor Linders and Jan Nordström Abstract A characterisation theorem for best uniform wavenumber aroximations
More informationAnalysis of High-Altitude Ionization Gauge Measurements Using the Direct Simulation Monte Carlo Method
th AIAA Thermohysics Conference June - July, Portland, Oregon AIAA - Analysis of High-Altitude Ionization Gauge Measurements Using the Direct Simulation Monte Carlo Method Quanhua Sun *, Chunei Cai, and
More informationState Estimation with ARMarkov Models
Deartment of Mechanical and Aerosace Engineering Technical Reort No. 3046, October 1998. Princeton University, Princeton, NJ. State Estimation with ARMarkov Models Ryoung K. Lim 1 Columbia University,
More informationRANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES
RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES AARON ZWIEBACH Abstract. In this aer we will analyze research that has been recently done in the field of discrete
More informationPreconditioning techniques for Newton s method for the incompressible Navier Stokes equations
Preconditioning techniques for Newton s method for the incomressible Navier Stokes equations H. C. ELMAN 1, D. LOGHIN 2 and A. J. WATHEN 3 1 Deartment of Comuter Science, University of Maryland, College
More informationTime Domain Calculation of Vortex Induced Vibration of Long-Span Bridges by Using a Reduced-order Modeling Technique
2017 2nd International Conference on Industrial Aerodynamics (ICIA 2017) ISBN: 978-1-60595-481-3 Time Domain Calculation of Vortex Induced Vibration of Long-San Bridges by Using a Reduced-order Modeling
More informationCHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules
CHAPTER-II Control Charts for Fraction Nonconforming using m-of-m Runs Rules. Introduction: The is widely used in industry to monitor the number of fraction nonconforming units. A nonconforming unit is
More informationJournal of Computational and Nonlinear Dynamics JANUARY 2010, Vol. 5 / Copyright 2010 by ASME
A. L. Schwab Laboratory for Engineering Mechanics, Delft University of Technology, Mekelweg 2, NL-2628 CD Delft, The Netherlands e-mail: a.l.schwab@tudelft.nl J. P. Meijaard Laboratory of Mechanical Automation
More informationPhase Equilibrium Calculations by Equation of State v2
Bulletin of Research Center for Comuting and Multimedia Studies, Hosei University, 7 (3) Published online (htt://hdl.handle.net/4/89) Phase Equilibrium Calculations by Equation of State v Yosuke KATAOKA
More informationSAMARIS. Draft report
SAMARIS Work Package 5 - Performance based secifications Selection and evaluation of models for rediction of ermanent deformations of unbound granular materials in road avements Draft reort Pierre Hornych
More informationChapter 1 Fundamentals
Chater Fundamentals. Overview of Thermodynamics Industrial Revolution brought in large scale automation of many tedious tasks which were earlier being erformed through manual or animal labour. Inventors
More informationOptimal Design of Truss Structures Using a Neutrosophic Number Optimization Model under an Indeterminate Environment
Neutrosohic Sets and Systems Vol 14 016 93 University of New Mexico Otimal Design of Truss Structures Using a Neutrosohic Number Otimization Model under an Indeterminate Environment Wenzhong Jiang & Jun
More informationON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES CONTROL PROBLEMS
Electronic Transactions on Numerical Analysis. Volume 44,. 53 72, 25. Coyright c 25,. ISSN 68 963. ETNA ON THE DEVELOPMENT OF PARAMETER-ROBUST PRECONDITIONERS AND COMMUTATOR ARGUMENTS FOR SOLVING STOKES
More informationComputational models of diamond anvil cell compression
UDC 519.6 Computational models of diamond anvil cell compression A. I. Kondrat yev Independent Researcher, 5944 St. Alban Road, Pensacola, Florida 32503, USA Abstract. Diamond anvil cells (DAC) are extensively
More informationCalculation of gravity due to a vertical cylinder using a spherical harmonic series and numerical integration
CSIRO PUBISHING Exloration Geohysics htt://dx.doi.org/.7/eg43 Calculation of gravity due to a vertical cylinder using a sherical harmonic series and numerical integration Sung-Ho Na,3 Hyoungrea Rim,3,4
More informationFor q 0; 1; : : : ; `? 1, we have m 0; 1; : : : ; q? 1. The set fh j(x) : j 0; 1; ; : : : ; `? 1g forms a basis for the tness functions dened on the i
Comuting with Haar Functions Sami Khuri Deartment of Mathematics and Comuter Science San Jose State University One Washington Square San Jose, CA 9519-0103, USA khuri@juiter.sjsu.edu Fax: (40)94-500 Keywords:
More informationCovariance Matrix Estimation for Reinforcement Learning
Covariance Matrix Estimation for Reinforcement Learning Tomer Lancewicki Deartment of Electrical Engineering and Comuter Science University of Tennessee Knoxville, TN 37996 tlancewi@utk.edu Itamar Arel
More informationANALYSIS OF ULTRA LOW CYCLE FATIGUE PROBLEMS WITH THE BARCELONA PLASTIC DAMAGE MODEL
XII International Conerence on Comutational Plasticity. Fundamentals and Alications COMPLAS XII E. Oñate, D.R.J. Owen, D. Peric and B. Suárez (Eds) ANALYSIS OF ULTRA LOW CYCLE FATIGUE PROBLEMS WITH THE
More informationDimensional perturbation theory for Regge poles
Dimensional erturbation theory for Regge oles Timothy C. Germann Deartment of Chemistry, University of California, Berkeley, California 94720 Sabre Kais Deartment of Chemistry, Purdue University, West
More informationStudy on Characteristics of Sound Absorption of Underwater Visco-elastic Coated Compound Structures
Vol. 3, No. Modern Alied Science Study on Characteristics of Sound Absortion of Underwater Visco-elastic Coated Comound Structures Zhihong Liu & Meiing Sheng College of Marine Northwestern Polytechnical
More informationPressure-sensitivity Effects on Toughness Measurements of Compact Tension Specimens for Strain-hardening Solids
American Journal of Alied Sciences (9): 19-195, 5 ISSN 1546-939 5 Science Publications Pressure-sensitivity Effects on Toughness Measurements of Comact Tension Secimens for Strain-hardening Solids Abdulhamid
More informationEstimation of the large covariance matrix with two-step monotone missing data
Estimation of the large covariance matrix with two-ste monotone missing data Masashi Hyodo, Nobumichi Shutoh 2, Takashi Seo, and Tatjana Pavlenko 3 Deartment of Mathematical Information Science, Tokyo
More informationCFD Modelling of Mass Transfer and Interfacial Phenomena on Single Droplets
Euroean Symosium on Comuter Arded Aided Process Engineering 15 L. Puigjaner and A. Esuña (Editors) 2005 Elsevier Science B.V. All rights reserved. CFD Modelling of Mass Transfer and Interfacial Phenomena
More informationThe directivity of the forced radiation of sound from panels and openings including the shadow zone
The directivity of the forced radiation of sound from anels and oenings including the shadow zone J. Davy RMIT University, Alied Physics, GPO Box 476V, 3001 Melbourne, Victoria, Australia john.davy@rmit.edu.au
More informationAnalysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia
roceedings World Geothermal Congress 00 Bali, Indonesia, 5-9 Aril 00 Analysis of ressure Transient Resonse for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia Jorge A.
More information1 Properties of Spherical Harmonics
Proerties of Sherical Harmonics. Reetition In the lecture the sherical harmonics Y m were introduced as the eigenfunctions of angular momentum oerators lˆz and lˆ2 in sherical coordinates. We found that
More informationKEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS
4 th International Conference on Earthquake Geotechnical Engineering June 2-28, 27 KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS Misko CUBRINOVSKI 1, Hayden BOWEN 1 ABSTRACT Two methods for analysis
More informationPhysical based Schottky barrier diode modeling for THz applications
Downloaded from orbit.dtu.dk on: Jan 6, 18 Physical based Schottky barrier diode modeling THz alications Yan, Lei; Krozer, iktor; Michaelsen, Rasmus Schandorh; Durhuus, Torsten; Johansen, Tom Keinicke
More information