24小时热门版块排行榜    

查看: 1923  |  回复: 7

i小木虫

木虫 (正式写手)

[求助] 关于泛函分析,有限元(finite element method)计算的一个问题,简单也不简单。 已有1人参与

最近一直头疼一个小问题,由于一些数学方面没基础,所以压根没看懂是怎么个计算法, 所以请版里大神能指导一二。可以在这里面讨论。
具体是这样的,已知一个长方体,在底面固定不变的情况下,什么时候表面积最小,如何求解,并得到对应的函数。
这里有一篇文献,做出来了,大概说了下计算方法,但是露珠实在看不懂,求大神帮解读。
图1. 这一步还理解,做了个大概近似;
关于泛函分析,有限元(finite element method)计算的一个问题,简单也不简单。
图2. 第二步,需要泛函分析,这就不明白了,一来泛函分析没基础,二来,怎么突然冒出了delta U, g 和偏导omega?
关于泛函分析,有限元(finite element method)计算的一个问题,简单也不简单。-1
图3.第三步,越来越糊涂了,拉格朗日乘数法?不理解啊,是为了干个啥?
关于泛函分析,有限元(finite element method)计算的一个问题,简单也不简单。-2
图4. 最后,得到了希望的模拟结果图,没说通过什么软件,大家都认为是matlab。
关于泛函分析,有限元(finite element method)计算的一个问题,简单也不简单。-3

可能这个问题有些复杂和繁琐,1.理论计算过程如何理解,比如应用到正六角柱,即上下底面都是正六边形,固定底面也可以使用该方法吗?2.如何在软件中操作,比如MATLAB,总不能他里面有Δu,我直接也写个就行?

欢迎大家能参与讨论,如果问题解决露珠保证追加金币!!!

问的问题可能有些已知量没交代清楚,具体在下面这个文献里面。
参考文献:Liu, X., et al. (2015). Fabrication and surface profile simulation of sapphire microlens array. Selected Proceedings of the Photoelectronic Technology Committee Conferences held August-October 2014, International Society for Optics and Photonics.
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gold2007

捐助贵宾 (正式写手)

从原泛函极小值问题导出微分方程的边值问题,是通过求泛函的变分导数得到的。参见图片上大学四年级课本的推导。精确的方程叫做极小曲面方程,作了线性近似就是Laplace方程。
关于泛函分析,有限元(finite element method)计算的一个问题,简单也不简单。-4


关于泛函分析,有限元(finite element method)计算的一个问题,简单也不简单。-5



发自小木虫IOS客户端
2楼2016-06-05 16:27:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gold2007

捐助贵宾 (正式写手)

在V等于常数的约束条件下求极值,正如第三步那样,引入乘子去掉约束;原理参见任一本高等数学,其严格证明需参考专著。

发自小木虫IOS客户端

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

3楼2016-06-05 16:45:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gold2007

捐助贵宾 (正式写手)

引用回帖:
3楼: Originally posted by gold2007 at 2016-06-05 16:45:10
在V等于常数的约束条件下求极值,正如第三步那样,引入乘子去掉约束;原理参见任一本高等数学,其严格证明需参考专著。

因为函数空间是无限维的,超出了高等数学教材范围。

发自小木虫IOS客户端
4楼2016-06-05 16:46:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

peterflyer

木虫之王 (文学泰斗)

peterflyer


【答案】应助回帖

感谢参与,应助指数 +1
第二步系由受约束的二维泛函的变分法中的欧拉方程而推得,最后欧拉方程转化为了一个调和微分方程。求解该定解问题即可得到体积一定时使表面积最小化的曲面形状函数。
第三步的来由是受约束问题目标泛函的拉格朗日乘子法。这可再去复习一下高数求极值那一相关章节。
此处的g是何物,本人也不清楚。希望楼主更详细地提供题目详情。
5楼2016-06-05 18:03:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

i小木虫

木虫 (正式写手)

送红花一朵
引用回帖:
3楼: Originally posted by gold2007 at 2016-06-05 16:45:10
在V等于常数的约束条件下求极值,正如第三步那样,引入乘子去掉约束;原理参见任一本高等数学,其严格证明需参考专著。

谢谢你,突然发现知道的太少,还不是很清楚如何能在Matlab上操作,我在下面在学习下。不知道为什么金币发不过去的,先送一束鲜花给你。
6楼2016-06-06 09:43:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

i小木虫

木虫 (正式写手)

引用回帖:
5楼: Originally posted by peterflyer at 2016-06-05 18:03:34
第二步系由受约束的二维泛函的变分法中的欧拉方程而推得,最后欧拉方程转化为了一个调和微分方程。求解该定解问题即可得到体积一定时使表面积最小化的曲面形状函数。
第三步的来由是受约束问题目标泛函的拉格朗日乘 ...

应该是一楼公式2.6中的边界条件,u▏偏导Ω=φ(x,y)
7楼2016-06-06 09:46:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gold2007

捐助贵宾 (正式写手)

引用回帖:
6楼: Originally posted by i小木虫 at 2016-06-06 09:43:56
谢谢你,突然发现知道的太少,还不是很清楚如何能在Matlab上操作,我在下面在学习下。不知道为什么金币发不过去的,先送一束鲜花给你。...

不客气,相互学习,三人行必有我师。
我不在意金币。

发自小木虫IOS客户端
8楼2016-06-06 10:06:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 i小木虫 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:48 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 xvJS1BYDxvBS 2026-10-31 8/400 2026-11-01 05:43 by YXM1FDCSE8Bn
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[公派出国] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 05:41 by YXM1FDCSE8Bn
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 05:22 by YXM1FDCSE8Bn
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 1uaZt2kAZj2S 2026-10-30 14/700 2026-11-01 05:17 by YXM1FDCSE8Bn
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 1uaZt2kAZj2S 2026-10-31 6/300 2026-11-01 03:56 by YXM1FDCSE8Bn
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-30 10/500 2026-11-01 03:51 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 15/750 2026-11-01 03:51 by YXM1FDCSE8Bn
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 1uaZt2kAZj2S 2026-10-30 13/650 2026-11-01 03:46 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 X5bHxOULxFqb 2026-10-30 16/800 2026-11-01 03:38 by YXM1FDCSE8Bn
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 X5bHxOULxFqb 2026-10-30 20/1000 2026-11-01 03:33 by YXM1FDCSE8Bn
[论文投稿] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 X5bHxOULxFqb 2026-10-30 25/1250 2026-11-01 03:30 by YXM1FDCSE8Bn
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 1uaZt2kAZj2S 2026-10-31 5/250 2026-11-01 02:25 by YXM1FDCSE8Bn
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 X5bHxOULxFqb 2026-10-30 12/600 2026-11-01 00:15 by YXM1FDCSE8Bn
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 X5bHxOULxFqb 2026-10-30 13/650 2026-10-31 17:12 by hKe7U1kad5Sg
[基金申请] 论文的挂名怎么评判贡献?基金申请中,排名第一的挂名文章 +3 biily 2026-10-25 3/150 2026-10-30 18:33 by 旁观2020
[硕博家园] 申博自荐 2篇一区TOP+1篇EI+2公开发明专利 2027博士申请 +3 小小阿光 2026-10-27 3/150 2026-10-29 23:56 by 起风了的快乐
[文学芳草园] 闲聊 +3 myrtle 2026-10-25 6/300 2026-10-29 15:28 by myrtle
信息提示
请填处理意见