hello
I am new to mechanical modeling of tissue. So, please forgive me if I say something obviously wrong.
what I want to do is to compute an "internal pressure" distribution for soft-tissue under compression, specifically, breast under mammography compression. I wrote a simple mass-spring program to comptue the deformation of a breast by supplying Hooke coefficients for each spring. With one step compression, the program comptue the steady-state node positions by solving the force-equilibrium equation at each node with an Euler's method.
now, from the steady-state, how can I extract such an "internal pressure"? I can compute the summation of all force amplitudes at each node, or the total force amplitudes devided by spring number connected to this node (i.e. averaged force per spring). But I can not do vector summation, because they are all zeros for internal nodes. How would the (averaged) total force amplitudes correlate to the internal pressure?
a friend suggested me to compute a stress for each spring, and display the distribution of these stresses (z-axis component). But I found this is spring-connection-topology dependent, the plot looks funny and does not make sense.
I am wondering if anyone can give me some suggestions to compute this distribution .
thanks
Qianqian