SfePy

Download

sfepy@python.org

January 2024

3 participants
2 discussions

Question about using sfepy
by Hongwei Wu Jan. 17, 2024

Jan. 17, 2024

Hi, everyone. I am starting to solve a chemomechanical problem, but I don't know whether the fully coupled equations can be solved in sfepy or not. Could you have a look at them? The equation are as followed: Governing Equations ∇∙σ+ρb=ρu ̈ Where σ is the stress tensor, ρ is the mass density, b is the body force per unit mass, u is the displacement vector. ∂c/∂t+∇∙j=∂r/∂t, where j=-(D*c_0)/RT ∇Ψ Where c is the concentration, D is the diffusion coefficient, c0 is the reference concentration, R is Boltzmann constant (gas constant), T is the absolute temperature, Ψ is the chemical potential. Constitutive relations σ=2Gε+(λe+bc)I Ψ=be+gc Where b=-β(2G+3λ), β is the chemical dilation coefficient, e is the volumetric strain and g represents the scalar chemistry modulus.

2 1

pseudo-2D Geometry to current flow
by Steve Spicklemire Jan. 4, 2024

Jan. 4, 2024

Hello SfePy folks, It’s been a long time (e.g., 10+ years) since I tinkered with SFEPY, but I dug out my old calculations and with a bit of tweaking, they still work. Yeah! Thank you for keeping on for so long with this excellent tool for python! I’m hoping you can lend me some guidance for a new calculation. I’m interested in using SfePy to answer questions about the behavior of PCB traces. For example, I’d like to be able to extract a trace from a Gerber file and convert it into a mesh that SfePy can grok. I can figure out the file conversion, I think. If nothing else I can fake it using the gmsh “geo” format. However, then I’d like to compute the current flow through the mesh with a fixed voltage difference between two boundaries while maintaining a zero current boundary condition on all the other boundaries. Is that a thing that could be done? Any suggestions on how to go about setting up the boundary conditions? Thank you, -Steve

3 10