FRP SPH(Smoothed Partcle Hydrodynamcs) (3) MPS(Movng Partcle sem-mplct) (4) MPS (5) (6) (7, 8) (9) (10) (11) Tama and Koshzuka LSMPS(Least Squares Mov
|
|
- ありおき ありたけ
- 7 years ago
- Views:
Transcription
1 Transactons of JSCES, Paper No Ansotropc Hgh Vscosty Flud Analyss Usng a Partcle Method for Evaluatng CFRTP Press Moldng Process Ryosaku SHINO, Tasuku TAMAI, Sech KOSHIZUKA, Akra MAKI and Takesh ISHIKAWA 1 ( ) 2 ( ) In ths paper, a new ansotropc hgh vscosty model usng a partcle method s proposed to analyze the carbon fber renforced thermoplastcs (CFRTP) press moldng process. Ansotropy of vscosty s consdered wth fber orentaton vector, and LSMPS method s used to mprove relablty and usablty. Analyss of a smple cube pressng process demonstrates the capablty of the present method to analyze ansotropc hgh vscosty flows wthout vtatng the ncompressblty. Key Words: Partcle Method, LSMPS Method, CFRTP, Ansotropc Vscous Flud, Hgh Vscosty 1. (Fber Renforced Plastcs : FRP) FRP FRP FRP (Carbon Fber Renforced Plastcs : CFRP) (Carbon Fber Renforced Thermoplastcs : CFRTP) CFRP CFRTP CFRTP FRP , , c Manuscrpt receved, October 04, 2015; fnal revson, Aprl 05, 2016; publshed, June 30, Copyrght c 2016 by the Japan Socety for Computatonal Engneerng and Scence. FRP FRP FRP FRP FRP (1) FRP CFRP (2) FRP
2 FRP SPH(Smoothed Partcle Hydrodynamcs) (3) MPS(Movng Partcle sem-mplct) (4) MPS (5) (6) (7, 8) (9) (10) (11) Tama and Koshzuka LSMPS(Least Squares Movng Partcle Sem-Implct) (12) LSMPS MPS SPH (12) 2. LSMPS (2) LSMPS 2.1 FRP Lee (13) Verweyst and Tucker (1) Lee (Hele-Shaw ) Hele-Shaw 2 Hele-Shaw (13) Verweyst and Tucker (14) 2 4 Verweyst and Tucker (2) Young consoldaton 2 (15) 2 FRP Du Dt = C 2 u 1 P + f (1) ρ u t P ρ f C Fg.1 r q 動粘度 p 繊維 Y 動粘度 Fber orentaton vector and knematc vscosty. Fg.1 p = {p x, p y, p z }, q = {q x, q y, q z }, r = {r x, r y, r z } Z X
3 Fg.1 X, Y, Z ) Fg.1 p, q, r ϕ ϕ p x q x r x ϕ = p y q y r y p z q z r z FRP (2) FRP ν hor ν ver ν ver < ν hor C 1 ν hor 0 0 (2) C = ϕ T 0 ν ver 0 ϕ (3) 0 0 ν ver 2.2 MPS (4) 10 m 2 /s x = x k + tu k (4) 3ũ 4u k + u k 1 = C 2 ũ 1 2 t ρ Pk + f (5) 2 ϕ k+1 = 3ρ 2 t ũ (6) 3u k+1 3ũ 2 t = 1 ρ ϕk+1 (7) P k+1 = P k + ϕ k+1 (8) x k+1 = x k + t ( u k+1 + u k) 2 (9) x x t u ũ ϕ k+1 ϕ k+1 = P k+1 P k (10) 2.3 LSMPS (5) 2 (6) (7) Tama and Koshzuka Standard LSMPS scheme Type-A (12) 1 1 ( ) T ϕ ϕ ϕ [( ) = H [1] r x, y, z s M [1] 1 ] b [1] r 1 s 0 0 H [1] = 0 r 1 s r 1 s [ M [1] = w ( ) ( ) ( )] x j, r e p [1] x j p [1] x j b [1] = j Λ j Λ [ w ( x j, r e ) p [1] ( x j (11) (12) (13) ) (ϕ j ϕ ) ] (14) x j = x j x (15) p [1] (x) = ( x, y, z ) T (16) Λ = { j 0 x j < re } (17) ϕ H rs (12) Scalng matrx M (13) Moment matrx b (14) p (x) (16) r e dlaton parameter scalng parameter w Tama and Koshzuka (12) ( ) 2 x 1 (0 x < r e ) w (x, r e ) = r e, 0, (r e x ) (18) (5) 1 (6) Standard LSMPS scheme Type- A (12) 2 2 [ H ]9 3 = ( 2 ϕ x 2, M = b = 2 ϕ y 2, ) 2 T ϕ [( ) = H z 2 M 1 ] b 2r 2 s r 2 s r 2 s j Λ j Λ [ w ( x j, r e ) p [ w ( x j, r e ) p ( x j ( x j ) ( )] p x j (19) (20) (21) ) (ϕ j ϕ ) ] (22) p (x) = ( x, y, z, x 2, y 2, z 2, xy, yz, zx ) T (23)
4 2.4 Posson (6) LSMPS ϕ k+1 = 0 (24) Drchlet ϕ k+1 n = 0 (25) Neumann (25) Neumann Tama and koshzuka Constrant LSMPS scheme type-a (12) ϕ ϕ k+1 n = g N (x) (26) (26) 2 ( ) 2 ϕ 2 ϕ 2 T ϕ [( ) = H x 2, y 2, z 2 M + N 1 ( )] b + c (27) (10) ( 3 ) (32) ( ũ x ũ y,ũ z ) 3 (31) (32) (32) 仮想粒子 l 0 n Support 半径 M H rs b (21) (22) N c Fg.2 自由表面粒子 Vrtual partcle of the surface partcle. N = w (0, r e ) p N (x ) p N (x ) (28) c = w (0, r e ) p N (x ) g N (x) (29) p N (x ) = ( n x, n y, n z, 0, 0, 0, 0, 0, 0 ) (30) n = {n x, n y, n z } 2.5 (5) ũ = u wall (31) [ p k I + ρc ũ ] n = 0 (32) p k k n ( (24) (8) p k ) MPS (16) ũ x ũ y ũ z ( ũ x ũ y ũ z ) ũ x ũ y ũ z Fg.2 (5) (32) l 0 1 M (32) A S d ũ d + A S ũ = r S (33) ũ A S (32) ũ d A S d (32) r S ũ d ũ d = A 1 S d r S A 1 S d A S ũ (34) (5) A Id ũ d + A I ũ = r I (35) A I (5) A Id (5) r I (34) (35) ũ d ( ) AI A Id A 1 S d A S ũ = r I A Id A 1 S d r S (36)
5 2.6 LSMPS (24) (32) MPS ( ) (4) FRP (17) (18) (19) (20) 2.7 (32) MPS l 0 x +l 0, l 0, y +l 0, l 0, z +l 0, l 0 Nomura (21) n j Λ x j w ( ) x j, r e n = j Λ x j w ( ) (37) x j, r e Nomura 3. y θ[degree] p r x z p:fber drecton q z y x Fg.3 Table1 Parameter ν ver [m 2 /s] l z V[m/s] Press T[second] h[m] l y l x Problem geometry Parameter of analyss Value ν hor [m 2 /s] ν ver 10, 50, 100, 500, 1000 θ [degree] 0 l x [m] l y [m] l z [m] h [m] T [s] ρ [kg/m 3 ] V [m/s] l 0 [m] r e [m] 3.1 l 0 [m] 1.3 l 0 t [s] Fg.3 Fg.3 xy y θ[degree] l x [m] l y [m] l z [m] V[m/s] h[m] T Table 1 l 0 T
6 Measurement (mm) Pressure (Pa) Pressure [Pa] Vscosty Rato 10 Vscosty Rato 50 Vscosty Rato 100 Vscosty Rato 500 Vscosty Rato 1000 Theoretcal Value Smulaton Tme [s] Fg.7 Maxmum pressure of the lamnate durng pressed X Y Z X Fber Drecton Fg.4 Analytcal result (ν hor = ν ver 1000, t = T ) Y Vscosty Rato Fber Drecton Fg.5 Transverse plane of the analytcal result (t = T) Fg.6 60 Wdth 40 Length 20 Intal Length and Wdth Double Value of the Intal Measurement Vscosty Rato Measurement of the lamnate after pressed 3.2 Fg.4 t = T[s] 1000 Fg.5 t = T[s] xy ( z = h Vt/2 ) Fg.6 t = T[s] (x y ) 60.5mm () 121mm ( 2 ) P (x) (22) P (x) = 6µḣ h 3 ( L 2 x 2) (38) µ ḣ ( ) L x (38) Fg.7 (38) µ = ν ver ρ ḣ = V h = l z Vt L = l z L x / (2h) x = 0
7 (a) t = 0.0 (b) t = 0.5 (c) t = 1.0 (d) t = 1.5 (e) t = 2.0 (f) t = 2.5 Pressure [Pa] Z Fber Drecton (g) ColorMap X (h) Axes and Fver Drecton Fg.8 Coronal plane of the analytcal result.(ν hor = ν ver 1000)
8 Error Volume [m^3] Z Y Fg.9 X Mesh constructed wth partcles model (ν hor = ν ver 1000, t = T ) 8.0E E E E-05 Volume Intal Volume 0.0E Vsocsty Rato 4.6E E E E E E-02 Fg.10 Estmated volume of the lamnate Error 4.0E Vscosty Rato Fg.11 Error of the volume value ( V V nt /V nt ) Fg.8 t = 0.0 t = 0.5 t = 1.0 t = 1.5 t = 2.0 t = 2.5 xz z x y MPS Fg.4 Fg.8 Fg.4 Fg.8 LSMPS ( ) Delaunay Delaunay TetGen (23) 1000 Fg.9 Fg.10 Fg.11 V nt V V error V error = V V nt V nt (39) Fg.11 5 % 4. LSMPS
9 Standard LSMPS scheme Type-A (12) (1) Verweyst, B. E. and Tucker C. L.III, Fber Suspensons n Complex Geometres: Flow/Orentaton Couplng, The Canadan Journal of Chemcal Engneerng, 80-6, 2002, (2) Ishkawa, T. and Ogasawara, H., and Tomoka, M., Evaluaton of flowablty of thermoplastc carbon fber compostes for compresson moldngs, Proceedng of ECCM16 - European Conference on Composte Materals, Sevlle: Span, June (3) Lucy, L. B., A numercal approach to the testng of the fsson hypothess, Astronomcal Journal, 82, 1977, pp (4) Koshzuka, S. and Oka, Y., Movng-Partcle Sem-Implct Method for Fragmentaton of Incompressble Flud, Nuclear Scence and Engneerng, 123-3, 1996, pp (5) Xong, J., Koshzuka, S. and Saka, M., Numercal analyss of droplet mpngement usng the movng partcle sem-mplct method, Journal of Nuclear Scence and Technology, 47-3, 2010, pp (6) Shbata, K. and Koshzuka, S., Numercal analyss of shppng water mpact on a deck usng a partcle method, Ocean Engneerng, , 2007, pp (7), Smulaton-based 4DCT, Medcal Imagng Technology, 28-4, pp , (8),, Medcal Imagng Technology, 31-4, 2013, pp (9) Khayyer, A. and Gotoh, H., Modfed Movng Partcle Sem-mplct methods for the predcton of 2D wave mpact pressure, Coastal Engneerng, 56-4, 2009, pp (10), MPS,, 2014, Paper No (11),,, Movng Partcle Semmplct,, 2015, Paper No (12) Tama, T. and Koshzuka, S., Least squares movng partcle sem-mplct method, Computatonal Partcle Mechancs, 1-3, 2014, pp (13) Lee, C.-C. and Lee, F. Folgar and C. L. Tucker, Smulaton of Compresson Moldng for Fber-Renforced Thermosettng Polymers, Journal of Engneerng for Industry, 106-2, 1984, pp (14) Advan, S. G. and Tucker, C. L., The Use of Tensors to Descrbe and Predct Fber Orentaton n Short Fber Compostes, Journal of Rheology, 31-8, (15) Young, W.B., Resn Flow Analyss n the Consoldaton of Mult-Drectonal Lamnated Compostes, Polymer Compostes, 16-3, 1995, pp (16) Sun, X., Saka, M., Shbata, K., Tochg, Y., and Fujwara, H. Numercal modelng on the dscharged flud flow from a glass melter by a Lagrangan approach Nuclear Engneerng and Desgn, 248, 2012, pp (17), Khayyer Abbas,,, B2( ), 65-1, 2009, pp (18),,, MPS, B2( ), 68-1, 2012, pp (19) Lee, E.-S., Moulnec, C., Xu, R., Voleau, D., Laurence, D. and Stansbyc, P., Comparsons of weakly compressble and truly ncompressble algorthms for the SPH mesh free partcle method, Journal of Computatonal Physcs, , 2008, pp (20),,,, October (21) Nomura, K., Koshzuka, S., Oka, Y. and Obata, H., Numercal Analyss of Droplet Breakup Behavor usng Partcle Method, Journal of Nuclear Scence and Technology, 38-12, 2001, pp (22) 4, ( ), 2004.
10 (23) Hang, S., TetGen, a Delaunay-Based Qualty Tetrahedral Mesh Generator, Transactons on Mathematcal Software, 41-2, 2015, No.11.
粒子法による流れの数値解析
21 2002 230 239. Numercal Analyss of Flow usng Partcle Method Sech KOSHIZUKA 1 1 2 Los Alamos PAF Partcle-and-Force MAC Marker-and- Cell MAC PIC Partcle-n-Cell 319-1188 2-22 E-mal: kosh@utnl.jp PIC Los
More informationIPSJ SIG Technical Report Vol.2017-HPC-158 No /3/9 OpenACC MPS 1,a) 1 Moving Particle Semi-implicit (MPS) MPS MPS OpenACC GPU 2 4 GPU NVIDIA K2
OpenACC MPS 1,a) 1 Movng Partcle Sem-mplct (MPS) MPS MPS OpenACC GPU 2 4 GPU NVIDIA K20c GTX1080 P100(PCIe) P100(NVlnk) 5 OpenACC 3.5 3 Fortran 29.0 74.5 GPU 1. MPS [1] 1 MPS MPS CUDA GPU [2] [3] [4] OpenACC
More information> σ, σ j, j σ j, σ j j σ σ j σ j (t) = σ (t ) σ j (t) = σ () j(t ) n j σ, σ j R lm σ = σ j, j V (8) t σ R σ d R lm σ = σ d V (9) t Fg.. Communcaton ln
IIC-- Dstrbuted Cooperatve Atttude Control for Multple Rgd Bodes wth Communcaton Delay Yoshhro achbana, oru Namerkawa (Keo Unversty) Abstract hs paper descrbes dstrbuted cooperatve atttude consensus and
More information3 1, 1, 1, 1 3D Environment Measurement Using Binocular Stereo and Motion Stereo by Mobile Robot with Omnidirectional Stereo Camera Shinichi GOTO 1, R
3 1, 1, 1, 1 3D Envronment Measurement Usng Bnocular Stereo and Moton Stereo by Moble Robot wth Omndrectonal Stereo Camera Shnch GOTO 1, Ryosuke KAWANISHI 1, Atsush YAMASHITA 1 and Toru KANEKO 1 1 Department
More information橡jttc2.PDF
1 ( ) 1 GA GA GA MOGA (Multple-Objectve Genetc Algorthm) GA GA GA MOGA GA GA MOGA GA GA 3.1MOGA ( ) x x j f = f, f, 1 2 L, f q x x j x j f ( x ) f ( x ) f ( x ) f ( x ) L f ( x ) f ( x ) ( ) ( ) 1 1 j
More information非線形長波モデルと流体粒子法による津波シミュレータの開発 I_ m ρ v p h g a b a 2h b r ab a b Fang W r ab h 5 Wendland 1995 q= r ab /h a d W r ab h
土木学会論文集 B2( 海岸工学 ) Vol. 70, No. 2, 2014, I_016-I_020 非線形長波モデルと流体粒子法による津波シミュレータの開発 Development of a Tsunami Simulator Integrating the Smoothed-Particle Hydrodynamics Method and the Nonlinear Shallow Water
More informationA Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical
A Precise Calculation Method of the Gradient Operator in Numerical Computation with the MPS Tsunakiyo IRIBE and Eizo NAKAZA A highly precise numerical calculation method of the gradient as a differential
More information10_4.dvi
Vol.44, No.1, 1/7 28 Synchronzed Control for Blateral Teleoperaton wth Dfferent Confguratons and Communcaton Delays Hsanosuke Kawada,KoueYoshda and Toru Namerkawa Ths paper addresses the problem of the
More information第62巻 第1号 平成24年4月/石こうを用いた木材ペレット
Bulletin of Japan Association for Fire Science and Engineering Vol. 62. No. 1 (2012) Development of Two-Dimensional Simple Simulation Model and Evaluation of Discharge Ability for Water Discharge of Firefighting
More informationKENZOU Karman) x
KENZO 8 8 31 8 1 3 4 5 6 Karman) 7 3 8 x 8 1 1.1.............................. 3 1............................................. 5 1.3................................... 5 1.4 /.........................
More information7 OpenFOAM 6) OpenFOAM (Fujitsu PRIMERGY BX9, TFLOPS) Fluent 8) ( ) 9, 1) 11 13) OpenFOAM - realizable k-ε 1) Launder-Gibson 15) OpenFOAM 1.6 CFD ( )
71 特集 オープンソースの大きな流れ Nonlinear Sloshing Analysis in a Three-dimensional Rectangular Pool Ken UZAWA, The Center for Computational Sciences and E-systems, Japan Atomic Energy Agency 1 1.1 ( ) (RIST) (ORNL/RSICC)
More informationTitle 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue
Title 混合体モデルに基づく圧縮性流体と移動する固体の熱連成計算手法 Author(s) 鳥生, 大祐 ; 牛島, 省 Citation 土木学会論文集 A2( 応用力学 ) = Journal of Japan Civil Engineers, Ser. A2 (2017), 73 Issue Date 2017 URL http://hdl.handle.net/2433/229150 Right
More informationmain.dvi
THE INSTITUTE OF ELECTRONICS, INFORMATION AND COMMUNICATION ENGINEERS TECHNICAL REPORT OF IEICE. ( ) Estmaton and Analyss of Topc Models n Tme Seres Japanese / Chnese News and Blogs Shuo HU,LyZHENG, Yusuke
More informationJFE.dvi
,, Department of Civil Engineering, Chuo University Kasuga 1-13-27, Bunkyo-ku, Tokyo 112 8551, JAPAN E-mail : atsu1005@kc.chuo-u.ac.jp E-mail : kawa@civil.chuo-u.ac.jp SATO KOGYO CO., LTD. 12-20, Nihonbashi-Honcho
More informationIPSJ SIG Techncal Report 2. RangeBased RangeFree. 2.1 Rangebased RangeBased TDOA(Tme Dfference Of Arrval) TOA(Tme Of Arrval) TDOA TDOA Actve Bat 2) Cr
IPSJ SIG Techncal Report 1 2 2 (SOM, Self Organzng Maps) 7). Self-Organzng Localzaton for Wreless Sensor Networks on Ansotropc Topology Yuto Takashma, 1 Naotosh Adach 2 and Yasuhsa Takzawa 2 On wreless
More information混相流 Vol.31 No.4
論 文 Movng Partcle Sem-mplct 法による ノズル内キャビテーションと液体噴流のシミュレーション * Numercal Smulaton of Cavtaton n a Nozzle and Lqud Jet Usng Movng Partcle Sem-mplct Method 関 根 章 裕 ** SEKINE Akhro 吉 村 一 樹 *** YOSHIMURA Kazuk
More informationk m m d2 x i dt 2 = f i = kx i (i = 1, 2, 3 or x, y, z) f i σ ij x i e ij = 2.1 Hooke s law and elastic constants (a) x i (2.1) k m σ A σ σ σ σ f i x
k m m d2 x i dt 2 = f i = kx i (i = 1, 2, 3 or x, y, z) f i ij x i e ij = 2.1 Hooke s law and elastic constants (a) x i (2.1) k m A f i x i B e e e e 0 e* e e (2.1) e (b) A e = 0 B = 0 (c) (2.1) (d) e
More information変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy,
変 位 変位とは 物体中のある点が変形後に 別の点に異動したときの位置の変化で あり ベクトル量である 変位には 物体の変形の他に剛体運動 剛体変位 が含まれている 剛体変位 P(x, y, z) 平行移動と回転 P! (x + u, y + v, z + w) Q(x + d x, y + dy, z + dz) Q! (x + d x + u + du, y + dy + v + dv, z +
More information修士論文
SAW 14 2 M3622 i 1 1 1-1 1 1-2 2 1-3 2 2 3 2-1 3 2-2 5 2-3 7 2-3-1 7 2-3-2 2-3-3 SAW 12 3 13 3-1 13 3-2 14 4 SAW 19 4-1 19 4-2 21 4-2-1 21 4-2-2 22 4-3 24 4-4 35 5 SAW 36 5-1 Wedge 36 5-1-1 SAW 36 5-1-2
More information64 3 g=9.85 m/s 2 g=9.791 m/s 2 36, km ( ) 1 () 2 () m/s : : a) b) kg/m kg/m k
63 3 Section 3.1 g 3.1 3.1: : 64 3 g=9.85 m/s 2 g=9.791 m/s 2 36, km ( ) 1 () 2 () 3 9.8 m/s 2 3.2 3.2: : a) b) 5 15 4 1 1. 1 3 14. 1 3 kg/m 3 2 3.3 1 3 5.8 1 3 kg/m 3 3 2.65 1 3 kg/m 3 4 6 m 3.1. 65 5
More informationR927清水信彦様.indd
Special Issue CFRP 455-8502 9 1 Development Status of Carbon Fiber Reinforced Plastics Nobuhiko SHIMIZU Automotive Center, Toray Industries, Inc., 9-1 Oe-cho, Minato-ku, Nagoya, Aichi 455-8502 Received
More information1 n 1 1 2 2 3 3 3.1............................ 3 3.2............................. 6 3.2.1.............. 6 3.2.2................. 7 3.2.3........................... 10 4 11 4.1..........................
More informationタイトル
Flud Flow Smulton wth Cellulr Automt 00N2100008J 2002 225 1. cellulr utomton n t+ 1 t t = f( r, L, + r (1 t t+ 1 f t r t +1 Prllel Vrtul Mchne Messge-Pssng Interfce 1 2. 2. 1 t t 0 t = 1 = 1 = t = 2 2
More informationOptical Lenses CCD Camera Laser Sheet Wind Turbine with med Diffuser Pitot Tube PC Fig.1 Experimental facility. Transparent Diffuser Double Pulsed Nd:
*1 *2 *3 PIV Measurement of Field of the Wind Turbine with a med Diffuser Kazuhiko TOSHIMITSU *4, Koutarou NISHIKAWA and Yuji OHYA *4 Department of Mechanical Engineering, Matsue National Collage of Technology,
More informationkeisoku01.dvi
2.,, Mon, 2006, 401, SAGA, JAPAN Dept. of Mechanical Engineering, Saga Univ., JAPAN 4 Mon, 2006, 401, SAGA, JAPAN Dept. of Mechanical Engineering, Saga Univ., JAPAN 5 Mon, 2006, 401, SAGA, JAPAN Dept.
More informationII A A441 : October 02, 2014 Version : Kawahira, Tomoki TA (Kondo, Hirotaka )
II 214-1 : October 2, 214 Version : 1.1 Kawahira, Tomoki TA (Kondo, Hirotaka ) http://www.math.nagoya-u.ac.jp/~kawahira/courses/14w-biseki.html pdf 1 2 1 9 1 16 1 23 1 3 11 6 11 13 11 2 11 27 12 4 12 11
More information( ; ) C. H. Scholz, The Mechanics of Earthquakes and Faulting : - ( ) σ = σ t sin 2π(r a) λ dσ d(r a) =
1 9 8 1 1 1 ; 1 11 16 C. H. Scholz, The Mechanics of Earthquakes and Faulting 1. 1.1 1.1.1 : - σ = σ t sin πr a λ dσ dr a = E a = π λ σ πr a t cos λ 1 r a/λ 1 cos 1 E: σ t = Eλ πa a λ E/π γ : λ/ 3 γ =
More informationO x y z O ( O ) O (O ) 3 x y z O O x v t = t = 0 ( 1 ) O t = 0 c t r = ct P (x, y, z) r 2 = x 2 + y 2 + z 2 (t, x, y, z) (ct) 2 x 2 y 2 z 2 = 0
9 O y O ( O ) O (O ) 3 y O O v t = t = 0 ( ) O t = 0 t r = t P (, y, ) r = + y + (t,, y, ) (t) y = 0 () ( )O O t (t ) y = 0 () (t) y = (t ) y = 0 (3) O O v O O v O O O y y O O v P(, y,, t) t (, y,, t )
More information.5 z = a + b + c n.6 = a sin t y = b cos t dy d a e e b e + e c e e e + e 3 s36 3 a + y = a, b > b 3 s363.7 y = + 3 y = + 3 s364.8 cos a 3 s365.9 y =,
[ ] IC. r, θ r, θ π, y y = 3 3 = r cos θ r sin θ D D = {, y ; y }, y D r, θ ep y yddy D D 9 s96. d y dt + 3dy + y = cos t dt t = y = e π + e π +. t = π y =.9 s6.3 d y d + dy d + y = y =, dy d = 3 a, b
More informationhttp://www.ns.kogakuin.ac.jp/~ft13389/lecture/physics1a2b/ pdf I 1 1 1.1 ( ) 1. 30 m µm 2. 20 cm km 3. 10 m 2 cm 2 4. 5 cm 3 km 3 5. 1 6. 1 7. 1 1.2 ( ) 1. 1 m + 10 cm 2. 1 hr + 6400 sec 3. 3.0 10 5 kg
More informationIPSJ SIG Technical Report Vol.2014-ARC-213 No.24 Vol.2014-HPC-147 No /12/10 GPU 1,a) 1,b) 1,c) 1,d) GPU GPU Structure Of Array Array Of
GPU 1,a) 1,b) 1,c) 1,d) GPU 1 GPU Structure Of Array Array Of Structure 1. MPS(Moving Particle Semi-Implicit) [1] SPH(Smoothed Particle Hydrodynamics) [] DEM(Distinct Element Method)[] [] 1 Tokyo Institute
More informationDPA,, ShareLog 3) 4) 2.2 Strino Strino STRain-based user Interface with tacticle of elastic Natural ObjectsStrino 1 Strino ) PC Log-Log (2007 6)
1 2 1 3 Experimental Evaluation of Convenient Strain Measurement Using a Magnet for Digital Public Art Junghyun Kim, 1 Makoto Iida, 2 Takeshi Naemura 1 and Hiroyuki Ota 3 We present a basic technology
More information数値計算:有限要素法
( ) 1 / 61 1 2 3 4 ( ) 2 / 61 ( ) 3 / 61 P(0) P(x) u(x) P(L) f P(0) P(x) P(L) ( ) 4 / 61 L P(x) E(x) A(x) x P(x) P(x) u(x) P(x) u(x) (0 x L) ( ) 5 / 61 u(x) 0 L x ( ) 6 / 61 P(0) P(L) f d dx ( EA du dx
More informationac b 0 r = r a 0 b 0 y 0 cy 0 ac b 0 f(, y) = a + by + cy ac b = 0 1 ac b = 0 z = f(, y) f(, y) 1 a, b, c 0 a 0 f(, y) = a ( ( + b ) ) a y ac b + a y
01 4 17 1.. y f(, y) = a + by + cy + p + qy + r a, b, c 0 y b b 1 z = f(, y) z = a + by + cy z = p + qy + r (, y) z = p + qy + r 1 y = + + 1 y = y = + 1 6 + + 1 ( = + 1 ) + 7 4 16 y y y + = O O O y = y
More informationTERG
Dscusson Paper No. 268 小標本特性に優れたパネル単位根検定 千木良弘朗 山本拓 2011 年 7 月 TOHOKU ECONOMICS RESEARCH GROUP GRADUATE SCHOOL OF ECONOMICS AND MANAGEMENT TOHOKU UNIVERSITY KAWAUCHI, AOBA-KU, SENDAI, 980-8576 JAPAN Λ z
More informationIPSJ SIG Techncal Report 歌声データベース 歌声の波形 スペクトル抽出 基本周波数抽出 HMM メルケプストラム ラベル HMM の学習 対数基本周波数 c 学習部 コンテキスト依存モデル c ( 合成部 楽譜 ラベル変換 ラベル... メルケプストラム パラメータ生成 ML
IPSJ SIG Techncal Report HMM HMM HMM HMM 60 Vbrato Modelng for HMM-based Sngng Voce Synthess Tomohko Yamada, Satoru Muto, Yoshhko ankaku, Shnj Sako and Kech Tokuda HMM-based sngng voce synthess can automatcally
More information, 3, STUDY ON IMPORTANCE OF OPTIMIZED GRID STRUCTURE IN GENERAL COORDINATE SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO
, 3, 2012 9 STUDY ON IMPORTANCE OF OPTIMIZED GRID STRUCTURE IN GENERAL COORDINATE SYSTEM 1 2 Hiroyasu YASUDA and Tsuyoshi HOSHINO 1 950-2181 2 8050 2 950-2181 2 8050 Numerical computation of river flows
More information立命館21_松本先生.indd
2143-552010 1 2 () 1 2 3 456 78- Key Words 1 3 12007 2 3 508 19992008 1 23 4 43 2120107 4 2008 1992 2003 2005 1989 2008 4 2 2-1 10 4 4 6 5 4 5 2946 1 155 11 41 44 45 4 2-2 2003 1 21 2 1 3 4 5 6 7 3 2120107
More information立命館20_服部先生.indd
20 93-105 2010 2 () ' - 1 ( ) ( ' 2005) Key Words 2 1967 93 20 2010 1 94 1 ' 2002 2 2 1996 1996 1999 2 2 2 1993 1999 4 1 2 1985 1989 1986 1994 4 1 2 1 1 4 2 4 4 1 4 1966 4 10-1970 20 1993-1972 95 20 2010
More information1996 2000 2004 1984 2005 7150 000 9 500 9 4 13 10 95 11 11 12 20002004 9 70
14 2006 1 Key Words 2002 3 1 2 3 3 1 2 3 1969 1987 69 1996 2000 2004 1984 2005 7150 000 9 500 9 4 13 10 95 11 11 12 20002004 9 70 14 2006 1 15 71 72 1 22 6 32 9 200 6 3 1 2 2000 10 1 2003 10 2005 6 5 4
More information立命館16_坂下.indd
1669-792008 1 - ' 85- -108 ' Key words 1 2 2003 69 1620082 5 3 1990 1997 4 5 2001 1307 6 7 1 1 1 1996 100 2 1997 71 3 1998 71 4 1999 96 5 2000 95 6 2001 145 7 2002 191 8 2003 174 9 2004 120 10 2005 122
More information立命館人間科学研究No.10
1 77 5 Key words 1 23 3 11417 14310045 20022004 2 2003 20022005 20022004 2 10200511 3 2003 1152003 59 1995 3 32003 19932002 20032003 2005 20052005 4 1997 2000521986 2001 42001 3 1981 6 1 7 5 1000248 1632647
More information立命館21_川端先生.indd
21 119-132 2010 ( ) ' Key Words 119 21 2010 7 1962 2001 2001 2007 1982 1988 1997 2007 1997 1998 1863 1880 1 1998 1998 2001 1599 120 121 1599 1695 8 1695 1714 4 1714 1715 5 1715 100 1812 9 1812 1864 2001
More information立命館14_前田.indd
1499-1122007 1 ) ) ) ) ) (1) -- ) ) ) ) ) ) ( ) ) ' ) ) ) ) - 1) 2)' 3) Key words 19811994 1721 99 1420073 100 20012004 2005 2004 2002 33 34 10 1987 45 20002003 2002 1 1 2000 1 1 2001 2 200341 12004 2
More information立命館17_坂下.indd
1793-1052008 1 () -- -- - Key words 1 93 1720088 2 3 2003 15 1996 50 3040 2 3 4050 50 10 1980 1995 1950 1958 1968 1972 94 95 1987 4 4 70 3 3 1 2000 2001 4 1720088 96 2001 2003 2 1978 1990 1997 130 2 3
More information立命館人間科学研究No.10
49 00 7 Key words 980 995 50 0005 90 997 99 990 994 99 99 99990 99 996 988988 994 99 995 995 995 994 984 988 986 997 997 00 995 00 5 7 5 5 999 997 997 998 6 998 999 997 997 998 0040000 994 996 000 00 5
More information立命館19_椎原他.indd
191-132009 1 ( ) ' Key Words 1 11 12007 2007 200520062 201120062 7558 4009 1 1920098 2007 2 2000 028 2005 1999 12 1999 13 1968 2 3 '1992 2007 2001 20052001 1977 2005 2005 2007 21 21 22 461927 3 13 1920098
More information立命館人間科学研究No.10
61 10 1990 Key words 1 102005 11 62 2000 1 1920 10 1892 63 1 1 19 100 1914 100 1 1 20 1 2 102005 11 64 1946 21 4 1947 1949 1947 22 1 3 1956 1959 1958 2 1964 65 2 10 1975 7080 70 2 1 1987 1990 40 1989 1989
More information立命館19_徳田.indd
1991-1022009 1 2 () )--) 28 2827 1 2 Key Words 1 (1721 2000 20012001 20082009 91 1920098 92 20042004 2004 2001 12 2005 20082009 2005 3 1997 200820096 2007 2 20082009 93 20012001 12 2008 2009 19611966 1
More information176 B B.1: ( ) ( ) ( ) (2 2 ) ( ) ( ) ( ) (quantitative nondestructive evaluation:qnde) (1) X X X X CT(computed tomography)
B 1) B.1 B.1.1 ( ) B.1 1 50 100 m B.1.2 (nondestructive testing:ndt) (nondestructive inspection:ndi) (nondestructive evaluation:nde) 175 176 B B.1: ( ) ( ) ( ) (2 2 ) ( ) ( ) ( ) (quantitative nondestructive
More information( ) x y f(x, y) = ax
013 4 16 5 54 (03-5465-7040) nkiyono@mail.ecc.u-okyo.ac.jp hp://lecure.ecc.u-okyo.ac.jp/~nkiyono/inde.hml 1.. y f(, y) = a + by + cy + p + qy + r a, b, c 0 y b b 1 z = f(, y) z = a + by + cy z = p + qy
More information2009 2010 2 23 (MHD ) GFV (Galium Field Visualizer) GFV OpenGL GFV GFV GFV 1 1 2 2 2.1.................... 2 2.2................................. 2 2.3...................... 3 3 6 3.1 GFV....................
More informationVRSJ-SIG-MR_okada_79dce8c8.pdf
THE INSTITUTE OF ELECTRONICS, INFORMATION AND COMMUNICATION ENGINEERS TECHNICAL REPORT OF IEICE. 630-0192 8916-5 E-mail: {kaduya-o,takafumi-t,goshiro,uranishi,miyazaki,kato}@is.naist.jp,.,,.,,,.,,., CG.,,,
More information7) FOE (Maxmum A Posteror : MAP) MAP 2. ( ) F (1) 2) p(f ; θ) = 1 Z all clques = 1 Z exp [ φ(f ; θ) all clques λ(f ; θ) ] (1) F f φ( ) λ( ) θ Z
1 2 2 (MRF) Natural Image pror model usng adaptve mult-varate Gaussan dstrbuton Ketaro Yamauch, 1 Masayuk Tanaka 2 and Masatosh Okutom 2 Many pror models are proposed whch model a dgtal mage by a Markov
More informationwakate2005-tobari-ver2.ppt
8 005 3 17 Outlne 1. Introducton. 3. 4. Introducton Image of MUSES-C on engne JAXA b = 4ε 9 q m V d 3 a large-scale) MHD MPD chokng VASIMR Super-Alfvénc flow plasma detachment HITOP(HIgh densty TOhoku
More informationMicrosoft Word - 予稿原稿2(中口).doc
社団法人電子情報通信学会 THE INSTITUTE OF ELECTRONICS, INFORMATION AND COMMUNICATION ENGINEERS 信学技報 TECHNICAL REPORT OF IEICE 複数の Level Set を用いた CT 画像からの膵臓抽出と術前シミュレーションシステムへの応用 中口俊哉 奥井雅博 津村徳道 Ken Museth 三宅洋一 千葉大学工学部
More informationA Higher Weissenberg Number Analysis of Die-swell Flow of Viscoelastic Fluids Using a Decoupled Finite Element Method Iwata, Shuichi * 1/Aragaki, Tsut
A Higher Weissenberg Number Analysis of Die-swell Flow of Viscoelastic Fluids Using a Decoupled Finite Element Method Iwata, Shuichi * 1/Aragaki, Tsutomu * 1/Mori, Hideki * 1 Ishikawa, Satoshi * 1/Shin,
More information[Ver. 0.2] 1 2 3 4 5 6 7 1 1.1 1.2 1.3 1.4 1.5 1 1.1 1 1.2 1. (elasticity) 2. (plasticity) 3. (strength) 4. 5. (toughness) 6. 1 1.2 1. (elasticity) } 1 1.2 2. (plasticity), 1 1.2 3. (strength) a < b F
More informationIPSJ SIG Technical Report Vol.2013-MUS-100 No /9/1 1,a) [1 9] 1 3 MIDI MIDI (Interonset interval) 1 NTT NTT Communication Science Labo
1,a) 1 1 1. [1 9] 1 3 MIDI MIDI (Interonset nterval) 1 N N Communcaton Scence Laboratores, Atsug, Kanagawa 43 0198, Japan he Insttute of Statstcal Mathematcs, achkawa, okyo 190 856 a) ohsh.yasunor@lab.ntt.co.jp
More information1
PalmGauss SC PGSC-5G Instruction Manual PalmGauss SC PGSC-5G Version 1.01 PalmGauss SC PGSC5G 1.... 3 2.... 3 3.... 3 3.1... 3 3.2... 3 3.3 PalmGauss... 4 3.4... 4 3.4.1 (Fig. 4)... 4 3.4.2 (Fig. 5)...
More information[ 30 p. 1-8 (2012)] / ** *** Numerical Analysis of Metal Transfer Phenomena - critical condition between globular and spray transfer mode - by KADOTA
[ 30 p. 1-8 (2012)] / ** *** Numerical Analysis of Metal Transfer Phenomena - critical condition between globular and spray transfer mode - by KADOTA Keiji and HIRATA Yoshinori Metal transfer modes in
More informationReport2009_Softsensor.dvi
/ http://www-pse.cheme.kyoto-u.ac.jp/ kano/ 2009 03 Copyright c 2009 Manabu Kano. All rights reserved. / 1 1 soft-sensor NIR Process Analytical Technology PAT NIR / 2 2 2.1 x m (m =1, 2,,M) 1 y M y = a
More informationシリコン結晶化過程の分子動力学
1-65 13 16 96154 1 4 1.1 5 1.1.1 5 1.1. 5 1.1.3 5 1.1.4 CVD 7 1. 8 1..1 8 1.. SPE 10 1.3 11 1.1 13. 14..1 Tesoff 14.. Lennad-Jones 17.3 18.4 19.5 0.6 Langevn.7 3.7.1 SPE 3.7. 4.7.3 6 3 SPE 8 3.1 9 3. [111]
More information21 2 26 i 1 1 1.1............................ 1 1.2............................ 3 2 9 2.1................... 9 2.2.......... 9 2.3................... 11 2.4....................... 12 3 15 3.1..........
More information1 M = (M, g) m Riemann N = (N, h) n Riemann M N C f : M N f df : T M T N M T M f N T N M f 1 T N T M f 1 T N C X, Y Γ(T M) M C T M f 1 T N M Levi-Civi
1 Surveys in Geometry 1980 2 6, 7 Harmonic Map Plateau Eells-Sampson [5] Siu [19, 20] Kähler 6 Reports on Global Analysis [15] Sacks- Uhlenbeck [18] Siu-Yau [21] Frankel Siu Yau Frankel [13] 1 Surveys
More information9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L) du (L) = f (9.3) dx (9.) P
9 (Finite Element Method; FEM) 9. 9. P(0) P(x) u(x) (a) P(L) f P(0) P(x) (b) 9. P(L) 9. 05 L x P(x) P(0) P(x) u(x) u(x) (0 < = x < = L) P(x) E(x) A(x) P(L) f ( d EA du ) = 0 (9.) dx dx u(0) = 0 (9.2) E(L)A(L)
More information空気の屈折率変調を光学的に検出する超指向性マイクロホン
23 2 1M36268 2 2 4 5 6 7 8 13 15 2 21 2 23 2 2 3 32 34 38 38 54 57 62 63 1-1 ( 1) ( 2) 1-1 a ( sinθ ) 2J D ( θ ) = 1 (1-1) kaka sinθ ( 3) 1-2 1 Back face hole Amplifier Diaphragm Equiphase wave surface
More informationThe Engneer Grapples wth Nonlnear Problems
The Engneer Grapples wth Nonlnear Problems Hstory of Computatonal Engneerng CAE Computer Aded Engneerng NASTRAN Factory Automaton Computatonal Mechancs IT Informaton Technology What s Computatonal Mechancs?
More informationEvidence for jet structure in hadron product by e+e-
G. Hanson et al. Phys. Rev. Lett. 5 (1975) 1609 Physcs Colloquum July 7th, 008 Evdence for Jet Structure n Hadron Producton by e + e - Annhlaton Contents: 1. Introducton. Exerment at SLAC. Analyss 4. Results
More information2. 2 P M A 2 F = mmg AP AP 2 AP (G > : ) AP/ AP A P P j M j F = n j=1 mm j G AP j AP j 2 AP j 3 P ψ(p) j ψ(p j ) j (P j j ) A F = n j=1 mgψ(p j ) j AP
1. 1 213 1 6 1 3 1: ( ) 2: 3: SF 1 2 3 1: 3 2 A m 2. 2 P M A 2 F = mmg AP AP 2 AP (G > : ) AP/ AP A P P j M j F = n j=1 mm j G AP j AP j 2 AP j 3 P ψ(p) j ψ(p j ) j (P j j ) A F = n j=1 mgψ(p j ) j AP
More informationSpacecraft Propulsion Using Solar Energy Spacecraft with Magnetic Field Light from the Sun Solar Wind Thrust Mirror Solar Sail Thrust production by li
2004.3.28 物理学会シンポジウム 磁気プラズマセイル の可能性と 深宇宙探査への挑戦 宇宙航空研究開発機構 船木一幸 Spacecraft Propulsion Using Solar Energy Spacecraft with Magnetic Field Light from the Sun Solar Wind Thrust Mirror Solar Sail Thrust production
More information1.... 3 2.... 5 3.... 8 3.1.... 9 3.2.... 13 3.3.... 18 4.... 21 4.1.... 21 4.2.... 23 4.3.... 26 5.... 32... 33 1... 35 2... 39 1.... 39 2.... 43 2
JILPT Dscusson Paper Seres 12-2 212 3 ( 22 6 18 ) 1 1.... 3 2.... 5 3.... 8 3.1.... 9 3.2.... 13 3.3.... 18 4.... 21 4.1.... 21 4.2.... 23 4.3.... 26 5.... 32... 33 1... 35 2... 39 1.... 39 2.... 43 2
More information第5章 偏微分方程式の境界値問題
October 5, 2018 1 / 113 4 ( ) 2 / 113 Poisson 5.1 Poisson ( A.7.1) Poisson Poisson 1 (A.6 ) Γ p p N u D Γ D b 5.1.1: = Γ D Γ N 3 / 113 Poisson 5.1.1 d {2, 3} Lipschitz (A.5 ) Γ D Γ N = \ Γ D Γ p Γ N Γ
More informationAuerbach and Kotlikoff(1987) (1987) (1988) 4 (2004) 5 Diamond(1965) Auerbach and Kotlikoff(1987) 1 ( ) ,
,, 2010 8 24 2010 9 14 A B C A (B Negishi(1960) (C) ( 22 3 27 ) E-mail:fujii@econ.kobe-u.ac.jp E-mail:082e527e@stu.kobe-u.ac.jp E-mail:iritani@econ.kobe-u.ac.jp 1 1 1 2 3 Auerbach and Kotlikoff(1987) (1987)
More information平成12年度
1 1-1 (1) 150[ml] 500[ml/] Cerebral Ventricle Brain 1-1 2 ( ) 1-1 1-2 0.20.5[mm] 13 14[mm] 1-2 3 ( ) (2) 4 2-1 (cerebral ventricle) (peritoneum) R O p O Cerebral Ventricle Valve Brain R o R i P i Peritoneum
More informationUntitled
II 14 14-7-8 8/4 II (http://www.damp.tottori-u.ac.jp/~ooshida/edu/fluid/) [ (3.4)] Navier Stokes [ 6/ ] Navier Stokes 3 [ ] Reynolds [ (4.6), (45.8)] [ p.186] Navier Stokes I 1 balance law t (ρv i )+ j
More information28 Horizontal angle correction using straight line detection in an equirectangular image
28 Horizontal angle correction using straight line detection in an equirectangular image 1170283 2017 3 1 2 i Abstract Horizontal angle correction using straight line detection in an equirectangular image
More information14 2 5
14 2 5 i ii Surface Reconstruction from Point Cloud of Human Body in Arbitrary Postures Isao MORO Abstract We propose a method for surface reconstruction from point cloud of human body in arbitrary postures.
More informationJKR Point loading of an elastic half-space 2 3 Pressure applied to a circular region Boussinesq, n =
JKR 17 9 15 1 Point loading of an elastic half-space Pressure applied to a circular region 4.1 Boussinesq, n = 1.............................. 4. Hertz, n = 1.................................. 6 4 Hertz
More information土木学会構造工学論文集(2011.3)
Vol.57A (11 3 ) RC Consecutve falling-weight impact test of large-scale RC girders under specified total input-impact energy * ** *** **** ***** Norimitsu Kishi, Hisashi Konno, Satoru Yamaguchi, Hiroshi
More informationJournal of Textile Engineering, Vol.53, No.5, pp
ORIGINAL PAPER Journal of Textile Engineering (2007), Vol.53, No.5, 203-210 2007 The Textile Machinery Society of Japan Analyzing the Path and the Tension of a Yarn under a False-twist Process Using a
More informationJ. Jpn. Inst. Light Met. 65(6): 224-228 (2015)
65 62015 224 228 ** Journal of The Japan Institute of Light Metals, Vol. 65, No. 6 (2015), 224 228 2015 The Japan Institute of Light Metals Investigation of heat flow behavior on die-casting core pin with
More informationuntitled
Stacking sequence optimization of composite wing using fractal branch and bound method Orbiting Plane : HOPE-X (JAXA) Fractal Branch and Bound Method (FBBM) Fractal structure of design space 5 y V a 9º
More information空力騒音シミュレータの開発
41 COSMOS-V, an Aerodynamic Noise Simulator Nariaki Horinouchi COSMOS-V COSMOS-V COSMOS-V 3 The present and future computational problems of the aerodynamic noise analysis using COSMOS-V, our in-house
More informationII No.01 [n/2] [1]H n (x) H n (x) = ( 1) r n! r!(n 2r)! (2x)n 2r. r=0 [2]H n (x) n,, H n ( x) = ( 1) n H n (x). [3] H n (x) = ( 1) n dn x2 e dx n e x2
II No.1 [n/] [1]H n x) H n x) = 1) r n! r!n r)! x)n r r= []H n x) n,, H n x) = 1) n H n x) [3] H n x) = 1) n dn x e dx n e x [4] H n+1 x) = xh n x) nh n 1 x) ) d dx x H n x) = H n+1 x) d dx H nx) = nh
More information