POLYHEDRON

所属分类:matlab编程
开发工具:matlab
文件大小:76KB
下载次数:10
上传日期:2014-04-22 09:57:45
上 传 者伍海军
说明:   采用线积分均匀任意形状的多面体源的全部重力张量的分析计算
(Analytical computation of the full gravity tensor of a homogeneous arbitrarily shaped polyhedral source using line integrals)

文件列表:
POLYHEDRON (0, 2012-03-24)
POLYHEDRON\xyzposnew (136, 2010-06-11)
POLYHEDRON\polaut2_out (13467, 2010-06-11)
POLYHEDRON\wuerfel.f (13477, 2012-03-24)
POLYHEDRON\additionaltables.pdf (67756, 2011-06-21)
POLYHEDRON\topoaut (66, 2008-07-09)
POLYHEDRON\wuerfel_dat (280, 2010-06-11)
POLYHEDRON\dataut (16, 1998-12-10)
POLYHEDRON\wuerfel2_out (18586, 2010-06-11)
POLYHEDRON\polyhedron.f (26189, 2012-03-24)
POLYHEDRON\ptw_dat (250, 1999-07-31)

The appended electronic supplement contains following 9 files: 1. polyhedron.f This file uses following input files 2. xyzposnew 3. topoaut 4. dataut and writes to the output file 5. polaut2_out Furthermore appended is the code for the computation of the gravity tensor of a right rectangular prism according to Mader (1951). The code is given by 6. wuerfel.f which uses the input files 7. wuerfel_dat 8. ptw_dat and writes to the output file 9. wuerfel2_out Note that the appended-type files 'polaut2_out' and 'wuerfel2_out' contain already 28 discrete cases as mentioned in the manuscript, where for each case file 'xyzposnew' expressing the exact geometric setting has also been added manually. For running or testing code wuerfel.f for any rectangular prismatic source you may need to remove the comments from lines 341, 391, 425 and 460 which refer to MSIMSL and read USE MSIMSL Module MSIMSL is specific to some Fortran compilers, e.g. Visual Compaq, and is necessary for computing the hyperbolic arcsine function. Most of the widely used compilers (gfortran, gfortran44 and ifort) already include all of the mathematical functions used in wuerfel.f and for these compilers the USE MSIMSL statement is redundant. For running or testing code polyhedron.f for any polyhedral source following steps are to be taken A. Insertion or 'save as' in a file named 'xyzposnew' of all vertex coordinates with respect to the computation point as X, Y, Z (the origin is located at the computation point). B. Definition/Computation of the topology matrix 'topoaut' of this source where vertex numbers should comply with the order of the vertices in file 'xyzposnew'. C. Definition of the varying number of segments building each of the polyhedral faces in file 'dataut'. File 'dataut' can be built simply by counting the vertices of each line (record) of file 'topoaut'. D. Manual editing of parameters 'nop' (number of planes) and 'nov' (number of vertices) in Line 127 of polyhedron.f with the actual total number of faces and vertices that build the source. Line 127 refers currently to the source body of Figure 2 of the paper and reads parameter(nop=7,nov=8,noed=100) Note that parameter 'noed' has been set to a conveniently large number that expresses the varying number of vertices building each face. If any of the polyhedral faces consists of more than 100 vertices (or, equally, segments) this parameter has to be changed accordingly as well. E. Compile and run 'polyhedron.f'. The results will be appended to the end of file 'polaut2_out' If you do not change anything to the present input files and simply compile and run 'polyhedron.f' and 'wuerfel.f' you should get an exact copy of the last 10 lines of files 'polaut2_out' and 'wuerfel2_out' appended at the end of these files. Dimitris Tsoulis, August 21, 2011

近期下载者

相关文件


收藏者