An efficient multi-harmonic method is proposed for studying the effects of mistuning on resonant features of bladed disks with blade-to-blade dry friction damping. This method is able to predict accurately the forced ...An efficient multi-harmonic method is proposed for studying the effects of mistuning on resonant features of bladed disks with blade-to-blade dry friction damping. This method is able to predict accurately the forced response of bladed disks in frequency domain, which is validated by numerical integration method in time domain. The resonant features of both tuned and mistuned systems are investigated by using this method under various system coupling strengths, viscous dampings, and dry friction darnpings, etc. The results demonstrate that the proposed multi-harmonic method is very efficient for studying the mistuning effects on the resonant response of bladed disks with blade-to-blade dry friction damping, especially considering the combined effects of various system parameters.展开更多
In this paper, the exact analytical and numerical solutions for rotating variable-thickness annular disk are presented. The inner and outer edges of the rotating variable-thickness annular disk are considered to have ...In this paper, the exact analytical and numerical solutions for rotating variable-thickness annular disk are presented. The inner and outer edges of the rotating variable-thickness annular disk are considered to have free boundary conditions. Two different annular disks for the radially varying thickness are given. The numerical Runge-Kutta solution as well as the exact analytical solution is available for the first disk while the exact analytical solution is not available for the second annular disk. Both exact and numerical results for stress function, stresses, strains and radial displacement will be investigated for the first annular disk of variable thickness. The accuracy of the present numerical solution is discussed and its ability of use for the second rotating variable-thickness annular disk is investigated. Finally, the distributions of stress function, displacement, strains, and stresses will be presented. The appropriate comparisons and discussions are made at the same angular velocity.展开更多
More than a decade of dedicated experimental work on the collisional physics of protoplanetary dust has brought us to a point at which the growth of dust aggregates can - for the first time - be self-consistently and ...More than a decade of dedicated experimental work on the collisional physics of protoplanetary dust has brought us to a point at which the growth of dust aggregates can - for the first time - be self-consistently and reliably modeled. In this article, the emergent collision model for protoplanetery dust aggregates, as well as the numerical model for the evolution of dust aggregates in protoplanetary disks, is reviewed. It turns out that, after a brief period of rapid collisional growth of fluffy dust aggregates to sizes of a few centimeters, the protoplanetary dust particles are subject to bouncing collisions, in which their porosity is considerably decreased. The model results also show that low-velocity fragmentation can reduce the final mass of the dust aggregates but that it does not trigger a new growth mode as discussed previously. According to the current stage of our model, the direct formation of kilometer-sized planetesimals by collisional sticking seems unlikely, implying that collective effects, such as the streaming instability and the gravitational instability in dust-enhanced regions of the protoplanetary disk, are the best candidates for the processes leading to planetesimals.展开更多
In this paper, the analytical and numerical solutions for rotating variable-thickness solid disk and numerical solution for rotating variable-thickness annular disk are presented. The outer edge of the solid disk and ...In this paper, the analytical and numerical solutions for rotating variable-thickness solid disk and numerical solution for rotating variable-thickness annular disk are presented. The outer edge of the solid disk and the inner and outer edges of the annular disk are considered to have clamped boundary conditions. Two different cases for the radially varying thickness of the solid and annular disks are given. The numerical solution as well as the analytical solution is available for the first case of the solid disk while the analytical solution is not available for the second case of the annular disk. Both analytical and numerical results for displacement and stresses will be investigated for the first case of radially varying thickness. The accuracy of the present numerical solution is discussed and its ability of use for the second case of radially varying thickness is investigated. Finally, the distributions of displacement and stresses will be presented and the appropriate comparisons and discussions are made at the same angular velocity.展开更多
The steady, laminar, incompressible and two dimensional micropolar flow between two porous disks was investigated using optimal homotopy asymptotic method(OHAM) and fourth order Runge–Kutta numerical method. Comparis...The steady, laminar, incompressible and two dimensional micropolar flow between two porous disks was investigated using optimal homotopy asymptotic method(OHAM) and fourth order Runge–Kutta numerical method. Comparison between OHAM and numerical method shows that OHAM is an exact and high efficient method for solving these kinds of problems. The results are presented to study the velocity and rotation profiles for different physical parameters such as Reynolds number, vortex viscosity parameter, spin gradient viscosity and microinertia density parameter. As an important outcome, the magnitude of the microrotation increases with an increase in the values of injection velocity while it decreases by increasing the values of suction velocity.展开更多
Stress intensity factors (SIFs) for the cracked circular disks under different distributing surface tractions are evaluated with the scaled boundary finite element method (SBFEM). In the SBFEM, the analytical adva...Stress intensity factors (SIFs) for the cracked circular disks under different distributing surface tractions are evaluated with the scaled boundary finite element method (SBFEM). In the SBFEM, the analytical advantage of the solution in the radial direction allows SIFs to be directly determined from its definition, therefore no special crack-tip treatment is necessary. Furthermore anisotropic material behavior can be treated easily. Different distributions of surface tractions are considered for the center and double-edge-cracked disks. The benchmark examples are modeled and an excellent agreement between the results in the present study and those in published literature is found. It shows that SBFEM is effective and possesses high accuracy. The SIFs of the cracked orthotropic material circular disks subjected to different surface tractions are also evaluated. The technique of substructure is applied to handle the multiple cracks problem.展开更多
We study the interaction between supermassive binary black holes in an elliptical orbit and their surrounding disk with a gap. The gap in the disk is a low density region formed due to the tidal effects of the less ma...We study the interaction between supermassive binary black holes in an elliptical orbit and their surrounding disk with a gap. The gap in the disk is a low density region formed due to the tidal effects of the less massive black hole. The binary we have investigated has a sub-parsec separation and is coplanar with the disk. We find that the maximum variation of the surface density in the gap reaches 50% during an orbital period. However, in other regions of the disk, the density variation is much less than 1%. Furthermore, we calculate the corresponding variation of spectral energy distribution within a period, but little variation is found. The reason for these results is that the viscosity timescale of the disk at the binary radius is much longer than the orbital period of the binary.展开更多
We present a statistical method to derive the stellar density profiles of the Milky Way from spectroscopic survey data, taking into account selection effects. We assume the selection function, which can be altered dur...We present a statistical method to derive the stellar density profiles of the Milky Way from spectroscopic survey data, taking into account selection effects. We assume the selection function, which can be altered during observations and data reductions, of the spectroscopic survey is based on photometric colors and magnitude. Then the underlying selection function for a line-of-sight can be recovered well by comparing the distribution of the spectroscopic stars in a color-magnitude plane with that of the photometric dataset. Subsequently, the stellar density profile along a line-of-sight can be derived from the spectroscopically measured stellar density profile multiplied by the selection function. The method is validated using Galaxia mock data with two different selection functions. We demonstrate that the derived stellar density profiles reconstruct the true ones well not only for the full set of targets, but also for sub-populations selected from the full dataset. Finally, the method is applied to map the density pro- files for the Galactic disk and halo, using the LAMOST RGB stars. The Galactic disk extends to about R = 19 kpc, where the disk still contributes about 10% to the total stellar surface density. Beyond this radius, the disk smoothly transitions to the halo without any truncation, bending or breaking. Moreover, no over-density corresponding to the Monoceros ring is found in the Galactic anti-center direction. The disk shows moderate north-south asymmetry at radii larger than 12 kpc. On the other hand, the R-Z tomographic map directly shows that the stellar halo is substantially oblate within a Galactocentric radius of 20 kpc and gradually becomes nearly spherical beyond 30 kpc.展开更多
This paper describes the robust optimum design which combines the geometrical optimization method proposed by Hashimoto and statistical method. Recently, 2.5″ hard disk drives (HDDs) are widely used for mobile device...This paper describes the robust optimum design which combines the geometrical optimization method proposed by Hashimoto and statistical method. Recently, 2.5″ hard disk drives (HDDs) are widely used for mobile devices such as laptops, video cameras and car navigation systems. In mobile applications, high durability towards external vibrations and shocks are essentials to the bearings of HDD spindle motor. In addition, the bearing characteristics are influenced by manufacturing error because of small size of the bearings of HDD. In this paper, the geometrical optimization is carried out to maximize the bearing stiffness using sequential quadratic programming to improve vibration characteristics. Additionally, the bearing stiffness is analyzed considering dimensional tolerance of the bearing using statistical method. The dimensional tolerance is assumed to distribute according to the Gaussian distribution, and then the bearing stiffness is estimated by combining the expectation and standard deviation. As a result, in the robust optimum design, new groove geometry of bearing can be obtained in which the bearing stiffness is four times higher than the stiffness of conventional spiral groove bearing. Moreover, the bearing has lower variability compared with the result of optimum design neglecting dimensional tolerance.展开更多
The unsteady, laminar, incompressible, and two-dimensional flow of a mi- cropolar fluid between two orthogonally moving porous coaxial disks is considered. The extension of von Karman's similarity transformations is ...The unsteady, laminar, incompressible, and two-dimensional flow of a mi- cropolar fluid between two orthogonally moving porous coaxial disks is considered. The extension of von Karman's similarity transformations is used to reduce the governing partial differential equations (PDEs) to a set of non-linear coupled ordinary differential equations (ODEs) in the dimensionless form. The analytical solutions are obtained by employing the homotopy analysis method (HAM). The effects of various physical param- eters such as the expansion ratio and the permeability Reynolds number on the velocity fields are discussed in detail.展开更多
Reliability analysis methods based on the linear damage accumulation law (LDAL) and load-life interference model are studied in this paper. According to the equal probability rule, the equivalent loads are derived, an...Reliability analysis methods based on the linear damage accumulation law (LDAL) and load-life interference model are studied in this paper. According to the equal probability rule, the equivalent loads are derived, and the reliability analysis method based on load-life interference model and recurrence formula is constructed. In conjunction with finite element analysis (FEA) program, the reliability of an aero engine turbine disk under low cycle fatigue (LCF) condition has been analyzed. The results show the turbine disk is safety and the above reliability analysis methods are feasible.展开更多
This paper studies the flow structural interaction (FSI) within a hard disk drive (HDD) through the use of a novel coupling method. The interaction studied was the fluid induced vibration in the HDD. A two step coupli...This paper studies the flow structural interaction (FSI) within a hard disk drive (HDD) through the use of a novel coupling method. The interaction studied was the fluid induced vibration in the HDD. A two step coupling approach was used, where the fluid and structural components were solved sequentially. The result obtained was a ratio of 0.65 between the vibration amplitudes of a fixed head stack assembly (HSA) and a moving HSA. The ratio was next applied on a real 3.5 inch HDD, to allow the parameter to be further improved upon. A new benchmark index of 0.69 was developed from this. This parameter may allow future researchers to model the out of plane vibrations of a HSA easily, saving precious time. A 31% more accurate simulation of FSI within 3.5 inch HDD at 15000 rpm is achieved by the use of this new coupling method and benchmark index.展开更多
This paper presents the analytical and numerical solutions for a rotating variable-thickness solid disk. The outer edge of the solid disk is considered to have free boundary conditions. The governing equation is deriv...This paper presents the analytical and numerical solutions for a rotating variable-thickness solid disk. The outer edge of the solid disk is considered to have free boundary conditions. The governing equation is derived from the basic equations of the rotating solid disk and it is solved analytically or numerically using finite difference algorithm. Both analytical and numerical results for the distributions of stress function and stresses of variable-thickness solid disks are obtained. Finally, the distributions of stress function and stresses are presented and the appropriate comparisons and discussions are made at the same angular velocity.展开更多
On the basis of the data of hydraulic fracturing stress measurement and mechanical property of the drilledcore from the Maoping 800 m bore hole, the failure mechanism of disk core in the bore hole was analysed. Inthis...On the basis of the data of hydraulic fracturing stress measurement and mechanical property of the drilledcore from the Maoping 800 m bore hole, the failure mechanism of disk core in the bore hole was analysed. Inthis paper, it is considered that the disk core is caused by the declination of principal axis of crustal stress under low crustal stress value. According to the section property of disk core and stress analysis, the form of disk core is tensile failure.展开更多
Fused deposition modelling(FDM), a widely used rapid prototyping process, is a promising technique in manufacturing engineering. In this work, a method for characterizing elastic constants of FDM-fabricated materials ...Fused deposition modelling(FDM), a widely used rapid prototyping process, is a promising technique in manufacturing engineering. In this work, a method for characterizing elastic constants of FDM-fabricated materials is proposed. First of all, according to the manufacturing process of FDM, orthotropic constitutive model is used to describe the mechanical behavior. Then the virtual fields method(VFM) is applied to characterize all the mechanical parameters(Q, Q, Q, Q) using the full-field strain,which is measured by digital image correlation(DIC). Since the principal axis of the FDM-fabricated structure is sometimes unknown due to the complexity of the manufacturing process, a disk in diametrical compression is used as the load configuration so that the loading angle can be changed conveniently. To verify the feasibility of the proposed method, finite element method(FEM) simulation is conducted to obtain the strain field of the disk. The simulation results show that higher accuracy can be achieved when the loading angle is close to 30?. Finally, a disk fabricated by FDM was used for the experiment. By rotating the disk, several tests with different loading angles were conducted. To determine the position of the principal axis in each test, two groups of parameters(Q, Q, Q, Q) are calculated by two different groups of virtual fields. Then the corresponding loading angle can be determined by minimizing the deviation between two groups of the parameters. After that, the four constants(Q, Q, Q, Q) were determined from the test with an angle of 27?.展开更多
基金National Nature Science Foundation of China (NO.50275121)
文摘An efficient multi-harmonic method is proposed for studying the effects of mistuning on resonant features of bladed disks with blade-to-blade dry friction damping. This method is able to predict accurately the forced response of bladed disks in frequency domain, which is validated by numerical integration method in time domain. The resonant features of both tuned and mistuned systems are investigated by using this method under various system coupling strengths, viscous dampings, and dry friction darnpings, etc. The results demonstrate that the proposed multi-harmonic method is very efficient for studying the mistuning effects on the resonant response of bladed disks with blade-to-blade dry friction damping, especially considering the combined effects of various system parameters.
文摘In this paper, the exact analytical and numerical solutions for rotating variable-thickness annular disk are presented. The inner and outer edges of the rotating variable-thickness annular disk are considered to have free boundary conditions. Two different annular disks for the radially varying thickness are given. The numerical Runge-Kutta solution as well as the exact analytical solution is available for the first disk while the exact analytical solution is not available for the second annular disk. Both exact and numerical results for stress function, stresses, strains and radial displacement will be investigated for the first annular disk of variable thickness. The accuracy of the present numerical solution is discussed and its ability of use for the second rotating variable-thickness annular disk is investigated. Finally, the distributions of stress function, displacement, strains, and stresses will be presented. The appropriate comparisons and discussions are made at the same angular velocity.
基金funded by the German Space Agency (DLR) under grant Nos. 50WM0336, 50WM0636 and 50WM0936the Deutsche Forschungsgemeinschaft (DFG) under grant No. Bl298/7-1
文摘More than a decade of dedicated experimental work on the collisional physics of protoplanetary dust has brought us to a point at which the growth of dust aggregates can - for the first time - be self-consistently and reliably modeled. In this article, the emergent collision model for protoplanetery dust aggregates, as well as the numerical model for the evolution of dust aggregates in protoplanetary disks, is reviewed. It turns out that, after a brief period of rapid collisional growth of fluffy dust aggregates to sizes of a few centimeters, the protoplanetary dust particles are subject to bouncing collisions, in which their porosity is considerably decreased. The model results also show that low-velocity fragmentation can reduce the final mass of the dust aggregates but that it does not trigger a new growth mode as discussed previously. According to the current stage of our model, the direct formation of kilometer-sized planetesimals by collisional sticking seems unlikely, implying that collective effects, such as the streaming instability and the gravitational instability in dust-enhanced regions of the protoplanetary disk, are the best candidates for the processes leading to planetesimals.
文摘In this paper, the analytical and numerical solutions for rotating variable-thickness solid disk and numerical solution for rotating variable-thickness annular disk are presented. The outer edge of the solid disk and the inner and outer edges of the annular disk are considered to have clamped boundary conditions. Two different cases for the radially varying thickness of the solid and annular disks are given. The numerical solution as well as the analytical solution is available for the first case of the solid disk while the analytical solution is not available for the second case of the annular disk. Both analytical and numerical results for displacement and stresses will be investigated for the first case of radially varying thickness. The accuracy of the present numerical solution is discussed and its ability of use for the second case of radially varying thickness is investigated. Finally, the distributions of displacement and stresses will be presented and the appropriate comparisons and discussions are made at the same angular velocity.
文摘The steady, laminar, incompressible and two dimensional micropolar flow between two porous disks was investigated using optimal homotopy asymptotic method(OHAM) and fourth order Runge–Kutta numerical method. Comparison between OHAM and numerical method shows that OHAM is an exact and high efficient method for solving these kinds of problems. The results are presented to study the velocity and rotation profiles for different physical parameters such as Reynolds number, vortex viscosity parameter, spin gradient viscosity and microinertia density parameter. As an important outcome, the magnitude of the microrotation increases with an increase in the values of injection velocity while it decreases by increasing the values of suction velocity.
基金financially supported by the National Natural Youth Foundation of China (Grant Nos. 51109134,51009019, 11102118 and 51208310)the Liaoning Province Education Administration Foundation (Grant No. L2010413)+1 种基金the China Postdoctoral Science Foundation (Grant No. 2011M500557)the Natural Science Foundation of Liaoning Province (Grant No.20102164)
文摘Stress intensity factors (SIFs) for the cracked circular disks under different distributing surface tractions are evaluated with the scaled boundary finite element method (SBFEM). In the SBFEM, the analytical advantage of the solution in the radial direction allows SIFs to be directly determined from its definition, therefore no special crack-tip treatment is necessary. Furthermore anisotropic material behavior can be treated easily. Different distributions of surface tractions are considered for the center and double-edge-cracked disks. The benchmark examples are modeled and an excellent agreement between the results in the present study and those in published literature is found. It shows that SBFEM is effective and possesses high accuracy. The SIFs of the cracked orthotropic material circular disks subjected to different surface tractions are also evaluated. The technique of substructure is applied to handle the multiple cracks problem.
基金Supported by the National Natural Science Foundation of China
文摘We study the interaction between supermassive binary black holes in an elliptical orbit and their surrounding disk with a gap. The gap in the disk is a low density region formed due to the tidal effects of the less massive black hole. The binary we have investigated has a sub-parsec separation and is coplanar with the disk. We find that the maximum variation of the surface density in the gap reaches 50% during an orbital period. However, in other regions of the disk, the density variation is much less than 1%. Furthermore, we calculate the corresponding variation of spectral energy distribution within a period, but little variation is found. The reason for these results is that the viscosity timescale of the disk at the binary radius is much longer than the orbital period of the binary.
基金supported by the Strategic Priority Research Program“The Emergence of Cosmological Structures”of the Chinese Academy of Sciences(Grant No.XDB09000000)the National Key Basic Research Program of China(2014CB845700)+1 种基金the National Natural Science Foundation of China(Grant Nos.11373032 and 11333003)a National Major Scientific Project built by the Chinese Academy of Sciences.Funding for the project has been provided by the project has been provided by the National Development and Reform Commission
文摘We present a statistical method to derive the stellar density profiles of the Milky Way from spectroscopic survey data, taking into account selection effects. We assume the selection function, which can be altered during observations and data reductions, of the spectroscopic survey is based on photometric colors and magnitude. Then the underlying selection function for a line-of-sight can be recovered well by comparing the distribution of the spectroscopic stars in a color-magnitude plane with that of the photometric dataset. Subsequently, the stellar density profile along a line-of-sight can be derived from the spectroscopically measured stellar density profile multiplied by the selection function. The method is validated using Galaxia mock data with two different selection functions. We demonstrate that the derived stellar density profiles reconstruct the true ones well not only for the full set of targets, but also for sub-populations selected from the full dataset. Finally, the method is applied to map the density pro- files for the Galactic disk and halo, using the LAMOST RGB stars. The Galactic disk extends to about R = 19 kpc, where the disk still contributes about 10% to the total stellar surface density. Beyond this radius, the disk smoothly transitions to the halo without any truncation, bending or breaking. Moreover, no over-density corresponding to the Monoceros ring is found in the Galactic anti-center direction. The disk shows moderate north-south asymmetry at radii larger than 12 kpc. On the other hand, the R-Z tomographic map directly shows that the stellar halo is substantially oblate within a Galactocentric radius of 20 kpc and gradually becomes nearly spherical beyond 30 kpc.
文摘This paper describes the robust optimum design which combines the geometrical optimization method proposed by Hashimoto and statistical method. Recently, 2.5″ hard disk drives (HDDs) are widely used for mobile devices such as laptops, video cameras and car navigation systems. In mobile applications, high durability towards external vibrations and shocks are essentials to the bearings of HDD spindle motor. In addition, the bearing characteristics are influenced by manufacturing error because of small size of the bearings of HDD. In this paper, the geometrical optimization is carried out to maximize the bearing stiffness using sequential quadratic programming to improve vibration characteristics. Additionally, the bearing stiffness is analyzed considering dimensional tolerance of the bearing using statistical method. The dimensional tolerance is assumed to distribute according to the Gaussian distribution, and then the bearing stiffness is estimated by combining the expectation and standard deviation. As a result, in the robust optimum design, new groove geometry of bearing can be obtained in which the bearing stiffness is four times higher than the stiffness of conventional spiral groove bearing. Moreover, the bearing has lower variability compared with the result of optimum design neglecting dimensional tolerance.
基金Project supported by the National Natural Science Foundation of China (Nos. 51004013,50936003,51174028,and 50905013)the Research Foundation of Engineering Research Institute of University of Science and Technology Beijing (No. Yj2011-015)the Fundamental Research Funds for the Central Universities (No. T-RF-TP-12-108A)
文摘The unsteady, laminar, incompressible, and two-dimensional flow of a mi- cropolar fluid between two orthogonally moving porous coaxial disks is considered. The extension of von Karman's similarity transformations is used to reduce the governing partial differential equations (PDEs) to a set of non-linear coupled ordinary differential equations (ODEs) in the dimensionless form. The analytical solutions are obtained by employing the homotopy analysis method (HAM). The effects of various physical param- eters such as the expansion ratio and the permeability Reynolds number on the velocity fields are discussed in detail.
基金Supports provided by Aviation Basic Science Foundation(00B53010)Aerospace Science Foundation(N3CH0502)Shaanxi Province Natural Science Foundation(N3CS0501)are gratefully appreciated.
文摘Reliability analysis methods based on the linear damage accumulation law (LDAL) and load-life interference model are studied in this paper. According to the equal probability rule, the equivalent loads are derived, and the reliability analysis method based on load-life interference model and recurrence formula is constructed. In conjunction with finite element analysis (FEA) program, the reliability of an aero engine turbine disk under low cycle fatigue (LCF) condition has been analyzed. The results show the turbine disk is safety and the above reliability analysis methods are feasible.
文摘This paper studies the flow structural interaction (FSI) within a hard disk drive (HDD) through the use of a novel coupling method. The interaction studied was the fluid induced vibration in the HDD. A two step coupling approach was used, where the fluid and structural components were solved sequentially. The result obtained was a ratio of 0.65 between the vibration amplitudes of a fixed head stack assembly (HSA) and a moving HSA. The ratio was next applied on a real 3.5 inch HDD, to allow the parameter to be further improved upon. A new benchmark index of 0.69 was developed from this. This parameter may allow future researchers to model the out of plane vibrations of a HSA easily, saving precious time. A 31% more accurate simulation of FSI within 3.5 inch HDD at 15000 rpm is achieved by the use of this new coupling method and benchmark index.
文摘This paper presents the analytical and numerical solutions for a rotating variable-thickness solid disk. The outer edge of the solid disk is considered to have free boundary conditions. The governing equation is derived from the basic equations of the rotating solid disk and it is solved analytically or numerically using finite difference algorithm. Both analytical and numerical results for the distributions of stress function and stresses of variable-thickness solid disks are obtained. Finally, the distributions of stress function and stresses are presented and the appropriate comparisons and discussions are made at the same angular velocity.
文摘On the basis of the data of hydraulic fracturing stress measurement and mechanical property of the drilledcore from the Maoping 800 m bore hole, the failure mechanism of disk core in the bore hole was analysed. Inthis paper, it is considered that the disk core is caused by the declination of principal axis of crustal stress under low crustal stress value. According to the section property of disk core and stress analysis, the form of disk core is tensile failure.
基金the financial support from the National Natural Science Foundation of China (Grants 11672153, 11232008, and 11227801)
文摘Fused deposition modelling(FDM), a widely used rapid prototyping process, is a promising technique in manufacturing engineering. In this work, a method for characterizing elastic constants of FDM-fabricated materials is proposed. First of all, according to the manufacturing process of FDM, orthotropic constitutive model is used to describe the mechanical behavior. Then the virtual fields method(VFM) is applied to characterize all the mechanical parameters(Q, Q, Q, Q) using the full-field strain,which is measured by digital image correlation(DIC). Since the principal axis of the FDM-fabricated structure is sometimes unknown due to the complexity of the manufacturing process, a disk in diametrical compression is used as the load configuration so that the loading angle can be changed conveniently. To verify the feasibility of the proposed method, finite element method(FEM) simulation is conducted to obtain the strain field of the disk. The simulation results show that higher accuracy can be achieved when the loading angle is close to 30?. Finally, a disk fabricated by FDM was used for the experiment. By rotating the disk, several tests with different loading angles were conducted. To determine the position of the principal axis in each test, two groups of parameters(Q, Q, Q, Q) are calculated by two different groups of virtual fields. Then the corresponding loading angle can be determined by minimizing the deviation between two groups of the parameters. After that, the four constants(Q, Q, Q, Q) were determined from the test with an angle of 27?.