Redian新闻
>
这样数值积分怎么搞定?
avatar
这样数值积分怎么搞定?# Computation - 科学计算
k*g
1
在三角域上积分,函数是piecewise polynomial的,怎搞?
avatar
w*d
2
you can work out the exact integral.
But if you really want to use quadrature, check Ronald Cools' homepage.

【在 k****g 的大作中提到】
: 在三角域上积分,函数是piecewise polynomial的,怎搞?
avatar
c*e
3
FEM里面的东西?

【在 k****g 的大作中提到】
: 在三角域上积分,函数是piecewise polynomial的,怎搞?
avatar
c*m
4
去查有限元手册
基本就是取几个点,数值积分
如果只有piecewise poly,可以用重心坐标

【在 k****g 的大作中提到】
: 在三角域上积分,函数是piecewise polynomial的,怎搞?
avatar
k*g
5
还是你牛。。一眼看出来了
any suggestion?

【在 c*******e 的大作中提到】
: FEM里面的东西?
avatar
k*g
6
啥是重心坐标啊?

【在 c******m 的大作中提到】
: 去查有限元手册
: 基本就是取几个点,数值积分
: 如果只有piecewise poly,可以用重心坐标

avatar
k*g
7
在FEM里这个exact integral不好求吧。
我现在正在尝试用Ronald的package,可惜是fortran的,:(

【在 w**d 的大作中提到】
: you can work out the exact integral.
: But if you really want to use quadrature, check Ronald Cools' homepage.

avatar
c*e
8
很久不用差不多都忘了,回头查一下看看能不能解决.

【在 k****g 的大作中提到】
: 还是你牛。。一眼看出来了
: any suggestion?

avatar
c*m
9
你还是去查手册吧
一般讲有限元编程的书都有
这些东西都是早就有标准做法的

【在 k****g 的大作中提到】
: 啥是重心坐标啊?
avatar
r*u
10
1.Using triangle coordinates to integrate the piecewise function by hand or
symbolic calculate software.
2.But I guess is your problem is not very special, there are a lot of standard
solution available on finite element analysis textbook or somewhere else. Try
to search keywords “triangle element” “FEM” “finite element analysis”

【在 c*******e 的大作中提到】
: 很久不用差不多都忘了,回头查一下看看能不能解决.
avatar
k*g
11
可是在每个单元里不是整个一块的polynomial,是piecewise的
所以用一般的quadrature方法都得不到精确解的。

【在 c******m 的大作中提到】
: 你还是去查手册吧
: 一般讲有限元编程的书都有
: 这些东西都是早就有标准做法的

avatar
r*u
12
You can integrate over n small 整个一块的polynomial and calculate the sum of n
results.
Or slice that piecewise function and make element finer to make sliced
function over every elements 整个一块的polynomial.
I remember lame Numerical Manifold Method deals with a similar case.

【在 k****g 的大作中提到】
: 可是在每个单元里不是整个一块的polynomial,是piecewise的
: 所以用一般的quadrature方法都得不到精确解的。

相关阅读
logo
联系我们隐私协议©2024 redian.news
Redian新闻
Redian.news刊载任何文章,不代表同意其说法或描述,仅为提供更多信息,也不构成任何建议。文章信息的合法性及真实性由其作者负责,与Redian.news及其运营公司无关。欢迎投稿,如发现稿件侵权,或作者不愿在本网发表文章,请版权拥有者通知本网处理。