πŸ”
brainless
  • What is this?
  • Responsibility
  • Changelog
  • meta
    • Sharing
      • Inspirations
      • Workflows
      • Social Media
    • Geography
      • Life
      • Death
        • Family Death
    • Research
      • Project Index
      • 3D Printing
      • Photogrammetry
      • Drone Building
    • External Websites
    • [unreleased]
      • [Template]
      • [TEMP] CL V0.0.5 or whatever
      • Skincare
      • Travel
      • Working and Staying Busy
      • Stride
      • Funeral Playlist
      • Notes and Ideas
      • Boredom
      • Four Noble Truths of "Thermo"
      • Respect
      • Work
  • STE[A]M
    • [guide]
    • Science
      • Materials Modeling
        • Syllabus Description
        • Lecture Slides
        • Student Notes
        • Assignments
          • index
          • Carbon Nanotubes
    • Technology
      • Computer Science
        • Commands
      • Photogrammetry
      • Quantum Computing
      • Computers
      • Programs
        • Matlab and Octave
        • Audacity
        • Google Chrome
          • Websites
            • Google Suite Sites
            • Github
              • Version Control
            • Product Hunt
            • Twitter
            • Youtube
              • Channels
            • Vimeo
          • Extensions
            • Dark Reader
            • Vimium
        • [miscellaneous]
          • Octave
          • PureRef
          • git
          • gnu stow
          • mermaid.js
        • Excel
        • Blender
        • LaTeX
        • Sublime: Text Editor
        • Spotify
        • VLC Media Player
      • Android and iOS
      • Operating Systems
        • macOS
          • mackup
        • Unix
          • folder structure
        • Windows
          • App Installation
          • Meshroom
          • Corsair Utility Engine
      • 3D Printing
    • Engineering
      • Accreditation
        • Fundamentals of Engineering
        • Professional Engineering
      • Continuum Mechanics
        • Fluid Mechanics
          • Incompressible Flow
            • corona final
            • indexhw4
            • indexhw3
            • index
            • hw2
            • hw1
          • Syllabus Description
          • Lecture Slides
          • Student Notes
            • Dynamic or Kinematic Viscosity
          • Assignments
            • Homeworks
              • Homework 1
              • Homework 2
              • Homework 3
              • Homework 4
              • Homework 5
            • Vortex Project
        • Solid Mechanics
          • Syllabus Description
          • Lecture Slides
          • Student Notes
          • Assignments
        • Incompresible Flow
          • Syllabus Description
          • Lecture Slides
          • Student Notes
          • Assignments
      • Experimental Mechanics
        • Syllabus Description
        • Lecture Slides
        • Student Notes
        • Assignments
      • Finite Element Methods
        • Intro to Finite Elements
          • Syllabus Description
          • Lecture Slides
          • Student Notes
          • Assignments
        • Fundamentals of FEM
          • Syllabus Description
            • index
          • Lecture Slides
          • Student Notes
          • Assignments
            • Project
              • index
              • Untitled
            • Homework 1
            • Homework 4
            • index
      • Heat Transfer
        • Syllabus Description
        • Lecture Slides
        • Student Notes
        • Assignments
          • homework
            • hw10
            • q9
            • q8
            • q7
            • hw7
            • hw6
            • q5
            • q3
            • 1 ec
          • Discussions
            • d11
            • d10
            • d9
            • d8
            • d6
            • d4
            • d3
          • Project Notes
      • Machine Dynamics
        • Syllabus Description
        • Lecture Slides
        • Student Notes
        • Assignments
    • Art
      • Color Theory
      • Origami
        • FolderMath
          • Surveying Origami Math
          • Represent a Folded Object
          • Creating a Crease Pattern
          • Making the Folds
          • Simulating Folding Origami
          • List of Resources
            • Codes
            • Papers, Programs, and Inspirations
    • Mathematics
      • Complex Numbers
        • What is i^i?
      • Analytic Hierarchy Process
      • Probability
      • Conway's Game of Life
      • Metallic Numbers
      • Cauchy's Formula for Repeated Integration
      • Wavelet Transform
      • Laplace Tidal Equation
      • Alternating Summation of Ones
      • Constants
      • Bad Maths
      • Calculus
        • Syllabus Description
        • Miscellaneous
  • Thoughts
    • Marksmanship
      • Archery
    • Schooling
    • ...and Ideas?
      • Perceived Time and Learning
      • Content Comprehension
    • Comics and Games
      • Rubik's Cube
      • Dungeons and Dragons
      • Beyond-All-Reason
      • Sekiro: Shadows Die Twice
      • Super Smash Bros
        • Project M
        • Project +
      • League of Legends
      • Satisfactory
    • Literature and Art
      • Books
      • Reading is Hard
      • Various Words and Phrases
      • Poems
      • Interviews
      • Quotes
        • Phrases
      • Jokes
      • ASCII Art
    • Shows and Films
      • Cowboy Bebop
      • My Hero Academia
      • Sword of the Stranger
    • Working and Life Balance
  • Projects
Powered by GitBook
On this page
  • Problem 1
  • Symmetry Condition
  • Bilinearity Condition
  • Problem 2
  • Part a
  • Part b
  • Part c
  • Part d
  • Problem 3
  1. STE[A]M
  2. Engineering
  3. Finite Element Methods
  4. Fundamentals of FEM
  5. Assignments

Homework 1

aaayyy

PreviousUntitledNextHomework 4

Last updated 4 years ago

Chris Nkinthorn

Finite Element Mechanics

Due Sept 17, 2019

Problem 1

Does the given function for $a(v, u)$ satisfy the symmetry and bilinearity conditions. Note $\kappa$ and $\lambda$ that and are given functions.

a ( v , u ) = \int _ { - 1 } ^ { 1 } \left( v _ {, x } \kappa u _ {, x } + v \lambda u \right) d x \label{prob1}

Symmetry Condition

Statement of this condition involves showing that the function $a$ can accept the functions $v$ and $u$ in either order. This is done by substituting each function in place of the other. Then employing algebra to form the arrangement of $[\ref{prob1}]$. So, as the arguments $u$ and $v$ are treated similarly: when multiplying $\lambda$, $u$ and $v$ are directly substituted for the other. When multiplying by $\kappa $, both functions have have derivatives taken with respect to $x$ :

\begin{align} a(u,v) = a(v,u) \quad \longrightarrow \quad \int _ { - 1 } ^ { 1 } \left( v _ {, x } \kappa u _ {, x } + v \lambda u \right) d x &= \int _ { - 1 } ^ { 1 } \left( \overbrace{ u _ {, x } \kappa v _ {, x } + u \lambda v }^{\text{flip w/ associative prop.}} \right) d x \label{prob1b} \\ & = \int _ { - 1 } ^ { 1 } \left( v _ {, x } \kappa u _ {, x } + v \lambda u \right) d x \quad \boxed{\checkmark} \label{prob1c} \end{align}.

Comparison of the lefthand side of Eqn. $[\ref{prob1b}]$ to the expression of Eqn. $[\ref{prob1c}]$, we prove that $a(u,v) = a(v,u) $.

Bilinearity Condition

When testing if the function $a(u,v)$ obeys bilinearity, proceed with the following the substitution the function $w$ in place of $u$. Now, with this variable $u$ is free, substitute in place of $v$, two scaled functions expressed as $c_1 u$ and $c_2 v$. The variables $c_1$ and $c_2$ are scaling coefficients. No offset coefficient is needed thanks to the assumption of positive definiteness.

Let’s attempt this:

By applying the linearity of the process integration to each term, we can group terms by the scaling coefficients. This separates into two terms which are each in analogous form.

Problem 2

For the following functions indicate if it is positive-definite, positive semi-definite, or neither ($\alpha$ is a scalar). Be sure to explain your answer.

Part a

Above is a function which is positive definite because the integrand is also positive definite. The integrand is the sum of two terms which are only ever positive. This is definite because the expression will only ever be zero when arguments are zero, thanks to the fixed interval $[-1,1]$

Part b

{ \left( \int _ { - 1 } ^ { 1 } \left( u ^ { 2 } + u _ { x x } ^ { 2 } \right) d x \right) ^ { 3 } } \label{prob2b}

The function above is neither positive definite nor positive semi-definite as the expression given is raised to an odd power, which means that the expression may be either positive or negative.

Part c

This is a positive semidefinite function as the variable $\alpha$ is an unspecified scalar, ranging anywhere along the real number line. If this condition is fulfilled, then the function will then definite-ly be positive.

Part d

Finally, the last component is another positive definite function. This follows similar logic as $[\ref{prob2b}]$ but with an even power (the squaring of the integral) restricting the output to $\mathbb{R}^+$.

Problem 3

In on the order of 1/2 page, indicate why you are interesting in finite element methods (including the application you would apply them to) and what you hope to get out of this course.

Finite Element Methods is fundamental to computational modeling. Having spent many hours in the SCOREC conference room listening to everyone talk to the Sandia researchers and having put it off from last year, it’s time that I get to peek under the mathematical hood of how modeling works. Perhaps not build up all the way to the multiphysics processing, similar to some modeling done in Albany; still, I would like to convert some FEM pseudo-code into a viable working code for some given set of conditions. That’s a skill which will prove valuable in the future, no matter what specific discipline I end up in. The stretch goal for this course is to be able to do some of this from memory: having the ability to derive these from wrote is a pretty good indication of mastery.

a(c1u+c2v,w)=c1a(u,w)+c2a(v,w)⟢=βˆ«βˆ’11((c1u+c2v),xΞΊw,x+(c1u+c2v)Ξ»w)dx=βˆ«βˆ’11((u,xΞΊw,x+uΞ»w)+c2(v,xΞΊw,x+vΞ»w))dx=βˆ«βˆ’11(c1(u,xΞΊw,x+uΞ»w))dx+βˆ«βˆ’11(c2(v,xΞΊw,x+vΞ»w))dx=c1βˆ«βˆ’11(u,xΞΊw,x+uΞ»w)dx⏟a(u,w)+c2βˆ«βˆ’11(v,xΞΊw,x+vΞ»w)dx⏟a(v,w).=c1a(u,w)+c2a(v,w)βœ“\begin{align} a \left( c _ { 1 } u + c _ { 2 } v , w \right) = c _ { 1 } a ( u , w ) + c _ { 2 } a ( v , w ) \quad \longrightarrow \quad & = \int _ { - 1 } ^ { 1 } \left( \left( c _ { 1 } u + c _ { 2 } v \right) _ {, x } \kappa w _ {, x } + \left( c _ { 1 } u + c _ { 2 } v \right) \lambda w \right) d x \\ &= \int _ { - 1 } ^ { 1 } \left( \left( u _ {, x } \kappa w _ {, x } + u \lambda w \right) + c _ { 2 } \left( v _ {, x } \kappa w _ {, x } + v \lambda w \right) \right) d x \\ &= \int _ { - 1 } ^ { 1 } \left( c _ { 1 } \left( u _ {, x } \kappa w _ {, x } + u \lambda w \right) \right) dx + \int_{-1}^1 \left( c _ { 2 } \left( v _ {, x } \kappa w _ {, x } + v \lambda w \right) \right) d x \\ &= c _ { 1 } \underbrace{ \int _ { - 1 } ^ { 1 } \left( u _ {, x } \kappa w _ {, x } + u \lambda w \right) dx }_{a(u,w)} + c _ { 2 } \underbrace{ \int_{-1}^1 \left( v _ {, x } \kappa w _ {, x } + v \lambda w \right) d x }_{a(v,w)}\quad. \\ & = c _ { 1 } a ( u , w ) + c _ { 2 } a ( v , w ) \quad \boxed{\checkmark} \end{align}a(c1​u+c2​v,w)=c1​a(u,w)+c2​a(v,w)βŸΆβ€‹=βˆ«βˆ’11​((c1​u+c2​v),x​κw,x​+(c1​u+c2​v)Ξ»w)dx=βˆ«βˆ’11​((u,x​κw,x​+uΞ»w)+c2​(v,x​κw,x​+vΞ»w))dx=βˆ«βˆ’11​(c1​(u,x​κw,x​+uΞ»w))dx+βˆ«βˆ’11​(c2​(v,x​κw,x​+vΞ»w))dx=c1​a(u,w)βˆ«βˆ’11​(u,x​κw,x​+uΞ»w)dx​​+c2​a(v,w)βˆ«βˆ’11​(v,x​κw,x​+vΞ»w)dx​​.=c1​a(u,w)+c2​a(v,w)βœ“β€‹β€‹β€‹
βˆ«βˆ’11(ux2+uxx2)dx\int _ { - 1 } ^ { 1 } \left( u _ { x } ^ { 2 } + u _ { x x } ^ { 2 } \right) d xβˆ«βˆ’11​(ux2​+uxx2​)dx
∫01Ξ±u2dx{ \int _ { 0 } ^ { 1 } \alpha u ^ { 2 } d x }∫01​αu2dx
(∫01Ξ±u2dx)2{ \left( \int _ { 0 } ^ { 1 } \alpha u ^ { 2 } d x \right) ^ { 2 } }(∫01​αu2dx)2
153KB
mane.6660.hw1.pdf
pdf
Homework 1 Prompt