Warping function for 3D beams
Dear All,
I am modeling 3D beams with arbitrary cross sections and for that purpose I need to calculate the warping function numerically. The warping function ψ(y,z) is defined as a solution of the following Laplace equation with Neumann boundary conditions:
∂2ψ/∂y2 + ∂2ψ/∂z2 = 0 in Ω,
∂ψ/∂n = z*ny - y*nz on Γ,
where Ω is the cross sectional area, and Γ is its contour. This system can be solved by different numerical methods, for example, finite element method, boundary element method, finite differences...
Numerically I obtain that
Hi All!
My problem is quite simple:
I need to change (add arguments) to a inp file, then submit the inp file
for analysis and finally to
read some results from the odb file.
I choose to use python (not the PDE provided by abaqus CAE but the
general python GUI). I’m sure this choice has consequences, perhaps the
import abaqus won’t work? This could bypassed if I run my code using abaqus python in the command window.