UDF-in-fluent

利用PHOTOSHOP软件可以量测出左右像片各控制点的像素坐标，通过公式可以计算出其像点坐标。再根据给定的控制点的空间坐标，利用直线线性变换解法计算出11个DLT参数即l的系数值以及内方位元素然后解算出外方位元素。另外一种方式是用单像空间后方交会的公式计算出外方位元素。利用单像空间后方交会的来结算影像的外方位元素的基本思想是：以单幅影像为基础，从该影像所覆盖地面范围内若干控制点的已知地面坐标和相应点的像坐标量测值出发，根据共线条件方程，解求该影像在航空摄影时刻的外方位元素。首先建立单文档的程序，建立一个对话框类，在这个类中进行像素坐标、控制点坐标文件的打开与读取，然后利用DLT公式进行外方位元素的解算。 数据检核部分用前方交会算法来进行，首先根据ID找到左右片的像点的同名点，然后计算变换矩阵及系数，根据给定的内方位元素，利用公式解算出像点的控制点坐标，然后和给定的控制点坐标相比较即可进行精度评定。 (PHOTOSHOP software can be measured using a photo of the control points around the pixel coordinates can be calculated by the formula sit out the image point Standard. According to the given coordinates of control points, the use of straight-line method to calculate the linear transformation of 11 DLT parameters that the coefficient l Value and calculate the orientation elements of exterior orientation elements and solutions. Another way is to use a single formula as the calculation of space resection Orientation elements out. Use as space resection alone to settle the image exterior orientation is the basic idea: to single video As the basis of images from the ground covered by a number of control points within the known ground coordinates and image coordinates of corresponding points in the measured values , The realities of collinearity equation, find the image in the time of aerial photography exterior orientation elements. First, the process to establish a single document)

自己写了一个二维Euler方程的间断有限元程序 上次发了一个三角形单元的程序 因为不是曲边单元 所以在圆柱后面容易形成涡 现在把程序改为曲边四边形单元了 没有涡出现 单元是8节点四边形单元 节点编号顺序是 1 5 2 6 3 7 4 8 也就是四个角点依次 是1 2 3 4 然后是边的中点编号 5 6 7 8. 时间推进采用 Runge-Kutta 方法 数值通量采用全局Lax-Friedrichs通量 仍然不能捕捉激波 因为没有做重构或者加人工粘性 等这个做出来了 再发一次。 程序没有进行优化 比如说内存的消耗没有优化 比如直边单元的边界积分仍然采 用了曲边的积分方法 增加了计算量 比如面积分、线积分都是采用的是Gauss- Legendre-Lobatto积分 积分精度会比一般的Gauss-Legendre积分精度低一阶 等 等问题。 二维的 纯属交流性质 就没有考虑这些问题 ^_^ 如果物面全部是直边 那么只要改变一个参数N 就可以获得不同的计算精度 且具 有谱精度 因为单元的节点是Gauss-Legendre-Lobatto积分点。 其实就是谱元法 （物面是曲边的情况我不清楚是不是也可以通过提高基函数的阶数 也就是增加N 来提高计算精度）(Wrote a two-dimensional Euler equations with discontinuous finite element program Last made ​ ​ a triangular element of the program, not curved edge unit is so easy to form a vortex in the cylinder behind the Program to curved edge quadrilateral element vortices appear Unit is the order of 8-node quadrilateral element node number is 15,263,748 which is the four corners of the points in turn Is 1234 and then the side of the midpoint of the number 5678. Time promote the use of Runge-Kutta method Numerical flux of the overall situation of Lax-Friedrichs, flux Still can not capture the shock wave did not do the reconstruction or artificial viscosity do it Zaifayici. The program is not optimized for example, memory consumption is not optimized such as straight-edge boundary integral of the unit is still mining Integral method to increase the amount of computation such as surface integral with a curved edge, the line integral using the Gauss- The Legendr)

说明：  c++编程bp神经网络，用c++语言编写的bp神经网络(C++ programming bp neural network, bp neural network written in c++ language.)

四杆桁架结构的有限元分析程序，在VS2008软件中使用C++语言编写，输入为四杆桁架结构的节点坐标，材料参数，边界条件，输出为各节点的应力和应变。(Four-bar truss structure finite element analysis program, the use of C++ written in VS2008 software, input node coordinate the four-bar truss structure, material parameters, boundary conditions, the output of the stress and strain of each node.)

四阶龙格库塔法函数更新四元数程序代码,采用VC编写(Fourth-order Runge-Kutta method function is updated quaternion code)

几种计算阶乘、阶乘之和的算法，由于计算机存储问题，一般阶乘不可能算得太大，只是一种方法。(Several computing factorial, factorial, and the algorithm, due to computer storage problems, the general factorial can not be considered too great, is only a way.)

湍流模型re计算流体力学软件udf 可以用于修改UDF文件(turbulence model FLUENT)

涡格法计算机翼的涡强，对不同剖面形状，拱度等进行了比较分析。并且与书籍上对比(Vortex lattice method the computer wing of the vortex is strong, the shape of different profiles, such as camber, a comparative analysis.)

共轭梯度法求解凸二次规划，代码由C语言开发，老外编写的程序，很值得一读。(Conjugate gradient method for solving convex quadratic programming, code generated by C language development, foreigners prepared program, it is worth reading.)

利用abaqus软件建模，然后将inp文件导入flac软件中。(ABAQUS software is used to model, and then the InP file is imported into the FLAC software.)