# Beam bending using Gridap

**URL:** https://discourse.julialang.org/t/beam-bending-using-gridap/112723
**Category:** New to Julia
**Created:** [April 9, 2024, 12:40pm UTC](https://discourse.julialang.org/t/beam-bending-using-gridap/112723 "2024-04-09T12:40:59Z")
**Posts on this page:** 1
**Page:** 1

<div class="post-metadata">

### Author: ![Mishal](https://sea2.discourse-cdn.com/julialang/user_avatar/discourse.julialang.org/mishal/32/214462_2.png) [@Mishal](https://discourse.julialang.org/u/Mishal)
#### Post date: [April 9, 2024, 12:40pm UTC](https://discourse.julialang.org/t/beam-bending-using-gridap/112723/1 "2024-04-09T12:40:59Z")

</div>

Hi Guys

I would appreciate any advice on the below code. I am struggling to get it up and running. It is cantilevered beam which I am trying to solve.

```julia
using Gridap

𝒯 = CartesianDiscreteModel((0,1),(5)) 
writevtk(model,"model") 

refFE = ReferenceFE(lagrangian,Float64,1)

V = TestFESpace(𝒯,refFE,dirichlet_tags=[1])

Q = FESpace(model, ReferenceFE(lagrangian,Float64,1))

uD = (0.1)

U = TrialFESpace(V,uD)
P = TrialFESpace(Q)

E = 200e9
I = 1e6

w = 10
# RHS(u,v) = ∫( w*v )dΩ
# LHS(u,v) = E*I

RHS(u,v) = ∫( w*v )dΩ # This is incorrect, need to define properly
LHS(u,v) = E*I * ∫(∂(∂(u)) * ∂(∂(v))) # This is also incorrect, defining the weak form properly

op = AffineFEOperator(RHS,LHS,U,V,uD)

```

![weak form](https://global.discourse-cdn.com/julialang/original/3X/d/5/d573b46fedaa47ce4afe56631935be924b75308a.jpeg)