- 积分
- 0
- 注册时间
- 2009-1-17
- 仿真币
-
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2009-2-19 10:26:18
|
显示全部楼层
来自 美国
本帖最后由 eric613 于 2009-2-20 09:57 编辑
先占位,等待进一步说明我的要求.
是这样的, 我有一个GIS的算例, 目的是要通过分析其内部电场的分布情况, 来对它的结构进行优化. 如果直接用有限元法计算的话, 因为网格剖分后系数矩阵的阶数很大, 达到30万阶左右, 这样求解线性代数方程组的计算量是很大的. 那么如果每次更改设计变量后, 都要进行如此大量的计算的话, 就很不方便了.
有一种方法叫做缩减基(reduced basis), 能降低系数矩阵的维数, 这样计算量就大大降低了. 这个reduced basis方法的核心就是: 在第一次求得线性代数方程的解后, 通过对系数矩阵进行求导, 得到缩减基矩阵. 这里是对设计变量进行求导, 所以要用到系数矩阵每个元素关于设计变量的表达式, 而不是一个直接的数值, 大家知道, 网格剖分后得到的系数矩阵其实是每个节点的坐标的累加, 而每个节点的坐标是和设计变量有关的(在本例中是这样的).
基本过程就是这样. 所以必须要有网格剖分的源程序, 才能进一步编程, 来实现对系数矩阵的求导, 从而得到缩减基矩阵, 只要得到缩减基矩阵, 剩下的工作就是解代数方程组的问题了! |
|