2021-5-17 有限元方法从0开始学习笔记 第一天

第一天

从今天开始从0学习有限元方法, 由于记纸质笔记太麻烦, 字也比较丑, 无法静下心来, 因此选择在博客上记录笔记, 每天能看多少记录多少, 虽然速度可能会慢一点, 但希望能够静下心来, 坚持学习.
参考资料
Brenner S, Scott R. The mathematical theory of finite element methods[M]. Springer Science & Business Media, 2007.

(0.1 节)

弱形式

考虑边值问题:
{ − d 2 u d x 2 = f ∈ ( 0 , 1 ) u ( 0 ) = 0 ,   u ′ ( 1 ) = 0 . \begin{cases} -\dfrac{d^{2}u}{dx^{2}} = f \in (0, 1) \\ u(0) = 0, \ u^{'}(1) = 0 \end{cases}. dx2d2u=f(0,1)u(0)=0, u(1)=0.
定义函数空间:
V = H ⋄ 1 ( 0 , 1 ) = { u ∈ H 1 ( 0 , 1 ) ∣ u ( 0 ) = 0 } . V = H^{1}_{\diamond}(0, 1) = \left\lbrace \left. u \in H^{1}(0, 1) \right| u(0) = 0 \right\rbrace. V=H1(0,1)={ uH1(0,1)u(0)=0}.
下面假设需要在空间 V V V 上找解 u u u ,将上面两边同事乘以 v ∈ V v \in V vV 并且进行分部积分, 得到如下弱问题:
find  u ∈ V   s . t . ,   a ( u , v ) = ( f , v )   ∀ v ∈ V .   ( WP ) \text{find} \ u \in V \ s.t., \ a(u, v) = (f, v) \ \forall v \in V. \ \color{red}{(\text{WP}}) find uV s.t., a(u,v)=(f,v) vV. (WP)
其中
a ( u , v ) = ∫ 0 1 u ′ ( x ) v ′ ( x ) d x , ( f , v ) = ∫ 0 1 f ( x ) v ( x ) d x . a(u, v) = \int_{0}^{1}u'(x)v'(x)dx, \\ (f,v) = \int_{0}^{1} f(x)v(x) dx. a(u,v)=01u(x)v(x)dx,(f,v)=01f(x)v(x)dx.
下面给出以下定理:
定理: 假设在问题 ( WP ) \color{red}{(\text{WP})} (WP)中有 f ∈ C 0 ( [ 0 , 1 ] ) f \in C^{0}([0, 1]) fC0([0,1]) 以及 u ∈ C 2 ( [ 0 , 1 ] ) u \in C^{2}([0, 1]) uC2([0,1])则, u u u 是原来边值问题的解.
证明:
v ∈ V ∩ C 1 ( [ 0 , 1 ] ) v \in V \cap C^{1}([0, 1]) vVC1([0,1])(为了保证能够在连续意义下进行分部积分.), 对问题 ( WP ) \color{red}{(\text{WP})} (WP)进行分部积分, 得:
( f , v ) = a ( u , v ) = ∫ 0 1 − u ′ ′ v d x + u ′ ( 1 ) v ( 1 ) . (f, v) = a(u, v) = \int_{0}^{1} - u^{''}v dx + u^{'}(1)v(1). (f,v)=a(u,v)=01uvdx+u(1)v(1).
代入右端边界条件, 移项得到:
( f + u ′ ′ , v ) = 0 ,   ∀   v ∈ V ∩ C 1 ( [ 0 , 1 ] ) & v ( 1 ) = 0. (f + u^{''}, v) = 0, \ \forall \ v \in V \cap C^{1}\left([0, 1]\right) \& v(1) = 0. (f+u,v)=0,  vVC1([0,1])&v(1)=0.
下面证明 w = f + u

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值