Hello I want to solve the following simple problem in 2D solid mechanics: A rectangular uniform mesh with a rectangular bar embedded inside at any arbitrary angle. The bar is represented by a density field for each element. The density field takes a value of 1 where the bar is, and 1e-5 (or a near-zero value to avoid singularities) where the domain is empty. I will simply apply a uniaxial tension to the bar. No dirichlet conditions are required. I am writing here to see if this has done before with the libMesh framework. The geometry is currently defined. I am going to apply the uniaxial forces by having a condition on the elements near the geometrical boundary. Even if that worked, I would still have to remove the null spaces from the matrix. How could I calculate this for a bar embedded in a uniform mesh?
I imagine that this problem has been solved before, hopefully in libMesh, any help will be appreciated. Thanks Miguel ------------------------------------------------------------------------------ Find and fix application performance issues faster with Applications Manager Applications Manager provides deep performance insights into multiple tiers of your business applications. It resolves application problems quickly and reduces your MTTR. Get your free trial! https://ad.doubleclick.net/ddm/clk/302982198;130105516;z _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
