24小时热门版块排行榜    

北京石油化工学院2025年硕士研究生招生接受调剂公告
查看: 1366  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 华南师范大学生物光学材料课题组招收调剂(统考数学,07,08大类不限专业) +6 llss0711 2025-04-03 12/600 2025-04-04 23:23 by 某十七
[考研] 0703化学352求调剂 +10 yqwerfdfjj 2025-04-02 11/550 2025-04-04 23:15 by 河大薇薇
[考研] 求调剂 +3 王世宁 2025-04-04 4/200 2025-04-04 22:46 by 诗与自由
[考研] 材料352找调剂 +3 沧海丶一粟 2025-04-04 3/150 2025-04-04 22:07 by 王元明
[考研] 求调剂 +5 cdw987 2025-04-04 6/300 2025-04-04 20:12 by yanenwang
[考研] 一志愿985 生物与医药调剂 +3 荔枝蛋糕 2025-04-03 6/300 2025-04-04 10:12 by 知更鸟6816
[考研] 材料与化工求调剂 +12 cs0106 2025-04-03 14/700 2025-04-04 09:58 by hongcg1977
[考研] 293求调剂 +17 mao1202 2025-04-02 17/850 2025-04-04 05:33 by 呼大怪
[考研] 求调剂 +9 蓝桉羽都 2025-04-02 10/500 2025-04-03 20:37 by asdhh1991
[考研] 273求调剂 +10 攀登T1 2025-04-03 10/500 2025-04-03 16:48 by mini小九
[考研] 315求调剂(11408) +5 19556298861 2025-04-03 5/250 2025-04-03 14:59 by 鱼遇雨欲语于余�
[考研] 求调剂 一志愿南理工材料专硕 307分 +5 QYJ1235 2025-03-30 5/250 2025-04-03 13:41 by Ade.
[考研] 本科211,11408 321分,求调剂 +5 无用一书生 2025-04-02 9/450 2025-04-03 07:57 by 无用一书生
[考研] 材料考研336分求调剂 +7 张天傲121 2025-04-02 7/350 2025-04-03 04:20 by alexanderp
[考研] 求调剂 +12 糖糖糖lu 2025-03-31 13/650 2025-04-03 00:55 by 诗与自由
[考研] 299求调剂 +7 LLLLcx618 2025-03-31 7/350 2025-04-02 23:26 by 诗与自由
[论文投稿] 求一个Ti2AlNb小论文辅导的 +3 reveve 2025-03-30 3/150 2025-04-02 11:28 by cftgs12
[考研] 070304(物理化学)309求调剂 +7 T&J_Tom 2025-04-01 9/450 2025-04-02 09:45 by kxxmaqq
[考研] 313求调剂 +4 Pj. 2025-03-30 4/200 2025-03-30 17:35 by jakeatcug
[考研] 土木工程学硕调剂 +4 哈哈哈哈哈哈哈� 2025-03-29 4/200 2025-03-29 21:34 by jiangkaiyu
信息提示
请填处理意见
婢跺秷鐦拫鍐ㄥ婢跺秷鐦拫鍐ㄥ
婢跺秷鐦拫鍐ㄥ娴滃瞼娣惍锟�閼惧嘲褰囨径宥堢槸 & 鐠嬪啫澧忛幐鍥у础閿涳拷瀵邦喕淇婇幍顐$閹碉拷
閸熷棗濮熼崥鍫滅稊閸熷棗濮熼崥鍫滅稊
閸熷棗濮熼崥鍫滅稊娴滃瞼娣惍锟�閸熷棗濮熼崥鍫滅稊閸溿劏顕�瀵邦喕淇婇幍顐$閹碉拷(婢跺洦鏁為敍姘櫌閸斺€虫値娴o拷)
鐎涳附婀抽獮鑼舵彛鐎涳附婀抽獮鑼舵彛
鐎涳附婀抽獮鑼舵彛娴滃瞼娣惍锟�瀵邦喕淇婇幍顐$閹碉拷娴滃棜袙閺堚偓閺傛澘顒熼張顖氬叡鐠愶拷
鐠佺儤鏋冩潏鍛嚤鐠佺儤鏋冩潏鍛嚤
鐠佺儤鏋冩潏鍛嚤娴滃瞼娣惍锟�瀵邦喕淇婇幍顐$閹碉拷濞h濮為懓浣哥瑎閿涘苯鍘ょ拹鐟版尒鐠囷拷
鐢牐绁幏娑樺珯鐢牐绁幏娑樺珯