main.dvi

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "main.dvi"

Transcription

1 Copyright c 2008 Tech Science Press CMES, vol.28, no.2, pp , 2008 Meshless Generalized Finite Difference Method and Human Carotid Atherosclerotic Plaque Progression Simulation Using Multi-Year MRI Patient-Tracking Data Chun Yang 1, Dalin Tang 2, Chun Yuan 3, William Kerwin 2,FeiLiu 3, Gador Canton 3 Thomas S. Hatsukami 3,4, Satya Atluri 5 Abstract: Atherosclerotic plaque rupture and progression have been the focus of intensive investigations in recent years. Plaque rupture is closely related to most severe cardiovascular syndromes such as heart attack and stroke. A computational procedure based on meshless generalized finite difference (MGFD) method and serial magnetic resonance imaging (MRI) data was introduced to quantify patient-specific carotid atherosclerotic plaque growth functions and simulate plaque progression. Participating patients were scanned three times (T 1, T 2,andT 3,atintervals of about 18 months) to obtain plaque progression data. Vessel wall thickness (WT) changes were used as the measure for plaque progression. Since there was insufficient data with the current technology to quantify individual plaque component growth, the whole plaque was assumed to be uniform, homogeneous, hyperelastic, isotropic and nearly incompressible. The linear elastic model was used. The 2D plaque model was discretized and solved using a meshless generalized finite difference (GFD) method. Starting from the T 2 plaque geometry, plaque progression was simulated by solving the solid model and adjusting wall thickness using plaque growth functions iteratively until T 3 is reached. Numerically simulated plaque progression agreed very well with actual 1 School of Math, Beijing Normal University, Beijing, China 2 Corresponding author, Math Dept, WPI, Worcester, MA USA 3 Deparment of Radiology, University of Washington, Seattle, WA Division of Vascular Surgery, University of Washington, Seattle, WA Center for Aerospace Research & Education, University of California, Irvine, CA 92612, USA. plaque geometry at T 3 given by MRI data. We believe this is the first time plaque progression simulation based on multi-year patient-tracking data was reported. Serial MRI-based progression simulation adds time dimension to plaque vulnerability assessment and will improve prediction accuracy for potential plaque rupture risk. Keyword: meshless, generalized finite difference, artery, plaque progression, plaque rupture, atherosclerosis. 1 Introduction Cardiovascular disease (CVD) has remained the number one cause of death in America since % of 45 year olds and 80% of those 75 and older have CVD to some degree. A large number of the fatal clinical events are caused by rupture of a vulnerable atherosclerotic plaque. [Fuster (1998); Fuster et al. (1990); Naghavi et al. (2003a, 2003b)]. Many victims of the disease who are apparently healthy die suddenly without prior symptoms. While major advancements in the treatment of CVD continue, progress has been very limited in early detection and treatment of at risk individuals. Our long-term goal is to develop non-invasive methods to assess plaque vulnerability and predict possible rupture before the fatal event actually happens. There has been considerable effort investigating mechanisms governing atherosclerotic plaque progression and rupture [Friedman, Bargeron, Deters, Hutchins and Mark (1987); Friedman and Giddens (2005); Giddens, Zarins, Glagov, S. (1993); Ku, Giddens, Zarins and Glagov (1985); Scotti et al. (2005); Yuan, Mitsumori, Beach,

2 96 Copyright c 2008 Tech Science Press CMES, vol.28, no.2, pp , 2008 and Maravilla (2001)]. Most efforts were focused on fluid dynamics side since it has been well accepted that atherosclerosis initiation and progression correlate positively with low and oscillating flow wall shear stresses. However, this low and oscillating shear stress hypothesis cannot explain why moderate and advanced plaques continue to grow under elevated flow shear stress conditions [Tang et al. (2005)]. Our recent results using serial MRI patient-tracking data and computational models ( data points/patient) indicated that 18 out of 21 patients studied showed significant negative correlation between plaque progression measured by wall thickness increase and plaque wall (structure) stress. However, computational models using patient-specific plaque progression data to simulate plaque growth and predict future plaque rupture risk are lacking in the current literature. In this paper, a computational procedure based on meshless generalized finite difference (MGFD) method and serial MRI data is introduced to quantify patient-specific carotid atherosclerotic plaque growth functions and simulate plaque progression. By adding time dimension into our play, plaque vulnerability assessment and clinical decisions can be based on multi-time MRI scans and simulated virtual plaque progression. Seeing is believing. With validation, our procedure can be implemented in clinical applications and will lead to considerable improvement in prediction accuracy for potential plaque rupture risk and possible prevention of fatal cardiovascular events. Because of the complexity of plaque geometry and structure, a meshless GFD method is introduced in this paper to avoid frequent mesh updates and adjustments. Computational modeling for engineering applications with meshless methods have made considerable advances in recent years [Ahrem, Beckert and Wendland (2006); Atluri (2004, 2005); Atluri, Yagawa, and Cruse (1995); Bathe (1996, 2002); Ling, Atluri (2006); Shu, Ding, and Yeo (2005); Wu, Shen, Tao (2007)]. A series of meshless local Petrov- Galerkin (MLPG) methods were introduced to solved 3-dimensional elasto-static and dynamical problems [Han and Atluti (2004a, 2004b, 2007); Han, Liu, Rajendran, Atluri (2006)] and nonlinear problems with large deformation and rotations [Han, Rajendran and Atluri (2004)]. A mixed approach was introduced to improve the MLPG method using finite volume method [Atluri, Han and Rajendran (2004)] and finite difference method [Atluri, Liu, and Han, (2006a, 2006b); Hu, Young, Fan (2008)]. A new meshless interpolation scheme for MLPG method was developed by Ma [Ma (2008)]. Analysis of structure with material interfaces was performed by Masuda and Noguchi [Masuda, Noguchi (2006)]. Perko and Sarler studied weight function shape parameter optimization in meshless methods for non-uniform grids [Perko and Sarler (2007)]. Remeshing and refining with moving finite elements were investigated by Wacher and Givoli [Wacher, Givoli (2006)]. Numerical methods were also developed to solve problems with free and moving boundaries [Zohouri, Pirooz, and Esmaeily (2005); Mai-Duy and Tran-Cong, (2004)]. GFD methods have been used in many engineering applications and in our previous papers where irregular geometries and free-moving boundaries are involved [Kleiber (1998); Liszka and Orkisz (1980); Tang, Chen, Yang, Kobayashi and Ku (2002); Tang, Yang, Kobayashi and Ku (2001)]. One advantage of using GFD is that generalized finite difference schemes can be derived for userselected irregular grid points which can be freely adjusted to accommodate plaque deformation and growth. The MGFD method introduced in this paper uses grid points from the local support of each nodal point so that theoretical MLPG framework can be applied [Atluri (2004)]. Details are given in the following sections. 2 Models and methods Due to the complexity of the problem, we start from 2D models in this paper to get some insight for further full 3D investigations. Patientspecific plaque progression data was acquired by serial MRI (scanning patients multiple times with time span at about 18 months). Corresponding slices were matched using carotid bifurcation as the registration point. 2D models were constructed for selected slices and solved by mesh-

3 Finite Difference Method and Plaque Progression Simulation 97 less GFD method to obtain the stress distributions in the plaque. Point-wise plaque growth functions were determined based on three time point stress and vessel wall thickness data. The growth functions were then used to simulate plaque progression starting from Time 2 plaque morphology. The simulated Time 3 plaque morphology was compared with actual Time 3 plaque geometry to validate our modeling method. Details are given below. 2.1 In vivo Serial MRI Data Acquisition Serial MRI data from one patient was provided by the University of Washington (UW) group using protocols approved by the University of Washington Institutional Review Board with informed consent obtained. Scan time intervals were about 18 months, subject to scheduling variations. MRI scans were conducted on a GE SIGNA 1.5-T whole body scanner using an established protocol outlined in the papers by Yuan and Kerwin et al. [Kerwin, Hooker, Spilker, Vicini, Ferguson, Hatsukami, and Yuan (2003); Yuan and Kerwin (2004)]. Upon completion of a review, an extensive report was generated and segmented contour lines for different plaque components for each slice were sent to Tang s group for model construction and further computational mechanical analysis. Details of the model construction process can be found from [Yang et al. (2007); Tang et al. (2008)]. Figure 1 shows 5 (selected from 12) MRI slices with 5 different weightings obtained from a human carotid plaque sample. Figure 2 gives the re-constructed 3D geometries of the plaque at three time points showing plaque progression. Figure 3 gives three-time point segmented contour plots of 5 selected slices. These slices were used for model construction, quantification of growth function, and validation of simulated plaque progressions. 2.2 The structure model Since there was insufficient data to quantify individual plaque component growth, the plaque was treated as a uniform material, which was assumed to be hyperelastic, isotropic, incompressible and homogeneous. The governing equations and cor- (a) MRI T1-Weighting (b) T2 (c) Proton Density (PD) (d) TOF (Time of Flight) (e) CET1 (Contrast Enhanced T1) Figure 1: Selected multi-weighting MRI slices (5 of 12 slices) of carotid plaque from a participating patient. Multi-weighting MRI techniques can better differentiate plaque components and provide more accurate plaque vulnerability assessments [Yuan and Kerwin (2004)]. responding initial and boundary conditions are given below [Fung (1994)]: ρu i,tt = σ ij, j, i, j = 1,2; sum over j, (1) ε ij =(u i, j +u j,i )/2, i, j = 1,2,3, (2) σ ij n j out_wall = 0, (3) σ ij n j Γ = p in (t) Γ, (4) u i t=0 = u i0, (5) u i,t t=0 = u i0 (6) where ρ is material density, u =(u 1, u 2 ) is the displacement vector, σ isstresstensor, ε is strain tensor, P in is the specified lumen pressure, Γ is vessel inner boundary, f, j stands for derivative of f with respect to the jth variable. Using a linear model, the strain-stress relationship

4 98 Copyright c 2008 Tech Science Press CMES, vol.28, no.2, pp , 2008 (a) T 1, Front (c) T 2, Front (e) T 3, Front (b) T 1, Back (d) T 2, Back (f) T 3, Back Figure 2: Re-constructed 3D geometry of a carotid plaque based on in vivo serial MRI data. Three time point data are shown. T 1, T 2 and T 3 refer to time points from here on, unless otherwise indicated. T 1 T 2 : 525 days; T 2 T 3 : 651 days. Red lumen; Cyan- outer wall; Yellow - necrotic core; Fresh red - hemorrhage in necrotic core); Purple - loose matrix; Dark blue calcification; Green - fibrous cap. (a) Time 1 (b) Time 2 (c) Time cm Figure 3: Segmented contour plots using CAS- CADE showing plaque components and plaque progression. 5 slices were selected with the bifurcation serving as the registration point. Data from three time points are shown. Red lipid; blue calcification (Ca). is [Fung (1994)]: σ 11 ε 11 σ 22 = C ε 22 (7) σ 12 γ 12 C = E 1 μ μ0 2 μ C 1 C 2 0 = C 2 C C 3 1 μ 0 2 (8) where E 0 is the Young s Modulus, μ 0 is the Poisson ratio, C 1, C 2, and C 3 are coefficients defined by (8) for convenience. Substituting (2) and (7)- (8) into (1), we have the displacement equations in scalar form: ρ 2 u 1 t 2 = C 2 u 1 2 u 1 1 x 2 +C 3 1 x 2 +(C 2 +C 3 ) 2 u 2, 2 x 1 x 2 (9) ρ 2 u 2 t 2 = C 2 u 2 2 u 2 3 x 2 +C 1 1 x 2 +(C 2 +C 3 ) 2 u 1, 2 x 1 x 2 (10) The boundary conditions are: ( u 1 u 2 u1 C 1 n 1 +C 2 n 1 +C 3 n 2 + u ) 2 = t 1, x 1 x 2 x 2 x 1 (11) ( u 1 u 2 u1 C 2 n 2 +C 1 n 2 +C 3 n 1 + u ) 2 = t 2, x 1 x 2 x 2 x 1 (12) where (n 1, n 2 ) is the normal direction of the boundary and (t 1,t 2 )=(P in n 1,P in n 2 ) is the fluid force applied at the inner boundary (lumen), and the outer boundary is treated as a free boundary, with (t 1,t 2 )= The meshless GFD method The advantage of MGFD method is that generalized finite difference schemes can be derived using arbitrarily distributed points (see Fig. 4). With

5 Finite Difference Method and Plaque Progression Simulation 99 MGFD, we will be able to use denser nodal point distributions where plaque has higher stress/strain concentration or critical morphological features. We will also be able to adjust, move, add or drop nodal points as needed. This leads to the desired flexibility in handling the complex geometry and plaque growth where not only the geometry changes, the total plaque area (volume if 3D) also changes. The GFD concept is explained by the following example. Fig. 4(a) gives a selected nodal point Pi, its round support and all surrounding points Z j, j = 1,...,Ni, which are used to derive the GFD scheme. To derive the second order GFD schemes for the first and second order derivatives of the unknown function f (x 1,x 2 ), we use the Taylor expansion of f at Pi. Omitting higher order terms, we have, f f j = f i +h j f x +k j 1 i x + h2 j 2 i 2 x 2 1 i +h j k j x 1 x + k2 j 2 i 2 +o(h 2 j +k2 j ) (13) i x 2 2 where the subscript i indicates the selected node Pi at which derivative GFD schemes are being derived, j = 1...Ni are the other nodes in the Ninode-star of Pi, f j is the function value at Zj, Zj are the neighboring points of Pi within the support, h j = x 1 j x 1i,k j = x 2 j x 2i. Ignore the higher order terms in (13), we have:. h 1 k 1 h h 1 k 1 k h h Ni k 2 Ni kni Ni 2 h Ni k 2 Ni 2.. = f x 1 f x 2 x 2 1 x 1 x 2 x 2 2 f 1 f i. f Ni f i i (14) Rewrite (14) as A f = g, and use the least-squares method to get the formula for f : f =(A T WA) 1 A T W g (15) (a) A Nodal-Star with Round Support x 2 Z1 Z6 Pi Z7 Z10 Z11 Z12 Z14 Z15 x 1 (b) Nodal Points for Slice 4, Time 2. Z2 Z3 Z4 Z5 Z8 Z9 Z13 Z14 Z16 (c) Nodal Points for Slice 4, Time 2, with Components. Lipid Figure 4: Meshless GFD scheme derivation and nodal point distributions. (a) A schematic plot illustrating the derivation process of meshless GFD schemes with round support; (b)-(c) Nodal points distributions on Slice 4, Time 2, with and without plaque components. where W is the weight matrix with w 1 0 W =... 0 w Ni Let w j = B =(b kl ) 5 Ni =(A T WA) 1 A T W Ca 1 (h 2 j +k2 j ) j = 1,2,...Ni, (16)

6 100 Copyright c 2008 Tech Science Press CMES, vol.28, no.2, pp , 2008 we have, f x 1 f x 2 x 2 1 x 1 x 2 x 2 2 = i ni j=1 ni j=1 ni j=1 ni j=1 ni j=1 b 1 j ( f j f i ) b 2 j ( f j f i ) b 3 j ( f j f i ) b 4 j ( f j f i ) b 5 j ( f j f i ) (17) Substituting (17) into (9)-(10), we get the discrete approximation for inner nodes. For boundary nodes, we do a similar work using 1 st order GFDM method. That is: f f j = f i +h j f x +k j 1 i x +o(h j +k j ) (18) 2 i h 1 k 1.. h Ni k Ni ( ) f x 1 f = x 2 i f 1 f i. (19) f Ni f i b 1 j ( f j f i ) (20) b 2 j ( f j f i ) ( ) f x 1 f = x 2 i ni j=1 ni j=1 Substituting (20) into (11)-(12), we get the discrete approximation for boundary nodes. Second order scheme for boundary nodes were also tested. However, it did not give good result. One order lower schemes for boundary nodes are common in computational schemes. Time derivatives are discretized using the 2 nd order center difference scheme: t 2 = f n+1 2 f n + f n 1 Δt 2 +O(Δt 2 ), (21) Assembling all the discretized equations and boundary conditions, we get the final linear system: K u = f, (22) The vector ( ) u1 u = u 2 is the displacement solution at time step (n+1). Equation (22) is solved by a sparse solver under Matlab environment. 2.4 The shrink-pressurize process and MGFD model validation When constructing in vivo MRI-based plaque models, a shrink-pressurize process needs to be used to (a) shrink the original in vivo MRI geometry to get the numerical starting geometry with zero lumen pressure and then (b) pressurize the reduced starting geometry to recover the original in vivo geometry with specified lumen pressure. Using Slice 4 (S4) at Time 2 (T 2 )as an example, Fig. 5 shows its original in vivo MRI geometry, the numerical starting geometry with 16.75% inner boundary shrinkage and 6.5% outer-boundary shrinkage, and the pressurized geometry obtained by solving the GFD plaque model. The outer boundary was reduced less so that the conservation law of mass (area for 2D models) is enforced. Young s modulus was set at E 0 = dyn/cm 2, based on our experimental data [Kobayashi, Tsunoda, Fukuzawa, Morikawa, Tang, Ku (2003); Tang et al. (2008); Tang, Yang, Zheng, Woodard, Saffitz, Petruccelli, Sicard and Yuan (2005)] and current literature [Fung (1993); Humphrey (2002)]. Patient-specific pressure = 136 mmhg (by arm) was used as the lumen pressure. A commercial finite element software package ADINA (ADINA R & D, Inc., Watertown, MA) was used to validate our MGFD model. ADINA has been validated by hundreds of realistic engineering and real life applications and is well accepted in the industry and research communities [Bathe (1996); Bathe (2002)]. We have been using ADINA in the past 10 years to construct and solve 2D/3D artery models which were validated by experimental measurements [Tang et al. (2005,2008)]. A finite element ADINA model was constructed using S4 geometry and following the same procedures using in [Tang et al. (2005)]. Figure 6 compares maximum principal stress (Stress-P 1 ) from MGFD and ADINA models and shows that results from both models had very good agreement (error < 3%).

7 Finite Difference Method and Plaque Progression Simulation 101 a) Original MRI Geometry (b) Numerical Starting Geometry c) In vivo MRI Compared with Starting Geometry (d) In vivo MRI Compared with Pressurized Geometry Figure 5: Comparison of (a) original in vivo MRI plaque geometry with (b) numerical starting geometry, inner boundary shrinkage, 16.75%, outerboundary shrinkage 6.5%, and (d) the pressurized geometry obtained by solving the GFD plaque model, Slice 4 at T 2 was used. 3.1 A piecewise equal-step method to define wall thickness Vessel wall thickness was selected as the measure for plaque progression. In our previous paper, the shortest distance method was used to determine vessel thickness, i.e., for a selected nodal point on the inner boundary (lumen), the shortest distance between that point and the out-boundary was defined at the vessel thickness at that lumen point. That led to uneven selection of nodal points from the out boundary as shown by Fig. 7(a). A piecewise equal-step method is introduced to fix the problem. The vessel is divided into several pieces (segments) according to its geometry (4 in the S4 case). For each piece, equal step is used for inner and outer boundaries respectively to choose equal number of nodal points. The corresponding points on the inner and out boundaries are paired and the distance between the paired points are defined as vessel wall thickness at the given lumen point (Fig. 7). This method is sufficient for the cases covered in this paper. (a) Stress-P 1 by GFD (b) Stress-P 1 by ADINA * Min=0.00 KPa * (a) Shortest Distance Method (b) Piecewise Equal-Step Method Δ Max=105.2 KPa Δ Max=108.6 KPa Min Universal Scale Max Figure 6: Comparison of maximum principal stress (Stress-P 1 ) by MGFD and ADINA for model validation. Figure 7: Piecewise equal-step method for determination of vessel wall thickness. Slice 4 at T2 is used for demonstration. The vessel wall was divided into 4 segments. 25 points were equally distributed on each segment. 3 Plaque growth function and progression simulation With the MGFD model validated, we can use it to derive plaque growth function and simulate plaque progression. 3.2 Quantifying plaque growth function (PGF) using serial MRI data Plaque growth function (PGF) is a function we use to determine the vessel wall thickness or nodal point displacement for every numerical time step based on current and past plaque geometry and

8 102 Copyright c 2008 Tech Science Press CMES, vol.28, no.2, pp , 2008 mechanical conditions. The following assumptions were made when deriving the plaque growth functions: a) Plaque growth depends on local plaque morphological and mechanical conditions; b) For simplicity, a single-point correspondence approach is used, i.e., plaque growth at a given nodal point is determined by information from the past and current time steps at the same point. In reality, influence from neighboring points should already be factored in by the plaque growth at the chosen nodal point; c) The first order time derivatives of the vessel wall thickness and stress conditions (which are functions of time, for every nodal point under consideration) should be included in PGF; d) The information at most current time step is most directly related to plaque growth in the nearest future. Two different plaque growth function forms were used in our derivation and simulation process. Three-time-point function. Data from three time steps T 1, T 2 and T 3 are used to fit wall thickness or displacement at T 3. The fitting function is: f T3 fit(i)=a i f T2 (i)+b i df dt (i) T2 +c i σ T2 (i)+d i dσ dt (i) (23) T2 where f is the wall-thickness or displacement function, i is the nodal point numbering index (100 points were chosen from the inner boundary), σ is the Stress-P 1 function, a i, b i, c i,and d i are coefficients determined by the least squares method to fit T 3 wall thickness or displacement data, and the time derivatives are calculated by dg dt (i)= T2 w g T 3 (i) g T2 (i) +(1 w) g T 2 (i) g T1 (i) T 3 T 2 T 2 T 1 (24) where g can be either f or σ, w is a weight function to be adjusted for better agreement when the growth function is used in progression simulation. When the growth function (23) is used in the plaque progression simulation code to adjust the inner and outer boundary at each numerical time step, the time derivatives are evaluated with T3 replaced by the current time corresponding to the numerical step: dg dt (i)= tn w g t n (i) g T2 (i) +(1 w) g T 2 (i) g T1 (i) t n T 2 T 2 T 1 (25) where t n indicates the current time in the numerical simulation. Two-time-point function Data from two time steps T 1 and T 2 are used to fit wall thickness or displacement at T 2. The fitting function is the same as that given by (23) except that the derivatives are calculated using two time points by: dg dt (i)= g T 2 (i) g T1 (i). (26) T2 T 2 T 1 In the simulation code, (25) will be used to calculate time derivatives, the same way as it is done for the three-point growth function. So the main difference between the three-point formula and the two-point formula is that the coefficients a i, b i, c i, and d i are determined using different time (two or three) data points. The formulas for the two growth functions in the simulation code are almost the same. 3.3 Plaque progression simulation Using S4 at T 2 and the plaque growth function determined in 3.2, the following procedure is used to simulate plaque progression: Step 1. Start from the original in vivo MRI geometry (S4 at T 2 ), use proper shrinkage to get the zero-pressure numerical starting geometry; Step 2. Discretize the geometry using the meshless GFD method, solve the model to get

9 Finite Difference Method and Plaque Progression Simulation 103 plaque geometry and stress/strain distributions under specified pressure conditions; (a) S4 at T 2 under 136 mmhg (b) S4 at T 3 under 136 mmhg (c) Overlapping (a) and (b) Step 3. Use the growth function to determine the plaque geometry for next numerical time step by adjusting the nodal points on inner and outer boundaries; Step 4. Adjust internal nodal points as needed; (d) Simulated S4 at Day 220 (blue). (e) Simulated S4 at Day 400 (blue). (f) Simulated S4 at Day 651 (T3). Step 5. Solve the plaque model using the updated plaque geometry; Step 6. Repeat Steps 3-5 till numerical time reaches T 3 (MRI scan time). Results obtained from the simulation code are presented in next section. 4 Results Slice 4 was used as the sample to demonstrate the simulation process. The procedure described in 3.3 was followed. Fig. 8 gives the starting, target and three simulated geometries, showing good agreement between simulated and actual plaque progression. Define the absolute and relative errors as: Absolute Error = ( WT_num(i) WT_T 3(i) ), (27) 100 Figure 8: Simulated plaque growth has good agreement with MRI data (error = 4.52%). Progression code starting time: T 2 ; ending time: T 3. S4 was used as the sample slice. The three-point growth function was used. S4 geometry at T 3 determined from MRI data was used as the benchmark data for validation. Lumen pressure: 136 mmhg. T 2 T 3 time span was 651 days. the two-point growth function (26). The errors are: Absolute Error = cm, (31) Relative Error = 11.4%, (32) It is clear that the three-point growth function gave more accurate predictions. Relative Error = [( WT_num(i) WT_T 3(i) )/(WT_T 3(i) ], 100 (28) (a) Simulated S4 at Day 220 (blue). (b) Simulated S4 at Day 440 (blue). (c) Simulated S4 at Day 651 (T 3 ). where 100=total number of nodal points selected from inner boundary, we have: Absolute Error = cm, (29) Relative Error = 4.52%, (30) for the case simulated (S4, 3-point formula). Figure 9 gives the simulated plaque geometries at Days 220, 440, and 651 (ending time, T 3 )using Figure 9: Simulated plaque geometries at three numerical time steps showing the prediction accuracy (error = 11.4%) using the two-point growth function was not as good as that from the threepoint growth function. Simulations were also conducted for S2, S3 and

10 104 Copyright c 2008 Tech Science Press CMES, vol.28, no.2, pp , 2008 S5. Results from S3 were similar to that of S4. S2 and S5 did not give good results because the original MRI geometries differ too drastically from T 1 to T 2 and from T 2 to T 3 and no clear trend could be quantified. 5 Discussion 5.1 The computational approach It is well-known that atherosclerotic plaque progression is a multi-faceted process involving not only mechanical factors, but also plaque type, component size and location, cell activities, blood conditions such as cholesterol level, diabetes, changes caused by medication such as statin, and other chemical conditions, inflammation and lumen surface condition. Investigations and findings from all the channels, modalities and disciplines could be integrated together to obtain better and more thorough understanding of the complicated atherosclerotic progression process. In stead of trying to identify the individual factors which contribute to plaque growth, our computational simulation approach takes the actual multiyear plaque progression data to simulate future growth, based on the assumptions that the trend that led to the current state would continue into the future. In other words, the trend would continue; the governing mechanisms (whether we know them or not) would remain the same; the contributing factors would continue to contribute the same way as they did in the past, with the changes taken into consideration by the terms included in the growth function. Our results from the 4 slices simulated indicated that our model worked well when the trend from T 2 to T 3 was more or less similar to that from T 1 and T 2. And the model did not work very well when T 2 geometry seemed to be out of place between T 1 and T 3. Our model was applicable to 50% of the cases considered. 5.2 Model limitations Clearly the current model is very limited and serves as initial demonstrations of the MGFD method and the potential significant contributions from the progression simulation model. The model needs to be extended to 3D. And the growth functions need to be adjusted to include the proper terms that can represent the plaque growth trend. Better understanding of the biological and mechanical factors will help us to better formulate the growth function. Fluid forces and blood conditions (cholesterol, lipid lowering medication factor) can be included for better accuracy of predictions. 5.3 Progression simulation and plaque vulnerability assessment The long term goal of the current study is that patient-specific quantitative plaque growth functions can be determined based on multiple annual MRI scans and used to simulate plaque progression. Future plaque vulnerability can be assessed based on predicted future plaque morphologies by the progression simulation models. We are adding the time dimension into plaque assessment technology to improve the predicting power and accuracy. If our studies are successful, annual MRI scans would be recommended to patients who are in their early-to-middle stages of atherosclerosis. Simulated plaque growth would be generated for their early diagnosis and proper treatment for prevention of serious or even fatal clinical cardiovascular events. 6 Conclusion We believe that this is the first time that human carotid atherosclerotic plaque progression was simulated based on patient-specific plaque morphology and point-wise plaque growth functions derived from multi-year MRI data. Our results indicated that our proposed progression simulation process was able to accurately predict future plaque morphology if the current progression trend was continued. The meshless GFD method worked well for the progression model. The predicted progression by the three-time-point growth function was considerably more accurate than that given by the two-time-point (4.52% vs. 11.4%). The current 2D model can be extended to 3D model with more terms added to the growth functions for better predictions. More case studies are needed to validate our findings. Accurate plaque

11 Finite Difference Method and Plaque Progression Simulation 105 progression simulation adds the time dimension to plaque vulnerability assessment strategies and should improve our predicting accuracies. Acknowledgement: This research was supported in part by NIH grant NIH/NIBIB, R01 EB as part of the NSF/NIH Collaborative Research in Computational Neuroscience Program, and in part by NSF grant DMS Drs Vasily Yarnykh, Baocheng Chu, and Fei Liu contributed in the in vivo MRI data acquisition and image processing and their efforts are happily acknowledged. References Ahrem, R.; Beckert, A.; Wendland, H. (2006): A meshless spatial coupling scheme for largescale fluid-structure-interaction problems, CMES: Computer Modeling in Engineering & Sciences, 12(2): Atluri,S.N.(2005): Methods of Computer Modeling in Engineering & the Sciences-Part I, Tech Science Press, Forsyth, GA. Atluri,S.N.(2004): The Meshless Local-Petrov- Galerkin Method for Domain & BIE Discretizations, Tech Science Press, Forsyth, GA. Atluri, S. N.; Han, Z. D.; Rajendran, A. M. (2004): A new implementation of the meshless finite volume method, through the MLPG Mixed approach, CMES: Computer Modeling in Engineering & Sciences, 6 (6): Atluri, S. N.; Liu, H. T.; Han, Z. D. (2006a): Meshless local Petrov-Galerkin (MLPG) mixed collocation method for elasticity problems, CMES: Computer Modeling in Engineering & Sciences, 14(3): Atluri,S.N.;Liu,H.T.;Han,Z.D.(2006b): Meshless local Petrov-Galerkin (MLPG) mixed finite difference method for solid mechanics, CMES: Computer Modeling in Engineering & Sciences, 15(1): Atluri, S. N.; Yagawa, G.; Cruse, T. A. (Editors), (1995): Computational Mechanics 95, Vols. I & II, Springer-Verlag, New York. Bathe, K. J. (1996): Finite Element Procedures, Prentice Hall, Inc., New Jersey. Bathe, K. J. (Ed.) (2002): Theory and Modeling Guide, Vol I & II: ADINA and ADINA-F, ADINA R & D, Inc., Watertown, MA. Friedman, M. H.; Bargeron, C. B.; Deters, O. J.;Hutchins,G.M.;Mark,F.F.(1987): Correlation between wall shear and intimal thickness at a coronary artery branch, Atherosclerosis, 68: Friedman, M. H.; Giddens, D. P. (2005): Blood flow in major blood vessels - modeling and experiments, Annals of Biomedical Engineering, 33(12): Fung, Y. C. (1994): A First Course in Continuum Mechanics, Third Edition, Englewood Cliffs, New Jersey. Fung, Y. C. (1993): Biomechanics: Mechanical properties of Living Tissues, Springer-Verlag, New York, pp. 68. Fuster, V. (1998): The Vulnerable Atherosclerotic Plaque: Understanding, Identification, and Modification, Editor: Fuster, V. Co-Editors: Cornhill, J. F.; Dinsmore, R. E.; Fallon, J. T.; Insull, W.; Libby, P.; Nissen, S.; Rosenfeld, M. E.; Wagner, W.D. AHA Monograph series, Futura Publishing, Armonk NY. Fuster V.; Stein, B.; Ambrose, J. A.; Badimon, L.; Badimon, J. J.; Chesebro, J. H. (1990): Atherosclerotic plaque rupture and thrombosis, evolving concept. Circulation. 82 Supplement II:II-47 II-59, Giddens, D. P.; Zarins, C. K.; Glagov, S. (1993): The role of fluid mechanics in the localization and detection of atherosclerosis. Journal of Biomechanical Engineering, 115: Han, Z. D.; Atluri, S. N. (2004a): Meshless local Petrov-Galerkin (MLPG) approaches for solving 3D problems in elasto-statics, CMES: Computer Modeling in Engineering & Sciences, 6(2): Han, Z.D.; Atluri S. N. (2004b): A meshless local Petrov-Galerkin (MLPG) approach for 3- dimensional elasto-dynamics, CMC: Computers, Materials & Continua, 1(2): Han, Z. D.; Atluri, S. N. (2007): A System-

12 106 Copyright c 2008 Tech Science Press CMES, vol.28, no.2, pp , 2008 atic Approach for the Development of Weakly Singular BIEs, CMES: Computer Modeling in Engineering & Sciences, 21(1): Han, Z. D.; Liu, H. T.; Rajendran, A. H.; Atluri,S.N.(2006): The Applications of Meshless Local Petrov-Galerkin (MLPG) Approaches in High-Speed Impact, Penetration and Perforation Problems, CMES: Computer Modeling in Engineering & Sciences, 14(2): Han, Z. D.; Rajendran, A. M.; Atluri, S. N. (2005): Meshless local Petrov-Galerkin (MLPG) approaches for solving nonlinear problems with large deformations and rotations, CMES: Computer Modeling in Engineering & Sciences, 10(1):1-12. Hu, S. P.; Young, D. L.; Fan, C. M. (2008): FDMFS for Diffusion Equation with Unsteady Forcing Function, CMES: Computer Modeling in Engineering & Sciences, 24(1):1-20. Humphrey, J. D. (2002): Cardiovascular Solid Mechanics, Springer-Verlag, New York. Kerwin, W.; Hooker, A.; Spilker, M.; Vicini, P.; Ferguson, M.; Hatsukami, T.; and Yuan, C. (2003): Quantitative magnetic resonance imaging analysis of neovasculature volume in carotid atherosclerotic plaque. Circulation. 107(6), Kleiber, M. (1998): Handbook of Computational Solid Mechanics, Springer-Verlag, New York. Kobayashi, S.; Tsunoda, D.; Fukuzawa, Y.; Morikawa, H.; Tang, D.; Ku, D. N. (2003): Flow and compression in arterial models of stenosis with lipid core, Proceedings of 2003 ASME Summer Bioengineering Conference, Miami,FL, Ku, D. N.; Giddens, D.P.; Zarins, C.K.; Glagov, S. (1985): Pulsatile flow and atherosclerosis in the human carotid bifurcation: positive correlation between plaque location and low and oscillating shear stress. Arteriosclerosis, 5: Ling,X.W.;Atluri,S.N.(2006): Stability analysis for inverse heat conduction problems, CMES: Computer Modeling in Engineering & Sciences, 13(3): Liszka, T.; Orkisz, J. (1980): The finite difference method at arbitrary irregular grids and its application in applied mechanics, Computers and Structures, Vol. 11, Ma,Q.W.(2008): A New Meshless Interpolation Scheme for MLPG_R Method, CMES: Computer Modeling in Engineering & Sciences, 23(2): Mai-Duy, N.; Tran-Cong, T. (2004): Boundary integral-based domain decomposition technique for solution of Navier-Stokes equations, CMES: Computer Modeling in Engineering & Sciences, 6(1): Masuda, S; Noguchi, H. (2006): Analysis of Structure with Material Interface by Meshfree Method, CMES: Computer Modeling in Engineering & Sciences, 11(3): Naghavi, M.; Libby, P.; Falk, E.; Casscells, S. W.; Litovsky, S.; Rumberger, J.; Badimon, J. J.; Stefanadis, C.; Moreno, P.; Pasterkamp, G.; Fayad, Z.; Stone, P. H.; Waxman, S.; Raggi P.; Madjid, M.; Zarrabi, A.; Burke, A.; Yuan, C.; Fitzgerald, P. J.; Siscovick, D. S.; de Korte, C. L.; Aikawa, M.; Juhani, K. E.; Airaksinen, K. E.; Assmann, G.; Becker, C. R.; Chesebro, J. H.; Farb A.; Galis, Z. S.; Jackson, C.; Jang, I. K.; Koenig, W.; Lodder, R. A.; March, K.; Demirovic, J.; Navab, M.; Priori, S. G.; Rekhter, M. D.; Bahr, R.; Grundy, S. M.; Mehran, R.; Colombo, A.; Boerwinkle, E.; Ballantyne, C.; Insull, W. Jr.; Schwartz, R. S.; Vogel, R.; Serruys, P. W.; Hansson, G. K.; Faxon, D. P.; Kaul, S.; Drexler, H.; Greenland, P.; Muller, J. E.; Virmani, R.; Ridker, P. M.; Zipes, D. P.; Shah, P. K.; Willerson, J. T. (2003a): From vulnerable plaque to vulnerable patient: a call for new definitions and risk assessment strategies: Part I. Circulation. 108(14): Naghavi, M.; (co-authors are the same as above). (2003b): From vulnerable plaque to vulnerable patient: a call for new definitions and risk assessment strategies: Part II. Circulation, 108(15): Perko, J.; Sarler, B. (2007): Weight Function Shape Parameter Optimization in Meshless Methods for Non-uniform Grids, CMES: Computer

13 Finite Difference Method and Plaque Progression Simulation 107 Modeling in Engineering & Sciences, 19(1): Scotti, C. M.; Shkolnik, A. D.; Muluk, S. C.; Finol, E. A. (2005): Fluid-structure interaction in abdominal aortic aneurysms: effects of asymmetry and wall thickness. Biomed Eng Online, vol. 4: Shu, C.; Ding, H.; Yeo, K. S. (2005): Computation of incompressible Navier-Stokes equations by local RBF-based differential quadrature method, CMES: Computer Modeling in Engineering & Sciences, 7(2): Tang, D.; Chen, X. K.; Yang, C.; Kobayashi, K. and Ku, D. N. (2002): A viscoelastic model and meshless GFD method for Blood Flow in collapsible Stenotic Arteries, Advances in Computational Engineering & Sciences, Chap. 11, Tech Science Press, International Conference on Computational Engineering and Sciences. Tang, D.; Yang, C.; Kobayashi, K. and Ku, D. N. (2001): A Generalized Finite Difference Method for 3-D Viscous Flow in Stenotic Tubes with Large Wall Deformation and Collapse, Applied Num. Math., 38: pp Tang, D.; Yang, C.; Mondal, S; Liu, F.; Canton,G.;Hatsukami,T.S.andYuan,C.(2008): A Negative Correlation between Human Carotid Atherosclerotic Plaque Progression and Plaque Wall Stress: In Vivo MRI-Based 2D/3D FSI Models, J. Biomechanics, 41(4) : Tang, D.; Yang, C.; Zheng, J.; Woodard, P. K.; Saffitz, J. E.; Petruccelli, J. D.; Sicard, G.A.; Yuan, C. (2005), Local maximal stress hypothesis and computational plaque vulnerability index for atherosclerotic plaque assessment, Annals of Biomedical Engineering, 33(12): Wacher, A.; Givoli, D. (2006): Remeshing and Refining with Moving Finite Elements. Application to Nonlinear Wave Problems, CMES: Computer Modeling in Engineering & Sciences, 15(3): Wu, X. H.; Shen, S. P.; Tao, W. Q. (2007): Meshless Local Petrov-Galerkin Collocation Method for Two-dimensional Heat Conduction Problems, CMES: Computer Modeling in Engineering & Sciences, 22(1): Yang, C.; Tang, D.; Yuan, C.; Hatsukami, T. S.; Zheng, J. and Woodard, P. K. (2007): In Vivo/Ex Vivo MRI-Based 3D Models with Fluid- Structure Interactions for Human Atherosclerotic Plaques Compared with Fluid/Wall-Only Models, CMES: Computer Modeling in Engineering and Sciences, 19(3): Yuan, C.; Kerwin, W. S. (2004): MRI of atherosclerosis. Journal of Magnetic Resonance Imaging. 19(6): Yuan, C.; Mitsumori, L. M.; Beach, K. W.; Maravilla, K. R. (2001): Special review: carotid atherosclerotic plaque: noninvasive MR characterization and identification of vulnerable lesions, Radiology, 221, Zohouri S.; Pirooz, M. D.; and Esmaeily, A. (2004): Predicting wave run-up using full ALE finite element approach considering moving boundary, CMES: Computer Modeling in Engineering & Sciences, 7(1):

14

Thesis for the Master degree in Engineering Research on Negative Pressure Wave Simulation and Signal Processing of Fluid-Conveying Pipeline Leak Candi

Thesis for the Master degree in Engineering Research on Negative Pressure Wave Simulation and Signal Processing of Fluid-Conveying Pipeline Leak Candi U17 10220 UDC624 Thesis for the Master degree in Engineering Research on Negative Pressure Wave Simulation and Signal Processing of Fluid-Conveying Pipeline Leak Candidate:Chen Hao Tutor: Xue Jinghong

More information

國家圖書館典藏電子全文

國家圖書館典藏電子全文 i ii Abstract The most important task in human resource management is to encourage and help employees to develop their potential so that they can fully contribute to the organization s goals. The main

More information

UDC Empirical Researches on Pricing of Corporate Bonds with Macro Factors 厦门大学博硕士论文摘要库

UDC Empirical Researches on Pricing of Corporate Bonds with Macro Factors 厦门大学博硕士论文摘要库 10384 15620071151397 UDC Empirical Researches on Pricing of Corporate Bonds with Macro Factors 2010 4 Duffee 1999 AAA Vasicek RMSE RMSE Abstract In order to investigate whether adding macro factors

More information

Fig. 1 1 The sketch for forced lead shear damper mm 45 mm 4 mm 200 mm 25 mm 2 mm mm Table 2 The energy dissip

Fig. 1 1 The sketch for forced lead shear damper mm 45 mm 4 mm 200 mm 25 mm 2 mm mm Table 2 The energy dissip * - 1 1 2 3 1. 100124 2. 100124 3. 210018 - ABAQUS - DOI 10. 13204 /j. gyjz201511033 EXPERIMENTAL STUDY AND THEORETICAL MODEL OF A NEW TYPE OF STEEL-LEAD DAMPING Shen Fei 1 Xue Suduo 1 Peng Lingyun 2 Ye

More information

Corporate Social Responsibility CSR CSR CSR 1 2 ~ CSR 6 CSR 7 CSR 8 CSR 9 10 ~ CSR 14 CSR CSR 2013 A A 23.

Corporate Social Responsibility CSR CSR CSR 1 2 ~ CSR 6 CSR 7 CSR 8 CSR 9 10 ~ CSR 14 CSR CSR 2013 A A 23. 24 3 Vol. 24 No. 3 2015 6 OPERATIONS RESEARCH AND MANAGEMENT SCIENCE Jun. 2015 1 2 2 1. 300071 2. 300071 Markowitz 10 F830. 59 A 1007-3221 2015 03-0275-13 Improvement of Portfolio Models Research An Empirical

More information

Time Estimation of Occurrence of Diabetes-Related Cardiovascular Complications by Ching-Yuan Hu A thesis submitted in partial fulfillment of the requi

Time Estimation of Occurrence of Diabetes-Related Cardiovascular Complications by Ching-Yuan Hu A thesis submitted in partial fulfillment of the requi Time Estimation of Occurrence of Diabetes-Related Cardiovascular Complications by Ching-Yuan Hu Master of Science 2011 Institute of Chinese Medical Sciences University of Macau Time Estimation of Occurrence

More information

東吳大學

東吳大學 律 律 論 論 療 行 The Study on Medical Practice and Coercion 林 年 律 律 論 論 療 行 The Study on Medical Practice and Coercion 林 年 i 讀 臨 療 留 館 讀 臨 律 六 礪 讀 不 冷 療 臨 年 裡 歷 練 禮 更 老 林 了 更 臨 不 吝 麗 老 劉 老 論 諸 見 了 年 金 歷 了 年

More information

Microsoft PowerPoint _代工實例-1

Microsoft PowerPoint _代工實例-1 4302 動態光散射儀 (Dynamic Light Scattering) 代工實例與結果解析 生醫暨非破壞性分析團隊 2016.10 updated Which Size to Measure? Diameter Many techniques make the useful and convenient assumption that every particle is a sphere. The

More information

( s y s t e m ) ( s t r e s s ) (stress model) ( s y s t e m ) [ ] [ 5 ] C o x [ 3 ] 1 [ 1, 2 ] [ 6-8 ] [ 9 ] Tw Fam Med Res 2003 Vol.1 No.1 23

( s y s t e m ) ( s t r e s s ) (stress model) ( s y s t e m ) [ ] [ 5 ] C o x [ 3 ] 1 [ 1, 2 ] [ 6-8 ] [ 9 ] Tw Fam Med Res 2003 Vol.1 No.1 23 1, 4 1, 2 3 1 4 1, 2 1977 E n g e l 1999 5 1 ( 1 ) ( 2 ) ( 3 ) ( 4 ) ( 5 ) 2003 1 2 2-32 stress model, stress analysis [ 1 ] ( p s y c h o g e n e s i s ) (holism) (biomedical model) ( r e d u c t i o

More information

θ 1 = φ n -n 2 2 n AR n φ i = 0 1 = a t - θ θ m a t-m 3 3 m MA m 1. 2 ρ k = R k /R 0 5 Akaike ρ k 1 AIC = n ln δ 2

θ 1 = φ n -n 2 2 n AR n φ i = 0 1 = a t - θ θ m a t-m 3 3 m MA m 1. 2 ρ k = R k /R 0 5 Akaike ρ k 1 AIC = n ln δ 2 35 2 2012 2 GEOMATICS & SPATIAL INFORMATION TECHNOLOGY Vol. 35 No. 2 Feb. 2012 1 2 3 4 1. 450008 2. 450005 3. 450008 4. 572000 20 J 101 20 ARMA TU196 B 1672-5867 2012 02-0213 - 04 Application of Time Series

More information

2 ( 自 然 科 学 版 ) 第 20 卷 波 ). 这 种 压 缩 波 空 气 必 然 有 一 部 分 要 绕 流 到 车 身 两 端 的 环 状 空 间 中, 形 成 与 列 车 运 行 方 向 相 反 的 空 气 流 动. 在 列 车 尾 部, 会 产 生 低 于 大 气 压 的 空 气 流

2 ( 自 然 科 学 版 ) 第 20 卷 波 ). 这 种 压 缩 波 空 气 必 然 有 一 部 分 要 绕 流 到 车 身 两 端 的 环 状 空 间 中, 形 成 与 列 车 运 行 方 向 相 反 的 空 气 流 动. 在 列 车 尾 部, 会 产 生 低 于 大 气 压 的 空 气 流 第 20 卷 第 3 期 2014 年 6 月 ( 自 然 科 学 版 ) JOURNAL OF SHANGHAI UNIVERSITY (NATURAL SCIENCE) Vol. 20 No. 3 June 2014 DOI: 10.3969/j.issn.1007-2861.2013.07.031 基 于 FLUENT 测 轨 道 交 通 隧 道 中 电 波 折 射 率 结 构 常 数 张 永

More information

穨423.PDF

穨423.PDF Chinese Journal of Science Education 2002,, 423-439 2002, 10(4), 423-439 1 2 1 1 1 2 90 8 10 91 4 9 91 8 22 ) NII 1995 7 14, 1999 1997 (Cooperative Remotely Accessible Learning CORAL) 424 (Collaborative

More information

2 199 Navier-Stokes { u t - v 2 u + u u + 1 p = f 1 ρ u = 0 ux y Case25CTA z= { u LC6 60% T n = 0 2 CAS Case25pre ux y z= 0 Case25post u = u x u

2 199 Navier-Stokes { u t - v 2 u + u u + 1 p = f 1 ρ u = 0 ux y Case25CTA z= { u LC6 60% T n = 0 2 CAS Case25pre ux y z= 0 Case25post u = u x u 42 2 Vol 42 No 2 2018 3 Journal of Jiangxi Normal UniversityNatural Science Mar 2018 1000-5862201802-0199-05 1 2 2 1* 2 2 1 3300222 518055 1 2 12 10 6 CPU 128 72 43% O 246 A DOI10 16357 /j cnki issn1000-5862

More information

(baking powder) 1 ( ) ( ) 1 10g g (two level design, D-optimal) 32 1/2 fraction Two Level Fractional Factorial Design D-Optimal D

(baking powder) 1 ( ) ( ) 1 10g g (two level design, D-optimal) 32 1/2 fraction Two Level Fractional Factorial Design D-Optimal D ( ) 4 1 1 1 145 1 110 1 (baking powder) 1 ( ) ( ) 1 10g 1 1 2.5g 1 1 1 1 60 10 (two level design, D-optimal) 32 1/2 fraction Two Level Fractional Factorial Design D-Optimal Design 1. 60 120 2. 3. 40 10

More information

1 引言

1 引言 P P 第 40 卷 Vol.40 第 7 期 No.7 计 算 机 工 程 Computer Engineering 014 年 7 月 July 014 开 发 研 究 与 工 程 应 用 文 章 编 号 :1000-348(014)07-081-05 文 献 标 识 码 :A 中 图 分 类 号 :TP391.41 摘 基 于 图 像 识 别 的 震 象 云 地 震 预 测 方 法 谢 庭,

More information

Microsoft Word - 刘 慧 板.doc

Microsoft Word - 刘  慧 板.doc 中 国 环 境 科 学 2012,32(5):933~941 China Environmental Science 系 统 动 力 学 在 空 港 区 域 规 划 环 境 影 响 评 价 中 的 应 用 刘 慧 1,2, 郭 怀 成 1*, 盛 虎 1, 都 小 尚 1,3, 李 娜 1 1, 杨 永 辉 (1. 北 京 大 学 环 境 科 学 与 工 程 学 院, 北 京 100871; 2.

More information

Improved Preimage Attacks on AES-like Hash Functions: Applications to Whirlpool and Grøstl

Improved Preimage Attacks on AES-like Hash Functions: Applications to Whirlpool and Grøstl SKLOIS (Pseudo) Preimage Attack on Reduced-Round Grøstl Hash Function and Others Shuang Wu, Dengguo Feng, Wenling Wu, Jian Guo, Le Dong, Jian Zou March 20, 2012 Institute. of Software, Chinese Academy

More information

#4 ~ #5 12 m m m 1. 5 m # m mm m Z4 Z5

#4 ~ #5 12 m m m 1. 5 m # m mm m Z4 Z5 2011 6 6 153 JOURNAL OF RAILWAY ENGINEERING SOCIETY Jun 2011 NO. 6 Ser. 153 1006-2106 2011 06-0014 - 07 300142 ABAQUS 4. 287 mm 6. 651 mm U455. 43 A Analysis of Impact of Shield Tunneling on Displacement

More information

10384 19020101152519 UDC Rayleigh Quasi-Rayleigh Method for computing eigenvalues of symmetric tensors 2 0 1 3 2 0 1 3 2 0 1 3 2013 , 1. 2. [4], [27].,. [6] E- ; [7], Z-. [15]. Ramara G. kolda [1, 2],

More information

Microsoft Word - 01李惠玲ok.doc

Microsoft Word - 01李惠玲ok.doc 康 寧 學 報 11:1-20(2009) 1 數 位 學 習 於 護 理 技 術 課 程 之 運 用 與 評 值 * 李 惠 玲 ** 高 清 華 *** 呂 莉 婷 摘 要 背 景 : 網 路 科 技 在 教 育 的 使 用 已 成 為 一 種 有 利 的 教 學 輔 助 工 具 網 路 教 學 的 特 性, 在 使 學 習 可 不 分 時 間 與 空 間 不 同 進 度 把 握 即 時 性 資

More information

Microsoft Word - TIP006SCH Uni-edit Writing Tip - Presentperfecttenseandpasttenseinyourintroduction readytopublish

Microsoft Word - TIP006SCH Uni-edit Writing Tip - Presentperfecttenseandpasttenseinyourintroduction readytopublish 我 难 度 : 高 级 对 们 现 不 在 知 仍 道 有 听 影 过 响 多 少 那 次 么 : 研 英 究 过 文 论 去 写 文 时 作 的 表 技 引 示 巧 言 事 : 部 情 引 分 发 言 该 生 使 在 中 用 过 去, 而 现 在 完 成 时 仅 表 示 事 情 发 生 在 过 去, 并 的 哪 现 种 在 时 完 态 成 呢 时? 和 难 过 道 去 不 时 相 关? 是 所 有

More information

THE APPLICATION OF ISOTOPE RATIO ANALYSIS BY INDUCTIVELY COUPLED PLASMA MASS SPECTROMETER A Dissertation Presented By Chaoyong YANG Supervisor: Prof.D

THE APPLICATION OF ISOTOPE RATIO ANALYSIS BY INDUCTIVELY COUPLED PLASMA MASS SPECTROMETER A Dissertation Presented By Chaoyong YANG Supervisor: Prof.D 10384 070302 9825042 UDC 2001.6. 2001.7. 20016 THE APPLICATION OF ISOTOPE RATIO ANALYSIS BY INDUCTIVELY COUPLED PLASMA MASS SPECTROMETER A Dissertation Presented By Chaoyong YANG Supervisor: Prof.Dr. Xiaoru

More information

度 身 體 活 動 量 ; 芬 蘭 幼 兒 呈 現 中 度 身 體 活 動 量 之 比 例 高 於 臺 灣 幼 兒 (5) 幼 兒 在 投 入 度 方 面 亦 達 顯 著 差 異 (χ²=185.35, p <.001), 芬 蘭 與 臺 灣 幼 兒 多 半 表 現 出 中 度 投 入 與 高 度

度 身 體 活 動 量 ; 芬 蘭 幼 兒 呈 現 中 度 身 體 活 動 量 之 比 例 高 於 臺 灣 幼 兒 (5) 幼 兒 在 投 入 度 方 面 亦 達 顯 著 差 異 (χ²=185.35, p <.001), 芬 蘭 與 臺 灣 幼 兒 多 半 表 現 出 中 度 投 入 與 高 度 臺 灣 與 芬 蘭 幼 兒 園 室 內 自 由 遊 戲 內 涵 之 探 討 林 昭 溶 毛 萬 儀 經 國 管 理 暨 健 康 學 院 幼 兒 保 育 系 副 教 授 joyce@ems.cku.edu.tw 吳 敏 而 國 家 教 育 研 究 院 研 究 員 rozwu@mail.naer.edu.tw wanyi@ems.cku.edu.tw 摘 要 自 由 遊 戲 被 視 為 是 幼 兒 的

More information

2008 Nankai Business Review 61

2008 Nankai Business Review 61 150 5 * 71272026 60 2008 Nankai Business Review 61 / 62 Nankai Business Review 63 64 Nankai Business Review 65 66 Nankai Business Review 67 68 Nankai Business Review 69 Mechanism of Luxury Brands Formation

More information

untitled

untitled Co-integration and VECM Yi-Nung Yang CYCU, Taiwan May, 2012 不 列 1 Learning objectives Integrated variables Co-integration Vector Error correction model (VECM) Engle-Granger 2-step co-integration test Johansen

More information

~ 10 2 P Y i t = my i t W Y i t 1000 PY i t Y t i W Y i t t i m Y i t t i 15 ~ 49 1 Y Y Y 15 ~ j j t j t = j P i t i = 15 P n i t n Y

~ 10 2 P Y i t = my i t W Y i t 1000 PY i t Y t i W Y i t t i m Y i t t i 15 ~ 49 1 Y Y Y 15 ~ j j t j t = j P i t i = 15 P n i t n Y * 35 4 2011 7 Vol. 35 No. 4 July 2011 3 Population Research 1950 ~ 1981 The Estimation Method and Its Application of Cohort Age - specific Fertility Rates Wang Gongzhou Hu Yaoling Abstract Based on the

More information

ENGG1410-F Tutorial 6

ENGG1410-F Tutorial 6 Jianwen Zhao Department of Computer Science and Engineering The Chinese University of Hong Kong 1/16 Problem 1. Matrix Diagonalization Diagonalize the following matrix: A = [ ] 1 2 4 3 2/16 Solution The

More information

: ( ),,

: ( ),, Case Discussion: ( ), RN, ScD, MPH : 24,,,,,,, ( ) ) ( ), : ( ),, randomized controlled trial (RCT). minimization program,,,, : : ; Apgar score > 7 ;, : ; ; ; ; (BT180/min; 2, 83,50 (60.2%).,

More information

报 告 1: 郑 斌 教 授, 美 国 俄 克 拉 荷 马 大 学 医 学 图 像 特 征 分 析 与 癌 症 风 险 评 估 方 法 摘 要 : 准 确 的 评 估 癌 症 近 期 发 病 风 险 和 预 后 或 者 治 疗 效 果 是 发 展 和 建 立 精 准 医 学 的 一 个 重 要 前

报 告 1: 郑 斌 教 授, 美 国 俄 克 拉 荷 马 大 学 医 学 图 像 特 征 分 析 与 癌 症 风 险 评 估 方 法 摘 要 : 准 确 的 评 估 癌 症 近 期 发 病 风 险 和 预 后 或 者 治 疗 效 果 是 发 展 和 建 立 精 准 医 学 的 一 个 重 要 前 东 北 大 学 中 荷 生 物 医 学 与 信 息 工 程 学 院 2016 年 度 生 物 医 学 与 信 息 工 程 论 坛 会 议 时 间 2016 年 6 月 8 日, 星 期 三,9:30 至 16:00 会 议 地 址 会 议 网 址 主 办 单 位 东 北 大 学 浑 南 校 区 沈 阳 市 浑 南 区 创 新 路 195 号 生 命 科 学 大 楼 B 座 619 报 告 厅 http://www.bmie.neu.edu.cn

More information

WTO

WTO 10384 X0115018 UDC MBA 2004 5 14 2004 6 1 WTO 2004 2006 7 2 Abstract According to the promise after our country enter into WTO, our country will open the readymade oil retail market in the end of 2004

More information

國立中山大學學位論文典藏.PDF

國立中山大學學位論文典藏.PDF I II III The Study of Factors to the Failure or Success of Applying to Holding International Sport Games Abstract For years, holding international sport games has been Taiwan s goal and we are on the way

More information

地質調査研究報告/Bulletin of the Geological Survey of Japan

地質調査研究報告/Bulletin of the Geological Survey of Japan Shigeru Suto, Takayuki Inomata, Hisashi Sasaki and Sakae Mukoyama (2007) Data base of the volcanic ash fall distribution map of Japan. Bull. Geol. Surv. Japan, vol. 58(9/10), p.261-321, 8 figs, 2 tables,

More information

PowerPoint Presentation

PowerPoint Presentation Decision analysis 量化決策分析方法專論 2011/5/26 1 Problem formulation- states of nature In the decision analysis, decision alternatives are referred to as chance events. The possible outcomes for a chance event

More information

20

20 37 92 19 40 19 20 21 1 7 22 1/5 6/30 5/3030 23 24 25 26 1 2 27 1 2 28 29 30 5 8 8 3 31 32 33 34 35 36 37 38 39 A Study Investigating Elementary School Students Concept of the Unit in Fraction in Northern

More information

Gassama Abdoul Gadiri University of Science and Technology of China A dissertation for master degree Ordinal Probit Regression Model and Application in Credit Rating for Users of Credit Card Author :

More information

國立高雄大學○○○○○○學系(研究所)(標楷體18號字

國立高雄大學○○○○○○學系(研究所)(標楷體18號字 國 立 高 雄 大 學 都 市 發 展 與 建 築 研 究 所 碩 士 論 文 高 雄 後 勁 地 區 傳 統 民 居 特 徵 之 研 究 The Study of The Characteristic of Traditional Residential Buildings of Houjing District in Kaohsiung 研 究 生 : 許 輝 隆 撰 指 導 教 授 : 陳 啟

More information

天 主 教 輔 仁 大 學 社 會 學 系 學 士 論 文 小 別 勝 新 婚? 久 別 要 離 婚? 影 響 遠 距 家 庭 婚 姻 感 情 因 素 之 探 討 Separate marital relations are getting better or getting worse? -Exp

天 主 教 輔 仁 大 學 社 會 學 系 學 士 論 文 小 別 勝 新 婚? 久 別 要 離 婚? 影 響 遠 距 家 庭 婚 姻 感 情 因 素 之 探 討 Separate marital relations are getting better or getting worse? -Exp 天 主 教 輔 仁 大 學 社 會 學 系 學 士 論 文 小 別 勝 新 婚? 久 別 要 離 婚? 影 響 遠 距 家 庭 婚 姻 感 情 因 素 之 探 討 Separate marital relations are getting better or getting worse? -Explore the impact of emotional factors couples do not

More information

從詩歌的鑒賞談生命價值的建構

從詩歌的鑒賞談生命價值的建構 Viktor E. Frankl (logotherapy) (will-to-meaning) (creative values) Ture (Good) (Beauty) (experiential values) (attitudinal values) 1 2 (logotherapy) (biological) (2) (psychological) (3) (noölogical) (4)

More information

Microsoft Word - 第四組心得.doc

Microsoft Word - 第四組心得.doc 徐 婉 真 這 四 天 的 綠 島 人 權 體 驗 營 令 我 印 象 深 刻, 尤 其 第 三 天 晚 上 吳 豪 人 教 授 的 那 堂 課, 他 讓 我 聽 到 不 同 於 以 往 的 正 義 之 聲 轉 型 正 義, 透 過 他 幽 默 熱 情 的 語 調 激 起 了 我 對 政 治 的 興 趣, 願 意 在 未 來 多 關 心 社 會 多 了 解 政 治 第 一 天 抵 達 綠 島 不 久,

More information

: 307, [], [2],,,, [3] (Response Surface Methodology, RSA),,, [4,5] Design-Expert 6.0,,,, [6] VPJ33 ph 3,, ph, OD, Design-Expert 6.0 Box-Behnken, VPJ3

: 307, [], [2],,,, [3] (Response Surface Methodology, RSA),,, [4,5] Design-Expert 6.0,,,, [6] VPJ33 ph 3,, ph, OD, Design-Expert 6.0 Box-Behnken, VPJ3 微 生 物 学 通 报 FEB 20, 2008, 35(2) 306~30 Microbiology 2008 by Institute of Microbiology, CAS tongbao@im.ac.cn 生 物 实 验 室 响 应 面 分 析 法 优 化 副 溶 血 性 弧 菌 生 长 条 件 刘 代 新 宁 喜 斌 * 张 继 伦 2 (. 200090) (2. 20202) 摘 要

More information

:1949, 1936, 1713 %, 63 % (, 1957, 5 ), :?,,,,,, (,1999, 329 ),,,,,,,,,, ( ) ; ( ), 1945,,,,,,,,, 100, 1952,,,,,, ,, :,,, 1928,,,,, (,1984, 109

:1949, 1936, 1713 %, 63 % (, 1957, 5 ), :?,,,,,, (,1999, 329 ),,,,,,,,,, ( ) ; ( ), 1945,,,,,,,,, 100, 1952,,,,,, ,, :,,, 1928,,,,, (,1984, 109 2006 9 1949 3 : 1949 2005, : 1949 1978, ; 1979 1997, ; 1998 2005,,, :,,, 1949, :, ;,,,, 50, 1952 1957 ; ; 60 ; 1978 ; 2003,,,,,,, 1953 1978 1953 1978,,,, 100,,,,, 3,, :100836, :wulijjs @263. net ;,, :

More information

学 校 编 码 :10384 分 类 号 密 级 学 号 :X2007155130 UDC 厦 门 怡 福 养 生 健 康 管 理 有 限 公 司 创 业 计 划 王 韬 指 导 教 师 姓 名 : 郭 霖 教 授 厦 门 大 学 硕 士 学 位 论 文 厦 门 怡 福 养 生 健 康 管 理 有 限 公 司 创 业 计 划 A Business Plan for Xiamen Eve Health

More information

07 353-408

07 353-408 編 號 :CCMP96-RD-006 以 德 菲 法 探 討 中 醫 臨 床 證 型 診 斷 基 準 - 六 淫 證 型 蘇 奕 彰 中 國 醫 藥 大 學 摘 要 德 菲 法 (Delphi Method) 是 廣 泛 被 使 用 的 研 究 方 法, 尤 其 適 用 於 新 的 研 究 領 域 (new research area) 及 探 索 性 研 究 (exploratory studies),

More information

附件1:

附件1: 附 件 1: 全 国 优 秀 教 育 硕 士 专 业 学 位 论 文 推 荐 表 单 位 名 称 : 西 南 大 学 论 文 题 目 填 表 日 期 :2014 年 4 月 30 日 数 学 小 组 合 作 学 习 的 课 堂 管 理 攻 硕 期 间 及 获 得 硕 士 学 位 后 一 年 内 获 得 与 硕 士 学 位 论 文 有 关 的 成 果 作 者 姓 名 论 文 答 辩 日 期 学 科 专

More information

1556 地 理 科 学 进 展 30 卷 他 关 于 农 村 住 房 结 构 与 抗 震 性 能 的 研 究, 则 多 是 从 工 程 抗 灾 的 角 度, 研 究 某 种 构 造 类 型 的 房 屋, 力 图 找 到 传 统 房 屋 的 结 构 失 误 和 新 建 房 屋 中 存 在 的 问 [

1556 地 理 科 学 进 展 30 卷 他 关 于 农 村 住 房 结 构 与 抗 震 性 能 的 研 究, 则 多 是 从 工 程 抗 灾 的 角 度, 研 究 某 种 构 造 类 型 的 房 屋, 力 图 找 到 传 统 房 屋 的 结 构 失 误 和 新 建 房 屋 中 存 在 的 问 [ 第 30 卷 第 12 期 2011 年 12 月 地 理 科 学 进 展 PROGRESS IN GEOGRAPHY Vol.30, No.12 Dec., 2012 中 国 农 村 地 区 住 房 结 构 的 区 域 差 异 及 其 影 响 因 素 1, 2, 李 方 一 3 1,, 高 晓 路 2 2, 王 英 杰 (1. 中 国 科 学 院 区 域 可 持 续 发 展 分 析 与 模 拟 重

More information

Abstract Today, the structures of domestic bus industry have been changed greatly. Many manufacturers enter into the field because of its lower thresh

Abstract Today, the structures of domestic bus industry have been changed greatly. Many manufacturers enter into the field because of its lower thresh SWOT 5 Abstract Today, the structures of domestic bus industry have been changed greatly. Many manufacturers enter into the field because of its lower threshold. All of these lead to aggravate drastically

More information

医学科研方法

医学科研方法 医 学 科 研 方 法 主 讲 苏 秀 兰 细 胞 生 物 学 教 授 第 一 讲 绪 论 一 概 念 : 在 医 学 科 学 研 究 活 动 中, 运 用 科 学 的 实 践 和 理 论 思 维 的 技 巧, 充 分 发 挥 自 己 的 智 慧, 去 寻 找 观 念 世 界 与 现 实 世 界 之 间 的 联 系 或 规 律, 从 而 发 现 新 现 象, 提 出 新 理 论 的 手 段 二 医

More information

JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION Vol. 31 No. 5 Oct /35 TU3521 P315.

JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION Vol. 31 No. 5 Oct /35 TU3521 P315. 31 5 2011 10 JOURNAL OF EARTHQUAKE ENGINEERING AND ENGINEERING VIBRATION Vol. 31 No. 5 Oct. 2011 1000-1301 2011 05-0075 - 09 510405 1 /35 TU3521 P315. 8 A Earthquake simulation shaking table test and analysis

More information

A Study on Grading and Sequencing of Senses of Grade-A Polysemous Adjectives in A Syllabus of Graded Vocabulary for Chinese Proficiency 2002 I II Abstract ublished in 1992, A Syllabus of Graded Vocabulary

More information

致 谢 本 人 自 2008 年 6 月 从 上 海 外 国 语 大 学 毕 业 之 后, 于 2010 年 3 月 再 次 进 入 上 外, 非 常 有 幸 成 为 汉 语 国 际 教 育 专 业 的 研 究 生 回 顾 三 年 以 来 的 学 习 和 生 活, 顿 时 感 觉 这 段 时 间 也

致 谢 本 人 自 2008 年 6 月 从 上 海 外 国 语 大 学 毕 业 之 后, 于 2010 年 3 月 再 次 进 入 上 外, 非 常 有 幸 成 为 汉 语 国 际 教 育 专 业 的 研 究 生 回 顾 三 年 以 来 的 学 习 和 生 活, 顿 时 感 觉 这 段 时 间 也 精 英 汉 语 和 新 实 用 汉 语 课 本 的 对 比 研 究 The Comparative Study of Jing Ying Chinese and The New Practical Chinese Textbook 专 业 : 届 别 : 姓 名 : 导 师 : 汉 语 国 际 教 育 2013 届 王 泉 玲 杨 金 华 1 致 谢 本 人 自 2008 年 6 月 从 上 海 外

More information

A VALIDATION STUDY OF THE ACHIEVEMENT TEST OF TEACHING CHINESE AS THE SECOND LANGUAGE by Chen Wei A Thesis Submitted to the Graduate School and Colleg

A VALIDATION STUDY OF THE ACHIEVEMENT TEST OF TEACHING CHINESE AS THE SECOND LANGUAGE by Chen Wei A Thesis Submitted to the Graduate School and Colleg 上 海 外 国 语 大 学 SHANGHAI INTERNATIONAL STUDIES UNIVERSITY 硕 士 学 位 论 文 MASTER DISSERTATION 学 院 国 际 文 化 交 流 学 院 专 业 汉 语 国 际 教 育 硕 士 题 目 届 别 2010 届 学 生 陈 炜 导 师 张 艳 莉 副 教 授 日 期 2010 年 4 月 A VALIDATION STUDY

More information

University of Science and Technology of China A dissertation for master s degree Research of e-learning style for public servants under the context of

University of Science and Technology of China A dissertation for master s degree Research of e-learning style for public servants under the context of 中 国 科 学 技 术 大 学 硕 士 学 位 论 文 新 媒 体 环 境 下 公 务 员 在 线 培 训 模 式 研 究 作 者 姓 名 : 学 科 专 业 : 导 师 姓 名 : 完 成 时 间 : 潘 琳 数 字 媒 体 周 荣 庭 教 授 二 一 二 年 五 月 University of Science and Technology of China A dissertation for

More information

The Development of Color Constancy and Calibration System

The Development of Color Constancy and Calibration System The Development of Color Constancy and Calibration System The Development of Color Constancy and Calibration System LabVIEW CCD BMP ii Abstract The modern technologies develop more and more faster, and

More information

BIBLID 0254-4466(2001)19:1 pp. 249-276 19 1 90 6 ** * ** 88 I 2000 8 249 250 19 1 251 1873-1929 1900 1 1902 1 35 1900 1960 7-12 252 19 1 2 3 2 1900 1902 3 2000 129-197 253 4 5 6 4 1902 1962 103 5 Joseph

More information

硕 士 学 位 论 文 论 文 题 目 : 北 岛 诗 歌 创 作 的 双 重 困 境 专 业 名 称 : 中 国 现 当 代 文 学 研 究 方 向 : 中 国 新 诗 研 究 论 文 作 者 : 奚 荣 荣 指 导 老 师 : 姜 玉 琴 2014 年 12 月

硕 士 学 位 论 文 论 文 题 目 : 北 岛 诗 歌 创 作 的 双 重 困 境 专 业 名 称 : 中 国 现 当 代 文 学 研 究 方 向 : 中 国 新 诗 研 究 论 文 作 者 : 奚 荣 荣 指 导 老 师 : 姜 玉 琴 2014 年 12 月 硕 士 学 位 论 文 论 文 题 目 : 北 岛 诗 歌 创 作 的 双 重 困 境 专 业 名 称 : 中 国 现 当 代 文 学 研 究 方 向 : 中 国 新 诗 研 究 论 文 作 者 : 奚 荣 荣 指 导 老 师 : 姜 玉 琴 2014 年 12 月 致 谢 文 学 是 我 们 人 类 宝 贵 的 精 神 财 富 两 年 半 的 硕 士 学 习 让 我 进 一 步 接 近 文 学,

More information

untitled

untitled 20 1 2010 10 Vol.20 Special 1 The Chinese Journal of Nonferrous Metals Oct. 2010 1004-0609(2010)S1-s0127-05 Ti-6Al-4V 1 2 2 (1. 710016 2., 710049) 500~1 000 20 Ti-6Al-4V(TC4) TC4 800 TC4 800 TC4 TC4 800

More information

BC04 Module_antenna__ doc

BC04 Module_antenna__ doc http://www.infobluetooth.com TEL:+86-23-68798999 Fax: +86-23-68889515 Page 1 of 10 http://www.infobluetooth.com TEL:+86-23-68798999 Fax: +86-23-68889515 Page 2 of 10 http://www.infobluetooth.com TEL:+86-23-68798999

More information

中國傳統醫學及養生保健學說中,與經脈及穴道有密不可分的關係

中國傳統醫學及養生保健學說中,與經脈及穴道有密不可分的關係 Journal of National Kaohsiung University of Applied Sciences The Study of Correlation between Time and Human Meridian Point Bioenergy Chien-Min Cheng Ching-Yi Sung 577 * ** 578 The study of correlation

More information

Chinese Journal of Applied Probability and Statistics Vol.25 No.4 Aug (,, ;,, ) (,, ) 应用概率统计 版权所有, Zhang (2002). λ q(t)

Chinese Journal of Applied Probability and Statistics Vol.25 No.4 Aug (,, ;,, ) (,, ) 应用概率统计 版权所有, Zhang (2002). λ q(t) 2009 8 Chinese Journal of Applied Probability and Statistics Vol.25 No.4 Aug. 2009,, 541004;,, 100124),, 100190), Zhang 2002). λ qt), Kolmogorov-Smirov, Berk and Jones 1979). λ qt).,,, λ qt),. λ qt) 1,.

More information

既 設 建 築 物 汙 水 處 理 設 施 生 化 需 氧 量 三 O 化 學 需 氧 量 一 OO 流 量 大 於 二 五 O 立 方 公 尺 / 日 懸 浮 固 體 三 O 大 腸 桿 菌 群 二 OO OOO 生 化 需 氧 量 五 O 既 設 建 築 物 指 中 華 民 國 九 化 學 需

既 設 建 築 物 汙 水 處 理 設 施 生 化 需 氧 量 三 O 化 學 需 氧 量 一 OO 流 量 大 於 二 五 O 立 方 公 尺 / 日 懸 浮 固 體 三 O 大 腸 桿 菌 群 二 OO OOO 生 化 需 氧 量 五 O 既 設 建 築 物 指 中 華 民 國 九 化 學 需 台 灣 基 層 透 析 協 會 本 會 宗 旨 : 關 懷 透 析 病 患, 提 升 透 析 品 質, 維 護 透 析 診 所 與 病 患 的 權 益, 最 新 消 息 : 協 助 政 府 制 定 相 關 醫 療 政 策 及 推 行 健 康 保 險 制 度 105 年 6 月 1 日, 腎 臟 醫 學 會 林 元 灝 秘 書 長 和 基 層 透 析 協 會 顏 大 翔 秘 書 長 至 劉 建 國 立

More information

2012 2 157 No. 2 2012 Jinan Journal Philosophy and Social Sciences Sum No. 157 100871 I127. 41 A 1000-5072 2012 02-0002 - 10 4 1 1 2 1540 1503 2 2011-06 - 27 1969 1 1985 19 2 3 1997 1998 3 2 34 2 3 1540

More information

08陈会广

08陈会广 第 34 卷 第 10 期 2012 年 10 月 2012,34(10):1871-1880 Resources Science Vol.34,No.10 Oct.,2012 文 章 编 号 :1007-7588(2012)10-1871-10 房 地 产 市 场 及 其 细 分 的 调 控 重 点 区 域 划 分 理 论 与 实 证 以 中 国 35 个 大 中 城 市 为 例 陈 会 广 1,

More information

/MPa / kg m - 3 /MPa /MPa 2. 1E ~ 56 ANSYS 6 Hz (a) 一阶垂向弯曲 (b) 一阶侧向弯曲 (c) 一阶扭转 (d) 二阶侧向弯曲 (e) 二阶垂向弯曲 (f) 弯扭组合 2 6 Hz

/MPa / kg m - 3 /MPa /MPa 2. 1E ~ 56 ANSYS 6 Hz (a) 一阶垂向弯曲 (b) 一阶侧向弯曲 (c) 一阶扭转 (d) 二阶侧向弯曲 (e) 二阶垂向弯曲 (f) 弯扭组合 2 6 Hz 31 3 Vol. 31 No. 3 218 9 Journal of Shijiazhuang Tiedao University Natural Science Edition Sep. 218 1 1 2 1 2 1 1. 543 2. 543 U462. 3 217-2 - 16 A 295-373 218 3-63 - 6 1-4 5-7 8-11 1 11 11 398 mm 86 mm

More information

[1-3] (Smile) [4] 808 nm (CW) W 1 50% 1 W 1 W Fig.1 Thermal design of semiconductor laser vertical stack ; Ansys 20 bar ; bar 2 25 Fig

[1-3] (Smile) [4] 808 nm (CW) W 1 50% 1 W 1 W Fig.1 Thermal design of semiconductor laser vertical stack ; Ansys 20 bar ; bar 2 25 Fig 40 6 2011 6 Vol.40 No.6 Infrared and Laser Engineering Jun. 2011 808 nm 2000 W 1 1 1 1 2 2 2 2 2 12 (1. 710119 2. 710119) : bar 808 nm bar 100 W 808 nm 20 bar 2 000 W bar LIV bar 808 nm : : TN248.4 TN365

More information

Microsoft PowerPoint - NCBA_Cattlemens_College_Darrh_B

Microsoft PowerPoint - NCBA_Cattlemens_College_Darrh_B Introduction to Genetics Darrh Bullock University of Kentucky The Model Trait = Genetics + Environment Genetics Additive Predictable effects that get passed from generation to generation Non-Additive Primarily

More information

~ a 3 h NCEP ~ 24 3 ~ ~ 8 9 ~ km m ~ 500 m 500 ~ 800 m 800 ~ m a 200

~ a 3 h NCEP ~ 24 3 ~ ~ 8 9 ~ km m ~ 500 m 500 ~ 800 m 800 ~ m a 200 31 1 2013 3 Journal of Arid Meteorology Vol. 31 No. 1 March 2013. J. 2013 31 1 62-69. doi 10. 11755 /j. issn. 1006-7639 2013-01 - 0062 211101 2000 ~ 2010 3 h 2. 5 2. 5 NCEP 2 5 ~ 7 02 11 05 ~ 08 2 ~ 8

More information

10384 X0115019 UDC (MBA) 2004 5 2004 6 2004 XTC An Research on Internationalization Strategy of XTC , XTC XTC XTC APT XTC XTC XTC XTC XTC XTC : Abstract Abstract Although it s well known that China s

More information

m m m ~ mm

m m m ~ mm 2011 10 10 157 JOURNAL OF RAILWAY ENGINEERING SOCIETY Oct 2011 NO. 10 Ser. 157 1006-2106 2011 10-0007 - 0124-05 710043 6 TBM TBM U455. 43 A Structural Calculation and Analysis of Transfer Node of Three

More information

<4D6963726F736F667420576F7264202D20B5DAC8FDB7BDBE57C9CFD6A7B8B6D6AEB7A8C2C98696EE7DCCBDBEBF2E646F63>

<4D6963726F736F667420576F7264202D20B5DAC8FDB7BDBE57C9CFD6A7B8B6D6AEB7A8C2C98696EE7DCCBDBEBF2E646F63> 題 目 : 第 三 方 網 上 支 付 之 法 律 問 題 探 究 Title:A study on legal issues of the third-party online payment 姓 名 Name 學 號 Student No. 學 院 Faculty 課 程 Program 專 業 Major 指 導 老 師 Supervisor 日 期 Date : 王 子 瑜 : 1209853J-LJ20-0021

More information

标题

标题 第 48 卷 第 6 期 2 1 6 年 6 月 哈 尔 滨 工 业 大 学 学 报 JOURNAL OF HARBIN INSTITUTE OF TECHNOLOGY Vol 48 No 6 Jun. 216 doi:1.11918 / j.issn.367 6234.216.6.1 大 跨 度 楼 盖 结 构 在 运 动 荷 载 下 的 振 动 性 能 杨 维 国 1, 马 伯 涛 2, 宋 毛

More information

Microsoft Word doc

Microsoft Word doc 中 考 英 语 科 考 试 标 准 及 试 卷 结 构 技 术 指 标 构 想 1 王 后 雄 童 祥 林 ( 华 中 师 范 大 学 考 试 研 究 院, 武 汉,430079, 湖 北 ) 提 要 : 本 文 从 结 构 模 式 内 容 要 素 能 力 要 素 题 型 要 素 难 度 要 素 分 数 要 素 时 限 要 素 等 方 面 细 致 分 析 了 中 考 英 语 科 试 卷 结 构 的

More information

by industrial structure evolution from 1952 to 2007 and its influence effect was first acceleration and then deceleration second the effects of indust

by industrial structure evolution from 1952 to 2007 and its influence effect was first acceleration and then deceleration second the effects of indust 2011 2 1 1 2 3 4 1. 100101 2. 100124 3. 100039 4. 650092 - - - 3 GDP U 20-30 60% 10% TK01 A 1002-9753 2011 02-0042 - 10 Analysis on Character and Potential of Energy Saving and Carbon Reducing by Structure

More information

: 23 S00017242 1 -----------------------------------------------------------------------------1 -----------------------------------------------------------------------------3 -------------------------------------------------------------------7

More information

% GIS / / Fig. 1 Characteristics of flood disaster variation in suburbs of Shang

% GIS / / Fig. 1 Characteristics of flood disaster variation in suburbs of Shang 20 6 2011 12 JOURNAL OF NATURAL DISASTERS Vol. 20 No. 6 Dec. 2011 1004-4574 2011 06-0094 - 05 200062 1949-1990 1949 1977 0. 8 0. 03345 0. 01243 30 100 P426. 616 A Risk analysis of flood disaster in Shanghai

More information

Microsoft PowerPoint - CH 04 Techniques of Circuit Analysis

Microsoft PowerPoint - CH 04 Techniques of Circuit Analysis Chap. 4 Techniques of Circuit Analysis Contents 4.1 Terminology 4.2 Introduction to the Node-Voltage Method 4.3 The Node-Voltage Method and Dependent Sources 4.4 The Node-Voltage Method: Some Special Cases

More information

405 急 性 心 肌 梗 死 是 临 床 较 为 常 见 的 心 血 管 疾 病, 病 情 危 急, 病 死 率 高 [1] 随 着 经 皮 冠 状 动 脉 介 入 治 疗 (percutaneous coronary intervention,pci) 技 术 在 急 性 心 肌 梗 死 急 诊

405 急 性 心 肌 梗 死 是 临 床 较 为 常 见 的 心 血 管 疾 病, 病 情 危 急, 病 死 率 高 [1] 随 着 经 皮 冠 状 动 脉 介 入 治 疗 (percutaneous coronary intervention,pci) 技 术 在 急 性 心 肌 梗 死 急 诊 404 冠 心 舒 通 胶 囊 预 防 急 性 心 肌 梗 死 PCI 术 后 再 狭 窄 的 临 床 研 究 李 静, 王 建 榜, 王 西 辉 西 安 医 学 院 第 二 附 属 医 院 心 内 科, 陕 西 西 安 710038 摘 要 : 目 的 研 究 冠 心 舒 通 胶 囊 预 防 急 性 心 肌 梗 死 经 皮 冠 状 动 脉 介 入 治 疗 (PCI) 术 后 再 狭 窄 的 临 床

More information

问 她! 我 们 把 这 只 手 机 举 起 来 借 着 它 的 光 看 到 了 我 老 婆 正 睁 着 双 眼 你 在 干 什 么 我 问, 我 开 始 想 她 至 少 是 闭 着 眼 睛 在 yun 酿 睡 意 的 我 睡 不 着 她 很 无 辜 地 看 着 我 我 问 她 yun 酿 的 yu

问 她! 我 们 把 这 只 手 机 举 起 来 借 着 它 的 光 看 到 了 我 老 婆 正 睁 着 双 眼 你 在 干 什 么 我 问, 我 开 始 想 她 至 少 是 闭 着 眼 睛 在 yun 酿 睡 意 的 我 睡 不 着 她 很 无 辜 地 看 着 我 我 问 她 yun 酿 的 yu 果 皮 云 写 作 NO.6: 响 水 不 滚 滚 水 不 响 时 间 :2011 年 12 月 25 日 主 编 : 乌 青 作 者 : 秦 留, 新 a, 伊 文 达, 乌 青, 张 墩 墩, 娜 娜, 女 斑 马 王, 马 其 顿 荒 原, 尼 码, 萨 尔 卡, 傀 儡 尫 仔, 东 成, 二 天, 老 马 迷 途, 曾 骞, 郑 在, 柚 子, 以 下 简 称 刘 某, 大 棋, 张 维,

More information

12-1b T Q235B ML15 Ca OH Table 1 Chemical composition of specimens % C Si Mn S P Cr Ni Fe

12-1b T Q235B ML15 Ca OH Table 1 Chemical composition of specimens % C Si Mn S P Cr Ni Fe * - - 100084 Q235B ML15 Ca OH 2 DOI 10. 13204 /j. gyjz201508023 STUDY OF GALVANIC CORROSION SENSITIVITY BETWEEN ANY COUPLE OF STUD WELDMENT OR BEAM Lu Xinying Li Yang Li Yuanjin Department of Civil Engineering

More information

Research of numerical simulation of high strength steel welding residual stress and fatigue life By Chen Song

Research of numerical simulation of high strength steel welding residual stress and fatigue life By Chen Song Research of numerical simulation of high strength steel welding residual stress and fatigue life By Chen Song I ABSTRACT They are very necessary and important to carry on the research on the welding residual

More information

2005 5,,,,,,,,,,,,,,,,, , , 2174, 7014 %, % 4, 1961, ,30, 30,, 4,1976,627,,,,, 3 (1993,12 ),, 2

2005 5,,,,,,,,,,,,,,,,, , , 2174, 7014 %, % 4, 1961, ,30, 30,, 4,1976,627,,,,, 3 (1993,12 ),, 2 3,,,,,, 1872,,,, 3 2004 ( 04BZS030),, 1 2005 5,,,,,,,,,,,,,,,,, 1928 716,1935 6 2682 1928 2 1935 6 1966, 2174, 7014 %, 94137 % 4, 1961, 59 1929,30, 30,, 4,1976,627,,,,, 3 (1993,12 ),, 2 , :,,,, :,,,,,,

More information

LaDefense Arch Petronas Towers 2009 CCTV MOMA Newmark Hahn Liu 8 Heredia - Zavoni Barranco 9 Heredia - Zavoni Leyva

LaDefense Arch Petronas Towers 2009 CCTV MOMA Newmark Hahn Liu 8 Heredia - Zavoni Barranco 9 Heredia - Zavoni Leyva 39 6 2011 12 Journal of Fuzhou University Natural Science Edition Vol 39 No 6 Dec 2011 DOI CNKI 35-1117 /N 20111220 0901 002 1000-2243 2011 06-0923 - 07 350108 105 m 14 69% TU311 3 A Seismic analysis of

More information

C doc

C doc No. C2004010 2004-7 No. C2004010 2004 7 16 1 No. C2004010 2004 7 16 2000 1990 1990 2000 ( ),, 2 / 2000 1990 1990 2000 ( ),, (1952-) 100871 It is not Appropriate to Remove the Birth Spacing Policy Now,

More information

厦 门 大 学 学 位 论 文 原 创 性 声 明 本 人 呈 交 的 学 位 论 文 是 本 人 在 导 师 指 导 下, 独 立 完 成 的 研 究 成 果 本 人 在 论 文 写 作 中 参 考 其 他 个 人 或 集 体 已 经 发 表 的 研 究 成 果, 均 在 文 中 以 适 当 方

厦 门 大 学 学 位 论 文 原 创 性 声 明 本 人 呈 交 的 学 位 论 文 是 本 人 在 导 师 指 导 下, 独 立 完 成 的 研 究 成 果 本 人 在 论 文 写 作 中 参 考 其 他 个 人 或 集 体 已 经 发 表 的 研 究 成 果, 均 在 文 中 以 适 当 方 学 校 编 码 :10384 分 类 号 密 级 学 号 : UDC 硕 士 学 位 论 文 浙 江 省 人 事 考 试 突 发 事 件 应 对 策 略 探 析 An Exploration of Zhejiang Province Personnel Examination Emergency Strategy 姜 海 峰 指 导 教 师 姓 名 : 王 玉 琼 教 授 专 业 名 称 : 公 共

More information

/3 CAD JPG GIS CAD GIS GIS 1 a CAD CAD CAD GIS GIS ArcGIS 9. x 10 1 b 1112 CAD GIS 1 c R2VArcscan CAD MapGIS CAD 1 d CAD U

/3 CAD JPG GIS CAD GIS GIS 1 a CAD CAD CAD GIS GIS ArcGIS 9. x 10 1 b 1112 CAD GIS 1 c R2VArcscan CAD MapGIS CAD 1 d CAD U 1006-3862 2010 05-0059 - 07 361005 1 GIS 2 3 What if 2. 0 1 2 3 4 GIS TU984. 11 A 1 Planning Support System MIS PSS 1989 1 90 23 4-7 GIS Planning Support GIS System SDSS PSS GIS GIS CAD GIS SDSS CAD CAD

More information

ABSTRACT ABSTRACT Based on analyzing public corporation in foreign countries, this paper studies basic theories of public legal establishment, with our country s reality in the social transferring period

More information

相 关 报 道 近 年 来 临 床 上 运 用 多 种 针 灸 疗 法 治 疗 肥 胖 及 其 引 起 的 并 发 症, 如 高 脂 血 症 糖 尿 病 高 血 压 及 痛 经 等, 均 获 得 满 意 疗 效. caused by sedentary lifestyle and genetic f

相 关 报 道 近 年 来 临 床 上 运 用 多 种 针 灸 疗 法 治 疗 肥 胖 及 其 引 起 的 并 发 症, 如 高 脂 血 症 糖 尿 病 高 血 压 及 痛 经 等, 均 获 得 满 意 疗 效. caused by sedentary lifestyle and genetic f 在 线 投 稿 : http://www.baishideng.com/wcjd/ch/index.aspx 帮 助 平 台 : http://www.wjgnet.com/esps/helpdesk.aspx DOI: 10.11569/wcjd.v24.i5.815 世 界 华 人 消 化 杂 志 2016 年 2 月 18 日 ; 24(5): 815-820 ISSN 1009-3079 (print)

More information

Revit Revit Revit BIM BIM 7-9 3D 1 BIM BIM 6 Revit 0 4D 1 2 Revit Revit 2. 1 Revit Revit Revit Revit 2 2 Autodesk Revit Aut

Revit Revit Revit BIM BIM 7-9 3D 1 BIM BIM 6 Revit 0 4D 1 2 Revit Revit 2. 1 Revit Revit Revit Revit 2 2 Autodesk Revit Aut 60 2 2016 2 RAILWAY STANDARD DESIGN Vol. 60 No. 2 Feb. 2016 1004-2954201602-0071-06 BIM 1 1 2 2 1 1. 7140992. 710054 BIM BIM 3D 4D nd BIM 1 3D 4D Revit BIM BIM U442. 5TP391. 72 A DOI10. 13238 /j. issn.

More information

2015年4月11日雅思阅读预测机经(新东方版)

2015年4月11日雅思阅读预测机经(新东方版) 剑 桥 雅 思 10 第 一 时 间 解 析 阅 读 部 分 1 剑 桥 雅 思 10 整 体 内 容 统 计 2 剑 桥 雅 思 10 话 题 类 型 从 以 上 统 计 可 以 看 出, 雅 思 阅 读 的 考 试 话 题 一 直 广 泛 多 样 而 题 型 则 稳 中 有 变 以 剑 桥 10 的 test 4 为 例 出 现 的 三 篇 文 章 分 别 是 自 然 类, 心 理 研 究 类,

More information

2015 1 330 GLOBAL EDUCATION Vol. 44 No1 2015 / 300387 / 300387 20 20 1926 20 1912 1912 1918 1919 34 1919 1921 1 1922 1923 1925 7 100 2 1922 3 1923 5 1923 1924 4 1927 1949 20 1. 20 5 35 2. 20 20 1927 1949

More information

Fig. 1 Frame calculation model 1 mm Table 1 Joints displacement mm

Fig. 1 Frame calculation model 1 mm Table 1 Joints displacement mm 33 2 2011 4 ol. 33 No. 2 Apr. 2011 1002-8412 2011 02-0104-08 1 1 1 2 361003 3. 361009 3 1. 361005 2. GB50023-2009 TU746. 3 A Study on Single-span RC Frame Reinforced with Steel Truss System Yuan Xing-ren

More information

1 ABSTRACT 3 5 13 14 14 15 15 16 17 18 18 18 18-21 28 33 34 46 47 53 54 2

1 ABSTRACT 3 5 13 14 14 15 15 16 17 18 18 18 18-21 28 33 34 46 47 53 54 2 : - S99017212 O O 1 1 ABSTRACT 3 5 13 14 14 15 15 16 17 18 18 18 18-21 28 33 34 46 47 53 54 2 - A L Aeq CNE 463 3-5 L Aeq L Aeq CNE WHO/ISH 1999 [49] 18.6kPa / 12.0kPa EpiInfo6.04 SPSS10.0 Logistic 10

More information

untitled

untitled and Due Diligence M&A in China Prelude and Due Diligence A Case For Proper A Gentleman s Agreement? 1 Respect for the Rule of Law in China mandatory under law? CRITICAL DOCUMENTS is driven by deal structure:

More information

作 主 动 追 求 知 识 获 取 技 能, 在 心 理 和 生 理 上 都 非 常 积 极 的 个 体 (Zimmerman & Pons, 1986) 在 此 期 间, 自 我 效 能 感 (self-efficacy) 自 我 控 制 (self-control) 自 我 管 理 (self-

作 主 动 追 求 知 识 获 取 技 能, 在 心 理 和 生 理 上 都 非 常 积 极 的 个 体 (Zimmerman & Pons, 1986) 在 此 期 间, 自 我 效 能 感 (self-efficacy) 自 我 控 制 (self-control) 自 我 管 理 (self- 网 络 环 境 中 英 语 学 习 者 自 我 调 控 的 个 体 特 征 研 究 A Qualitative Study of Chinese EFL Learners Online Self-regulation * 1,2 郑 春 萍, 穆 婕 2, 章 僖 格 2, 苑 仁 庆 2, 李 芒 1 北 京 师 范 大 学 教 育 学 部 2 北 京 邮 电 大 学 人 文 学 院 * zhengchunping@bupt.edu.cn

More information

A Study on the Relationships of the Co-construction Contract A Study on the Relationships of the Co-Construction Contract ( ) ABSTRACT Co-constructio in the real estate development, holds the quite

More information

<4D6963726F736F667420576F7264202D203338B4C12D42A448A4E5C3C0B34EC3FE2DAB65ABE1>

<4D6963726F736F667420576F7264202D203338B4C12D42A448A4E5C3C0B34EC3FE2DAB65ABE1> ϲ ฯ र ቑ ጯ 高雄師大學報 2015, 38, 63-93 高雄港港史館歷史變遷之研究 李文環 1 楊晴惠 2 摘 要 古老的建築物往往承載許多回憶 也能追溯某些歷史發展的軌跡 位於高雄市蓬 萊路三號 現為高雄港港史館的紅磚式建築 在高雄港三號碼頭作業區旁的一片倉庫 群中 格外搶眼 這棟建築建成於西元 1917 年 至今已將近百年 不僅躲過二戰戰 火無情轟炸 並保存至今 十分可貴 本文透過歷史考證

More information

Microsoft PowerPoint - ryz_030708_pwo.ppt

Microsoft PowerPoint - ryz_030708_pwo.ppt Long Term Recovery of Seven PWO Crystals Ren-yuan Zhu California Institute of Technology CMS ECAL Week, CERN Introduction 20 endcap and 5 barrel PWO crystals went through (1) thermal annealing at 200 o

More information

國家圖書館典藏電子全文

國家圖書館典藏電子全文 A Study on the Job Stress and the Ways of Coping for the Director of Elementary School in the Middle Area of Taiwan Abstract This study aims at probing the subject current status as related to stress and

More information

10384 27720071152270 UDC SHIBOR - Research on Dynamics of Short-term Shibor via Parametric and Nonparametric Models 2 0 1 0 0 5 2 0 1 0 0 5 2 0 1 0 0 5 2010 , 1. 2. Shibor 2006 10 8 2007 1 4 Shibor

More information

Microsoft PowerPoint - ATF2015.ppt [相容模式]

Microsoft PowerPoint - ATF2015.ppt [相容模式] Improving the Video Totalized Method of Stopwatch Calibration Samuel C.K. Ko, Aaron Y.K. Yan and Henry C.K. Ma The Government of Hong Kong Special Administrative Region (SCL) 31 Oct 2015 1 Contents Introduction

More information

<4D6963726F736F667420576F7264202D2035B171AB73B6CBA8ECAB73A6D3A4A3B6CBA158B3AFA46CA9F9BB50B169A445C4D6AABAB750B94AB8D6B9EFA4F1ACE3A873>

<4D6963726F736F667420576F7264202D2035B171AB73B6CBA8ECAB73A6D3A4A3B6CBA158B3AFA46CA9F9BB50B169A445C4D6AABAB750B94AB8D6B9EFA4F1ACE3A873> 中 正 漢 學 研 究 2012 年 第 一 期 ( 總 第 十 九 期 ) 2012 年 6 月 頁 111~134 國 立 中 正 大 學 中 國 文 學 系 111 從 哀 傷 到 哀 而 不 傷 : 陳 子 昂 與 張 九 齡 的 感 遇 詩 對 比 研 究 * 丁 涵 摘 要 在 中 國 古 典 文 學 語 境 中, 一 個 主 題 的 奠 立 往 往 需 要 歷 時 彌 久, 而 這 本

More information