24小时热门版块排行榜    

查看: 1226  |  回复: 6

lijunjiexm

铁虫 (小有名气)

[求助] 请教形函数的导数在任意四边形的积分,谢谢

计算形函数的导数在任意四边形的积分,
应用格林公式:

但是我不知道这个积分最终等于多少(假设为4节点四边形单元,i=1,2,3,4). 因为 是x,y的函数,所以上式的曲线积分我不知道怎么处理x,请大家多多指教。谢谢!
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuyx_78

金虫 (小有名气)

【答案】应助回帖


soliton923(金币+1): 谢谢参与讨论~~~ 2011-11-05 12:00:15
lijunjiexm(金币+5): 谢谢您的回复,请看我后面的回复,谢谢 2011-11-05 12:32:47
由于具体的形函的多样性和四边形的多多样性,在四边形的边界上y与x的关系确定,即y是x的函数,从而得到每条边上的Ni(x,y)=Ni(x,y(x))及dy,这样计算四边开边界上线积分可以转化为计算四个一元函数的定积分的和。但这样的话每条边界上都要计算,如果剖分单元在太多使计算量大副增加。因此在很多情况下使用用数值积分去直接计算重积分。
2楼2011-11-05 10:46:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Jackie2011

木虫 (正式写手)

【答案】应助回帖


soliton923(金币+1): 谢谢参与讨论~~~ 2011-11-05 12:00:22
如果能写出四条边的线段的参数方程,就可代入继续求曲线积分啊!
学好外语,不忘数学!
3楼2011-11-05 11:10:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijunjiexm

铁虫 (小有名气)

引用回帖:
3楼: Originally posted by Jackie2011 at 2011-11-05 11:10:10:
如果能写出四条边的线段的参数方程,就可代入继续求曲线积分啊!

举个例子,矩形四边形


计算边界kl上的积分(节点按顺时针编号,i,j,k,l所以相应的形函数要重排(对应逆时针编号的情形)), 我的计算结果在粉红色的部分(帮我看看这部分哪里有问题就行)


但是正确结果应该是上面的再乘以一个1/2 !
4楼2011-11-05 12:29:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuyx_78

金虫 (小有名气)

【答案】应助回帖

lijunjiexm(金币+25): 明白,谢谢,标准单位四边形长为2,我一直记为1,这是问题所在 2011-11-05 15:00:31
引用回帖:
4楼: Originally posted by lijunjiexm at 2011-11-05 12:29:36:
举个例子,矩形四边形


计算边界kl上的积分(节点按顺时针编号,i,j,k,l所以相应的形函数要重排(对应逆时针编号的情形)), 我的计算结果在粉红色的部分(帮 ...

你采用的形函数是标准备的单位正方形[-1,1]*[-1,1]上的,而图中的单元是一般的矩形,在图中边上\eta与y间的关系是\eta=2(y-c)/(y_k-y_l),其中c=(y_k+y_l)/2,于是在积分中dy=(y_l-y_k)/2* d\eta.
5楼2011-11-05 14:14:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lijunjiexm

铁虫 (小有名气)

引用回帖:
5楼: Originally posted by xuyx_78 at 2011-11-05 14:14:24:
你采用的形函数是标准备的单位正方形[-1,1]*[-1,1]上的,而图中的单元是一般的矩形,在图中边上\eta与y间的关系是\eta=2(y-c)/(y_k-y_l),其中c=(y_k+y_l)/2,于是在积分中dy=(y_l-y_k)/2* d\eta.

是否对任意凸四边形都可以呢?我试试看看,能不能得到类似的公式

» 本帖已获得的红花(最新10朵)

6楼2011-11-05 18:44:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xuyx_78

金虫 (小有名气)

送鲜花一朵
引用回帖:
6楼: Originally posted by lijunjiexm at 2011-11-05 18:44:29:
是否对任意凸四边形都可以呢?我试试看看,能不能得到类似的公式

对任意四边形单元,你可以参考朱明权《解任意四边形板弯曲问题的样条有限元》(1987)这篇文章。
7楼2011-11-05 19:59:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 lijunjiexm 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 评正教授需要两个国家级项目,有人用子课题糊弄,结果在评审前资格公示时被举报拿下了 +22 瞬息宇宙 2024-11-12 31/1550 2024-11-15 22:02 by 潇湘之迷
[考博] 中南大学 粉末冶金国家重点实验室 陈超教授课题组拟招收 1~2位博士研究生通知 +7 中南大学-金属材 2024-11-14 18/900 2024-11-15 20:22 by 中南大学-金属材
[论文投稿] 找不到审稿人 +6 吉祥如意141913 2024-11-13 6/300 2024-11-15 16:48 by LEyblyt
[论文投稿] 职称论文 +10 xiaoqi_66 2024-11-09 13/650 2024-11-15 16:45 by sugarmei
[论文投稿] 核心期刊那个阶段最难 +4 公共茅斯 2024-11-11 4/200 2024-11-15 16:44 by sugarmei
[论文投稿] 小论文影响大论文查重吗? +3 阿道夫深度 2024-11-14 4/200 2024-11-15 16:19 by sugarmei
[基金申请] 博后面上太难了 +14 不安定因素 2024-11-14 17/850 2024-11-15 16:05 by lupbao
[基金申请] 变性人申请基金咋算性别? +21 leicahfs 2024-11-11 21/1050 2024-11-15 14:54 by 跳鱼拨刺
[基金申请] 去年七月底入站的还能申请下一批吗? +4 brightwo 2024-11-14 4/200 2024-11-15 14:18 by gazi1111
[硕博家园] reject后小感 +6 sioc-sunj 2024-11-14 8/400 2024-11-15 11:50 by 畅21
[论文投稿] 论文返修状态变成了awaiting AE recommendation 10+3 猪小耍 2024-11-13 12/600 2024-11-15 08:41 by 北京莱茵润色
[硕博家园] 课题和专业一点都不相关怎么办 +7 imt4n9 2024-11-11 10/500 2024-11-15 01:55 by 成长虫
[论文投稿] OE返修遇到expired +3 隔壁老王来了 2024-11-14 7/350 2024-11-14 20:46 by 隔壁老王来了
[硕博家园] 研究生的生活该是什么样 +4 lqy0719 2024-11-14 4/200 2024-11-14 16:45 by 阿荣喝酒
[有机交流] 同一个反应回流情况不同 20+3 1853846 2024-11-12 3/150 2024-11-14 15:42 by 太阳谷
[教师之家] 处在人生职业的分水岭 +4 otani 2024-11-13 4/200 2024-11-14 14:17 by mddzwo
[论文投稿] 爱思唯尔投稿系统里的通讯作者可以和文章里的通讯作者标注不同吗 +7 Omnissiah 2024-11-10 7/350 2024-11-12 14:07 by holypower
[硕博家园] 同步辐射 +4 小张要加油努力 2024-11-10 7/350 2024-11-11 13:47 by mpdfwxgui
[论文投稿] renewable energy 添加作者 15+3 Hebauwww12 2024-11-09 4/200 2024-11-11 08:42 by 北京莱茵润色
[基金申请] 要求论文发表日期在项目执行期内,论文发表日期是在线日期还是见刊日期 +6 Jiangnanyu1 2024-11-09 7/350 2024-11-10 17:49 by jurkat.1640
信息提示
请填处理意见